Transcriptome analysis identified aberrant gene expression in pollen developmental pathways leading to CGMS in cotton (Gossypium hirsutum L.)

Genomic insights into CKX genes: key players in cotton fibre development and abiotic stress responses

Cytokinin oxidase/dehydrogenase (CKX), responsible for irreversible cytokinin degradation, also controls plant growth and development and response to abiotic stress. While the CKX gene has been studied in other plants extensively, its function in cotton is still unknown. Therefore, a genome-wide study to identify the CKX gene family in the four cotton species was conducted using transcriptomics, quantitative real-time PCR (qRT-PCR) and bioinformatics. As a result, in G. hirsutum and G. barbadense (the tetraploid cotton species), 87 and 96 CKX genes respectively and 62 genes each in G. arboreum and G. raimondii, were identified. Based on the evolutionary studies, the cotton CKX gene family has been divided into five distinct subfamilies. It was observed that CKX genes in cotton have conserved sequence logos and gene family expansion was due to segmental duplication or whole genome duplication (WGD). Collinearity and multiple synteny studies showed an expansion of gene families during evolution and purifying selection pressure has been exerted. G. hirsutum CKX genes displayed multiple exons/introns, uneven chromosomal distribution, conserved protein motifs, and cis-elements related to growth and stress in their promoter regions. Cis-elements related to resistance, physiological metabolism and hormonal regulation were identified within the promoter regions of the CKX genes. Expression analysis under different stress conditions (cold, heat, drought and salt) revealed different expression patterns in the different tissues. Through virus-induced gene silencing (VIGS), the GhCKX34A gene was found to improve cold resistance by modulating antioxidant-related activity. Since GhCKX29A is highly expressed during fibre development, we hypothesize that the increased expression of GhCKX29A in fibres has significant effects on fibre elongation. Consequently, these results contribute to our understanding of the involvement of GhCKXs in both fibre development and response to abiotic stress.

Dynamic roles of small RNAs and DNA methylation associated with heterosis in allotetraploid cotton (Gossypium hirsutum L.)

BACKGROUND

Heterosis is a complex phenomenon wherein the hybrids outperform their parents. Understanding the underlying molecular mechanism by which hybridization leads to higher yields in allopolyploid cotton is critical for effective breeding programs. Here, we integrated DNA methylation, transcriptomes, and small RNA profiles to comprehend the genetic and molecular basis of heterosis in allopolyploid cotton at three developmental stages.

RESULTS

Transcriptome analysis revealed that numerous DEGs responsive to phytohormones (auxin and salicylic acid) were drastically altered in F1 hybrid compared to the parental lines. DEGs involved in energy metabolism and plant growth were upregulated, whereas DEGs related to basal defense were downregulated. Differences in homoeologous gene expression in F1 hybrid were greatly reduced after hybridization, suggesting that higher levels of parental expression have a vital role in heterosis. Small RNAome and methylome studies showed that the degree of DNA methylation in hybrid is higher when compared to the parents. A substantial number of allele-specific expression genes were found to be strongly regulated by CG allele-specific methylation levels. The hybrid exhibited higher 24-nt-small RNA (siRNA) expression levels than the parents. The regions in the genome with increased levels of 24-nt-siRNA were chiefly related to genes and their flanking regulatory regions, demonstrating a possible effect of these molecules on gene expression. The transposable elements correlated with siRNA clusters in the F1 hybrid had higher methylation levels but lower expression levels, which suggest that these non-additively expressed siRNA clusters, reduced the activity of transposable elements through DNA methylation in the hybrid.

CONCLUSIONS

These multi-omics data provide insights into how changes in epigenetic mechanisms and gene expression patterns can lead to heterosis in allopolyploid cotton. This makes heterosis a viable tool in cotton breeding.

Transcriptome analysis provides insights into the stress response in cultivated peanut (Arachis hypogaea L.) subjected to drought-stress

BACKGROUND

Peanut (Arachis hypogaea L.) is one of the valuable oilseed crops grown in drought-prone areas worldwide. Drought severely limits peanut production and productivity significantly.

METHOD AND RESULTS

