Structural and Functional Evolution of Positively Selected Sites in Pine Glutathione S-Transferase Enzyme Family

Phylogenetic analyses have identified positive selection as an important driver of protein evolution, both structural and functional. However, the lack of appropriate combined functional and structural assays has generally hindered attempts to elucidate patterns of positively selected sites and their effects on enzyme activity and substrate specificity. In this study we investigated the evolutionary divergence of the glutathione S-transferase (GST) family in Pinus tabuliformis, a pine that is widely distributed from northern to central China, including cold temperate and drought-stressed regions. GSTs play important roles in plant stress tolerance and detoxification. We cloned 44 GST genes from P. tabuliformis and found that 26 of the 44 belong to the largest (Tau) class of GSTs and are differentially expressed across tissues and developmental stages. Substitution models identified five positively selected sites in the Tau GSTs. To examine the functional significance of these positively selected sites, we applied protein structural modeling and site-directed mutagenesis. We found that four of the five positively selected sites significantly affect the enzyme activity and specificity; thus their variation broadens the GST family substrate spectrum. In addition, positive selection has mainly acted on secondary substrate binding sites or sites close to (but not directly at) the primary substrate binding site; thus their variation enables the acquisition of new catalytic functions without compromising the protein primary biochemical properties. Our study sheds light on selective aspects of the functional and structural divergence of the GST family in pine and other organisms.     

Increased Expression and Protein Divergence in Duplicate Genes Is Associated with Morphological Diversification

The differentiation of both gene expression and protein function is thought to be important as a mechanism of the functionalization of duplicate genes. However, it has not been addressed whether expression or protein divergence of duplicate genes is greater in those genes that have undergone functionalization compared with those that have not. We examined a total of 492 paralogous gene pairs associated with morphological diversification in a plant model organism (Arabidopsis thaliana). Classifying these paralogous gene pairs into high, low, and no morphological diversification groups, based on knock-out data, we found that the divergence rate of both gene expression and protein sequences were significantly higher in either high or low morphological diversification groups compared with those in the no morphological diversification group. These results strongly suggest that the divergence of both expression and protein sequence are important sources for morphological diversification of duplicate genes. Although both mechanisms are not mutually exclusive, our analysis suggested that changes of expression pattern play the minor role (33%–41%) and that changes of protein sequence play the major role (59%–67%) in morphological diversification. Finally, we examined to what extent duplicate genes are associated with expression or protein divergence exerting morphological diversification at the whole-genome level. Interestingly, duplicate genes randomly chosen from A. thaliana had not experienced expression or protein divergence that resulted in morphological diversification. These results indicate that most duplicate genes have experienced minor functionalization.     

Micro-Plasticity of Genomes As Illustrated by the Evolution of Glutathione Transferases in 12 Drosophila Species

Glutathione transferases (GST) are an ancient superfamily comprising a large number of paralogous proteins in a single organism. This multiplicity of GSTs has allowed the copies to diverge for neofunctionalization with proposed roles ranging from detoxication and oxidative stress response to involvement in signal transduction cascades. We performed a comparative genomic analysis using FlyBase annotations and Drosophila melanogaster GST sequences as templates to further annotate the GST orthologs in the 12 Drosophila sequenced genomes. We found that GST genes in the Drosophila subgenera have undergone repeated local duplications followed by transposition, inversion, and micro-rearrangements of these copies. The colinearity and orientations of the orthologous GST genes appear to be unique in many of the species which suggests that genomic rearrangement events have occurred multiple times during speciation. The high micro-plasticity of the genomes appears to have a functional contribution utilized for evolution of this gene family.     

Functional divergence and catalytic properties of dehydroascorbate reductase family proteins from Populus tomentosa

