Genetic and epigenetic alteration among three homoeologous genes of a class E MADS box gene in hexaploid wheat

Genome-wide analysis of wheat xyloglucan endotransglucosylase/hydrolase (XTH) gene family revealed TaXTH17 involved in abiotic stress responses

BACKGROUND

Environmental stresses, including high salinity and drought, severely diminish wheat yield and quality globally. The xyloglucan endotransglucosylase/hydrolase (XTH) family represents a class of cell wall-modifying enzymes and plays important roles in plants growth, development and stress adaptation. However, systematic analyses of XTH family genes and their functions under salt and drought stresses have not been undertaken in wheat.

RESULTS

In this study, we identified a total of 135 XTH genes in wheat, which were clustered into three evolutionary groups. These TaXTHs were unevenly distributed on 21 chromosomes of wheat with a majority of TaXTHs located on homelogous groups 2, 3 and 7. Gene duplication analysis revealed that segmental and tandem duplication were the main reasons for the expansion of XTH family in wheat. Interaction network predictions indicated that TaXTHs could interact with multiple proteins, including three kinases, one methyltransferase and one gibberellin-regulated protein. The promoters of the TaXTH genes harbored various cis-acting elements related to stress and hormone responses. RNA-seq data analyses showed that some TaXTH genes were induced by salt and drought stresses. Furthermore, we verified that TaXTH17 was induced by abiotic stresses and phytohormone treatments, and demonstrated that TaXTH17 was localized in the secretory pathway and cell wall. Functional analyses conducted in heterologous expression systems and in wheat established that TaXTH17 plays a negative role in plant resistance to salt and drought.

CONCLUSIONS

We identified 135 XTH genes in wheat and conducted comprehensive analyses of their phylogenetic relationships, gene structures, conserved motifs, gene duplication events, chromosome locations, interaction networks, cis-acting elements and gene expression patterns. Furthermore, we provided solid evidence supporting the notion that TaXTH17 plays a negative role in plant resistance to salt and drought stresses. Collectively, our results provide valuable insights into understanding wheat XTHs, particularly their involvement in plant stress responses, and establish a foundation for further functional and mechanistic studies of TaXTHs.

LHP1-mediated epigenetic buffering of subgenome diversity and defense responses confers genome plasticity and adaptability in allopolyploid wheat

Polyploidization is a major driver of genome diversification and environmental adaptation. However, the merger of different genomes may result in genomic conflicts, raising a major question regarding how genetic diversity is interpreted and regulated to enable environmental plasticity. By analyzing the genome-wide binding of 191 trans-factors in allopolyploid wheat, we identified like heterochromatin protein 1 (LHP1) as a master regulator of subgenome-diversified genes. Transcriptomic and epigenomic analyses of LHP1 mutants reveal its role in buffering the expression of subgenome-diversified defense genes by controlling H3K27me3 homeostasis. Stripe rust infection releases latent subgenomic variations by eliminating H3K27me3-related repression. The simultaneous inactivation of LHP1 homoeologs by CRISPR-Cas9 confers robust stripe rust resistance in wheat seedlings. The conditional repression of subgenome-diversified defenses ensures developmental plasticity to external changes, while also promoting neutral-to-non-neutral selection transitions and adaptive evolution. These findings establish an LHP1-mediated buffering system at the intersection of genotypes, environments, and phenotypes in polyploid wheat. Manipulating the epigenetic buffering capacity offers a tool to harness cryptic subgenomic variations for crop improvement.

The TaTCP4/10-B1 cascade regulates awn elongation in wheat (Triticum aestivum L.)

Awns are important morphological markers for wheat and exert a strong physiological effect on wheat yield. The awn elongation suppressor B1 has recently been cloned through association and linkage analysis in wheat. However, the mechanism of awn inhibition centered around B1 remains to be clarified. Here, we identified an allelic variant in the coding region of B1 through analysis of re-sequencing data; this variant causes an amino acid substitution and premature termination, resulting in a long-awn phenotype. Transcriptome analysis indicated that B1 inhibited awn elongation by impeding cytokinin- and auxin-promoted cell division. Moreover, B1 directly repressed the expression of TaRAE2 and TaLks2, whose orthologs have been reported to promote awn development in rice or barley. More importantly, we found that TaTCP4 and TaTCP10 synergistically inhibited the expression of B1, and a G-to-A mutation in the B1 promoter attenuated its inhibition by TaTCP4/10. Taken together, our results reveal novel mechanisms of awn development and provide genetic resources for trait improvement in wheat.

Analysis of Genetic Diversity and Phylogenetic Relationships of Wheat (Triticum aestivum L.) Genotypes Using Phenological, Molecular and DNA Barcoding Markers

Wheat (Triticum aestivum L.) is a key food crop, accounting for approximately 765 million tons produced worldwide. The present study evaluated 16 wheat genotypes using 19 morphological and phenological traits, 16 molecular markers (Inter Simple Sequence Repeats and Start Codon Targeted; ISSR and SCoT) and rbcL and matK plastid gene barcoding. The 16 wheat genotypes showed significant genetic variation using the markers assayed. Cell plot of phenological parameters revealed significant differences among the 16-day-old seedlings of wheat genotypes at Z1.1 growth stage. Collectively, W2 genotype had the lowest shoot length (SL), length of first internodes (LFI) and leaf area (LA) values, while W8 genotype had the highest diameter of first internode (DFI) and LA values. Furthermore, W7 genotype had the maximum plant biomass (PB) and leaf width (LW) values. Geometric models grouped wheat kernels into "rounded" and "nearly elongated". Estimates of heritability (H2) for these morphological characters ranged from 4.93 to 100%. The highest H2 values were recorded for root number (RN) (100%) followed by SL (88.72%), LFI (88.30%), LA (87.76%) and Feret diameter (86.68%), while the lowest H2 value was recorded for DFI (4.93%). Furthermore, highly significant genotypic and phenotypic correlations were also observed among those traits. Reproducible fingerprinting profiles and high levels of polymorphism (PPB%) of SCoT (95.46%) and ISSR (82.41%) were recorded, indicating that they are effective tools for detecting genetic variation levels among wheat genotypes. The informativeness of markers were measured through estimation of polymorphic information content (PIC), resolving power (RP) and marker index (MI). The RP and PPB% of SCoT were significantly higher compared to those of ISSR. Comparatively, the two molecular markers were effective for studying genetic diversity among wheat genotypes, but SCoT markers were more informative. Moreover, based on the two chloroplast DNA regions (rbcL and matK), MatK was found to be more reliable for differentiating among T. aestivum genotypes. Taken together, using all the studied attributes, a clear taxonomic relationship can be used to identify T. aestivum species and improve their pragmatic production and development.

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.

Tracing the Evolution of the SEPALLATA Subfamily across Angiosperms Associated with Neo- and Sub-Functionalization for Reproductive and Agronomically Relevant Traits

SEPALLATA transcription factors (SEP TFs) have been extensively studied in angiosperms as pivotal components of virtually all the MADS-box tetrameric complex master regulators of floral organ identities. However, there are published reports that suggest that some SEP members also regulate earlier reproductive events, such as inflorescence meristem determinacy and inflorescence architecture, with potential for application in breeding programs in crops. The SEP subfamily underwent a quite complex pattern of duplications during the radiation of the angiosperms. Taking advantage of the many whole genomic sequences now available, we present a revised and expanded SEP phylogeny and link it to the known functions of previously characterized genes. This snapshot supports the evidence that the major SEP3 clade is highly specialized for the specification of the three innermost floral whorls, while its sister LOFSEP clade is functionally more versatile and has been recruited for diverse roles, such as the regulation of extra-floral bract formation and inflorescence determinacy and shape. This larger pool of angiosperm SEP genes confirms previous evidence that their evolution was driven by whole-genome duplications rather than small-scale duplication events. Our work may help to identify those SEP lineages that are the best candidates for the improvement of inflorescence traits, even in far distantly related crops.

A gain-of-function allele of a DREB transcription factor gene ameliorates drought tolerance in wheat

