Transcriptome profiling of wheat glumes in wild emmer, hulled landraces and modern cultivars

Kernel Transcriptome Profiles of Susceptible Wheat Genotypes in Response to Wheat Dwarf Bunt

Wheat dwarf bunt is caused by Tilletia controversa J. G. Kühn (TCK), which is a serious fungal diseases affecting kernels of wheat. In order to identify candidate genes involved in the abnormal development of kernels in wheat, we used RNA sequencing technology to analyze the transcriptome of the abnormal and healthy kernels of a susceptible variety (Yili053) at the mid-filling stage, late-filling stage, and maturity stage, respectively. The differentially expressed genes (DEGs) were analyzed, and there were 3930 DEGs, 28,422 DEGs, and 20,874 DEGs found at the mid-filling stage, late-filling stage, and maturity stage in Yili053, respectively. A total of 1592 DEGs (506 DEGs up-regulated) showed continuously differential expression in the three stages. Gene ontology analysis showed that these DEGs were related to biological regulation, metabolic processes, and the response to stimulus. Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that these DEGs play major roles in pathways including photosynthesis, carbon metabolism, carbon fixation in photosynthetic organisms, and glyoxylate and dicarboxylate metabolism. Moreover, we predicted that 13 MADS-MIKC transcription factors, which were continuously up-regulated, were crucial for regulating the maturation and senescence of eukaryotes. Some 21 genes related to the plant hormone signaling transduction pathway and 61 genes related to the response to stimulus were analyzed. A total of 26 of them were successful validated with a qPCR analysis. These genes were thought to be involved in the abnormal development of kernels infected by TCK. A transcriptomics analysis of wheat kernels in response to TCK will contribute to understanding the interaction of TCK and wheat, and may provide a basis for knowledge of molecular events in the abnormal development of kernels, which will be helpful for more efficient TCK management.

The Triticum ispahanicum elongated glume locus P2 maps to chromosome 6A and is associated with the ectopic expression of SVP-A1

We propose the MADS-box transcription factor SVP-A1 as a promising candidate gene for the elongated glume locus P2, which maps to chromosome 6A instead of the previously proposed chromosome 7B. In rice and wheat, glume and floral organ length are positively correlated with grain size, making them an important target to increase grain size and potentially yield. The wheat subspecies Triticum ispahanicum is known to develop elongated glumes and floral organs as well as long grains. These multiple phenotypic effects are controlled by the P2 locus, which was previously mapped to wheat chromosome 7B. Using three mapping populations, we show that the long glume locus P2 does not map to chromosome 7B, but instead maps to a 1.68 Mbp interval on chromosome 6A. Within this interval, we identified SVP-A1, a MADS box transcription factor which is the direct ortholog of the maize gene underlying the 'pod corn' Tunicate locus and is a paralog to the T. polonicum elongated glume P1 gene. In T. ispahanicum, we identified a unique allele which has a 482-bp deletion in the SVP-A1 promoter and is associated with ectopic and higher expression of SVP-A1 in the elongated glumes and floral organs. We used near-isogenic lines (NILs) to show that P2 has a consistent positive effect on the length of glume, lemma, palea, spike and grain. Based on the mapping data, natural variation, biological function of SVP genes in cereals and expression analyses, we propose the MADS-box transcription factor SVP-A1 as a promising candidate for P2.

Wheat Escapes Low Light Stress by Altering Pollination Types

Although low light stress seriously affects florets fertility and grain number during the reproductive period, crops can be fertilized by heterologous pollen to alleviate the reduction of grain number. However, wheat is strongly autogamous, how to change to outcross after low light remains unclear. To understand the mechanisms of this change process, an approach combined morphological, physiological, and transcriptomic analyses was performed under low light stress imposed at the young microspore stage the booting stage from tetrad to uni-nucleate microspores stage. The results showed that low light stress caused pollen abortion, and the unfertilized ovary is fertilized by heterologous pollen after floret opening. Compared to control, the opening angle of lemma and glume were increased by 11.6-48.6 and 48.4-78.5%, respectively. The outcross of stressed wheat compensated for the 2.1-18.0% of grain number loss. During this process, phytohormones played an important role. Jasmonic acid (JA) and methyl jasmonate (MeJA) levels in spikelets were increased. Meanwhile, lignin and cellulose content decreased, and genes associated with cell wall related GO terms were enriched. Among the differentially expressed genes (DEGs), were identified 88-710 transcription factors genes, of which some homologs in Arabidopsis are proposed to function in lignin and cellulose, influencing the glume and lemma opening. Our finding can provide new insight into a survival mechanism to set seeds through pollination way alteration in the absence of self-fertilization after the stress of adversity.

