GC content of plant genes is linked to past gene duplications

The frequency of G and C nucleotides in genomes varies from species to species, and sometimes even between different genes in the same genome. The monocot grasses have a bimodal distribution of genic GC content absent in dicots. We categorized plant genes from 5 dicots and 4 monocot grasses by synteny to related species and determined that syntenic genes have significantly higher GC content than non-syntenic genes at their 5`-end in the third position within codons for all 9 species. Lower GC content is correlated with gene duplication, as lack of synteny to distantly related genomes is associated with past interspersed gene duplications. Two mutation types can account for biased GC content, mutation of methylated C to T and gene conversion from A to G. Gene conversion involves non-reciprocal exchanges between homologous alleles and is not detectable when the alleles are identical or heterozygous for presence-absence variation, both likely situations for genes duplicated to new loci. Gene duplication can cause production of siRNA which can induce targeted methylation, elevating mC→T mutations. Recently duplicated plant genes are more frequently methylated and less likely to undergo gene conversion, each of these factors synergistically creating a mutational environment favoring AT nucleotides. The syntenic genes with high GC content in the grasses compose a subset that have undergone few duplications, or for which duplicate copies were purged by selection. We propose a “biased gene duplication / biased mutation” (BDBM) model that may explain the origin and trajectory of the observed link between duplication and genic GC bias. The BDBM model is supported by empirical data based on joint analyses of 9 angiosperm species with their genes categorized by duplication status, GC content, methylation levels and functional classes.

Genome-wide identification of the Liriodendron chinense WRKY gene family and its diverse roles in response to multiple abiotic stress

Background Liriodendron chinense (Lchi) is a tree species within the Magnoliaceae family and is considered a basal angiosperm. The too low or high temperature or soil drought will restrict its growth as the adverse environmental conditions, thus improving L. chinense abiotic tolerance was the key issues to study. WRKYs are a major family of plant transcription factors known to often be involved in biotic and abiotic stress responses. So far, it is still largely unknown if and how the LchiWRKY gene family is tied to regulating L. chinense stress responses. Therefore, studying the involvement of the WRKY gene family in abiotic stress regulation in L. chinense could be very informative in showing how this tree deals with such stressful conditions. Results In this research, we performed a genome-wide analysis of the Liriodendron chinense (Lchi) WRKY gene family, studying their classification relationships, gene structure, chromosomal locations, gene duplication, cis-element, and response to abiotic stress. The 44 members of the LchiWRKY gene family contain a significant amount of sequence diversity, with their lengths ranging from 525 bp to 40,981 bp. Using classification analysis, we divided the 44 LchiWRKY genes into three phylogenetic groups (I, II, II), with group II then being further divided into five subgroups (IIa, IIb, IIc, IId, IIe). Comparative phylogenetic analysis including the WRKY families from 17 plant species suggested that LchiWRKYs are closely related to the Magnolia Cinnamomum kanehirae WRKY family, and has fewer family members than higher plants. We found the LchiWRKYs to be evenly distributed across 15 chromosomes, with their duplication events suggesting that tandem duplication may have played a major role in LchiWRKY gene expansion model. A Ka/Ks analysis indicated that they mainly underwent purifying selection and distributed in the group IId. Motif analysis showed that LchiWRKYs contained 20 motifs, and different phylogenetic groups contained conserved motif. Gene ontology (GO) analysis showed that LchiWRKYs were mainly enriched in two categories, i.e., biological process and molecular function. Two group IIc members (LchiWRKY10 and LchiWRKY37) contain unique WRKY element sequence variants (WRKYGKK and WRKYGKS). Gene structure analysis showed that most LchiWRKYs possess 3 exons and two different types of introns: the R- and V-type which are both contained within the WRKY domain (WD). Additional promoter cis-element analysis indicated that 12 cis-elements that play different functions in environmental adaptability occur across all LchiWRKY groups. Heat, cold, and drought stress mainly induced the expression of group II and I LchiWRKYs, some of which had undergone gene duplication during evolution, and more than half of which had three exons. LchiWRKY33 mainly responded to cold stress and LchiWRKY25 mainly responded to heat stress, and LchiWRKY18 mainly responded to drought stress, which was almost 4-fold highly expressed, while 5 LchiWRKYs (LchiWRKY5, LchiWRKY23, LchiWRKY14, LchiWRKY27, and LchiWRKY36) responded equally three stresses with more than 6-fold expression. Subcellular localization analysis showed that all LchiWRKYs were localized in the nucleus, and subcellular localization experiments of LchiWRKY18 and 36 also showed that these two transcription factors were expressed in the nucleus. Conclusions This study shows that in Liriodendron chinense, several WRKY genes like LchiWRKY33, LchiWRKY25, and LchiWRKY18, respond to cold or heat or drought stress, suggesting that they may indeed play a role in regulating the tree’s response to such conditions. This information will prove a pivotal role in directing further studies on the function of the LchiWRKY gene family in abiotic stress response and provides a theoretical basis for popularizing afforestation in different regions of China.

Toxoplasma ceramide synthases: a curious case of gene duplication, divergence and key functionality

Toxoplasma gondii is an obligate, intracellular eukaryotic apicomplexan protozoan parasite that can cause foetal damage and abortion in both animals and humans. Sphingolipids have indispensable functions as signaling molecules and are essential and ubiquitous components of eukaryotic membranes that are both synthesized and scavenged by the Apicomplexa. Ceramide is the precursor for all sphingolipids, and here we report the identification, localisation and analyses of the Toxoplasma ceramide synthases Tg CerS1 and Tg CerS2 and, using a conditional gene regulation approach, establish their roles in pathogenicity and parasite fitness. Interestingly, we observed that whilst Tg CerS1 was a fully functional orthologue of the yeast Lag1p capable of catalysing the conversion of sphinganine to ceramide, in contrast Tg CerS2 was catalytically inactive. Furthermore, genomic deletion of Tg CerS1 using CRISPR/Cas-9 led to viable but slow growing parasites indicating its importance but not indispensability. In contrast, genomic knock out of Tg CerS2 was only accessible utilising the rapamycin-inducible Cre recombinase system. Surprisingly, the results demonstrated that this ‘pseudo’ ceramide synthase, Tg CerS2, has an even greater role in parasite fitness than its catalytically active orthologue (Tg CerS1). Phylogenetic analyses indicated that, as in humans and plants, the ceramide synthase isoforms found in Toxoplasma and other Apicomplexa arose through gene duplication. However, in the Apicomplexa the duplicated copy subsequently evolved into a non-functional ‘pseudo’ ceramide synthase. This arrangement is unique to the Apicomplexa and further illustrates the unusual biology that characterize these protozoan parasites, a feature that could potentially be exploited in the development of new antiprotozoals.

HSDFinder: A BLAST-Based Strategy for Identifying Highly Similar Duplicated Genes in Eukaryotic Genomes

Gene duplication is an important evolutionary mechanism capable of providing new genetic material for adaptive and nonadaptive evolution. However, bioinformatics tools for identifying duplicate genes are often limited to the detection of paralogs in multiple species or to specific types of gene duplicates, such as retrocopies. Here, we present a user-friendly, BLAST-based web tool, called HSDFinder, which can identify, annotate, categorize, and visualize highly similar duplicate genes (HSDs) in eukaryotic nuclear genomes. HSDFinder includes an online heatmap plotting option, allowing users to compare HSDs among different species and visualize the results in different Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway functional categories. The external software requirements are BLAST, InterProScan, and KEGG. The utility of HSDFinder was tested on various model eukaryotic species, including Chlamydomonas reinhardtii, Arabidopsis thaliana, Oryza sativa, and Zea mays as well as the psychrophilic green alga Chlamydomonas sp. UWO241, and was proven to be a practical and accurate tool for gene duplication analyses. The web tool is free to use at http://hsdfinder.com. Documentation and tutorials can be found via the GitHub: https://github.com/zx0223winner/HSDFinder.

Evolution of 14-3-3 Proteins in Angiosperm Plants: Recurring Gene Duplication and Loss

14-3-3 proteins are key regulatory factors in plants and are involved in a broad range of physiological processes. We addressed the evolutionary history of 14-3-3s from 46 angiosperm species, including basal angiosperm Amborella and major lineage of monocotyledons and eudicotyledons. Orthologs of Arabidopsis isoforms were detected. There were several rounds of duplication events in the evolutionary history of the 14-3-3 protein family in plants. At least four subfamilies (iota, epsilon, kappa, and psi) formed as a result of ancient duplication in a common ancestor of angiosperm plants. Recent duplication events followed by gene loss in plant lineage, among others Brassicaceae, Fabaceae, and Poaceae, further shaped the high diversity of 14-3-3 isoforms in plants. Coexpression data showed that 14-3-3 proteins formed different functional groups in different species. In some species, evolutionarily related groups of 14-3-3 proteins had coexpressed together under certain physiological conditions, whereas in other species, closely related isoforms expressed in the opposite manner. A possible explanation is that gene duplication and loss is accompanied by functional plasticity of 14-3-3 proteins.

Horizontal Transfer and Evolution of the Biosynthetic Gene Cluster for Benzoxazinoid

Benzoxazinoids are a class of protective and allelopathic plant secondary metabolites, first identified in maize (Zea mays) and subsequently shown to be encoded by a biosynthetic gene cluster (BGC), the Bx cluster. Data mining of mining 40 high-quality grass genomes identified complete Bx clusters (containing genes Bx1 to Bx5 and Bx8) in three genera (Zea, Echinochloa and Dichanthelium) in the Panicoideae and partial clusters in the Triticeae. The Bx cluster originated from gene duplication of native analogues of Bx genes and chromosomal translocation. An ancient Bx cluster including additional Bx genes (e.g., Bx6) is found in ancestral Panicoideae. The ancient Bx cluster was gained by the Triticeae ancestor via a horizontal transfer (HT) event from the ancestral Panicoideae and later separated into three parts on different chromosomes. Bx6 appears to have been under less constrained selection during evolution of the Panicoideae as evidenced by the fact that was translocated ∼1.31-Mb away from the Bx cluster in Z. mays, moved to other chromosomes in Echinochloa, and even lost in Dichanthelium. Further investigation indicated that intense selection and polyploidization shaped the evolutionary trajectory of the Bx cluster in the grass family. This study provides the first case of HT of BGCs among plants and sheds new insights on the evolution of BGCs.SignificanceBiosynthetic gene clustering and horizontal gene transfer are two evolutionary inventions for rapid adaption by organisms. Horizontal transfer of a gene cluster has been reported in fungi and bacteria, but not in plants up to now. By mining the genomes of 40 monocot species, we deciphered the organization of Bx gene cluster, a biosynthetic gene cluster for benzoxazinoids in grasses. We found that the Bx cluster was formed by gene duplication of native analogues of individual Bx genes and directional translocation. More importantly, the Bx cluster in Triticeae was inherited from the Panicoideae via horizontal transfer. Compared with the native analogues, Bx clusters in grasses show constrained purifying selection underscoring their significance in environmental adaption.

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-seq 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.

Muller’s ratchet of the Y chromosome with gene conversion

Muller’s ratchet is a process in which deleterious mutations are fixed irreversibly in the absence of recombination. The degeneration of the Y chromosome, and the gradual loss of its genes, can be explained by Muller’s ratchet. However, most theories consider single-copy genes, and may not be applicable to Y chromosomes, which have a number of duplicated genes in many species, which are probably undergoing concerted evolution by gene conversion. We developed a model of Muller’s ratchet to explore the evolution of the Y chromosome. The model assumes a non-recombining chromosome with both single-copy and duplicated genes. We used analytical and simulation approaches to obtain the rate of gene loss in this model, with special attention to the role of gene conversion. Homogenization by gene conversion makes both duplicated copies either mutated or intact. The former promotes the ratchet, and the latter retards, and we ask which of these counteracting forces dominates under which conditions. We found that the effect of gene conversion is complex, and depends upon the fitness effect of gene duplication. When duplication has no effect on fitness, gene conversion accelerates the ratchet of both single-copy and duplicated genes. If duplication has an additive fitness effect, the ratchet of single-copy genes is accelerated by gene duplication, regardless of the gene conversion rate, whereas gene conversion slows the degeneration of duplicated genes. Our results suggest that the evolution of the Y chromosome involves several parameters, including the fitness effect of gene duplication by increasing dosage and gene conversion rate.