In order to decipher the drought tolerance mechanism in peanut under drought stress, RNA sequencing was performed in TAG - 24 (drought tolerant genotype) and JL-24 (drought susceptible genotype). Approximately 51 million raw reads were generated from four different libraries of two genotypes subjected to drought stress exerted by 20% PEG 6000 stress and control conditions, of which ~ 41 million (80.87%) filtered reads were mapped to the Arachis hypogaea L. reference genome. The transcriptome analysis detected 1,629 differentially expressed genes (DEGs), 186 genes encoding transcription factors (TFs) and 30,199 SSR among the identified DEGs. Among the differentially expressed TF encoding genes, the highest number of genes were WRKY followed by bZIP, C2H2, and MYB during drought stress. The comparative analysis between the two genotypes revealed that TAG-24 exhibits activation of certain key genes and transcriptional factors that are involved in essential biological processes. Specifically, TAG-24 showed activation of genes involved in the plant hormone signaling pathway such as PYL9, Auxin response receptor gene, and ABA. Additionally, genes related to water deprivation such as LEA protein and those involved in combating oxidative damage such as Glutathione reductase were also found to be activated in TAG-24.

CONCLUSION

This genome-wide transcription map, therefore, provides a valuable tool for future transcript profiling under drought stress and enriches the genetic resources available for this important oilseed crop.

Comparative Transcriptome Profiling of CMS-D2 and CMS-D8 Systems Characterizes Fertility Restoration Genes Network in Upland Cotton

Resolving the genetic basis of fertility restoration for cytoplasmic male sterility (CMS) can improve the efficiency of three-line hybrid breeding. However, the genetic determinants of male fertility restoration in cotton are still largely unknown. This study comprehensively compared the full-length transcripts of CMS-D2 and CMS-D8 systems to identify potential genes linked with fertility restorer genes Rf₁ or Rf₂. Target comparative analysis revealed a higher percentage of differential genes in each restorer line as compared to their corresponding sterile and maintainer lines. An array of genes with specific expression in the restorer line of CMS-D2 had functional annotations related to floral development and pathway enrichments in various secondary metabolites, while specifically expressed genes in the CMS-D8 restorer line showed functional annotations related to anther development and pathway enrichment in the biosynthesis of secondary metabolites. Further analysis identified potentially key genes located in the target region of fertility restorer genes Rf₁ or Rf₂. In particular, Ghir_D05G032450 can be the candidate gene related to restorer gene Rf₁, and Ghir_D05G035690 can be the candidate gene associated with restorer gene Rf₂. Further gene expression validation with qRT-PCR confirmed the accuracy of our results. Our findings provide useful insights into decoding the potential regulatory network that retrieves pollen fertility in cotton and will help to further reveal the differences in the genetic basis of fertility restoration for two CMS systems.

Unravelling the treasure trove of drought-responsive genes in wild-type peanut through transcriptomics and physiological analyses of root

Peanut is one of the most valuable legumes, grown mainly in arid and semi-arid regions, where its production may be hindered by the lack of water. Therefore, breeding drought tolerant varieties is of great importance for peanut breeding programs around the world. Unlike cultivated peanuts, wild peanuts have greater genetic diversity and are an important source of alleles conferring tolerance/resistance to abiotic and biotic stresses. To decipher the transcriptome changes under drought stress, transcriptomics of roots of highly tolerant Arachis duranensis (ADU) and moderately susceptible A. stenosperma (AST) genotypes were performed. Transcriptome analysis revealed an aggregate of 1465 differentially expressed genes (DEGs), and among the identified DEGs, there were 366 single nucleotide polymorphisms (SNPs). Gene ontology and Mapman analyses revealed that the ADU genotype had a higher number of transcripts related to DNA methylation or demethylation, phytohormone signal transduction and flavonoid production, transcription factors, and responses to ethylene. The transcriptome analysis was endorsed by qRT-PCR, which showed a strong correlation value (R² = 0.96). Physio-biochemical analysis showed that the drought-tolerant plants produced more osmolytes, ROS phagocytes, and sugars, but less MDA, thus attenuating the effects of drought stress. In addition, three SNPs of the gene encoding transcription factor NFAY (Aradu.YE2F8), expansin alpha (Aradu.78HGD), and cytokinin dehydrogenase 1-like (Aradu.U999X) exhibited polymorphism in selected different genotypes. Such SNPs could be useful for the selection of drought-tolerant genotypes.