Drought is a major environmental factor limiting wheat production worldwide. However, the genetic components underlying wheat drought tolerance are largely unknown. Here, we identify a DREB transcription factor gene (TaDTG6-B) by genome-wide association study that is tightly associated with drought tolerance in wheat. Candidate gene association analysis revealed that a 26-bp deletion in the TaDTG6-B coding region induces a gain-of-function for TaDTG6-BDel574, which exhibits stronger transcriptional activation, protein interactions, and binding activity to dehydration-responsive elements (DRE)/CRT cis-elements than the TaDTG6-BIn574 encoded by the allele lacking the deletion, thus conferring greater drought tolerance in wheat seedlings harboring this variant. Knockdown of TaDTG6-BDel574 transcripts attenuated drought tolerance in transgenic wheat, whereas its overexpression resulted in enhanced drought tolerance without accompanying phenotypic abnormalities. Furthermore, the introgression of the TaDTG6-BDel574 elite allele into drought-sensitive cultivars improved their drought tolerance, thus providing a valuable genetic resource for wheat breeding. We also identified 268 putative target genes that are directly bound and transcriptionally regulated by TaDTG6-BDel574. Further analysis showed that TaDTG6-BDel574 positively regulates TaPIF1 transcription to enhance wheat drought tolerance. These results describe the genetic basis and accompanying mechanism driving phenotypic variation in wheat drought tolerance, and provide a novel genetic resource for crop breeding programs.

Distribution, expression of hexaploid wheat Fes1s and functional characterization of two TaFes1As in Arabidopsis

Hexaploid wheat is a major food crop and is sensitive to heat stress. It is necessary to discover genes related to thermotolerance in wheat. Fes1s is a class of nucleotide exchange factor of heat shock protein 70s, proven to be participated in heat response in human, yeast, and Arabidopsis. However, little is known about Fes1s in hexaploid wheat. In this study, we identified nine Fes1s in hexaploid wheat (TaFes1s) and found that they present as three triads. A phylogenetic relationship analysis revealed that these Fes1s grouped into Fes1A, Fes1B and Fes1C subclades, and Fes1As and Fes1Bs were divergent in monocots, but possibly not in dicots. The sequences, gene structures and protein motifs of TaFes1s homoeologues within a triad were highly conserved. Through cis-elements analysis including heat shock elements, and miRNA targets prediction, we found that regulation of three TaFes1s homoeologues may be different, while the expression patterns of three homoeologues were similar. The expression levels of TaFes1As were higher than those of TaFes1Bs and TaFes1Cs, and based on these expressions, TaFes1As were chosen for functional characterization. Intriguingly, neither TaFes1A-5A nor TaFes1A-5D could not rescue the thermotolerance defect of Arabidopsis fes1a mutants at seedling stage, but in the transgenic plants seed germination was accelerated under normal and heat stress condition. The functional characterization indicated that roles of Fes1As would be different in Arabidopsis and hexaploid wheat, and function retention of TaFes1As may occur during wheat evolution. In conclusion, our study comprehensively characterized the distribution and expression of Fes1s in hexaploid wheat and found that two TaFes1As could accelerate seed germination under normal and heat stress condition.

A transcriptomic-guided strategy used in identification of a wheat rust pathogen target and modification of the target enhanced host resistance to rust pathogens

Transcriptional reprogramming is an essential feature of plant immunity and is governed by transcription factors (TFs) and co-regulatory proteins associated with discrete transcriptional complexes. On the other hand, effector proteins from pathogens have been shown to hijack these vast repertoires of plant TFs. Our current knowledge of host genes' role (including TFs) involved in pathogen colonization is based on research employing model plants such as Arabidopsis and rice with minimal efforts in wheat rust interactions. In this study, we begun the research by identifying wheat genes that benefit rust pathogens during infection and editing those genes to provide wheat with passive resistance to rust. We identified the wheat MYC4 transcription factor (TF) located on chromosome 1B (TaMYC4-1B) as a rust pathogen target. The gene was upregulated only in susceptible lines in the presence of the pathogens. Down-regulation of TaMYC4-1B using barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) in the susceptible cultivar Chinese Spring enhanced its resistance to the stem rust pathogen. Knockout of the TaMYC4-1BL in Cadenza rendered new resistance to races of stem, leaf, and stripe rust pathogens. We developed new germplasm in wheat via modifications of the wheat TaMYC4-1BL transcription factor.

Flavonoid Biosynthesis Genes in Triticum aestivum L.: Methylation Patterns in Cis-Regulatory Regions of the Duplicated CHI and F3H Genes

Flavonoids are a diverse group of secondary plant metabolites that play an important role in the regulation of plant development and protection against stressors. The biosynthesis of flavonoids occurs through the activity of several enzymes, including chalcone isomerase (CHI) and flavanone 3-hydroxylase (F3H). A functional divergence between some copies of the structural TaCHI and TaF3H genes was previously shown in the allohexaploid bread wheat Triticum aestivum L. (BBAADD genome). We hypothesized that the specific nature of TaCHI and TaF3H expression may be induced by the methylation of the promoter. It was found that the predicted position of CpG islands in the promoter regions of the analyzed genes and the actual location of methylation sites did not match. We found for the first time that differences in the methylation status could affect the expression of TaCHI copies, but not the expression of TaF3Hs. At the same time, we revealed significant differences in the structure of the promoters of only the TaF3H genes, while the TaCHI promoters were highly homologous. We assume that the promoter structure in TaF3Hs primarily affects the change in the nature of gene expression. The data obtained are important for understanding the mechanisms that regulate the synthesis of flavonoids in allopolyploid wheat and show that differences in the structure of promoters have a key effect on gene expression.

Physio-molecular traits of contrasting bread wheat genotypes associated with 15N influx exhibiting homeolog expression bias in nitrate transporter genes under different external nitrate concentrations

The high affinity nitrate transport system is a potential target for improving nitrogen use efficiency of bread wheat growing either under optimal or limiting nitrate concentration. Nitrate uptake is one of the most important traits to take into account to improve nitrogen use efficiency in wheat (Triticum aestivum L.). In this study, we aimed to gain an insight into the regulation of NO₃^- -uptake and translocation systems in two contrasting wheat genotypes [K9107(K9) vs. Choti Lerma (CL)]. Different conditions, such as NO₃^--uptake rates, soil-types, N-free solid external media, and external NO₃^- levels at the seedling stage, were considered. We also studied the contribution of homeolog expression of five genes encoding two nitrate transporters in the root tissue, along with their overall transcript expression levels relative to specific external nitrate availability. We observed that K9107 had a higher ¹⁵N influx than Choti Lerma under both limiting as well as optimum external N conditions in vermiculite-perlite (i.e., N-free solid) medium, with the improved translocation efficiency in Choti Lerma. However, in different soil types, different levels of ¹⁵N-enrichment in both the genotypes were found. Our results also demonstrated that the partitioning of dry matter in root and shoot was different under these growing conditions. Moreover, K9107 showed significantly higher relative expression of TaNRT2.1 at the lowest and TaNPF6.1 and TaNPF6.2 at the highest external nitrate concentrations. We also observed genotype-specific and nitrate starvation-dependent homeolog expression bias in all five nitrate transporter genes. Our data suggest that K9107 had a higher NO₃^- influx capacity, involving different nitrate transporters, than Choti Lerma at the seedling stage.

Genome-Wide Analysis of HSP70s in Hexaploid Wheat: Tandem Duplication, Heat Response, and Regulation

HSP70s play crucial roles in plant growth and development, as well as in stress response. Knowledge of the distribution and heat response of HSP70s is important to understand heat adaptation and facilitate thermotolerance improvement in wheat. In this study, we comprehensively analyzed the distribution of HSP70s in hexaploid wheat (TaHSP70s) and its relatives, and we found an obvious expansion of TaHSP70s in the D genome of hexaploid wheat. Meanwhile, a large portion of tandem duplication events occurred in hexaploid wheat. Among the 84 identified TaHSP70s, more than 64% were present as homeologs. The expression profiles of TaHSP70s in triads tended to be expressed more in non-stressful and heat stress conditions. Intriguingly, many TaHSP70s were especially heat responsive. Tandem duplicated TaHSP70s also participated in heat response and growth development. Further HSE analysis revealed divergent distribution of HSEs in the promoter regions of TaHSP70 homeologs, which suggested a distinct heat regulatory mechanism. Our results indicated that the heat response of TaHSP70s may experience a different regulation, and this regulation, together with the expression of tandem duplicated TaHSP70s, may help hexaploid wheat to adapt to heat conditions.

