DNA sequence evidence for the segmental allotetraploid origin of maize

Accurate Inference of the Polyploid Continuum Using Forward-Time Simulations

Multiple rounds of whole-genome duplication (WGD) followed by diploidization have occurred throughout the evolutionary history of angiosperms. Much work has been done to model the genomic consequences and evolutionary significance of WGD. While researchers have historically modeled polyploids as either allopolyploids or autopolyploids, the variety of natural polyploids span a continuum of differentiation across multiple parameters, such as the extent of polysomic versus disomic inheritance, and the degree of genetic differentiation between the ancestral lineages. Here we present a forward-time polyploid genome evolution simulator called SpecKS. SpecKS models polyploid speciation as originating from a 2D continuum, whose dimensions account for both the level of genetic differentiation between the ancestral parental genomes, as well the time lag between ancestral speciation and their subsequent reunion in the derived polyploid. Using extensive simulations, we demonstrate that changes in initial conditions along either dimension of the 2D continuum deterministically affect the shape of the Ks histogram. Our findings indicate that the error in the common method of estimating WGD time from the Ks histogram peak scales with the degree of allopolyploidy, and we present an alternative, accurate estimation method that is independent of the degree of allopolyploidy. Lastly, we use SpecKS to derive tests that infer both the lag time between parental divergence and WGD time, and the diversity of the ancestral species, from an input Ks histogram. We apply the latter test to transcriptomic data from over 200 species across the plant kingdom, the results of which are concordant with the prevailing theory that the majority of angiosperm lineages are derived from diverse parental genomes and may be of allopolyploid origin.

Variant calling in polyploids for population and quantitative genetics

Advancements in genome assembly and sequencing technology have made whole genome sequence (WGS) data and reference genomes accessible to study polyploid species. Compared to popular reduced-representation sequencing approaches, the genome-wide coverage and greater marker density provided by WGS data can greatly improve our understanding of polyploid species and polyploid biology. However, biological features that make polyploid species interesting also pose challenges in read mapping, variant identification, and genotype estimation. Accounting for characteristics in variant calling like allelic dosage uncertainty, homology between subgenomes, and variance in chromosome inheritance mode can reduce errors. Here, I discuss the challenges of variant calling in polyploid WGS data and discuss where potential solutions can be integrated into a standard variant calling pipeline.

Physiology, genomics, and evolutionary aspects of desert plants

BACKGROUND

Despite the exposure to arid environmental conditions across the globe ultimately hampering the sustainability of the living organism, few plant species are equipped with several unique genotypic, biochemical, and physiological features to counter such harsh conditions. Physiologically, they have evolved with reduced leaf size, spines, waxy cuticles, thick leaves, succulent hydrenchyma, sclerophyll, chloroembryo, and photosynthesis in nonfoliar and other parts. At the biochemical level, they are evolved to perform efficient photosynthesis through Crassulacean acid metabolism (CAM) and C4 pathways with the formation of oxaloacetic acid (Hatch-Slack pathway) instead of the C3 pathway. Additionally, comparative genomics with existing data provides ample evidence of the xerophytic plants' positive selection to adapt to the arid environment. However, adding more high-throughput sequencing of xerophyte plant species is further required for a comparative genomic study toward trait discovery related to survival. Learning from the mechanism to survive in harsh conditions could pave the way to engineer crops for future sustainable agriculture.

AIM OF THE REVIEW

The distinct physiology of desert plants allows them to survive in harsh environments. However, the genomic composition also contributes significantly to this and requires great attention. This review emphasizes the physiological and genomic adaptation of desert plants. Other important parameters, such as desert biodiversity and photosynthetic strategy, are also discussed with recent progress in the field. Overall, this review discusses the different features of desert plants, which prepares them for harsh conditions intending to translate knowledge to engineer plant species for sustainable agriculture.

KEY SCIENTIFIC CONCEPTS OF REVIEW

This review comprehensively presents the physiology, molecular mechanism, and genomics of desert plants aimed towards engineering a sustainable crop.

Modeling transposable elements dynamics during polyploidization in plants

In this work we study the proliferation of transposable elements (TEs) and the epigenetic response of plants during the process of polyploidization. Through a deterministic model, expanding on our previous work on TE proliferation under epigenetic regulation, we study the long-term TE distribution and TE stability in the subgenomes of both autopolyploids and allopolyploids. We also explore different small-interfering RNA (siRNA) action modes on the subgenomes, including a model where siRNAs are not directed to specific genomes and one where siRNAs are directed - i.e. more active - in subgenomes with higher TE loads. In the autopolyploid case, we find long-term stable equilbria that tend to equilibrate the number of active TEs between subgenomes. In the allopolyploid case, directed siRNA action is fundamental to avoid a "winner takes all" outcome of the competition between the TEs in the different subgenomes. We also show that decaying oscillations in the number of TEs occur naturally in all cases, perhaps explaining some of the observed features of 'genomic shock' after hybridization events, and that the balance in the dynamics of the different types of siRNA is determinant for the synchronization of these oscillations.

Polyploid speciation in Zea (Poaceae): cytogenetic insights

MAIN CONCLUSION

The analysis of meiotic pairing affinities and genomic formulae in species and hybrids of Zea allowed us to speculate an evolutionary model to recreate the ancient polyploidization of maize and allied species. The meiotic pairing affinities and the genomic formulae analysis in Zea species and hybrids obtained in new and previous crosses, together with the molecular data known in the genus, allowed us to speculate an evolutionary model to attempt to recreate the ancient polyploidization process of Zea species. We propose that x = 5 semispecies are the ancestors of all modern species of the genus. The complex evolutionary process that originated the different taxa could be included hybridization between sympatric diploid ancestral semispecies (2n = 10) and recurrent duplication of the hybrid chromosome number, resulting in distinct auto- and allopolyploids. After the merger and doubling of independent genomes would have undergone cytological and genetical diploidization, implying revolutionary changes in genome organization and genic balance processes. Based on the meiotic behaviour of the 2n = 30 hybrids, that showed homoeology between the A subgenomes of all parental species, we propose that this subgenome A would be pivotal in all the species and would have conserved the rDNA sequences and the pairing regulator locus (PrZ). In the hypothetical model postulated here, the ancestral semispecies with the pivotal subgenome A would have had a wide geographic distribution, co-occurring and hybridizing with the semispecies harbouring B subgenomes, thus enabling sympatric speciation.

Demographic history inference and the polyploid continuum

Polyploidy is an important generator of evolutionary novelty across diverse groups in the Tree of Life, including many crops. However, the impact of whole-genome duplication depends on the mode of formation: doubling within a single lineage (autopolyploidy) versus doubling after hybridization between two different lineages (allopolyploidy). Researchers have historically treated these two scenarios as completely separate cases based on patterns of chromosome pairing, but these cases represent ideals on a continuum of chromosomal interactions among duplicated genomes. Understanding the history of polyploid species thus demands quantitative inferences of demographic history and rates of exchange between subgenomes. To meet this need, we developed diffusion models for genetic variation in polyploids with subgenomes that cannot be bioinformatically separated and with potentially variable inheritance patterns, implementing them in the dadi software. We validated our models using forward SLiM simulations and found that our inference approach is able to accurately infer evolutionary parameters (timing, bottleneck size) involved with the formation of auto- and allotetraploids, as well as exchange rates in segmental allotetraploids. We then applied our models to empirical data for allotetraploid shepherd's purse (Capsella bursa-pastoris), finding evidence for allelic exchange between the subgenomes. Taken together, our model provides a foundation for demographic modeling in polyploids using diffusion equations, which will help increase our understanding of the impact of demography and selection in polyploid lineages.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

