Transcriptome analysis reveals important candidate genes involved in grain-size formation at the stage of grain enlargement in common wheat cultivar "Bainong 4199"

Temporal Gene Expression Profiles From Pollination to Seed Maturity in Sorghum Provide Core Candidates for Engineering Seed Traits

Sorghum (Sorghum bicolor (L.) Moench) is a highly nutritional multipurpose millet crop. However, the genetic and molecular regulatory mechanisms governing sorghum grain development and the associated agronomic traits remain unexplored. In this study, we performed a comprehensive transcriptomic analysis of pistils collected 1-2 days before pollination, and developing seeds collected -2, 10, 20 and 30 days after pollination of S. bicolor variety M35-1. Out of 31 337 genes expressed in these stages, 12 804 were differentially expressed in the consecutive stages of seed development. These exhibited 10 dominant expression patterns correlated with the distinct pathways and gene functions. Functional analysis, based on the pathway mapping, transcription factor enrichment and orthology, delineated the key patterns associated with pollination, fertilization, early seed development, grain filling and seed maturation. Furthermore, colocalization with previously reported quantitative trait loci (QTLs) for grain weight/size revealed 48 differentially expressed genes mapping to these QTL regions. Comprehensive literature mining integrated with QTL mapping and expression data shortlisted 25, 17 and 8 core candidates for engineering grain size, starch and protein content, respectively.

Genomic analysis of Zhou8425B, a key founder parent, reveals its genetic contributions to elite agronomic traits in wheat breeding

High-quality genome information is essential for efficiently deciphering and improving crop traits. Here, we report a highly contiguous and accurate hexaploid genome assembly for the key wheat breeding parent Zhou8425B, an elite 1BL/1RS translocation line with durable adult plant resistance (APR) against yellow rust (YR) disease. By integrating HiFi and Hi-C sequencing reads, we have generated a 14.75-Gb genome assembly for Zhou8425B with a contig N50 of 70.94 and a scaffold N50 of 735.11 Mb. Comparisons with previously sequenced common wheat cultivars shed light on structural changes in the 1RS chromosome arm, which has been extensively used in wheat improvement. Interestingly, Zhou8425B 1RS carries more genes encoding AP2/ERF-ERF or B3 transcription factors than its counterparts in four previously sequenced wheat and rye genotypes. The Zhou8425B genome assembly aided in the fine mapping of a new APR locus (YrZH3BS) that confers resistance to YR disease and promotes grain yield under field conditions. Notably, pyramiding YrZH3BS with two previously characterized APR loci (YrZH22 and YrZH84) can further reduce YR severity and enhance grain yield, with the triple combination (YrZH3B + YrZH22 + YrZH84) having the greatest effect. Finally, the founder genotype effects of Zhou8425B were explored using publicly available genome resequencing data, which reveals the presence of important Zhou8425B genomic blocks in its derivative cultivars. Our data demonstrate the value of the Zhou8425B genome assembly for further study of the structural and functional characteristics of 1RS, the genetic basis of durable YR resistance, and founder genotype effects in wheat breeding. Our resources will facilitate the development of elite wheat cultivars through genomics-assisted breeding.

Starch and storage protein dynamics in the developing and matured grains of durum wheat and diploid progenitor species

Durum wheat, less immunogenically intolerant than bread wheat, originates from diploid progenitors known for nutritional quality and stress tolerance. Present study involves the analysis of major grain parameters, viz. size, weight, sugar, starch, and protein content of Triticum durum (AABB genome) and its diploid progenitors, Triticum monococcum (AA genome) and Aegilops speltoides (BB genome). Samples were collected during 2-5 weeks after anthesis (WAA), and at maturity. The investigation revealed that T. durum displayed the maximum grain size and weight. Expression analysis of Grain Weight 2 (GW2) and Glutamine Synthase (GS2), negative and positive regulators of grain weight and size, respectively, revealed higher GW2 expression in Ae. speltoides and higher GS2 expression in T. durum. Further we explored total starch, sugar and protein content, observing higher levels of starch and sugar in durum wheat while AA genome species exhibited higher protein content dominated by the fractions of albumin/globulin. HPLC profiling revealed unique sub-fractions in all three genome species. Additionally, a comparative transcriptome analysis also corroborated with the starch and protein content in the grains. This study provides valuable insights into the genetic and biochemical distinctions among durum wheat and its diploid progenitors, offering a foundation for their nutritional composition.

