Evolution of Gene Duplication in Plants

Assessment of low-coverage nanopore long read sequencing for SNP genotyping in doubled haploid canola (Brassica napus L.)

Despite the high accuracy of short read sequencing (SRS), there are still issues with attaining accurate single nucleotide polymorphism (SNP) genotypes at low sequencing coverage and in highly duplicated genomes due to misalignment. Long read sequencing (LRS) systems, including the Oxford Nanopore Technologies (ONT) minION, have become popular options for de novo genome assembly and structural variant characterisation. The current high error rate often requires substantial post-sequencing correction and would appear to prevent the adoption of this system for SNP genotyping, but nanopore sequencing errors are largely random. Using low coverage ONT minION sequencing for genotyping of pre-validated SNP loci was examined in 9 canola doubled haploids. The minION genotypes were compared to the Illumina sequences to determine the extent and nature of genotype discrepancies between the two systems. The significant increase in read length improved alignment to the genome and the absence of classical SRS biases results in a more even representation of the genome. Sequencing errors are present, primarily in the form of heterozygous genotypes, which can be removed in completely homozygous backgrounds but requires more advanced bioinformatics in heterozygous genomes. Developments in this technology are promising for routine genotyping in the future.

Rapid evolution of α-gliadin gene family revealed by analyzing Gli-2 locus regions of wild emmer wheat

α-Gliadins are a major group of gluten proteins in wheat flour that contribute to the end-use properties for food processing and contain major immunogenic epitopes that can cause serious health-related issues including celiac disease (CD). α-Gliadins are also the youngest group of gluten proteins and are encoded by a large gene family. The majority of the gene family members evolved independently in the A, B, and D genomes of different wheat species after their separation from a common ancestral species. To gain insights into the origin and evolution of these complex genes, the genomic regions of the Gli-2 loci encoding α-gliadins were characterized from the tetraploid wild emmer, a progenitor of hexaploid bread wheat that contributed the AABB genomes. Genomic sequences of Gli-2 locus regions for the wild emmer A and B genomes were first reconstructed using the genome sequence scaffolds along with optical genome maps. A total of 24 and 16 α-gliadin genes were identified for the A and B genome regions, respectively. α-Gliadin pseudogene frequencies of 86% for the A genome and 69% for the B genome were primarily caused by C to T substitutions in the highly abundant glutamine codons, resulting in the generation of premature stop codons. Comparison with the homologous regions from the hexaploid wheat cv. Chinese Spring indicated considerable sequence divergence of the two A genomes at the genomic level. In comparison, conserved regions between the two B genomes were identified that included α-gliadin pseudogenes containing shared nested TE insertions. Analyses of the genomic organization and phylogenetic tree reconstruction indicate that although orthologous gene pairs derived from speciation were present, large portions of α-gliadin genes were likely derived from differential gene duplications or deletions after the separation of the homologous wheat genomes ~ 0.5 MYA. The higher number of full-length intact α-gliadin genes in hexaploid wheat than that in wild emmer suggests that human selection through domestication might have an impact on α-gliadin evolution. Our study provides insights into the rapid and dynamic evolution of genomic regions harboring the α-gliadin genes in wheat.

Tissue-specific study across the stem reveals the chemistry and transcriptome dynamics of birch bark

Tree bark is a highly specialized array of tissues that plays important roles in plant protection and development. Bark tissues develop from two lateral meristems; the phellogen (cork cambium) produces the outermost stem-environment barrier called the periderm, while the vascular cambium contributes with phloem tissues. Although bark is diverse in terms of tissues, functions and species, it remains understudied at higher resolution. We dissected the stem of silver birch (Betula pendula) into eight major tissue types, and characterized these by a combined transcriptomics and metabolomics approach. We further analyzed the varying bark types within the Betulaceae family. The two meristems had a distinct contribution to the stem transcriptomic landscape. Furthermore, inter- and intraspecies analyses illustrated the unique molecular profile of the phellem. We identified multiple tissue-specific metabolic pathways, such as the mevalonate/betulin biosynthesis pathway, that displayed differential evolution within the Betulaceae. A detailed analysis of suberin and betulin biosynthesis pathways identified a set of underlying regulators and highlighted the important role of local, small-scale gene duplication events in the evolution of metabolic pathways. This work reveals the transcriptome and metabolic diversity among bark tissues and provides insights to its development and evolution, as well as its biotechnological applications.

Genetic Interactions and Molecular Evolution of the Duplicated Genes ICARUS2 and ICARUS1 Help Arabidopsis Plants Adapt to Different Ambient Temperatures

Understanding how plants adapt to ambient temperatures has become a major challenge prompted by global climate change. This has led to the identification of several genes regulating the thermal plasticity of plant growth and flowering time. However, the mechanisms accounting for the natural variation and evolution of such developmental plasticity remain mostly unknown. In this study, we determined that natural variation at ICARUS2 (ICA2), which interacts genetically with its homolog ICA1, alters growth and flowering time plasticity in relation to temperature in Arabidopsis (Arabidopsis thaliana). Transgenic analyses demonstrated multiple functional effects for ICA2 and supported the notion that structural polymorphisms in ICA2 likely underlie its natural variation. Two major ICA2 haplogroups carrying distinct functionally active alleles showed high frequency, strong geographic structure, and significant associations with climatic variables related to annual and daily fluctuations in temperature. Genome analyses across the plant phylogeny indicated that the prevalent plant ICA genes encoding two tRNAHis guanylyl transferase 1 units evolved ∼120 million years ago during the early divergence of mono- and dicotyledonous clades. In addition, ICA1/ICA2 duplication occurred specifically in the Camelineae tribe (Brassicaceae). Thus, ICA2 appears to be ubiquitous across plant evolution and likely contributes to climate adaptation through modifications of thermal developmental plasticity in Arabidopsis.

The m6A pathway protects the transcriptome integrity by restricting RNA chimera formation in plants

This study reveals that an m⁶A-assisted polyadenylation pathway comprising conserved m⁶A writer proteins and a plant-specific m⁶A reader contributes to transcriptome integrity in Arabidopsis thaliana by restricting RNA chimera formation at rearranged loci.

Global, segmental, and gene duplication–related processes are driving genome size and complexity in plants. Despite their evolutionary potentials, those processes can also have adverse effects on genome regulation, thus implying the existence of specialized corrective mechanisms. Here, we report that an N6-methyladenosine (m⁶A)–assisted polyadenylation (m-ASP) pathway ensures transcriptome integrity in Arabidopsis thaliana. Efficient m-ASP pathway activity requires the m⁶A methyltransferase-associated factor FIP37 and CPSF30L, an m⁶A reader corresponding to an YT512-B Homology Domain-containing protein (YTHDC)-type domain containing isoform of the 30-kD subunit of cleavage and polyadenylation specificity factor. Targets of the m-ASP pathway are enriched in recently rearranged gene pairs, displayed an atypical chromatin signature, and showed transcriptional readthrough and mRNA chimera formation in FIP37- and CPSF30L-deficient plants. Furthermore, we showed that the m-ASP pathway can also restrict the formation of chimeric gene/transposable-element transcript, suggesting a possible implication of this pathway in the control of transposable elements at specific locus. Taken together, our results point to selective recognition of 3′-UTR m⁶A as a safeguard mechanism ensuring transcriptome integrity at rearranged genomic loci in plants.

A chromosome-level sequence assembly reveals the structure of the Arabidopsis thaliana Nd-1 genome and its gene set

In addition to the BAC-based reference sequence of the accession Columbia-0 from the year 2000, several short read assemblies of THE plant model organism Arabidopsis thaliana were published during the last years. Also, a SMRT-based assembly of Landsberg erecta has been generated that identified translocation and inversion polymorphisms between two genotypes of the species. Here we provide a chromosome-arm level assembly of the A. thaliana accession Niederzenz-1 (AthNd-1_v2c) based on SMRT sequencing data. The best assembly comprises 69 nucleome sequences and displays a contig length of up to 16 Mbp. Compared to an earlier Illumina short read-based NGS assembly (AthNd-1_v1), a 75 fold increase in contiguity was observed for AthNd-1_v2c. To assign contig locations independent from the Col-0 gold standard reference sequence, we used genetic anchoring to generate a de novo assembly. In addition, we assembled the chondrome and plastome sequences. Detailed analyses of AthNd-1_v2c allowed reliable identification of large genomic rearrangements between A. thaliana accessions contributing to differences in the gene sets that distinguish the genotypes. One of the differences detected identified a gene that is lacking from the Col-0 gold standard sequence. This de novo assembly extends the known proportion of the A. thaliana pan-genome.

