Transcriptome analysis and identification of genes associated with oil accumulation in upland cotton

Integrative analyses of long and short-read RNA sequencing reveal the spliced isoform regulatory network of seedling growth dynamics in upland cotton

The polyploid genome of cotton has significantly increased the transcript complexity. Recent advances in full-length transcript sequencing are now widely used to characterize the complete landscape of transcriptional events. Such studies in cotton can help us to explore the genetic mechanisms of the cotton seedling growth. Through long-read single-molecule RNA sequencing, this study compared the transcriptomes of three yield contrasting genotypes of upland cotton. Our analysis identified different numbers of spliced isoforms from 31,166, 28,716, and 28,713 genes in SJ48, Z98, and DT8 cotton genotypes, respectively, most of which were novel compared to previous cotton reference transcriptomes, and showed significant differences in the number of exon structures and coding sequence length due to intron retention. Quantification of isoform expression revealed significant differences in expression in the root and leaf of each genotype. An array of key isoform target genes showed protein kinase or phosphorylation functions, and their protein interaction network contained most of the circadian oscillator proteins. Spliced isoforms from the GIGANTEA (GI) protien were differentially regulated in each genotype and might be expected to regulate translational activities, including the sequence and function of target proteins. In addition, these spliced isoforms generate diurnal expression profiles in cotton leaves, which may alter the transcriptional regulatory network of seedling growth. Silencing of the novel spliced GI isoform Gh_A02G0645_N17 significantly affected biomass traits, contributed to variable growth, and increased transcription of the early flowering pathway gene ELF in cotton. Our high-throughput hybrid sequencing results will be useful to dissect functional differences among spliced isoforms in the polyploid cotton genome.

Identification of key genes for triacylglycerol biosynthesis and storage in herbaceous peony (Paeonia lactifolra Pall.) seeds based on full-length transcriptome

BACKGROUND

The herbaceous peony (Paeonia lactiflora Pall.) is extensively cultivated in China due to its root being used as a traditional Chinese medicine known as 'Radix Paeoniae Alba'. In recent years, it has been discovered that its seeds incorporate abundant unsaturated fatty acids, thereby presenting a potential new oilseed plant. Surprisingly, little is known about the full-length transcriptome sequencing of Paeonia lactiflora, limiting research into its gene function and molecular mechanisms.

RESULTS

A total of 484,931 Reads of Inserts (ROI) sequences and 1,455,771 full-Length non-chimeric reads (FLNC) sequences were obtained for CDS prediction, TF analysis, SSR analysis and lncRNA identification. In addition, gene function annotation and gene structure analysis were performed. A total of 4905 transcripts were related to lipid metabolism biosynthesis pathway, belonging to 28 enzymes. We use these data to identify 10 oleosin (OLE) and 5 diacylglycerol acyltransferase (DGAT) gene members after de-redundancy. The analysis of physicochemical properties and secondary structure showed them similarity in gene family respectively. The phylogenetic analysis showed that the distribution of OLE and DGAT family members was roughly the same as that of Arabidopsis. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses revealed expression changes in different seed development stages, and showed a trend of increasing and then decreasing.

CONCLUSION

In summary, these results provide new insights into the molecular mechanism of triacylglycerol (TAG) biosynthesis and storage during the seedling stage in Paeonia lactiflora. It provides theoretical references for selecting and breeding oil varieties and understanding the functions of oil storage as well as lipid synthesis related genes in Paeonia lactiflora.

Genome-wide identification and expression-pattern analysis of sulfate transporter (SULTR) gene family in cotton under multiple abiotic stresses and fiber development

Sulfate transporter (SULTR) proteins are in charge of the transport and absorption on sulfate substances, and have been reported to play vital roles in the biological processes of plant growth and stress response. However, there were few reports of genome-wide identification and expression-pattern analysis of SULTRs in Hibiscus mutabilis. Gossypium genus is a ideal model for studying the allopolyploidy, therefore two diploid species (G. raimondii and G. arboreum) and two tetraploid species (G. hirsutum and G. barbadense) were chosen in this study to perform bioinformatic analyses, identifying 18, 18, 35, and 35 SULTR members, respectively. All the 106 cotton SULTR genes were utilized to construct the phylogenetic tree together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 8 Zea mays ones, which was divided into Group1-Group4. The clustering analyses of gene structures and 10 conserved motifs among the cotton SULTR genes showed the consistent evolutionary relationship with the phylogenetic tree, and the results of gene-duplication identification among the four representative Gossypium species indicated that genome-wide or segment duplication might make main contributions to the expansion of SULTR gene family in cotton. Having conducted the cis-regulatory element analysis in promoter region, we noticed that the existing salicylic acid (SA), jasmonic acid (JA), and abscisic acid (ABA) elements could have influences with expression levels of cotton SULTR genes. The expression patterns of GhSULTR genes were also investigated on the 7 different tissues or organs and the developing ovules and fibers, most of which were highly expressed in root, stem, sepal, receptacel, ovule at 10 DPA, and fiber at 20 and 25 DPA. In addition, more active regulatory were observed in GhSULTR genes responding to multiple abiotic stresses, and 12 highly expressed genes showed the similar expression patterns in the quantitative Real-time PCR experiments under cold, heat, salt, and drought treatments. These findings broaden our insight into the evolutionary relationships and expression patterns of the SULTR gene family in cotton, and provide the valuable information for further screening the vital candidate genes on trait improvement.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

