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.     

Comparison of the starch synthesis genes between maize and rice: copies, chromosome location and expression divergence

Gene duplication and divergence are important evolutionary processes. It has been suggested that a whole genome duplication (WGD) event occurred in the Gramineae, predating its divergence, and a second WGD occurred in maize during its evolution. In this study we compared the fate of the genes involved in the core pathway of starch biosynthesis following the ancient and second WGDs in maize and rice. In total, thirty starch synthesis genes were detected in the maize genome, which covered all the starch synthesis gene families encoded by 27 genes in rice. All of these genes, except ZmGBSSIIb and ZmBEIII, are anchored within large-scale synteny blocks of rice and maize chromosomes. Previous findings and our results indicate that two of the current copies of many starch synthesis genes (including AGPL, AGPS, GBSS, SSII, SSIII, and BEII) probably arose from the ancient WGD in the Gramineae and are still present in the maize and rice genome. Furthermore, two copies of at least six genes (AGPS1, SSIIb, SSIIIb, GBSSII, BEI, and ISA3) appear to have been retained in the maize genome after its second WGD, although complete coding regions were only detected among the duplicate sets of AGPS1, SSIIb, and SSIIIb. The expression patterns of the remaining duplicate sets of starch synthesis genes (AGPL1/2, AGPS1/2, SSIIa/b, SSIIIa/b, GBSSI/II, and BEIIa/b) differ in their expression and could be classified into two groups in maize. The first group is mainly expressed in the endosperm, whereas the second is expressed in other organs and the early endosperm development. The four duplicate sets of ZmGBSSII, ZmSSIIb, ZmSSIIIb and AGPS1, which arose from the second WGD diverged in gene structure and/or expression patterns in maize. These results indicated that some duplicated starch synthesis genes were remained, whereas others diverged in gene structure and/or expression pattern in maize. For most of the duplicated genes, one of the copies has disappeared in the maize genome after the WGD and the subsequent “diploidization”. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Functional Divergence in a Multi-gene Family Is a Key Evolutionary Innovation for Anaerobic Growth in Saccharomyces cerevisiae

The amplification and diversification of genes into large multi-gene families often mark key evolutionary innovations, but this process often creates genetic redundancy that hinders functional investigations. When the model budding yeast Saccharomyces cerevisiae transitions to anaerobic growth conditions, the cell massively induces the expression of seven serine/threonine-rich anaerobically-induced cell wall mannoproteins (anCWMPs): TIP1, TIR1, TIR2, TIR3, TIR4, DAN1, and DAN4. Here, we show that these genes likely derive evolutionarily from a single ancestral anCWMP locus, which was duplicated and translocated to new genomic contexts several times both prior to and following the budding yeast whole genome duplication (WGD) event. Based on synteny and their phylogeny, we separate the anCWMPs into four gene subfamilies. To resolve prior inconclusive genetic investigations of these genes, we constructed a set of combinatorial deletion mutants to determine their contributions toward anaerobic growth in S. cerevisiae. We found that two genes, TIR1 and TIR3, were together necessary and sufficient for the anCWMP contribution to anaerobic growth. Overexpressing either gene alone was insufficient for anaerobic growth, implying that they encode non-overlapping functional roles in the cell during anaerobic growth. We infer from the phylogeny of the anCWMP genes that these two important genes derive from an ancient duplication that predates the WGD event, whereas the TIR1 subfamily experienced gene family amplification after the WGD event. Taken together, the genetic and molecular evidence suggests that one key anCWMP gene duplication event, several auxiliary gene duplication events, and functional divergence underpin the evolution of anaerobic growth in budding yeasts.     

Comparative transcriptional profiling and evolutionary analysis of the GRAM domain family in eukaryotes

