Genomic asymmetry in allopolyploid plants: wheat as a model

Functional Characterization of Accessible Chromatin in Common Wheat

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

Morphological and genetic evidence suggest gene flow among native and naturalized mint species

PREMISE

Cultivation and naturalization of plants beyond their natural range can bring previously geographically isolated taxa together, increasing the opportunity for hybridization, the outcomes of which are not predictable. Here, we explored the phenotypic and genomic effects of interspecific gene flow following the widespread cultivation of Mentha spicata (spearmint), M. longifolia, and M. suaveolens.

METHODS

We morphologically evaluated 155 herbarium specimens of three Mentha species and sequenced the genomes of a subset of 93 specimens. We analyzed the whole genomes in a population and the phylogenetic framework and associated genomic classifications in conjunction with the morphological assessments.

RESULTS

The allopolyploid M. spicata, which likely evolved in cultivation, had altered trichome characters, that is possibly a product of human selection for a more palatable plant or a byproduct of selection for essential oils. There were signs of genetic admixture between mints, including allopolyploids, indicating that the reproductive barriers between Mentha species with differences in ploidy are likely incomplete. Still, despite gene flow between species, we found that genetic variants associated with the cultivated trichome morphology continue to segregate.

CONCLUSIONS

Although hybridization, allopolyploidization, and human selection during cultivation can increase species richness (e.g., by forming hybrid taxa), we showed that unless reproductive barriers are strong, these processes can also result in mixing of genes between species and the potential loss of natural biodiversity.

DNA methylation role in subgenome expression dominance of Juglans regia and its wild relative J. mandshurica

Subgenome expression dominance plays a crucial role in the environmental adaptation of polyploids. However, the epigenetic molecular mechanism underlying this process has not been thoroughly investigated, particularly in perennial woody plants. Persian walnut (Juglans regia) and its wild relative, Manchurian walnut (Juglans mandshurica), are woody plants of great economic importance and are both paleopolyploids that have undergone whole-genome duplication events. In this study, we explored the characteristics of subgenome expression dominance in these 2 Juglans species and examined its epigenetic basis. We divided their genomes into dominant subgenome (DS) and submissive subgenome (SS) and found that the DS-specific genes might play critical roles in biotic stress response or pathogen defense. We comprehensively elucidated the characteristics of biased gene expression, asymmetric DNA methylation, transposable elements (TEs), and alternative splicing (AS) events of homoeologous gene pairs between subgenomes. The results showed that biased expression genes (BEGs) in 2 Juglans species were mainly related to external stimuli response, while non-BEGs were related to complexes that might be involved in signal transduction. DS genes had higher expression and more AS events while having less DNA methylation and TEs than homoeologous genes from the SS in the 2 Juglans species. Further studies showed that DNA methylation might contribute to the biased expression of gene pairs by modifying LTR/TIR/nonTIR TEs and improving the AS efficiency of corresponding precursor mRNAs in a particular context. Our study contributes to understanding the epigenetic basis of subgenome expression dominance and the environmental adaptation of perennial woody plants.

New insights on the evolution of nucleolar dominance in newly resynthesized hexaploid wheat Triticum zhukovskyi

Nucleolar dominance (ND) is a widespread epigenetic phenomenon in hybridizations where nucleolus transcription fails at the nucleolus organizer region (NOR). However, the dynamics of NORs during the formation of Triticum zhukovskyi (GGAu Au Am Am ), another evolutionary branch of allohexaploid wheat, remains poorly understood. Here, we elucidated genetic and epigenetic changes occurring at the NOR loci within the Am , G, and D subgenomes during allopolyploidization by synthesizing hexaploid wheat GGAu Au Am Am and GGAu Au DD. In T. zhukovskyi, Au genome NORs from T. timopheevii (GGAu Au ) were lost, while the second incoming NORs from T. monococcum (Am Am ) were retained. Analysis of the synthesized T. zhukovskyi revealed that rRNA genes from the Am genome were silenced in F1 hybrids (GAu Am ) and remained inactive after genome doubling and subsequent self-pollinations. We observed increased DNA methylation accompanying the inactivation of NORs in the Am genome and found that silencing of NORs in the S1 generation could be reversed by a cytidine methylase inhibitor. Our findings provide insights into the ND process during the evolutionary period of T. zhukovskyi and highlight that inactive rDNA units may serve as a 'first reserve' in the form of R-loops, contributing to the successful evolution of T. zhukovskyi.

Wheat PHT1;9 acts as one candidate arsenate absorption transporter for phytoremediation

Arsenate (AsV) is one of the most common forms of arsenic (As) in environment and plant high-affinity phosphate transporters (PHT1s) are the primary plant AsV transporters. However, few PHT1s involved in AsV absorption have been identified in crops. In our previous study, TaPHT1;3, TaPHT1;6 and TaPHT1;9 were identified to function in phosphate absorption. Here, their AsV absorption capacities were evaluated using several experiments. Ectopic expression in yeast mutants indicated that TaPHT1;9 had the highest AsV absorption rates, followed by TaPHT1;6, while not for TaPHT1;3. Under AsV stress, further, BSMV-VIGS-mediated TaPHT1;9-silencing wheat plants exhibited higher AsV tolerance and lower As concentrations than TaPHT1;6-silenced plants, whereas TaPHT1;3-silencing plants had similar phenotype and AsV concentrations to control. These suggested that TaPHT1;9 and TaPHT1;6 possessed AsV absorption capacity with the former showing higher activities. Under hydroponic condition, furthermore, CRISPR-edited TaPHT1;9 wheat mutants showed the enhanced tolerance to AsV with decreased As distributions and concentrations, whereas TaPHT1;9 ectopic expression transgenic rice plants had the opposite results. Also, under AsV-contaminated soil condition, TaPHT1;9 transgenic rice plants exhibited depressed AsV tolerance with increased As concentrations in roots, straws and grains. Moreover, Pi addition alleviated the AsV toxicity. These suggested that TaPHT1;9 should be a candidate target gene for AsV phytoremediation.

Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis

BACKGROUND

Homoeologs are defined as homologous genes resulting from allopolyploidy. Bread wheat, Triticum aestivum, is an allohexaploid species with many homoeologs. Homoeolog expression bias, referring to the relative contribution of homoeologs to the transcriptome, is critical for determining the traits that influence wheat growth and development. Asymmetric transcription of homoeologs has been so far investigated in a tissue or organ-specific manner, which could be misleading due to a mixture of cell types.

RESULTS

Here, we perform single nuclei RNA sequencing and ATAC sequencing of wheat root to study the asymmetric gene transcription, reconstruct cell differentiation trajectories and cell-type-specific gene regulatory networks. We identify 22 cell types. We then reconstruct cell differentiation trajectories that suggest different origins between epidermis/cortex and endodermis, distinguishing bread wheat from Arabidopsis. We show that the ratio of asymmetrically transcribed triads varies greatly when analyzing at the single-cell level. Hub transcription factors determining cell type identity are also identified. In particular, we demonstrate that TaSPL14 participates in vasculature development by regulating the expression of BAM1. Combining single-cell transcription and chromatin accessibility data, we construct the pseudo-time regulatory network driving root hair differentiation. We find MYB3R4, REF6, HDG1, and GATAs as key regulators in this process.

