Species-specific partial gene duplication in Arabidopsis thaliana evolved novel phenotypic effects on morphological traits under strong positive selection

Gene duplication is increasingly recognized as an important mechanism for the origination of new genes, as revealed by comparative genomic analysis. However, how new duplicate genes contribute to phenotypic evolution remains largely unknown, especially in plants. Here, we identified the new gene EXOV, derived from a partial gene duplication of its parental gene EXOVL in Arabidopsis thaliana. EXOV is a species-specific gene that originated within the last 3.5 million years and shows strong signals of positive selection. Unexpectedly, RNA-sequencing analyses revealed that, despite its young age, EXOV has acquired many novel direct and indirect interactions in which the parental gene does not engage. This observation is consistent with the high, selection-driven substitution rate of its encoded protein, in contrast to the slowly evolving EXOVL, suggesting an important role for EXOV in phenotypic evolution. We observed significant differentiation of morphological changes for all phenotypes assessed in genome-edited and T-DNA insertional single mutants and in double T-DNA insertion mutants in EXOV and EXOVL. We discovered a substantial divergence of phenotypic effects by principal component analyses, suggesting neofunctionalization of the new gene. These results reveal a young gene that plays critical roles in biological processes that underlie morphological evolution in A. thaliana.

A young gene produced by species-specific partial gene duplication quickly evolved critical developmental functions in Arabidopsis.     

Asymmetric Functional Divergence of Young, Dispersed Gene Duplicates in Arabidopsis thaliana

One prediction of the classic Ohno model of gene duplication predicts that new genes form from the asymmetric functional divergence of a newly arisen, redundant duplicate locus. In order to understand the mechanisms which give rise to functional divergence of newly formed dispersed duplicates, we assessed the expression and molecular evolutionary divergence of a suite of 19 highly similar dispersed duplicates in Arabidopsis thaliana. These duplicates have a K _sil equal to or less than 5 % and are specific to the A. thaliana lineage; thus, they predictably represent some of the youngest duplicates in the A. thaliana genome. We found that the majority of young duplicate loci exhibit asymmetric expression patterns, with the daughter locus exhibiting reduced expression across all tissues analyzed relative to the progenitor locus or simply not expressed. Furthermore, daughter loci, on the whole, have significantly more nonsynonymous substitutions than the progenitor loci. We also identified four pairs of loci which exhibit significant (P < 0.05) evolutionary rate asymmetry, three of which exhibit elevated dN/dS in the duplicate copy. We suggest, based on these data, that functional diversification initially takes the form of asymmetric regulatory divergence that can be a direct consequence of the mode of duplication. The reduced and/or absence of expression in the daughter copy relaxes functional constraint on its protein coding sequence leading to the asymmetric accumulation of nonsynonymous mutations. Thus, our data both affirm Ohno’s prediction while explaining the mechanism by which functional divergence initially occurs following duplication for dispersed gene duplicates.     

CCCH-type zinc finger family in maize: genome-wide identification, classification and expression profiling under abscisic acid and drought treatments

Background

CCCH-type zinc finger proteins comprise a large protein family. Increasing evidence suggests that members of this family are RNA-binding proteins with regulatory functions in mRNA processing. Compared with those in animals, functions of CCCH-type zinc finger proteins involved in plant growth and development are poorly understood.

Methodology/Principal Findings

Here, we performed a genome-wide survey of CCCH-type zinc finger genes in maize (Zea mays L.) by describing the gene structure, phylogenetic relationships and chromosomal location of each family member. Promoter sequences and expression profiles of putative stress-responsive members were also investigated. A total of 68 CCCH genes (ZmC3H1-68) were identified in maize and divided into seven groups by phylogenetic analysis. These 68 genes were found to be unevenly distributed on 10 chromosomes with 15 segmental duplication events, suggesting that segmental duplication played a major role in expansion of the maize CCCH family. The Ka/Ks ratios suggested that the duplicated genes of the CCCH family mainly experienced purifying selection with limited functional divergence after duplication events. Twelve maize CCCH genes grouped with other known stress-responsive genes from Arabidopsis were found to contain putative stress-responsive cis-elements in their promoter regions. Seven of these genes chosen for further quantitative real-time PCR analysis showed differential expression patterns among five representative maize tissues and over time in response to abscisic acid and drought treatments.

Conclusions

The results presented in this study provide basic information on maize CCCH proteins and form the foundation for future functional studies of these proteins, especially for those members of which may play important roles in response to abiotic stresses.     

Genome-wide analysis of the maize superoxide dismutase (SOD) gene family reveals important roles in drought and salt responses

Superoxide dismutase proteins (SODs) are antioxidant enzymes with important roles in abiotic stress responses. The SOD gene family has been systematically analyzed in many plants; however, it is still poorly understood in maize. Here, a bioinformatics analysis of maize SOD gene family was conducted by describing gene structure, conserved motifs, phylogenetic relationships, gene duplications, promoter cis-elements and GO annotations. In total, 13 SOD genes were identified in maize and five members were involved in segmental duplication. Phylogenetic analysis indicated that SODs from maize and other plants comprised two groups, which could be further classified into different subgroups, with most members in the same subgroup having the same subcellular localization. The ZmSOD promoters contained 2-10 stress-responsive cis-elements with different distributions. Heatmap analysis indicated that ZmSODs were expressed in most of the detected tissues and organs. The expression patterns of ZmSODs were investigated under drought and salt treatments by qRT-PCR, and most members were responsive to drought or salt stress, especially some ZmSODs with significant expression changes were identified, such as ZmCSD2 and ZmMSD2, suggesting the important roles of ZmSODs in abiotic stress responses. Our results provide an important basis for further functional study of ZmSODs in future study.     