Epigenetics of wheat-rust interaction: an update

The outcome of different host-pathogen interactions is influenced by both genetic and epigenetic systems, which determine the response of plants to pathogens and vice versa. This review highlights key molecular mechanisms and conceptual advances involved in epigenetic research and the progress made in epigenetics of wheat-rust interactions. Epigenetics implies the heritable changes in the way of gene expression as a consequence of the modification of DNA bases, histone proteins, and/or non-coding-RNA biogenesis without disturbing the underlying nucleotide sequence. The changes occurring between DNA and its surrounding chromatin without altering its DNA sequence and leading to significant changes in the genome of any organism are called epigenetic changes. Epigenetics has already been used successfully to explain the mechanism of human pathogens and in the identification of pathogen-induced modifications within various host plants. Wheat rusts are one of the most vital fungal diseases throughout the major wheat-growing areas of the world. The epigenome in plant pathogens causing diseases such as wheat rusts is mysterious. The investigations of host and pathogen epigenetics in the wheat rusts system can offer a piece of suitable evidence for elucidation of the molecular basis of host-pathogen interaction. Besides, the information on the epigenetic regulation of the genes involved in resistance or pathogenicity will provide better insights into the complex resistance signaling pathways and could provide answers to certain key questions, such as whether epigenetic regulation of certain genes is imparting resistance to host in response of certain pathogen elicitors or not. In the last few years, there has been an upsurge in research on the host as well as pathogen epigenetics and its outcome in plant-pathogen interactions. This review summarizes the progress made in the areas related to the epigenetic control of host-pathogen interaction with particular emphasis on wheat rusts.

The wheat AGL6-like MADS-box gene is a master regulator for floral organ identity and a target for spikelet meristem development manipulation

The AGAMOUS-LIKE6 (AGL6)-like genes are ancient MADS-box genes and are functionally studied in a few model plants. The knowledge of these genes in wheat remains limited. Here, by studying a 'double homoeolog mutant' of the AGL6 gene in tetraploid wheat, we showed that AGL6 was required for the development of all four whorls of floral organs with dosage-dependent effect on floret fertility. Yeast two-hybrid analyses detected interactions of AGL6 with all classes of MADS-box proteins in the ABCDE model for floral organ development. AGL6 was found to interact with several additional proteins, including the G protein β and γ (DEP1) subunits. Analysis of the DEP1-B mutant showed a significant reduction in spikelet number per spike in tetraploid wheat, while overexpression of AGL6 in common wheat increased the spikelet number per spike and hence the grain number per spike. RNA-seq analysis identified the regulation of several meristem activity genes by AGL6, such as FUL2 and TaMADS55. Our work therefore extensively updated the wheat ABCDE model and proposed an alternative approach to improve wheat grain yield by manipulating the AGL6 gene.

Alternative Splicing of TaGS3 Differentially Regulates Grain Weight and Size in Bread Wheat

The heterotrimeric G-protein mediates growth and development by perceiving and transmitting signals in multiple organisms. Alternative splicing (AS), a vital process for regulating gene expression at the post-transcriptional level, plays a significant role in plant adaptation and evolution. Here, we identified five splicing variants of G_γ subunit gene TaGS3 (TaGS3.1 to TaGS3.5), which showed expression divergence during wheat polyploidization, and differential function in grain weight and size determination. TaGS3.1 overexpression significantly reduced grain weight by 5.89% and grain length by 5.04%, while TaGS3.2–3.4 overexpression did not significantly alter grain size compared to wild type. Overexpressing TaGS3.5 significantly increased the grain weight by 5.70% and grain length by 4.30%. Biochemical assays revealed that TaGS3 isoforms (TaGS3.1–3.4) with an intact OSR domain interact with WGB1 to form active G_βγ heterodimers that further interact with WGA1 to form inactive G_αβγ heterotrimers. Truncated isoforms TaGS3.2–3.4 , which lack the C-terminal Cys-rich region but have enhanced binding affinity to WGB1, antagonistically compete with TaGS3.1 to bind WGB1, while TaGS3.5 with an incomplete OSR domain does not interact with WGB1. Taking these observations together, we proposed that TaGS3 differentially regulates grain size via AS, providing a strategy by which the grain size is fine-tuned and regulated at the post-transcriptional level.

Evolution of the MLO gene families in octoploid strawberry (Fragaria ×ananassa) and progenitor diploid species identified potential genes for strawberry powdery mildew resistance

Powdery mildew (PM) caused by Podosphaera aphanis is a major fungal disease of cultivated strawberry. Mildew Resistance Locus O (MLO) is a gene family described for having conserved seven-transmembrane domains. Induced loss-of-function in specific MLO genes can confer durable and broad resistance against PM pathogens. However, the genomic structure and potential role of MLO genes for PM resistance have not been characterized yet in the octoploid cultivated strawberry. In the present study, MLO gene families were characterized in four diploid progenitor species (Fragaria vesca, F. iinumae, F. viridis, and F. nipponica) and octoploid cultivated (Fragaria ×ananassa) strawberry, and potential sources of MLO-mediated susceptibility were identified. Twenty MLO sequences were identified in F. vesca and 68 identified in F. ×ananassa. Phylogenetic analysis divided diploid and octoploid strawberry MLO genes into eight different clades, in which three FveMLO (MLO10, MLO17, and MLO20) and their twelve orthologs of FaMLO were grouped together with functionally characterized MLO genes conferring PM susceptibility. Copy number variations revealed differences in MLO composition among homoeologous chromosomes, supporting the distinct origin of each subgenome during the evolution of octoploid strawberry. Dissecting genomic sequence and structural variations in candidate FaMLO genes revealed their potential role associated with genetic controls and functionality in strawberry against PM pathogen. Furthermore, the gene expression profiling and RNAi silencing of putative FaMLO genes in response to the pathogen indicate the function in PM resistance. These results are a critical first step in understanding the function of strawberry MLO genes and will facilitate further genetic studies of PM resistance in cultivated strawberry.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Evolution of Homeologous Gene Expression in Polyploid Wheat

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

TaANR1-TaBG1 and TaWabi5-TaNRT2s/NARs Link ABA Metabolism and Nitrate Acquisition in Wheat Roots

Nitrate is the preferred form of nitrogen for most plants, acting both as a nutrient and a signaling molecule. However, the components and regulatory factors governing nitrate uptake in bread wheat (Triticum aestivum), one of the world's most important crop species, have remained unclear, largely due to the complexity of its hexaploid genome. Here, based on recently released whole-genome information for bread wheat, the high-affinity nitrate transporter2 (NRT2) and the nitrate-assimilation-related (NAR) gene family are characterized. We show that abscisic acid (ABA)- Glc ester deconjugation is stimulated in bread wheat roots by nitrate resupply following nitrate withdrawal, leading to enhanced root-tissue ABA accumulation, and that this enhancement, in turn, affects the expression of root-type NRT2/NAR genes. TaANR1 is shown to regulate nitrate-mediated ABA accumulation by directly activating TaBG1, while TaWabi5 is involved in ABA-mediated NO3 - induction of NRT2/NAR genes. Building on previous evidence establishing ABA involvement in the developmental response to high-nitrate stress, our study suggests that ABA also contributes to the optimization of nitrate uptake by regulating the expression of NRT2/NAR genes under limited nitrate supply, offering a new target for improvement of nitrate absorption in crops.

BnaMPK6 is a determinant of quantitative disease resistance against Sclerotinia sclerotiorum in oilseed rape