Dehydroascorbate reductase (DHAR) is a key enzyme in the ascorbate–glutathione cycle that maintains reduced pools of ascorbic acid and serves as an important antioxidant. In this study, to investigate functional divergence of plant DHAR family and catalytic characteristics of the glutathione binding site (G-site) residues of DHAR proteins, we cloned three DHAR genes (PtoDHAR1/2/3) from Populus tomentosa and predicted the G-site residues. PtoDHAR1 protein was localized in chloroplast, while PtoDHAR2/3 proteins showed cytosolic localizations. Three DHAR proteins showed different enzymatic activities, apparent kinetic characteristics, optimum T _m and pH profiles, indicating their functional divergence. Cys20, Lys8, Pro61, Asp72 and Ser73 of PtoDHAR2 were predicted as G-site residues based on their N-terminal amino acid sequence identity and the available crystal structures of glutathione S-transferases. The biochemical functions of these residues are examined in this study through site-directed mutagenesis. The aforementioned five residues are critical components of active sites that contribute to the enzyme’s catalytic activity. Cys20, Pro61 and Asp72 of PtoDHAR2 are also responsible for maintaining proper protein structure. This study provides new insights into the functional divergence of the plant DHAR family and biochemical properties of the G-site residues in plant DHAR proteins.     

Evolution of plant Ash1 SET genes: structural divergence and functional differentiation

Plant Ash1 SET proteins are involved in H3K36 methylation, and play a key role in plant reproductive development. Genes encoding Ash1 SET proteins constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. To investigate the evolutionary history and functional differentiation of the Ash1 SET gene family, we made a comprehensive evolutionary analysis of this gene family from eleven major representatives of green plants. A novel deep sister relationship grouping previously resolved II-1 and II-2 orthologous groups was identified. The absence of AWS domain in the group II-2 suggests that the independent losses of AWS domain have occurred during evolution. A diversity of gene structures in plant Ash1 SET gene family have been presented since the divergence of Physcomitrella patens (moss) from the other land plants. A small proportion of codons in SET domain regions were detected to be under positive selection along the branches ancestral to land plant and angiosperms, which may have allowed changes of substrate specificity among different evolutionary groups while maintaining the primary function of SET domains. Our predictive subcellular localization and comparative anatomical meta-expression analyses can assort with the structural divergences of Ash1 SET proteins.     

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.     

A Comprehensive Analysis of the Cupin Gene Family in Soybean (Glycine max)

Cupin superfamily of proteins, including germin and germin-like proteins (GLPs) from higher plants, is known to play crucial roles in plant development and defense. To date, no systematic analysis has been conducted in soybean (Glycine max) incorporating genome organization, gene structure, expression compendium. In this study, 69 putative Cupin genes were identified from the whole-genome of soybean, which were non-randomly distributed on 17 of the 20 chromosomes. These Gmcupin proteins were phylogenetically clustered into ten distinct subgroups among which the gene structures were highly conserved. Eighteen pairs (52.2%) of duplicate paralogous genes were preferentially retained in duplicated regions of the soybean genome. The distributions of GmCupin genes implied that long segmental duplications contributed significantly to the expansion of the GmCupin gene family. According to the RNA-seq data analysis, most of the Gmcupins were differentially expressed in tissue-specific expression pattern and the expression of some duplicate genes were partially redundant while others showed functional diversity, suggesting the Gmcupins have been retained by substantial subfunctionalization during soybean evolutionary processes. Selective analysis based on single nucleotide polymorphisms (SNPs) in cultivated and wild soybeans revealed sixteen Gmcupins had selected site(s), with all SNPs in Gmcupin10.3 and Gmcupin07.2 genes were selected sites, which implied these genes may have undergone strong selection effects during soybean domestication. Taken together, our results contribute to the functional characterization of Gmcupin genes in soybean.     

Aldehyde Dehydrogenase Gene Superfamily in Populus: Organization and Expression Divergence between Paralogous Gene Pairs