Genome-wide comparative analysis of the valine glutamine motif containing genes in four Ipomoea species

BACKGROUND

Genes with valine glutamine (VQ) motifs play an essential role in plant growth, development, and resistance to biotic and abiotic stresses. However, little information on the VQ genes in sweetpotato and other Ipomoea species is available.

RESULTS

This study identified 55, 58, 50 and 47 VQ genes from sweetpotato (I. batatas), I.triflida, I. triloba and I. nil, respectively. The phylogenetic analysis revealed that the VQ genes formed eight clades (I-VII), and the members in the same group exhibited similar exon-intron structure and conserved motifs distribution. The distribution of the VQ genes among the chromosomes of Ipomoea species was disproportional, with no VQ genes mapped on a few of each species' chromosomes. Duplication analysis suggested that segmental duplication significantly contributes to their expansion in sweetpotato, I.trifida, and I.triloba, while the segmental and tandem duplication contributions were comparable in I.nil. Cis-regulatory elements involved in stress responses, such as W-box, TGACG-motif, CGTCA-motif, ABRE, ARE, MBS, TCA-elements, LTR, and WUN-motif, were detected in the promoter regions of the VQ genes. A total of 30 orthologous groups were detected by syntenic analysis of the VQ genes. Based on the analysis of RNA-seq datasets, it was found that the VQ genes are expressed distinctly among different tissues and hormone or stress treatments. A total of 40 sweetpotato differentially expressed genes (DEGs) refer to biotic (sweetpotato stem nematodes and Ceratocystis fimbriata pathogen infection) or abiotic (cold, salt and drought) stress treatments were detected. Moreover, IbVQ8, IbVQ25 and IbVQ44 responded to the five stress treatments and were selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis.

CONCLUSIONS

Our study may provide new insights into the evolution of VQ genes in the four Ipomoea genomes and contribute to the future molecular breeding of sweetpotatoes.

Conserved noncoding sequences and de novo Mutator insertion alleles are imprinted in maize

Genomic imprinting is an epigenetic phenomenon in which differential allele expression occurs in a parent-of-origin-dependent manner. Imprinting in plants is tightly linked to transposable elements (TEs), and it has been hypothesized that genomic imprinting may be a consequence of demethylation of TEs. Here, we performed high-throughput sequencing of ribonucleic acids from four maize (Zea mays) endosperms that segregated newly silenced Mutator (Mu) transposons and identified 110 paternally expressed imprinted genes (PEGs) and 139 maternally expressed imprinted genes (MEGs). Additionally, two potentially novel paternally suppressed MEGs are associated with de novo Mu insertions. In addition, we find evidence for parent-of-origin effects on expression of 407 conserved noncoding sequences (CNSs) in maize endosperm. The imprinted CNSs are largely localized within genic regions and near genes, but the imprinting status of the CNSs are largely independent of their associated genes. Both imprinted CNSs and PEGs have been subject to relaxed selection. However, our data suggest that although MEGs were already subject to a higher mutation rate prior to their being imprinted, imprinting may be the cause of the relaxed selection of PEGs. In addition, although DNA methylation is lower in the maternal alleles of both the maternally and paternally expressed CNSs (mat and pat CNSs), the difference between the two alleles in H3K27me3 levels was only observed in pat CNSs. Together, our findings point to the importance of both transposons and CNSs in genomic imprinting in maize.

GENESPACE tracks regions of interest and gene copy number variation across multiple genomes

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence–absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to existing gene family and synteny programs, especially in polyploid, outbred, and other complex genomes.

The genome is the complete DNA sequence of an individual. It is a crucial foundation for many studies in medicine, agriculture, and conservation biology. Advances in genetics have made it possible to rapidly sequence, or read out, the genome of many organisms. For closely related species, scientists can then do detailed comparisons, revealing similar genes with a shared past or a common role, but comparing more distantly related organisms remains difficult.

One major challenge is that genes are often lost or duplicated over evolutionary time. One way to be more confident is to look at ‘synteny’, or how genes are organized or ordered within the genome. In some groups of species, synteny persists across millions of years of evolution. Combining sequence similarity with gene order could make comparisons between distantly related species more robust.

To do this, Lovell et al. developed GENESPACE, a software that links similarities between DNA sequences to the order of genes in a genome. This allows researchers to visualize and explore related DNA sequences and determine whether genes have been lost or duplicated. To demonstrate the value of GENESPACE, Lovell et al. explored evolution in vertebrates and flowering plants. The software was able to highlight the shared sequences between unique sex chromosomes in birds and mammals, and it was able to track the positions of genes important in the evolution of grass crops including maize, wheat, and rice.

Exploring the genetic code in this way could lead to a better understanding of the evolution of important sections of the genome. It might also allow scientists to find target genes for applications like crop improvement. Lovell et al. have designed the GENESPACE software to be easy for other scientists to use, allowing them to make graphics and perform analyses with few programming skills.

Genome-Wide Characterization of Nitrogenase Reductase (nifH) Genes in the Sweet Potato [Ipomoea batatas (L.) Lam] and Its Wild Ancestors

The sweet potato (Ipomoea batatas (L.) Lam.) is an important and widely grown crop, and the nitrogenase reductase (nifH) gene is the most widely sequenced marker gene used to identify nitrogen-fixing bacteria and archaea. There have been many examples of the isolation of the diazotrophic endophytes in sweet potatoes, and there has been no report on whether sweet potatoes and their wild ancestors harbored nifH genes. In this study, a comprehensive analysis of nifH genes has been conducted on these species by using bioinformatics and molecular biology methods. A total of 20, 19 and 17 nifH genes were identified for the first time in sweet potatoes, I. trifida and I. triloba, respectively. Based on a phylogenetic analysis, all of the nifH genes, except for g10233.t1, itf14g14040.t1 and itb14g15470.t1, were clustered into five independent clades: I, II, III, IV and V. The nifH genes clustered in the same phylogenetic branch showed a more similar distribution of conserved motifs and exons–introns than those of the other ones. All of the identified genes were further mapped on the 15 chromosomes of the sweet potato, I. trifida and I. triloba. No segmental duplication was detected in each genome of three Ipomoea species, and 0, 8 and 7 tandemly duplicated gene pairs were detected in the genome of the sweet potato, I. trifida and I. triloba, respectively. Synteny analysis between the three Ipomoea species revealed that there were 7, 7 and 8 syntenic gene pairs of nifH genes detected between the sweet potato and I. trifida, between the sweet potato and I. triloba and between I. trifida and I. triloba, respectively. All of the duplicated and syntenic nifH genes were subjected to purifying selection inside duplicated genomic elements during speciation, except for the tandemly duplicated gene pair itf11g07340.t2_itf11g07340.t3, which was subjected to positive selection. Different expression profiles were detected in the sweet potato, I. trifida and I. triloba. According to the above results, four nifH genes of the sweet potato (g950, g16683, g27094 and g33987) were selected for quantitative real-time polymerase chain reaction (qRT-PCR) analysis in two sweet potato cultivars (Eshu 15 and Long 9) under nitrogen deficiency (N0) and normal (N1) conditions. All of them were upregulated in the N1 treatment and were consistent with the analysis of the RNA-seq data. We hope that these results will provide new insights into the nifH genes in the sweet potato and its wild ancestors and will contribute to the molecular breeding of sweet potatoes in the future.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

Convergent evolution of polyploid genomes from across the eukaryotic tree of life