The GRAM domain was found in glucosyltransferases, myotubularins and other membrane-associated proteins. So far, functions for majority of these proteins are yet to be uncovered. In order to address the evolutionary and functional significance of this family members, we have performed a comprehensive investigation on their genome-wide identification, phylogenetic relationship and expression divergence in five different organisms representing monocot/dicot plants, vertebrate/invertebrate animals and yeast, namely, Oryza sativa, Arabidopsis thaliana, Mus musculus, Drosophila melanogaster and Saccharomyces cerevisiae, respectively. We have identified 65 members of GRAM domain family from these organisms. Our data revealed that this family was an ancient group and various organisms had evolved into different family sizes. Large-scale genome duplication and divergence in both expression patterns and functions were significantly contributed to the expansion and retention of this family. Mouse and Drosophila members showed higher divergences in their proteins as indicated by higher Ka/Ks ratios and possessed multiple domains in various combinations. However, in plants, their protein functions were possibly retained with a relatively low divergence as signified by lower Ka/Ks ratios and only one additional domain was combined during evolution. On the other hand, this family in all five organisms exhibited high divergence in their expression patterns both at tissue level and under various biotic and abiotic stresses. These highly divergent expression patterns unraveled the complexity of functions of GRAM domain family. Each member may play specialized roles in a specific tissue or stress condition and may function as regulators of environmental and hormonal signaling.     

The Evolutionary History of Sarco(endo)plasmic Calcium ATPase (SERCA)

Investigating the phylogenetic relationships within physiologically essential gene families across a broad range of taxa can reveal the key gene duplication events underlying their family expansion and is thus important to functional genomics studies. P-Type II ATPases represent a large family of ATP powered transporters that move ions across cellular membranes and includes Na⁺/K⁺ transporters, H⁺/K⁺ transporters, and plasma membrane Ca²⁺ pumps. Here, we examine the evolutionary history of one such transporter, the Sarco(endo)plasmic reticulum calcium ATPase (SERCA), which maintains calcium homeostasis in the cell by actively pumping Ca²⁺ into the sarco(endo)plasmic reticulum. Our protein-based phylogenetic analyses across Eukaryotes revealed two monophyletic clades of SERCA proteins, one containing animals, fungi, and plants, and the other consisting of plants and protists. Our analyses suggest that the three known SERCA proteins in vertebrates arose through two major gene duplication events after the divergence from tunicates, but before the separation of fishes and tetrapods. In plants, we recovered two SERCA clades, one being the sister group to Metazoa and the other to Apicomplexa clade, suggesting an ancient duplication in an early eukaryotic ancestor, followed by subsequent loss of one copy in Opisthokonta, the other in protists, and retention of both in plants. We also report relatively recent and independent gene duplication events within invertebrate taxa including tunicates and the leech Helobdella robusta. Thus, it appears that both ancient and recent gene duplication events have played an important role in the evolution of this ubiquitous gene family across the eukaryotic domain.     

Molecular evolution and functional divergence of the cytochrome P450 3 (CYP3) Family in Actinopterygii (ray-finned fish)

Background

The cytochrome P450 (CYP) superfamily is a multifunctional hemethiolate enzyme that is widely distributed from Bacteria to Eukarya. The CYP3 family contains mainly the four subfamilies CYP3A, CYP3B, CYP3C and CYP3D in vertebrates; however, only the Actinopterygii (ray-finned fish) have all four subfamilies and detailed understanding of the evolutionary relationship of Actinopterygii CYP3 family members would be valuable.

Methods and Findings

Phylogenetic relationships were constructed to trace the evolutionary history of the Actinopterygii CYP3 family genes. Selection analysis, relative rate tests and functional divergence analysis were combined to interpret the relationship of the site-specific evolution and functional divergence in the Actinopterygii CYP3 family. The results showed that the four CYP3 subfamilies in Actinopterygii might be formed by gene duplication. The first gene duplication event was responsible for divergence of the CYP3B/C clusters from ancient CYP3 before the origin of the Actinopterygii, which corresponded to the fish-specific whole genome duplication (WGD). Tandem repeat duplication in each of the homologue clusters produced stable CYP3B, CYP3C, CYP3A and CYP3D subfamilies. Acceleration of asymmetric evolutionary rates and purifying selection together were the main force for the production of new subfamilies and functional divergence in the new subset after gene duplication, whereas positive selection was detected only in the retained CYP3A subfamily. Furthermore, nearly half of the functional divergence sites appear to be related to substrate recognition, which suggests that site-specific evolution is closely related with functional divergence in the Actinopterygii CYP3 family.

Conclusions