Old and young duplicate genes reveal different responses to environmental changes in Arachis duranensis

Old and young duplicate genes have been reported in some organisms. However, little is known about the properties of old and young duplicate genes in Arachis. Here, we have identified old and young duplicate genes in Arachis duranensis, and analyzed the evolution, gene complexity, gene expression pattern, and functional divergence between old and young duplicate genes. Our results showed different evolutionary, gene complexity and gene expression patterns, as well as differing correlations between old and young duplicate genes. Gene ontology results showed that old duplicate genes play a crucial role in lipid and amino acid biosynthesis and the oxidation-reduction process and that young duplicate genes are preferentially involved in photosynthesis and response to biotic stimulus. Transcriptome data sets revealed that most old and young duplicate genes had asymmetric function, and only a few duplicate genes exhibited symmetric function under drought and nematode stress. We found that old duplicate genes are preferentially involved in lipid and amino acid metabolism and response to abiotic stress, while young duplicate genes are likely to participate in photosynthesis and response to biotic stress. This work provides a better understanding of the evolution and functional divergence of old and young duplicate genes in A. duranensis.

Genomic dissection and transcriptional profiling of Cysteine-rich receptor-like kinases in five cereals and functional characterization of TaCRK68-A

Cysteine-rich receptor-like kinases (CRK) constitute one of the largest subfamily of receptor-like kinases, which play crucial roles in plant development and stress response. In total, 43, 37, 36, 38 and 170 CRK genes including duplicated genes were identified in the genome of Brachypodium distachyon, Hordeum vulgare, Oryza sativa, Sorghum bicolor and Triticum aestivum, respectively. These CRK proteins were tightly clustered into four phylogenetic groups and exhibited close syntenic relationship among orthologous genes. Majority of CRK proteins contain a transmembrane domain for plasma membrane localization. The organization of exon/intron, domains and motifs were variably conserved. Tissue-specific expression suggested the involvement of certain CRK genes in plant development. Modulated expression revealed their specific stress-responsive functions. Co-expression and interaction analysis indicated their role in signaling. Ks value and divergence time analysis suggested duplication of TaCRK genes before the hybridization of T. aestivum sub-genomes. Expression comparison of duplicated TaCRK genes revealed functional retention, neofunctionalization or pseudo-functionalization. Recombinant expression of a stress-responsive gene TaCRK68-A in Escherichia coli and Saccharomyces cerevisiae displayed enhanced tolerance against heat, drought, cold and salinity stresses. The study suggested vital functions of CRKs during development and stresses, and provides the basis for functional characterization of each gene in future studies.

Evolutionary history and structure of nuclear matrix constituent proteins, the plant analogues of lamins

Nuclear matrix constituent proteins (NMCPs), the structural components of the plant lamina, are considered to be the analogues of lamins in plants based on numerous structural and functional similarities. Current phylogenetic knowledge suggests that, in contrast to lamins, which are widely distributed in eukaryotes, NMCPs are taxonomically restricted to Streptophyta. At present, most information about NMCPs comes from angiosperms, and virtually no data are available from more ancestral groups. In angiosperms, the NMCP family comprises two phylogenetic groups, NMCP1 and NMCP2, which evolved from the NMCP1 and NMCP2 progenitor genes. Based on sequence conservation and the presence of NMCP-specific domains, we determined the structure and number of NMCP genes present in different Streptophyta clades. We analysed 91 species of embryophytes and report additional NMCP sequences from mosses, liverworts, clubmosses, horsetail, ferns, gymnosperms, and Charophyta algae. Our results confirm an origin of NMCPs in Charophyta (the earliest diverging group of Streptophyta), resolve the number and structure of NMCPs in the different clades, and propose the emergence of additional NMCP homologues by whole-genome duplication events. Immunofluorescence microscopy demonstrated localization of a basal NMCP from the moss Physcomitrella patens at the nuclear envelope, suggesting a functional conservation for basal and more evolved NMCPs.

Fine-Tuning the Expression of Duplicate Genes by Translational Regulation in Arabidopsis and Maize