Stage specific comparative transcriptomic analysis to reveal gene networks regulating iron and zinc content in pearl millet [Pennisetum glaucum (L.) R. Br.]

Pearl millet is an important staple food crop of poor people and excels all other cereals due to its unique features of resilience to adverse climatic conditions. It is rich in micronutrients like iron and zinc and amenable for focused breeding for these micronutrients along with high yield. Hence, this is a key to alleviate malnutrition and ensure nutritional security. This study was conducted to identify and validate candidate genes governing grain iron and zinc content enabling the desired modifications in the genotypes. Transcriptome sequencing using ION S5 Next Generation Sequencer generated 43.5 million sequence reads resulting in 83,721 transcripts with N₅₀ of 597 bp and 84.35% of transcripts matched with the pearl millet genome assembly. The genotypes having high iron and zinc showed differential gene expression during different stages. Of which, 155 were up-regulated and 251 were down-regulated while during flowering stage and milking stage 349 and 378 transcripts were differentially expressed, respectively. Gene annotation and GO term showed the presence of transcripts involved in metabolic activities associated with uptake and transport of iron and zinc. Information generated will help in gaining insights into iron and zinc metabolism and develop genotypes with high yield, grain iron and zinc content.

Comparative transcriptome analysis of fiber and nonfiber tissues to identify the genes preferentially expressed in fiber development in Gossypium hirsutum

Cotton is an important natural fiber crop and economic crop worldwide. The quality of cotton fiber directly determines the quality of cotton textiles. Identifying cotton fiber development-related genes and exploring their biological functions will not only help to better understand the elongation and development mechanisms of cotton fibers but also provide a theoretical basis for the cultivation of new cotton varieties with excellent fiber quality. In this study, RNA sequencing technology was used to construct transcriptome databases for different nonfiber tissues (root, leaf, anther and stigma) and fiber developmental stages (7 days post-anthesis (DPA), 14 DPA, and 26 DPA) of upland cotton Coker 312. The sizes of the seven transcriptome databases constructed ranged from 4.43 to 5.20 Gb, corresponding to approximately twice the genome size of Gossypium hirsutum (2.5 Gb). Among the obtained clean reads, 83.32% to 88.22% could be compared to the upland cotton TM-1 reference genome. By analyzing the differential gene expression profiles of the transcriptome libraries of fiber and nonfiber tissues, we obtained 1205, 1135 and 937 genes with significantly upregulated expression at 7 DPA, 14 DPA and 26 DPA, respectively, and 124, 179 and 213 genes with significantly downregulated expression. Subsequently, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analyses were performed, which revealed that these genes were mainly involved in catalytic activity, carbohydrate metabolism, the cell membrane and organelles, signal transduction and other functions and metabolic pathways. Through gene annotation analysis, many transcription factors and genes related to fiber development were screened. Thirty-six genes were randomly selected from the significantly upregulated genes in fiber, and expression profile analysis was performed using qRT-PCR. The results were highly consistent with the gene expression profile analyzed by RNA-seq, and all of the genes were specifically or predominantly expressed in fiber. Therefore, our RNA sequencing-based comparative transcriptome analysis will lay a foundation for future research to provide new genetic resources for the genetic engineering of improved cotton fiber quality and for cultivating new transgenic cotton germplasms for fiber quality improvement.

Genic microsatellite marker characterization and development in little millet (Panicum sumatrense) using transcriptome sequencing

