Co-expression network and comparative transcriptome analysis for fiber initiation and elongation reveal genetic differences in two lines from upland cotton CCRI70 RIL population.
Ontology highlight
ABSTRACT:
SUBMITTER: Jiang X
PROVIDER: S-EPMC8308610 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundImproving the yield and fiber quality of upland cotton is a goal of plant breeders. However, increasing the yield and quality of cotton fibers is becoming more urgent. While the growing human population needs more cotton fiber, climate change is reducing the amount of land on which cotton can be planted, or making it difficult to ensure that water and other resources will be available in optimal quantities. The most logical means of improving yield and quality is understanding and manipulating the genes involved. Here, we used comparative transcriptomics to explore differences in gene expression between long- and short-fiber cotton lines to identify candidate genes useful for cotton improvement.ResultsLight and electron microscopy revealed that the initial fiber density was significantly greater in our short-fiber group (SFG) than in our long-fiber group (LFG). Compared with the SFG fibers, the LFG fibers were longer at all developmental stages. Comparison of the LFG and SFG transcriptomes revealed a total of 3538 differentially expressed genes (DEGs). Notably, at all three developmental stages examined, two expression patterns, consistently downregulated (profile 0) and consistently upregulated (profile 7), were identified, and both were significantly enriched in the SFG and LFG. Twenty-two DEGs known to be involved in fiber initiation were detected in profile 0, while 31 DEGs involved in fiber elongation were detected in profile 7. Functional annotation suggested that these DEGs, which included ERF1, TUA2, TUB1, and PER64, affect fiber elongation by participating in the ethylene response, microtubule synthesis, and/or the peroxidase (POD) catalytic pathway. qRT-PCR was used to confirm the RNA sequencing results for select genes.ConclusionsA comparison of SFG and LFG transcription profiles revealed modest but important differences in gene expression between the groups. Notably, our results confirm those of previous studies suggesting that genes involved in ethylene, tubulin, and POD pathways play important roles in fiber development. The 22 consistently downregulated DEGs involved in fiber initiation and the 31 consistently upregulated genes involved in fiber elongation are seemingly good candidate genes for improving fiber initiation and elongation in cotton.
| S-EPMC6683361 | biostudies-literature
Project description:BackgroundHow to develop new cotton varieties possessing high yield traits of Upland cotton and superior fiber quality traits of Sea Island cotton remains a key task for cotton breeders and researchers. While multiple attempts bring in little significant progresses, the development of Chromosome Segment Substitution Lines (CSSLs) from Gossypium barbadense in G. hirsutum background provided ideal materials for aforementioned breeding purposes in upland cotton improvement. Based on the excellent fiber performance and relatively clear chromosome substitution segments information identified by Simple Sequence Repeat (SSR) markers, two CSSLs, MBI9915 and MBI9749, together with the recurrent parent CCRI36 were chosen to conduct transcriptome sequencing during the development stages of fiber elongation and Secondary Cell Wall (SCW) synthesis (from 10DPA and 28DPA), aiming at revealing the mechanism of fiber development and the potential contribution of chromosome substitution segments from Sea Island cotton to fiber development of Upland cotton.ResultsIn total, 15 RNA-seq libraries were constructed and sequenced separately, generating 705.433 million clean reads with mean GC content of 45.13% and average Q30 of 90.26%. Through multiple comparisons between libraries, 1801 differentially expressed genes (DEGs) were identified, of which the 902 up-regulated DEGs were mainly involved in cell wall organization and response to oxidative stress and auxin, while the 898 down-regulated ones participated in translation, regulation of transcription, DNA-templated and cytoplasmic translation based on GO annotation and KEGG enrichment analysis. Subsequently, STEM software was performed to explicate the temporal expression pattern of DEGs. Two peroxidases and four flavonoid pathway-related genes were identified in the "oxidation-reduction process", which could play a role in fiber development and quality formation. Finally, the reliability of RNA-seq data was validated by quantitative real-time PCR of randomly selected 20 genes.ConclusionsThe present report focuses on the similarities and differences of transcriptome profiles between the two CSSLs and the recurrent parent CCRI36 and provides novel insights into the molecular mechanism of fiber development, and into further exploration of the feasible contribution of G. barbadense substitution segments to fiber quality formation, which will lay solid foundation for simultaneously improving fiber yield and quality of upland cotton through CSSLs.
| S-EPMC5591532 | biostudies-literature
Project description:To identify potential miRNAs involved in fiber development and elucidate their expression differences between G. barbadense and G. hirsutum, we constructed two small RNA libraries, Gb10 and Gh10, prepared from fibers of 3-79 (G. barbadense) and TM-1 (G. hirsutum) collected at 10 days post-anthesis (DPA). We identified 28 conserved miRNA families, including 24 that exactly match known plant miRNA families in miRBase. With MIREAP and newly developed software miRsearcher, 7 candidate-novel miRNAs were found. 5 candidate-novel miRNAs were expressed in both species, 2 candidate-novel miRNAs were expressed only in one species. Moreover, 4 miRNA families showed significant expression differences between sea-island cotton and upland cotton in 10 DPA fibers.
2010-12-03 | GSE25773 | GEO
Project description:To identify potential miRNAs involved in fiber development and elucidate their expression differences between G. barbadense and G. hirsutum, we constructed two small RNA libraries, Gb10 and Gh10, prepared from fibers of 3-79 (G. barbadense) and TM-1 (G. hirsutum) collected at 10 days post-anthesis (DPA). We identified 28 conserved miRNA families, including 24 that exactly match known plant miRNA families in miRBase. With MIREAP and newly developed software miRsearcher, 7 candidate-novel miRNAs were found. 5 candidate-novel miRNAs were expressed in both species, 2 candidate-novel miRNAs were expressed only in one species. Moreover, 4 miRNA families showed significant expression differences between sea-island cotton and upland cotton in 10 DPA fibers. two examples including 3-79 and TM-1 10 DPA fibers
2010-12-03 | E-GEOD-25773 | biostudies-arrayexpress
Project description:Histone methylation is one of the most significant epigenetic modifications in plants because it is responsible for regulating and controlling the expression of many important functional genes. Although there have been many studies concerning histone methylation, its specific functions and mechanisms in cotton fiber initiation remain unclear. To determine if the H3K27me3 modification had regulated cotton fiber initiation in upland cotton (Gossypium hirsutum), we performed a comparative analysis on the cleavage under targets and tagmentation (CUT&Tag) between the wild-type cultivar Xu142 and its fibreless mutant Xu142 fl. This research was carried out three time points during the initiation period of cotton fibers, which were -1, 0 and 1 days post-anthesis (DPA). We also used the H3K27me3 antibody in Xu142 and Xu142 fl to examine the differential deposition of H3K27me3. Combing these results with RNA-seq analysis, we discovered a negative correlation between transcription and H3K27me3. We identified 336 genes with opposing trends in H3K27me3 deposition and expression. These genes are involved in various metabolic pathways, including fatty acid biosynthesis, galactose metabolism, plant hormone signal transduction, plant-pathogen interaction, glycolysis-/-gluconeogenesis, phenylpropanoid biosynthesis, and cutin, suberin and wax biosynthesis. Further physiological research revealed that in vitro ovule culture with an RDS 3434 inhibitor could alters the level of H3K27me3 in ovules, thus influencing the process of cotton fiber development. Taken together, these results indicate that H3K27me3 plays an important role in the fiber initiation process of cotton, and demonstrates the H3K27me3 network involved in fiber initiation in cotton ovules. The study also provides a foundation for the future study of genes involved in fiber development.
2024-03-14 | GSE227870 | GEO
Project description:BackgroundCotton is the most essential textile crop worldwide, and phytohormones are critical for cotton fiber development. One example is the role of auxin in fiber initiation, but we know little molecular basis. MicroRNAs (miRNAs) have a significant function in cotton development; nevertheless their role in fiber initiation remains unclear. Here, exogenous IAA was applied to cotton plant before anthesis. Utilizing small RNA sequencing, the mechanism underlying miRNA-mediated regulation of fiber initiation under exogenous IAA treatment was investigated.ResultsWith exogenous IAA application, the endogenous IAA and GA contents of IAA treated (IT) ovules were higher than control (CK) ovules at the fiber initiation stage, while endogenous ABA content was lower in IT than CK. Using scanning electron microscopy, we found the fiber number and size were significantly promoted in IT at 0 DPA. Fiber quality analysis showed that fiber length, uniformity, strength, elongation, and micronaire of IT were higher than CK, though not statistically significant, while lint percent was significantly higher in IT. We generated six small RNA libraries using - 3, 0, and 3 DPA ovules of IT and CK, and identified 58 known miRNAs and 83 novel miRNAs together with the target genes. The differential expressed miRNAs number between IT and CK at - 3, 0, 3 DPA was 34, 16 and 24, respectively. Gene ontology and KEGG pathway enrichment analyses for the target genes of the miRNAs expressed in a differential manner showed that they were significantly enriched in 30 terms and 8 pathways. QRT-PCR for those identified miRNAs and the target genes related to phytohormones and fiber development was performed, and results suggested a potential role of these miRNAs in fiber initiation.ConclusionsThe exogenous IAA application affected the relative phytohormone contents in ovule and promoted fiber initiation in cotton. Identification and profiling of miRNAs and their targets at the fiber initiation stage provided insights for miRNAs' regulation function of fiber initiation. These findings not only shed light on the regulatory network of fiber growth but also offer clues for cotton fiber amelioration strategies in cotton.
| S-EPMC6537205 | biostudies-literature
Project description:BackgroundAppropriate plant architecture can improve the amount of cotton boll opening and allow increased planting density, thus increasing the level of cotton mechanical harvesting and cotton yields. The internodes of cotton fruiting branches are an important part of cotton plant architecture. Thus, studying the molecular mechanism of internode elongation in cotton fruiting branches is highly important.ResultsIn this study, we selected internodes of cotton fruiting branches at three different stages from two cultivars whose internode lengths differed significantly. A total of 76,331 genes were detected by transcriptome sequencing. By KEGG pathway analysis, we found that DEGs were significantly enriched in the plant hormone signal transduction pathway. The transcriptional data and qRT-PCR results showed that members of the GH3 gene family, which are involved in auxin signal transduction, and CKX enzymes, which can reduce the level of CKs, were highly expressed in the cultivar XLZ77, which has relatively short internodes. Genes related to ethylene synthase (ACS), EIN2/3 and ERF in the ethylene signal transduction pathway and genes related to JAR1, COI1 and MYC2 in the JA signal transduction pathway were also highly expressed in XLZ77. Plant hormone determination results showed that the IAA and CK contents significantly decreased in cultivar XLZ77 compared with those in cultivar L28, while the ACC (the precursor of ethylene) and JA contents significantly increased. GO enrichment analysis revealed that the GO categories associated with promoting cell elongation, such as cell division, the cell cycle process and cell wall organization, were significantly enriched, and related genes were highly expressed in L28. However, genes related to the sphingolipid metabolic process and lignin biosynthetic process, whose expression can affect cell elongation, were highly expressed in XLZ77. In addition, 2067 TFs were differentially expressed. The WRKY, ERF and bHLH TF families were the top three largest families whose members were active in the two varieties, and the expression levels of most of the genes encoding these TFs were upregulated in XLZ77.ConclusionsAuxin and CK are positive regulators of internode elongation in cotton branches. In contrast, ethylene and JA may act as negative regulators of internode elongation in cotton branches. Furthermore, the WRKY, ERF and bHLH TFs were identified as important inhibitors of internode elongation in cotton. In XLZ77(a short-internode variety), the mass synthesis of ethylene and amino acid conjugation of auxin led to the inhibition of plant cell elongation, while an increase in JA content and degradation of CKs led to a slow rate of cell division, which eventually resulted in a phenotype that presented relatively short internodes on the fruiting branches. The results of this study not only provide gene resources for the genetic improvement of cotton plant architecture but also lay a foundation for improved understanding of the molecular mechanism of the internode elongation of cotton branches.
| S-EPMC6781417 | biostudies-literature
Project description:BackgroundIn upland cotton (Gossypium hirsutum L.), genotypes with the same mature fiber length (FL) might possess different genes and exhibit differential expression of genes related to fiber elongation at different fiber developmental stages. However, there is a lack of information on the genetic variation influencing fiber length and its quantitative trait loci (QTLs) during the fiber elongation stage. In this study, a subset of upland cotton accessions was selected based on a previous GWAS conducted in China and grown in multiple environments to determine the dynamic fiber length at 10, 15, 20, and 25 days post-anthesis (DPA) and maturity. The germplasm lines were genotyped with the Cotton 63 K Illumina single-nucleotide polymorphism (SNP) array for GWAS.ResultsA total of 25, 38, 57, 89 and 88 SNPs showed significant correlations with fiber length at 10, 15, 20 and 25 DPA and maturity, respectively. In addition, 60 more promising SNPs were detected in at least two tests and two FL developmental time points, and 20 SNPs were located within the confidence intervals of QTLs identified in previous studies. The fastest fiber-length growth rates were obtained at 10 to 15 DPA in 69 upland cotton lines and at 15 to 20 DPA in 14 upland cotton accessions, and 10 SNPs showed significant correlations with the fiber-length growth rate. A combined transcriptome and qRT-PCR analysis revealed that two genes (D10G1008 and D13G2037) showed differential expression between two long-fiber genotypes and two short-fiber genotypes.ConclusionsThis study provides important new insights into the genetic basis of the time-dependent fiber-length trait and reveals candidate SNPs and genes for improving fiber length in upland cotton.
| S-EPMC6282333 | biostudies-literature
Project description:Fiber strength is the key trait that determines fiber quality in cotton, and it is closely related to secondary cell wall synthesis. To understand the mechanism underlying fiber strength, we compared fiber transcriptomes from different G. barbadense chromosome introgression lines (CSILs) that had higher fiber strengths than their recipient, G. hirsutum acc. TM-1. A total of 18,288 differentially expressed genes (DEGs) were detected between CSIL-35431 and CSIL-31010, two CSILs with stronger fiber and TM-1 during secondary cell wall synthesis. Functional classification and enrichment analysis revealed that these DEGs were enriched for secondary cell wall biogenesis, glucuronoxylan biosynthesis, cellulose biosynthesis, sugar-mediated signaling pathways, and fatty acid biosynthesis. Pathway analysis showed that these DEGs participated in starch and sucrose metabolism (328 genes), glycolysis/gluconeogenesis (122 genes), phenylpropanoid biosynthesis (101 genes), and oxidative phosphorylation (87 genes), etc. Moreover, the expression of MYB- and NAC-type transcription factor genes were also dramatically different between the CSILs and TM-1. Being different to those of CSIL-31134, CSIL-35431 and CSIL-31010, there were many genes for fatty acid degradation and biosynthesis, and also for carbohydrate metabolism that were down-regulated in CSIL-35368. Metabolic pathway analysis in the CSILs showed that different pathways were changed, and some changes at the same developmental stage in some pathways. Our results extended our understanding that carbonhydrate metabolic pathway and secondary cell wall biosynthesis can affect the fiber strength and suggested more genes and/or pathways be related to complex fiber strength formation process.
| S-EPMC3998979 | biostudies-literature
Project description:Two fiber tissues harvested 10 days post anthesis from upland cotton trees grown under the same green house conditions except for different seasons of the year were used for RNA extraction. Small RNA molecules under 30 bases were amplified and isolated from an agarose gel. The purified DNA was used directly for cluster generation and sequencing analysis using the Illumina Genome Analyzer according to the manufacturer's instructions. The 35nt sequence tags from sequencing went through data cleaning first, which included getting rid of the low-quality tags and several kinds of contaminants from the 35nt tags. All clean tag sequences with copy numbers were then summarized into a fasta format file.
2011-07-01 | GSE27697 | GEO