Plant genomes are extensively shaped by various types of gene duplication. However, in this active area of investigation, the vast majority of studies focus on the sequence and transcription of duplicate genes, leaving open the question of how translational regulation impacts the expression and evolution of duplicate genes. We explored this issue by analyzing the ribo- and mRNA-seq data sets across six tissue types and stress conditions in Arabidopsis thaliana and maize (Zea mays). We dissected the relative contributions of transcriptional and translational regulation to the divergence in the abundance of ribosome footprint (RF) for different types of duplicate genes. We found that the divergence in RF abundance was largely programmed at the transcription level and that translational regulation plays more of a modulatory role. Intriguingly, translational regulation is characterized by its strong directionality, with the divergence in translational efficiency (TE) globally counteracting the divergence in mRNA abundance, indicating partial buffering of the transcriptional divergence between paralogs by translational regulation. Divergence in TE was associated with several sequence features. The faster-evolving copy in a duplicate pair was more likely to show lower RF abundance, which possibly results from relaxed purifying selection compared with its paralog. A considerable proportion of duplicates displayed differential TE across tissue types and stress conditions, most of which were enriched in photosynthesis, energy production, and translation-related processes. Additionally, we constructed a database TDPDG-DB (http://www.plantdupribo.tk), providing an online platform for data exploration. Overall, our study illustrates the roles of translational regulation in fine-tuning duplicate gene expression in plants.

Rapid functional divergence after small-scale gene duplication in grasses

BACKGROUND

Gene duplication has played an important role in the evolution and domestication of flowering plants. Yet little is known about how plant duplicate genes evolve and are retained over long timescales, particularly those arising from small-scale duplication (SSD) rather than whole-genome duplication (WGD) events.

RESULTS

We address this question in the Poaceae (grass) family by analyzing gene expression data from nine tissues of Brachypodium distachyon, Oryza sativa japonica (rice), and Sorghum bicolor (sorghum). Consistent with theoretical predictions, expression profiles of most grass genes are conserved after SSD, suggesting that functional conservation is the primary outcome of SSD in grasses. However, we also uncover support for widespread functional divergence, much of which occurs asymmetrically via the process of neofunctionalization. Moreover, neofunctionalization preferentially targets younger (child) duplicate gene copies, is associated with RNA-mediated duplication, and occurs quickly after duplication. Further analysis reveals that functional divergence of SSD-derived genes is positively correlated with both sequence divergence and tissue specificity in all three grass species, and particularly with anther expression in B. distachyon.

CONCLUSIONS

Our results suggest that SSD-derived grass genes often undergo rapid functional divergence that may be driven by natural selection on male-specific phenotypes. These observations are consistent with those in several animal species, suggesting that duplicate genes take similar evolutionary trajectories in plants and animals.

Evolution of Cold Acclimation and Its Role in Niche Transition in the Temperate Grass Subfamily Pooideae

The grass subfamily Pooideae dominates the grass floras in cold temperate regions and has evolved complex physiological adaptations to cope with extreme environmental conditions like frost, winter, and seasonality. One such adaptation is cold acclimation, wherein plants increase their frost tolerance in response to gradually falling temperatures and shorter days in the autumn. However, understanding how complex traits like cold acclimation evolve remains a major challenge in evolutionary biology. Here, we investigated the evolution of cold acclimation in Pooideae and found that a phylogenetically diverse set of Pooideae species displayed cold acclimation capacity. However, comparing differential gene expression after cold treatment in transcriptomes of five phylogenetically diverse species revealed widespread species-specific responses of genes with conserved sequences. Furthermore, we studied the correlation between gene family size and number of cold-responsive genes as well as between selection pressure on coding sequences of genes and their cold responsiveness. We saw evidence of protein-coding and regulatory sequence evolution as well as the origin of novel genes and functions contributing toward evolution of a cold response in Pooideae. Our results reflect that selection pressure resulting from global cooling must have acted on already diverged lineages. Nevertheless, conservation of cold-induced gene expression of certain genes indicates that the Pooideae ancestor may have possessed some molecular machinery to mitigate cold stress. Evolution of adaptations to seasonally cold climates is regarded as particularly difficult. How Pooideae evolved to transition from tropical to temperate biomes sheds light on how complex traits evolve in the light of climate changes.

Ancient Polyploidy and Genome Evolution in Palms

Mechanisms of genome evolution are fundamental to our understanding of adaptation and the generation and maintenance of biodiversity, yet genome dynamics are still poorly characterized in many clades. Strong correlations between variation in genomic attributes and species diversity across the plant tree of life suggest that polyploidy or other mechanisms of genome size change confer selective advantages due to the introduction of genomic novelty. Palms (order Arecales, family Arecaceae) are diverse, widespread, and dominant in tropical ecosystems, yet little is known about genome evolution in this ecologically and economically important clade. Here, we take a phylogenetic comparative approach to investigate palm genome dynamics using genomic and transcriptomic data in combination with a recent, densely sampled, phylogenetic tree. We find conclusive evidence of a paleopolyploid event shared by the ancestor of palms but not with the sister clade, Dasypogonales. We find evidence of incremental chromosome number change in the palms as opposed to one of recurrent polyploidy. We find strong phylogenetic signal in chromosome number, but no signal in genome size, and further no correlation between the two when correcting for phylogenetic relationships. Palms thus add to a growing number of diverse, ecologically successful clades with evidence of whole-genome duplication, sister to a species-poor clade with no evidence of such an event. Disentangling the causes of genome size variation in palms moves us closer to understanding the genomic conditions facilitating adaptive radiation and ecological dominance in an evolutionarily successful, emblematic tropical clade.

Genome-Wide Identification, Comprehensive Gene Feature, Evolution, and Expression Analysis of Plant Metal Tolerance Proteins in Tobacco Under Heavy Metal Toxicity

Plant metal tolerance proteins (MTPs) comprise a family of membrane divalent cation transporters that play essential roles in plant mineral nutrition maintenance and heavy metal stresses resistance. However, the evolutionary relationships and biological functions of MTP family in tobacco remain unclear. In the present study, 26, 13, and 12 MTPs in three main Nicotiana species (N. tabacum, N. sylvestris, and N. tomentosiformis) were identified and designated, respectively. The phylogenetic relationships, gene structures, chromosome distributions, conserved motifs, and domains of NtMTPs were systematic analyzed. According to the phylogenetic features, 26 NtMTPs were classified into three major substrate-specific groups that were Zn-cation diffusion facilitators (CDFs), Zn/Fe-CDFs, and Mn-CDFs, and seven primary groups (1, 5, 6, 7, 8, 9, and 12). All of the NtMTPs contained a modified signature sequence and the cation_efflux domain, whereas some of them also harbored the ZT_dimer. Evolutionary analysis showed that NtMTP family of N. tabacum originated from its parental genome of N. sylvestris and N. tomentosiformis, and further underwent gene loss and expanded via one segmental duplication event. Moreover, the prediction of cis-acting elements (CREs) and the microRNA target sites of NtMTP genes suggested the diverse and complex regulatory mechanisms that control NtMTPs gene expression. Expression profile analysis derived from transcriptome data and quantitative real-time reverse transcription-PCR (qRT-PCR) analysis showed that the tissue expression patterns of NtMTPs in the same group were similar but varied among groups. Besides, under heavy metal toxicity, NtMTP genes exhibited various responses in either tobacco leaves or roots. 19 and 15 NtMTPs were found to response to at least one metal ion treatment in leaves and roots, respectively. In addition, NtMTP8.1, NtMTP8.4, and NtMTP11.1 exhibited Mn transport abilities in yeast cells. These results provided a perspective on the evolution of MTP genes in tobacco and were helpful for further functional characterization of NtMTP genes.

Genome-Wide Identification of the Aux/IAA Family Genes (MdIAA) and Functional Analysis of MdIAA18 for Apple Tree Ideotype

The Aux/IAA (auxin/indole-3-acetic acid) gene family is one of the early auxin-responsive gene families, which play a central role in auxin response. Few reports are involved in Aux/IAA genes in fruit trees, especially in apple (Malus × domestica Borkh.). A total of 33 MdIAA members were identified, of which 27 members contained four conserved domains, whereas the others lost one or two conserved domains. Several cis-elements in promoters of MdIAAs were predicted responsive to hormones and abiotic stress. Tissue-specific expression patterns of MdIAAs in different apple tree ideotypes were investigated by quantitative real-time PCR. A large number of MdIAAs were highly expressed in leaf buds and reproductive organs, and MdIAAs clustered in same group showed similar expression profiles. Overexpression of MdIAA18 in Arabidopsis resulted in compact phenotype. These results indicated that MdIAA genes may be involved in vegetative and reproductive growth of apple. Taken together, the results provide useful clues to reveal the function of MdIAAs in apple and control apple tree architecture by manipulation of MdIAAs.

Evolution of TOR-SnRK dynamics in green plants and its integration with phytohormone signaling networks

The target of rapamycin (TOR)-sucrose non-fermenting 1 (SNF1)-related protein kinase 1 (SnRK1) signaling is an ancient regulatory mechanism that originated in eukaryotes to regulate nutrient-dependent growth. Although the TOR-SnRK1 signaling cascade shows highly conserved functions among eukaryotes, studies in the past two decades have identified many important plant-specific innovations in this pathway. Plants also possess SnRK2 and SnRK3 kinases, which originated from the ancient SnRK1-related kinases and have specialized roles in controlling growth, stress responses and nutrient homeostasis in plants. Recently, an integrative picture has started to emerge in which different SnRKs and TOR kinase are highly interconnected to control nutrient and stress responses of plants. Further, these kinases are intimately involved with phytohormone signaling networks that originated at different stages of plant evolution. In this review, we highlight the evolution and divergence of TOR-SnRK signaling components in plants and their communication with each other as well as phytohormone signaling to fine-tune growth and stress responses in plants.

Comprehensive analysis of Lon proteases in plants highlights independent gene duplication events

The degradation of damaged proteins is essential for cell viability. Lon is a highly conserved ATP-dependent serine-lysine protease that maintains proteostasis. We performed a comparative genome-wide analysis to determine the evolutionary history of Lon proteases. Prokaryotes and unicellular eukaryotes retained a single Lon copy, whereas multicellular eukaryotes acquired a peroxisomal copy, in addition to the mitochondrial gene, to sustain the evolution of higher order organ structures. Land plants developed small Lon gene families. Despite the Lon2 peroxisomal paralog, Lon genes triplicated in the Arabidopsis lineage through sequential evolutionary events including whole-genome and tandem duplications. The retention of Lon1, Lon4, and Lon3 triplicates relied on their differential and even contrasting expression patterns, distinct subcellular targeting mechanisms, and functional divergence. Lon1 seems similar to the pre-duplication ancestral gene unit, whereas the duplication of Lon3 and Lon4 is evolutionarily recent. In the wider context of plant evolution, papaya is the only genome with a single ancestral Lon1-type gene. The evolutionary trend among plants is to acquire Lon copies with ambiguous pre-sequences for dual-targeting to mitochondria and chloroplasts, and a substrate recognition domain that deviates from the ancestral Lon1 type. Lon genes constitute a paradigm of dynamic evolution contributing to understanding the functional fate of gene duplicates.

Adaptive and degenerative evolution of the S-Phase Kinase-Associated Protein 1-Like family in Arabidopsis thaliana

Genome sequencing has uncovered tremendous sequence variation within and between species. In plants, in addition to large variations in genome size, a great deal of sequence polymorphism is also evident in several large multi-gene families, including those involved in the ubiquitin-26S proteasome protein degradation system. However, the biological function of this sequence variation is yet not clear. In this work, we explicitly demonstrated a single origin of retroposed Arabidopsis Skp1-Like (ASK) genes using an improved phylogenetic analysis. Taking advantage of the 1,001 genomes project, we here provide several lines of polymorphism evidence showing both adaptive and degenerative evolutionary processes in ASK genes. Yeast two-hybrid quantitative interaction assays further suggested that recent neutral changes in the ASK2 coding sequence weakened its interactions with some F-box proteins. The trend that highly polymorphic upstream regions of ASK1 yield high levels of expression implied negative expression regulation of ASK1 by an as-yet-unknown transcriptional suppression mechanism, which may contribute to the polymorphic roles of Skp1-CUL1-F-box complexes. Taken together, this study provides new evolutionary evidence to guide future functional genomic studies of SCF-mediated protein ubiquitylation.

Genomic insights into HSFs as candidate genes for high-temperature stress adaptation and gene editing with minimal off-target effects in flax

Flax (Linum usitatissimum) is a cool season crop commercially cultivated for seed oil and stem fibre production. A comprehensive characterization of the heat shock factor (HSF) candidate genes in flax can accelerate genetic improvement and adaptive breeding for high temperature stress tolerance. We report the genome-wide identification of 34 putative HSF genes from the flax genome, which we mapped on 14 of the 15 chromosomes. Through comparative homology analysis, we classified these genes into three broad groups, and sub-groups. The arrangement of HSF-specific protein motifs, DNA-binding domain (DBD) and hydrophobic heptad repeat (HR-A/B), and exon-intron boundaries substantiated the phylogenetic separation of these genes. Orthologous relationships and evolutionary analysis revealed that the co-evolution of the LusHSF genes was due to recent genome duplication events. Digital and RT-qPCR analyses provided significant evidence of the differential expression of the LusHSF genes in various tissues, at various developmental stages, and in response to high-temperature stress. The co-localization of diverse cis-acting elements in the promoters of the LusHSF genes further emphasized their regulatory roles in the abiotic stress response. We further confirmed DNA-binding sites on the LusHSF proteins and designed guide RNA sequences for gene editing with minimal off-target effects. These results will hasten functional investigations of LusHSFs or assist in devising genome engineering strategies to develop high-temperature stress tolerant flax cultivars.

Landscape of gene transposition-duplication within the Brassicaceae family

We developed the CLfinder-OrthNet pipeline that detects co-linearity among multiple closely related genomes, finds orthologous gene groups, and encodes the evolutionary history of each orthologue group into a representative network (OrthNet). Using a search based on network topology, we identified 1,394 OrthNets that included gene transposition-duplication (tr-d) events, out of 17,432 identified in six Brassicaceae genomes. Occurrences of tr-d shared by subsets of Brassicaceae genomes mirrored the divergence times between the genomes and their repeat contents. The majority of tr-d events resulted in truncated open reading frames (ORFs) in the duplicated loci. However, the duplicates with complete ORFs were significantly more frequent than expected from random events. These were derived from older tr-d events and had a higher chance of being expressed. We also found an enrichment of tr-d events with complete loss of intergenic sequence conservation between the original and duplicated loci. Finally, we identified tr-d events uniquely found in two extremophytes among the six Brassicaceae genomes, including tr-d of SALT TOLERANCE 32 and ZINC TRANSPORTER 3 that relate to their adaptive evolution. CLfinder-OrthNet provides a flexible toolkit to compare gene order, visualize evolutionary paths among orthologues as networks, and identify gene loci that share an evolutionary history.

Copy Number Variation in Domestication

Domesticated plants have long served as excellent models for studying evolution. Many genes and mutations underlying important domestication traits have been identified, and most causal mutations appear to be SNPs. Copy number variation (CNV) is an important source of genetic variation that has been largely neglected in studies of domestication. Ongoing work demonstrates the importance of CNVs as a source of genetic variation during domestication, and during the diversification of domesticated taxa. Here, we review how CNVs contribute to evolutionary processes underlying domestication, and review examples of domestication traits caused by CNVs. We draw from examples in plant species, but also highlight cases in animal systems that could illuminate the roles of CNVs in the domestication process.

Variation of gene expression in plants is influenced by gene architecture and structural properties of promoters

In higher eukaryotes, gene architecture and structural properties of promoters have emerged as significant factors influencing variation in number of transcripts (expression level) and specificity of gene expression in a tissue (expression breadth), which eventually shape the phenotype. In this study, transcriptome data of different tissue types at various developmental stages of A. thaliana, O. sativa, S. bicolor and Z. mays have been used to understand the relationship between properties of gene components and its expression. Our findings indicate that in plants, among all gene architecture and structural properties of promoters, compactness of genes in terms of intron content is significantly linked to gene expression level and breadth, whereas in human an exactly opposite scenario is seen. In plants, for the first time we have carried out a quantitative estimation of effect of a particular trait on expression level and breadth, by using multiple regression analysis and it confirms that intron content of primary transcript (as %) is a powerful determinant of expression breadth. Similarly, further regression analysis revealed that among structural properties of the promoters, stability is negatively linked to expression breadth, while DNase1 sensitivity strongly governs gene expression breadth in monocots and gene expression level in dicots. In addition, promoter regions of tissue specific genes are found to be enriched with TATA box and Y-patch motifs. Finally, multi copy orthologous genes in plants are found to be longer, highly regulated and tissue specific.

Genome-wide identification and expression analysis of the AAAP family in Medicago truncatula

The amino acid/auxin permease (AAAP) gene family plays an important role in the long-distance amino acid transport pathway and takes part in various stages of plant growth and development. However, little is known about the AAAP gene family in Medicago truncatula. Here, we identified 86 putative MtAAAP family members using genome sequence information. Based on phylogenetic analysis, these MtAAAP genes were categorized into eight distinct subfamilies. The MtAAAP genes were mapped on 8 chromosomes and duplication events appeared widely, with 19 and 21 pairs of MtAAAP genes showing segment and tandem duplication events, respectively. Ratio of Ka/Ks indicated that duplicated genes underwent purifying selection. Analysis of RNA-seq data showed that MtAAAP genes exhibited specific expression patterns among different tissues and abiotic stress, indicating that MtAAAP members were involved in plant developmental regulation and stress responses. Expression patterns of 16 MtAAAP genes under abiotic stress were verified by qRT-PCR. The present study provides a foundation for the functional analysis of MtAAAPs in developmental regulation and stress responses.

Multispecies genome-wide analysis defines the MAP3K gene family in Gossypium hirsutum and reveals conserved family expansions

Background

Gene families are sets of structurally and evolutionarily related genes – in one or multiple species – that typically share a conserved biological function. As such, the identification and subsequent analyses of entire gene families are widely employed in the fields of evolutionary and functional genomics of both well established and newly sequenced plant genomes. Currently, plant gene families are typically identified using one of two major ways: 1) HMM-profile based searches using models built on Arabidopsis thaliana genes or 2) coding sequence homology searches using curated databases. Integrated databases containing functionally annotated genes and gene families have been developed for model organisms and several important crops; however, a comprehensive methodology for gene family annotation is currently lacking, preventing automated annotation of newly sequenced genomes.