Gene Expression Analysis of Different Organs and Identification of AP2 Transcription Factors in Flax (Linum usitatissimum L.)

Flax (Linum usitatissimum L.) is an important oilseed crop widely cultivated for its oil and fiber. This study conducted transcriptome analysis to analyze the gene expression profiles of roots, leaves, stamens, pistils, and fruits in the flax cultivar Longya10. A total of 43,471 genes were detected in the RNA-seq data, with 34,497 genes showing differential expression levels between different organs. Gene expression patterns varied across different organs, with differences observed in expression-regulating genes within specific organs. However, 23,448 genes were found to be commonly expressed across all organs. Further analysis revealed organ-specific gene expressions, with 236, 690, 544, 909, and 1212 genes identified in pistils, fruits, leaves, roots, and stamens, respectively. Gene Ontology (GO) enrichment analysis was performed on these organ-specific genes, and significant enrichment was observed in various biological processes, cellular components, and molecular functions, providing new insights for the specific growth patterns of flax organs. Furthermore, we investigated the expression differences of AP2 transcription factors in various tissues and organs of Longya10. We identified 96 AP2 genes that were differentially expressed in different organs and annotated them into various biological pathways. Our results suggest that AP2 transcription factors may play important roles in regulating the growth and development of flax organs including stress response. In summary, our study provides a comprehensive analysis of gene expression patterns in different organs and tissues of flax plant and identifies potential critical regulators of flax organ growth and development. These findings contribute to a better understanding of the molecular mechanisms underlying flax organ development and may have important implications for the genetic improvement of flax crops.

Evaluation of altered starch mutants and identification of candidate genes responsible for starch variation in wheat

BACKGROUND

Induction of mutation through chemical mutagenesis is a novel approach for preparing diverse germplasm. Introduction of functional alleles in the starch biosynthetic genes help in the improvement of the quality and yield of cereals.

RESULTS

In the present study, a set of 350 stable mutant lines were used to evaluate dynamic variation of the total starch contents. A megazyme kits were used for measuring the total starch content, resistant starch, amylose, and amylopectin content. Analysis of variance showed significant variation (p < 0.05) in starch content within the population. Furthermore, two high starch mutants (JE0173 and JE0218) and two low starch mutants (JE0089 and JE0418) were selected for studying different traits. A multiple comparison test showed that significant variation in all physiological and morphological traits, with respect to the parent variety (J411) in 2019-2020 and 2020-2021. The quantitative expression of starch metabolic genes revealed that eleven genes of JE0173 and twelve genes of JE0218 had consistent expression in high starch mutant lines. Similarly, in low starch mutant lines, eleven genes of JE0089 and thirteen genes of JE0418 had consistent expression in all stages of seed development. An additional two candidate genes showed over-expression (PHO1, PUL) in the high starch mutant lines, indicating that other starch metabolic genes may also contribute to the starch biosynthesis. The overexpression of SSII, SSIII and SBEI in JE0173 may be due to presence of missense mutations in these genes and SSI also showed overexpression which may be due to 3-primer_UTR variant. These mutations can affect the other starch related genes and help to increase the starch content in this mutant line (JE0173).

CONCLUSIONS

This study screened a large scale of mutant population and identified mutants, could provide useful genetic resources for the study of starch biosynthesis and genetic improvement of wheat in the future. Further study will help to understand new genes which are responsible for the fluctuation of total starch.

Transcriptome Analysis Reveals Potential Mechanism in Storage Protein Trafficking within Developing Grains of Common Wheat