Sclerotinia sclerotiorum causes a devastating disease in oilseed rape (Brassica napus), resulting in major economic losses. Resistance response of B. napus against S. sclerotiorum exhibits a typical quantitative disease resistance (QDR) characteristic, but the molecular determinants of this QDR are largely unknown. In this study, we isolated a B. napus mitogen-activated protein kinase gene, BnaMPK6, and found that BnaMPK6 expression is highly responsive to infection by S. sclerotiorum and treatment with salicylic acid (SA) or jasmonic acid (JA). Moreover, overexpression (OE) of BnaMPK6 significantly enhances resistance to S. sclerotiorum, whereas RNAi in BnaMPK6 significantly reduces this resistance. These results showed that BnaMPK6 plays an important role in defense to S. sclerotiorum. Furthermore, expression of defense genes associated with SA-, JA- and ethylene (ET)-mediated signaling was investigated in BnaMPK6-RNAi, WT and BnaMPK6-OE plants after S. sclerotiorum infection, and consequently, it was indicated that the activation of ET signaling by BnaMPK6 may play a role in the defense. Further, four BnaMPK6-encoding homologous loci were mapped in the B. napus genome. Using the allele analysis and expression analysis on the four loci, we demonstrated that the locus BnaA03.MPK6 makes an important contribution to QDR against S. sclerotiorum. Our data indicated that BnaMPK6 is a previously unknown determinant of QDR against S. sclerotiorum in B. napus.

Gene Expression Changes During the Allo-/Deallopolyploidization Process of Brassica napus

Gene expression changes due to allopolyploidization have been extensively studied in plants over the past few decades. Nearly all these studies focused on comparing the changes before and after genome merger. In this study, we used the uniquely restituted Brassica rapa (RBR, A_eA_e, 2n = 20) obtained from Brassica napus (A_nA_nC_nC_n, 2n = 38) to analyze the gene expression changes and its potential mechanism during the process of allo-/deallopolyploidization. RNA-seq-based transcriptome profiling identified a large number of differentially expressed genes (DEGs) between RBR and natural B. rapa (A_rA_r), suggesting potential effects of allopolyploidization/domestication of AA component of B. napus at the tetrapolyploid level. Meanwhile, it was revealed that up to 20% of gene expressions were immediately altered when compared with those in the A_n-subgenome. Interestingly, one fifth of these changes are in fact indicative of the recovery of antecedent gene expression alternations occurring since the origin of B. napus and showed association with homoeologous expression bias between A_n and C_n subgenomes. Enrichment of distinct gene ontology (GO) categories of the above sets of genes further indicated potential functional cooperation of the A_n and C_n subgenome of B. napus. Whole genome methylation analysis revealed a small number of DEGs were identified in the differentially methylated regions.

Conservation and divergence of ancestral AGAMOUS/SEEDSTICK subfamily genes from the basal angiosperm Magnolia wufengensis

AGAMOUS/SEEDSTICK (AG/STK) subfamily genes play crucial roles in the reproductive development of plants. However, most of our current knowledge of AG/STK subfamily genes is restricted to core eudicots and grasses, and the knowledge of ancestral exon-intron structures, expression patterns, protein-protein interaction patterns and functions of AG/STK subfamily genes remains unclear. To determine these, we isolated AG/STK subfamily genes (MawuAG1, MawuAG2 and MawuSTK) from a woody basal angiosperm Magnolia wufengensis (Magnoliaceae). MawuSTK arose from the gene duplication event occurring before the diversification of extant angiosperms, and MawuAG1 and MawuAG2 may result from a gene duplication event occurring before the divergence of Magnoliaceae and Lauraceae. Gene duplication led to apparent diversification in their expression and interaction patterns. It revealed that expression in both stamens and carpels likely represents the ancestral expression profiles of AG lineage genes, and expression of STK-like genes in stamens may have been lost soon after the appearance of the STK lineage. Moreover, AG/STK subfamily proteins may have immediately established interactions with the SEPALLATA (SEP) subfamily proteins following the emergence of the SEP subfamily; however, their interactions with the APETALA1/FRUITFULL subfamily proteins or themselves differ from those found in monocots and basal and core eudicots. MawuAG1 plays highly conserved roles in the determinacy of stamen, carpel and ovule identity, while gene duplication contributed to the functional diversification of MawuAG2 and MawuSTK. In addition, we investigated the evolutionary history of exon-intron structural changes of the AG/STK subfamily, and a novel splice-acceptor mode (GUU-AU) and the convergent evolution of N-terminal extension in the euAG and PLE subclades were revealed for the first time. These results further advance our understanding of ancestral AG/STK subfamily genes in terms of phylogeny, exon-intron structures, expression and interaction patterns, and functions, and provide strong evidence for the significance of gene duplication in the expansion and evolution of the AG/STK subfamily.

Genome-wide analysis of MIKC-type MADS-box genes in wheat: pervasive duplications, functional conservation and putative neofunctionalization

Wheat (Triticum aestivum) is one of the most important crops worldwide. Given a growing global population coupled with increasingly challenging cultivation conditions, facilitating wheat breeding by fine-tuning important traits is of great importance. MADS-box genes are prime candidates for this, as they are involved in virtually all aspects of plant development. Here, we present a detailed overview of phylogeny and expression of 201 wheat MIKC-type MADS-box genes. Homoeolog retention is significantly above the average genome-wide retention rate for wheat genes, indicating that many MIKC-type homoeologs are functionally important and not redundant. Gene expression is generally in agreement with the expected subfamily-specific expression pattern, indicating broad conservation of function of MIKC-type genes during wheat evolution. We also found extensive expansion of some MIKC-type subfamilies, especially those potentially involved in adaptation to different environmental conditions like flowering time genes. Duplications are especially prominent in distal telomeric regions. A number of MIKC-type genes show novel expression patterns and respond, for example, to biotic stress, pointing towards neofunctionalization. We speculate that conserved, duplicated and neofunctionalized MIKC-type genes may have played an important role in the adaptation of wheat to a diversity of conditions, hence contributing to the importance of wheat as a global staple food.

Molecular Characterization of Auxin Efflux Carrier- ABCB1 in hexaploid wheat

Auxin is an important phytohormone that regulates response, differentiation, and development of plant cell, tissue, and organs. Along with its local production, long-distance transport coordinated by the efflux/influx membrane transporters is instrumental in plant development and architecture. In the present study, we cloned and characterized a wheat (Triticum aestivum) auxin efflux carrier ABCB1. The TaABCB1 was physically localized to the proximal 15% of the short arm of wheat homoeologous group 7 chromosomes. Size of the Chinese spring (CS) homoeologs genomic copies ranged from 5.3–6.2 kb with the 7A copy being the largest due to novel insertions in its third intron. The three homoeologous copies share 95–97% sequence similarity at the nucleotide, 98–99% amino acid, and overall Q-score of 0.98 at 3-D structure level. Though detected in all analyzed tissues, TaABCB1 predominantly expressed in the meristematic tissues likely due to the presence of meristem-specific activation regulatory element identified in the promoter region. RNAi plants of TaABCB1 gene resulted in reduced plant height and increased seed width. Promoter analysis revealed several responsive elements detected in the promoter region including that for different hormones as auxin, gibberellic acid, jasmonic acid and abscisic acid, light, and circadian regulated elements.

Overexpression of paralogues of the wheat expansin gene TaEXPA8 improves low-temperature tolerance in Arabidopsis

Low temperature is one of the important factors limiting wheat yield in cold regions. Expansins are nonenzymatic proteins that loosen cell walls and play important roles in diverse biological processes related to cell wall modification, including development and stress tolerance. Many studies have shown that expansins are involved in resistance to various abiotic stresses, such as heat and drought. However, the role of expansins in response to low-temperature stress remains unclear. Based on our previous transcriptome data of a winter wheat cultivar Dongnongdongmai 2 (DN2), we found that one of the expansin genes, TaEXPA8, was significantly induced by low temperature, indicating a role for TaEXPA8 in cold resistance. In this study, the paralogous TaEXPA8 genes TaEXPA8-A, TaEXPA8-B and TaEXPA8-D were cloned by RT-PCR. These three genes were then transformed into Arabidopsis by the floral dip method. Expression patterns of TaEXPA8 genes in different tissues and in response to several abiotic stresses and hormones were detected by quantitative real-time PCR (qRT-PCR). The results showed that TaEXPA8-A and TaEXPA8-B were expressed mainly in roots, while TaEXPA8-D was expressed predominantly in flowers. TaEXPA8 genes were induced by low-temperature and drought. The overexpression of TaEXPA8-B and TaEXPA8-D enhanced low-temperature resistance and had increased superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activity and soluble protein, MDA and proline content. In summary, our study suggested that the expansins TaEXPA8-B and TaEXPA8-D are involved in the response to low temperature and possibly play a role in cold resistance by activating the protective enzyme system.