CONCLUSIONS

Our findings reveal the transcriptional landscape of root organization and asymmetric gene transcription at single-cell resolution in polyploid wheat.

Molecular evidence for adaptive evolution of drought tolerance in wild cereals

The considerable drought tolerance of wild cereal crop progenitors has diminished during domestication in the pursuit of higher productivity. Regaining this trait in cereal crops is essential for global food security but requires novel genetic insight. Here, we assessed the molecular evidence for natural variation of drought tolerance in wild barley (Hordeum spontaneum), wild emmer wheat (Triticum dicoccoides), and Brachypodium species collected from dry and moist habitats at Evolution Canyon, Israel (ECI). We report that prevailing moist vs dry conditions have differentially shaped the stomatal and photosynthetic traits of these wild cereals in their respective habitats. We present the genomic and transcriptomic evidence accounting for differences, including co-expression gene modules, correlated with physiological traits, and selective sweeps, driven by the xeric site conditions on the African Slope (AS) at ECI. Co-expression gene module 'circadian rhythm' was linked to significant drought-induced delay in flowering time in Brachypodium stacei genotypes. African Slope-specific differentially expressed genes are important in barley drought tolerance, verified by silencing Disease-Related Nonspecific Lipid Transfer 1 (DRN1), Nonphotochemical Quenching 4 (NPQ4), and Brassinosteroid-Responsive Ring-H1 (BRH1). Our results provide new genetic information for the breeding of resilient wheat and barley in a changing global climate with increasingly frequent drought events.

Subgenome-biased expression and functional diversification of a Na+/H+ antiporter homoeologs in salt tolerance of polyploid wheat

Common wheat (Triticum aestivum, BBAADD) is an allohexaploid species combines the D genome from Ae. tauschii and with the AB genomes from tetraploid wheat (Triticum turgidum). Compared with tetraploid wheat, hexaploid wheat has wide-ranging adaptability to environmental adversity such as salt stress. However, little is known about the molecular basis underlying this trait. The plasma membrane Na⁺/H⁺ transporter Salt Overly Sensitive 1 (SOS1) is a key determinant of salt tolerance in plants. Here we show that the upregulation of TaSOS1 expression is positively correlated with salt tolerance variation in polyploid wheat. Furthermore, both transcriptional analysis and GUS staining on transgenic plants indicated TaSOS1-A and TaSOS1-B exhibited higher basal expression in roots and leaves in normal conditions and further up-regulated under salt stress; while TaSOS1-D showed markedly lower expression in roots and leaves under normal conditions, but significant up-regulated in roots but not leaves under salt stress. Moreover, transgenic studies in Arabidopsis demonstrate that three TaSOS1 homoeologs display different contribution to salt tolerance and TaSOS1-D plays the prominent role in salt stress. Our findings provide insights into the subgenomic homoeologs variation potential to broad adaptability of natural polyploidy wheat, which might effective for genetic improvement of salinity tolerance in wheat and other crops.

Genome-Wide Association Study in Bread Wheat Identifies Genomic Regions Associated with Grain Yield and Quality under Contrasting Water Availability

The objectives of this study were to identify genetic loci in the bread wheat genome that would influence yield stability and quality under water stress, and to identify accessions that can be recommended for cultivation in dry and hot regions. We performed a genome-wide association study (GWAS) using a panel of 232 wheat accessions spanning diverse ecogeographic regions. Plants were evaluated in the Israeli Northern Negev, under two environments: water-limited (D; 250 mm) and well-watered (W; 450 mm) conditions; they were genotyped with ~71,500 SNPs derived from exome capture sequencing. Of the 14 phenotypic traits evaluated, 12 had significantly lower values under D compared to W conditions, while the values for two traits were higher under D. High heritability (H2 = 0.5-0.9) was observed for grain yield, spike weight, number of grains per spike, peduncle length, and plant height. Days to heading and grain yield could be partitioned based on accession origins. GWAS identified 154 marker-trait associations (MTAs) for yield and quality-related traits, 82 under D and 72 under W, and identified potential candidate genes. We identified 24 accessions showing high and/or stable yields under D conditions that can be recommended for cultivation in regions under the threat of global climate change.

Effects of homoeologous exchange on gene expression and alternative splicing in a newly formed allotetraploid wheat

Homoeologous exchange (HE) is a major mechanism generating post-polyploidization genetic variation with important evolutionary consequences. However, the direct impacts of HE on gene expression and transcript diversity in allopolyploids without the intertwined evolutionary processes remain to be fully understood. Here, we analyzed high-throughput RNA-seq data of young leaves from plant groups of a synthetic allotetraploid wheat (AADD), which contained variable numbers of HEs. We aimed to investigate if and to which extent HE directly impacts gene expression and alternative splicing (AS). We found that HE impacts expression of genes located within HE regions primarily via a cis-acting dosage effect, which led to significant changes in the total expression level of homoeologous gene pairs, especially for homoeologs whose original expression was biased. In parallel, HE also influences expression of a large number of genes residing in non-HE regions by trans-regulation leading to convergent expression of homoeologs. Intriguingly, when taking the original relative homoeolog expression states into account, homoeolog pairs under trans-effect are more prone to manifesting a convergent response to the HEs whereas those under cis-regulation tended to show further exacerbated subgenome-biased expression. Moreover, HE-induced quantitative, largely individual-specific, changes of AS events were detected. Similar to homoeologous expression, homoeo-AS events under trans-effect were more responsive to HE. HE therefore exerts multifaceted immediate effects on gene expression and, to a less extent, on individualized transcript diversity in nascent allopolyploidy.

Quick mapping and characterization of a co-located kernel length and thousand-kernel weight-related QTL in wheat

A co-located KL and TKW-related QTL with no negative effect on PH and AD was rapidly identified using BSA and wheat 660 K SNP array. Its effect was validated in a panel of 218 wheat accessions. Kernel length (KL) and thousand-kernel weight (TKW) of wheat (Triticum aestivum L.) contribute significantly to kernel yield. In the present study, a recombinant inbred line (RIL) population derived from the cross between the wheat line S849-8 with larger kernels and more spikelets per spike and the line SY95-71 was developed. Further, of both the bulked segregant analysis (BSA) and the wheat 660 K single nucleotide polymorphism (SNP) array were used to rapidly identify genomic regions for kernel-related traits from this RIL population. Kompetitive Allele Specific PCR markers were further developed in the SNP-enriched region on the 2D chromosome to construct a genetic map. Both QKL.sicau-SSY-2D for KL and QTKW.sicau-SSY-2D for TKW were identified at multiple environments on chromosome arm 2DL. These two QTLs explained 9.68-23.02% and 6.73-18.32% of the phenotypic variation, respectively. The effects of this co-located QTL were successfully verified in a natural population consisting of 218 Sichuan wheat accessions. Interestingly, the major QTL was significantly and positively correlated with spike length, but did not negatively affect spikelet number per spike (SNS), plant height, or anthesis date. These results indicated that it is possible to synchronously improve kernel weight and SNS by using this QTL. Additionally, several genes associated with kernel development and filling rate were predicted and sequenced in the QTL-containing physical intervals of reference genomes of 'Chinese spring' and Aegilops tauschii. Collectively, these results provide a QTL with great breeding potential and its linked markers which should be helpful for fine mapping and molecular breeding.