Results

This paper proposes a combined measure of homology identification, motif conservation, phylogenomic and integrated gene expression analyses to define gene family structures in multiple plant species. The MAP3K gene families in seven plant species, including two currently unexamined species Gossypium hirsutum, and Zostera marina, were characterized to reveal new insights into their collective function and evolution and demonstrate the effectiveness of our novel methodology.

Conclusion

Compared with recent reports, this methodology performs significantly better for the identification and analysis of gene family members in several monocots/dicots, diploid as well as polyploid plant species.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2624-9) contains supplementary material, which is available to authorized users.

Expression and regulatory asymmetry of retained Arabidopsis thaliana transcription factor genes derived from whole genome duplication

Background

Transcription factors (TFs) play a key role in regulating plant development and response to environmental stimuli. While most genes revert to single copy after whole genome duplication (WGD) event, transcription factors are retained at a significantly higher rate. Little is known about how TF duplicates have diverged in their expression and regulation, the answer to which may contribute to a better understanding of the elevated retention rate among TFs.

Results

Here we assessed what features may explain differences in the retention of TF duplicates and other genes using Arabidopsis thaliana as a model. We integrated 34 expression, sequence, and conservation features to build a linear model for predicting the extent of duplicate retention following WGD events among TFs and 19 groups of genes with other functions. We found that TFs was the least well predicted, demonstrating the features of TFs are substantially deviated from duplicate genes in other function groups. Consistent with this, the evolution of TF expression patterns and cis-regulatory cites favors the partitioning of ancestral states among the resulting duplicates: one “ancestral” TF duplicate retains most ancestral expression and cis-regulatory sites, while the “non-ancestral” duplicate is enriched for novel regulatory sites. By modeling the retention of ancestral expression and cis-regulatory states in duplicate pairs using a system of differential equations, we found that TF duplicate pairs in a partitioned state are preferentially maintained.

