GhLTPG1, a cotton GPI-anchored lipid transfer protein, regulates the transport of phosphatidylinositol monophosphates and cotton fiber elongation

An aldehyde dehydrogenase gene, GhALDH7B4_A06, positively regulates fiber strength in upland cotton (Gossypium hirsutum L.)

High fiber strength (FS) premium cotton has significant market demand. Consequently, enhancing FS is a major objective in breeding quality cotton. However, there is a notable lack of known functionally applicable genes that can be targeted for breeding. To address this issue, our study used specific length-amplified fragment sequencing combined with bulk segregant analysis to study FS trait in an F2 population. Subsequently, we integrated these results with previous quantitative trait locus mapping results regarding fiber quality, which used simple sequence repeat markers in F2, F2:3, and recombinant inbred line populations. We identified a stable quantitative trait locus qFSA06 associated with FS located on chromosome A06 (90.74-90.83 Mb). Within this interval, we cloned a gene, GhALDH7B4_A06, which harbored a critical mutation site in coding sequences that is distinct in the two parents of the tested cotton line. In the paternal parent Ji228, the gene is normal and referred to as GhALDH7B4_A06O; however, there is a nonsense mutation in the maternal parent Ji567 that results in premature termination of protein translation, and this gene is designated as truncated GhALDH7B4_A06S. Validation using recombinant inbred lines and gene expression analysis revealed that this mutation site is correlated with cotton FS. Virus-induced gene silencing of GhALDH7B4 in cotton caused significant decreases in FS and fiber micronaire. Conversely, GhALDH7B4_A06O overexpression in Arabidopsis boosted cell wall component contents in the stem. The findings of our study provide a candidate gene for improving cotton fiber quality through molecular breeding.

Characterization and expression analysis of GLABRA3 (GL3) genes in cotton: insights into trichome development and hormonal regulation

BACKGROUND

GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) genes encode a typical helix-loop-helix (bHLH) transcription factors that primarily regulate trichome branching and root hair development, DNA endoreduplication, trichoblast size, and stomatal formation. The functions of GL3 genes in cotton crop have been poorly characterized. In this study, we performed comprehensive genome-wide scans for GL3 and EGL3 homologs to enhance our comprehension of their potential roles in trichome and fiber development in cotton crop.

METHODS AND RESULTS

Our findings paraded that Gossypium hirsutum and G. barbadense have 6 GL3s each, unevenly distributed on 4 chromosomes whereas, G. arboreum, and G. raimondii have 3 GL3s each, unevenly distributed on 2 chromosomes. Gh_A08G2088 and Gb_A09G2187, despite having the same bHLH domain as the other GL3 genes, were excluded due to remarkable short sequences and limited number of motifs, indicating a lack of potential functional activity. The phylogenetic analysis categorized remaining 16 GL3s into three subfamilies (Group I-III) closely related to A. thaliana. The 16 GL3s have complete bHLH domain, encompassing 590-631 amino acids, with molecular weights (MWs) ranging from 65.92 to 71.36 kDa. Within each subfamily GL3s depicted shared similar gene structures and motifs, indicating conserved characteristics within respective groups. Promoter region analysis revealed 27 cis-acting elements, these elements were responsive to salicylic acid, abscisic acid (ABA), methyl jasmonate (MeJA), and gibberellin. The expression of GL3 genes was analyzed across 12 tissues in both G. barbadense and G. hirsutum using the publicly available RNA-seq data. Among GL3s, Gb_D11G0219, Gb_D11G0214, and Gb_D08G2182, were identified as relatively highly expressed across different tissues, consequently selected for hormone treatment and expression validation in G. barbadense. RT-qPCR results demonstrated significant alterations in the expression levels of Gb_D11G0219 and Gb_D11G0214 following MeJA, GA, and ABA treatment. Subcellular localization prediction revealed that most GL3 proteins were predominantly expressed in the nucleus, while a few were localized in the cytoplasm and chloroplasts.

CONCLUSIONS

In summary, this study lays the foundation for subsequent functional validation of GL3 genes by identifying hormonal regulation patterns and probable sites of action in cotton trichome formation and fiber development. The results stipulate a rationale to elucidate the roles and regulatory mechanisms of GL3 genes in the intricate process of cotton fibre and trichome development.

LIPID TRANSFER PROTEIN4 regulates cotton ceramide content and activates fiber cell elongation