Cross-kingdom transcriptomic trends in the evolution of hybrid gene expression

Hybridization is a route to speciation that occurs widely across the eukaryote tree of life. The success of allopolyploids (hybrid species with increased ploidy) and homoploid hybrids (with unchanged ploidy) is well documented. However, their formation and establishment is not straightforward, with a suite of near‐instantaneous and longer term biological repercussions faced by the new species. Central to these challenges is the rewiring of gene regulatory networks following the merger of distinct genomes inherited from both parental species. Research on the evolution of hybrid gene expression has largely involved studies on a single hybrid species or a few gene families. Here, we present the first standardized transcriptome‐wide study exploring the fates of genes following hybridization across three kingdoms: animals, plants and fungi. Within each kingdom, we pair an allopolyploid system with a closely related homoploid hybrid to decouple the influence of increased ploidy from genome merger. Genome merger, not changes in ploidy, has the greatest effect on posthybridization expression patterns across all study systems. Strikingly, we find that differentially expressed genes in parent species preferentially switch to more similar expression in hybrids across all kingdoms, likely as a consequence of regulatory trans‐acting cross‐talk within the hybrid nucleus. We also highlight the prevalence of gene loss or silencing among extremely differentially expressed genes in hybrid species across all kingdoms. These shared patterns suggest that the evolutionary process of hybridization leads to common high‐level expression outcomes, regardless of the particular species or kingdom.

Evolution and origin of bread wheat

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500-9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat's large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.

The Potential Role of RNA Structure in Crop Molecular Breeding

The continually growing human population creates a concomitantly increasing demand for nutritious crops with high yields. Advances in high throughput sequencing technologies have revealed the genetic architecture of major crops. This includes extensive information enabling comprehensive genetic markers for breeding selection, new gene discoveries, and novel gene regulatory strategies for crop editing. RNA structure is an important type of genetic feature, essential for post-transcriptional regulation of gene expression. Here, we summarize recent advances in genome-wide RNA structure studies in crops and review the associated RNA structure-mediated regulation of gene expression. We also discuss the functional importance of those single nucleotide variations that induce large RNA structure disparities. Lastly, we discuss the potential role of RNA structure in crop molecular breeding.

Rubiscosome gene expression is balanced across the hexaploid wheat genome

Functional and active Rubisco is essential for CO₂ fixation and is a primary target for engineering approaches to increasing crop yields. However, the assembly and maintenance of active Rubisco are dependent on the coordinated biosynthesis of at least 11 nuclear-encoded proteins, termed the 'Rubiscosome'. Using publicly available gene expression data for wheat (Triticum aestivum L.), we show that the expression of Rubiscosome genes is balanced across the three closely related subgenomes that form the allohexaploid genome. Each subgenome contains a near complete set of homoeologous genes and contributes equally to overall expression, both under optimal and under heat stress conditions. The expression of the wheat thermo-tolerant Rubisco activase isoform 1β increases under heat stress and remains balanced across the subgenomes, albeit with a slight shift towards greater contribution from the D subgenome. The findings show that the gene copies in all three subgenomes need to be accounted for when designing strategies for crop improvement.

Over-expression of TaDWF4 increases wheat productivity under low and sufficient nitrogen through enhanced carbon assimilation

There is a strong pressure to reduce nitrogen (N) fertilizer inputs while maintaining or increasing current cereal crop yields. We show that overexpression of TaDWF4-B, the dominant shoot expressed homoeologue of OsDWF4, in wheat can increase plant productivity by up to 105% under a range of N levels on marginal soils, resulting in increased N use efficiency (NUE). We show that a two to four-fold increase in TaDWF4 transcript levels enhances the responsiveness of genes regulated by N. The productivity increases seen were primarily due to the maintenance of photosystem II operating efficiency and carbon assimilation in plants when grown under limiting N conditions and not an overall increase in photosynthesis capacity. The increased biomass production and yield per plant in TaDWF4 OE lines could be linked to modified carbon partitioning and changes in expression pattern of the growth regulator Target Of Rapamycin, offering a route towards breeding for sustained yield and lower N inputs.

In wheat, overexpression of TaDWF4 overrides normal nutrient sensing allowing for increased biomass when grown under limiting nutrient conditions. This maintenance of growth is associated with modified carbon partitioning and changes in expression of growth regulator TaTOR, offering a route towards breeding for sustained yields with lower nitrogen inputs.

GWAS for main effects and epistatic interactions for grain morphology traits in wheat

In the present study in wheat, GWAS was conducted for identification of marker trait associations (MTAs) for the following six grain morphology traits: (1) grain cross-sectional area (GCSA), (2) grain perimeter (GP), (3) grain length (GL), (4) grain width (GWid), (5) grain length-width ratio (GLWR) and (6) grain form-density (GFD). The data were recorded on a subset of spring wheat reference set (SWRS) comprising 225 diverse genotypes, which were genotyped using 10,904 SNPs and phenotyped for two consecutive years (2017-2018, 2018-2019). GWAS was conducted using five different models including two single-locus models (CMLM, SUPER), one multi-locus model (FarmCPU), one multi-trait model (mvLMM) and a model for Q x Q epistatic interactions. False discovery rate (FDR) [P value -log10(p) ≥ 5] and Bonferroni correction [P value -log10(p) ≥ 6] (corrected p value < 0.05) were applied to eliminate false positives due to multiple testing. This exercise gave 88 main effect and 29 epistatic MTAs after FDR and 13 main effect and 6 epistatic MTAs after Bonferroni corrections. MTAs obtained after Bonferroni corrections were further utilized for identification of 55 candidate genes (CGs). In silico expression analysis of CGs in different tissues at different parts of the seed at different developmental stages was also carried out. MTAs and CGs identified during the present study are useful addition to available resources for MAS to supplement wheat breeding programmes after due validation and also for future strategic basic research.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01164-w.

Transcriptome Analysis of Developing Wheat Grains at Rapid Expanding Phase Reveals Dynamic Gene Expression Patterns

Simple Summary

Understanding the regulatory mechanism underlying grain development is essential for wheat improvement. The early grain expanding phase boasts critical biological events like embryogenesis and initiation of grain filling. RNA sequencing analysis of this developmental stage revealed dynamic expressions of genes related to cell division, starch biosynthesis, and hormone biosynthesis. An unbalanced expression among triads may play critical roles as shown by multiple enriched metabolic pathways. Our work demonstrated complex regulation mechanisms in early grain development and provided useful information for future wheat improvement.

Abstract