Aldehyde dehydrogenases (ALDHs) constitute a superfamily of NAD(P)⁺-dependent enzymes that catalyze the irreversible oxidation of a wide range of reactive aldehydes to their corresponding nontoxic carboxylic acids. ALDHs have been studied in many organisms from bacteria to mammals; however, no systematic analyses incorporating genome organization, gene structure, expression profiles, and cis-acting elements have been conducted in the model tree species Populus trichocarpa thus far. In this study, a comprehensive analysis of the Populus ALDH gene superfamily was performed. A total of 26 Populus ALDH genes were found to be distributed across 12 chromosomes. Genomic organization analysis indicated that purifying selection may have played a pivotal role in the retention and maintenance of PtALDH gene families. The exon-intron organizations of PtALDHs were highly conserved within the same family, suggesting that the members of the same family also may have conserved functionalities. Microarray data and qRT-PCR analysis indicated that most PtALDHs had distinct tissue-specific expression patterns. The specificity of cis-acting elements in the promoter regions of the PtALDHs and the divergence of expression patterns between nine paralogous PtALDH gene pairs suggested that gene duplications may have freed the duplicate genes from the functional constraints. The expression levels of some ALDHs were up- or down-regulated by various abiotic stresses, implying that the products of these genes may be involved in the adaptation of Populus to abiotic stresses. Overall, the data obtained from our investigation contribute to a better understanding of the complexity of the Populus ALDH gene superfamily and provide insights into the function and evolution of ALDH gene families in vascular plants.     

Relationship between gene duplicability and diversifiability in the topology of biochemical networks

Background

Selective gene duplicability, the extensive expansion of a small number of gene families, is universal. Quantitatively, the number of genes (P_(K)) with K duplicates in a genome decreases precipitously as K increases, and often follows a power law (P_(k)∝k^-α). Functional diversification, either neo- or sub-functionalization, is a major evolution route for duplicate genes.

Results

Using three lines of genomic datasets, we studied the relationship between gene duplicability and diversifiability in the topology of biochemical networks. First, we explored scenario where two pathways in the biochemical networks antagonize each other. Synthetic knockout of respective genes for the two pathways rescues the phenotypic defects of each individual knockout. We identified duplicate gene pairs with sufficient divergences that represent this antagonism relationship in the yeast S. cerevisiae. Such pairs overwhelmingly belong to large gene families, thus tend to have high duplicability. Second, we used distances between proteins of duplicate genes in the protein interaction network as a metric of their diversification. The higher a gene’s duplicate count, the further the proteins of this gene and its duplicates drift away from one another in the networks, which is especially true for genetically antagonizing duplicate genes. Third, we computed a sequence-homology-based clustering coefficient to quantify sequence diversifiability among duplicate genes – the lower the coefficient, the more the sequences have diverged. Duplicate count (K) of a gene is negatively correlated to the clustering coefficient of its duplicates, suggesting that gene duplicability is related to the extent of sequence divergence within the duplicate gene family.

Conclusion

Thus, a positive correlation exists between gene diversifiability and duplicability in the context of biochemical networks – an improvement of our understanding of gene duplicability.     

Identification of a Novel Gig2 Gene Family Specific to Non-Amniote Vertebrates

Gig2 (grass carp reovirus (GCRV)-induced gene 2) is first identified as a novel fish interferon (IFN)-stimulated gene (ISG). Overexpression of a zebrafish Gig2 gene can protect cultured fish cells from virus infection. In the present study, we identify a novel gene family that is comprised of genes homologous to the previously characterized Gig2. EST/GSS search and in silico cloning identify 190 Gig2 homologous genes in 51 vertebrate species ranged from lampreys to amphibians. Further large-scale search of vertebrate and invertebrate genome databases indicate that Gig2 gene family is specific to non-amniotes including lampreys, sharks/rays, ray-finned fishes and amphibians. Phylogenetic analysis and synteny analysis reveal lineage-specific expansion of Gig2 gene family and also provide valuable evidence for the fish-specific genome duplication (FSGD) hypothesis. Although Gig2 family proteins exhibit no significant sequence similarity to any known proteins, a typical Gig2 protein appears to consist of two conserved parts: an N-terminus that bears very low homology to the catalytic domains of poly(ADP-ribose) polymerases (PARPs), and a novel C-terminal domain that is unique to this gene family. Expression profiling of zebrafish Gig2 family genes shows that some duplicate pairs have diverged in function via acquisition of novel spatial and/or temporal expression under stresses. The specificity of this gene family to non-amniotes might contribute to a large extent to distinct physiology in non-amniote vertebrates.     