Cell elongation is a fundamental process for plant growth and development. Studies have shown lipid metabolism plays important role in cell elongation; however, the related functional mechanisms remain largely unknown. Here, we report that cotton (Gossypium hirsutum) LIPID TRANSFER PROTEIN4 (GhLTP4) promotes fiber cell elongation via elevating ceramides (Cers) content and activating auxin-responsive pathways. GhLTP4 was preferentially expressed in elongating fibers. Over-expression and down-regulation of GhLTP4 led to longer and shorter fiber cells, respectively. Cers were greatly enriched in GhLTP4-overexpressing lines and decreased dramatically in GhLTP4 down-regulating lines. Moreover, auxin content and transcript levels of indole-3-acetic acid (IAA)-responsive genes were significantly increased in GhLTP4-overexpressing cotton fibers. Exogenous application of Cers promoted fiber elongation, while NPA (N-1-naphthalic acid, a polar auxin transport inhibitor) counteracted the promoting effect, suggesting that IAA functions downstream of Cers in regulating fiber elongation. Furthermore, we identified a basic helix-loop-helix transcription factor, GhbHLH105, that binds to the E-box element in the GhLTP4 promoter region and promotes the expression of GhLTP4. Suppression of GhbHLH105 in cotton reduced the transcripts level of GhLTP4, resulting in smaller cotton bolls and decreased fiber length. These results provide insights into the complex interactions between lipids and auxin-signaling pathways to promote plant cell elongation.

Domestication over Speciation in Allopolyploid Cotton Species: A Stronger Transcriptomic Pull

Cotton has been domesticated independently four times for its fiber, but the genomic targets of selection during each domestication event are mostly unknown. Comparative analysis of the transcriptome during cotton fiber development in wild and cultivated materials holds promise for revealing how independent domestications led to the superficially similar modern cotton fiber phenotype in upland (G. hirsutum) and Pima (G. barbadense) cotton cultivars. Here we examined the fiber transcriptomes of both wild and domesticated G. hirsutum and G. barbadense to compare the effects of speciation versus domestication, performing differential gene expression analysis and coexpression network analysis at four developmental timepoints (5, 10, 15, or 20 days after flowering) spanning primary and secondary wall synthesis. These analyses revealed extensive differential expression between species, timepoints, domestication states, and particularly the intersection of domestication and species. Differential expression was higher when comparing domesticated accessions of the two species than between the wild, indicating that domestication had a greater impact on the transcriptome than speciation. Network analysis showed significant interspecific differences in coexpression network topology, module membership, and connectivity. Despite these differences, some modules or module functions were subject to parallel domestication in both species. Taken together, these results indicate that independent domestication led G. hirsutum and G. barbadense down unique pathways but that it also leveraged similar modules of coexpression to arrive at similar domesticated phenotypes.

Recent progression and future perspectives in cotton genomic breeding

Upland cotton is an important global cash crop for its long seed fibers and high edible oil and protein content. Progress in cotton genomics promotes the advancement of cotton genetics, evolutionary studies, functional genetics, and breeding, and has ushered cotton research and breeding into a new era. Here, we summarize high-impact genomics studies for cotton from the last 10 years. The diploid Gossypium arboreum and allotetraploid Gossypium hirsutum are the main focus of most genetic and genomic studies. We next review recent progress in cotton molecular biology and genetics, which builds on cotton genome sequencing efforts, population studies, and functional genomics, to provide insights into the mechanisms shaping abiotic and biotic stress tolerance, plant architecture, seed oil content, and fiber development. We also suggest the application of novel technologies and strategies to facilitate genome-based crop breeding. Explosive growth in the amount of novel genomic data, identified genes, gene modules, and pathways is now enabling researchers to utilize multidisciplinary genomics-enabled breeding strategies to cultivate "super cotton", synergistically improving multiple traits. These strategies must rise to meet urgent demands for a sustainable cotton industry.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes for Fuzz Development in Gossypium hirsutum

Fuzzless Gossypium hirsutum mutants are ideal materials for investigating cotton fiber initiation and development. In this study, we used the fuzzless G. hirsutum mutant Xinluzao 50 FLM as the research material and combined it with other fuzzless materials for verification by RNA sequencing to explore the gene expression patterns and differences between genes in upland cotton during the fuzz period. A gene ontology (GO) enrichment analysis showed that differentially expressed genes (DEGs) were mainly enriched in the metabolic process, microtubule binding, and other pathways. A weighted gene co-expression network analysis (WGCNA) showed that two modules of Xinluzao 50 and Xinluzao 50 FLM and four modules of CSS386 and Sicala V-2 were highly correlated with fuzz. We selected the hub gene with the highest KME value among the six modules and constructed an interaction network. In addition, we selected some genes with high KME values from the six modules that were highly associated with fuzz in the four materials and found 19 common differential genes produced by the four materials. These 19 genes are likely involved in the formation of fuzz in upland cotton. Several hub genes belong to the arabinogalactan protein and GDSL lipase, which play important roles in fiber development. According to the differences in expression level, 4 genes were selected from the 19 genes and tested for their expression level in some fuzzless materials. The modules, hub genes, and common genes identified in this study can provide new insights into the formation of fiber and fuzz, and provide a reference for molecular design breeding for the genetic improvement of cotton fiber.

