Gene duplication is a major mechanism to create new genes. After gene duplication, some duplicated genes undergo functionalization, whereas others largely maintain redundant functions. Duplicated genes comprise various degrees of functional diversification in plants. However, the evolutionary fate of high and low diversified duplicates is unclear at genomic scale. To infer high and low diversified duplicates in Arabidopsis thaliana genome, we generated a prediction method for predicting whether a pair of duplicate genes was subjected to high or low diversification based on the phenotypes of knock-out mutants. Among 4,017 pairs of recently duplicated A. thaliana genes, 1,052 and 600 are high and low diversified duplicate pairs, respectively. The predictions were validated based on the phenotypes of generated knock-down transgenic plants. We determined that the high diversified duplicates resulting from tandem duplications tend to have lineage-specific functions, whereas the low diversified duplicates produced by whole-genome duplications are related to essential signaling pathways. To assess the evolutionary impact of high and low diversified duplicates in closely related species, we compared the retention rates and selection pressures on the orthologs of A. thaliana duplicates in two closely related species. Interestingly, high diversified duplicates resulting from tandem duplications tend to be retained in multiple lineages under positive selection. Low diversified duplicates by whole-genome duplications tend to be retained in multiple lineages under purifying selection. Taken together, the functional diversities determined by different duplication mechanisms had distinct effects on plant evolution.

BackgroundAmong human T cell leukemia virus type 1 (HTLV-1)-infected individuals, there is an association between HTLV-1 tax subgroups (subgroup-A or subgroup-B) and the risk of HAM/TSP in the Japanese population. To investigate the role of HTLV-1 subgroups in viral pathogenesis, we studied the functional difference in the subgroup-specific viral transcriptional regulators Tax and HBZ using microarray analysis, reporter gene assays, and evaluation of viral-host protein–protein interaction.Results(1) Transcriptional changes in Jurkat Tet-On human T-cells that express each subgroup of Tax or HBZ protein under the control of an inducible promoter revealed different target gene profiles; (2) the number of differentially regulated genes induced by HBZ was 2–3 times higher than that induced by Tax; (3) Tax and HBZ induced the expression of different classes of non-coding RNAs (ncRNAs); (4) the chemokine CXCL10, which has been proposed as a prognostic biomarker for HAM/TSP, was more efficiently induced by subgroup-A Tax (Tax-A) than subgroup-B Tax (Tax-B), in vitro as well as in unmanipulated (ex vivo) PBMCs obtained from HAM/TSP patients; (5) reporter gene assays indicated that although transient Tax expression in an HTLV-1-negative human T-cell line activated the CXCL10 gene promoter through the NF-κB pathway, there was no difference in the ability of each subgroup of Tax to activate the CXCL10 promoter; however, (6) chromatin immunoprecipitation assays showed that the ternary complex containing Tax-A is more efficiently recruited onto the promoter region of CXCL10, which contains two NF-κB binding sites, than that containing Tax-B.ConclusionsOur results indicate that different HTLV-1 subgroups are characterized by different patterns of host gene expression. Differential expression of pathogenesis-related genes by subgroup-specific Tax or HBZ may be associated with the onset of HAM/TSP.Electronic supplementary materialThe online version of this article (10.1186/s12977-018-0454-x) contains supplementary material, which is available to authorized users.

Lineage-specific gene duplications contribute to a large variation in specialized metabolites among different plant species. There is also considerable variability in the specialized metabolites within a single plant species. However, it is unclear whether copy number variations (CNVs) derived from gene duplication events contribute to the diversity of specialized metabolites within species. We identified metabolome quantitative trait genes (mQTGs) associated with quantitative metabolite variations and examined the relationship between mQTGs and CNVs. We obtained 1,335 specialized metabolite signals from 53 worldwide A. thaliana accessions using liquid chromatography-quadrupole time-of-flight mass spectrometry. In this study, genes associated with specialized metabolites were inferred by either a generally authorized genome-wide association study (GWAS) approach or a novel analysis of the association between gene expression and metabolite accumulation. Genes qualified by both analyses are defined to be mQTGs. The integrated method enabled us to detect mQTGs with a low false positive rate (=5.71 × 10-4). We also identified 5,654 genes associated with 1,335 specialized metabolites. Of these genes, 4.4% were affected by CNVs, which was more than expected (χ2 test: P < 0.01). This result suggests that CNVs contribute to variations in specialized metabolites within a species. To assess the contribution of CNVs to adaptive evolution in A. thaliana, we examined the selective sweeps around the mQTGs. We observed that the mQTGs with CNVs tended to undergo selective sweeps. These observations imply that variations in specialized metabolites caused by CNVs contribute to the adaptive evolution of A. thaliana.

DNA methylation is an important factor regulating gene expression in organisms. However, whether DNA methylation plays a key role in adaptive evolution is unknown. Here, we show evidence of naturally selected DNA methylation in Arabidopsis thaliana. In comparison with single nucleotide polymorphisms, three types of methylation—methylated CGs (mCGs), mCHGs, and mCHHs—contributed highly to variable gene expression levels among an A. thaliana population. Such variably expressed genes largely affect a large variation of specialized metabolic quantities. Among the three types of methylations, only mCGs located in promoter regions of genes associated with specialized metabolites show a selective sweep signature in the A. thaliana population. Thus, naturally selected mCGs appear to be key mutations that cause the expressional diversity associated with specialized metabolites during plant evolution.

There is considerable diversity in the specialized metabolites within a single plant species (intra-species) and among different plant species (inter-species). The functional divergence associated with gene duplications largely contributes to the inter-species diversity in the specialized metabolites, whereas the intra-species diversity is due to gene dosage changes via gene duplications [i.e., copy number variants (CNVs)] at the intra-species level of evolution. This is because CNVs are thought to undergo associated with less functional divergence at the intra-species level of evolution. However, functional divergence caused by CNVs may induce specialized metabolite diversity at the intra-species and inter-species levels of evolution. We herein discuss the functional divergence of CNVs in metabolic quantitative trait genes (mQTGs). We focused on 5,654 previously identified mQTGs in 270 Arabidopsis thaliana accessions. The ratio of nonsynonymous to synonymous variations tends to be higher for mQTGs with CNVs than for mQTGs without CNVs within A. thaliana accessions, suggesting that CNVs are responsible for the functional divergence among mQTGs at the intra-species level of evolution. To evaluate the contribution of CNVs to inter-species diversity, we calculated the ratio of nonsynonymous to synonymous substitutions in the Arabidopsis lineage. The ratio tends to be higher for the mQTGs with CNVs than for the mQTGs without CNVs. Additionally, we determined that mQTGs with CNVs are subject to positive selection in the Arabidopsis lineage. Our data suggest that CNVs are closely related to functional divergence contributing to adaptations via the production of diverse specialized metabolites at the intra-species and inter-species levels of evolution.