Genome-Wide Survey and Expression Analysis of Chlamydomonas reinhardtii U-box E3 Ubiquitin Ligases (CrPUBs) Reveal a Functional Lipid Metabolism Module

E3 ubiquitin ligases determine the substrate specificity of ubiquitination. Plant U-box (PUB) E3 ligases, with a typical 70-amino acid U-box domain, participate in plant developmental processes and environmental responses. Thus far, 64 PUB proteins have been identified in Arabidopsis and 77 PUB proteins have been identified in Oryza. However, detailed studies on U-box genes in the model microalgae Chlamydomonas reinhardtii are lacking. Here, we present a comprehensive analysis of the genes encoding U-box family proteins in C. reinhardtii. Following BLASTP analysis, 30 full-length U-box genes were identified in the C. reinhardtii genome sequence. Bioinformatics analyses of CrPUB genes were performed to characterize the phylogenetic relationships, chromosomal locations and gene structures of each member. The 30 identified CrPUB proteins are clustered into 3 distinct subfamilies, and the genes for these proteins are unevenly distributed among 14 chromosomes. Furthermore, the quantitative real-time RT-PCR or semi-quantitative RT-PCR analysis of 30 CrPUB mRNA abundances under nitrogen starvation showed that 18 CrPUB genes were induced by N starvation and that 7 genes were repressed in the N-poor environment. We selected five CrPUB genes exhibiting marked changes in expression under N-free conditions for further analysis in RNAi experiments and examined the oil content of these gene-silenced transgenic strains. The silencing of CrPUB5 and CrPUB14, which are typically down-regulated under N starvation, induced 9.8%-45.0% and 14.4%-61.8% lipid accumulation, respectively. In contrast, the silencing of CrPUB11, CrPUB23 and CrPUB28, which are markedly up-regulated under N-free conditions, decreased the lipid content by 5.5%-27.8%, 8.1%-27.3% and 6.6%-27.9%, respectively. These results provide a useful reference for the identification and functional analysis of this gene family and fundamental information for microalgae lipid metabolism research.     

Consequences of Whole-Genome Triplication as Revealed by Comparative Genomic Analyses of the Wild Radish Raphanus raphanistrum and Three Other Brassicaceae Species

Polyploidization events are frequent among flowering plants, and the duplicate genes produced via such events contribute significantly to plant evolution. We sequenced the genome of wild radish (Raphanus raphanistrum), a Brassicaceae species that experienced a whole-genome triplication event prior to diverging from Brassica rapa. Despite substantial gene gains in these two species compared with Arabidopsis thaliana and Arabidopsis lyrata, ∼70% of the orthologous groups experienced gene losses in R. raphanistrum and B. rapa, with most of the losses occurring prior to their divergence. The retained duplicates show substantial divergence in sequence and expression. Based on comparison of A. thaliana and R. raphanistrum ortholog floral expression levels, retained radish duplicates diverged primarily via maintenance of ancestral expression level in one copy and reduction of expression level in others. In addition, retained duplicates differed significantly from genes that reverted to singleton state in function, sequence composition, expression patterns, network connectivity, and rates of evolution. Using these properties, we established a statistical learning model for predicting whether a duplicate would be retained postpolyploidization. Overall, our study provides new insights into the processes of plant duplicate loss, retention, and functional divergence and highlights the need for further understanding factors controlling duplicate gene fate.     

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.     

Divergence in structure and function of tau class glutathione transferase from Pinus tabulaeformis,P. yunnanensis and P. densata