Little millet is a climate-resilient and high-nutrient value plant. The lack of molecular markers severely limits the adoption of modern genomic approaches in millet breeding studies. Here the transcriptome of three samples were sequenced. A total of 4443 genic-SSR motifs were identified in 30,220 unigene sequences. SSRs were found at a rate of 12.25 percent, with an average of one SSR locus per 10 kb. Among different repeat motifs, tri-nucleotide repeat (66.67) was the most abundant one, followed by di- (27.39P), and tetra- (3.83P) repeats. CDS contained fewer motifs with the majority of tri-nucleotides, while 3' and 5' UTR carry more motifs but have shorter repeats. Functional annotation of unigenes containing microsatellites, revealed that most of them were linked to metabolism, gene expression regulation, and response to environmental stresses. Fifty primers were randomly chosen and validated in five little millet and 20 minor millet genotypes; 48% showed polymorphism, with a high transferability (70%) rate. Identified microsatellites can be a noteworthy resource for future research into QTL-based breeding, genetic resource conservation, MAS selection, and evolutionary genetics.

Heat stress response mechanisms in pollen development

Being rooted in place, plants are faced with the challenge of responding to unfavourable local conditions. One such condition, heat stress, contributes massively to crop losses globally. Heatwaves are predicted to increase, and it is of vital importance to generate crops that are tolerant to not only heat stress but also to several other abiotic stresses (e.g. drought stress, salinity stress) to ensure that global food security is protected. A better understanding of the molecular mechanisms that underlie the temperature stress response in pollen will be a significant step towards developing effective breeding strategies for high and stable production in crop plants. While most studies have focused on the vegetative phase of plant growth to understand heat stress tolerance, it is the reproductive phase that requires more attention as it is more sensitive to elevated temperatures. Every phase of reproductive development is affected by environmental challenges, including pollen and ovule development, pollen tube growth, male-female cross-talk, fertilization, and embryo development. In this review we summarize how pollen is affected by heat stress and the molecular mechanisms employed during the stress period, as revealed by classical and -omics experiments.

MicroRNAs Involved in Regulatory Cytoplasmic Male Sterility by Analysis RNA-seq and Small RNA-seq in Soybean

Cytoplasmic male sterility (CMS) is an important plant characteristic for exploiting heterosis to enhance crop traits during breeding. However, the CMS regulatory network remains unclear in plants, even though researchers have attempted to isolate genes associated with CMS. In this study, we performed high-throughput sequencing and degradome analyses to identify microRNAs (miRNAs) and their targets in a soybean CMS line (JLCMS9A) and its maintainer line (JLCMS9B). Additionally, the differentially expressed genes during reproductive development were identified using RNA-seq data. A total of 280 miRNAs matched soybean miRNA sequences in miRBase, including mature miRNAs and pre-miRNAs. Of the 280 miRNAs, 30, 23, and 21 belonged to the miR166, miR156, and miR171 families, respectively. Moreover, 410 novel low-abundant miRNAs were identified in the JLCMS9A and JLCMS9B flower buds. Furthermore, 303 and 462 target genes unique to JLCMS9A and JLCMS9B, respectively, as well as 782 common targets were predicted based on the degradome analysis. Target genes differentially expressed between the CMS line and the maintainer line were revealed by an RNA-seq analysis. Moreover, all target genes were annotated with diverse functions related to biological processes, cellular components, and molecular functions, including transcriptional regulation, the nucleus, meristem maintenance, meristem initiation, cell differentiation, auxin-activated signaling, plant ovule development, and anther development. Finally, a network was built based on the interactions. Analyses of the miRNA, degradome, and transcriptome datasets generated in this study provided a comprehensive overview of the reproductive development of a CMS soybean line. The data presented herein represent useful information for soybean hybrid breeding. Furthermore, the study results indicate that miRNAs might contribute to the soybean CMS regulatory network by modulating the expression of CMS-related genes. These findings lay the foundation for future studies on the molecular mechanisms underlying soybean CMS.

Genome-wide comparative transcriptome analysis of the A4-CMS line ICPA 2043 and its maintainer ICPB 2043 during the floral bud development of pigeonpea