The split of fish-specific CYP3 subfamilies was related to the fish-specific WGD, and site-specific acceleration of asymmetric evolutionary rates and purifying selection was the main force for the origin of the new subfamilies and functional divergence in the new subset after gene duplication. Site-specific evolution in substrate recognition was related to functional divergence in the Actinopterygii CYP3 family.     

Small auxin upregulated RNA （SAUR） gene family in maize： Identification,evolution, and its phylogenetic comparison with Arabidopsis,rice, and sorghum

Small auxin-up RNAs （.SAURs） are the early auxin- responsive genes represented by a large multigene family in plants. Here, we identified 79 SAUR gene family members from maize （Zea mays subsp, mays） by a reiterative database search and manual annotation. Phylogenetic analysis indicated that the SAUR proteins from Arabidopsis, rice, sorghum, and maize had divided into 16 groups. These genes were non-randomly distributed across the maize chromosomes, and segmental duplication and tandem duplication contributed to the expansion of the maize .SAUR gene family. Synteny analysis established ortholos~J relationships and functional linkages between SAUR genes in maize and sorghum genomes. We also found that the auxin-responsive elements were conserved in the upstream sequences of maize SAUR members. Selection analyses identified some significant site-specific constraints acted on most SAUR paralogs. Expression profiles based on microarray data have provided insights into the possible functional divergence among members of the .SAUR gene family. Quantitative real-time PCR analysis indicated that some of the 10 randomly selected ZmSAUR genes could be induced at least in maize shoot or root tissue tested. The results reveal a comprehensive overview of the maize .SAUR gene family and may pave the way for deciphering their function during pJant development.     

Molecular Properties and Functional Divergence of the Dehydroascorbate Reductase Gene Family in Lower and Higher Plants

Dehydroascorbate reductase (DHAR), which reduces oxidized ascorbate, is important for maintaining an appropriate ascorbate redox state in plant cells. To date, genome-wide molecular characterization of DHARs has only been conducted in bryophytes (Physcomitrella patens) and eudicots (e.g. Arabidopsis thaliana). In this study, to gain a general understanding of the molecular properties and functional divergence of the DHARs in land plants, we further conducted a comprehensive analysis of DHARs from the lycophyte Selaginella moellendorffii, gymnosperm Picea abies and monocot Zea mays. DHARs were present as a small gene family in all of the land plants we examined, with gene numbers ranging from two to four. All the plants contained cytosolic and chloroplastic DHARs, indicating dehydroascorbate (DHA) can be directly reduced in the cytoplasm and chloroplast by DHARs in all the plants. A novel vacuolar DHAR was found in Z. mays, indicating DHA may also be reduced in the vacuole by DHARs in Z. mays. The DHARs within each species showed extensive functional divergence in their gene structures, subcellular localizations, and enzymatic characteristics. This study provides new insights into the molecular characteristics and functional divergence of DHARs in land plants.     

Comparative Genome Analysis of Scutellaria baicalensis and Scutellaria barbata Reveals the Evolution of Active Flavonoid Biosynthesis

Scutellaria baicalensis (S. baicalensis) and Scutellaria barbata (S. barbata) are common medicinal plants of the Lamiaceae family. Both produce specific flavonoid compounds, including baicalein, scutellarein, norwogonin, and wogonin, as well as their glycosides, which exhibit antioxidant and antitumor activities. Here, we report chromosome-level genome assemblies of S. baicalensis and S. barbata with quantitative chromosomal variation (2n = 18 and 2n = 26, respectively). The divergence of S. baicalensis and S. barbata occurred far earlier than previously reported, and a whole-genome duplication (WGD) event was identified. The insertion of long terminal repeat elements after speciation might be responsible for the observed chromosomal expansion and rearrangement. Comparative genome analysis of the congeneric species revealed the species-specific evolution of chrysin and apigenin biosynthetic genes, such as the S. baicalensis-specific tandem duplication of genes encoding phenylalanine ammonia lyase and chalcone synthase, and the S. barbata-specific duplication of genes encoding 4-CoA ligase. In addition, the paralogous duplication, colinearity, and expression diversity of CYP82D subfamily members revealed the functional divergence of genes encoding flavone hydroxylase between S. baicalensis and S. barbata. Analyzing these Scutellaria genomes reveals the common and species-specific evolution of flavone biosynthetic genes. Thus, these findings would facilitate the development of molecular breeding and studies of biosynthesis and regulation of bioactive compounds.     