By modeling the homoeologous gene losses that occurred in 50 genomes deriving from ten distinct polyploidy events, we show that the evolutionary forces acting on polyploids are remarkably similar, regardless of whether they occur in flowering plants, ciliates, fishes, or yeasts. We show that many of the events show a relative rate of duplicate gene loss before the first postpolyploidy speciation that is significantly higher than in later phases of their evolution. The relatively weak selective constraint experienced by the single-copy genes these losses produced leads us to suggest that most of the purely selectively neutral duplicate gene losses occur in the immediate postpolyploid period. Nearly all of the events show strong evidence of biases in the duplicate losses, consistent with them being allopolyploidies, with 2 distinct progenitors contributing to the modern species. We also find ongoing and extensive reciprocal gene losses (alternative losses of duplicated ancestral genes) between these genomes. With the exception of a handful of closely related taxa, all of these polyploid organisms are separated from each other by tens to thousands of reciprocal gene losses. As a result, it is very unlikely that viable diploid hybrid species could form between these taxa, since matings between such hybrids would tend to produce offspring lacking essential genes. It is, therefore, possible that the relatively high frequency of recurrent polyploidies in some lineages may be due to the ability of new polyploidies to bypass reciprocal gene loss barriers.

Genome-wide identification and characterization of NBS-encoding genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.)

The most predominant type of resistance (R) genes contain nucleotide-binding sites and leucine-rich repeat (NBS-LRR) domains, characterization of which is helpful for plant resistance improvement. However, the NBS genes of Ipomoea trifida (H.B.K.) remain insufficient to date. In this study, a genome-wide analysis of the NBS-encoding gene in I. trifida (H.B.K.) was carried out. A total of 442 NBS encoding genes were identified, amounting to 1.37% of the total genes of I. trifida (H.B.K.). Based on the analysis of the domains, the identified ItfNBS genes were further classified into seven groups: CNL, NL, CN, N, TNL, TN, and RNL. Phylogenetic analysis showed that the I. trifida NBS genes clustered into three independent clades: RNL, TNL, and CNL. Chromosome location analysis revealed that the distribution of ItfNBS genes in chromosomes was uneven, with a number ranging from 3 to 45. Multiple stress-related regulatory elements were detected in the promoters of the NBS-encoding genes, and their expression profiles were obtained. The qRT-PCR analysis revealed that IbNBS10, IbNBS20, IbNBS258, and IbNBS88 responded to stem nematode infection. These results provide critical proof for further characterization and analysis of NBS-encoding genes with important functions.

High-quality chromosome-scale de novo assembly of the Paspalum notatum 'Flugge' genome

BACKGROUND

Paspalum notatum 'Flugge' is a diploid with 20 chromosomes (2n = 20) multi-purpose subtropical herb native to South America and has a high ecological significance. It is currently widely planted in tropical and subtropical regions. Despite the gene pool of P. notatum 'Flugge' being unearthed to a large extent in the past decade, no details about the genomic information of relevant species in Paspalum have been reported. In this study, the complete genome information of P. notatum was established and annotated through sequencing and de novo assembly of its genome.

RESULTS

The latest PacBio third-generation HiFi assembly and sequencing revealed that the genome size of P. notatum 'Flugge' is 541 M. The assembly result is the higher index among the genomes of the gramineous family published so far, with a contig N50 = 52Mbp, scaffold N50 = 49Mbp, and BUSCOs = 98.1%, accounting for 98.5% of the estimated genome. Genome annotation revealed 36,511 high-confidence gene models, thus providing an important resource for future molecular breeding and evolutionary research. A comparison of the genome annotation results of P. notatum 'Flugge' with other closely related species revealed that it had a close relationship with Zea mays but not close compared to Brachypodium distachyon, Setaria viridis, Oryza sativa, Puccinellia tenuiflora, Echinochloa crusgalli. An analysis of the expansion and contraction of gene families suggested that P. notatum 'Flugge' contains gene families associated with environmental resistance, increased reproductive ability, and molecular evolution, which explained its excellent agronomic traits.

CONCLUSION

This study is the first to report the high-quality chromosome-scale-based genome of P. notatum 'Flugge' assembled using the latest PacBio third-generation HiFi sequencing reads. The study provides an excellent genetic resource bank for gramineous crops and invaluable perspectives regarding the evolution of gramineous plants.

Peripheral structures in unlabelled trees and the accumulation of subgenomes in the evolution of polyploids

Many angiosperms have undergone some series of polyploidization events over the course of their evolutionary history. In these genomes, especially those resulting from multiple autopolyploidization, it may be relatively easy to recognize all the ξ sets of n homeologous chromosomes, but it is much harder, if not impossible, to partition these chromosomes into n subgenomes, each representing one distinct genomic component of ξ chromosomes making up the original polyploid. Thus, if we wish to infer the polyploidization history of the genome, we could make use of all the gene trees inferred from the genes in one set of homeologous chromosomes to construct a consensus tree, but there is no evident way of combining the trees from the ξ different sets, because we have no labelling of the chromosomes that is known to be consistent across these sets. We suggest here that lacking a consistent leaf-labelling, the topological structure of the trees may display sufficient resemblance so that a higher level consensus could be revealing of evolutionary history. This would be especially true of the peripheral structures of the tree, likely representing events that occurred more recently and have thus been less obscured by subsequent evolutionary processes. Here, we present a statistical test to assess whether the subgenomes in a polyploid genome could have been added one at a time. The null hypothesis is that the accumulation of chromosomes follows a stochastic process in which transition from one generation to the next is through randomly choosing an edge, and then subdividing this edge in order to link the new internal vertex to a new external vertex. We analyze the probability distributions of a number of peripheral tree substructures, namely leaf- or terminal-pairs, triples and quadruples, arising from this stochastic process, in terms of some exact recurrences. We propose some conjectures regarding the asymptotic behaviours of these distributions. Applying our analysis to a sugarcane genome, we demonstrate that it is unlikely that the accumulation of subgenomes has occurred one at a time.

Genomic and Meiotic Changes Accompanying Polyploidization

Hybridization and polyploidy have been considered as significant evolutionary forces in adaptation and speciation, especially among plants. Interspecific gene flow generates novel genetic variants adaptable to different environments, but it is also a gene introgression mechanism in crops to increase their agronomical yield. An estimate of 9% of interspecific hybridization has been reported although the frequency varies among taxa. Homoploid hybrid speciation is rare compared to allopolyploidy. Chromosome doubling after hybridization is the result of cellular defects produced mainly during meiosis. Unreduced gametes, which are formed at an average frequency of 2.52% across species, are the result of altered spindle organization or orientation, disturbed kinetochore functioning, abnormal cytokinesis, or loss of any meiotic division. Meiotic changes and their genetic basis, leading to the cytological diploidization of allopolyploids, are just beginning to be understood especially in wheat. However, the nature and mode of action of homoeologous recombination suppressor genes are poorly understood in other allopolyploids. The merger of two independent genomes causes a deep modification of their architecture, gene expression, and molecular interactions leading to the phenotype. We provide an overview of genomic changes and transcriptomic modifications that particularly occur at the early stages of allopolyploid formation.

Characterization of Nucleotide Binding -Site-Encoding Genes in Sweetpotato, Ipomoea batatas(L.) Lam., and Their Response to Biotic and Abiotic Stresses