Overexpression of NtabDOG1L promotes plant growth and enhances drought tolerance in Nicotiana tabacum

Drought is one of the major environmental stresses limiting crop growth and production. It is very important to exploit and utilize drought-tolerance genes to improve crop drought-resistance. In this study, we identified two homoeologs of a Nicotiana tabacum (Ntab) DELAY OF GERMINATION (DOG) 1 like gene, named as NtabDOG1L-T and NtabDOG1L-S, respectively. The NtabDOG1L genes were preferentially expressed in roots and their expression levels were induced by polyethylene glycol, high salt, cold, and abscisic acid treatments. Subcellular localization results indicated that NtabDOG1L-T was localized in the nucleus, cytoplasm and cell membrane. Overexpression of NtabDOG1L-T in tobacco resulted in roots growth enhancement in transgenic plants. Furthermore, overexpression of NtabDOG1L-T enhanced drought stress tolerance in transgenic tobacco. The transgenic tobacco lines exhibited lower leaf water loss and electrolyte leakage, lower content of malondialdehyde and reactive oxygen species (ROS), and higher antioxidant enzymes activities after drought treatment when compared with wild type (WT) plants. In addition, the expression levels of several genes encoding key antioxidant enzymes and drought-related proteins were higher in the transgenic plants than in the WT plants under drought stress. Taken together, our results showed that NtabDOG1L functions as a novel regulator that improves plant growth and drought tolerance in tobacco.

Genetic and Molecular Control of Floral Organ Identity in Cereals

Grasses represent a major family of monocots comprising mostly cereals. When compared to their eudicot counterparts, cereals show a remarkable morphological diversity. Understanding the molecular basis of floral organ identity and inflorescence development is crucial to gain insight into the grain development for yield improvement purposes in cereals, however, the exact genetic mechanism of floral organogenesis remains elusive due to their complex inflorescence architecture. Extensive molecular analyses of Arabidopsis and other plant genera and species have established the ABCDE floral organ identity model. According to this model, hierarchical combinatorial activities of A, B, C, D, and E classes of homeotic genes regulate the identity of different floral organs with partial conservation and partial diversification between eudicots and cereals. Here, we review the developmental role of A, B, C, D, and E gene classes and explore the recent advances in understanding the floral development and subsequent organ specification in major cereals with reference to model plants. Furthermore, we discuss the evolutionary relationships among known floral organ identity genes. This comparative overview of floral developmental genes and associated regulatory factors, within and between species, will provide a thorough understanding of underlying complex genetic and molecular control of flower development and floral organ identity, which can be helpful to devise innovative strategies for grain yield improvement in cereals.

Genome-wide investigation of the MADS gene family and dehulling genes in tartary buckwheat (Fagopyrum tataricum)

Genome-wide identification, expression analysis and potential functional characterization of previously uncharacterized MADS family of tartary buckwheat, emphasized the importance of this gene family in plant growth and development. The MADS transcription factor is a key regulatory factor in the development of most plants. The MADS gene in plants controls all aspects of tissue and organ growth and reproduction and can be used to regulate plant seed cracking. However, there has been little research on the MADS genes of tartary buckwheat (Fagopyrum tataricum), which is an important edible and medicinal crop. The recently published whole genome sequence of tartary buckwheat allows us to study the tissue and expression profiles of the MADS gene in tartary buckwheat at a genome-wide level. In this study, 65 MADS genes of tartary buckwheat were identified and renamed according to the chromosomal distribution of the FtMADS genes. Here, we provide a complete overview of the gene structure, gene expression, genomic mapping, protein motif organization, and phylogenetic relationships of each member of the gene family. According to the phylogenetic relationship of MADS genes, the transcription factor family was divided into two subfamilies, the M subfamily (28 genes) and the MIKC subfamily (37 genes). The results showed that the FtMADS genes belonged to related sister pairs and the chromosomal map showed that the replication of FtMADSs was related to the replication of chromosome blocks. In different tissues and at different fruit development stages, the FtMADS genes obtained by real-time quantitative PCR (RT-qPCR) showed obvious expression patterns. A comprehensive analysis of the MADS genes in tartary buckwheat was conducted. Through systematic analysis, the potential genes that may regulate the growth and development of tartary buckwheat and the genes that may regulate the easy dehulling of tartary buckwheat fruit were screened, which laid a solid foundation for improving the quality of tartary buckwheat.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The genome-wide transcriptional consequences of the nullisomic-tetrasomic stocks for homoeologous group 7 in bread wheat

Background

Hexaploid bread wheat (Triticum aestivum L) arose by two polyploidisation events from three diploid species with homoeologous genomes. Nullisomic-tetrasomic (nulli-tetra or NT) lines are aneuploid wheat plants lacking two and adding two of six homoeologous chromosomes. These plants can grow normally, but with significantly morphological variations because the adding two chromosomes or the remaining four chromosomes compensate for those absent. Despite these interesting phenomena, detailed molecular mechanisms underlying dosage deletion and compensation in these useful genetic materials have not been determined.

Results

By sequencing the transcriptomes of leaves in two-week-old seedlings, we showed that the profiles of differentially expressed genes between NT stocks for homoeologous group 7 and the parent hexaploid Chinese Spring (CS) occurred throughout the whole genome with a subgenome and chromosome preference. The deletion effect of nulli-chromosomes was compensated partly by the tetra-chromosomes via the dose level of expressed genes, according to the types of homoeologous genes. The functions of differentially regulated genes primarily focused on carbon metabolic process, photosynthesis process, hormone metabolism, and responding to stimulus, and etc., which might be related to the defective phenotypes that included reductions in plant height, flag leaf length, spikelet number, and kernels per spike.

Conclusions

The perturbation of the expression levels of transcriptional genes among the NT stocks for homoeologous group 7 demonstrated the gene dosage effect of the subgenome at the genome-wide level. The gene dosage deletion and compensation can be used as a model to elucidate the functions of the subgenomes in modern polyploid plants.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5421-3) contains supplementary material, which is available to authorized users.

Genome-Wide Identification and Expression Profiling of the TCP Family Genes in Spike and Grain Development of Wheat (Triticum aestivum L.)

The TCP family genes are plant-specific transcription factors and play important roles in plant development. TCPs have been evolutionarily and functionally studied in several plants. Although common wheat (Triticum aestivum L.) is a major staple crop worldwide, no systematic analysis of TCPs in this important crop has been conducted. Here, we performed a genome-wide survey in wheat and found 66 TCP genes that belonged to 22 homoeologous groups. We then mapped these genes on wheat chromosomes and found that several TCP genes were duplicated in wheat including the ortholog of the maize TEOSINTE BRANCHED 1. Expression study using both RT-PCR and in situ hybridization assay showed that most wheat TCP genes were expressed throughout development of young spike and immature seed. Cis-acting element survey along promoter regions suggests that subfunctionalization may have occurred for homoeologous genes. Moreover, protein-protein interaction experiments of three TCP proteins showed that they can form either homodimers or heterodimers. Finally, we characterized two TaTCP9 mutants from tetraploid wheat. Each of these two mutant lines contained a premature stop codon in the A subgenome homoeolog that was dominantly expressed over the B subgenome homoeolog. We observed that mutation caused increased spike and grain lengths. Together, our analysis of the wheat TCP gene family provides a start point for further functional study of these important transcription factors in wheat.

Overexpression of Lilium formosanumMADS-box (LFMADS) Causing Floral Defects While Promoting Flowering in Arabidopsis thaliana, Whereas Only Affecting Floral Transition Time in Nicotiana tabacum