Individual lipid transfer proteins from Tanacetum parthenium show different specificity for extracellular accumulation of sesquiterpenes

A highly specialized function for individual LTPs for different products from the same terpenoid biosynthesis pathway is described and the function of an LTP GPI anchor is studied. Sequiterpenes produced in glandular trichomes of the medicinal plant Tanacetum parthenium (feverfew) accumulate in the subcuticular extracellular space. Transport of these compounds over the plasma membrane is presumably by specialized membrane transporters, but it is still not clear how these hydrophobic compounds are subsequently transported over the hydrophilic cell wall. Here we identified eight so-called non-specific Lipid transfer proteins (nsLTPs) genes that are expressed in feverfew trichomes. A putative function of these eight nsLTPs in transport of the lipophilic sesquiterpene lactones produced in feverfew trichomes, was tested in an in-planta transport assay using transient expression in Nicotiana benthamiana. Of eight feverfew nsLTP candidate genes analyzed, two (TpLTP1 and TpLTP2) can specifically improve extracellular accumulation of the sesquiterpene costunolide, while one nsLTP (TpLTP3) shows high specificity towards export of parthenolide. The specificity of the nsLTPs was also tested in an assay that test for the exclusion capacity of the nsLTP for influx of extracellular substrates. In such assay, TpLTP3 was identified as most effective in blocking influx of both costunolide and parthenolide, when these substrates are infiltrated into the apoplast. The TpLTP3 is special in having a GPI-anchor domain, which is essential for the export activity of TpLTP3. However, addition of the TpLTP3 GPI-anchor domain to TpLTP1 resulted in loss of TpLTP1 export activity. These novel export and exclusion assays thus provide new means to test functionality of plant nsLTPs.

Lipid transfer proteins involved in plant-pathogen interactions and their molecular mechanisms

Nonspecific lipid transfer proteins (LTPs) are small, cysteine-rich proteins that play numerous functional roles in plant growth and development, including cutin wax formation, pollen tube adhesion, cell expansion, seed development, germination, and adaptation to changing environmental conditions. LTPs contain eight conserved cysteine residues and a hydrophobic cavity that provides a wide variety of lipid-binding specificities. As members of the pathogenesis-related protein 14 family (PR14), many LTPs inhibit fungal or bacterial growth, and act as positive regulators in plant disease resistance. Over the past decade, these essential immunity-related roles of LTPs in plant immune processes have been documented in a growing body of literature. In this review, we summarize the roles of LTPs in plant-pathogen interactions, emphasizing the underlying molecular mechanisms in plant immune responses and specific LTP functions.

Analyzing Quantitative Trait Loci for Fiber Quality and Yield-Related Traits From a Recombinant Inbred Line Population With Gossypium hirsutum Race palmeri as One Parent

Fiber quality and yield-related traits are important agronomic traits in cotton breeding. To detect the genetic basis of fiber quality and yield related traits, a recombinant inbred line (RIL) population consisting of 182 lines was established from a cross between Gossypium hirsutum cultivar CCRI35 and G. hirsutum race palmeri accession TX-832. The RIL population was deeply genotyped using SLAF-seq and was phenotyped in six environments. A high-density genetic linkage map with 15,765 SNP markers and 153 SSR markers was constructed, with an average distance of 0.30 cM between adjacent markers. A total of 210 fiber quality quantitative trait loci (QTLs) and 73 yield-related QTLs were identified. Of the detected QTLs, 62 fiber quality QTLs and 10 yield-related QTLs were stable across multiple environments. Twelve and twenty QTL clusters were detected on the At and Dt subgenome, respectively. Twenty-three major QTL clusters were further validated through associated analysis and five candidate genes of four stable fiber quality QTLs were identified. This study revealed elite loci influencing fiber quality and yield and significant phenotypic selection regions during G. hirsutum domestication, and set a stage for future utilization of molecular marker assisted breeding in cotton breeding programs.

Novel insights into water-deficit-responsive mRNAs and lncRNAs during fiber development in Gossypium hirsutum