Sweetpotato, Ipomoea batatas (L.) Lam., is an important and widely grown crop, yet its production is affected severely by biotic and abiotic stresses. The nucleotide binding site (NBS)-encoding genes have been shown to improve stress tolerance in several plant species. However, the characterization of NBS-encoding genes in sweetpotato is not well-documented to date. In this study, a comprehensive analysis of NBS-encoding genes has been conducted on this species by using bioinformatics and molecular biology methods. A total of 315 NBS-encoding genes were identified, and 260 of them contained all essential conserved domains while 55 genes were truncated. Based on domain architectures, the 260 NBS-encoding genes were grouped into 6 distinct categories. Phylogenetic analysis grouped these genes into 3 classes: TIR, CC (I), and CC (II). Chromosome location analysis revealed that the distribution of NBS-encoding genes in chromosomes was uneven, with a number ranging from 1 to 34. Multiple stress-related regulatory elements were detected in the promoters, and the NBS-encoding genes' expression profiles under biotic and abiotic stresses were obtained. According to the bioinformatics analysis, 9 genes were selected for RT-qPCR analysis. The results revealed that IbNBS75, IbNBS219, and IbNBS256 respond to stem nematode infection; Ib-NBS240, IbNBS90, and IbNBS80 respond to cold stress, while IbNBS208, IbNBS71, and IbNBS159 respond to 30% PEG treatment. We hope these results will provide new insights into the evolution of NBS-encoding genes in the sweetpotato genome and contribute to the molecular breeding of sweetpotato in the future.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

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.

Differential Regulation of Maize and Sorghum Orthologs in Response to the Fungal Pathogen Exserohilum turcicum

Pathogens that infect more than one host offer an opportunity to study how resistance mechanisms have evolved across different species. Exserohilum turcicum infects both maize and sorghum and the isolates are host-specific, offering a unique system to examine both compatible and incompatible interactions. We conducted transcriptional analysis of maize and sorghum in response to maize-specific and sorghum-specific E. turcicum isolates and identified functionally related co-expressed modules. Maize had a more robust transcriptional response than sorghum. E. turcicum responsive genes were enriched in core orthologs in both crops, but only up to 16% of core orthologs showed conserved expression patterns. Most changes in gene expression for the core orthologs, including hub genes, were lineage specific, suggesting a role for regulatory divergent evolution. We identified several defense-related shared differentially expressed (DE) orthologs with conserved expression patterns between the two crops, suggesting a role for parallel evolution of those genes in both crops. Many of the differentially expressed genes (DEGs) during the incompatible interaction were related to quantitative disease resistance (QDR). This work offers insights into how different hosts with relatively recent divergence interact with a common pathogen. Our results are important for developing resistance to this critical pathogen and understanding the evolution of host-pathogen interactions.

The genome of the warm-season turfgrass African bermudagrass (Cynodon transvaalensis)

Cynodon species can be used for multiple purposes and have high economic and ecological significance. However, the genetic basis of the favorable agronomic traits of Cynodon species is poorly understood, partially due to the limited availability of genomic resources. In this study, we report a chromosome-scale genome assembly of a diploid Cynodon species, C. transvaalensis, obtained by combining Illumina and Nanopore sequencing, BioNano, and Hi-C. The assembly contains 282 scaffolds (~423.42 Mb, N50 = 5.37 Mb), which cover ~93.2% of the estimated genome of C. transvaalensis (~454.4 Mb). Furthermore, 90.48% of the scaffolds (~383.08 Mb) were anchored to nine pseudomolecules, of which the largest was 60.78 Mb in length. Evolutionary analysis along with transcriptome comparison provided a preliminary genomic basis for the adaptation of this species to tropical and/or subtropical climates, typically with dry summers. The genomic resources generated in this study will not only facilitate evolutionary studies of the Chloridoideae subfamily, in particular, the Cynodonteae tribe, but also facilitate functional genomic research and genetic breeding in Cynodon species for new leading turfgrass cultivars in the future.

Structural variation at the maize WUSCHEL1 locus alters stem cell organization in inflorescences

Structural variation in plant genomes is a significant driver of phenotypic variability in traits important for the domestication and productivity of crop species. Among these are traits that depend on functional meristems, populations of stem cells maintained by the CLAVATA-WUSCHEL (CLV-WUS) negative feedback-loop that controls the expression of the WUS homeobox transcription factor. WUS function and impact on maize development and yield remain largely unexplored. Here we show that the maize dominant Barren inflorescence3 (Bif3) mutant harbors a tandem duplicated copy of the ZmWUS1 gene, ZmWUS1-B, whose novel promoter enhances transcription in a ring-like pattern. Overexpression of ZmWUS1-B is due to multimerized binding sites for type-B RESPONSE REGULATORs (RRs), key transcription factors in cytokinin signaling. Hypersensitivity to cytokinin causes stem cell overproliferation and major rearrangements of Bif3 inflorescence meristems, leading to the formation of ball-shaped ears and severely affecting productivity. These findings establish ZmWUS1 as an essential meristem size regulator in maize and highlight the striking effect of cis-regulatory variation on a key developmental program.

The WUSCHEL transcription factor promotes plant stem cell proliferation. Here the authors show that the maize Bif3 mutant contains a duplication of the ZmWUS1 locus leading to cytokinin hypersensitivity and overproliferation at the shoot meristem demonstrating the role of WUSCHEL in maize and how structural variation can impact plant morphology.

Nucleoside Metabolism Is Induced in Common Bean During Early Seedling Development

Nucleoside hydrolases (NSH; nucleosidases) catalyze the cleavage of nucleosides into ribose and free nucleobases. These enzymes have been postulated as key elements controlling the ratio between nucleotide salvage and degradation. Moreover, they play a pivotal role in ureidic legumes by providing the substrate for the synthesis of ureides. Furthermore, nucleotide metabolism has a crucial role during germination and early seedling development, since the developing seedlings require high amount of nucleotide simultaneously to the mobilization of nutrient in cotyledons. In this study, we have cloned two nucleosidases genes from Phaseolus vulgaris, PvNSH1 and PvNSH2, expressed them as recombinant proteins, and characterized their catalytic activities. Both enzymes showed a broad range of substrate affinity; however, PvNSH1 exhibited the highest activity with uridine, followed by xanthosine, whereas PvNSH2 hydrolyses preferentially xanthosine and shows low activity with uridine. The study of the regulation of nucleosidases during germination and early postgerminative development indicated that nucleosidases are induced in cotyledons and embryonic axes just after the radicle emergence, coincident with the induction of nucleases activity and the synthesis of ureides in the embryonic axes, with no remarkable differences in the level of expression of both nucleosidase genes. In addition, nucleosides and nucleobase levels were determined as well in cotyledons and embryonic axes. Our results suggest that PvNSH1 and PvNSH2 play an important role in the mobilization of nutrients during this crucial stage of plant development.

An updated census of the maize TIFY family

The TIFY gene family is a plant-specific gene family encoding a group of proteins characterized by its namesake, the conservative TIFY domain and members can be organized into four subfamilies: ZML, TIFY, PPD and JAZ (Jasmonate ZIM-domain protein) by presence of additional conserved domains. The TIFY gene family is intensively explored in several model and agriculturally important crop species and here, yet the composition of the TIFY family of maize has remained unresolved. This study increases the number of maize TIFY family members known by 40%, bringing the total to 47 including 38 JAZ, 5 TIFY, and 4 ZML genes. The majority of the newly identified genes were belonging to the JAZ subfamily, six of which had aberrant TIFY domains, suggesting loss JAZ-JAZ or JAZ-NINJA interactions. Six JAZ genes were found to have truncated Jas domain or an altered degron motif, suggesting resistance to classical JAZ degradation. In addition, seven membranes were found to have an LxLxL-type EAR motif which allows them to recruit TPL/TPP co-repressors directly without association to NINJA. Expression analysis revealed that ZmJAZ14 was specifically expressed in the seeds and ZmJAZ19 and 22 in the anthers, while the majority of other ZmJAZs were generally highly expressed across diverse tissue types. Additionally, ZmJAZ genes were highly responsive to wounding and JA treatment. This study provides a comprehensive update of the maize TIFY/JAZ gene family paving the way for functional, physiological, and ecological analysis.

