Origin and evolution of genes related to ABA metabolism and its signaling pathways

Since plants cannot move to avoid stress, they have sophisticated acclimation mechanisms against a variety of abiotic stresses. The phytohormone abscisic acid (ABA) plays essential roles in abiotic stress tolerances in land plants. Therefore, it is interesting to address the evolutionary origins of ABA metabolism and its signaling pathways in land plants. Here, we focused on 48 ABA-related Arabidopsis thaliana genes with 11 protein functions, and generated 11 orthologous clusters of ABA-related genes from A. thaliana, Arabidopsis lyrata, Populus trichocarpa, Oryza sativa, Selaginella moellendorffii, and Physcomitrella patens. Phylogenetic analyses suggested that the common ancestor of these six species possessed most of the key protein functions of ABA-related genes. In two species (A. thaliana and O. sativa), duplicate genes related to ABA signaling pathways contribute to the expression variation in different organs or stress responses. In particular, there is significant expansion of gene families related to ABA in evolutionary periods associated with morphological divergence. Taken together, these results suggest that expansion of the gene families related to ABA signaling pathways may have contributed to the sophisticated stress tolerance mechanisms of higher land plants.     

Organization and structural evolution of four multigene families in Arabidopsis thaliana: AtLCAD,AtLGT, AtMYST and AtHD-GL2

The Arabidopsis Genome Initiative has released up to now more than 80% of the genome sequence of Arabidopsis thaliana. About 70% of the identified genes have at least one paralogue. In order to understand the biological function of individual genes, it is essential to study the structure, expression and organization of the entire multigene family. A systematic analysis of multigene families, made possible by the amount of genomic sequence data available, provides important clues for the understanding of genome evolution and plasticity. In this paper, four multigene families of A. thaliana are characterized, namely LCAD, HD-GL2, LGT and MYST. Members of HD-GL2 and LCAD have already been reported in plants. The LGT genes specify proteins containing motifs of glycosyl transferase. No plant genes similar to the LGT genes have been reported to date. The novel MYST family, most likely plant-specific, encodes proteins with no identified function. Sequencing and in silico analysis led to the characterization of 29 novel genes belonging to these four gene families. The organization, structure and evolution of all the members of the four families are discussed, as well as their chromosome location. Expression data of some of the paralogues of each family are also presented.     

Patterns of univariate and multivariate plasticity to elevated carbon dioxide in six European populations of Arabidopsis thaliana

The impact of elevated carbon dioxide on plants is a growing concern in evolutionary ecology and global change biology. Characterizing patterns of phenotypic integration and multivariate plasticity to elevated carbon dioxide can provide insights into ecological and evolutionary dynamics in future human‐altered environments. Here, we examined univariate and multivariate responses to carbon enrichment in six functional traits among six European accessions of Arabidopsis thaliana. We detected phenotypic plasticity in both univariate and multivariate phenotypes, but did not find significant variation in plasticity (genotype by environment interactions) within or among accessions. Eigenvector, eigenvalue variance, and common principal components analyses showed that elevated carbon dioxide altered patterns of trait covariance, reduced the strength of phenotypic integration, and decreased population‐level differentiation in the multivariate phenotype. Our data suggest that future carbon dioxide conditions may influence evolutionary dynamics in natural populations of A. thaliana.     

Acylated monogalactosyl diacylglycerol: prevalence in the plant kingdom and identification of an enzyme catalyzing galactolipid head group acylation in Arabidopsis thaliana

The lipid phase of the thylakoid membrane is mainly composed of the galactolipids mono‐ and digalactosyl diacylglycerol (MGDG and DGDG, respectively). It has been known since the late 1960s that MGDG can be acylated with a third fatty acid to the galactose head group (acyl‐MGDG) in plant leaf homogenates. In certain brassicaceous plants like Arabidopsis thaliana, the acyl‐MGDG frequently incorporates oxidized fatty acids in the form of the jasmonic acid precursor 12‐oxo‐phytodienoic acid (OPDA). In the present study we further investigated the distribution of acylated and OPDA‐containing galactolipids in the plant kingdom. While acyl‐MGDG was found to be ubiquitous in green tissue of plants ranging from non‐vascular plants to angiosperms, OPDA‐containing galactolipids were only present in plants from a few genera. A candidate protein responsible for the acyl transfer was identified in Avena sativa (oat) leaf tissue using biochemical fractionation and proteomics. Knockout of the orthologous gene in A. thaliana resulted in an almost total elimination of the ability to form both non‐oxidized and OPDA‐containing acyl‐MGDG. In addition, heterologous expression of the A. thaliana gene in E. coli demonstrated that the protein catalyzed acylation of MGDG. We thus demonstrate that a phylogenetically conserved enzyme is responsible for the accumulation of acyl‐MGDG in A. thaliana. The activity of this enzyme in vivo is strongly enhanced by freezing damage and the hypersensitive response.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.     