Grain development, as a vital process in the crop’s life cycle, is crucial for determining crop quality and yield. The wheat grain expanding phase is the early process involving the rapid morphological changes and initiation of grain filling. However, little is known about the molecular basis of grain development at this stage. Here, we provide a time-series transcriptome profile of developing wheat grain at 0, 2, 4, 6, 8, and 10 days after pollination of the wheat landrace Chinese Spring. A total of 26,892 differentially expressed genes, including 1468 transcription factors, were found between adjacent time points. Co-expression cluster analysis and Gene Ontology enrichment revealed dynamic expressions of cell division and starch biosynthesis related structural genes and transcription factors. Moreover, diverse, differential and drastically varied expression trends of the key genes related to hormone metabolism were identified. Furthermore, ~30% of triads showed unbalanced expression patterns enriching for genes in multiple pivotal metabolic pathways. Hormone metabolism related genes, such as YUC10 (YUCCA flavin-containing monooxygenase 10), AOS2 (allene oxide synthase 2), CYP90D2 (cytochrome P450 90D2), and CKX1 (cytokinin dehydrogenase 1), were dominantly contributed by A or D homoeologs of the triads. Our study provided a systematic picture of transcriptional regulation of wheat grains at the early grain expanding phase which should deepen our understanding of wheat grain development and help in wheat yield improvement.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Chloroplast Phylogenomic Analyses Reveal a Maternal Hybridization Event Leading to the Formation of Cultivated Peanuts

Peanuts (Arachis hypogaea L.) offer numerous healthy benefits, and the production of peanuts has a prominent role in global food security. As a result, it is in the interest of society to improve the productivity and quality of peanuts with transgenic means. However, the lack of a robust phylogeny of cultivated and wild peanut species has limited the utilization of genetic resources in peanut molecular breeding. In this study, a total of 33 complete peanut plastomes were sequenced, analyzed and used for phylogenetic analyses. Our results suggest that sect. Arachis can be subdivided into two lineages. All the cultivated species are contained in Lineage I with AABB and AA are the two predominant genome types present, while species in Lineage II possess diverse genome types, including BB, KK, GG, etc. Phylogenetic studies also indicate that all allotetraploid cultivated peanut species have been derived from a possible maternal hybridization event with one of the diploid Arachis duranensis accessions being a potential AA sub-genome ancestor. In addition, Arachis monticola, a tetraploid wild species, is placed in the same group with all the cultivated peanuts, and it may represent a transitional species, which has been through the recent hybridization event. This research could facilitate a better understanding of the taxonomic status of various Arachis species/accessions and the evolutionary relationship among them, and assists in the correct and efficient use of germplasm resources in breeding efforts to improve peanuts for the benefit of human beings.

Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry

BACKGROUND

Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry.

RESULTS

Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat.

CONCLUSIONS

We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

Chromosome-scale assembly and evolution of the tetraploid Salvia splendens (Lamiaceae) genome

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

Histone acetyltransferase TaHAG1 acts as a crucial regulator to strengthen salt tolerance of hexaploid wheat

Polyploidy occurs prevalently and plays an important role during plant speciation and evolution. This phenomenon suggests polyploidy could develop novel features that enable them to adapt wider range of environmental conditions compared with diploid progenitors. Bread wheat (Triticum aestivum L., BBAADD) is a typical allohexaploid species and generally exhibits greater salt tolerance than its tetraploid wheat progenitor (BBAA). However, little is known about the underlying molecular basis and the regulatory pathway of this trait. Here, we show that the histone acetyltransferase TaHAG1 acts as a crucial regulator to strengthen salt tolerance of hexaploid wheat. Salinity-induced TaHAG1 expression was associated with tolerance variation in polyploidy wheat. Overexpression, silencing, and CRISPR-mediated knockout of TaHAG1 validated the role of TaHAG1 in salinity tolerance of wheat. TaHAG1 contributed to salt tolerance by modulating reactive oxygen species (ROS) production and signal specificity. Moreover, TaHAG1 directly targeted a subset of genes that are responsible for hydrogen peroxide production, and enrichment of TaHAG1 triggered increased H3 acetylation and transcriptional upregulation of these loci under salt stress. In addition, we found the salinity-induced TaHAG1-mediated ROS production pathway is involved in salt tolerance difference of wheat accessions with varying ploidy. Our findings provide insight into the molecular mechanism of how an epigenetic regulatory factor facilitates adaptability of polyploidy wheat and highlights this epigenetic modulator as a strategy for salt tolerance breeding in bread wheat.

Reconstruction of proto-vertebrate, proto-cyclostome and proto-gnathostome genomes provides new insights into early vertebrate evolution

Ancient polyploidization events have had a lasting impact on vertebrate genome structure, organization and function. Some key questions regarding the number of ancient polyploidization events and their timing in relation to the cyclostome-gnathostome divergence have remained contentious. Here we generate de novo long-read-based chromosome-scale genome assemblies for the Japanese lamprey and elephant shark. Using these and other representative genomes and developing algorithms for the probabilistic macrosynteny model, we reconstruct high-resolution proto-vertebrate, proto-cyclostome and proto-gnathostome genomes. Our reconstructions resolve key questions regarding the early evolutionary history of vertebrates. First, cyclostomes diverged from the lineage leading to gnathostomes after a shared tetraploidization (1R) but before a gnathostome-specific tetraploidization (2R). Second, the cyclostome lineage experienced an additional hexaploidization. Third, 2R in the gnathostome lineage was an allotetraploidization event, and biased gene loss from one of the subgenomes shaped the gnathostome genome by giving rise to remarkably conserved microchromosomes. Thus, our reconstructions reveal the major evolutionary events and offer new insights into the origin and evolution of vertebrate genomes.

The endosperm-specific transcription factor TaNAC019 regulates glutenin and starch accumulation and its elite allele improves wheat grain quality

In wheat (Triticum aestivum L.), breeding efforts have focused intensively on improving grain yield and quality. For quality, the content and composition of seed storage proteins (SSPs) determine the elasticity of wheat dough and flour processing quality. Moreover, starch levels in seeds are associated with yield. However, little is known about the mechanisms that coordinate SSP and starch accumulation in wheat. In this study, we explored the role of the endosperm-specific NAC transcription factor TaNAC019 in coordinating SSP and starch accumulation. TaNAC019 binds to the promoters of TaGlu-1 loci, encoding high molecular weight glutenin (HMW-GS), and of starch metabolism genes. Triple knock-out mutants of all three TaNAC019 homoeologs exhibited reduced transcript levels for all SSP types and genes involved in starch metabolism, leading to lower gluten and starch contents, and in flour processing quality parameters. TaNAC019 directly activated the expression of HMW-GS genes by binding to a specific motif in their promoters and interacting with the TaGlu-1 regulator TaGAMyb. TaNAC019 also indirectly regulated the expression of TaSPA, an ortholog of maize Opaque2 that activates SSP accumulation. Therefore, TaNAC019 regulation of starch- and SSP-related genes has key roles in wheat grain quality. Finally, we identified an elite allele (TaNAC019-BI) associated with flour processing quality, providing a candidate gene for breeding wheat with improved quality.

Conservation and trans-regulation of histone modification in the A and B subgenomes of polyploid wheat during domestication and ploidy transition

BACKGROUND

Polyploidy has played a prominent role in the evolution of plants and many other eukaryotic lineages. However, how polyploid genomes adapt to the abrupt presence of two or more sets of chromosomes via genome regulation remains poorly understood. Here, we analyzed genome-wide histone modification and gene expression profiles in relation to domestication and ploidy transition in the A and B subgenomes of polyploid wheat.