Crop breeding - From experience-based selection to precision design

Crops are the foundation of human society, not only by providing needed nutrition, but also by feeding livestock and serving as raw materials for industry. Cereal crops, which supply most of our calories, have been supporting humans for thousands of years. However food security is facing many challenges nowadays, including growing populations, water shortage, and increased incidence of biotic and abiotic stresses. According to statistical data from the Food and Agriculture Organization of the United Nations (FAO, http://www.fao.org/), the people suffering severe food insecurity increased from 7.9 % in 2015 to 9.7 % in 2019 and the number of people exposed to moderate or severe food insecurity have increased by 400 million over the same time period. Although there are many ways to cope with these challenges, crop breeding remains the most crucial and direct manner. With the development of molecular genetics, the speed of cloning genetic variations underlying corresponding phenotypes of agricultural importance is considerably more rapid. As a consequence breeding methods have evolved from phenotype-based to genome-based selection. In the future, knowledge-driven crop design, which integrates multi-omics data to reveal the connections between genotypes and phenotypes and to build selection models, will undoubtedly become the most efficient way to shape plants, to improve crops, and to ensure food security.

Unbiased Subgenome Evolution in Allotetraploid Species of Ephedra and Its Implications for the Evolution of Large Genomes in Gymnosperms

The evolutionary dynamics of polyploid genomes and consequences of polyploidy have been studied extensively in angiosperms but very rarely in gymnosperms. The gymnospermous genus Ephedra is characterized by a high frequency of polyploidy, and thus provides an ideal system to investigate the evolutionary mode of allopolyploid genomes and test whether subgenome dominance has occurred in gymnosperms. Here, we sequenced transcriptomes of two allotetraploid species of Ephedra and their putative diploid progenitors, identified expressed homeologs, and analyzed alternative splicing and homeolog expression based on PacBio Iso-Seq and Illumina RNA-seq data. We found that the two subgenomes of the allotetraploids had similar numbers of expressed homeologs, similar percentages of homeologs with dominant expression, and approximately equal numbers of isoforms with alternative splicing, showing an unbiased subgenome evolution as in a few polyploid angiosperms, with a divergence of the two subgenomes at ∼8 Ma. In addition, the nuclear DNA content of the allotetraploid species is almost equal to the sum of two putative progenitors, suggesting limited genome restructuring after allotetraploid speciation. The allopolyploid species of Ephedra might have undergone slow diploidization, and the unbiased subgenome evolution implies that the formation of large genomes in gymnosperms could be attributed to even and slow fractionation following polyploidization.

Profiling Alternative 3' Untranslated Regions in Sorghum using RNA-seq Data

Sorghum is an important crop widely used for food, feed, and fuel. Transcriptome-wide studies of 3′ untranslated regions (3′UTR) using regular RNA-seq remain scarce in sorghum, while transcriptomes have been characterized extensively using Illumina short-read sequencing platforms for many sorghum varieties under various conditions or developmental contexts. 3′UTR is a critical regulatory component of genes, controlling the translation, transport, and stability of messenger RNAs. In the present study, we profiled the alternative 3′UTRs at the transcriptome level in three genetically related but phenotypically contrasting lines of sorghum: Rio, BTx406, and R9188. A total of 1,197 transcripts with alternative 3′UTRs were detected using RNA-seq data. Their categorization identified 612 high-confidence alternative 3′UTRs. Importantly, the high-confidence alternative 3′UTR genes significantly overlapped with the genesets that are associated with RNA N⁶-methyladenosine (m⁶A) modification, suggesting a clear indication between alternative 3′UTR and m⁶A methylation in sorghum. Moreover, taking advantage of sorghum genetics, we provided evidence of genotype specificity of alternative 3′UTR usage. In summary, our work exemplifies a transcriptome-wide profiling of alternative 3′UTRs using regular RNA-seq data in non-model crops and gains insights into alternative 3′UTRs and their genotype specificity.

Evolution and Domestication Footprints Uncovered from the Genomes of Coix

Coix lacryma-jobi, a plant species closely related to Zea and Sorghum, is an important food and medicinal crop in Asia. However, no reference genome of this species has been reported, and its exact phylogeny within the Andropogoneae remains unresolved. Here, we generated a high-quality genome assembly of coix comprising ∼1.73 Gb with 44 485 predicted protein-coding genes. We found coix to be a typical diploid plant with an overall 1-to-1 syntenic relationship with the Sorghum genome, despite its drastic genome expansion (∼2.3-fold) due mainly to the activity of transposable elements. Phylogenetic analysis revealed that coix diverged with sorghum ∼10.41 million years ago, which was ∼1.49 million years later than the divergence between sorghum and maize. Resequencing of 27 additional coix accessions revealed that they could be unambiguously separated into wild relatives and cultivars, and suggested that coix experienced a strong genetic bottleneck, resulting in the loss of about half of the genetic diversity during domestication, even though many traits have remained undomesticated. Our data not only provide novel comparative genomic and evolutionary insights into the Andropogoneae lineage, but also an important resource that will greatly benefit molecular breeding of this important crop.

Identification, characterization, and functional prediction of circular RNAs in maize

Circular RNAs (circRNAs) are a new type of intracellular regulator that have been widely identified in animals and plants by high-throughput sequencing. However, there are still few functional studies on circRNAs in plants. To better understand maize circRNAs and their potential functions, we identified 1199 circRNAs in maize from RiboMinus RNA-Seq transcriptome data, and found distinct features of splicing site selection bias, longer flanking introns, and miniature inverted-repeat transposable element (MITE) insertions in flanking introns in maize circRNAs compared to other plant circRNAs. In total, 31 and 36 orthologous circRNAs were identified in rice and maize, respectively, but the orthologous parental genes could not produce orthologous circRNAs, mostly because of long-sequence insertions/deletions at flanking introns and approximately 24.3% of them contained MITE sequences. The majority of maize circRNAs showed high diversity of expression under different treatments and/or in different genetic backgrounds, implying that circRNAs could be involved in various regulatory networks. Twenty-six ecircRNAs were predicted to contain one or more target mimics, and 229 circRNAs had high coding potential, indicating that circRNAs could perform peptide-encoding functions in plants. These results will broaden understanding of the roles of circRNAs in plants and support further functional work on maize.

A two-step mutation process in the double WS1 homologs drives the evolution of burley tobacco, a special chlorophyll-deficient mutant with abnormal chloroplast development

MAIN CONCLUSION

The functional homologs WS1A and WS1B, identified by map-based cloning, control the burley character by affecting chloroplast development in tobacco, contributing to gene isolation and genetic improvement in polyploid crops. Burley represents a special type of tobacco (Nicotiana tabacum L.) cultivar that is characterized by a white stem with a high degree of chlorophyll deficiency. Although important progress in the research of burley tobacco has been made, the molecular mechanisms underlying this character remain unclear. Here, on the basis of our previous genetic analyses and preliminary mapping results, we isolated the White Stem 1A (WS1A) and WS1B genes using a map-based cloning approach. WS1A and WS1B are functional homologs with completely identical biological functions and highly similar expression patterns that control the burley character in tobacco. WS1A and WS1B are derived from Nicotiana sylvestris and Nicotiana tomentosiformis, the diploid ancestors of Nicotiana tabacum, respectively. The two genes encode zinc metalloproteases of the M50 family that are highly homologous to the Ethylene-dependent Gravitropism-deficient and Yellow-green 1 (EGY1) protein of Arabidopsis and the Lutescent 2 (L2) protein of tomato. Transmission electron microscopic examinations indicated that WS1A and WS1B are involved in the development of chloroplasts by controlling the formation of thylakoid membranes, very similar to that observed for EGY1 and L2. The genotyping of historical tobacco varieties revealed that a two-step mutation process occurred in WS1A and WS1B during the evolution of burley tobacco. We also discussed the strategy for gene map-based cloning in polyploid plants with complex genomes. This study will facilitate the identification of agronomically important genes in tobacco and other polyploid crops and provide insights into crop improvement via molecular approaches.

Genome-Wide Analysis of TCP Family Genes in Zea mays L. Identified a Role for ZmTCP42 in Drought Tolerance

The Teosinte-branched 1/Cycloidea/Proliferating (TCP) plant-specific transcription factors (TFs) have been demonstrated to play a fundamental role in plant development and organ patterning. However, it remains unknown whether or not the TCP gene family plays a role in conferring a tolerance to drought stress in maize, which is a major constraint to maize production. In this study, we identified 46 ZmTCP genes in the maize genome and systematically analyzed their phylogenetic relationships and synteny with rice, sorghum, and ArabidopsisTCP genes. Expression analysis of the 46 ZmTCP genes in different tissues and under drought conditions, suggests their involvement in maize response to drought stress. Importantly, genetic variations in ZmTCP32 and ZmTCP42 are significantly associated with drought tolerance at the seedling stage. RT-qPCR results suggest that ZmTCP32 and ZmTCP42 RNA levels are both induced by ABA, drought, and polyethylene glycol treatments. Based on the significant association between the genetic variation of ZmTCP42 and drought tolerance, and the inducible expression of ZmTCP42 by drought stress, we selected ZmTCP42, to investigate its function in drought response. We found that overexpression of ZmTCP42 in Arabidopsis led to a hypersensitivity to ABA in seed germination and enhanced drought tolerance, validating its function in drought tolerance. These results suggested that ZmTCP42 functions as an important TCP TF in maize, which plays a positive role in drought tolerance.

Glycome and Proteome Components of Golgi Membranes Are Common between Two Angiosperms with Distinct Cell-Wall Structures

The plant endoplasmic reticulum-Golgi apparatus is the site of synthesis, assembly, and trafficking of all noncellulosic polysaccharides, proteoglycans, and proteins destined for the cell wall. As grass species make cell walls distinct from those of dicots and noncommelinid monocots, it has been assumed that the differences in cell-wall composition stem from differences in biosynthetic capacities of their respective Golgi. However, immunosorbence-based screens and carbohydrate linkage analysis of polysaccharides in Golgi membranes, enriched by flotation centrifugation from etiolated coleoptiles of maize (Zea mays) and leaves of Arabidopsis (Arabidopsis thaliana), showed that arabinogalactan-proteins and arabinans represent substantial portions of the Golgi-resident polysaccharides not typically found in high abundance in cell walls of either species. Further, hemicelluloses accumulated in Golgi at levels that contrasted with those found in their respective cell walls, with xyloglucans enriched in maize Golgi, and xylans enriched in Arabidopsis. Consistent with this finding, maize Golgi membranes isolated by flotation centrifugation and enriched further by free-flow electrophoresis, yielded >200 proteins known to function in the biosynthesis and metabolism of cell-wall polysaccharides common to all angiosperms, and not just those specific to cell-wall type. We propose that the distinctive compositions of grass primary cell walls compared with other angiosperms result from differential gating or metabolism of secreted polysaccharides post-Golgi by an as-yet unknown mechanism, and not necessarily by differential expression of genes encoding specific synthase complexes.

Evolution of the Cell Wall Gene Families of Grasses

Grasses and related commelinid monocot species synthesize cell walls distinct in composition from other angiosperm species. With few exceptions, the genomes of all angiosperms contain the genes that encode the enzymes for synthesis of all cell-wall polysaccharide, phenylpropanoid, and protein constituents known in vascular plants. RNA-seq analysis of transcripts expressed during development of the upper and lower internodes of maize (Zea mays) stem captured the expression of cell-wall-related genes associated with primary or secondary wall formation. High levels of transcript abundances were not confined to genes associated with the distinct walls of grasses but also of those associated with xyloglucan and pectin synthesis. Combined with proteomics data to confirm that expressed genes are translated, we propose that the distinctive cell-wall composition of grasses results from sorting downstream from their sites of synthesis in the Golgi apparatus and hydrolysis of the uncharacteristic polysaccharides and not from differential expression of synthases of grass-specific polysaccharides.

Widespread ancient whole-genome duplications in Malpighiales coincide with Eocene global climatic upheaval

Whole-genome duplications (WGDs) are widespread and prevalent in vascular plants and frequently coincide with major episodes of global and climatic upheaval, including the mass extinction at the Cretaceous-Tertiary boundary (c. 65 Ma) and during more recent periods of global aridification in the Miocene (c. 10-5 Ma). Here, we explore WGDs in the diverse flowering plant clade Malpighiales. Using transcriptomes and complete genomes from 42 species, we applied a multipronged phylogenomic pipeline to identify, locate, and determine the age of WGDs in Malpighiales using three means of inference: distributions of synonymous substitutions per synonymous site (Ks ) among paralogs, phylogenomic (gene tree) reconciliation, and a likelihood-based gene-count method. We conservatively identify 22 ancient WGDs, widely distributed across Malpighiales subclades. Importantly, these events are clustered around the Eocene-Paleocene transition (c. 54 Ma), during which time the planet was warmer and wetter than any period in the Cenozoic. These results establish that the Eocene Climatic Optimum likely represents a previously unrecognized period of prolific WGDs in plants, and lends further support to the hypothesis that polyploidization promotes adaptation and enhances plant survival during episodes of global change, especially for tropical organisms like Malpighiales, which have tight thermal tolerances.

A Malvaceae mystery: A mallow maelstrom of genome multiplications and maybe misleading methods?

Previous research suggests that Gossypium has undergone a 5- to 6-fold multiplication following its divergence from Theobroma. However, the number of events, or where they occurred in the Malvaceae phylogeny remains unknown. We analyzed transcriptomic and genomic data from representatives of eight of the nine Malvaceae subfamilies. Phylogenetic analysis of nuclear data placed Dombeya (Dombeyoideae) as sister to the rest of Malvadendrina clade, but the plastid DNA tree strongly supported Durio (Helicteroideae) in this position. Intraspecific Ks plots indicated that all sampled taxa, except Theobroma (Byttnerioideae), Corchorus (Grewioideae), and Dombeya (Dombeyoideae), have experienced whole genome multiplications (WGMs). Quartet analysis suggested WGMs were shared by Malvoideae-Bombacoideae and Sterculioideae-Tilioideae, but did not resolve whether these are shared with each other or Helicteroideae (Durio). Gene tree reconciliation and Bayesian concordance analysis suggested a complex history. Alternative hypotheses are suggested, each involving two independent autotetraploid and one allopolyploid event. They differ in that one entails an allopolyploid origin for the Durio lineage, whereas the other invokes an allopolyploid origin for Malvoideae-Bombacoideae. We highlight the need for more genomic information in the Malvaceae and improved methods to resolve complex evolutionary histories that may include allopolyploidy, incomplete lineage sorting, and variable rates of gene and genome evolution.

Transcriptome and miRNAs analyses enhance our understanding of the evolutionary advantages of polyploidy

Polyploid organisms have more than two sets of chromosomes, including autopolyploid via intraspecific genome doubling, and allopolyploid via merging genomes of distinct species by hybridization. Polyploid organisms are widespread in plants, indicating that polyploidy has some evolutionary advantages over its diploid ancestor. Actually, polyploidy is always tightly associated with hybrid vigor and adaptation to adverse environmental conditions. However, why polyploidy can develop such advantages is poorly known. MicroRNAs (miRNAs) are endogenous ∼21 nt small RNAs which can play important regulatory roles in animals and plants by targeting mRNAs for cleavage or translational repression. MicroRNAs are essential for cell development, differentiation, signal transduction, and show an adaptive response to biotic and abiotic stresses. Environmental stresses cause plants to over- or under-express certain miRNAs or synthesize new miRNAs to cope with stress. We have here reviewed our current knowledge on the molecular mechanisms, which can account for the evolutionary advantages of polyploidy over its diploid ancestor from genome-wide gene expression and microRNAs expression perspectives.

An AGAMOUS-like factor is associated with the origin of two domesticated varieties in Cymbidium sinense (Orchidaceae)

Cymbidium has been artificially domesticated for centuries in Asia, which produced numerous cultivated varieties. Flowers with stamenoid tepals or those with multiple tepals have been found in different species of Cymbidium; however, the molecular basis controlling the formation of these phenotypes is still largely unknown. Previous work demonstrated that AGAMOUS/AG lineage MADS genes function in floral meristem determinacy as well as in reproductive organs development in both dicots and monocots, indicating a possible relationship with the origin of two flower varieties in Cymbidium. Here, we characterized and analyzed two AG lineage paralogues, CsAG1 and CsAG2, from Cymbidium sinense, both of which were highly expressed in the gynostemium column of a standard C. sinense. Interestingly, we detected ectopic expression of CsAG1 rather than CsAG2 in all floral organs of a stamenoid-tepal variety and significant down-regulation of CsAG1 in a variety with multiple tepals. Over-expression of CsAG1 in wild type Arabidopsis resulted in petal-to-stamen homeotic conversion, suggesting a conserved C-function of CsAG1 in the development of Cymbidium flower. Altogether, our results supported a hypothesis that disruption of a single AG-like factor would be associated with the formation of two domesticated varieties in C. sinense.

Mysterious meiotic behavior of autopolyploid and allopolyploid maize

This study was aimed to investigate the stability of chromosomes during meiosis in autopolyploid and allopolyploid maize, as well as to determine an association of chromosomes between maize (Zea mays ssp. mays Linnaeus, 1753) and Z. perennis (Hitchcock, 1922) Reeves & Mangelsdor, 1942, by producing a series of autopolyploid and allopolyploid maize hybrids. The intra-genomic and inter-genomic meiotic pairings in these polyploids were quantified and compared using dual-color genomic in-situ hybridization. The results demonstrated higher level of chromosome stability in allopolyploid maize during meiosis as compared to autopolyploid maize. In addition, the meiotic behavior of Z. perennis was relatively more stable as compared to the allopolyploid maize. Moreover, ten chromosomes of "A" subgenome in maize were homologous to twenty chromosomes of Z. perennis genome with a higher pairing frequency and little evolutionary differentiation. At the same time, little evolutionary differentiation has been shown by chromosomes of "A" subgenome in maize, while chromosomes of "B" subgenome, had a lower pairing frequency and higher evolutionary differentiation. Furthermore, 5IM + 5IIPP + 5IIIMPP and 5IIMM + 5IIPP + 5IVMMPP were observed in allotriploids and allotetraploids respectively, whereas homoeologous chromosomes were found between the "A" and "B" genome of maize and Z. perennis.

Allele phasing is critical to revealing a shared allopolyploid origin of Medicago arborea and M. strasseri (Fabaceae)

BACKGROUND

Whole genome duplication plays a central role in plant evolution. There are two main classes of polyploid formation: autopolyploids which arise within one species by doubling of similar homologous genomes; in contrast, allopolyploidy (hybrid polyploidy) arise via hybridization and subsequent doubling of nonhomologous (homoeologous) genomes. The distinction between polyploid origins can be made using gene phylogenies, if alleles from each genome can be correctly retrieved. We examined whether two closely related tetraploid Mediterranean shrubs (Medicago arborea and M. strasseri) have an allopolyploid origin - a question that has remained unsolved despite substantial previous research. We sequenced and analyzed ten low-copy nuclear genes from these and related species, phasing all alleles. To test the efficacy of allele phasing on the ability to recover the evolutionary origin of polyploids, we compared these results to analyses using unphased sequences.

RESULTS

In eight of the gene trees the alleles inferred from the tetraploids formed two clades, in a non-sister relationship. Each of these clades was more closely related to alleles sampled from other species of Medicago, a pattern typical of allopolyploids. However, we also observed that alleles from one of the remaining genes formed two clades that were sister to one another, as is expected for autopolyploids. Trees inferred from unphased sequences were very different, with the tetraploids often placed in poorly supported and different positions compared to results obtained using phased alleles.

CONCLUSIONS

The complex phylogenetic history of M. arborea and M. strasseri is explained predominantly by shared allotetraploidy. We also observed that an increase in woodiness is correlated with polyploidy in this group of species and present a new possibility that woodiness could be a transgressive phenotype. Correctly phased homoeologues are likely to be critical for inferring the hybrid origin of allopolyploid species, when most genes retain more than one homoeologue. Ignoring homoeologous variation by merging the homoeologues can obscure the signal of hybrid polyploid origins and produce inaccurate results.

Cytological diploidization of paleopolyploid genus Zea: Divergence between homoeologous chromosomes or activity of pairing regulator genes?

Cytological diploidization process is different in autopolyploid and allopolyploid species. Colchicine applied at the onset of meiosis suppresses the effect of pairing regulator genes resulting multivalents formation in bivalent-forming species. Colchicine treated maizes (4x = 2n = 20, A_mA_mB_mB_m) showed up to 5IV, suggesting pairing between chromosomes from genomes homoeologous A_m and B_m. In untreated individuals of the alloautooctoploid Zea perennis (8x = 2n = 40, A_pA_pA_p´A_p´B_p1B_p1B_p2B_p2) the most frequent configuration was 5IV+10II (formed by A and B genomes, respectively). The colchicine treated Z. perennis show up to 10IV revealing higher affinity within genomes A and B, but any homology among them. These results suggest the presence of a paring regulator locus (PrZ) in maize and Z. perennis, whose expression is suppressed by colchicine. It could be postulated that in Z. perennis, PrZ would affect independently the genomes A and B, being relevant the threshold of homology, the fidelity of pairing in each genomes and the ploidy level. Cytological analysis of the treated hexaploid hybrids (6x = 2n = 30), with Z. perennis as a parental, strongly suggests that PrZ is less effective in only one doses. This conclusion was reinforced by the homoeologous pairing observed in untreated dihaploid maizes, which showed up to 5II. Meiotic behaviour of individuals treated with different doses of colchicine allowed to postulate that PrZ affect the homoeologous association by controlling entire genomes (A_m or B_m) rather than individual chromosomes. Based on cytological and statistical results it is possible to propose that the cytological diploidization in Zea species occurs by restriction of pairing between homoeologous chromosomes or by genetical divergence of the homoeologous chromosomes, as was observed in untreated Z. mays ssp. parviglumis. These are independent but complementary systems and could be acting jointly in the same nucleus.

The limited role of differential fractionation in genome content variation and function in maize (Zea mays L.) inbred lines

Maize is a diverse paleotetraploid species with considerable presence/absence variation and copy number variation. One mechanism through which presence/absence variation can arise is differential fractionation. Fractionation refers to the loss of duplicate gene pairs from one of the maize subgenomes during diploidization. Differential fractionation refers to non-shared gene loss events between individuals following a whole-genome duplication event. We investigated the prevalence of presence/absence variation resulting from differential fractionation in the syntenic portion of the genome using two whole-genome de novo assemblies of the inbred lines B73 and PH207. Between these two genomes, syntenic genes were highly conserved with less than 1% of syntenic genes being subject to differential fractionation. The few variably fractionated syntenic genes that were identified are unlikely to contribute to functional phenotypic variation, as there is a significant depletion of these genes in annotated gene sets. In further comparisons of 60 diverse inbred lines, non-syntenic genes were six times more likely to be variable than syntenic genes, suggesting that comparisons among additional genome assemblies are not likely to result in the discovery of large-scale presence/absence variation among syntenic genes.

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.

Genome Halving with an Outgroup

Some genomes are known to have incurred a genome doubling (tetraploidization) event in their evolutionary history, and this is reflected today in patterns of duplicated segments scattered throughout their chromosomes. These duplications may be used as data to “halve” the genome, i.e. to reconstruct the ancestral genome at the moment of tetraploidization, but the solution is often highly non-unique. To resolve this problem, we adapt the genome halving algorithm of El-Mabrouk and Sankoff to take account of an external reference genome. We apply this to reconstruct the tetraploid ancestor of maize, using either rice or sorghum as the reference.

Competitive Ability of Maize Pollen Grains Requires Paralogous Serine Threonine Protein Kinases STK1 and STK2

serine threonine kinase1 (stk1) and serine threonine kinase2 (stk2) are closely related maize paralogous genes predicted to encode serine/threonine protein kinases. Pollen mutated in stk1 or stk2 competes poorly with normal pollen, pointing to a defect in pollen tube germination or growth. Both genes are expressed in pollen, but not in most other tissues. In germination media, STK1 and STK2 fluorescent fusion proteins localize to the plasma membrane of the vegetative cell. RNA-seq experiments identified 534 differentially expressed genes in stk1 mutant pollen relative to wild type. Gene ontology (GO) molecular functional analysis uncovered several differentially expressed genes with putative ribosome initiation and elongation functions, suggesting that stk1 might affect ribosome function. Of the two paralogs, stk1 may play a more important role in pollen development than stk2, as stk2 mutations have a smaller pollen transmission effect. However, stk2 does act as an enhancer of stk1 because the double mutant combination is only infrequently pollen-transmitted in double heterozygotes. We conclude that the stk paralogs play an essential role in pollen development.

Maize YABBY Genes drooping leaf1 and drooping leaf2 Regulate Plant Architecture

Leaf architecture directly influences canopy structure, consequentially affecting yield. We discovered a maize (Zea mays) mutant with aberrant leaf architecture, which we named drooping leaf1 (drl1). Pleiotropic mutations in drl1 affect leaf length and width, leaf angle, and internode length and diameter. These phenotypes are enhanced by natural variation at the drl2 enhancer locus, including reduced expression of the drl2-Mo17 allele in the Mo17 inbred. A second drl2 allele, produced by transposon mutagenesis, interacted synergistically with drl1 mutants and reduced drl2 transcript levels. The drl genes are required for proper leaf patterning, development and cell proliferation of leaf support tissues, and for restricting auricle expansion at the midrib. The paralogous loci encode maize CRABS CLAW co-orthologs in the YABBY family of transcriptional regulators. The drl genes are coexpressed in incipient and emergent leaf primordia at the shoot apex, but not in the vegetative meristem or stem. Genome-wide association studies using maize NAM-RIL (nested association mapping-recombinant inbred line) populations indicated that the drl loci reside within quantitative trait locus regions for leaf angle, leaf width, and internode length and identified rare single nucleotide polymorphisms with large phenotypic effects for the latter two traits. This study demonstrates that drl genes control the development of key agronomic traits in maize.

Diverse genome organization following 13 independent mesopolyploid events in Brassicaceae contrasts with convergent patterns of gene retention

Hybridization and polyploidy followed by genome-wide diploidization had a significant impact on the diversification of land plants. The ancient At-α whole-genome duplication (WGD) preceded the diversification of crucifers (Brassicaceae). Some genera and tribes also experienced younger, mesopolyploid WGDs concealed by subsequent genome diploidization. Here we tested if multiple base chromosome numbers originated due to genome diploidization after independent mesopolyploid WGDs and how diploidization affected post-polyploid gene retention. Sixteen species representing 10 Brassicaceae tribes were analyzed by comparative chromosome painting and/or whole-transcriptome analysis of gene age distributions and phylogenetic analyses of gene duplications. Overall, we found evidence for at least 13 independent mesopolyploidies followed by different degrees of diploidization across the Brassicaceae. New mesotetraploid events were uncovered for the tribes Anastaticeae, Iberideae and Schizopetaleae, and mesohexaploid WGDs for Cochlearieae and Physarieae. In contrast, we found convergent patterns of gene retention and loss among these independent WGDs. Our combined analyses of genomic data for Brassicaceae indicate that extant chromosome number variation in many plant groups, and especially monophyletic taxa with multiple base chromosome numbers, can result from clade-specific genome duplications followed by diploidization. Our observation of parallel gene retention and loss across multiple independent WGDs provides one of the first multi-species tests of the predictability of patterns of post-polyploid genome evolution.

Natural variation for gene expression responses to abiotic stress in maize

Plants respond to abiotic stress through a variety of physiological, biochemical, and transcriptional mechanisms. Many genes exhibit altered levels of expression in response to abiotic stress, which requires concerted action of both cis- and trans-regulatory features. In order to study the variability in transcriptome response to abiotic stress, RNA sequencing was performed using 14-day-old maize seedlings of inbreds B73, Mo17, Oh43, PH207 and B37 under control, cold and heat conditions. Large numbers of genes that responded differentially to stress between parental inbred lines were identified. RNA sequencing was also performed on similar tissues of the F₁ hybrids produced by crossing B73 and each of the three other inbred lines. By evaluating allele-specific transcript abundance in the F₁ hybrids, we were able to measure the abundance of cis- and trans-regulatory variation between genotypes for both steady-state and stress-responsive expression differences. Although examples of trans-regulatory variation were observed, cis-regulatory variation was more common for both steady-state and stress-responsive expression differences. The genes with cis-allelic variation for response to cold or heat stress provided an opportunity to study the basis for regulatory diversity.

Evolution of meiotic recombination genes in maize and teosinte

Background

Meiotic recombination is a major source of genetic variation in eukaryotes. The role of recombination in evolution is recognized but little is known about how evolutionary forces affect the recombination pathway itself. Although the recombination pathway is fundamentally conserved across different species, genetic variation in recombination components and outcomes has been observed. Theoretical predictions and empirical studies suggest that changes in the recombination pathway are likely to provide adaptive abilities to populations experiencing directional or strong selection pressures, such as those occurring during species domestication. We hypothesized that adaptive changes in recombination may be associated with adaptive evolution patterns of genes involved in meiotic recombination.

Results

To examine how maize evolution and domestication affected meiotic recombination genes, we studied patterns of sequence polymorphism and divergence in eleven genes controlling key steps in the meiotic recombination pathway in a diverse set of maize inbred lines and several accessions of teosinte, the wild ancestor of maize. We discovered that, even though the recombination genes generally exhibited high sequence conservation expected in a pathway controlling a key cellular process, they showed substantial levels and diverse patterns of sequence polymorphism. Among others, we found differences in sequence polymorphism patterns between tropical and temperate maize germplasms. Several recombination genes displayed patterns of polymorphism indicative of adaptive evolution.

Conclusions

Despite their ancient origin and overall sequence conservation, meiotic recombination genes can exhibit extensive and complex patterns of molecular evolution. Changes in these genes could affect the functioning of the recombination pathway, and may have contributed to the successful domestication of maize and its expansion to new cultivation areas.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3486-z) contains supplementary material, which is available to authorized users.