BACKGROUND

The fiber yield and quality of cotton are greatly and periodically affected by water deficit. However, the molecular mechanism of the water deficit response in cotton fiber cells has not been fully elucidated.

RESULTS

In this study, water deficit caused a significant reduction in fiber length, strength, and elongation rate but a dramatic increase in micronaire value. To explore genome-wide transcriptional changes, fibers from cotton plants subjected to water deficit (WD) and normal irrigation (NI) during fiber development were analyzed by transcriptome sequencing. Analysis showed that 3427 mRNAs and 1021 long noncoding RNAs (lncRNAs) from fibers were differentially expressed between WD and NI plants. The maximum number of differentially expressed genes (DEGs) and lncRNAs (DERs) was identified in fibers at the secondary cell wall biosynthesis stage, suggesting that this is a critical period in response to water deficit. Twelve genes in cotton fiber were differentially and persistently expressed at ≥ five time points, suggesting that these genes are involved in both fiber development and the water-deficit response and could potentially be used in breeding to improve cotton resistance to drought stress. A total of 540 DEGs were predicted to be potentially regulated by DERs by analysis of coexpression and genomic colocation, accounting for approximately 15.76% of all DEGs. Four DERs, potentially acting as target mimics for microRNAs (miRNAs), indirectly regulated their corresponding DEGs in response to water deficit.

CONCLUSIONS

This work provides a comprehensive transcriptome analysis of fiber cells and a set of protein-coding genes and lncRNAs implicated in the cotton response to water deficit, significantly affecting fiber quality during the fiber development stage.

Transcriptional loss of domestication-driven cytoskeletal GhPRF1 gene causes defective floral and fiber development in cotton (Gossypium)

Constitutive- and fiber-specific RNAi of GhPRF1 gene illustrated strong correlation between domestication-driven profilin genes and floral/fiber architecture in cotton. During morpho-transformation of short-fuzz of wild cotton into the elongating spinnable fibers under the millennia of human selection, actin-polymerizing cytoskeletal profilin genes had undergone significant sequence alterations and spatiotemporal shift in their transcription levels. To comprehend the expression dynamics of profilin genes with their phenotypic implications, transgenic expression modulation of cotton profilin 1 (GhPRF1) gene was performed in the constitutive- and fiber-specific manner in Coker 310FR cotton cultivar. The constitutive GhPRF1-RNAi lines (35S:GhPRF1-RNAi) exhibited distorted 'monadelphous' staminal-tube, reduced pollen-viability and poorly developed fibers, whereas floral and fiber development of fiber-specific GhPRF1-RNAi lines showed no abnormalities. Moreover, the fiber-specific GhPRF1 overexpression lines (FBP7:GhPRF1-Ox) showed increased emergence of fiber-initials on the ovule surface, on the contrary to no fiber-initials in fiber-specific RNAi lines (FBP7:GhPRF1-RNAi). Interestingly, the average seed weight and fiber weight of FBP7:GhPRF1-Ox lines increased > 60% and > 38%, respectively, compared with FBP7:GhPRF1-RNAi lines and untransformed control seeds. On a molecular basis, the aberrant floral and fiber development of 35S:GhPRF1-RNAi lines was largely associated with sugar metabolism and hormone-signaling mechanisms. These observations illustrated the strong correlation between domestication-driven GhPRF genes, and floral/fiber development in cotton. Also, the enhanced agronomic traits in GhPRF1-Ox lines of cotton empowered us to recognize their imperative roles, and their future deployment for the sustainable cotton crop improvement.

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.

Comparative Transcriptomics of Non-Embryogenic and Embryogenic Callus in Semi-Recalcitrant and Non-Recalcitrant Upland Cotton Lines

Somatic embryogenesis-mediated plant regeneration is essential for the genetic manipulation of agronomically important traits in upland cotton. Genotype specific recalcitrance to regeneration is a primary challenge in deploying genome editing and incorporating useful transgenes into elite cotton germplasm. In this study, transcriptomes of a semi-recalcitrant cotton (Gossypium hirsutum L.) genotype ‘Coker312’ were analyzed at two critical stages of somatic embryogenesis that include non-embryogenic callus (NEC) and embryogenic callus (EC) cells, and the results were compared to a non-recalcitrant genotype ‘Jin668’. We discovered 305 differentially expressed genes in Coker312, whereas, in Jin668, about 6-fold more genes (2155) were differentially expressed. A total of 154 differentially expressed genes were common between the two genotypes. Gene enrichment analysis of the upregulated genes identified functional categories, such as lipid transport, embryo development, regulation of transcription, sugar transport, and vitamin biosynthesis, among others. In Coker312 EC cells, five major transcription factors were highly upregulated: LEAFY COTYLEDON 1 (LEC1), WUS-related homeobox 5 (WOX5), ABSCISIC ACID INSENSITIVE3 (ABI3), FUSCA3 (FUS3), and WRKY2. In Jin668, LEC1, BABY BOOM (BBM), FUS3, and AGAMOUS-LIKE15 (AGL15) were highly expressed in EC cells. We also found that gene expression of these embryogenesis genes was typically higher in Jin668 when compared to Coker312. We conclude that significant differences in the expression of the above genes between Coker312 and Jin668 may be a critical factor affecting the regenerative ability of these genotypes.