Phylogenomic analysis of gene co‐expression networks reveals the evolution of functional modules

Molecular evolutionary studies correlate genomic and phylogenetic information with the emergence of new traits of organisms. These traits are, however, the consequence of dynamic gene networks composed of functional modules, which might not be captured by genomic analyses. Here, we established a method that combines large‐scale genomic and phylogenetic data with gene co‐expression networks to extensively study the evolutionary make‐up of modules in the moss Physcomitrella patens, and in the angiosperms Arabidopsis thaliana and Oryza sativa (rice). We first show that younger genes are less annotated than older genes. By mapping genomic data onto the co‐expression networks, we found that genes from the same evolutionary period tend to be connected, whereas old and young genes tend to be disconnected. Consequently, the analysis revealed modules that emerged at a specific time in plant evolution. To uncover the evolutionary relationships of the modules that are conserved across the plant kingdom, we added phylogenetic information that revealed duplication and speciation events on the module level. This combined analysis revealed an independent duplication of cell wall modules in bryophytes and angiosperms, suggesting a parallel evolution of cell wall pathways in land plants. We provide an online tool allowing plant researchers to perform these analyses at http://www.gene2function.de .     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

GABI‐DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana

Owing to duplication events in its progenitor, more than 90% of the genes in the Arabidopsis thaliana genome are members of multigene families. A set of 2108 gene families, each consisting of precisely two unlinked paralogous genes, was identified in the nuclear genome of A. thaliana on the basis of sequence similarity. A systematic method for the creation of double knock‐out lines for such gene pairs, designated as DUPLO lines, was established and 200 lines are now publicly available. Their initial phenotypic characterisation led to the identification of seven lines with defects that emerge only in the adult stage. A further six lines display seedling lethality and 23 lines were lethal before germination. Another 14 lines are known to show phenotypes under non‐standard conditions or at the molecular level. Knock‐out of gene pairs with very similar coding sequences or expression profiles is more likely to produce a mutant phenotype than inactivation of gene pairs with dissimilar profiles or sequences. High coding sequence similarity and highly similar expression profiles are only weakly correlated, implying that promoter and coding regions of these gene pairs display different degrees of diversification.     

SMG1 is an ancient nonsense‐mediated mRNA decay effector

Nonsense‐mediated mRNA decay (NMD) is a eukaryotic process that targets selected mRNAs for destruction, for both quality control and gene regulatory purposes. SMG1, the core kinase of the NMD machinery in animals, phosphorylates the highly conserved UPF1 effector protein to activate NMD. However, SMG1 is missing from the genomes of fungi and the model flowering plant Arabidopsis thaliana, leading to the conclusion that SMG1 is animal‐specific and questioning the mechanistic conservation of the pathway. Here we show that SMG1 is not animal‐specific, by identifying SMG1 in a range of eukaryotes, including all examined green plants with the exception of A. thaliana. Knockout of SMG1 by homologous recombination in the basal land plant Physcomitrella patens reveals that SMG1 has a conserved role in the NMD pathway across kingdoms. SMG1 has been lost at various points during the evolution of eukaryotes from multiple lineages, including an early loss in the fungal lineage and a very recent observable gene loss in A. thaliana. These findings suggest that the SMG1 kinase functioned in the NMD pathway of the last common eukaryotic ancestor.     

Predicting genome-wide redundancy using machine learning

Background  