RESULTS

We found that epigenetic modification patterns by two typical euchromatin histone markers, H3K4me3 and H3K27me3, for the great majority of homoeologous triad genes in A and B subgenomes were highly conserved between wild and domesticated tetraploid wheats and remained stable in the process of ploidy transitions from hexaploid to extracted tetraploid and then back to resynthesized hexaploid. However, a subset of genes was differentially modified during tetraploid and hexaploid wheat domestication and in response to ploidy transitions, and these genes were enriched for particular gene ontology (GO) terms. The extracted tetraploid wheat manifested higher overall histone modification levels than its hexaploid donor, and which were reversible and restored to normal levels in the resynthesized hexaploid. Further, while H3K4me3 marks were distally distributed along each chromosome and significantly correlated with subgenome expression as expected, H3K27me3 marks showed only a weak distal bias and did not show a significant correlation with gene expression.

CONCLUSIONS

Our results reveal overall high stability of histone modification patterns in the A and B subgenomes of polyploid wheat during domestication and in the process of ploidy transitions. However, modification levels of a subset of functionally relevant genes in the A and B genomes were trans-regulated by the D genome in hexaploid wheat.

New insights into homoeologous copy number variations in the hexaploid wheat genome

Bread wheat is an allohexaploid species originating from two successive and recent rounds of hybridization between three diploid species that were very similar in terms of chromosome number, genome size, TE content, gene content and synteny. As a result, it has long been considered that most of the genes were in three pairs of homoeologous copies. However, these so-called triads represent only one half of wheat genes, while the remaining half belong to homoeologous groups with various number of copies across subgenomes. In this study, we examined and compared the distribution, conservation, function, expression and epigenetic profiles of triads with homoeologous groups having undergone a deletion (dyads) or a duplication (tetrads) in one subgenome. We show that dyads and tetrads are mostly located in distal regions and have lower expression level and breadth than triads. Moreover, they are enriched in functions related to adaptation and more associated with the repressive H3K27me3 modification. Altogether, these results suggest that triads mainly correspond to housekeeping genes and are part of the core genome, while dyads and tetrads belong to the Triticeae dispensable genome. In addition, by comparing the different categories of dyads and tetrads, we hypothesize that, unlike most of the allopolyploid species, subgenome dominance and biased fractionation are absent in hexaploid wheat. Differences observed between the three subgenomes are more likely related to two successive and ongoing waves of post-polyploid diploidization, that had impacted A and B more significantly than D, as a result of the evolutionary history of hexaploid wheat.

Functional characterization and regulatory mechanism of wheat CPK34 kinase in response to drought stress

BACKGROUND

Drought is one of the most adverse environmental factors limiting crop productions and it is important to identify key genetic determinants for food safety. Calcium-dependent protein kinases (CPKs) are known to be involved in plant growth, development, and environmental stresses. However, biological functions and regulatory mechanisms of many plant CPKs have not been explored. In our previous study, abundance of the wheat CPK34 (TaCPK34) protein was remarkably upregulated in wheat plants suffering from drought stress, inferring that it could be involved in this stress. Therefore, here we further detected its function and mechanism in response to drought stress.

RESULTS

Transcripts of the TaCPK34 gene were significantly induced after PEG-stimulated water deficiency (20% PEG6000) or 100 μM abscisic acid (ABA) treatments. The TaCPK34 gene was transiently silenced in wheat genome by using barley stripe mosaic virus-induced silencing (BSMV-VIGS) method. After 14 days of drought stress, the transiently TaCPK34-silenced wheat seedlings showed more sensitivity compared with control, and the plant biomasses and relative water contents significantly decreased, whereas soluble sugar and MDA contents increased. The iTRAQ-based quantitative proteomics was employed to measure the protein expression profiles in leaves of the transiently TaCPK34-silenced wheat plants after drought stress. There were 6103 proteins identified, of these, 51 proteins exhibited significantly altered abundance, they were involved in diverse function. And sequence analysis on the promoters of genes, which encoded the above identified proteins, indicated that some promoters harbored some ABA-responsive elements. We determined the interactions between TaCPK34 and three identified proteins by using bimolecular fluorescent complementation (BiFC) method and our data indicated that TaCPK34directly interacted with the glutathione S-transferase 1 and prx113, respectively.

CONCLUSIONS

Our study suggested that the TaCPK34 gene played positive roles in wheat response to drought stress through directly or indirectly regulating the expression of ABA-dependent manner genes, which were encoding identified proteins from iTRAQ-based quantitative proteomics. And it could be used as one potential gene to develop crop cultivars with improved drought tolerance.

Current Progress in Understanding and Recovering the Wheat Genes Lost in Evolution and Domestication

The modern cultivated wheat has passed a long evolution involving origin of wild emmer (WEM), development of cultivated emmer, formation of spelt wheat and finally establishment of modern bread wheat and durum wheat. During this evolutionary process, rapid alterations and sporadic changes in wheat genome took place, due to hybridization, polyploidization, domestication, and mutation. This has resulted in some modifications and a high level of gene loss. As a result, the modern cultivated wheat does not contain all genes of their progenitors. These lost genes are novel for modern wheat improvement. Exploring wild progenitor for genetic variation of important traits is directly beneficial for wheat breeding. WEM wheat (Triticum dicoccoides) is a great genetic resource with huge diversity for traits. Few genes and quantitative trait loci (QTL) for agronomic, quantitative, biotic and abiotic stress-related traits have already been mapped from WEM. This resource can be utilized for modern wheat improvement by integrating identified genes or QTLs through breeding.

Development of RNA-seq-based molecular markers for characterizing Thinopyrum bessarabicum and Secale introgressions in wheat

Sequence-based markers have added a new dimension in the efficiency of identifying alien introgressions in wheat. Expressed sequence tag-sequence tagged sites (EST-STS) markers have proved useful in tracing alien chromatin. In this study, we report the development of Thinopyrum bessarabicum- and Secale anatolicum-specific EST-STS markers and their application in tracing respective alien chromatin introgressions in wheat. The parental lines, Chinese Spring (CS), ISR991.1 (CS/Th. bessarabicum amphidiploid), and ISR1049.2 (CS/Secale anatolicum amphidiploid), were used as core experimental materials. Using comparative analysis of RNA-Seq data, 10 903 and 10 660 candidate sequences specific to Th. bessarabicum and S. anatolicum, respectively, were assembled and identified. To validate the genome specificity of these candidate sequences, 68 and 64 EST-STS markers were developed from randomly selected candidate sequences of Th. bessarabicum and S. anatolicum, respectively, and tested on sets of alien addition lines. Fifty-five and 53 markers for Th. bessarabicum and S. anatolicum chromatin, respectively, were assigned to chromosomal location(s), covering all seven chromosomes. Approximately 83% of S. anatolicum-specific markers were transferable to S. cereale. The genome-specific candidate sequences identified and the EST-STS markers developed will be valuable resources for exploitation of Th. bessarabicum and Secale species diversity in wheat and triticale breeding.

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.