Rising Atmospheric Temperature Impact on Wheat and Thermotolerance Strategies

Temperature across the globe is increasing continuously at the rate of 0.15–0.17 °C per decade since the industrial revolution. It is influencing agricultural crop productivity. Therefore, thermotolerance strategies are needed to have sustainability in crop yield under higher temperature. However, improving thermotolerance in the crop is a challenging task for crop scientists. Therefore, this review work was conducted with the aim of providing information on the wheat response in three research areas, i.e., physiology, breeding, and advances in genetics, which could assist the researchers in improving thermotolerance. The optimum temperature for wheat growth at the heading, anthesis, and grain filling duration is 16 ± 2.3 °C, 23 ± 1.75 °C, and 26 ± 1.53 °C, respectively. The high temperature adversely influences the crop phenology, growth, and development. The pre-anthesis high temperature retards the pollen viability, seed formation, and embryo development. The post-anthesis high temperature declines the starch granules accumulation, stem reserve carbohydrates, and translocation of photosynthates into grains. A high temperature above 40 °C inhibits the photosynthesis by damaging the photosystem-II, electron transport chain, and photosystem-I. Our review work highlighted that genotypes which can maintain a higher accumulation of proline, glycine betaine, expression of heat shock proteins, stay green and antioxidant enzymes activity viz., catalase, peroxidase, super oxide dismutase, and glutathione reductase can tolerate high temperature efficiently through sustaining cellular physiology. Similarly, the pre-anthesis acclimation with heat treatment, inorganic fertilizer such as nitrogen, potassium nitrate and potassium chloride, mulches with rice husk, early sowing, presoaking of a 6.6 mM solution of thiourea, foliar application of 50 ppm dithiothreitol, 10 mg per kg of silicon at heading and zinc ameliorate the crop against the high temperature. Finally, it has been suggested that modern genomics and omics techniques should be used to develop thermotolerance in wheat.

A novel model for the X-chromosome inactivation association on survival data

The X-linked genetic association is overlooked in most of the genetic studies because of the complexity of X-chromosome inactivation process. In fact, the biological process of the gene at the same locus can vary across different subjects. Besides, the skewness of X-chromosome inactivation is inherently subject-specific (even tissue-specific within subjects) and cannot be accurately represented by a population-level parameter. To tackle this issue, a new model is proposed to incorporate the X-linked genetic association into right-censored survival data. The novel model can present that the X-linked genes on different subjects may go through different biological processes via a mixed distribution. Further, a random effect is employed to describe the uncertainty of the biological process for every subject. The proposed method can derive the probability for the escape of X-chromosome inactivation and derive the unbiased estimates of the model parameters. The Legendre-Gauss Quadrature method is used to approximate the integration over the random effect. Finite sample performance of this method is examined via extensive simulation studies. An application is illustrated with the implementation on a cancer genetic study with right-censored survival data.

Comparative transcriptome profiling of multi-ovary wheat under heterogeneous cytoplasm suppression