A Comparative Genome Analysis of PME and PMEI Families Reveals the Evolution of Pectin Metabolism in Plant Cell Walls

Pectins are fundamental polysaccharides in the plant primary cell wall. Pectins are synthesized and secreted to cell walls as highly methyl-esterified polymers and then demethyl-esterified by pectin methylesterases (PMEs), which are spatially regulated by pectin methylesterase inhibitors (PMEIs). Although PME and PMEI genes are pivotal in plant cell wall formation, few studies have focused on the evolutionary patterns of the PME and PMEI gene families. In this study, the gene origin, evolution, and expression diversity of these two families were systematically analyzed using 11 representative species, including algae, bryophytes, lycophytes and flowering land plants. The results show that 1) for the two subfamilies (PME and proPME) of PME, the origin of the PME subfamily is consistent with the appearance of pectins in early charophyte cell walls, 2) Whole genome duplication (WGD) and tandem duplication contribute to the expansion of proPME and PMEI families in land plants, 3) Evidence of selection pressure shows that the proPME and PMEI families have rapidly evolved, particularly the PMEI family in vascular plants, and 4) Comparative expression profile analysis of the two families indicates that the eudicot Arabidopsis and monocot rice have different expression patterns. In addition, the gene structure and sequence analyses show that the origin of the PMEI domain may be derived from the neofunctionalization of the pro domain after WGD. This study will advance the evolutionary understanding of the PME and PMEI families and plant cell wall development.     

Expansion and evolutionary patterns of GDSL-type esterases/lipases in Rosaceae genomes

GDSL-type esterase/lipase (GELP) is mainly characterized by a conserved GDSL domain at N terminus, and is widely found in all living species, both prokaryotes and eukaryotes. GELP gene family consists of a wide range of members playing important roles in plant physiological processes, such as development, stress responses, and functional divergences. In our study, 597 GELP genes were identified from six Rosaceae genomes (i.e., Fragaria vesca, Prunus persica, Prunus avium, Prunus mume, Pyrus bretschneideri, and Malus domestica) by a comprehensive analysis. All GELP genes were further divided into ten subfamilies based on phylogenetic tree analysis. Subfamily D and subfamily E are the two largest subfamilies. Microcollinearity analysis suggested that WGD/segmental events contribute to the expansion of the GELP gene family in M. domestica and P. bretschneideri compared to F. vesca, P. persica, P. avium, and P. mume. Some PbGELPs were expressed during the fruit development of P. bretschneideri and pollen tubes, indicating their activity in these tissues. The expression divergence of PbGELP duplication gene pairs suggests that many mutations were allowed during evolution, although the structure of GELP genes was highly conserved. The current study results provided the feasibility to understand the expansion and evolution patterns of GELP in Rosaceae genomes, and highlight the function during P. bretschneideri fruits and pollen tubes development.     

Comparative genomic analysis of the Sm gene family in rice and maize

Sm proteins are a group of ubiquitous ring-shaped oligomers that function in multiple aspects of RNA metabolism. However, until this study, no comprehensive study incorporating phylogeny, chromosomal location, gene organization, adaptive evolution, expression profiling and functional networks has been reported for rice and maize. In this study, twenty-five and thirty-three Sm genes have been identified in rice and maize, respectively. Phylogenetic analyses identified eighteen gene groups. Results by gene locations indicated that segmental duplication contributes to the expansion of this gene family in rice and maize. Gene organization and motif compositions of the Sm members are highly conserved in each group, indicative of their functional conservation. Expression profiles have provided insights into the possible functional divergence among members of the Sm gene family. Adaptive evolution analyses suggested that purifying selection was the main force driving Sm evolution, but some critical sites might be responsible for functional divergence. In addition, four hundred and seventy-nine interactions were identified by functional network analyses, and most of which were associated with binding, cellular macromolecule biosynthesis, pre-mRNA processing and transferase activity. Overall, the data contribute to a better understanding of the complexity of Sm gene family in rice and maize and will provide a solid foundation for future functional studies.