Advances about the Roles of Membranes in Cotton Fiber Development

Cotton fiber is an extremely elongated single cell derived from the ovule epidermis and is an ideal model for studying cell development. The plasma membrane is tremendously expanded and accompanied by the coordination of various physiological and biochemical activities on the membrane, one of the three major systems of a eukaryotic cell. This review compiles the recent progress and advances for the roles of the membrane in cotton fiber development: the functions of membrane lipids, especially the fatty acids, sphingolipids, and phytosterols; membrane channels, including aquaporins, the ATP-binding cassette (ABC) transporters, vacuolar invertase, and plasmodesmata; and the regulation mechanism of membrane proteins, such as membrane binding enzymes, annexins, and receptor-like kinases.

Differential expression of SlKLUH controlling fruit and seed weight is associated with changes in lipid metabolism and photosynthesis-related genes

The sizes of plant organs such as fruit and seed are crucial yield components. Tomato KLUH underlies the locus fw3.2, an important regulator of fruit and seed weight. However, the mechanism by which the expression levels of KLUH affect organ size is poorly understood. We found that higher expression of SlKLUH increased cell proliferation in the pericarp within 5 d post-anthesis in tomato near-isogenic lines. Differential gene expression analyses showed that lower expression of SlKLUH was associated with increased expression of genes involved in lipid metabolism. Lipidomic analysis revealed that repression of SlKLUH mainly increased the contents of certain non-phosphorus glycerolipids and phospholipids and decreased the contents of four unknown lipids. Co-expression network analyses revealed that lipid metabolism was possibly associated with but not directly controlled by SlKLUH, and that this gene instead controls photosynthesis-related processes. In addition, many transcription factors putatively involved in the KLUH pathway were identified. Collectively, we show that SlKLUH regulates fruit and seed weight which is associated with altered lipid metabolism. The results expand our understanding of fruit and seed weight regulation and offer a valuable resource for functional studies of candidate genes putatively involved in regulation of organ size in tomato and other crops.

Genetic variants and underlying mechanisms influencing variance heterogeneity in maize

Traditional genetic studies focus on identifying genetic variants associated with the mean difference in a quantitative trait. Because genetic variants also influence phenotypic variation via heterogeneity, we conducted a variance-heterogeneity genome-wide association study to examine the contribution of variance heterogeneity to oil-related quantitative traits. We identified 79 unique variance-controlling single nucleotide polymorphisms (vSNPs) from the sequences of 77 candidate variance-heterogeneity genes for 21 oil-related traits using the Levene test (P < 1.0 × 10^-5 ). About 30% of the candidate genes encode enzymes that work in lipid metabolic pathways, most of which define clear expression variance quantitative trait loci. Of the vSNPs specifically associated with the genetic variance heterogeneity of oil concentration, 89% can be explained by additional linked mean-effects genetic variants. Furthermore, we demonstrated that gene × gene interactions play important roles in the formation of variance heterogeneity for fatty acid compositional traits. The interaction pattern was validated for one gene pair (GRMZM2G035341 and GRMZM2G152328) using yeast two-hybrid and bimolecular fluorescent complementation analyses. Our findings have implications for uncovering the genetic basis of hidden additive genetic effects and epistatic interaction effects, and we indicate opportunities to stabilize efficient breeding and selection of high-oil maize (Zea mays L.).

Genome-wide identification of GhAAI genes reveals that GhAAI66 triggers a phase transition to induce early flowering