DUOII is a multi-ovary wheat line with two or three pistils and three stamens in each floret. The multi-ovary trait of DUOII is controlled by a dominant gene, whose expression can be suppressed by the heterogeneous cytoplasm of TeZhiI (TZI), a line with the nucleus of common wheat and the cytoplasm of Aegilops. DUOII (♀) × TZI (♂) shows multi-ovary trait, while TZI (♀) × DUOII (♂) shows mono-ovary. Observing the developmental process, we found that the critical stage of additional pistil primordium development was when the young spikes were 2–6 mm long. To elucidate the molecular mechanisms that are responsible for the heterogeneous cytoplasmic suppression of the multi-ovary gene, we RNA-sequenced the entire transcriptome of 2–6 mm long young spikes obtained from the reciprocal crosses between DUOII and TZI. A total of 600 differentially expressed genes (DEGs) was identified. Functional annotation of these DEGs showed that the heterogeneous cytoplasmic suppression of additional pistil development mainly involved four pathways, i.e., chloroplast metabolism, DNA replication and repair, hormone signal transduction, and trehalose-6-phosphate in the primordium development stage, which cooperated to modulate the multi-ovary gene expression under heterogeneous cytoplasmic suppression.

Btr1-A Induces Grain Shattering and Affects Spike Morphology and Yield-Related Traits in Wheat

Spike brittleness represents an important domestication trait in crops. Although the brittle rachis of wild wheat was cloned, however, the molecular mechanism underlying spike brittleness is yet to be elucidated. Here, we identified a single dominant brittle rachis gene Br-Ab on chromosome arm 3AbS using an F2 population of diploid wheat and designated Btr1-Ab. Sequence analysis of the Btr1-A gene in 40 diploid wheat accessions, 80 tetraploid wheat accessions and 38 hexaploid wheat accessions showed that two independent mutations (Ala119Thr for diploid and Gly97* for polyploids) in the Btr1-A coding region resulting in the nonbrittle rachis allele. Overexpression of Btr1-Ab in nonbrittle hexaploid wheat led to brittle rachis in transgenic plants. RNA-Seq analysis revealed that Btr1-A represses the expression of cell wall biosynthesis genes during wheat rachis development. In addition, we found that Btr1-A can modify spike morphology and reduce threshability, grain size and thousand grain weight in transgenic wheat. These results demonstrated that Btr1-A reduces cell wall synthesis in rachis nodes, resulting in natural spikelet shattering, and that the transition from Btr1-A to btr1-A during wheat domestication had profound effects on evolution of spike morphology and yield-related traits.

Transcriptomic response of durum wheat to cold stress at reproductive stage

Understanding the genetic basis of cold tolerance is a key step towards obtaining new and improved crop varieties. Current geographical distribution of durum wheat in Argentina exposes the plants to frost damage when spikes have already emerged. Biochemical pathways involved in cold tolerance are known to be early activated at above freezing temperatures. In this study we reported the transcriptome of CBW0101 spring durum wheat by merging data from untreated control and cold (5 °C) treated plant samples at reproductive stage. A total of 128,804 unigenes were predicted. Near 62% of the unigenes were annotated in at least one database. In total 876 unigenes were differentially expressed (DEGs), 562 were up-regulated and 314 down-regulated in treated samples. DEGs are involved in many critical processes including, photosynthetic activity, lipid and carbohydrate synthesis and accumulation of amino acids and seed proteins. Twenty-eight transcription factors (TFs) belonging to 14 families resulted differentially expressed from which eight families comprised of only TFs induced by cold. We also found 31 differentially expressed Long non-coding RNAs (lncRNAs), most of them up-regulated in treated plants. Two of these lncRNAs could operate via microRNAs (miRNAs) target mimic. Our results suggest a reprogramming of expression patterns in CBW0101 that affects a number of genes that is closer to the number reported in winter genotypes. These observations could partially explain its moderate tolerance (low proportion of frost-damaged spikes) when exposed to freezing days in the field.

Domestication reduces alternative splicing expression variations in sorghum