Conclusions

These TF duplicates with asymmetrically partitioned ancestral states are likely maintained because one copy retains ancestral functions while the other, at least in some cases, acquires novel cis-regulatory sites that may be important for novel, adaptive traits.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1398-z) contains supplementary material, which is available to authorized users.

Genome-Wide Identification, Expression Profile, and Alternative Splicing Analysis of the Brassinosteroid-Signaling Kinase (BSK) Family Genes in Arabidopsis

Brassinosteroids (BRs) are steroid hormones essential for different biological processes, ranging from growth to environmental adaptation in plants. The plant brassinosteroid-signaling kinase (BSK) proteins belong to a family of receptor-like cytoplasmic kinases, which have been reported to play an important role in BR signal transduction. However, the knowledge of BSK genes in plants is still quite limited. In the present study, a total of 143 BSK proteins were identified by a genome-wide search in 17 plant species. A phylogenetic analysis showed that the BSK gene originated in embryophytes, with no BSK found in green algae, and these BSK genes were divided into six groups by comparison with orthologs/paralogs. A further study using comparative analyses of gene structure, expression patterns and alternative splicing of BSK genes in Arabidopsis revealed that all BSK proteins shared similar protein structure with some exception and post-translation modifications including sumolyation and ubiquitination. An expression profile analysis showed that most Arabidopsis BSK genes were constitutively expressed in different tissues; of these, several BSK genes were significantly expressed in response to some hormones or abiotic stresses. Furthermore, reverse transcription-polymerase chain reaction (RT-PCR) assays showed that BSK5, BSK7, and BSK9 underwent alternative splicing in specific stress induced and tissue-dependent patterns. Collectively, these results lay the foundation for further functional analyses of these genes in plants.

A Ma10 gene encoding P-type ATPase is involved in fruit organic acid accumulation in apple

Acidity is one of the main determinants of fruit organoleptic quality. Here, comparative transcriptome analysis was conducted between two cultivars that showed a significant difference in fruit acidity, but contained homozygous non-functional alleles at the major gene Ma1 locus controlling apple fruit acidity. A candidate gene for fruit acidity, designated M10, was identified. The M10 gene encodes a P-type proton pump, P3A -ATPase, which facilitates malate uptake into the vacuole. The Ma10 gene is significantly associated with fruit malate content, accounting for ~7.5% of the observed phenotypic variation in apple germplasm. Subcellular localization assay showed that the Ma10 is targeted to the tonoplast. Overexpression of the Ma10 gene can complement the defect in proton transport of the mutant YAK2 yeast strain and enhance the accumulation of malic acid in apple callus. Moreover, its ectopic expression in tomato induces a decrease in fruit pH. These results suggest that the Ma10 gene has the capacity for proton pumping and plays an important role in fruit vacuolar acidification in apple. Our study provides useful knowledge towards comprehensive understanding of the complex mechanism regulating apple fruit acidity.

Molecular evolution of the plant ECERIFERUM1 and ECERIFERUM3 genes involved in aliphatic hydrocarbon production

The Arabidopsis ECERIFERUM1 (CER1) protein is a decarbonylase that converts fatty acid metabolites into alkanes. Alkanes are components of waxes in the plant cuticle, a waterproof barrier serving to protect land plants from both biotic and abiotic stimuli. CER1 enzymes can be used to produce alternative and sustainable hydrocarbons in eukaryotic systems. In this report we identified 193 CER1 and 128 CER3 sequences from 56 land plants respectively. CER1 and CER3 proteins have high amino acid similarity and both are involved in alkane synthesis in Arabidopsis. The common homologues of CER1 and CER3 genes were identified in three species of chlorophytes, which may be one of the earliest plant taxa that possess CER1 and CER3 genes. To facilitate potential applications, the 3-dimensional structure and conserved motifs of CER1 proteins were also characterized. CER1 and CER3 proteins are structurally similar, but CER1 proteins have more conserved histidine-containing motifs common to fatty acid hydroxylases and stearoyl-CoA desaturases. There was no significant loss or gain of protein motifs after ancient and recent duplications, suggesting that varied properties of CER1 proteins may be associated with less-conserved regions. Among 56 land plants, the codon-based assessments of selection modes revealed that neither entire proteins nor individual amino acids of CER1 proteins were significantly subjected to positive selection, indicating that CER1 proteins are highly conserved throughout evolution.

Angiosperm-Wide and Family-Level Analyses of AP2/ERF Genes Reveal Differential Retention and Sequence Divergence After Whole-Genome Duplication

Plants are immobile and often face stressful environmental conditions, prompting the evolution of genes regulating environmental responses. Such evolution is achieved largely through gene duplication and subsequent divergence. One of the most important gene families involved in regulating plant environmental responses and development is the AP2/ERF superfamily; however, the evolutionary history of these genes is unclear across angiosperms and in major angiosperm families adapted to various ecological niches. Specifically, the impact on gene copy number of whole-genome duplication events occurring around the time of the origins of several plant families is unknown. Here, we present the first angiosperm-wide comparative study of AP2/ERF genes, identifying 75 Angiosperm OrthoGroups (AOGs), each derived from an ancestral angiosperm gene copy. Among these AOGs, 21 retain duplicates with increased copy number in many angiosperm lineages, while the remaining 54 AOGs tend to maintain low copy number. Further analyses of multiple species in the Brassicaceae family indicated that family-specific duplicates experienced differential selective pressures in coding regions, with some paralogs showing signs of positive selection. Further, cis regulatory elements also exhibit extensive divergence between duplicates in Arabidopsis. Moreover, comparison of expression levels suggested that AP2/ERF genes with frequently retained duplicates are enriched for broad expression patterns, offering increased opportunities for functional diversification via changes in expression patterns, and providing a mechanism for repeated duplicate retention in some AOGs. Our results represent the most comprehensive evolutionary history of the AP2/ERF gene family, and support the hypothesis that AP2/ERF genes with broader expression patterns are more likely to be retained as duplicates than those with narrower expression profiles, which could lead to a higher chance of duplicate gene subfunctionalization. The greater tendency of some AOGs to retain duplicates, allowing expression and functional divergence, may facilitate the evolution of complex signaling networks in response to new environmental conditions.

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

High-throughput sequencing data and the impact of plant gene annotation quality

The use of draft genomes of different species and re-sequencing of accessions and populations are now common tools for plant biology research. The de novo assembled draft genomes make it possible to identify pivotal divergence points in the plant lineage and provide an opportunity to investigate the genomic basis and timing of biological innovations by inferring orthologs between species. Furthermore, re-sequencing facilitates the mapping and subsequent molecular characterization of causative loci for traits, such as those for plant stress tolerance and development. In both cases high-quality gene annotation-the identification of protein-coding regions, gene promoters, and 5'- and 3'-untranslated regions-is critical for investigation of gene function. Annotations are constantly improving but automated gene annotations still require manual curation and experimental validation. This is particularly important for genes with large introns, genes located in regions rich with transposable elements or repeats, large gene families, and segmentally duplicated genes. In this opinion paper, we highlight the impact of annotation quality on evolutionary analyses, genome-wide association studies, and the identification of orthologous genes in plants. Furthermore, we predict that incorporating accurate information from manual curation into databases will dramatically improve the performance of automated gene predictors.

Myc-like transcriptional factors in wheat: structural and functional organization of the subfamily I members

BACKGROUND