Gene duplication can lead to genetic redundancy, which masks the function of mutated genes in genetic analyses. Methods to increase sensitivity in identifying genetic redundancy can improve the efficiency of reverse genetics and lend insights into the evolutionary outcomes of gene duplication. Machine learning techniques are well suited to classifying gene family members into redundant and non-redundant gene pairs in model species where sufficient genetic and genomic data is available, such as Arabidopsis thaliana, the test case used here.     

Biological implications of the occurrence of 32 members of the XTH (xyloglucan endotransglucosylase/hydrolase) family of proteins in the bryophyte Physcomitrella patens

This comprehensive overview of the xyloglucan endotransglucosylase/hydrolase (XTH) family of genes and proteins in bryophytes, based on research using genomic resources that are newly available for the moss Physcomitrella patens, provides new insights into plant evolution. In angiosperms, the XTH genes are found in large multi‐gene families, probably reflecting the diverse roles of individual XTHs in various cell types. As there are fewer cell types in P. patens than in angiosperms such as Arabidopsis and rice, it is tempting to deduce that there are fewer XTH family genes in bryophytes. However, the present study unexpectedly identified as many as 32 genes that potentially encode XTH family proteins in the genome of P. patens, constituting a fairly large multi‐gene family that is comparable in size with those of Arabidopsis and rice. In situ localization of xyloglucan endotransglucosylase activity in this moss indicates that some P. patens XTH proteins exhibit biochemical functions similar to those found in angiosperms, and that their expression profiles are tissue‐dependent. However, comparison of structural features of families of XTH genes between P. patens and angiosperms demonstrated the existence of several bryophyte‐specific XTH genes with distinct structural and functional features that are not found in angiosperms. These bryophyte‐specific XTH genes might have evolved to meet morphological and functional needs specific to the bryophyte. These findings raise interesting questions about the biological implications of the XTH family of proteins in non‐seed plants.     

Evolutionary rate variation among genes involved in galactomannan biosynthesis in Coffea canephora

The endosperm cell walls of mature coffee seeds accumulate large amounts of mannan storage polysaccharides, which serve as nutrient reserve for embryo and contribute to beverage quality. Our study investigated the evolutionary patterns of key galactomannan (GM) biosynthesis genes using d_N/d_S ratio, synteny, and phylogenetic analysis and detected heterogeneity in rate of evolution among gene copies. Selection ratio index revealed evidence of positive selection in the branch editing gene Coffea canephora alpha (α) galactosidase (Cc‐alpha Gal) at Cc11_g15950 copy (ω = 1.12), whereas strong purifying selection on deleterious mutations was observed in the Coffea canephora uridine diphosphate (UDP)‐glucose 4′‐epimerase (Cc‐UG4E) and Coffea canephora mannose‐1P guanylytransferase (Cc‐MGT) genes controlling the crucial nucleotide carbon sugar building blocks flux in the pathway. Relatively low sequence diversity and strong syntenic linkages were detected in all GM pathway genes except in Cc‐alpha Gal, which suggests a correlation between selection pressure and nucleotide diversity or synteny analysis. In addition, phylogenetic analysis revealed independent evolution or expansion of GM pathway genes in different plant species, with no obvious inferable clustering patterns according to either gene family or congruent with evolutionary plants lineages tested due to high dynamic nature and specific biochemical cell wall modification requirements. Altogether, our study shows a significant high rate of evolutionary variation among GM pathway genes in the diploid C. canephora and demonstrates the inherent variation in evolution of gene copies and their potential role in understanding selection rates in a homogenously connected metabolic pathway.     

Analysis of promoter activity of members of the <Emphasis Type="Italic">PECTATE LYASE-LIKE (PLL)</Emphasis> gene family in cell separation in Arabidopsis

Background  

Pectate lyases depolymerize pectins by catalyzing the eliminative cleavage of α-1,4-linked galacturonic acid. Pectate lyase-like (PLL) genes make up among the largest and most complex families in plants, but their cellular and organismal roles have not been well characterized, and the activity of these genes has been assessed only at the level of entire organs or plant parts, potentially obscuring important sub-organ or cell-type-specific activities. As a first step to understand the potential functional diversity of PLL genes in plants and specificity of individual genes, we utilized a reporter gene approach to document the spatial and temporal promoter activity for 23 of the 26 members of the Arabidopsis thaliana (Arabidopsis) PLL gene family throughout development, focusing on processes involving cell separation.