Plants undergo a phase transition from vegetative to reproductive development that triggers floral induction. Genes containing an AAI (α-amylase inhibitor) domain form a large gene family, but there have been no comprehensive analyses of this gene family in any plant species. Here, we identified 336 AAI genes from nine plant species including122 AAI genes in cotton (Gossypium hirsutum). The AAI gene family has evolutionarily conserved amino acid residues throughout the plant kingdom. Phylogenetic analysis classified AAI genes into five major clades with significant polyploidization and showing effects of genome duplication. Our study identified 42 paralogous and 216 orthologous gene pairs resulting from segmental and whole-genome duplication, respectively, demonstrating significant contributions of gene duplication to expansion of the cotton AAI gene family. Further, GhAAI66 was preferentially expressed in flower tissue and as responses to phytohormone treatments. Ectopic expression of GhAAI66 in Arabidopsis and silencing in cotton revealed that GhAAI66 triggers a phase transition to induce early flowering. Further, GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis of RNA sequencing data and qRT-PCR (quantitative reverse transcription-PCR) analysis indicated that GhAAI66 integrates multiple flower signaling pathways including gibberellin, jasmonic acid, and floral integrators to trigger an early flowering cascade in Arabidopsis. Therefore, characterization of the AAI family provides invaluable insights for improving cotton breeding.

Updates on molecular mechanisms in the development of branched trichome in Arabidopsis and nonbranched in cotton

Trichomes are specialized epidermal cells and a vital plant organ that protect plants from various harms and provide valuable resources for plant development and use. Some key genes related to trichomes have been identified in the model plant Arabidopsis thaliana through glabrous mutants and gene cloning, and the hub MYB-bHLH-WD40, consisting of several factors including GLABRA1 (GL1), GL3, TRANSPARENT TESTA GLABRA1 (TTG1), and ENHANCER OF GLABRA3 (EGL3), has been established. Subsequently, some upstream transcription factors, phytohormones and epigenetic modification factors have also been studied in depth. In cotton, a very important fibre and oil crop globally, in addition to the key MYB-like factors, more important regulators and potential molecular mechanisms (e.g. epigenetic modifiers, distinct metabolic pathways) are being exploited during different fibre developmental stages. This occurs due to increased cotton research, resulting in the discovery of more complex regulation mechanisms from the allotetraploid genome of cotton. In addition, some conservative as well as specific mediators are involved in trichome development in other species. This study summarizes molecular mechanisms in trichome development and provides a detailed comparison of the similarities and differences between Arabidopsis and cotton, analyses the possible reasons for the discrepancy in identification of regulators, and raises future questions and foci for understanding trichome development in more detail.

Long non-coding RNAs and their potential functions in Ligon-lintless-1 mutant cotton during fiber development

Background

Long non-coding RNAs (LncRNAs) are part of genes, which are not translated into proteins and play a vital role in plant growth and development. Nevertheless, the presence of LncRNAs and how they functions in Ligon-lintless-1 mutant during the early cessation of cotton fiber development are still not well understood. In order to investigate the function of LncRNAs in cotton fiber development, it is necessary and important to identify LncRNAs and their potential roles in fiber cell development.

Results

In this work, we identified 18,333 LncRNAs, with the proportion of long intergenic noncoding RNAs (LincRNAs) (91.5%) and anti-sense LncRNAs (8.5%), all transcribed from Ligon-lintless-1 (Li1) and wild-type (WT). Expression differences were detected between Ligon-lintless-1 and wild-type at 0 and 8 DPA (day post anthesis). Pathway analysis and Gene Ontology based on differentially expressed LncRNAs on target genes, indicated fatty acid biosynthesis and fatty acid elongation being integral to lack of fiber in mutant cotton. The result of RNA-seq and RT-qPCR clearly singles out two potential LncRNAs, LNC_001237 and LNC_017085, to be highly down-regulated in the mutant cotton. The two LncRNAs were found to be destabilized or repressed by ghr-miR2950. Both RNA-seq analysis and RT-qPCR results in Ligon-lintless-1 mutant and wild-type may provide strong evidence of LNC_001237, LNC_017085 and ghr-miR2950 being integral molecular elements participating in various pathways of cotton fiber development.

Conclusion

The results of this study provide fundamental evidence for the better understanding of LncRNAs regulatory role in the molecular pathways governing cotton fiber development. Further research on designing and transforming LncRNAs will help not only in the understanding of their functions but will also in the improvement of fiber quality.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5978-5) contains supplementary material, which is available to authorized users.

Involvement of OpsLTP1 from Opuntia streptacantha in abiotic stress adaptation and lipid metabolism