Gluten proteins are the major storage protein fraction in the mature wheat grain. They are restricted to the starchy endosperm, which defines the viscoelastic properties of wheat dough. The synthesis of these storage proteins is controlled by the endoplasmic reticulum (ER) and is directed into the vacuole via the Golgi apparatus. In the present study, transcriptome analysis was used to explore the potential mechanism within critical stages of grain development of wheat cultivar "Shaannong 33" and its sister line used as the control (CK). Samples were collected at 10 DPA (days after anthesis), 14 DPA, 20 DPA, and 30 DPA for transcriptomic analysis. The comparative transcriptome analysis identified that a total of 18,875 genes were differentially expressed genes (DEGs) between grains of four groups "T10 vs. CK10, T14 vs. CK14, T20 vs. CK20, and T30 vs. CK30", including 2824 up-regulated and 5423 down-regulated genes in T30 vs. CK30. Further, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment highlighted the maximum number of genes regulating protein processing in the endoplasmic reticulum (ER) during grain enlargement stages (10-20 DPA). In addition, KEGG database analysis reported 1362 and 788 DEGs involved in translation, ribosomal structure, biogenesis, flavonoid biosynthesis pathway and intracellular trafficking, secretion, and vesicular transport through protein processing within ER pathway (ko04141). Notably, consistent with the higher expression of intercellular storage protein trafficking genes at the initial 10 DPA, there was relatively low expression at later stages. Expression levels of nine randomly selected genes were verified by qRT-PCR, which were consistent with the transcriptome data. These data suggested that the initial stages of "cell division" played a significant role in protein quality control within the ER, thus maintaining the protein quality characteristics at grain maturity. Furthermore, our data suggested that the protein synthesis, folding, and trafficking pathways directed by a different number of genes during the grain enlargement stage contributed to the observed high-quality characteristics of gluten protein in Shaannong 33 (Triticum aestivum L.).

Gene Expression in the Developing Seed of Wild and Domesticated Rice

The composition and nutritional properties of rice are the product of the expression of genes in the developing seed. RNA-Seq was used to investigate the level of gene expression at different stages of seed development in domesticated rice (Oryza sativa ssp. japonica var. Nipponbare) and two Australian wild taxa from the primary gene pool of rice (Oryza meridionalis and Oryza rufipogon type taxa). Transcriptome profiling of all coding sequences in the genome revealed that genes were significantly differentially expressed at different stages of seed development in both wild and domesticated rice. Differentially expressed genes were associated with metabolism, transcriptional regulation, nucleic acid processing, and signal transduction with the highest number of being linked to protein synthesis and starch/sucrose metabolism. The level of gene expression associated with domestication traits, starch and sucrose metabolism, and seed storage proteins were highest at the early stage (5 days post anthesis (DPA)) to the middle stage (15 DPA) and declined late in seed development in both wild and domesticated rice. However, in contrast, black hull colour (Bh4) gene was significantly expressed throughout seed development. A substantial number of novel transcripts (38) corresponding to domestication genes, starch and sucrose metabolism, and seed storage proteins were identified. The patterns of gene expression revealed in this study define the timing of metabolic processes associated with seed development and may be used to explain differences in rice grain quality and nutritional value.

Major Genomic Regions for Wheat Grain Weight as Revealed by QTL Linkage Mapping and Meta-Analysis

Grain weight is a key determinant for grain yield potential in wheat, which is highly governed by a type of quantitative genetic basis. The identification of major quantitative trait locus (QTL) and functional genes are urgently required for molecular improvements in wheat grain yield. In this study, major genomic regions and putative candidate genes for thousand grain weight (TGW) were revealed by integrative approaches with QTL linkage mapping, meta-analysis and transcriptome evaluation. Forty-five TGW QTLs were detected using a set of recombinant inbred lines, explaining 1.76-12.87% of the phenotypic variation. Of these, ten stable QTLs were identified across more than four environments. Meta-QTL (MQTL) analysis were performed on 394 initial TGW QTLs available from previous studies and the present study, where 274 loci were finally refined into 67 MQTLs. The average confidence interval of these MQTLs was 3.73-fold less than that of initial QTLs. A total of 134 putative candidate genes were mined within MQTL regions by combined analysis of transcriptomic and omics data. Some key putative candidate genes similar to those reported early for grain development and grain weight formation were further discussed. This finding will provide a better understanding of the genetic determinants of TGW and will be useful for marker-assisted selection of high yield in wheat breeding.

Transcriptome analysis reveals important candidate gene families related to oligosaccharides biosynthesis in Morindaofficinalis