The Formosa lily (Lilium formosanum) is one of the most common horticultural species in Taiwan. To explore gene regulation involved in this species, we used transcriptome analysis to generate PH-FB (mixed floral buds) and PH-LF (mature leaves) datasets. Combination of the PH-FB and PH-LF constructed a de novo assembly of the ALL dataset, including 18,041 contigs and 23,807 unigenes by Nr, GO, COG, and KEGG databases. The differential gene expression (DGE) analysis revealed 9937 genes were upregulated while 10,383 genes were downregulated in the developing floral buds compared to mature leaves. Seven putative genes (LFMADS1 to 7) encoding floral organ identity proteins were selected for further analysis. LFMADS1-6 genes were specifically expressed in the floral organ, while LFMADS7 in the floral buds and mature leaves. Phylogenetic analysis revealed that LFMADS1-3 is classified into B-class, LFMADS4 into C-class, LFMADS5 into D-class, and LFMADS6-7 into E-class, respectively. LFMADS-GFP fusion proteins appeared to localize in the nucleus, supporting their roles as transcription factors (TFs). Overexpression of the LFMADS2, LFMADS4, and LFMADS6 genes in Arabidopsis resulted in early flowering and floral defect, however, only early flowering in transgenic tobacco was observed. Highly expressed floral integrator genes, including AtFT, AtLFY, and AtFUL in transgenic Arabidopsis and NtFUL and NtSOC1 in transgenic tobacco, resulted in early flowering phenotype through qRT-PCR analysis. Yeast two-hybrid analysis suggested that LFMADSs may form higher order complexes with the B-, C-, D, and/or E-class proteins to determine the floral organ identity. Furthermore, E-class LFMADS proteins may function as a glue to mediate and strengthen the protein-protein interactions. Therefore, our de novo datasets would provide information for investigating other differentially expressed candidate transcripts. In addition, functional conservation of LFMADSs appears to be vital in floral transition and floral organ identity.

Dissecting the role of MADS-box genes in monocot floral development and diversity

Many monocot plants have high social and economic value. These include grasses such as rice (Oryza sativa), wheat (Triticum aestivum), and barley (Hordeum vulgare), which produce soft commodities for many food and beverage industries, and ornamental flowers such ase lily (Lilium longiflorum) and orchid (Oncidium Gower Ramsey), which represent an important component of international flower markets. There is constant pressure to improve the development and diversity of these species, with a significant emphasis on flower development, and this is particularly relevant considering the impact of changing environments on reproduction and thus yield. MADS-box proteins are a family of transcription factors that contain a conserved 60 amino acid MADS-box motif. In plants, attention has been devoted to characterization of this family due to their roles in inflorescence and flower development, which holds promise for the modification of floral architecture for plant breeding. This has been explored in diverse angiosperms, but particularly the dicot model Arabidopsis thaliana. The focus of this review is on the less well characterized roles of the MADS-box proteins in monocot flower development and how changes in MADS-box proteins throughout evolution may have contributed to creating a diverse range of flowers. Examining these changes within the monocots can identify the importance of certain genes and pinpoint those which might be useful in future crop improvement and breeding strategies.

Epigenetic perspectives on the evolution and domestication of polyploid plant and crops

Polyploidy or whole genome duplication (WGD) is a prominent feature for genome evolution of some animals and all flowering plants, including many important crops such as wheat, cotton, and canola. In autopolyploids, genome duplication often perturbs dosage regulation on biological networks. In allopolyploids, interspecific hybridization could induce genetic and epigenetic changes, the effects of which could be amplified by genome doubling (ploidy changes). Albeit the importance of genetic changes, some epigenetic changes can be stabilized and transmitted as epialleles into the progeny, which are subject to natural selection, adaptation, and domestication. Here we review recent advances for general and specific roles of epigenetic changes in the evolution of flowering plants and domestication of agricultural crops.

Transcriptome analysis reveals potential mechanisms for different grain size between natural and resynthesized allohexaploid wheats with near-identical AABB genomes

BACKGROUND

Common wheat is a typical allohexaploid species (AABBDD) derived from the interspecific crossing between allotetraploid wheat (AABB) and Aegilops tauschii (DD). Wide variation in grain size and shape observed among Aegilops tauschii can be retained in synthetic allohexaploid wheats, but the underlying mechanism remains enigmatic. Here, the natural and resynthesized allohexaploid wheats with near-identical AB genomes and different D genomes (TAA10 and XX329) were employed for analysis.

RESULTS

Significant differences in grain size and weight between TAA10 and XX329 were observed at the early stages of development, which could be mainly attributed to the higher growth rates of the pericarp and endosperm cells in XX329 compared to TAA10. Furthermore, comparative transcriptome analysis identified that 8891 of 69,711 unigenes (12.75%) were differentially expressed between grains at 6 days after pollination (DAP) of TAA10 and XX329, including 5314 up-regulated and 3577 down-regulated genes in XX329 compared to TAA10. The MapMan functional annotation and enrichment analysis revealed that the differentially expressed genes were significantly enriched in categories of cell wall, carbohydrate and hormone metabolism. Notably, consistent with the up-regulation of sucrose synthase genes in resynthesized relative to natural allohexaploid wheat, the resynthesized allohexaploid wheat accumulated much higher contents of glucose and fructose in 6-DAP grains than those of the natural allohexaploid wheat.

CONCLUSIONS

These data indicated that the genetic variation of the D genome induced drastic alterations of gene expression in grains of the natural and resynthesized allohexaploid wheats, which may contribute to the observed differences in grain size and weight.

The Wheat MYB Transcription Factor TaMYB31 Is Involved in Drought Stress Responses in Arabidopsis

Drought is one of the major environmental stresses limiting crop growth and production. MYB family transcription factors play crucial roles in response to abiotic stresses. Previous studies found that TaMYB31 is transcriptionally induced by drought stress. However, the biological functions of TaMYB31 in drought stress responses remained unknown. In this study, three TaMYB31 homoeologous genes from hexaploid wheat, designated TaMYB31-A, TaMYB31-B, and TaMYB31-D, were cloned and characterized. Expression analysis showed that TaMYB31 genes have different tissue expression patterns, and TaMYB31-B has relatively high expression levels in most tested tissues. All the three homoeologs were up-regulated by polyethylene glycol (PEG) 6000 and abscisic acid (ABA) treatments. Subcellular localization analyses revealed that TaMYB31 is localized to the nucleus. Ectopic expression of the TaMYB31-B gene in Arabidopsis affected plants growth and enhanced drought tolerance. In addition, seed germination and seedling root growth of TaMYB31-B transgenic plants were more sensitive to exogenous ABA treatment compared to wild type control. RNA-seq analysis indicated that TaMYB31 functions through up-regulation of wax biosynthesis genes and drought-responsive genes. These results provide evidence that TaMYB31 acts as a positive regulator of drought resistance, and justify its potential application in genetic modification of crop drought tolerance.

Homoeologous exchanges cause extensive dosage-dependent gene expression changes in an allopolyploid crop

Structural variation is a major source of genetic diversity and an important substrate for selection. In allopolyploids, homoeologous exchanges (i.e. between the constituent subgenomes) are a very frequent type of structural variant. However, their direct impact on gene content and gene expression had not been determined. Here, we used a tissue-specific mRNA-Seq dataset to measure the consequences of homoeologous exchanges (HE) on gene expression in Brassica napus, a representative allotetraploid crop. We demonstrate that expression changes are proportional to the change in gene copy number triggered by the HEs. Thus, when homoeologous gene pairs have unbalanced transcriptional contributions before the HE, duplication of one copy does not accurately compensate for loss of the other and combined homoeologue expression also changes. These effects are, however, mitigated over time. This study sheds light on the origins, timing and functional consequences of homeologous exchanges in allopolyploids. It demonstrates that the interplay between new structural variation and the resulting impacts on gene expression, influences allopolyploid genome evolution.

BnDGAT1s Function Similarly in Oil Deposition and Are Expressed with Uniform Patterns in Tissues of Brassica napus

As an allotetraploid oilcrop, Brassica napus contains four duplicated Acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) genes, which catalyze one of the rate-limiting steps in triacylglycerol (TAG) biosynthesis in plants. While all four BnDGAT1s have been expressed functionally in yeast, their expression patterns in different germplasms and tissues and also consequent contribution to seed oil accumulation in planta remain to be elucidated. In this study, the coding regions of the four BnDGAT1s were expressed in an Arabidopsis dgat1 mutant. All four BnDGAT1s showed similar effects on oil content and fatty acid composition, a result which is different from that observed in previous studies of their expression in yeast. Expression patterns of BnDGAT1s were analyzed in developing seeds of 34 B. napus inbred lines and in different tissues of 14 lines. Different expression patterns were observed for the four BnDGAT1s, which suggests that they express independently or randomly in different germplasm sources. Higher expression of BnDGAT1s was correlated with higher seed oil content lines. Tissue-specific analyses showed that the BnDGAT1s were expressed in a uniform pattern in different tissues. Our results suggest that it is important to maintain expression of the four BnDGAT1s for maximum return on oil content.