A Brassica rapa Linkage Map of EST-based SNP Markers for Identification of Candidate Genes Controlling Flowering Time and Leaf Morphological Traits

For identification of genes responsible for varietal differences in flowering time and leaf morphological traits, we constructed a linkage map of Brassica rapa DNA markers including 170 EST-based markers, 12 SSR markers, and 59 BAC sequence-based markers, of which 151 are single nucleotide polymorphism (SNP) markers. By BLASTN, 223 markers were shown to have homologous regions in Arabidopsis thaliana, and these homologous loci covered nearly the whole genome of A. thaliana. Synteny analysis between B. rapa and A. thaliana revealed 33 large syntenic regions. Three quantitative trait loci (QTLs) for flowering time were detected. BrFLC1 and BrFLC2 were linked to the QTLs for bolting time, budding time, and flowering time. Three SNPs in the promoter, which may be the cause of low expression of BrFLC2 in the early-flowering parental line, were identified. For leaf lobe depth and leaf hairiness, one major QTL corresponding to a syntenic region containing GIBBERELLIN 20 OXIDASE 3 and one major QTL containing BrGL1, respectively, were detected. Analysis of nucleotide sequences and expression of these genes suggested possible involvement of these genes in leaf morphological traits.     

Isolation and Activity Analysis of a Seed-Abundant soyAP1 Gene Promoter from Soybean

The soybean aspartic proteinase gene soyAP1 has previously been shown to be expressed specifically in soybean seeds. To investigate the expression pattern and active cis-elements of the soyAP1 promoter, the 1,650-bp 5??-upstream genomic DNA fragment named PS-552 was isolated by PCR walking. Sequence analysis revealed that this fragment contains a series of motifs related to seed-specific promoters and some pollen-expressed elements. Stable expression in transgenic Arabidopsis thaliana showed that the PS-552 promoter can regulate beta-glucuronidase gene accumulation in mature seeds at much higher levels than other tissues, especially vegetative tissues, and exhibits similar activity to the 35S promoter in mature seeds. These results show that the PS-552 promoter is a highly active promoter controlling downstream gene expression, mainly in mature seeds. The 5??-end deletion studies of PS-552 showed that the cis-elements of CAAACAC, AACA, E-box, and CCAA play a role in increasing the seed-specific activity. The proportion of mature seed activity and flower activity was increased as the deletion fragment lengthened, indicating that seed cis-elements possibly lessen or suppress the effect of pollen-expressed elements, increasing the activity of PS-552 in mature seeds.     

Late Embryogenesis Abundant (<Emphasis Type="Italic">LEA</Emphasis>) Gene Family in Maize: Identification,Evolution, and Expression Profiles

Late embryogenesis abundant (LEA) proteins are identified as a large and highly diverse group of polypeptides accumulating in response to cellular dehydration in many organisms. However, there are only very limited reports of this protein family in maize until this study. In the present paper, we identified 32 LEA genes in maize. A total of 83 LEA proteins including 51 members in Arabidopsis and 32 putative members in maize were classified into nine groups. Gene organization and motif compositions of the LEA members are highly conserved in each of the groups, indicative of their functional conservation. The predicted ZmLEA genes were non-random distributed across chromosomes, and transposition event and segmental duplication contributed to the expansion of the LEA gene family in maize. Some abiotic stress-responsive cis-elements were also found in the promoters of ZmLEA genes. Microarray expression analyses revealed different accumulation patterns of ZmLEA family members. Moreover, some members of ZmLEAs were regulated under IAA and some abiotic stresses. This study will provide comprehensive information for maize LEA gene family and may pave the way for deciphering their functions in further studies.     

Evolution of Stress-Regulated Gene Expression in Duplicate Genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >~0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments.     

Genetic Control of Individual Differences in Gene-Specific Methylation in Human Brain

We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island “shores” and at further distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites. cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations—the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both cis- and trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.     

Allele-specific DNA methylation analyses associated with siRNAs in Arabidopsis hybrids

Accumulating evidence has suggested that epigenetic marks including DNA methylation,small RNA and histone modification may involve hybrid vigor in plants.However,knowledge about how epigenetic marks in hybrids regulate gene expression is still limited.Based on genome-wide DNA methylation landscapes of Arabidopsis thaliana Ler and C24 ecotypes and their reciprocal F1 hybrids which were obtained in our previous work,we analyzed allele-specific DNA methylation and distinguished cis-and trans-regulated DNA methylation in hybrids.Our study indicated that both cis-and trans-regulated DNA methylation played roles in hybrids,when cis-regulation played a major role in CG methylation and trans-regulation played major roles in CHG and CHH methylation.In addition,we observed correlations between trans-regulated DNA methylation and siRNA densities.Enriched siRNA regions were significantly concurrent with highly trans-regulated DNA methylation regions.Our results illustrated DNA methylation regulation patterns integrated with siRNAs in Arabidopsis hybrids,and shed light on understanding the mechanism of epigenetic reprogramming for hybrid vigor.