Plant lipid transfer proteins (LTPs) exhibit the ability to transfer lipids between membranes in vitro, and have been implicated in diverse physiological processes associated to plant growth, reproduction, development, biotic and abiotic stress responses. However, their mode of action is not yet fully understood. To explore the functions of the OpsLTP1 gene encoding a LTP from cactus pear Opuntia streptacantha Lem., we generated transgenic Arabidopsis thaliana (L.) Heynh. plants to overexpress OpsLTP1 and contrasted our results with the loss-of-function mutant ltp3 from A. thaliana under abiotic stress conditions. The ltp3 mutant seeds showed impaired germination under salt and osmotic treatments, in contrast to OpsLTP1 overexpressing lines that displayed significant increases in germination rate. Moreover, stress recovery assays showed that ltp3 mutant seedlings were more sensitive to salt and osmotic treatments than wild-type plants suggesting that AtLTP3 is required for stress-induced responses, while the OpsLTP1 overexpressing line showed no significant differences. In addition, OpsLTP1 overexpressing and ltp3 mutant seeds stored lower amount of total lipids compared with wild-type seeds, showing changes primarily on 16C and 18C fatty acids. However, ltp3 mutant also lead changes in lipid profile and no over concrete lipids which may suggest a compensatory activation of other LTPs. Interestingly, linoleic acid (18:2ω6) was consistently increased in neutral, galactoglycerolipids and phosphoglycerolipids of OpsLTP1 overexpressing line indicating a role of OpsLTP1 in the modulation of lipid composition in A. thaliana.

Systematic Analysis of Cotton Non-specific Lipid Transfer Protein Family Revealed a Special Group That Is Involved in Fiber Elongation

Non-specific lipid transfer proteins (nsLTPs) had been previously isolated from cotton fiber but their functions were unclear so far. Bioinformatic analysis of the tetraploid cotton genome database identified 138 nsLTP genes, falling into the 11 groups as reported previously. Different from Arabidopsis, cacao, and other crops, cotton type XI genes were considerably expanded and diverged earlier on chromosome At11, Dt11, and Dt08. Corresponding to the type XI genes, the type XI proteins (GhLtpXIs) all contained an extra N-terminal cap resulting in larger molecular weight. The research revealed that the expression of type XI genes was dramatically increased in fibers of tetraploid cotton compared with the two diploid progenitors. High-level of GhLtpXIs expression was observed in long-fibered cotton cultivars during fiber elongation. Ectopic expression of GhLtpXIs in Arabidopsis significantly enhanced trichome length, suggesting that GhLtpXIs promoted fiber elongation. Overall, the findings of this research provide insights into phenotypic evolution of Gossypium species and regulatory mechanism of nsLTPs during fiber development. HIGHLIGHT A specific group, type XI nsLTPs, was identified with predominant expression in elongating fibers of Gossypium hirsutum based on evolutionary, transcriptional, and functional analyses.

Genome-wide association study discovered genetic variation and candidate genes of fibre quality traits in Gossypium hirsutum L

Genetic improvement of fibre quality is one of the main breeding goals for the upland cotton, Gossypium hirsutum, but there are difficulties with precise selection of traits. Therefore, it is important to improve the understanding of the genetic basis of phenotypic variation. In this study, we conducted phenotyping and genetic variation analyses of 719 diverse accessions of upland cotton based on multiple environment tests and a recently developed Cotton 63K Illumina Infinium SNP array and performed a genome-wide association study (GWAS) of fibre quality traits. A total of 10 511 polymorphic SNPs distributed in 26 chromosomes were screened across the cotton germplasms, and forty-six significant SNPs associated with five fibre quality traits were detected. These significant SNPs were scattered over 15 chromosomes and were involved in 612 unique candidate genes, many related to polysaccharide biosynthesis, signal transduction and protein translocation. Two major haplotypes for fibre length and strength were identified on chromosomes Dt11 and At07. Furthermore, by combining GWAS and transcriptome analysis, we identified 163 and 120 fibre developmental genes related to length and strength, respectively, of which a number of novel genes and 19 promising genes were screened. These results provide new insight into the genetic basis of fibre quality in G. hirsutum and provide candidate SNPs and genes to accelerate the improvement of upland cotton.

Proteomic Analysis of Differences in Fiber Development between Wild and Cultivated Gossypium hirsutum L