Myc-like regulatory factors carrying the basic helix-loop-helix (bHLH) domain belong to a large superfamily of transcriptional factors (TFs) present in all eukaryotic kingdoms. In plants, the representatives of this superfamily regulate diverse biological processes including growth and development as well as response to various stresses. As members of the regulatory MBW complexes, they participate in biosynthesis of flavonoids. In wheat, only one member (TaMyc1) of the Myc-like TFs family has been studied, while structural and functional organization of further members remained uncharacterized. From two Myc-subfamilies described recently in the genomes of Triticeae tribe species, we investigated thoroughly the members of the subfamily I which includes the TaMyc1 gene.

RESULTS

Comparison of the promoter regions of the Myc subfamily I members in wheat suggested their division into two groups (likely homoeologous sets): TaMyc-1 (TaMyc-A1/TaMyc1, TaMyc-B1, TaMyc-D1) and TaMyc-2 (TaMyc-A2 and TaMyc-D2). It was demonstrated that the TaMyc-D1 copy has lost its functionality due to the frame shift mutation. The study of functional features of the other four copies suggested some of them to be involved in the biosynthesis of anthocyanins. In particular, TaMyc-B1 is assumed to be a co-regulator of the gene TaC1-A1 (encoding R2R3-Myb factor) in the MBW regulatory complex activating anthocyanin synthesis in wheat coleoptile. The mRNA levels of the TaMyc-A1, TaMyc-B1, TaMyc-A2 and TaMyc-D2 genes increased significantly in wheat seedlings exposed to osmotic stress. Salinity stress induced expression of TaMyc-B1 and TaMyc-A2, while TaMyc-A1 was repressed.

CONCLUSIONS

The features of the structural and functional organization of the members of subfamily I of Myc-like TFs in wheat were determined. Myc-like co-regulator (TaMyc-B1) of anthocyanin synthesis in wheat coleoptile was described for the first time. The Myc-encoding genes presumably involved in response to drought and salinity were determined in wheat. The results obtained are important for further manipulations with Myc genes, aimed on increasing wheat adaptability.

Molecular Evolution of the Vacuolar Iron Transporter (VIT) Family Genes in 14 Plant Species

The vacuolar iron transporter (VIT) proteins are involved in the storage and transport of iron. However, the evolution of this gene family in plants is unknown. In this study, I first identified 114 VIT genes in 14 plant species and classified these genes into seven groups by phylogenetic analysis. Conserved gene organization and motif distribution implied conserved function in each group. I also found that tandem duplication, segmental duplication and transposition contributed to the expansion of this gene family. Additionally, several positive selection sites were identified. Divergent expression patterns of soybean VIT genes were further investigated in different development stages and under iron stress. Functional network analysis exhibited 211 physical or functional interactions. The results will provide the basis for further functional studies of the VIT genes in plants.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

Evolution of the modular, disordered stress proteins known as dehydrins

Dehydrins, plant proteins that are upregulated during dehydration stress conditions, have modular sequences that can contain three conserved motifs (the Y-, S-, and K-segments). The presence and order of these motifs are used to classify dehydrins into one of five architectures: Kn, SKn, KnS, YnKn, and YnSKn, where the subscript n describes the number of copies of that motif. In this study, an architectural and phylogenetic analysis was performed on 426 dehydrin sequences that were identified in 53 angiosperm and 3 gymnosperm genomes. It was found that angiosperms contained all five architectures, while gymnosperms only contained Kn and SKn dehydrins. This suggests that the ancestral dehydrin in spermatophytes was either Kn or SKn, and the Y-segment containing dehydrins first arose in angiosperms. A high-level split between the YnSKn dehydrins from either the Kn or SKn dehydrins could not be confidently identified, however, two lower level architectural divisions appear to have occurred after different duplication events. The first likely occurred after a whole genome duplication, resulting in the duplication of a Y3SK2 dehydrin; the duplicate subsequently lost an S- and K- segment to become a Y3K1 dehydrin. The second split occurred after a tandem duplication of a Y1SK2 dehydrin, where the duplicate lost both the Y- and S- segment and gained four K-segments, resulting in a K6 dehydrin. We suggest that the newly arisen Y3K1 dehydrin is possibly on its way to pseudogenization, while the newly arisen K6 dehydrin developed a novel function in cold protection.

Genome-wide identification and comparative evolutionary analysis of the Dof transcription factor family in physic nut and castor bean

DNA-binding with one finger (Dof) proteins comprise a plant-specific transcription factor family involved in plant growth, development and stress responses. This study presents a genome-wide comparison of Dof family genes in physic nut (Jatropha curcas) and castor bean (Ricinus communis), two Euphorbiaceae plants that have not experienced any recent whole-genome duplication. A total of 25 or 24 Dof genes were identified from physic nut and castor genomes, respectively, where JcDof genes are distributed across nine out of 11 chromosomes. Phylogenetic analysis assigned these genes into nine groups representing four subfamilies, and 24 orthologous groups were also proposed based on comparison of physic nut, castor, Arabidopsis and rice Dofs. Conserved microsynteny was observed between physic nut and castor Dof-coding scaffolds, which allowed anchoring of 23 RcDof genes to nine physic nut chromosomes. In contrast to how no recent duplicate was present in castor, two tandem duplications and one gene loss were found in the Dof gene family of physic nut. Global transcriptome profiling revealed diverse patterns of Jc/RcDof genes over various tissues, and key Dof genes involved in flower development and stress response were also identified in physic nut. These findings provide valuable information for further studies of Dof genes in physic nut and castor.

Genome Sequence of Jaltomata Addresses Rapid Reproductive Trait Evolution and Enhances Comparative Genomics in the Hyper-Diverse Solanaceae

Within the economically important plant family Solanaceae, Jaltomata is a rapidly evolving genus that has extensive diversity in flower size and shape, as well as fruit and nectar color, among its ∼80 species. Here, we report the whole-genome sequencing, assembly, and annotation, of one representative species (Jaltomata sinuosa) from this genus. Combining PacBio long reads (25×) and Illumina short reads (148×) achieved an assembly of ∼1.45 Gb, spanning ∼96% of the estimated genome. Ninety-six percent of curated single-copy orthologs in plants were detected in the assembly, supporting a high level of completeness of the genome. Similar to other Solanaceous species, repetitive elements made up a large fraction (∼80%) of the genome, with the most recently active element, Gypsy, expanding across the genome in the last 1–2 Myr. Computational gene prediction, in conjunction with a merged transcriptome data set from 11 tissues, identified 34,725 protein-coding genes. Comparative phylogenetic analyses with six other sequenced Solanaceae species determined that Jaltomata is most likely sister to Solanum, although a large fraction of gene trees supported a conflicting bipartition consistent with substantial introgression between Jaltomata and Capsicum after these species split. We also identified gene family dynamics specific to Jaltomata, including expansion of gene families potentially involved in novel reproductive trait development, and loss of gene families that accompanied the loss of self-incompatibility. This high-quality genome will facilitate studies of phenotypic diversification in this rapidly radiating group and provide a new point of comparison for broader analyses of genomic evolution across the Solanaceae.

Dissecting the Genomic Diversification of Late Embryogenesis Abundant (LEA) Protein Gene Families in Plants

Late embryogenesis abundant (LEA) proteins include eight multigene families that are expressed in response to water loss during seed maturation and in vegetative tissues of desiccation tolerant species. To elucidate LEA proteins evolution and diversification, we performed a comprehensive synteny and phylogenetic analyses of the eight gene families across 60 complete plant genomes. Our integrated comparative genomic approach revealed that synteny conservation and diversification contributed to LEA family expansion and functional diversification in plants. We provide examples that: 1) the genomic diversification of the Dehydrin family contributed to differential evolution of amino acid sequences, protein biochemical properties, and gene expression patterns, and led to the appearance of a novel functional motif in angiosperms; 2) ancient genomic diversification contributed to the evolution of distinct intrinsically disordered regions of LEA_1 proteins; 3) recurrent tandem-duplications contributed to the large expansion of LEA_2; and 4) dynamic synteny diversification played a role on the evolution of LEA_4 and its function on plant desiccation tolerance. Taken together, these results show that multiple evolutionary mechanisms have not only led to genomic diversification but also to structural and functional plasticity among LEA proteins which have jointly contributed to the adaptation of plants to water-limiting environments.

