Many or most genes in Arabidopsis transposed after the origin of the order Brassicales

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Detection of colinear blocks and synteny and evolutionary analyses based on utilization of MCScanX

As different taxa evolve, gene order often changes slowly enough that chromosomal 'blocks' with conserved gene orders (synteny) are discernible. The MCScanX toolkit ( https://github.com/wyp1125/MCScanX ) was published in 2012 as freely available software for the detection of such 'colinear blocks' and subsequent synteny and evolutionary analyses based on genome-wide gene location and protein sequence information. Owing to its simplicity and high efficiency for colinear block detection, MCScanX provides a powerful tool for conducting diverse synteny and evolutionary analyses. Moreover, the detection of colinear blocks has been embraced as an integral step for pangenome graph construction. Here, new application trends of MCScanX are explored, striving to better connect this increasingly used tool to other tools and accelerate insight generation from exponentially growing sequence data. We provide a detailed protocol that covers how to install MCScanX on diverse platforms, tune parameters, prepare input files from data from the National Center for Biotechnology Information, run MCScanX and its visualization and evolutionary analysis tools, and connect MCScanX with external tools, including MCScanX-transposed, Circos and SynVisio. This protocol is easily implemented by users with minimal computational background and is adaptable to new data of interest to them. The data and utility programs for this protocol can be obtained from http://bdx-consulting.com/mcscanx-protocol .

Genome-Wide Analysis of MYB Gene Family in Chrysanthemum ×morifolium Provides Insights into Flower Color Regulation

MYBs constitute the second largest transcription factor (TF) superfamily in flowering plants with substantial structural and functional diversity, which have been brought into focus because they affect flower colors by regulating anthocyanin biosynthesis. Up to now, the genomic data of several Chrysanthemum species have been released, which provides us with abundant genomic resources for revealing the evolution of the MYB gene family in Chrysanthemum species. In the present study, comparative analyses of the MYB gene family in six representative species, including C. lavandulifolium, C. seticuspe, C. ×morifolium, Helianthus annuus, Lactuca sativa, and Arabidopsis thaliana, were performed. A total of 1104 MYBs, which were classified into four subfamilies and 35 lineages, were identified in the three Chrysanthemum species (C. lavandulifolium, C. seticuspe, and C. ×morifolium). We found that whole-genome duplication and tandem duplication are the main duplication mechanisms that drove the occurrence of duplicates in CmMYBs (particularly in the R2R3-MYB subfamily) during the evolution of the cultivated chrysanthemums. Sequence structure and selective pressure analyses of the MYB gene family revealed that some of R2R3-MYBs were subjected to positive selection, which are mostly located on the distal telomere segments of the chromosomes and contain motifs 7 and 8. In addition, the gene expression analysis of CmMYBs in different organs and at various capitulum developmental stages of C. ×morifolium indicated that CmMYBS2, CmMYB96, and CmMYB109 might be the negative regulators for anthocyanin biosynthesis. Our results provide the phylogenetic context for research on the genetic and functional evolution of the MYB gene family in Chrysanthemum species and deepen our understanding of the regulatory mechanism of MYB TFs on the flower color of C. ×morifolium.

Comparative Analysis Reveals Different Evolutionary Fates and Biological Functions in Wheat Duplicated Genes (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is a staple food crop that provides 20% of total human calorie consumption. Gene duplication has been considered to play an important role in evolution by providing new genetic resources. However, the evolutionary fates and biological functions of the duplicated genes in wheat remain to be elucidated. In this study, the resulting data showed that the duplicated genes evolved faster with shorter gene lengths, higher codon usage bias, lower expression levels, and higher tissue specificity when compared to non-duplicated genes. Our analysis further revealed functions of duplicated genes in various biological processes with significant enrichment to environmental stresses. In addition, duplicated genes derived from dispersed, proximal, tandem, transposed, and whole-genome duplication differed in abundance, evolutionary rate, gene compactness, expression pattern, and genetic diversity. Tandem and proximal duplicates experienced stronger selective pressure and showed a more compact gene structure with diverse expression profiles than other duplication modes. Moreover, genes derived from different duplication modes showed an asymmetrical evolutionary pattern for wheat A, B, and D subgenomes. Several candidate duplication hotspots associated with wheat domestication or polyploidization were characterized as potential targets for wheat molecular breeding. Our comprehensive analysis revealed the evolutionary trajectory of duplicated genes and laid the foundation for future functional studies on wheat.

A high-quality Bougainvillea genome provides new insights into evolutionary history and pigment biosynthetic pathways in the Caryophyllales

Bougainvillea is a perennial ornamental shrub that is highly regarded in ornamental horticulture around the world. However, the absence of genome data limits our understanding of the pathways involved in bract coloration and breeding. Here, we report a chromosome-level assembly of the giga-genome of Bougainvillea × buttiana 'Mrs Butt', a cultivar thought to be the origin of many other Bougainvillea cultivars. The assembled genome is ~5 Gb with a scaffold N50 of 151 756 278 bp and contains 86 572 genes which have undergone recent whole-genome duplication. We confirmed that multiple rounds of whole-genome multiplication have occurred in the evolutionary history of the Caryophyllales, reconstructed the relationship in the Caryophyllales at whole genome level, and found discordance between species and gene trees as the result of complex introgression events. We investigated betalain and anthocyanin biosynthetic pathways and found instances of independent evolutionary innovations in the nine different Caryophyllales species. To explore the potential formation mechanism of diverse bract colors in Bougainvillea, we analyzed the genes involved in betalain and anthocyanin biosynthesis and found extremely low expression of ANS and DFR genes in all cultivars, which may limit anthocyanin biosynthesis. Our findings indicate that the expression pattern of the betalain biosynthetic pathway did not directly correlate with bract color, and a higher expression level in the betalain biosynthetic pathway is required for colored bracts. This improved understanding of the correlation between gene expression and bract color allows plant breeding outcomes to be predicted with greater certainty.

Genome-wide analysis of the serine carboxypeptidase-like protein family reveals Ga09G1039 is involved in fiber elongation in cotton

The Gossypium is a model genus for understanding polyploidy and the evolutionary pattern of inheritance. This study aimed to investigate the characteristics of SCPLs in different cotton species and their role in fiber development. A total of 891 genes from one typical monocot and ten dicot species were naturally divided into three classes based on phylogenetic analysis. The SCPL gene family in cotton has undergone intense purifying selection with some functional variation. Segmental duplication and whole genome duplication were shown to be the two main reasons for the increase in the number of genes during cotton evolution. The identification of Gh_SCPL genes exhibiting differential expression in particular tissues or response to environmental stimuli provides a new measure for the in-depth characterization of selected genes of importance. Ga09G1039 was involved in the developmental process of fibers and ovules, and it is significantly different from proteins from other cotton species in terms of phylogenetic, gene structure, conserved protein motifs and tertiary structure. Overexpression of Ga09G1039 significantly increased the length of stem trichomes. Ga09G1039 may be a serine carboxypeptidase protein with hydrolase activity, according to functional region, prokaryotic expression, and western blotting analysis. The results provide a comprehensive overview of the genetic basis of SCPLs in Gossypium and further our knowledge in understanding the key aspects of SCPLs in cotton with their potential role in fiber development and stress resistance.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Genomic imbalance modulates transposable element expression in maize

Genomic imbalance refers to the more severe phenotypic consequences of changing part of a chromosome compared with the whole genome set. Previous genome imbalance studies in maize have identified prevalent inverse modulation of genes on the unvaried chromosomes (trans) with both the addition or subtraction of chromosome arms. Transposable elements (TEs) comprise a substantial fraction of the genome, and their reaction to genomic imbalance is therefore of interest. Here, we analyzed TE expression using RNA-seq data of aneuploidy and ploidy series and found that most aneuploidies showed an inverse modulation of TEs, but reductions in monosomy and increases in disomy and trisomy were also common. By contrast, the ploidy series showed little TE modulation. The modulation of TEs and genes in the same experimental group were compared, and TEs showed greater modulation than genes, especially in disomy. Class I and II TEs were differentially modulated in most aneuploidies, and some superfamilies in each TE class also showed differential modulation. Finally, the significantly upregulated TEs in three disomies (TB-7Lb, TB9Lc, and TB-10L19) did not increase the proportion of adjacent gene expression when compared with non-differentially expressed TEs, indicating that modulations of TEs do not compound the effect on genes. These results suggest that the prevalent inverse TE modulation in aneuploidy results from stoichiometric upset of the regulatory machinery used by TEs, similar to the response of core genes to genomic imbalance.

MicroRNAs play regulatory roles in genomic balance

Classic genetics studies found that genomic imbalance caused by changing the dosage of part of the genome (aneuploidy) has more detrimental effects than altering the dosage of the whole genome (ploidy). Previous analysis revealed global modulation of gene expression triggered by aneuploidy across various species, including maize (Zea mays), Arabidopsis, yeast, mammals, etc. Plant microRNAs (miRNAs) are a class of 20- to 24-nt endogenous small noncoding RNAs that carry out post-transcriptional gene expression regulation. That miRNAs and their putative targets are preferentially retained as duplicates after whole-genome duplication, as are many transcription factors and signaling components, indicates miRNAs are likely to be dosage-sensitive and potentially involved in genomic balance networks. This review addresses the following questions regarding the role of miRNAs in genomic imbalance. (1) How do aneuploidy and polyploidy impact the expression of miRNAs? (2) Do miRNAs play a regulatory role in modulating the expression of their targets under genomic imbalance?

The characteristics of mRNA m6A methylomes in allopolyploid Brassica napus and its diploid progenitors

Genome duplication events, comprising whole-genome duplication and single-gene duplication, produce a complex genomic context leading to multiple levels of genetic changes. However, the characteristics of m⁶A modification, the most widespread internal eukaryotic mRNA modification, in polyploid species are still poorly understood. This study revealed the characteristics of m⁶A methylomes within the early formation and following the evolution of allopolyploid Brassica napus. We found a complex relationship between m⁶A modification abundance and gene expression level depending on the degree of enrichment or presence/absence of m⁶A modification. Overall, the m⁶A genes had lower gene expression levels than the non-m⁶A genes. Allopolyploidization may change the expression divergence of duplicated gene pairs with identical m⁶A patterns and diverged m⁶A patterns. Compared with duplicated genes, singletons with a higher evolutionary rate exhibited higher m⁶A modification. Five kinds of duplicated genes exhibited distinct distributions of m⁶A modifications in transcripts and gene expression level. In particular, tandem duplication-derived genes showed unique m⁶A modification enrichment around the transcript start site. Active histone modifications (H3K27ac and H3K4me3) but not DNA methylation were enriched around genes of m⁶A peaks. These findings provide a new understanding of the features of m ⁶A modification and gene expression regulation in allopolyploid plants with sophisticated genomic architecture.

Excision and reinsertion of Ac macrotransposons in maize

Eukaryotic Macrotransposons (MTns) can be formed by 2 nearby elements flanking a segment of host DNA. The maize Ac transposon can form Ac::MTns, but little is known about Ac::MTn transposition activities. Here, we studied 3 Ac::MTns at the maize p1 locus, each of which is composed of a segment of maize p1 genomic DNA (up to 15 kb) bounded by a fractured Ac element (fAc, 2039 bp), and a full-length Ac element in direct orientation. The resulting Ac::MTns are of 16, 16.5, and 22 kb total length. From these 3 Ac::MTns, we identified 10 independent cases of macrotransposition, and observed similar features of transposition between Ac::MTn and standard Ac/Ds, including characteristic excision footprints and insertion target site duplications. Nine out of the 10 Ac::MTn reinsertion targets were genetically linked to the donor sites, another similarity with Ac/Ds standard transposition. We also identified a MTn-like structure in the maize B73 reference genome and 5 NAM founder lines. The MTn in diverse lines is flanked by target site duplications, confirming the historic occurrence of MTn transposition during genome evolution. Our results show that Ac::MTns are capable of mobilizing segments of DNA long enough to include a typical full-length plant gene and in theory could erode gene colinearity in syntenic regions during plant genome evolution.

Transcription factor retention through multiple polyploidization steps in wheat

Whole-genome duplication is widespread in plant evolutionary history and is followed by nonrandom gene loss to return to a diploid state. Across multiple angiosperm species, the retained genes tend to be dosage-sensitive regulatory genes such as transcription factors, yet data for younger polyploid species is sparse. Here, we analyzed the retention, expression, and genetic variation in transcription factors in the recent allohexaploid bread wheat (Triticum aestivum L.). By comparing diploid, tetraploid, and hexaploid wheat, we found that, following each of two hybridization and whole-genome duplication events, the proportion of transcription factors in the genome increased. Transcription factors were preferentially retained over other genes as homoeologous groups in tetraploid and hexaploid wheat. Across cultivars, transcription factor homoeologs contained fewer deleterious missense mutations than nontranscription factors, suggesting that transcription factors are maintained as three functional homoeologs in hexaploid wheat populations. Transcription factor homoeologs were more strongly coexpressed than nontranscription factors, indicating conservation of function between homoeologs. We found that the B3, MADS-M-type, and NAC transcription factor families were less likely to have three homoeologs present than other families, which was associated with low expression levels and high levels of tandem duplication. Together, our results show that transcription factors are preferentially retained in polyploid wheat genomes although there is variation between families. Knocking out one transcription factor homoeolog to alter gene dosage, using TILLING or CRISPR, could generate new phenotypes for wheat breeding.

The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

Interplay between gene expression and gene architecture as a consequence of gene and genome duplications: evidence from metabolic genes of Arabidopsis thaliana

Gene and genome duplications have been widespread during the evolution of flowering plant which resulted in the increment of biological complexity as well as creation of plasticity of a genome helping the species to adapt to changing environments. Duplicated genes with higher evolutionary rates can act as a mechanism of generating novel functions in secondary metabolism. In this study, we explored duplication as a potential factor governing the expression heterogeneity and gene architecture of Primary Metabolic Genes (PMGs) and Secondary Metabolic Genes (SMGs) of Arabidopsis thaliana. It is remarkable that different types of duplication processes controlled gene expression and tissue specificity differently in PMGs and SMGs. A complex relationship exists between gene architecture and expression patterns of primary and secondary metabolic genes. Our study reflects, expression heterogeneity and gene structure variation of primary and secondary metabolism in Arabidopsis thaliana are partly results of duplication events of different origins. Our study suggests that duplication has differential effect on PMGs and SMGs regarding expression pattern by controlling gene structure, epigenetic modifications, multifunctionality and subcellular compartmentalization. This study provides an insight into the evolution of metabolism in plants in the light of gene and genome scale duplication.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01188-2.

Gene duplication and rate variation in the evolution of plastid ACCase and Clp genes in angiosperms

While the chloroplast (plastid) is known for its role in photosynthesis, it is also involved in many other metabolic pathways essential for plant survival. As such, plastids contain an extensive suite of enzymes required for non-photosynthetic processes. The evolution of the associated genes has been especially dynamic in flowering plants (angiosperms), including examples of gene duplication and extensive rate variation. We examined the role of ongoing gene duplication in two key plastid enzymes, the acetyl-CoA carboxylase (ACCase) and the caseinolytic protease (Clp), responsible for fatty acid biosynthesis and protein turnover, respectively. In plants, there are two ACCase complexes-a homomeric version present in the cytosol and a heteromeric version present in the plastid. Duplications of the nuclear-encoded homomeric ACCase gene and retargeting of one resultant protein to the plastid have been previously reported in multiple species. We find that these retargeted homomeric ACCase proteins exhibit elevated rates of sequence evolution, consistent with neofunctionalization and/or relaxation of selection. The plastid Clp complex catalytic core is composed of nine paralogous proteins that arose via ancient gene duplication in the cyanobacterial/plastid lineage. We show that further gene duplication occurred more recently in the nuclear-encoded core subunits of this complex, yielding additional paralogs in many species of angiosperms. Moreover, in six of eight cases, subunits that have undergone recent duplication display increased rates of sequence evolution relative to those that have remained single copy. We also compared substitution patterns between pairs of Clp core paralogs to gain insight into post-duplication evolutionary routes. These results show that gene duplication and rate variation continue to shape the plastid proteome.

One Hundred Years of Gene Balance: How Stoichiometric Issues Affect Gene Expression, Genome Evolution, and Quantitative Traits

A century ago experiments with the flowering plant Datura stramonium and the fruit fly Drosophila melanogaster revealed that adding an extra chromosome to a karyotype was much more detrimental than adding a whole set of chromosomes. This phenomenon was referred to as gene balance and has been recapitulated across eukaryotic species. Here, we retrace some developments in this field. Molecular studies suggest that the basis of balance involves stoichiometric relationships of multi-component interactions. This concept has implication for the mechanisms controlling gene expression, genome evolution, sex chromosome evolution/dosage compensation, speciation mechanisms, and the underlying genetics of quantitative traits.

Recurrent mutations promote widespread structural and functional divergence of MULE-derived genes in plants

Transposable element (TE)-derived genes are increasingly recognized as major sources conferring essential traits in agriculturally important crops but underlying evolutionary mechanisms remain obscure. We updated previous annotations and constructed 18,744 FAR-RED IMPAIRED RESPONSE1 (FAR1) genes, a transcription factor family derived from Mutator-like elements (MULEs), from 80 plant species, including 15,546 genes omitted in previous annotations. In-depth sequence comparison of the updated gene repertoire revealed that FAR1 genes underwent continuous structural divergence via frameshift and nonsense mutations that caused premature translation termination or specific domain truncations. CRISPR/Cas9-based genome editing and transcriptome analysis determined a novel gene involved in fertility-regulating transcription of rice pollen, denoting the functional capacity of our re-annotated gene models especially in monocots which had the highest copy numbers. Genomic evidence showed that the functional gene adapted by obtaining a shortened form through a frameshift mutation caused by a tandem duplication of a 79-bp sequence resulting in premature translation termination. Our findings provide improved resources for comprehensive studies of FAR1 genes with beneficial agricultural traits and unveil novel evolutionary mechanisms generating structural divergence and subsequent adaptation of TE-derived genes in plants.

Genomic imbalance determines positive and negative modulation of gene expression in diploid maize

Genomic imbalance caused by changing the dosage of individual chromosomes (aneuploidy) has a more detrimental effect than varying the dosage of complete sets of chromosomes (ploidy). We examined the impact of both increased and decreased dosage of 15 distal and 1 interstitial chromosomal regions via RNA-seq of maize (Zea mays) mature leaf tissue to reveal new aspects of genomic imbalance. The results indicate that significant changes in gene expression in aneuploids occur both on the varied chromosome (cis) and the remainder of the genome (trans), with a wider spread of modulation compared with the whole-ploidy series of haploid to tetraploid. In general, cis genes in aneuploids range from a gene-dosage effect to dosage compensation, whereas for trans genes the most common effect is an inverse correlation in that expression is modulated toward the opposite direction of the varied chromosomal dosage, although positive modulations also occur. Furthermore, this analysis revealed the existence of increased and decreased effects in which the expression of many genes under genome imbalance are modulated toward the same direction regardless of increased or decreased chromosomal dosage, which is predicted from kinetic considerations of multicomponent molecular interactions. The findings provide novel insights into understanding mechanistic aspects of gene regulation.

Predominantly inverse modulation of gene expression in genomically unbalanced disomic haploid maize

The phenotypic consequences of the addition or subtraction of part of a chromosome is more severe than changing the dosage of the whole genome. By crossing diploid trisomies to a haploid inducer, we identified 17 distal segmental haploid disomies that cover ∼80% of the maize genome. Disomic haploids provide a level of genomic imbalance that is not ordinarily achievable in multicellular eukaryotes, allowing the impact to be stronger and more easily studied. Transcriptome size estimates revealed that a few disomies inversely modulate most of the transcriptome. Based on RNA sequencing, the expression levels of genes located on the varied chromosome arms (cis) in disomies ranged from being proportional to chromosomal dosage (dosage effect) to showing dosage compensation with no expression change with dosage. For genes not located on the varied chromosome arm (trans), an obvious trans-acting effect can be observed, with the majority showing a decreased modulation (inverse effect). The extent of dosage compensation of varied cis genes correlates with the extent of trans inverse effects across the 17 genomic regions studied. The results also have implications for the role of stoichiometry in gene expression, the control of quantitative traits, and the evolution of dosage-sensitive genes.

Starch Synthesis-Related Genes (SSRG) Evolution in the Genus Oryza

Cooking quality is an important attribute in Common/Asian rice (Oryzasativa L.) varieties, being highly dependent on grain starch composition. This composition is known to be highly dependent on a cultivar’s genetics, but the way in which their genes express different phenotypes is not well understood. Further analysis of variation of grain quality genes using new information obtained from the wild relatives of rice should provide important insights into the evolution and potential use of these genetic resources. All analyses were conducted using bioinformatics approaches. The analysis of the protein sequences of grain quality genes across the Oryza suggest that the deletion/mutation of amino acids in active sites result in variations that can negatively affect specific steps of starch biosynthesis in the endosperm. On the other hand, the complete deletion of some genes in the wild species may not affect the amylose content. Here we present new insights for Starch Synthesis-Related Genes (SSRGs) evolution from starch-specific rice phenotypes.

Asymmetric epigenome maps of subgenomes reveal imbalanced transcription and distinct evolutionary trends in Brassica napus

The complexity of the epigenome landscape and transcriptional regulation is significantly increased during plant polyploidization, which drives genome evolution and contributes to the increased adaptability to diverse environments. However, a comprehensive epigenome map of Brassica napus is still unavailable. In this study, we performed integrative analysis of five histone modifications, RNA polymerase II occupancy, DNA methylation, and transcriptomes in two B. napus lines (2063A and B409), and established global maps of regulatory elements, chromatin states, and their dynamics for the whole genome (including the An and Cn subgenomes) in four tissue types (young leaf, flower bud, silique, and root) of these two lines. Approximately 65.8% of the genome was annotated with different epigenomic signals. Compared with the Cn subgenome, the An subgenome possesses a higher level of active epigenetic marks and lower level of repressive epigenetic marks. Genes from subgenome-unique regions contribute to the major differences between the An and Cn subgenomes. Asymmetric histone modifications between homeologous gene pairs reflect their biased expression patterns. We identified a novel bivalent chromatin state (with H3K4me1 and H3K27me3) in B. napus that is associated with tissue-specific gene expression. Furthermore, we observed that different types of duplicated genes have discrepant patterns of histone modification and DNA methylation levels. Collectively, our findings provide a valuable epigenetic resource for allopolyploid plants.

Evolution of Conserved Noncoding Sequences in Arabidopsis thaliana

Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence-absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.

A large genome with chromosome-scale assembly sheds light on the evolutionary success of a true toad (Bufo gargarizans)

We present a high-quality genome assembly for the Asiatic toad (Bufo gargarizans) and explore the evolution of several large gene families in amphibians. With a large genome assembly size of 4.55 Gb, the chromosome-scale assembly includes 747 scaffolds with an N50 of 539.8 Mb and 1.79% gaps. Long terminal repeats (LTRs) constitute a high proportion of the genome and their expansion is a key contributor to the inflated genome size in this species. This is very different from other small amphibian genomes, but similar to that of the enormous axolotl genome. The genome retains a large number of duplicated genes, with tandem (TD) and proximal duplications (PD) the predominant mode of duplication. A total of 122 gene families have undergone significant expansion and were mainly enriched in sensory perception of smell and bitter taste. The CYP2C subfamily, which plays an important role in metabolic detoxification, specifically expanded via TD and PD in the Asiatic toad and the cane toad (true toads). Most of Na⁺ /K⁺ -ATPase genes experienced accelerated evolution along Bufonid lineages and two amino acid sites involving toad-toxin resistance were found to experience positive selection. We also revealed a dynamic evolution of olfactory and vomeronasal receptor gene families which was likely driven by the water-to-land transition. The high-quality genome of the Asiatic toad will provide a solid foundation to understand the genetic basis of its many biological processes.

Single-parent expression drives dynamic gene expression complementation in maize hybrids

Single-parent expression (SPE) is defined as gene expression in only one of the two parents. SPE can arise from differential expression between parental alleles, termed non-presence/absence (non-PAV) SPE, or from the physical absence of a gene in one parent, termed PAV SPE. We used transcriptome data of diverse Zea mays (maize) inbreds and hybrids, including 401 samples from five different tissues, to test for differences between these types of SPE genes. Although commonly observed, SPE is highly genotype and tissue specific. A positive correlation was observed between the genetic distance of the two inbred parents and the number of SPE genes identified. Regulatory analysis showed that PAV SPE and non-PAV SPE genes are mainly regulated by cis effects, with a small fraction under trans regulation. Polymorphic transposable element insertions in promoter sequences contributed to the high level of cis regulation for PAV SPE and non-PAV SPE genes. PAV SPE genes were more frequently expressed in hybrids than non-PAV SPE genes. The expression of parentally silent alleles in hybrids of non-PAV SPE genes was relatively rare but occurred in most hybrids. Non-PAV SPE genes with expression of the silent allele in hybrids are more likely to exhibit above high parent expression level than hybrids that do not express the silent allele, leading to non-additive expression. This study provides a comprehensive understanding of the nature of non-PAV SPE and PAV SPE genes and their roles in gene expression complementation in maize hybrids.

The Gene Balance Hypothesis: Epigenetics and Dosage Effects in Plants

Dosage effects in plants are caused by changes in the copy number of chromosomes, segments of chromosomes, or multiples of individual genes. Genes often exhibit a dosage effect in which the amount of product is closely correlated with the number of copies present. However, when larger segments of chromosomes are varied, there are trans-acting effects across the genome that are unleashed that modulate gene expression in cascading effects. These appear to be mediated by the stoichiometric relationship of gene regulatory machineries. There are both positive and negative modulations of target gene expression, but the latter is the plurality effect. When this inverse effect is combined with a dosage effect, compensation for a gene can occur in which its expression is similar to the normal diploid regardless of the change in chromosomal dosage. In contrast, changing the whole genome in a polyploidy series has fewer relative effects as the stoichiometric relationship is not disrupted. Together, these observations suggest that the stoichiometry of gene regulation is important as a reflection of the mode of assembly of the individual subunits involved in the effective regulatory macromolecular complexes. This principle has implications for gene expression mechanisms, quantitative trait genetics, and the evolution of genes depending on the mode of duplication, either segmentally or via whole-genome duplication.

Genes and gene models, an important distinction

Genome sequencing has fundamentally changed how plant biologists think about genes. All or nearly all genes can ultimately be associated with a gene model. However, many gene models appear to play little or no role in the traits of an organism. A range of structural, molecular, population and evolutionary features all show a separation between genes with known phenotypes and the overall set of annotated gene models. These different features could be combined to develop models to distinguish the genes that determine the traits of plants from the subset gene other annotated gene models which are unlikely to play a role in doing so. Efforts to identify the subset of annotated gene models likely involved in specifying the characteristics of plants would help aid a wide range of researchers.

Plant pan-genomes are the new reference

Recent years have seen a surge in plant genome sequencing projects and the comparison of multiple related individuals. The high degree of genomic variation observed led to the realization that single reference genomes do not represent the diversity within a species, and led to the expansion of the pan-genome concept. Pan-genomes represent the genomic diversity of a species and includes core genes, found in all individuals, as well as variable genes, which are absent in some individuals. Variable gene annotations often show similarities across plant species, with genes for biotic and abiotic stress commonly enriched within variable gene groups. Here we review the growth of pan-genomics in plants, explore the origins of gene presence and absence variation, and show how pan-genomes can support plant breeding and evolution studies.

Genes Encoding Recognition of the Cladosporium fulvum Effector Protein Ecp5 Are Encoded at Several Loci in the Tomato Genome

The molecular interactions between tomato and Cladosporium fulvum have been an important model for molecular plant pathology. Complex genetic loci on tomato chromosomes 1 and 6 harbor genes for resistance to Cladosporium fulvum, encoding receptor like-proteins that perceive distinct Cladosporium fulvum effectors and trigger plant defenses. Here, we report classical mapping strategies for loci in tomato accessions that respond to Cladosporium fulvum effector Ecp5, which is very sequence-monomorphic. We screened 139 wild tomato accessions for an Ecp5-induced hypersensitive response, and in five accessions, the Ecp5-induced hypersensitive response segregated as a monogenic trait, mapping to distinct loci in the tomato genome. We identified at least three loci on chromosomes 1, 7 and 12 that harbor distinct Cf-Ecp5 genes in four different accessions. Our mapping showed that the Cf-Ecp5 in Solanum pimpinellifolium G1.1161 is located at the Milky Way locus. The Cf-Ecp5 in Solanum pimpinellifolium LA0722 was mapped to the bottom arm of chromosome 7, while the Cf-Ecp5 genes in Solanum lycopersicum Ontario 7522 and Solanum pimpinellifolium LA2852 were mapped to the same locus on the top arm of chromosome 12. Bi-parental crosses between accessions carrying distinct Cf-Ecp5 genes revealed putative genetically unlinked suppressors of the Ecp5-induced hypersensitive response. Our mapping also showed that Cf-11 is located on chromosome 11, close to the Cf-3 locus. The Ecp5-induced hypersensitive response is widely distributed within tomato species and is variable in strength. This novel example of convergent evolution could be used for choosing different functional Cf-Ecp5 genes according to individual plant breeding needs.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Robust non-syntenic gene expression patterns in diverse maize hybrids during root development

Distantly related maize (Zea mays L.) inbred lines exhibit an exceptional degree of structural genomic diversity, which is probably unique among plants. This study systematically investigated the developmental and genotype-dependent regulation of the primary root transcriptomes of a genetically diverse panel of maize F1-hybrids and their parental inbred lines. While we observed substantial transcriptomic changes during primary root development, we demonstrated that hybrid-associated gene expression patterns, including differential, non-additive, and allele-specific transcriptome profiles, are particularly robust to these developmental fluctuations. For instance, differentially expressed genes with preferential expression in hybrids were highly conserved during development in comparison to their parental counterparts. Similarly, in hybrids a major proportion of non-additively expressed genes with expression levels between the parental values were particularly conserved during development. Importantly, in these expression patterns non-syntenic genes that evolved after the separation of the maize and sorghum lineages were systemically enriched. Furthermore, non-syntenic genes were substantially linked to the conservation of all surveyed gene expression patterns during primary root development. Among all F1-hybrids, between ~40% of the non-syntenic genes with unexpected allelic expression ratios and ~60% of the non-syntenic differentially and non-additively expressed genes were conserved and therefore robust to developmental changes. Hence, the enrichment of non-syntenic genes during primary root development might be involved in the developmental adaptation of maize roots and thus the superior performance of hybrids.

Molecular evolution of CRH and CRHR subfamily before the evolutionary origin of vertebrate

Corticotropin-releasing hormone (CRH) is well-cited for its important role in governing the stress responses via neuroendocrine system in vertebrates. After the identification of homologs of CRH receptor (CRHR) in both deuterostome and arthropod lineages, it was suggested that the ancestral homolog of CRH-CRHR molecular system is present in the bilaterian. However, homolog sequences from arthropods differ considerably from vertebrate CRH-like peptide sequences. Due to the significant difference between the biological system, as well as the gene regulatory network, of protostome and that of vertebrate, physiological studies on the protostomes may not provide important insight into the evolutionary history of vertebrate CRH system, while tunicate and amphioxus, two close relatives to vertebrate, which have diverged before two rounds of whole genome duplication (2WGDs) do. Given the identification of amphioxus CRH-like peptide by our group, this review aims to reexamine the current hypotheses on the evolution of CRH subfamily. It is generally accepted that paralogs of CRH and CRHR have been produced through 2WGDs, which occurred during the early vertebrate evolution. The identification of a single crh-like gene in amphioxi and tunicates by in silico search and the presence of two paralogons with a total of 5 crh-like genes in gnathostomes has shown that an additional duplication event might have happened to the ancestral crh-like gene before 2WGDs. On the other hand, the evolution of crhr gene subfamily appears to be mainly influenced by 2WGDs and only two receptor genes have been retained in the genomes of jawed vertebrates.

Alignment of Rutaceae Genomes Reveals Lower Genome Fractionation Level Than Eudicot Genomes Affected by Extra Polyploidization

Owing to their nutritional and commercial values, the genomes of several citrus plants have been sequenced, and the genome of one close relative in the Rutaceae family, atalantia (Atalantia buxifolia), has also been sequenced. Here, we show a family-level comparative analysis of Rutaceae genomes. By using grape as the outgroup and checking cross-genome gene collinearity, we systematically performed a hierarchical and event-related alignment of Rutaceae genomes, and produced a gene list defining homologous regions based on ancestral polyploidization or speciation. We characterized genome fractionation resulting from gene loss or relocation, and found that erosion of gene collinearity could largely be described by a geometric distribution. Moreover, we found that well-assembled Rutaceae genomes retained significantly more genes (65-82%) than other eudicots affected by recursive polyploidization. Additionally, we showed divergent evolutionary rates among Rutaceae plants, with sweet orange evolving faster than others, and by performing evolutionary rate correction, re-dated major evolutionary events during their evolution. We deduced that the divergence between the Rutaceae family and grape occurred about 81.15-91.74 million years ago (mya), while the split between citrus and atalantia plants occurred <10 mya. In addition, we showed that polyploidization led to a copy number expansion of key gene families contributing to the biosynthesis of vitamin C. Overall, the present effort provides an important comparative genomics resource and lays a foundation to understand the evolution and functional innovation of Rutaceae genomes.

A Dead Gene Walking: Convergent Degeneration of a Clade of MADS-Box Genes in Crucifers

Genes are "born," and eventually they "die." These processes shape the phenotypic evolution of organisms and are hence of great biological interest. If genes die in plants, they generally do so quite rapidly. Here, we describe the fate of GOA-like genes that evolve in a dramatically different manner. GOA-like genes belong to the subfamily of Bsister genes of MIKC-type MADS-box genes. Typical MIKC-type genes encode conserved transcription factors controlling plant development. We show that ABS-like genes, a clade of Bsister genes, are indeed highly conserved in crucifers (Brassicaceae) maintaining the ancestral function of Bsister genes in ovule and seed development. In contrast, their closest paralogs, the GOA-like genes, have been undergoing convergent gene death in Brassicaceae. Intriguingly, erosion of GOA-like genes occurred after millions of years of coexistence with ABS-like genes. We thus describe Delayed Convergent Asymmetric Degeneration, a so far neglected but possibly frequent pattern of duplicate gene evolution that does not fit classical scenarios. Delayed Convergent Asymmetric Degeneration of GOA-like genes may have been initiated by a reduction in the expression of an ancestral GOA-like gene in the stem group of Brassicaceae and driven by dosage subfunctionalization. Our findings have profound implications for gene annotations in genomics, interpreting patterns of gene evolution and using genes in phylogeny reconstructions of species.

Parallelism and convergence in post-domestication adaptation in cereal grasses

The selection of desirable traits in crops during domestication has been well studied. Many crops share a suite of modified phenotypic characteristics collectively known as the domestication syndrome. In this sense, crops have convergently evolved. Previous work has demonstrated that, at least in some instances, convergence for domestication traits has been achieved through parallel molecular means. However, both demography and selection during domestication may have placed limits on evolutionary potential and reduced opportunities for convergent adaptation during post-domestication migration to new environments. Here we review current knowledge regarding trait convergence in the cereal grasses and consider whether the complexity and dynamism of cereal genomes (e.g., transposable elements, polyploidy, genome size) helped these species overcome potential limitations owing to domestication and achieve broad subsequent adaptation, in many cases through parallel means. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

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.

Genome-wide identification and functional analysis of the splicing component SYF2/NTC31/p29 across different plant species

This study systematically identifies plant SYF2/NTC31/p29 genes from 62 plant species by a combinatory bioinformatics approach, revealing the importance of this gene family in phylogenetics, duplication, transcriptional, and post-transcriptional regulation. Alternative splicing is a post-transcriptional regulatory mechanism, which is critical for plant development and stress responses. The entire process is strictly attenuated by a complex of splicing-related proteins, designated splicing factors. Human p29, also referred to as synthetic lethal with cdc forty 2 (SYF2) or the NineTeen complex 31 (NTC31), is a core protein found in the NTC complex of humans and yeast. This splicing factor participates in a variety of biological processes, including DNA damage repair, control of the cell cycle, splicing, and tumorigenesis. However, its function in plants has been seldom reported. Thus, we have systematically identified 89 putative plant SYF2s from 62 plant species among the deposited entries in the Phytozome database. The phylogenetic relationships and evolutionary history among these plant SYF2s were carefully examined. The results revealed that plant SYF2s exhibited distinct patterns regarding their gene structure, promoter sequences, and expression levels, suggesting their functional diversity in response to developmental cues or stress treatments. Although local duplication events, such as tandem duplication and retrotransposition, were found among several plant species, most of the plant species contained only one copy of SYF2, suggesting the existence of additional mechanisms to confer duplication resistance. Further investigation using the model dicot and monocot representatives Arabidopsis and rice SYF2s indicated that the splicing pattern and resulting protein isoforms might play an alternative role in the functional diversity.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

Different Modes of Gene Duplication Show Divergent Evolutionary Patterns and Contribute Differently to the Expansion of Gene Families Involved in Important Fruit Traits in Pear (Pyrus bretschneideri)

Pear is an important fruit crop of the Rosaceae family and has experienced two rounds of ancient whole-genome duplications (WGDs). However, whether different types of gene duplications evolved differently after duplication remains unclear in the pear genome. In this study, we identified the different modes of gene duplication in pear. Duplicate genes derived from WGD, tandem, proximal, retrotransposed, DNA-based transposed or dispersed duplications differ in genomic distribution, gene features, selection pressure, expression divergence, regulatory divergence and biological roles. Widespread sequence, expression and regulatory divergence have occurred between duplicate genes over the 30-45 million years of evolution after the recent genome duplication in pear. The retrotransposed genes show relatively higher expression and regulatory divergence than other gene duplication modes. In contrast, WGD genes underwent a slower sequence divergence and may be influenced by abundant gene conversion events. Moreover, the different classes of duplicate genes exhibited biased functional roles. We also investigated the evolution and expansion patterns of the gene families involved in sugar and organic acid metabolism pathways, which are closely related to the fruit quality and taste in pear. Single-gene duplications largely account for the extensive expansion of gene families involved in the sorbitol metabolism pathway in pear. Gene family expansion was also detected in the sucrose metabolism pathway and tricarboxylic acid cycle pathways. Thus, this study provides insights into the evolutionary fates of duplicated genes.

Single-Parent Expression Is a General Mechanism Driving Extensive Complementation of Non-syntenic Genes in Maize Hybrids

Maize (Zea mays L.) displays an exceptional degree of structural genomic diversity [1, 2]. In addition, variation in gene expression further contributes to the extraordinary phenotypic diversity and plasticity of maize. This study provides a systematic investigation on how distantly related homozygous maize inbred lines affect the transcriptomic plasticity of their highly heterozygous F₁ hybrids. The classical dominance model of heterosis explains the superiority of hybrid plants by the complementation of deleterious parental alleles by superior alleles of the second parent at many loci [3]. Genes active in one inbred line but inactive in another represent an extreme instance of allelic diversity defined as single-parent expression [4]. We observed on average ∼1,000 such genes in all inbred line combinations during primary root development. These genes consistently displayed expression complementation (i.e., activity) in their hybrid progeny. Consequently, extreme expression complementation is a general mechanism that results on average in ∼600 additionally active genes and their encoded biological functions in hybrids. The modern maize genome is complemented by a set of non-syntenic genes, which emerged after the separation of the maize and sorghum lineages and lack syntenic orthologs in any other grass species [5]. We demonstrated that non-syntenic genes are the driving force of gene expression complementation in hybrids. Among those, the highly diversified families of bZIP and bHLH transcription factors [6] are systematically overrepresented. In summary, extreme gene expression complementation extensively shapes the transcriptomic plasticity of maize hybrids and might therefore be one factor controlling the developmental plasticity of hybrids.

Extensive gene content variation in the Brachypodium distachyon pan-genome correlates with population structure

While prokaryotic pan-genomes have been shown to contain many more genes than any individual organism, the prevalence and functional significance of differentially present genes in eukaryotes remains poorly understood. Whole-genome de novo assembly and annotation of 54 lines of the grass Brachypodium distachyon yield a pan-genome containing nearly twice the number of genes found in any individual genome. Genes present in all lines are enriched for essential biological functions, while genes present in only some lines are enriched for conditionally beneficial functions (e.g., defense and development), display faster evolutionary rates, lie closer to transposable elements and are less likely to be syntenic with orthologous genes in other grasses. Our data suggest that differentially present genes contribute substantially to phenotypic variation within a eukaryote species, these genes have a major influence in population genetics, and transposable elements play a key role in pan-genome evolution.

The role of differential gene content in the evolution and function of eukaryotic genomes remains poorly explored. Here the authors assemble and annotate the Brachypodium distachyon pan-genome consisting of 54 diverse lines and reveal the differential present genes as a major driver of phenotypic variation.

Genome-wide characterization and expression profiling of NAC transcription factor genes under abiotic stresses in radish (Raphanus sativus L.)

NAC (NAM, no apical meristem; ATAF, Arabidopsis transcription activation factor and CUC, cup-shaped cotyledon) proteins are among the largest transcription factor (TF) families playing fundamental biological processes, including cell expansion and differentiation, and hormone signaling in response to biotic and abiotic stresses. In this study, 172 RsNACs comprising 17 membrane-bound members were identified from the whole radish genome. In total, 98 RsNAC genes were non-uniformly distributed across the nine radish chromosomes. In silico analysis revealed that expression patterns of several NAC genes were tissue-specific such as a preferential expression in roots and leaves. In addition, 21 representative NAC genes were selected to investigate their responses to heavy metals (HMs), salt, heat, drought and abscisic acid (ABA) stresses using real-time polymerase chain reaction (RT-qPCR). As a result, differential expressions among these genes were identified where RsNAC023 and RsNAC080 genes responded positively to all stresses except ABA, while RsNAC145 responded more actively to salt, heat and drought stresses compared with other genes. The results provides more valuable information and robust candidate genes for future functional analysis for improving abiotic stress tolerances in radish.

Patterns and Consequences of Subgenome Differentiation Provide Insights into the Nature of Paleopolyploidy in Plants

Polyploidy is an important feature of plant genomes, but the nature of many polyploidization events remains to be elucidated. Here, we demonstrate that the evolutionary fates of the subgenomes in maize (Zea mays) and soybean (Glycine max) have followed different trajectories. One subgenome has been subject to relaxed selection, lower levels of gene expression, higher rates of transposable element accumulation, more small interfering RNAs and DNA methylation around genes, and higher rates of gene loss in maize, whereas none of these features were observed in soybean. Nevertheless, individual gene pairs exhibit differentiation with respect to these features in both species. In addition, we observed a higher number of chromosomal rearrangements and higher frequency of retention of duplicated genes in soybean than in maize. Furthermore, soybean "singletons" were found to be more frequently tandemly duplicated than "duplicates" in soybean, which may, to some extent, counteract the genome imbalance caused by gene loss. We propose that unlike in maize, in which two subgenomes were distinct prior to the allotetraploidization event and thus experienced global differences in selective constraints, in soybean, the two subgenomes were far less distinct prior to polyploidization, such that individual gene pairs, rather than subgenomes, experienced stochastic differences over longer periods of time, resulting in retention of the majority of duplicates.

Reciprocally Retained Genes in the Angiosperm Lineage Show the Hallmarks of Dosage Balance Sensitivity

In several organisms, particular functional categories of genes, such as regulatory and complex-forming genes, are preferentially retained after whole-genome multiplications but rarely duplicate through small-scale duplication, a pattern referred to as reciprocal retention. This peculiar duplication behavior is hypothesized to stem from constraints on the dosage balance between the genes concerned and their interaction context. However, the evidence for a relationship between reciprocal retention and dosage balance sensitivity remains fragmentary. Here, we identified which gene families are most strongly reciprocally retained in the angiosperm lineage and studied their functional and evolutionary characteristics. Reciprocally retained gene families exhibit stronger sequence divergence constraints and lower rates of functional and expression divergence than other gene families, suggesting that dosage balance sensitivity is a general characteristic of reciprocally retained genes. Gene families functioning in regulatory and signaling processes are much more strongly represented at the top of the reciprocal retention ranking than those functioning in multiprotein complexes, suggesting that regulatory imbalances may lead to stronger fitness effects than classical stoichiometric protein complex imbalances. Finally, reciprocally retained duplicates are often subject to dosage balance constraints for prolonged evolutionary times, which may have repercussions for the ease with which genome multiplications can engender evolutionary innovation.

Maize RNA PolIV affects the expression of genes with nearby TE insertions and has a genome-wide repressive impact on transcription

BACKGROUND

RNA-directed DNA methylation (RdDM) is a plant-specific epigenetic process that relies on the RNA polymerase IV (Pol IV) for the production of 24 nucleotide small interfering RNAs (siRNA) that guide the cytosine methylation and silencing of genes and transposons. Zea mays RPD1/RMR6 gene encodes the largest subunit of Pol IV and is required for normal plant development, paramutation, transcriptional repression of certain transposable elements (TEs) and transcriptional regulation of specific alleles.

RESULTS

In this study we applied a total RNA-Seq approach to compare the B73 and rpd1/rmr6 leaf transcriptomes. Although previous studies indicated that loss of siRNAs production in RdDM mutants provokes a strong loss of CHH DNA methylation but not massive gene or TEs transcriptional activation in both Arabidopsis and maize, our total RNA-Seq analysis of rpd1/rmr6 transcriptome reveals that loss of Pol IV activity causes a global increase in the transcribed fraction of the maize genome. Our results point to the genes with nearby TE insertions as being the most strongly affected by Pol IV-mediated gene silencing. TEs modulation of nearby gene expression is linked to alternative methylation profiles on gene flanking regions, and these profiles are strictly dependent on specific characteristics of the TE member inserted. Although Pol IV is essential for the biogenesis of siRNAs, the genes with associated siRNA loci are less affected by the pol IV mutation.

CONCLUSIONS

This deep and integrated analysis of gene expression, TEs distribution, smallRNA targeting and DNA methylation levels, reveals that loss of Pol IV activity globally affects genome regulation, pointing at TEs as modulator of nearby gene expression and indicating the existence of multiple level epigenetic silencing mechanisms. Our results also suggest a predominant role of the Pol IV-mediated RdDM pathway in genome dominance regulation, and subgenome stability and evolution in maize.

Recruitment of pre-existing networks during the evolution of C4 photosynthesis

During C4 photosynthesis, CO2 is concentrated around the enzyme RuBisCO. The net effect is to reduce photorespiration while increasing water and nitrogen use efficiencies. Species that use C4 photosynthesis have evolved independently from their C3 ancestors on more than 60 occasions. Along with mimicry and the camera-like eye, the C4 pathway therefore represents a remarkable example of the repeated evolution of a highly complex trait. In this review, we provide evidence that the polyphyletic evolution of C4 photosynthesis is built upon pre-existing metabolic and genetic networks. For example, cells around veins of C3 species show similarities to those of the C4 bundle sheath in terms of C4 acid decarboxylase activity and also the photosynthetic electron transport chain. Enzymes of C4 photosynthesis function together in gluconeogenesis during early seedling growth of C3Arabidopsis thaliana Furthermore, multiple C4 genes appear to be under control of both light and chloroplast signals in the ancestral C3 state. We, therefore, hypothesize that relatively minor rewiring of pre-existing genetic and metabolic networks has facilitated the recurrent evolution of this trait. Understanding how these changes are likely to have occurred could inform attempts to install C4 traits into C3 crops.This article is part of the themed issue 'Enhancing photosynthesis in crop plants: targets for improvement'.

Phylogenomic Synteny Network Analysis of MADS-Box Transcription Factor Genes Reveals Lineage-Specific Transpositions, Ancient Tandem Duplications, and Deep Positional Conservation

Conserved genomic context provides critical information for comparative evolutionary analysis. With the increase in numbers of sequenced plant genomes, synteny analysis can provide new insights into gene family evolution. Here, we exploit a network analysis approach to organize and interpret massive pairwise syntenic relationships. Specifically, we analyzed synteny networks of the MADS-box transcription factor gene family using 51 completed plant genomes. In combination with phylogenetic profiling, several novel evolutionary patterns were inferred and visualized from synteny network clusters. We found lineage-specific clusters that derive from transposition events for the regulators of floral development (APETALA3 and PI) and flowering time (FLC) in the Brassicales and for the regulators of root development (AGL17) in Poales. We also identified two large gene clusters that jointly encompass many key phenotypic regulatory Type II MADS-box gene clades (SEP1, SQUA, TM8, SEP3, FLC, AGL6, and TM3). Gene clustering and gene trees support the idea that these genes are derived from an ancient tandem gene duplication that likely predates the radiation of the seed plants and then expanded by subsequent polyploidy events. We also identified angiosperm-wide conservation of synteny of several other less studied clades. Combined, these findings provide new hypotheses for the genomic origins, biological conservation, and divergence of MADS-box gene family members.

ARGOS8 variants generated by CRISPR-Cas9 improve maize grain yield under field drought stress conditions

Maize ARGOS8 is a negative regulator of ethylene responses. A previous study has shown that transgenic plants constitutively overexpressing ARGOS8 have reduced ethylene sensitivity and improved grain yield under drought stress conditions. To explore the targeted use of ARGOS8 native expression variation in drought-tolerant breeding, a diverse set of over 400 maize inbreds was examined for ARGOS8 mRNA expression, but the expression levels in all lines were less than that created in the original ARGOS8 transgenic events. We then employed a CRISPR-Cas-enabled advanced breeding technology to generate novel variants of ARGOS8. The native maize GOS2 promoter, which confers a moderate level of constitutive expression, was inserted into the 5'-untranslated region of the native ARGOS8 gene or was used to replace the native promoter of ARGOS8. Precise genomic DNA modification at the ARGOS8 locus was verified by PCR and sequencing. The ARGOS8 variants had elevated levels of ARGOS8 transcripts relative to the native allele and these transcripts were detectable in all the tissues tested, which was the expected results using the GOS2 promoter. A field study showed that compared to the WT, the ARGOS8 variants increased grain yield by five bushels per acre under flowering stress conditions and had no yield loss under well-watered conditions. These results demonstrate the utility of the CRISPR-Cas9 system in generating novel allelic variation for breeding drought-tolerant crops.

Stability of Single-Parent Gene Expression Complementation in Maize Hybrids upon Water Deficit Stress

Heterosis is the superior performance of F1 hybrids compared with their homozygous, genetically distinct parents. In this study, we monitored the transcriptomic divergence of the maize (Zea mays) inbred lines B73 and Mo17 and their reciprocal F1 hybrid progeny in primary roots under control and water deficit conditions simulated by polyethylene glycol treatment. Single-parent expression (SPE) of genes is an extreme instance of gene expression complementation, in which genes are active in only one of two parents but are expressed in both reciprocal hybrids. In this study, 1,997 genes only expressed in B73 and 2,024 genes only expressed in Mo17 displayed SPE complementation under control and water deficit conditions. As a consequence, the number of active genes in hybrids exceeded the number of active genes in the parental inbred lines significantly independent of treatment. SPE patterns were substantially more stable to expression changes by water deficit treatment than other genotype-specific expression profiles. While, on average, 75% of all SPE patterns were not altered in response to polyethylene glycol treatment, only 17% of the remaining genotype-specific expression patterns were not changed by water deficit. Nonsyntenic genes that lack syntenic orthologs in other grass species, and thus evolved late in the grass lineage, were significantly overrepresented among SPE genes. Hence, the significant overrepresentation of nonsyntenic genes among SPE patterns and their stability under water limitation might suggest a function of these genes during the early developmental manifestation of heterosis under fluctuating environmental conditions in hybrid progeny of the inbred lines B73 and Mo17.