Illegitimate Recombination Between Homeologous Genes in Wheat Genome

Polyploidies produce a large number of duplicated regions and genes in genomes, which have a long-term impact and stimulate genetic innovation. The high similarity between homeologous chromosomes, forming different subgenomes, or homologous regions after genome repatterning, may permit illegitimate DNA recombination. Here, based on gene colinearity, we aligned the (sub)genomes of common wheat (Triticum aestivum, AABBDD genotype) and its relatives, including Triticum urartu (AA), Aegilops tauschii (DD), and T. turgidum ssp. dicoccoides (AABB) to detect the homeologous (paralogous or orthologous) colinear genes within and between (sub)genomes. Besides, we inferred more ancient paralogous regions produced by a much ancient grass-common tetraploidization. By comparing the sequence similarity between paralogous and orthologous genes, we assumed abnormality in the topology of constructed gene trees, which could be explained by gene conversion as a result of illegitimate recombination. We found large numbers of inferred converted genes (>2,000 gene pairs) suggested long-lasting genome instability of the hexaploid plant, and preferential donor roles by DD genes. Though illegitimate recombination was much restricted, duplicated genes produced by an ancient whole-genome duplication, which occurred millions of years ago, also showed evidence of likely gene conversion. As to biological function, we found that ~40% catalytic genes in colinearity, including those involved in starch biosynthesis, were likely affected by gene conversion. The present study will contribute to understanding the functional and structural innovation of the common wheat genome.

Independent evolution of transcript abundance and gene regulatory dynamics

Changes in gene expression drive novel phenotypes, raising interest in how gene expression evolves. In contrast to the static genome, cells modulate gene expression in response to changing environments. Previous comparative studies focused on specific conditions, describing interspecies variation in expression levels, but providing limited information about variation across different conditions. To close this gap, we profiled mRNA levels of two related yeast species in hundreds of conditions and used coexpression analysis to distinguish variation in the dynamic pattern of gene expression from variation in expression levels. The majority of genes whose expression varied between the species maintained a conserved dynamic pattern. Cases of diverged dynamic pattern correspond to genes that were induced under distinct subsets of conditions in the two species. Profiling the interspecific hybrid allowed us to distinguish between genes with predominantly cis- or trans-regulatory variation. We find that trans-varying alleles are dominantly inherited, and that cis-variations are often complemented by variations in trans Based on these results, we suggest that gene expression diverges primarily through changes in expression levels, but does not alter the pattern by which these levels are dynamically regulated.

Allele mining of TaGRF-2D gene 5'-UTR in Triticum aestivum and Aegilops tauschii genotypes

The low diversity of the D-subgenome of bread wheat requires the involvement of new alleles for breeding. In grasses, the allelic state of Growth Regulating Factor (GRF) gene is correlated with nitrogen uptake. In this study, we characterized the sequence of TaGRF-2D and assessed its diversity in bread wheat and goatgrass Aegilops tauschii (genome DD). In silico analysis was performed for reference sequence searching, primer pairs design and sequence assembly. The gene sequence was obtained using Illumina and Sanger sequencing. The complete sequences of TaGRF-2D were obtained for 18 varieties of wheat. The polymorphism in the presence/absence of two GCAGCC repeats in 5' UTR was revealed and the GRF-2D-SSR marker was developed. Our results showed that the alleles 5' UTR-250 and 5' UTR-238 were present in wheat varieties, 5' UTR-250 was presented in the majority of wheat varieties. In Ae. tauschii ssp. strangulata (likely donor of the D-subgenome of polyploid wheat), most accessions carried the 5' UTR-250 allele, whilst most Ae. tauschii ssp. tauschii have 5' UTR-244. The developed GRF-2D-SSR marker can be used to study the genetic diversity of wheat and Ae. tauschii.

A high-throughput BAC end analysis protocol (BAC-anchor) for profiling genome assembly and physical mapping

Traditional approaches for sequencing insertion ends of bacterial artificial chromosome (BAC) libraries are laborious and expensive, which are currently some of the bottlenecks limiting a better understanding of the genomic features of auto- or allopolyploid species. Here, we developed a highly efficient and low-cost BAC end analysis protocol, named BAC-anchor, to identify paired-end reads containing large internal gaps. Our approach mainly focused on the identification of high-throughput sequencing reads carrying restriction enzyme cutting sites and searching for large internal gaps based on the mapping locations of both ends of the reads. We sequenced and analysed eight libraries containing over 3 200 000 BAC end clones derived from the BAC library of the tetraploid potato cultivar C88 digested with two restriction enzymes, Cla I and Mlu I. About 25% of the BAC end reads carrying cutting sites generated a 60-100 kb internal gap in the potato DM reference genome, which was consistent with the mapping results of Sanger sequencing of the BAC end clones and indicated large differences between autotetraploid and haploid genotypes in potato. A total of 5341 Cla I- and 165 Mlu I-derived unique reads were distributed on different chromosomes of the DM reference genome and could be used to establish a physical map of target regions and assemble the C88 genome. The reads that matched different chromosomes are especially significant for the further assembly of complex polyploid genomes. Our study provides an example of analysing high-coverage BAC end libraries with low sequencing cost and is a resource for further genome sequencing studies.

Genomic resources for energy cane breeding in the post genomics era

Sugarcane is one of the most sustainable energy crops among cultivated crops presenting the highest tonnage of cultivated plants. Its high productivity of sugar, bioethanol and bioelectricity make it a promising green alternative to petroleum. Furthermore, the myriad of products that can be derived from sugarcane biomass has been driving breeding programs towards varieties with a higher yield of fiber and a more vigorous and sustainable performance: the energy cane. Here we provide an overview of the energy cane including plant description, breeding efforts, types, and end-uses. In addition, we describe recently published genomic resources for the development of this crop, discuss current knowledge of cell wall metabolism, bioinformatic tools and databases available for the community.

The subgenomes show asymmetric expression of alleles in hybrid lineages of Megalobrama amblycephala × Culter alburnus

Hybridization drives rapid speciation by shaping novel genotypic and phenotypic profiles. Genomic incompatibility and transcriptome shock have been observed in hybrids, although this is rarer in animals than in plants. Using the newly sequenced genomes of the blunt snout bream (Megalobrama amblycephala [BSB]) and the topmouth culter (Culter alburnus [TC]), we focused on the sequence variation and gene expression changes in the reciprocal intergeneric hybrid lineages (F₁-F₃) of BSB × TC. A genome-wide transcriptional analysis identified 145-974 expressed recombinant genes in the successive generations of hybrid fish, suggesting the rapid emergence of allelic variation following hybridization. Some gradual changes of gene expression with additive and dominance effects and various cis and trans regulations were observed from F₁ to F₃ in the two hybrid lineages. These asymmetric patterns of gene expression represent the alternative strategies for counteracting deleterious effects of the subgenomes and improving adaptability of novel hybrids. Furthermore, we identified positive selection and additive expression patterns in transforming growth factor, beta 1b (tgfb1b), which may account for the morphological variations of the pharyngeal jaw in the two hybrid lineages. Our current findings provide insights into the evolution of vertebrate genomes immediately following hybridization.