Morinda officinalis How (MO) is one of the best-known traditional herbs and is widely cultivated in subtropical and tropical areas for many years, especially in southern China. Oligosaccharides are the major constituents in the roots of MO, which is well known for its therapeutic effects with anti-depression, anti-osteoporosis, memory-enhancing, ect. To date, the main gene families that regulate the biosynthetic pathway of MO oligosaccharides metabolism yet have been published. In our study, six cDNA libraries generated from six plants of MO were sequenced utilizing an Illumina HiSeq 4000 platform. Corresponding totals of more than 132.60 million clean reads were obtained from the six libraries and assembled into 25,812 unigenes with an average length of 1288 bp. Moreover, 6036 unigenes were found to be allocated to 26 pathways maps using several public databases, and 2538 differential expression genes (DEGs) were screened. Among them, 25 genes from three families were selected as the mainly candidate genes related to MO oligosaccharides biosynthesis. Then, the expression patterns of six DEGs closely related to MO oligosaccharides biosynthesis were verified by quantitative real-time PCR (qRT-PCR). Besides, the MO was clustered more closely to Coffea arabica of Rubiaceae. In summary, the transcriptomic analysis was used to investigate the differences in expression genes of oligosaccharides biosynthesis, with the notable outcome that several key gene families were closely linked to oligosaccharides biosynthesis.

OsBRKq1, Related Grain Size Mapping, and Identification of Grain Shape Based on QTL Mapping in Rice

The world population is growing rapidly, and food shortage remains a critical issue. Quantitative trait locus (QTL) mapping is a statistical analytical method that uses both phenotypic and genotypic data. The purpose of QTL mapping is to determine the exact gene location for various complex traits. Increasing grain weight is a way to increase yield in rice. Genes related to grain size were mapped using the Samgang/Nagdong double haploid (SNDH) populations. Grain sizes were diversely distributed in SNDH 113 populations, and OsBRKq1 was detected on chromosome 1 in an analysis of QTL mapping that used 1000 grain weight, grain length, and grain width. OsBRKq1 exhibited high sequence similarity with the brassinosteroid leucine-rich repeat-receptor kinases of Arabidopsis thaliana and Zea mays. It was also predicted to have a similar function because of its high homology. OsBRKq1 interacts with various grain-size control genes. Among the SNDH populations, the analysis of the relative expression level during the panicle formation stage of OsBRKq1 in panicles of SNDH117, which has the largest grain size, and SNDH6, which has the smallest grain size, the relative expression level was significantly increased in SNDH117 panicles. SNDH populations have been advancing generations for 10 years; various genetic traits have been fixed and are currently being used as bridging parents. Therefore, the stable expression level of OsBRKq1 was confirmed via QTL mapping. In the future, OsBRKq1 can be effectively used to increase the yield of rice and solve food problems by increasing the size of seeds.

Genome-wide transcriptome profiling indicates the putative mechanism underlying enhanced grain size in a wheat mutant

Grain size is an important trait for crops. The endogenous hormones brassinosteroids (BRs) play key roles in grain size and mass. In this study, we identified an ethyl methylsulfonate (EMS) mutant wheat line, SM482gs, with increased grain size, 1000-grain weight, and protein content, but decreased starch content, compared with the levels in the wild type (WT). Comparative transcriptomic analysis of SM482gs and WT at four developmental stages [9, 15, 20, and 25 days post-anthesis (DPA)] revealed a total of 264, 267, 771, and 1038 differentially expressed genes (DEGs) at these stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) database analysis showed that some DEGs from the comparison at 15 DPA were involved in the pathway of "brassinosteroid biosynthesis," and eight genes involved in BR biosynthesis and signal transduction were significantly upregulated in SM482gs during at least one stage. This indicated that the enhanced BR signaling in SM482gs might have contributed to its increased grain size via network interactions. The expression of seed storage protein (SSP)-encoding genes in SM482gs was upregulated, mostly at 15 and 20 DPA, while most of the starch synthetase genes showed lower expression in SM482gs at all stages, compared with that in WT. The expression patterns of starch synthase genes and seed storage protein-encoding genes paralleled the decreased level of starch and increased storage protein content of SM482gs, which might be related to the increased seed weight and wrinkled phenotype.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-020-02579-6.

Nanopore RNA Sequencing Revealed Long Non-Coding and LTR Retrotransposon-Related RNAs Expressed at Early Stages of Triticale SEED Development