Messages From the Past: New Insights in Plant Lectin Evolution

Lectins are a large and diverse class of proteins, found in all kingdoms of life. Plants are known to express different types of carbohydrate-binding proteins, each containing at least one particular lectin domain which enables them to specifically recognize and bind carbohydrate structures. The group of plant lectins is heterogeneous in terms of structure, biological activity and function. Lectins control various aspects of plant development and defense. Some lectins facilitate recognition of exogenous danger signals or play a role in endogenous signaling pathways, while others are considered as storage proteins or involved in symbiotic relationships. In this study, we revisit the origin of the different plant lectin families in view of the recently reshaped tree of life. Due to new genomic sampling of previously unknown microbial lineages, the tree of life has expanded and was reshaped multiple times. In addition, more plant genomes especially from basal Phragmoplastophyta, bryophytes, and Salviniales (e.g., Chara braunii, Marchantia polymorpha, Physcomitrella patens, Azolla filiculoides, and Salvinia cucullata) have been analyzed, and annotated genome sequences have become accessible. We searched 38 plant genome sequences including core eudicots, monocots, gymnosperms, fern, lycophytes, bryophytes, charophytes, chlorophytes, glaucophytes, and rhodophytes for lectin motifs, performed an extensive comparative analysis of lectin domain architectures, and determined the phylogenetic and evolutionary history of lectins in the plant lineage. In conclusion, we describe the conservation of particular domains in plant lectin sequences obtained from algae to higher plants. The strong conservation of several lectin motifs highlights their significance for plants.

Robust predictions of specialized metabolism genes through machine learning

Specialized metabolites are critical for plant–environment interactions, e.g., attracting pollinators or defending against herbivores, and are important sources of plant-based pharmaceuticals. However, it is unclear what proportion of enzyme-encoding genes play a role in specialized metabolism (SM) as opposed to general metabolism (GM) in any plant species. This is because of the diversity of specialized metabolites and the considerable number of incompletely characterized pathways responsible for their production. In addition, SM gene ancestors frequently played roles in GM. We evaluate features distinguishing SM and GM genes and build a computational model that accurately predicts SM genes. Our predictions provide candidates for experimental studies, and our modeling approach can be applied to other species that produce medicinally or industrially useful compounds.

Plant specialized metabolism (SM) enzymes produce lineage-specific metabolites with important ecological, evolutionary, and biotechnological implications. Using Arabidopsis thaliana as a model, we identified distinguishing characteristics of SM and GM (general metabolism, traditionally referred to as primary metabolism) genes through a detailed study of features including duplication pattern, sequence conservation, transcription, protein domain content, and gene network properties. Analysis of multiple sets of benchmark genes revealed that SM genes tend to be tandemly duplicated, coexpressed with their paralogs, narrowly expressed at lower levels, less conserved, and less well connected in gene networks relative to GM genes. Although the values of each of these features significantly differed between SM and GM genes, any single feature was ineffective at predicting SM from GM genes. Using machine learning methods to integrate all features, a prediction model was established with a true positive rate of 87% and a true negative rate of 71%. In addition, 86% of known SM genes not used to create the machine learning model were predicted. We also demonstrated that the model could be further improved when we distinguished between SM, GM, and junction genes responsible for reactions shared by SM and GM pathways, indicating that topological considerations may further improve the SM prediction model. Application of the prediction model led to the identification of 1,220 A. thaliana genes with previously unknown functions, each assigned a confidence measure called an SM score, providing a global estimate of SM gene content in a plant genome.

Genomic and evolutionary aspects of chloroplast tRNA in monocot plants

BACKGROUND

Chloroplasts are one of the most indispensable organelles that make life forms on the earth possible by their capacity to photosynthesize. These organelles possess a circular genome with a number of coding genes responsible for self-regulation. tRNAs are an important evolutionary-conserved gene family that are responsible for protein translation. However, within the chloroplast genome, tRNA machinery are poorly understood.

RESULTS

In the present study, the chloroplast genome of six monocot plants, Oryza nivara (NC_005973), Oryza sativa (NC_001320), Sachharum officinarum (NC_006084), Sorghum bicolor (NC_008602), Triticum aestivum (NC_002762), and Zea mays (NC_001666) were downloaded and analyzed to identify tRNA sequences. Further analysis of the tRNA sequences in the chloroplast genomes of the monocot plants resulted in the identification of several novel features. The length of tRNAs in the chloroplast genome of the monocot plants ranged from 59 to 155 nucleotides. Pair-wise sequence alignment revealed the presence of a conserved A-C-x-U-A-x-U-A-x-U-x5-U-A-A nucleotide consensus sequence. In addition, the tRNAs in chloroplast genomes of the monocot plants also contain 21-28 anti-codons against 61 sense codons in the genome. They also contain a group I intron and a C-A-U anti-codon for tRNAIle, which is a common anti-codon of tRNAMet. Evolutionary analysis indicates that tRNAs in the chloroplast genome have evolved from multiple common ancestors, and tRNAMet appears to be the ancestral tRNA that underwent duplication and diversification to give rise to other tRNAs.

CONCLUSION

The results obtained from the study of chloroplast tRNA will greatly help to increase our understanding of tRNA biology at a new level. Functional studies of the reported novel aspects of the chloroplast tRNA of the monocot plants will greatly help to decipher their roles in diverse cellular processes.

Transcriptomic signature reveals mechanism of flower bud distortion in witches'-broom disease of soybean (Glycine max)

BACKGROUND

Soybean (Glycine max L. Merril) crop is major source of edible oil and protein for human and animals besides its various industrial uses including biofuels. Phytoplasma induced floral bud distortion syndrome (FBD), also known as witches' broom syndrome (WBS) has been one of the major biotic stresses adversely affecting its productivity. Transcriptomic approach can be used for knowledge discovery of this disease manifestation by morpho-physiological key pathways.

RESULTS

We report transcriptomic study using Illumina HiSeq NGS data of FBD in soybean, revealing 17,454 differentially expressed genes, 5561 transcription factors, 139 pathways and 176,029 genic region putative markers single sequence repeats, single nucleotide polymorphism and Insertion Deletion. Roles of PmbA, Zn-dependent protease, SAP family and auxin responsive system are described revealing mechanism of flower bud distortion having abnormalities in pollen, stigma development. Validation of 10 randomly selected genes was done by qPCR. Our findings describe the basic mechanism of FBD disease, right from sensing of phytoplasma infection by host plant triggering molecular signalling leading to mobilization of carbohydrate and protein, phyllody, abnormal pollen development, improved colonization of insect in host plants to spread the disease. Study reveals how phytoplasma hijacks metabolic machinery of soybean manifesting FBD.

CONCLUSIONS

This is the first report of transcriptomic signature of FBD or WBS disease of soybean revealing morphological and metabolic changes which attracts insect for spread of disease. All the genic region putative markers may be used as genomic resource for variety improvement and new agro-chemical development for disease control to enhance soybean productivity.

Genome-wide identification and expression analyses of WRKY transcription factor family members from chickpea (Cicer arietinum L.) reveal their role in abiotic stress-responses

BACKGROUND

WRKY proteins play a vital role in the regulation of several imperative plant metabolic processes and pathways, especially under biotic and abiotic stresses. Although WRKY genes have been characterized in various major crop plants, their identification and characterization in pulse legumes is still in its infancy. Chickpea (Cicer arietinum L.) is the most important pulse legume grown in arid and semi-arid tropics.

OBJECTIVE

In silico identification and characterization of WRKY transcription factor-encoding genes in chickpea genome.

METHODS

For this purpose, a systematic genome-wide analysis was carried out to identify the non-redundant WRKY transcription factors in the chickpea genome.

RESULTS