Transcriptional reference map of hormone responses in wheat spikes

Background

Phytohormones are key regulators of plant growth, development, and signalling networks involved in responses to diverse biotic and abiotic stresses. Transcriptional reference maps of hormone responses have been reported for several model plant species such as Arabidopsis thaliana, Oryza sativa, and Brachypodium distachyon. However, because of species differences and the complexity of the wheat genome, these transcriptome data are not appropriate reference material for wheat studies.

Results

We comprehensively analysed the transcriptomic responses in wheat spikes to seven phytohormones, including indole acetic acid (IAA), gibberellic acid (GA), abscisic acid (ABA), ethylene (ET), cytokinin (CK), salicylic acid (SA), and methyl jasmonic acid (MeJA). A total of 3386 genes were differentially expressed at 24 h after the hormone treatments. Furthermore, 22.7% of these genes exhibited overlapping transcriptional responses for at least two hormones, implying there is crosstalk among phytohormones. We subsequently identified genes with expression levels that were significantly and differentially induced by a specific phytohormone (i.e., hormone-specific responses). The data for these hormone-responsive genes were then compared with the transcriptome data for wheat spikes exposed to biotic (Fusarium head blight) and abiotic (water deficit) stresses.

Conclusion

Our data were used to develop a transcriptional reference map of hormone responses in wheat spikes.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5726-x) contains supplementary material, which is available to authorized users.

Broadening the bread wheat D genome

Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. In this review, we discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. The D genome of allohexaploid bread wheat (Triticum aestivum, 2n = AABBDD) is the least diverse of the three wheat genomes and is unarguably less diverse than that of diploid progenitor Aegilops tauschii (2n = DD). Useful genetic variation and phenotypic traits also exist within each of the wheat group species containing a copy of the D genome: allopolyploid Aegilops species Ae. cylindrica (2n = D^cD^cC^cC^c), Ae. crassa 4x (2n = D¹D¹X^crX^cr), Ae. crassa 6x (2n = D¹D¹X^crX^crD^crD^cr), Ae. ventricosa (2n = D^vD^vN^vN^v), Ae. vavilovii (2n = D¹D¹X^crX^crS^vS^v) and Ae. juvenalis (2n = D¹D¹X^crX^crU^jU^j). Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. Some of these D genomes appear to be modified relative to the bread wheat and Ae. tauschii D genomes, and others present in the allopolyploids may also contain useful variation as a result of adaptation to an allopolyploid, multi-genome environment. We summarise the genetic relationships, karyotypic variation and phenotypic traits known to be present in each of the D genome species that could be of relevance for bread wheat improvement and discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. Better understanding of factors controlling chromosome inheritance and recombination in wheat group interspecific hybrids, as well as effective utilisation of new and developing genetics and genomics technologies, have great potential to improve the agronomic potential of the bread wheat D genome.

The Polymorphisms of Oligonucleotide Probes in Wheat Cultivars Determined by ND-FISH

Non-denaturing fluorescence in situ hybridization (ND-FISH) has been used to distinguish wheat chromosomes and to detect alien chromosomes in the wheat genome. In this study, five different oligonucleotide probes were used with ND-FISH to examine 21 wheat cultivars and lines. These oligonucleotide probes distinguished 42 wheat chromosomes and also detected rye chromatin in the wheat genome. Moreover, the signal patterns of the oligonucleotide probes Oligo-pTa535-1 and Oligo-pSc119.2-1 showed high polymorphism in the wheat chromosomes. A total of 17.6% of the A group chromosomes, 25.9% of the B group chromosomes and 8.9% of the D group chromosomes showed obvious mutations when they were compared to the standard ND-FISH signal patterns, and most of them were Oligo-pSc119.2-1 mutants. The results suggested that these polymorphisms could be induced by the crossing of wheat cultivars. The results provided more information for the further application of oligonucleotide probes and ND-FISH.

Prediction of Subgenome Additive and Interaction Effects in Allohexaploid Wheat

Whole genome duplications have played an important role in the evolution of angiosperms. These events often occur through hybridization between closely related species, resulting in an allopolyploid with multiple subgenomes. With the availability of affordable genotyping and a reference genome to locate markers, breeders of allopolyploids now have the opportunity to manipulate subgenomes independently. This also presents a unique opportunity to investigate epistatic interactions between homeologous orthologs across subgenomes. We present a statistical framework for partitioning genetic variance to the subgenomes of an allopolyploid, predicting breeding values for each subgenome, and determining the importance of inter-genomic epistasis. We demonstrate using an allohexaploid wheat breeding population evaluated in Ithaca, NY and an important wheat dataset from CIMMYT previously shown to demonstrate non-additive genetic variance. Subgenome covariance matrices were constructed and used to calculate subgenome interaction covariance matrices for variance component estimation and genomic prediction. We propose a method to extract population structure from all subgenomes at once before covariances are calculated to reduce collinearity between subgenome estimates. Variance parameter estimation was shown to be reliable for additive subgenome effects, but was less reliable for subgenome interaction components. Predictive ability was equivalent to current genomic prediction methods. Including only inter-genomic interactions resulted in the same increase in accuracy as modeling all pairwise marker interactions. Thus, we provide a new tool for breeders of allopolyploid crops to characterize the genetic architecture of existing populations, determine breeding goals, and develop new strategies for selection of additive effects and fixation of inter-genomic epistasis.

Homeologous Epistasis in Wheat: The Search for an Immortal Hybrid

Hybridization between related species results in the formation of an allopolyploid with multiple subgenomes. These subgenomes will each contain complete, yet evolutionarily divergent, sets of genes. Like a diploid hybrid, allopolyploids will have two versions, or homeoalleles, for every gene. Partial functional redundancy between homeologous genes should result in a deviation from additivity. These epistatic interactions between homeoalleles are analogous to dominance effects, but are fixed across subgenomes through self pollination. An allopolyploid can be viewed as an immortalized hybrid, with the opportunity to select and fix favorable homeoallelic interactions within inbred varieties. We present a subfunctionalization epistasis model to estimate the degree of functional redundancy between homeoallelic loci and a statistical framework to determine their importance within a population. We provide an example using the homeologous dwarfing genes of allohexaploid wheat, Rht-1, and search for genome-wide patterns indicative of homeoallelic subfunctionalization in a breeding population. Using the IWGSC RefSeq v1.0 sequence, 23,796 homeoallelic gene sets were identified and anchored to the nearest DNA marker to form 10,172 homeologous marker sets. Interaction predictors constructed from products of marker scores were used to fit the homeologous main and interaction effects, as well as estimate whole genome genetic values. Some traits displayed a pattern indicative of homeoallelic subfunctionalization, while other traits showed a less clear pattern or were not affected. Using genomic prediction accuracy to evaluate importance of marker interactions, we show that homeologous interactions explain a portion of the nonadditive genetic signal, but are less important than other epistatic interactions.

Sympatric and allopatric niche shift of endemic Gypsophila (Caryophyllaceae) taxa in the Iberian Peninsula