The intergenic space of plant genomes encodes many functionally important yet unexplored RNAs. The genomic loci encoding these RNAs are often considered “junk”, DNA as they are frequently associated with repeat-rich regions of the genome. The latter makes the annotations of these loci and the assembly of the corresponding transcripts using short RNAseq reads particularly challenging. Here, using long-read Nanopore direct RNA sequencing, we aimed to identify these “junk” RNA molecules, including long non-coding RNAs (lncRNAs) and transposon-derived transcripts expressed during early stages (10 days post anthesis) of seed development of triticale (AABBRR, 2n = 6x = 42), an interspecific hybrid between wheat and rye. Altogether, we found 796 lncRNAs and 20 LTR retrotransposon-related transcripts (RTE-RNAs) expressed at this stage, with most of them being previously unannotated and located in the intergenic as well as intronic regions. Sequence analysis of the lncRNAs provide evidence for the frequent exonization of Class I (retrotransposons) and class II (DNA transposons) transposon sequences and suggest direct influence of “junk” DNA on the structure and origin of lncRNAs. We show that the expression patterns of lncRNAs and RTE-related transcripts have high stage specificity. In turn, almost half of the lncRNAs located in Genomes A and B have the highest expression levels at 10–30 days post anthesis in wheat. Detailed analysis of the protein-coding potential of the RTE-RNAs showed that 75% of them carry open reading frames (ORFs) for a diverse set of GAG proteins, the main component of virus-like particles of LTR retrotransposons. We further experimentally demonstrated that some RTE-RNAs originate from autonomous LTR retrotransposons with ongoing transposition activity during early stages of triticale seed development. Overall, our results provide a framework for further exploration of the newly discovered lncRNAs and RTE-RNAs in functional and genome-wide association studies in triticale and wheat. Our study also demonstrates that Nanopore direct RNA sequencing is an indispensable tool for the elucidation of lncRNA and retrotransposon transcripts.

Metabolic and transcriptomic analysis of two Cucurbita moschata germplasms throughout fruit development

BACKGROUND

Pumpkins (Cucurbita moschata; Cucurbitaceae) are valued for their fruits and seeds and are rich in nutrients. Carotenoids and sugar contents, as main feature of pumpkin pulp, are used to determine the fruit quality.

RESULTS

Two pumpkin germplasms, CMO-X and CMO-E, were analyzed regarding the essential quality traits such as dry weight, soluble solids, organic acids, carotenoids and sugar contents. For the comparison of fruit development in these two germplasms, fruit transcriptome was analyzed at 5 different developmental stages from 0 d to 40 d in a time course manner. Putative pathways for carotenoids biosynthesis and sucrose metabolism were developed in C. moschata fruit and homologs were identified for each key gene involved in the pathways. Gene expression data was found consistent with the accumulation of metabolites across developmental stages and also between two germplasms. PSY, PDS, ZEP, CRTISO and SUS, SPS, HK, FK were found highly correlated with the accumulation of carotenoids and sucrose metabolites, respectively, at different growth stages of C. moschata as shown by whole transcriptomic analysis. The results of qRT-PCR analysis further confirmed the association of these genes.

CONCLUSION

Developmental regulation of the genes associated with the metabolite accumulation can be considered as an important factor for the determination of C. moschata fruit quality. This research will facilitate the investigation of metabolic profiles in other cultivars.

Transcriptome Analysis Reveals Important Candidate Genes Related to Nutrient Reservoir, Carbohydrate Metabolism, and Defence Proteins during Grain Development of Hexaploid Bread Wheat and Its Diploid Progenitors

Wheat grain development after anthesis is an important biological process, in which major components of seeds are synthesised, and these components are further required for germination and seed vigour. We have made a comparative RNA-Seq analysis between hexaploid wheat and its individual diploid progenitors to know the major differentially expressed genes (DEGs) involved during grain development. Two libraries from each species were generated with an average of 55.63, 55.23, 68.13, and 103.81 million reads, resulting in 79.3K, 113.7K, 90.6K, and 121.3K numbers of transcripts in AA, BB, DD, and AABBDD genome species respectively. Number of expressed genes in hexaploid wheat was not proportional to its genome size, but marginally higher than that of its diploid progenitors. However, to capture all the transcripts in hexaploid wheat, sufficiently higher number of reads was required. Functional analysis of DEGs, in all the three comparisons, showed their predominance in three major classes of genes during grain development, i.e., nutrient reservoirs, carbohydrate metabolism, and defence proteins; some of them were subsequently validated through real time quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR). Further, developmental stage-specific gene expression showed most of the defence protein genes expressed during initial developmental stages in hexaploid contrary to the diploids at later stages. Genes related to carbohydrates anabolism expressed during early stages, whereas catabolism genes expressed at later stages in all the species. However, no trend was observed in case of different nutrient reservoirs gene expression. This data could be used to study the comparative gene expression among the three diploid species and homeologue-specific expression in hexaploid.