Domestication is known to strongly reduce genomic diversity through population bottlenecks. The resulting loss of polymorphism has been thoroughly documented in numerous cultivated species. Here we investigate the impact of domestication on the diversity of alternative transcript expressions using RNAseq data obtained on cultivated and wild sorghum accessions (ten accessions for each pool). In that aim, we focus on genes expressing two isoforms in sorghum and estimate the ratio between expression levels of those isoforms in each accession. Noticeably, for a given gene, one isoform can either be overexpressed or underexpressed in some wild accessions, whereas in the cultivated accessions, the balance between the two isoforms of the same gene appears to be much more homogenous. Indeed, we observe in sorghum significantly more variation in isoform expression balance among wild accessions than among domesticated accessions. The possibility exists that the loss of nucleotide diversity due to domestication could affect regulatory elements, controlling transcription or degradation of these isoforms. Impact on the isoform expression balance is discussed. As far as we know, this is the first time that the impact of domestication on transcript isoform balance has been studied at the genomic scale. This could pave the way towards the identification of key domestication genes with finely tuned isoform expressions in domesticated accessions while being highly variable in their wild relatives.

Novel insights into the effect of nitrogen on storage protein biosynthesis and protein body development in wheat caryopsis

Molecular and cytological mechanisms concerning the effects of nitrogen on wheat (Triticum aestivum L.) storage protein biosynthesis and protein body development remain largely elusive. We used transcriptome sequencing, proteomics techniques, and light microscopy to investigate these issues. In total, 2585 differentially expressed genes (DEGs) and 57 differentially expressed proteins (DEPs) were found 7 days after anthesis (DAA), and 2456 DEGs and 64 DEPs were detected 18 DAA after nitrogen treatment. Gene ontology terms related to protein biosynthesis processes enriched these numbers by 678 and 582 DEGs at 7 and 18 DAA, respectively. Further, 25 Kyoto Encyclopedia of Genes and Genomes pathways were involved in protein biosynthesis at both 7 and 18 DAA. DEPs related to storage protein biosynthesis contained gliadin and glutenin subunits, most of which were up-regulated after nitrogen treatment. Quantitative real-time PCR analysis indicated that some gliadin and glutenin subunit encoding genes were differentially expressed at 18 DAA. Structural observation revealed that wheat endosperm accumulated more and larger protein bodies after nitrogen treatment. Collectively, our findings suggest that nitrogen treatment enhances storage protein content, endosperm protein body quantity, and partial processing quality by altering the expression levels of certain genes involved in protein biosynthesis pathways and storage protein expression at the proteomics level.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.)

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.

Evolution and Adaptation of Wild Emmer Wheat Populations to Biotic and Abiotic Stresses

The genetic bottlenecks associated with plant domestication and subsequent selection in man-made agroecosystems have limited the genetic diversity of modern crops and increased their vulnerability to environmental stresses. Wild emmer wheat, the tetraploid progenitor of domesticated wheat, distributed along a wide range of ecogeographical conditions in the Fertile Crescent, has valuable "left behind" adaptive diversity to multiple diseases and environmental stresses. The biotic and abiotic stress responses are conferred by series of genes and quantitative trait loci (QTLs) that control complex resistance pathways. The study of genetic diversity, genomic organization, expression profiles, protein structure and function of biotic and abiotic stress-resistance genes, and QTLs could shed light on the evolutionary history and adaptation mechanisms of wild emmer populations for their natural habitats. The continuous evolution and adaptation of wild emmer to the changing environment provide novel solutions that can contribute to safeguarding food for the rapidly growing human population.

Molecular markers for tolerance of European ash (Fraxinus excelsior) to dieback disease identified using Associative Transcriptomics

Tree disease epidemics are a global problem, impacting food security, biodiversity and national economies. The potential for conservation and breeding in trees is hampered by complex genomes and long lifecycles, with most species lacking genomic resources. The European Ash tree Fraxinus excelsior is being devastated by the fungal pathogen Hymenoscyphus fraxineus, which causes ash dieback disease. Taking this system as an example and utilizing Associative Transcriptomics for the first time in a plant pathology study, we discovered gene sequence and gene expression variants across a genetic diversity panel scored for disease symptoms and identified markers strongly associated with canopy damage in infected trees. Using these markers we predicted phenotypes in a test panel of unrelated trees, successfully identifying individuals with a low level of susceptibility to the disease. Co-expression analysis suggested that pre-priming of defence responses may underlie reduced susceptibility to ash dieback.