Several species of the Gypsophila genus are endemic to the Iberian Peninsula, including gypsophytes of particular ecological, evolutionary and biochemical interest, and taxa that have undergone both sympatric and allopatric genetic differentiation. The niche shift among these taxa has been assessed using ecological niche modelling and ordination techniques, adopting a niche overlap approach to compare the similarity and equivalency of the ecological niches. We used the Maximum Entropy method to study the potential distribution of these taxa in different eras: the Last Glacial Maximum (LGM), the Mid Holocene and the current conditions. We present evidence of niche shift during the speciation of G. bermejoi, with a strong niche overlap between the parental taxa (G. struthium subsp. struthium and G. tomentosa), yet both overlap much more weakly with the hybrid species. This phenomenon may be explained by genetic and epigenetic interactions, and it has been described in other species. We also studied the sister subspecies G. struthium subsp. struthium and G. struthium subsp. hispanica, with mostly allopatric distributions and with the Iberian System mountain range acting as a geographical barrier. The Iberian System and other mountain ranges may have favored differences in the climatic conditions on either side of the mountain range, which is consistent with an incipient process of bioclimatic ecological speciation. These results seem to indicate that niche shift can occur over very different timespans. In the case of G. bermejoi, speciation may have produced significant niche shifting in one or two generations due to its alloploid nature. By contrast, G. struthium subsp. struthium and G. struthium subsp. hispanica seem to have undergone a more gradual process of allopatric genetic differentiation driven by bioclimatic factors. Both these processes are relatively recent and they will have been strongly influenced by the climate change at the end of LGM.

The non-random patterns of genetic variation induced by asymmetric somatic hybridization in wheat

BACKGROUND

Asymmetric somatic hybridization is an efficient crop breeding approach by introducing several exogenous chromatin fragments, which leads to genomic shock and therefore induces genome-wide genetic variation. However, the fundamental question concerning the genetic variation such as whether it occurs randomly and suffers from selection pressure remains unknown.

RESULTS

Here, we explored this issue by comparing expressed sequence tags of a common wheat cultivar and its asymmetric somatic hybrid line. Both nucleotide substitutions and indels (insertions and deletions) had lower frequencies in coding sequences than in un-translated regions. The frequencies of nucleotide substitutions and indels were both comparable between chromosomes with and without introgressed fragments. Nucleotide substitutions distributed unevenly and were preferential to indel-flanking sequences, and the frequency of nucleotide substitutions at 5'-flanking sequences of indels was obviously higher in chromosomes with introgressed fragments than in those without exogenous fragment. Nucleotide substitutions and indels both had various frequencies among seven groups of allelic chromosomes, and the frequencies of nucleotide substitutions were strongly negatively correlative to those of indels. Among three sets of genomes, the frequencies of nucleotide substitutions and indels were both heterogeneous, and the frequencies of nucleotide substitutions exhibited drastically positive correlation to those of indels.

CONCLUSIONS

Our work demonstrates that the genetic variation induced by asymmetric somatic hybridization is attributed to both whole genomic shock and local chromosomal shock, which is a predetermined and non-random genetic event being closely associated with selection pressure. Asymmetric somatic hybrids provide a worthwhile model to further investigate the nature of genomic shock induced genetic variation.

Three genomes differentially contribute to the seedling lateral root number in allohexaploid wheat: evidence from phenotype evolution and gene expression

Common wheat is an allohexaploid (BBAADD) that originated from the hybridization and polyploidization of the diploid Aegilops tauschii (DD) with the allotetraploid Triticum turgidum (BBAA). Phenotypic changes often arise with the formation and evolution of allopolyploid wheat, but little is known about the evolution of root traits in different wheat species with varying ploidy levels. Here, we reported that the lateral root number on the primary root (LRNPR) of synthetic and natural allohexaploid wheats (BBAADD) is significantly higher than that of their allotetraploid (BBAA) and diploid (AA and SS) progenitors, but is much lower than that of their diploid (DD) progenitors. The expression of the wheat gene TaLBD16, an ortholog of the Arabidopsis LATERAL ORGAN BOUNDARIES-DOMAIN16/ASYMMETRIC LEAVES2-LIKE18 (LBD16), which is involved in lateral root development in Arabidopsis, was positively correlated with the LRNPR in diploid and allopolyploid wheats. In natural and synthetic allohexaploid wheats, the transcript of the TaLBD16 from the D genome (TaLBD16-D) was relatively more abundant compared with TaLBD16-A and TaLBD16-B. Consistent with the observed variation in LRNPR, the divergence in the expression of TaLBD16 homoeologous genes occurred before the formation of polyploidy wheat. Collectively, our observations indicate that the D genome played a crucial role in the increased lateral root number of allohexaploid wheats compared with their allotetraploid progenitors, and that TaLBD16-D was one of the key genes involved in the formation of lateral root number during wheat evolution.

The transcriptional landscape of polyploid wheat

The coordinated expression of highly related homoeologous genes in polyploid species underlies the phenotypes of many of the world's major crops. Here we combine extensive gene expression datasets to produce a comprehensive, genome-wide analysis of homoeolog expression patterns in hexaploid bread wheat. Bias in homoeolog expression varies between tissues, with ~30% of wheat homoeologs showing nonbalanced expression. We found expression asymmetries along wheat chromosomes, with homoeologs showing the largest inter-tissue, inter-cultivar, and coding sequence variation, most often located in high-recombination distal ends of chromosomes. These transcriptionally dynamic genes potentially represent the first steps toward neo- or subfunctionalization of wheat homoeologs. Coexpression networks reveal extensive coordination of homoeologs throughout development and, alongside a detailed expression atlas, provide a framework to target candidate genes underpinning agronomic traits in wheat.

Remapping of the stripe rust resistance gene Yr10 in common wheat

Yr10 is an important gene to control wheat stripe rust, and the search for Yr10 needs to be continued. Wheat stripe rust or yellow rust is a devastating fungal disease caused by Puccinia striiformis f. sp. tritici (Pst). Host disease resistance offers a primary source for controlling wheat stripe rust. The stripe rust resistance gene Yr10 confers the race-specific resistance to most tested Pst races in China including CYR29. Early studies proposed that Yr10 was a nucleotide-binding site, leucine-rich repeat gene archived as GenBank accession AF149112 (hereafter designated the Yr10 candidate gene or Yr10 _CG ). In this study, we revealed that 15 Chinese wheat cultivars positive for Yr10 _CG are susceptible to CYR29. We then expressed the Yr10 _CG cDNA in the common wheat 'Bobwhite'. The Yr10 _CG -cDNA positive transgenic plants were also susceptible to CYR29. Thus, it is highly unlikely that Yr10 _CG corresponds to the Yr10 resistance gene. Using the Yr10 donor 'Moro' and the Pst-susceptible wheat 'Huixianhong', we generated two F₃ populations that displayed a single Mendelian segregation on the Yr10 gene, and used them to remap the Yr10 gene. Six markers were placed in the Yr10 region, with the Yr10 _CG gene now mapping about 1.2-cM proximal to the Yr10 locus and the Xsdauw79 marker is completely linked to the Yr10 locus. Apparently, the Yr10 gene has not yet been identified. Fine mapping and positional cloning of Yr10 is important for gene pyramiding for stripe rust resistance in wheat.