We have computationally identified 70 WRKY-encoding non-redundant genes which were randomly distributed on all the chickpea chromosomes except chromosome 8. The evolutionary phylogenetic analysis classified the WRKY proteins into three major groups (I, II and III) and seven sub-groups (IN, IC, IIa, IIb, IIc, IId and IIe). The gene structure analysis revealed the presence of 2-7 introns among the family members. Along with the presence of absolutely conserved signatory WRKY domain, 19 different domains were also found to be conserved in a group-specific manner. Insights of gene duplication analysis revealed the predominant role of segmental duplications for the expansion of WRKY genes in chickpea. Purifying selection seems to be operated during the evolution and expansion of paralogous WRKY genes. The transcriptome data-based in silico expression analysis revealed the differential expression of CarWRKY genes in root and shoot tissues under salt, drought, and cold stress conditions. Moreover, some of these genes showed identical expression pattern under these stresses, revealing the possibility of involvement of these genes in conserved abiotic stress-response pathways.

CONCLUSION

This genome-wide computational analysis will serve as a base to accelerate the functional characterization of WRKY TFs especially under biotic and abiotic stresses.

Effect of water-deficit on tassel development in maize

Maize often exhibits asynchronous pollination under abiotic and biotic stress conditions; however, the molecular basis of this developmental deficiency has not been elucidated. Tassel development is a key process affecting the anthesis-silking interval (ASI) in maize. In this study, we showed that pollen shedding was delayed and ASI was significantly increased in B73 and Chang7-2 inbred lines under water deficit conditions, which resulted in longer barren tip length and decreased yields under both controlled and field conditions. Comparative transcriptome analysis performed on immature tassels derived from plants grown under well-watered and water deficit conditions identified 1931 and 1713 differentially expressed genes (DEGs) in B73 and Chang7-2, respectively. Further, 28 differentially co-expressed transcription factors were identified across both lines. Collectively, we demonstrated that the molecular regulation of tassel development is associated with water deficit stress at early vegetative stage in maize. This finding extends our understanding of the molecular basis of maize tassel development during abiotic stress.

Different divergence events for three pairs of PEBPs in Gossypium as implied by evolutionary analysis

INTRODUCTION

The phosphatidylethanolamine-binding protein (PEBP) gene family plays a crucial role in seed germination, reproductive transformation, and other important developmental processes in plants, but its distribution in Gossypium genomes or species, evolutionary properties, and the fates of multiple duplicated genes remain unclear.

OBJECTIVES

The primary objectives of this study were to elucidate the distribution and characteristics of PEBP genes in Gossypium, as well as the evolutionary pattern of duplication and deletion, and functional differentiation of PEBPs in plants.

METHODS

Using the PEBP protein sequences in Arabidopsis thaliana as queries, blast alignment was carried out for the identification of PEBP genes in four sequenced cotton species. Using the primers designed according to the PEBP genome sequences, PEBP genes were cloned from 15 representative genomes of Gossypium genus, and the gene structure, CDS sequence, protein sequence and properties were predicted and phylogenetic analysis was performed. Taking PEBP proteins of grape as reference, grouping of orthologous gene, analysis of phylogeny and divergence of PEBPs in nine species were conducted to reconstruct the evolutionary pattern of PEBP genes in plants.

RESULTS

We identified and cloned 160 PEBPs from 15 cotton species, and the phylogenetic analysis showed that the genes could be classified into the following three subfamilies: MFT-like, FT-like and TFL1-like. There were eight single orthologous group (OG) members in each diploid and 16 double OG members in each tetraploid. An analysis of the expression and selective pressure indicated that expression divergence and strong purification selection within the same OG presented in the PEBP gene family.

CONCLUSION

An evolutionary pattern of duplication and deletion of the PEBP family in the evolutionary history of Gossypium was suggested, and three pairs of genes resulted from different whole-genome duplication events.

Loss of Phosphoethanolamine N-Methyltransferases Abolishes Phosphatidylcholine Synthesis and Is Lethal

Plants use several pathways to synthesize phosphatidylcholine (PC), the major phospholipid of eukaryotic cells. PC has important structural and signaling roles. One pathway plants use for synthesis is the phospho-base methylation pathway, which forms the head-group phosphocholine through the triple methylation of phosphoethanolamine (PEA) catalyzed by phosphoethanolamine N-methyltransferases (PEAMTs). Our understanding of that pathway and its physiological importance remains limited. We recently reported that disruption of Arabidopsis thaliana PEAMT1/NMT1 and PEAMT3/NMT3 induces severe PC deficiency leading to dwarfism and impaired development. However, the double nmt1 nmt3 knock-out mutant is viable. Here, we show that this is enabled by residual PEAMT activity through a third family member, NMT2. The triple nmt1 nmt2 nmt3 knock-out mutant cannot synthesize PC from PEA and is lethal. This shows that, unlike mammals and yeast, Arabidopsis cannot form PC from phosphatidyl ethanolamine (PE), and demonstrates that methylation of PEA is the sole, and vital, entry point to PC synthesis. We further show that Arabidopsis has evolved an expanded family of four nonredundant PEAMTs through gene duplication and alternate use of the NMT2 promoter. NMT2 encodes two PEAMT variants, which greatly differ in their ability to perform the initial phospho-base methylation of PEA. Five amino acids at the N terminus of PEAMTs are shown to each be critical for the catalysis of that step committing to PC synthesis. As a whole, these findings open new avenues for enzymatic engineering and the exploration of ways to better tune phosphocholine and PC synthesis to environmental conditions for improved plant performance.

Beyond Arabidopsis: Differential UV-B Response Mediated by UVR8 in Diverse Species

Ultraviolet-B radiation (UV-B, 280-315 nm) is an important environmental signal that regulates growth and development in plants. Two dose-dependent UV-B response pathways were described in plants: a specific one, mediated by UVR8 (the specific UV-B receptor) and an unspecific one, activated by the oxidative damage produced by radiation. The constitutively expressed receptor appears inactive as a dimer, with the two monomers dissociating upon UV-B irradiation. The monomer then interacts with COP1, an ubiquitin ligase, hindering its ability to poly-ubiquitinate transcriptional factor HY5, thus averting its degradation and activating the photomorphogenic response. HY5 induces the synthesis of proteins RUP1 and RUP2, which interact with UVR8, releasing COP1, and inducing the re-dimerization of UVR8. This mechanism has been thoroughly characterized in Arabidopsis, where studies have demonstrated that the UVR8 receptor is key in UV-B response. Although Arabidopsis importance as a model plant many mechanisms described in this specie differ in other plants. In this paper, we review the latest information regarding UV-B response mediated by UVR8 in different species, focusing on the differences reported compared to Arabidopsis. For instance, UVR8 is not only induced by UV-B but also by other agents that are expressed differentially in diverse tissues. Also, in some of the species analyzed, proteins with low homology to RUP1 and RUP2 were detected. We also discuss how UVR8 is involved in other developmental and stress processes unrelated to UV-B. We conclude that the receptor is highly versatile, showing differences among species.

Comparative Genomics Analysis Provides New Insight Into Molecular Basis of Stomatal Movement in Kalanchoë fedtschenkoi

CO₂ uptake and water loss in plants are regulated by microscopic pores on the surface of leaves, called stomata. This enablement of gas exchange by the opening and closing of stomata is one of the most essential processes in plant photosynthesis and transpiration, affecting water-use efficiency (WUE) and thus drought susceptibility. In plant species with crassulacean acid metabolism (CAM) photosynthesis, diel stomatal movement pattern is inverted relative to C₃ and C₄ photosynthesis species, resulting in much higher WUE and drought tolerance. However, little is known about the molecular basis of stomatal movement in CAM species. The goal of this study is to identify candidate genes that could play a role in stomatal movement in an obligate CAM species, Kalanchoë fedtschenkoi. By way of a text-mining approach, proteins were identified in various plant species, spanning C₃, C₄, and CAM photosynthetic types, which are orthologous to proteins known to be involved in stomatal movement. A comparative analysis of diel time-course gene expression data was performed between K. fedtschenkoi and two C₃ species (i.e., Arabidopsis thaliana and Solanum lycopersicum) to identify differential gene expression between the dusk and dawn phases of the 24-h cycle. A rescheduled catalase gene known to be involved in stomatal movement was identified, suggesting a role for H₂O₂ in CAM-like stomatal movement. Overall, these results provide new insights into the molecular regulation of stomatal movement in CAM plants, facilitating genetic improvement of drought resistance in agricultural crops through manipulation of stomata-related genes.