Genomic and co-expression network analyses reveal candidate genes for oil accumulation based on an introgression population in Upland cotton (Gossypium hirsutum)

KEY MESSAGE

Integrated QTL mapping and WGCNA condense the potential gene regulatory network involved in oil accumulation. A glycosyl hydrolases gene (GhHSD1) for oil biosynthesis was confirmed in Arabidopsis, which will provide useful knowledge to understand the functional mechanism of oil biosynthesis in cotton. Cotton is an economical source of edible oil for the food industry. The genetic mechanism that regulates oil biosynthesis in cottonseeds is essential for the genetic enhancement of oil content (OC). To explore the functional genomics of OC, this study utilized an interspecific backcross inbred line population to dissect the quantitative trait locus (QTL) interlinked with OC. In total, nine OC QTLs were identified, four of which were novel, and each QTL explained 3.62-34.73% of the phenotypic variation of OC. The comprehensive transcript profiling of developing cottonseeds revealed 3,646 core genes differentially expressed in both inbred parents. Functional enrichment analysis determined 43 genes were annotated with oil biosynthesis processes. Implementation of weighted gene co-expression network analysis showed that 803 differential genes had a significant correlation with the OC phenotype. Further integrated analysis identified seven important genes located in OC QTLs. Of which, the GhHSD1 gene located in stable QTL qOC-Dt3-1 exhibited the highest functional linkages with the other network genes. Phylogenetic analysis showed significant evolutionary differences in the HSD1 sequences between oilseed- and starch- crops. Furthermore, the overexpression of GhHSD1 in Arabidopsis yielded almost 6.78% higher seed oil. This study not only uncovers important genetic loci for oil accumulation in cottonseed, but also provides a set of new candidate genes that potentially influence the oil biosynthesis pathway in cottonseed.

Omics-driven exploration and mining of key functional genes for the improvement of food and fiber crops

The deployment of omics technologies has obtained an incredible boost over the past few decades with the advances in next-generation sequencing (NGS) technologies, innovative bioinformatics tools, and the deluge of available biological information. The major omics technologies in the limelight are genomics, transcriptomics, proteomics, metabolomics, and phenomics. These biotechnological advances have modernized crop breeding and opened new horizons for developing crop varieties with improved traits. The genomes of several crop species are sequenced, and a huge number of genes associated with crucial economic traits have been identified. These identified genes not only provide insights into the understanding of regulatory mechanisms of crop traits but also decipher practical grounds to assist in the molecular breeding of crops. This review discusses the potential of omics technologies for the acquisition of biological information and mining of the genes associated with important agronomic traits in important food and fiber crops, such as wheat, rice, maize, potato, tomato, cassava, and cotton. Different functional genomics approaches for the validation of these important genes are also highlighted. Furthermore, a list of genes discovered by employing omics approaches is being represented as potential targets for genetic modifications by the latest genome engineering methods for the development of climate-resilient crops that would in turn provide great impetus to secure global food security.

Oil candidate genes in seeds of cotton (Gossypium hirsutum L.) and functional validation of GhPXN1

BACKGROUND

Cottonseed oil is a promising edible plant oil with abundant unsaturated fatty acids. However, few studies have been conducted to explore the characteristics of cottonseed oil. The molecular mechanism of cottonseed oil accumulation remains unclear.

RESULTS

In the present study, we conducted comparative transcriptome and weighted gene co-expression network (WGCNA) analysis for two G. hirsutum materials with significant difference in cottonseed oil content. Results showed that, between the high oil genotype 6053 (H6053) and the low oil genotype 2052 (L2052), a total of 412, 507, 1,121, 1,953, and 2,019 differentially expressed genes (DEGs) were detected at 10, 15, 20, 25, and 30 DPA, respectively. Remarkably, a large number of the down-regulated DEGs were enriched in the phenylalanine metabolic processes. Investigation into the dynamic changes of expression profiling of genes associated with both phenylalanine metabolism and oil biosynthesis has shed light on a significant competitive relationship in substrate allocation during cottonseed development. Additionally, the WGCNA analysis of all DEGs identified eight distinct modules, one of which includes GhPXN1, a gene closely associated with oil accumulation. Through phylogenetic analysis, we hypothesized that GhPXN1 in G. hirsutum might have been introgressed from G. arboreum. Overexpression of the GhPXN1 gene in tobacco leaf suggested a significant reduction in oil content compared to the empty-vector transformants. Furthermore, ten other crucial oil candidate genes identified in this study were also validated using quantitative real-time PCR (qRT-PCR).