Cytoplasmic male sterility (CMS) offers a unique system to understand cytoplasmic nuclear crosstalk, and is also employed for exploitation of hybrid vigor in various crops. Pigeonpea A4-CMS, a predominant source of male sterility, is being used for efficient hybrid seed production. The molecular mechanisms of CMS trait remain poorly studied in pigeonpea. We performed genome-wide transcriptome profiling of A4-CMS line ICPA 2043 and its isogenic maintainer ICPB 2043 at two different stages of floral bud development (stage S1 and stage S2). Consistent with the evidences from some other crops, we also observed significant difference in the expression levels of genes in the later stage, i.e., stage S2. Differential expression was observed for 143 and 55 genes within the two stages of ICPA 2043 and ICPB 2043, respectively. We obtained only 10 differentially expressed genes (DEGs) between the stage S1 of the two genotypes, whereas expression change was significant for 582 genes in the case of stage S2. The qRT-PCR assay of randomly selected six genes supported the differential expression of genes between ICPA 2043 and ICPB 2043. Further, GO and KEGG pathway mapping suggested a possible compromise in key bioprocesses during flower and pollen development. Besides providing novel insights into the functional genomics of CMS trait, our results were in strong agreement with the gene expression atlas of pigeonpea that implicated various candidate genes like sucrose-proton symporter 2 and an uncharacterized protein along with pectate lyase, pectinesterase inhibitors, L-ascorbate oxidase homolog, ATPase, β-galactosidase, polygalacturonase, and aldose 1-epimerase for pollen development of pigeonpea. The dataset presented here provides a rich genomic resource to improve understanding of CMS trait and its deployment in heterosis breeding in pigeonpea.

Uncloaking lncRNA-meditated gene expression as a potential regulator of CMS in cotton (Gossypium hirsutum L.)

Cytoplasmic male sterility is a well-proven mechanism for cotton hybrid production. Long non-coding RNAs belong to a class of transcriptional regulators that function in multiple biological processes. The cDNA libraries from the flower buds of the cotton CGMS, it's restorer (Rf) and maintainer lines were sequenced using high throughput NGS technique. A total of 1531 lncRNAs showed significant differential expression patterns between these three lines. Functional analysis of the co-expression network of lncRNA-mRNA using gene ontology vouchsafes that, lncRNAs play a crucial role in cytoplasmic male sterility and fertility restoration through pollen development, INO80 complex, development of anther wall tapetum, chromatin remodeling, and histone modification. Additionally, 94 lncRNAs were identified as putative precursors of 49 miRNAs. qRT-PCR affirms the concordance of expression pattern to RNA-seq data. These findings divulge the lncRNA driven miRNA-mediated regulation of gene expression profiling superintended for a better understanding of the CMS mechanisms of cotton.

Comparative RNA-Seq profiling of a resistant and susceptible peanut (Arachis hypogaea) genotypes in response to leaf rust infection caused by Puccinia arachidis

The goal of this study was to identify differentially expressed genes (DEGs) responsible for peanut plant (Arachis hypogaea) defence against Puccinia arachidis (causative agent of rust disease). Genes were identified using a high-throughput RNA-sequencing strategy. In total, 86,380,930 reads were generated from RNA-Seq data of two peanut genotypes, JL-24 (susceptible), and GPBD-4 (resistant). Gene Ontology (GO) and KEGG analysis of DEGs revealed essential genes and their pathways responsible for defence response to P. arachidis. DEGs uniquely upregulated in resistant genotype included pathogenesis-related (PR) proteins, MLO such as protein, ethylene-responsive factor, thaumatin, and F-box, whereas, other genes down-regulated in susceptible genotype were Caffeate O-methyltransferase, beta-glucosidase, and transcription factors (WRKY, bZIP, MYB). Moreover, various genes, such as Chitinase, Cytochrome P450, Glutathione S-transferase, and R genes such as NBS-LRR were highly up-regulated in the resistant genotype, indicating their involvement in the plant defence mechanism. RNA-Seq analysis data were validated by RT-qPCR using 15 primer sets derived from DEGs producing high correlation value (R ² = 0.82). A total of 4511 EST-SSRs were identified from the unigenes, which can be useful in evaluating genetic diversity among genotypes, QTL mapping, and plant variety improvement through marker-assisted breeding. These findings will help to understand the molecular defence mechanisms of the peanut plant in response to P. arachidis infection.