Features of Ppd-B1 expression regulation and their impact on the flowering time of wheat near-isogenic lines

BACKGROUND

Photoperiod insensitive Ppd-1a alleles determine early flowering of wheat. Increased expression of homoeologous Ppd-D1a and Ppd-A1a result from deletions in the promoter region, and elevated expression of Ppd-B1a is determined by an increased copy number.

RESULTS

In this study, using bread wheat cultivars Sonora and PSL2, which contrast in flowering time, and near-isogenic lines resulting from their cross, "Ppd-m" and "Ppd-w" with Ppd-B1a introgressed from Sonora, we investigated the putative factors that influence Ppd-B1a expression. By analyzing the Ppd-B1a three distinct copies, we identified an indel and the two SNPs, which distinguished the investigated allele from other alleles with a copy number variation. We studied the expression of the Ppd-A1, Ppd-B1a, and Ppd-D1 genes along with genes that are involved in light perception (PhyA, PhyB, PhyC) and the flowering initiation (Vrn-1, TaFT1) and discussed their interactions. Expression of Ppd-B1a in the "Ppd-m" line, which flowered four days earlier than "Ppd-w", was significantly higher. We found PhyC to be up-regulated in lines with Ppd-B1a alleles. Expression of PhyC was higher in "Ppd-m". Microsatellite genotyping demonstrated that in the line "Ppd-m", there is an introgression in the pericentromeric region of chromosome 5B from the early flowering parental Sonora, while the "Ppd-w" does not have this introgression. FHY3/FAR1 is known to be located in this region. Expression of the transcription factor FHY3/FAR1 was higher in the "Ppd-m" line than in "Ppd-w", suggesting that FHY3/FAR1 is important for the wheat flowering time and may cause earlier flowering of "Ppd-m" as compared to "Ppd-w".

CONCLUSIONS

We propose that there is a positive bidirectional regulation of Ppd-B1a and PhyC with an FHY3/FAR1 contribution. The bidirectional regulation can be proposed for Ppd-A1a and Ppd-D1a. Using in silico analysis, we demonstrated that the specificity of the Ppd-B1 regulation compared to that of homoeologous genes involves not only a copy number variation but also distinct regulatory elements.

Altered expression of the TaRSL2 gene contributed to variation in root hair length during allopolyploid wheat evolution

Altered expression of the TaRSL2 gene was positively correlated with variation in root hair length during allopolyploid wheat evolution, and overexpression of TaRSL2 in Arabidopsis increases root hair length. Root hairs aid nutrient and water uptake and anchor the plant in the soil. Allopolyploid wheats display significant growth vigor in terms of root hair length compared to their diploid progenitors, but little is known about the molecular basis of variation in root hair length during wheat allopolyploidization. Here, we isolated three orthologs of the Arabidopsis root hair gene ROOT HAIR DEFECTIVE SIX-LIKE 2 (AtRSL2) in allohexaploid wheat, designated TaRSL2-4A, TaRSL2-4B and TaRSL2-4D. The deduced polypeptides of these three TaRSL2 homoeologous genes shared high similarity, and a conserved basic helix-loop-helix (bHLH) domain was present in their C-terminal regions. Notably, the expression of TaRSL2 was positively correlated with root hair length of wheat accessions with different ploidy levels. Moreover, ectopic overexpression of TaRSL2-4D in Arabidopsis could increase root hair length. We found that the transcript levels of TaRSL2 homoeologous genes dynamically changed during allopolyploid wheat evolution, implicating the complexity of the underlying molecular mechanism. Collectively, we propose that altered expression of the TaRSL2 gene contributed to variation in root hair length in allopolyploid wheats.

Down-regulation of BnDA1, whose gene locus is associated with the seeds weight, improves the seeds weight and organ size in Brassica napus

Brassica napus L. is an important oil crop worldwide and is the main raw material for biofuel. Seed weight and seed size are the main contributors to seed yield. DA1 (DA means big in Chinese) is an ubiquitin receptor and negatively regulates seed size. Down-regulation of AtDA1 in Arabidopsis leads to larger seeds and organs by increasing cell proliferation in integuments. In this study, BnDA1 was down-regulated in B. napus by over expressed of AtDA1^R358K , which is a functional deficiency of DA1 with an arginine-to-lysine mutation at the 358th amino acid. The results showed that the biomass and size of the seeds, cotyledons, leaves, flowers and siliques of transgenic plants all increased significantly. In particular, the 1000 seed weight increased 21.23% and the seed yield per plant increased 13.22% in field condition. The transgenic plants had no negative traits related to yield. The candidate gene association analysis demonstrated that the BnDA1 locus was contributed to the seeds weight. Therefore, our study showed that regulation of DA1 in B. napus can increase the seed yield and biomass, and DA1 is a promising target for crop improvement.

Chromatin modification contributes to the expression divergence of three TaGS2 homoeologs in hexaploid wheat

Plastic glutamine synthetase (GS2) is responsible for ammonium assimilation. The reason that TaGS2 homoeologs in hexaploid wheat experience different selection pressures in the breeding process remains unclear. TaGS2 were minimally expressed in roots but predominantly expressed in leaves, and TaGS2-B had higher expression than TaGS2-A and TaGS2-D. ChIP assays revealed that the activation of TaGS2-B expression in leaves was correlated with increased H3K4 trimethylation. The transcriptional silencing of TaGS2 in roots was correlated with greater cytosine methylation and less H3K4 trimethylation. Micrococcal nuclease and DNase I accessibility experiments indicated that the promoter region was more resistant to digestion in roots than leaves, which indicated that the closed nucleosome conformation of the promoter region was important to the transcription initiation for the spatial-temporal expression of TaGS2. In contrast, the transcribed regions possess different nuclease accessibilities of three TaGS2 homoeologs in the same tissue, suggesting that nucleosome conformation of the transcribed region was part of the fine adjustment of TaGS2 homoeologs. This study provides evidence that histone modification, DNA methylation and nuclease accessibility coordinated the control of the transcription of TaGS2 homoeologs. Our results provided important evidence that TaGS2-B experienced the strongest selection pressures during the breeding process.

Molecular and Functional Characterization of Wheat ARGOS Genes Influencing Plant Growth and Stress Tolerance

Auxin Regulated Gene involved in Organ Size (ARGOS) is significantly and positively associated with organ size and is involved in abiotic stress responses in plants. However, no studies on wheat ARGOS genes have been reported to date. In the present study, three TaARGOS homoeologous genes were isolated and located on chromosomes 4A, 4B, and 4D of bread wheat, all of which are highly conserved in wheat and its wild relatives. Comparisons of gene expression in different tissues demonstrated that the TaARGOSs were mainly expressed in the stem. Furthermore, the TaARGOS transcripts were significantly induced by drought, salinity, and various phytohormones. Transient expression of the TaARGOS-D protein in wheat protoplasts showed that TaARGOS-D localized to the endoplasmic reticulum. Moreover, overexpression of TaARGOS-D in Arabidopsis resulted in an enhanced germination rate, larger rosette diameter, increased rosette leaf area, and higher silique number than in wild-type (WT) plants. The roles of TaARGOS-D in the control of plant growth were further studied via RNA-seq, and it was found that 105 genes were differentially expressed; most of these genes were involved in 'developmental processes.' Interestingly, we also found that overexpression of TaARGOS-D in Arabidopsis improved drought and salinity tolerance and insensitivity to ABA relative to that in WT plants. Taken together, these results demonstrate that the TaARGOSs are involved in seed germination, seedling growth, and abiotic stress tolerance.

Differences in meristem size and expression of branching genes are associated with variation in panicle phenotype in wild and domesticated African rice

Background