CONCLUSIONS

Overall, this study enhances our comprehension of the molecular mechanisms underlying cottonseed oil accumulation.

Association Analysis of Tiller-Related Traits with EST-SSR Markers in Psathyrostachys juncea

Psathyrostachys juncea is a long-lived perennial Gramineae grass with dense basal tillers and soft leaves. It is used widely in cold and dry areas of Eurasia and North America to establish grazing pasture and is even used as an ideal plant for revegetation and ecological restoration. Plant architecture, especially tillering traits, is critical for bunch grasses in breeding programs, and these traits in plants are mostly quantitative traits. In this study, the genetic diversity, population structure, and linkage disequilibrium of 480 individual lines were analyzed using 127 pairs of the EST-SSR marker, and a significant association between ten plant-architecture-related traits of P. juncea and molecular markers was found. The results of the genetic diversity analysis showed that the number of observed alleles was 1.957, the number of effective alleles was 1.682, Shannon's information index was 0.554, observed heterozygosity was 0.353, expected heterozygosity was 0.379, and the polymorphism information content was 0.300. A total of 480 individual lines were clustered into five groups based on population genetic structure, principal coordinate analysis, and unweighted pair group method with arithmetic mean analysis (UPGMA). The linkage disequilibrium coefficient (r2) was between 0.00 and 0.68, with an average of 0.04, which indicated a relatively low level of linkage disequilibrium among loci. The results of the association analysis revealed 55 significant marker-trait associations (MTA). Moreover, nine SSR markers were associated with multiple traits. This study provides tools with promising applications in the molecular selection and breeding of P. juncea germplasm.

Genetic and Morpho-Physiological Differences among Transgenic and No-Transgenic Cotton Cultivars

Three carbon-chain extension genes associated with fatty acid synthesis in upland cotton (Gossypium hirsutum), namely GhKAR, GhHAD, and GhENR, play important roles in oil accumulation in cotton seeds. In the present study, these three genes were cloned and characterized. The expression patterns of GhKAR, GhHAD, and GhENR in the high seed oil content cultivars 10H1014 and 10H1041 differed somewhat compared with those of 10H1007 and 2074B with low seed oil content at different stages of seed development. GhKAR showed all three cultivars showed higher transcript levels than that of 2074B at 10-, 40-, and 45-days post anthesis (DPA). The expression pattern of GhHAD showed a lower transcript level than that of 2074B at both 10 and 30 DPA but a higher transcript level than that of 2074B at 40 DPA. GhENR showed a lower transcript level than that of 2074B at both 15 and 30 DPA. The highest transcript levels of GhKAR and GhENR were detected at 15 DPA in 10H1007, 10H1014, and 10H1041 compared with 2074B. From 5 to 45 DPA cotton seed, the oil content accumulated continuously in the developing seed. Oil accumulation reached a peak between 40 DPA and 45 DPA and slightly decreased in mature seed. In addition, GhKAR and GhENR showed different expression patterns in fiber and ovule development processes, in which they showed high expression levels at 20 DPA during the fiber elongation stage, but their expression level peaked at 15 DPA during ovule development processes. These two genes showed the lowest expression levels at the late seed maturation stage, while GhHAD showed a peak of 10 DPA in fiber development. Compared to 2074B, the oil contents of GhKAR and GhENR overexpression lines increased 1.05~1.08 folds. These results indicated that GhHAD, GhENR, and GhKAR were involved in both seed oil synthesis and fiber elongation with dual biological functions in cotton.

Genome-Wide Identification and Expression Analyses of the Cotton AGO Genes and Their Potential Roles in Fiber Development and Stress Response

Argonaute proteins (AGOs) are indispensable components of RNA silencing. However, systematic characterization of the AGO genes have not been completed in cotton until now. In this study, cotton AGO genes were identified and analyzed with respect to their evolution and expression profile during biotic and abiotic stresses. We identified 14 GaAGO, 14 GrAGO, and 28 GhAGO genes in the genomes of Gossypium arboreum, Gossypium raimondii, and Gossypium hirsutum. Cotton AGO proteins were classified into four subgroups. Structural and functional conservation were observed in the same subgroups based on the analysis of the gene structure and conserved domains. Twenty-four duplicated gene pairs were identified in GhAGO genes, and all of them exhibited strong purifying selection during evolution. Moreover, RNA-seq analysis showed that most of the GhAGO genes exhibit high expression levels in the fiber initiation and elongation processes. Furthermore, the expression profiles of GhAGO genes tested by quantitative real-time polymerase chain reaction (qPCR) demonstrated that they were sensitive to Verticillium wilt infection and salt and drought stresses. Overall, our results will pave the way for further functional investigation of the cotton AGO gene family, which may be involved in fiber development and stress response.