Upland cotton (Gossypium hirsutum L.) is one of the world's most important fiber crops, accounting for more than 90% of all cotton production. While their wild progenitors have relatively short and coarse, often tan-colored fibers, modern cotton cultivars possess longer, finer, stronger, and whiter fiber. In this study, the wild and cultivated cottons (YU-3 and TM-1) selected show significant differences on fibers at 10 days postanthesis (DPA), 20 DPA, and mature stages at the morphological level. To explore the effects of domestication, reveal molecular mechanisms underlying these phenotypic differences, and better inform our efforts to further enhance cotton fiber quality, isobaric tags for relative and absolute protein quantification-facilitated proteomic methods were performed on developing fibers. There were 6990 proteins identified; among them, 336 were defined as differentially expressed proteins between fibers of wild versus domesticated cotton. The down- or up-regulated proteins in wild cotton were involved in phenylpropanoid biosynthesis, zeatin biosynthesis, fatty acid elongation, and other processes. Association analysis between transcriptome and proteome showed positive correlations between transcripts and proteins at both 10 DPA and 20 DPA. Differences in proteomics have been verified at the mRNA level by quantitative real-time polymerase chain reaction and have been validated at the physiological and biochemical levels by POD (peroxidase) activity assays and ZA (zeatin) content estimates. This work corroborates the major pathways involved in cotton fiber development and demonstrates that POD activity and zeatin content have a great potential related to fiber elongation and thickening.

Glycosylphosphatidylinositol (GPI) Modification Serves as a Primary Plasmodesmal Sorting Signal

Plasmodesmata (Pd) are membranous channels that serve as a major conduit for cell-to-cell communication in plants. The Pd-associated β-1,3-glucanase (BG_pap) and CALLOSE BINDING PROTEIN1 (PDCB1) were identified as key regulators of Pd conductivity. Both are predicted glycosylphosphatidylinositol-anchored proteins (GPI-APs) carrying a conserved GPI modification signal. However, the subcellular targeting mechanism of these proteins is unknown, particularly in the context of other GPI-APs not associated with Pd Here, we conducted a comparative analysis of the subcellular targeting of the two Pd-resident and two unrelated non-Pd GPI-APs in Arabidopsis (Arabidopsis thaliana). We show that GPI modification is necessary and sufficient for delivering both BG_pap and PDCB1 to Pd Moreover, the GPI modification signal from both Pd- and non-Pd GPI-APs is able to target a reporter protein to Pd, likely to plasma membrane microdomains enriched at Pd As such, the GPI modification serves as a primary Pd sorting signal in plant cells. Interestingly, the ectodomain, a region that carries the functional domain in GPI-APs, in Pd-resident proteins further enhances Pd accumulation. However, in non-Pd GPI-APs, the ectodomain overrides the Pd targeting function of the GPI signal and determines a specific GPI-dependent non-Pd localization of these proteins at the plasma membrane and cell wall. Domain-swap analysis showed that the non-Pd localization is also dominant over the Pd-enhancing function mediated by a Pd ectodomain. In conclusion, our results indicate that segregation between Pd- and non-Pd GPI-APs occurs prior to Pd targeting, providing, to our knowledge, the first evidence of the mechanism of GPI-AP sorting in plants.

GhMCS1, the Cotton Orthologue of Human GRIM-19, Is a Subunit of Mitochondrial Complex I and Associated with Cotton Fibre Growth

GRIM-19 (Gene associated with Retinoid-Interferon-induced Mortality 19) is a subunit of mitochondrial respiratory complex I in mammalian systems, and it has been demonstrated to be a multifunctional protein involved in the cell cycle, cell motility and innate immunity. However, little is known about the molecular functions of its homologues in plants. Here, we characterised GhMCS1, an orthologue of human GRIM-19 from cotton (Gossypium hirsutum L.), and found that it was essential for maintaining complex integrity and mitochondrial function in cotton. GhMCS1 was detected in various cotton tissues, with high levels expressed in developing fibres and flowers and lower levels in leaves, roots and ovules. In fibres at different developmental stages, GhMCS1 expression peaked at 5-15 days post anthesis (dpa) and then decreased at 20 dpa and diminished at 25 dpa. By Western blot analysis, GhMCS1 was observed to be localised to the mitochondria of cotton leaves and to colocalise with complex I. In Arabidopsis, GhMCS1 overexpression enhanced the assembly of complex I and thus respiratory activity, whereas the GhMCS1 homologue (At1g04630) knockdown mutants showed significantly decreased respiratory activities. Furthermore, the mutants presented with some phenotypic changes, such as smaller whole-plant architecture, poorly developed seeds and fewer trichomes. More importantly, in the cotton fibres, both the GhMCS1 transcript and protein levels were correlated with respiratory activity and fibre developmental phase. Our results suggest that GhMCS1, a functional ortholog of the human GRIM-19, is an essential subunit of mitochondrial complex I and is involved in cotton fibre development. The present data may deepen our knowledge on the potential roles of mitochondria in fibre morphogenesis.