Expression profile analysis of cotton fiber secondary cell wall thickening stage.

To determine the genes associated with the fiber strength trait in cotton, three different cotton cultivars were selected: Sea Island cotton (Xinhai 32, with hyper-long fibers labeled as HL), and upland cotton (17-24, with long fibers labeled as L, and 62-33, with short fibers labeled as S). These cultivars were chosen to assess fiber samples with varying qualities. RNA-seq technology was used to analyze the expression profiles of cotton fibers at the secondary cell wall (SCW) thickening stage (20, 25, and 30 days post-anthesis (DPA)). The results showed that a large number of differentially expressed genes (DEGs) were obtained from the three assessed cotton cultivars at different stages of SCW development. For instance, at 20 DPA, Sea Island cotton (HL) had 6,215 and 5,364 DEGs compared to upland cotton 17-24 (L) and 62-33 (S), respectively. Meanwhile, there were 1,236 DEGs between two upland cotton cultivars, 17-24 (L) and 62-33 (S). Gene Ontology (GO) term enrichment identified 42 functions, including 20 biological processes, 11 cellular components, and 11 molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified several pathways involved in SCW synthesis and thickening, such as glycolysis/gluconeogenesis, galactose metabolism, propanoate metabolism, biosynthesis of unsaturated fatty acids pathway, valine, leucine and isoleucine degradation, fatty acid elongation pathways, and plant hormone signal transduction. Through the identification of shared DEGs, 46 DEGs were found to exhibit considerable expressional differences at different fiber stages from the three cotton cultivars. These shared DEGs have functions including REDOX enzymes, binding proteins, hydrolases (such as GDSL thioesterase), transferases, metalloproteins (cytochromatin-like genes), kinases, carbohydrates, and transcription factors (MYB and WRKY). Therefore, RT-qPCR was performed to verify the expression levels of nine of the 46 identified DEGs, an approach which demonstrated the reliability of RNA-seq data. Our results provided valuable molecular resources for clarifying the cell biology of SCW biosynthesis during fiber development in cotton.

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.

Analysis of PAT1 subfamily members in the GRAS family of upland cotton and functional characterization of GhSCL13-2A in Verticillium dahliae resistance.

KEY MESSAGE: GhSCL13-2A, a member of the PAT1 subfamily in the GRAS family, positively regulates cotton resistance to Verticillium dahliae by mediating the jasmonic acid and salicylic acid signaling pathways and accumulation of reactive oxygen species. Verticillium wilt (VW) is a devastating disease of upland cotton (Gossypium hirsutum) that is primarily caused by the soil-borne fungus Verticillium dahliae. Scarecrow-like (SCL) proteins are known to be involved in plant abiotic and biotic stress responses, but their roles in cotton defense responses are still unclear. In this study, a total of 25 GhPAT1 subfamily members in the GRAS family were identified in upland cotton. Gene organization and protein domain analysis showed that GhPAT1 members were highly conserved. GhPAT1 genes were widely expressed in various tissues and at multiple developmental stages, and they were responsive to jasmonic acid (JA), salicylic acid (SA), and ethylene (ET) signals. Furthermore, GhSCL13-2A was induced by V. dahliae infection. V. dahliae resistance was enhanced in Arabidopsis thaliana by ectopic overexpression of GhSCL13-2A, whereas cotton GhSCL13-2A knockdowns showed increased susceptibility. Levels of reactive oxygen species (ROS) and JA were also increased and SA content was decreased in GhSCL13-2A knockdowns. At the gene expression level, PR genes and SA signaling marker genes were down-regulated and JA signaling marker genes were upregulated in GhSCL13-2A knockdowns. GhSCL13-2A was shown to be localized to the cell membrane and the nucleus. Yeast two-hybrid and luciferase complementation assays indicated that GhSCL13-2A interacted with GhERF5. In Arabidopsis, V. dahliae resistance was enhanced by GhERF5 overexpression; in cotton, resistance was reduced in GhERF5 knockdowns. This study revealed a positive role of GhSCL13-2A in V. dahliae resistance, establishing it as a strong candidate gene for future breeding of V. dahliae-resistant cotton cultivars.

Identification, evolutionary analysis and functional diversification of RAV gene family in cotton (G. hirsutum L.).

MAIN CONCLUSION: Genome wide analysis, expression pattern analysis, and functional characterization of RAV genes highlight their roles in roots, stem development and hormonal response. RAV (Related to ABI3 and VP1) gene family members have been involved in tissues/organs growth and hormone signaling in various plant species. Here, we identified 247 RAVs from 12 different species with 33 RAV genes from G. hirsutum. Phylogenetic analysis classified RAV genes into four distinct groups. Analysis of gene structure showed that most GhRAVs lack introns. Motif distribution pattern and protein sequence logos indicated that GhRAV genes were highly conserved during the process of evolution. Promotor cis-acting elements revealed that promotor regions of GhRAV genes encode numerous elements related to plant growth, abiotic stresses and phytohormones. Chromosomal location information showed uneven distribution of 33 GhRAV genes on different chromosomes. Collinearity analysis identified 628 and 52 orthologous/ paralogous gene pairs in G. hirsutum and G. barbadense, respectively. Ka/Ks values indicated that GhRAV and GbRAV genes underwent strong purifying selection pressure. Selecton model and codon model selection revealed that GhRAV amino acids were under purifying selection and adaptive evolution exists among GhRAV proteins. Three dimensional structure of GhRAVs indicated the presence of numerous alpha helix and beta-barrels. Expression level revealed that some GhRAV genes exhibited high expression in roots (GhRAV3, GhRAV4, GhRAV11, GhRAV18, GhRAV20 and GhRAV30) and stem (GhRAV3 and GhRAV18), indicating their potential role in roots and stem development. GhRAV genes can be regulated by phytohormonal stresses (BL, JA and IAA). Our study provides a reference for future studies related to the functional analysis of GhRAVs in cotton.