The African rice Oryza glaberrima was domesticated from its wild relative Oryza barthii about 3000 years ago. During the domestication process, panicle complexity changed from a panicle with low complexity in O. barthii, to a highly branched panicle carrying more seeds in O. glaberrima. To understand the basis of this differential panicle development between the two species, we conducted morphological and molecular analyses of early panicle development.

Results

Using X-ray tomography, we analyzed the morphological basis of early developmental stages of panicle development. We uncovered evidence for a wider rachis meristem in O. glaberrima than in O. barthii. At the molecular level, spatial and temporal expression profiles of orthologs of O. sativa genes related to meristem activity and meristem fate control were obtained using in situ hybridization and qRT-PCR. Despite highly conserved spatial expression patterns between O. glaberrima and O. barthii, differences in the expression levels of these early acting genes were detected.

Conclusion

The higher complexity of the O. glaberrima panicle compared to that of its wild relative O. barthii is associated with a wider rachis meristem and a modification of expression of branching-related genes. Our study indicates that the expression of genes in the miR156/miR529/SPL and TAW1 pathways, along with that of their target genes, is altered from the unbranched stage of development. This suggests that differences in panicle complexity between the two African rice species result from early alterations to gene expression during reproductive development.

Electronic supplementary material

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

Polyploidy: Pitfalls and paths to a paradigm

Investigators have long searched for a polyploidy paradigm-rules or principles that might be common following polyploidization (whole-genome duplication, WGD). Here we attempt to integrate what is known across the more thoroughly investigated polyploid systems on topics ranging from genetics to ecology. We found that while certain rules may govern gene retention and loss, systems vary in the prevalence of gene silencing vs. homeolog loss, chromosomal change, the presence of a dominant genome (in allopolyploids), and the relative importance of hybridization vs. genome doubling per se. In some lineages, aspects of polyploidization are repeated across multiple origins, but in other species multiple origins behave more stochastically in terms of genetic and phenotypic change. Our investigation also reveals that the path to synthesis is hindered by numerous gaps in our knowledge of even the best-known systems. Particularly concerning is the absence of linkage between genotype and phenotype. Moreover, most recent studies have focused on the genetic and genomic attributes of polyploidy, but rarely is there an ecological or physiological context. To promote a path to a polyploidy paradigm (or paradigms), we propose a major community goal over the next 10-20 yr to fill the gaps in our knowledge of well-studied polyploids. Before a meaningful synthesis is possible, more complete data sets are needed for comparison-systems that include comparable genetic, genomic, chromosomal, proteomic, as well as morphological, physiological, and ecological data. Also needed are more natural evolutionary model systems, as most of what we know about polyploidy continues to come from a few crop and genetic models, systems that often lack the ecological context inherent in natural systems and necessary for understanding the drivers of biodiversity.

Creating Order from Chaos: Epigenome Dynamics in Plants with Complex Genomes

Flowering plants have strikingly distinct genomes, although they contain a similar suite of expressed genes. The diversity of genome structures and organization is largely due to variation in transposable elements (TEs) and whole-genome duplication (WGD) events. We review evidence that chromatin modifications and epigenetic regulation are intimately associated with TEs and likely play a role in mediating the effects of WGDs. We hypothesize that the current structure of a genome is the result of various TE bursts and WGDs and it is likely that the silencing mechanisms and the chromatin structure of a genome have been shaped by these events. This suggests that the specific mechanisms targeting chromatin modifications and epigenomic patterns may vary among different species. Many crop species have likely evolved chromatin-based mechanisms to tolerate silenced TEs near actively expressed genes. These interactions of heterochromatin and euchromatin are likely to have important roles in modulating gene expression and variability within species.

Altered expression of TaRSL4 gene by genome interplay shapes root hair length in allopolyploid wheat

Polyploidy is a major driving force in plant evolution and speciation. Phenotypic changes often arise with the formation, natural selection and domestication of polyploid plants. However, little is known about the consequence of hybridization and polyploidization on root hair development. Here, we report that root hair length of synthetic and natural allopolyploid wheats is significantly longer than those of their diploid progenitors, whereas no difference is observed between allohexaploid and allotetraploid wheats. The expression of wheat gene TaRSL4, an orthologue of AtRSL4 controlling the root hair development in Arabidopsis, was positively correlated with the root hair length in diploid and allotetraploid wheats. Moreover, transcript abundance of TaRSL4 homoeologue from A genome (TaRSL4-A) was much higher than those of other genomes in natural allopolyploid wheat. Notably, increased root hair length by overexpression of the TaRSL4-A in wheat led to enhanced shoot fresh biomass under nutrient-poor conditions. Our observations indicate that increased root hair length in allohexaploid wheat originated in the allotetraploid progenitors and altered expression of TaRSL4 gene by genome interplay shapes root hair length in allopolyploid wheat.

Homoeologous copy-specific expression patterns of MADS-box genes for floral formation in allopolyploid wheat

The consensus model for floral organ formation in higher plants, the so-called ABCDE model, proposes that floral whorl-specific combinations of class A, B, C, D, and E genes specify floral organ identity. Class A, B, C, D and E genes encode MADS-box transcription factors; the single exception being the class A gene APETALA2. Bread wheat (Triticum aestivum) is a hexaploid species with a genome constitution AABBDD; the hexaploid originated from a cross between tetraploid T. turgidum (AABB) and diploid Aegilops tauschii (DD). Tetraploid wheat is thought to have originated from a cross between the diploid species T. urartu (AA) and Ae. speltoides (BB). Consequently, the hexaploid wheat genome contains triplicated homoeologous copies (homoeologs) of each gene derived from the different ancestral diploid species. In this study, we examined the expression patterns of homoeologs of class B, C and D MADS-box genes during floral development. For the class B gene wheat PISTILLATA2 (WPI2), the homoeologs from the A and D genomes were expressed, while expression of the B genome homoeolog was suppressed. For the class C gene wheat AGAMOUS1 (WAG1), the homoeologs on the A and B genomes were expressed, while expression of the D genome homoeolog was suppressed. For the class D gene wheat SEEDSTICK (WSTK), the B genome homoeolog was preferentially expressed. These differential patterns of homoeolog expression were consistently observed among different hexaploid wheat varieties and synthetic hexaploid wheat lines developed by artificial crosses between tetraploid wheat and Ae. tauschii. These results suggest that homoeolog-specific regulation of the floral MADS-box genes occurs in allopolyploid wheat.

Temporal transcriptome profiling reveals expression partitioning of homeologous genes contributing to heat and drought acclimation in wheat (Triticum aestivum L.)

BACKGROUND

Hexaploid wheat (Triticum aestivum) is a globally important crop. Heat, drought and their combination dramatically reduce wheat yield and quality, but the molecular mechanisms underlying wheat tolerance to extreme environments, especially stress combination, are largely unknown. As an allohexaploid, wheat consists of three closely related subgenomes (A, B, and D), and was reported to show improved tolerance to stress conditions compared to tetraploid. But so far very little is known about how wheat coordinates the expression of homeologous genes to cope with various environmental constraints on the whole-genome level.

RESULTS

To explore the transcriptional response of wheat to the individual and combined stress, we performed high-throughput transcriptome sequencing of seedlings under normal condition and subjected to drought stress (DS), heat stress (HS) and their combination (HD) for 1 h and 6 h, and presented global gene expression reprograms in response to these three stresses. Gene Ontology (GO) enrichment analysis of DS, HS and HD responsive genes revealed an overlap and complexity of functional pathways between each other. Moreover, 4,375 wheat transcription factors were identified on a whole-genome scale based on the released scaffold information by IWGSC, and 1,328 were responsive to stress treatments. Then, the regulatory network analysis of HSFs and DREBs implicated they were both involved in the regulation of DS, HS and HD response and indicated a cross-talk between heat and drought stress. Finally, approximately 68.4 % of homeologous genes were found to exhibit expression partitioning in response to DS, HS or HD, which was further confirmed by using quantitative RT-PCR and Nullisomic-Tetrasomic lines.

CONCLUSIONS

A large proportion of wheat homeologs exhibited expression partitioning under normal and abiotic stresses, which possibly contributes to the wide adaptability and distribution of hexaploid wheat in response to various environmental constraints.