Glutathione transferases (GSTs) are a family of enzymes that play important roles in stress tolerance and detoxification in plants. The plant GSTs are divided into four classes (phi, tau, zeta and theta), among which tau is the most numerously represented. To date, studies on GSTs in plants have focused largely on crop species. There is extremely little information on the molecular characteristics of GSTs in gymnosperms. Generalization on GST characteristics unique to gymnosperms and the patterns of GST evolution in plants cannot be made before more members of the gene family in conifers are described. In this study we report three new GSTs from Pinus tabulaeformis, Pinus densata and Pinus yunnanensis. Structural and phylogenetic analyses placed these three GSTs in tau class. The tau GST class is subdivided into three clades and this subdivision seems an ancient event that may have pre-dated the gymnosperm and angiosperm split. Sequence analysis revealed a highly conserved N-terminal domain in contrast to a highly variable C-terminal domain. Mutations even outside the critical glutathione-binding site in N-terminal domain can have pronounced effect on GST catalytic property. Thus, sequence similarity does not parallel functional specificity. The high diversity in C-terminal domain determines a wide range of substrate selectivity and specificity among tau GSTs. Thus the a few conserved residues in C-terminal domain seem essential to maintain the structure of the domain and the protein dimer. More extensive data on GST family organization and a thorough gene-by-gene analysis in conifers are needed to advance our understanding of the true diversity and evolution of GST in structure and function in plants.     

hunchback and Ikaros-like zinc finger genes control reproductive system development in Caenorhabditis elegans

Here we provide evidence for a C2H2 zinc finger gene family with similarity to Ikaros and hunchback. The founding member of this family is Caenorhabditis elegans ehn-3, which has important and poorly understood functions in somatic gonad development. We examined the expression and function of four additional hunchback/Ikaros-like (HIL) genes in C. elegans reproductive system development. Two genes, ehn-3 and R08E3.4, are expressed in somatic gonadal precursors (SGPs) and have overlapping functions in their development. In ehn-3; R08E3.4 double mutants, we find defects in the generation of distal tip cells, anchor cells, and spermatheca; three of the five tissues derived from the SGPs. We provide in vivo evidence that C. elegans HIL proteins have functionally distinct zinc finger domains, with specificity residing in the N-terminal set of four zinc fingers and a likely protein-protein interaction domain provided by the C-terminal pair of zinc fingers. In addition, we find that a chimeric human Ikaros protein containing the N-terminal zinc fingers of EHN-3 functions in C. elegans. Together, these results lend support to the idea that the C. elegans HIL genes and Ikaros have similar functional domains. We propose that hunchback, Ikaros, and the HIL genes arose from a common ancestor that was present prior to the divergence of protostomes and deuterostomes.     

Divergence of Gene Body DNA Methylation and Evolution of Plant Duplicate Genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.     

Characterization of a new rhamnogalacturonan acetyl esterase from Bacillus halodurans C-125 with a new putative carbohydrate binding domain

BH1115 is a gene from Bacillus halodurans strain C-125 that hypothetically encodes a rhamnogalacturonan acetyl esterase (RGAE) of the CE-12 family. As confirmation, this gene was cloned, and the product was expressed in Escherichia coli strain Rosetta (DE3) cells and purified. The enzyme obtained was monomeric, with a molecular mass of 45 kDa, and exhibited alkaliphilic properties. A study of the inhibition of the activity by some modulators confirmed that the catalytic triad for the esterase activity was Ser-His-Asp. This enzyme also presents broad substrate specificity and is active toward 7-aminocephalosporanic acid, cephalosporin C, p-nitrophenyl acetate, β-naphthyl acetate, glucose pentaacetate, and acetylated xylan. Moreover, RGAE from B. halodurans achieves a synergistic effect with xylanase A toward acetylated xylan. As a member of the SGNH family, it does not adopt the common α/β hydrolase fold. The homology between the folds of RGAE from Aspergillus aculeatus and the hypothetical YxiM precursor from Bacillus subtilis, which both belong to the SGNH family, illustrates the divergence of such proteins from a common ancestor. Furthermore, the enzyme possesses a putative substrate binding region at the N terminus of the protein which has never been described to date for any RGAE.