The Ligon lintless -2 Short Fiber Mutation Is Located within a Terminal Deletion of Chromosome 18 in Cotton

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Systematically and Comprehensively Understanding the Regulation of Cotton Fiber Initiation: A Review

Cotton fibers provide an important source of raw materials for the textile industry worldwide. Cotton fiber is a kind of single cell that differentiates from the epidermis of the ovule and provides a perfect research model for the differentiation and elongation of plant cells. Cotton fiber initiation is the first stage throughout the entire developmental process. The number of fiber cell initials on the seed ovule epidermis decides the final fiber yield. Thus, it is of great significance to clarify the mechanism underlying cotton fiber initiation. Fiber cell initiation is controlled by complex and interrelated regulatory networks. Plant phytohormones, transcription factors, sugar signals, small signal molecules, functional genes, non-coding RNAs, and histone modification play important roles during this process. Here, we not only summarize the different kinds of factors involved in fiber cell initiation but also discuss the mechanisms of these factors that act together to regulate cotton fiber initiation. Our aim is to synthesize a systematic and comprehensive review of different factors during fiber initiation that will provide the basics for further illustrating these mechanisms and offer theoretical guidance for improving fiber yield in future molecular breeding work.

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development

Introduction

Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels.

Methods

Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes.

Results

Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers.

Discussion

Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

N6 -Methyladenosine mRNA modification regulates transcripts stability associated with cotton fiber elongation

N6 -Methyladenosine (m6 A) is the most abundant methylation modification in eukaryotic mRNA. The discovery of the dynamic and reversible regulatory mechanism of m6 A has greatly promoted the development of m6 A-led epitranscriptomics. However, the characterization of m6 A in cotton fiber is still unknown. Here, we reveal the potential link between m6 A modification and cotton fiber elongation by parallel m6 A-immunoprecipitation-sequencing (m6 A-seq) and RNA-seq analysis of fibers from the short fiber mutants Ligonliness-2 (Li2 ) and wild-type (WT). This study demonstrated a higher level of m6 A in the Li2 mutant, with the enrichment of m6 A modifications in the stop codon, 3'-untranslated region and coding sequence regions than in WT cotton. In the correlation analysis between genes containing differential m6 A modifications and differentially expressed genes, we identified several genes that could potentially regulate fiber elongation, including cytoskeleton, microtubule binding, cell wall and transcription factors (TFs). We further confirmed that the methylation of m6 A affected the mRNA stability of these fiber elongation-related genes including the TF GhMYB44, which showed the highest expression level in the RNA-seq data and m6 A methylation in the m6 A-seq data. Next, the overexpression of GhMYB44 reduces fiber elongation, whereas the silencing of GhMYB44 produces longer fibers. In summary, these results uncover that m6 A methylation regulated the expression of genes related to fiber development by affecting mRNA's stability, ultimately affecting cotton fiber elongation.

Characterization of chromatin accessibility and gene expression reveal the key genes involved in cotton fiber elongation

Cotton (Gossypium hirsutum L.) is an important economic crop, and cotton fiber is one of the longest plant cells, which provides an ideal model for the study of cell elongation and secondary cell wall synthesis. Cotton fiber length is regulated by a variety of transcription factors (TF) and their target genes; however, the mechanism of fiber elongation mediated by transcriptional regulatory networks is still unclear to a large extent. Here, we used a comparative assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) assay and RNA-seq analysis to identify fiber elongation transcription factors and genes using the short-fiber mutant ligon linless-2 (Li2 ) and wild type (WT). A total of 499 differential target genes were identified and GO analysis shows that differential genes are mainly involved in plant secondary wall synthesis and microtubule-binding processes. Analysis of the genomic regions preferentially accessible (Peak) has identified a number of overrepresented TF-binding motifs, highlighting sets of TFs that are important for cotton fiber development. Using ATAC-seq and RNA-seq data, we have constructed a functional regulatory network of each TF regulatory target gene and also the network pattern of TF regulating differential target genes. Further, to obtain the genes related to fiber length, the differential target genes were combined with FLGWAS data to identify the genes highly related to fiber length. Our work provides new insights into cotton fiber elongation.

Phenomics and transcriptomics analyses reveal deposition of suberin and lignin in the short fiber cell walls produced from a wild cotton species and two mutants

Fiber length is one of the major properties determining the quality and commercial value of cotton. To understand the mechanisms regulating fiber length, genetic variations of cotton species and mutants producing short fibers have been compared with cultivated cottons generating long and normal fibers. However, their phenomic variation other than fiber length has not been well characterized. Therefore, we compared physical and chemical properties of the short fibers with the long fibers. Fiber characteristics were compared in two sets: 1) wild diploid Gossypium raimondii Ulbrich (short fibers) with cultivated diploid G. arboreum L and tetraploid G. hirsutum L. (long fibers); 2) G. hirsutum short fiber mutants, Ligon-lintless 1 (Li1) and 2 (Li2) with their near isogenic line (NIL), DP-5690 (long fibers). Chemical analyses showed that the short fibers commonly consisted of greater non-cellulosic components, including lignin and suberin, than the long fibers. Transcriptomic analyses also identified up-regulation of the genes related to suberin and lignin biosynthesis in the short fibers. Our results may provide insight on how high levels of suberin and lignin in cell walls can affect cotton fiber length. The approaches combining phenomic and transcriptomic analyses of multiple sets of cotton fibers sharing a common phenotype would facilitate identifying genes and common pathways that significantly influence cotton fiber properties.

A deletion/duplication in the Ligon lintless-2 locus induces siRNAs that inhibit cotton fiber cell elongation

Most cultivated cotton (Gossypium hirsutum L.) varieties have two types of seed fibers: short fuzz fiber strongly adhered to the seed coat, and long lint fiber used in the textile industry. The Ligon lintless-2 (Li2) cotton mutant has a normal vegetative phenotype but produces very short lint fiber on the seeds. The Li2 mutation is controlled by a single dominant gene. We discovered a large structural rearrangement at the end of chromosome D13 in the Li2 mutant based on whole-genome sequencing and genetic mapping of segregating populations. The rearrangement contains a 177-kb deletion and a 221-kb duplication positioned as a tandem inverted repeat. The gene Gh_D13G2437 is located at the junction of the inverted repeat in the duplicated region. During transcription such structure spontaneously forms self-complementary hairpin RNA of Gh_D13G2437 followed by production of small interfering RNA (siRNA). Gh_D13G2437 encodes a Ran-Binding Protein 1 (RanBP1) that preferentially expresses during cotton fiber elongation. The abundance of siRNA produced from Gh_D13G2437 reciprocally corresponds with the abundance of highly homologous (68%-98% amino acid sequence identity) RanBP1 family transcripts during fiber elongation, resulting in a shorter fiber phenotype in the Li2. Overexpression of Gh_D13G2437 in the Li2 mutant recovered the long lint fiber phenotype. Taken together, our findings revealed that siRNA-induced silencing of a family of RanBP1s inhibit elongation of cotton fiber cells in the Li2 mutant.

Analysis of transcriptome data and quantitative trait loci enables the identification of candidate genes responsible for fiber strength in Gossypium barbadense

Gossypium barbadense possesses a superior fiber quality because of its fiber length and strength. An in-depth analysis of the underlying genetic mechanism could aid in filling the gap in research regarding fiber strength and could provide helpful information for Gossypium barbadense breeding. Three quantitative trait loci related to fiber strength were identified from a Gossypium barbadense recombinant inbred line (PimaS-7 × 5917) for further analysis. RNA sequencing was performed in the fiber tissues of PimaS-7 × 5917 0–35 days postanthesis. Four specific modules closely related to the secondary wall-thickening stage were obtained using the weighted gene coexpression network analysis. In total, 55 genes were identified as differentially expressed from 4 specific modules. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes were used for enrichment analysis, and Gbar_D11G032910, Gbar_D08G020540, Gbar_D08G013370, Gbar_D11G033670, and Gbar_D11G029020 were found to regulate fiber strength by playing a role in the composition of structural constituents of cytoskeleton and microtubules during fiber development. Quantitative real-time PCR results confirmed the accuracy of the transcriptome data. This study provides a quick strategy for exploring candidate genes and provides new insights for improving fiber strength in cotton.

Dynamic roles and intricate mechanisms of ethylene in epidermal hair development in Arabidopsis and cotton

Ethylene affects many aspects of plant growth and development, including root hairs and trichomes growth in Arabidopsis, as well as fiber development in cotton, though the underlying mechanism is unclear. In this article, we update the research progress associated with the main genes in ethylene biosynthesis and signaling pathway, and we propose a clear ethylene pathway based on genome-wide identification of homologues in cotton. Expression pattern analysis using transcriptome data revealed that some candidate genes may contribute to cotton fiber development through the ethylene pathway. Moreover, we systematically summarized the effects of ethylene on the development of epidermal hair and the underlying regulatory mechanisms in Arabidopsis. Based on the knowledge of ethylene-promoted cell differentiation, elongation, and development in different tissues or plants, we advised a possible regulatory network for cotton fiber development with ethylene as the hub. Importantly, we emphasized the roles of ethylene as an important node in regulating cotton vegetative growth, and stress resistance, and suggested utilizing multiple methods to subtly modify ethylene synthesis or signaling in a tissue or spatiotemporal-specific manner to clarify its exact effect on architecture, adaptability of the plant, and fiber development, paving the way for basic research and genetic improvement of the cotton crop.

Genome-wide identification and expression analysis of GL2-interacting-repressor (GIR) genes during cotton fiber and fuzz development

GL2-interacting-repressor (GIR) family members may contribute to fiber/fuzz formation via a newly discovered unique pathway in Gossypium arboreum. There are similarities between cotton fiber development and the formation of trichomes and root hairs. The GL2-interacting-repressors (GIRs) are crucial regulators of root hair and trichome formation. The GaFzl gene, annotated as GaGIR1, is negatively associated with trichome development and fuzz initiation. However, there is relatively little available information regarding the other GIR genes in cotton, especially regarding their effects on cotton fiber development. In this study, 21 GIR family genes were identified in the diploid cotton species Gossypium arboreum; these genes were divided into three groups. The GIR genes were characterized in terms of their phylogenetic relationships, structures, chromosomal distribution and evolutionary dynamics. These GIR genes were revealed to be unequally distributed on 12 chromosomes in the diploid cotton genome, with no GIR gene detected on Ga06. The cis-acting elements in the promoter regions were predicted to be responsive to light, phytohormones, defense activities and stress. The transcriptomic data and qRT-PCR results revealed that most GIR genes were not differentially expressed between the wild-type control and the fuzzless mutant line. Moreover, 14 of 21 family genes were expressed at high levels, indicating these genes may play important roles during fiber development and fuzz formation. Furthermore, Ga01G0231 was predominantly expressed in root samples, suggestive of a role in root hair formation rather than in fuzz initiation and development. The results of this study have enhanced our understanding of the GIR genes and their potential utility for improving cotton fiber through breeding.

Quantitative Trait Loci and Transcriptome Analysis Reveal Genetic Basis of Fiber Quality Traits in CCRI70 RIL Population of Gossypium hirsutum

Upland cotton (Gossypium hirsutum) is widely planted around the world for its natural fiber, and producing high-quality fiber is essential for the textile industry. CCRI70 is a hybrid cotton plant harboring superior yield and fiber quality, whose recombinant inbred line (RIL) population was developed from two upland cotton varieties (sGK156 and 901-001) and were used here to investigate the source of high-quality related alleles. Based on the material of the whole population, a high-density genetic map was constructed using specific locus-amplified fragment sequencing (SLAF-seq). It contained 24,425 single nucleotide polymorphism (SNP) markers, spanning a distance of 4,850.47 centimorgans (cM) over 26 chromosomes with an average marker interval of 0.20 cM. In evaluating three fiber quality traits in nine environments to detect multiple environments stable quantitative trait loci (QTLs), we found 289 QTLs, of which 36 of them were stable QTLs and 18 were novel. Based on the transcriptome analysis for two parents and two RILs, 24,941 unique differentially expressed genes (DEGs) were identified, 473 of which were promising genes. For the fiber strength (FS) QTLs, 320 DEGs were identified, suggesting that pectin synthesis, phenylpropanoid biosynthesis, and plant hormone signaling pathways could influence FS, and several transcription factors may regulate fiber development, such as GAE6, C4H, OMT1, AFR18, EIN3, bZIP44, and GAI. Notably, the marker D13_56413025 in qFS-chr18-4 provides a potential basis for enhancing fiber quality of upland cotton via marker-assisted breeding and gene cloning of important fiber quality traits.

Co-expression network and comparative transcriptome analysis for fiber initiation and elongation reveal genetic differences in two lines from upland cotton CCRI70 RIL population

Upland cotton is the most widely planted for natural fiber around the world, and either lint percentage (LP) or fiber length (FL) is the crucial component tremendously affecting cotton yield and fiber quality, respectively. In this study, two lines MBZ70-053 and MBZ70-236 derived from G. hirsutum CCRI70 recombinant inbred line (RIL) population presenting different phenotypes in LP and FL traits were chosen to conduct RNA sequencing on ovule and fiber samples, aiming at exploring the differences of molecular and genetic mechanisms during cotton fiber initiation and elongation stages. As a result, 249/128, 369/206, 4296/1198 and 3547/2129 up-/down- regulated differentially expressed genes (DGEs) in L2 were obtained at -3, 0, 5 and 10 days post-anthesis (DPA), respectively. Seven gene expression profiles were discriminated using Short Time-series Expression Miner (STEM) analysis; seven modules and hub genes were identified using weighted gene co-expression network analysis. The DEGs were mainly enriched into energetic metabolism and accumulating as well as auxin signaling pathway in initiation and elongation stages, respectively. Meanwhile, 29 hub genes were identified as 14-3-3ω , TBL35, GhACS, PME3, GAMMA-TIP, PUM-7, etc., where the DEGs and hub genes revealed the genetic and molecular mechanisms and differences during cotton fiber development.

Wheat heat tolerance is impaired by heightened deletions in the distal end of 4AL chromosomal arm

Heat stress (HS) causes substantial damages to worldwide crop production. As a cool season crop, wheat (Triticum aestivum) is sensitive to HS-induced damages. To support the genetic improvement of wheat HS tolerance (HST), we conducted fine mapping of TaHST1, a locus required for maintaining wheat vegetative and reproductive growth under elevated temperatures. TaHST1 was mapped to the distal terminus of 4AL chromosome arm using genetic populations derived from two BC₆ F₆ breeding lines showing tolerance (E6015-4T) or sensitivity (E6015-3S) to HS. The 4AL region carrying TaHST1 locus was approximately 0.949 Mbp and contained the last 19 high confidence genes of 4AL according to wheat reference genome sequence. Resequencing of E6015-3S and E6015-4T and haplotype analysis of 3087 worldwide wheat accessions revealed heightened deletion polymorphisms in the distal 0.949 Mbp region of 4AL, which was confirmed by the finding of frequent gene losses in this region in eight genome-sequenced hexaploid wheat cultivars. The great majority (86.36%) of the 3087 lines displayed different degrees of nucleotide sequence deletions, with only 13.64% of them resembling E6015-4T in this region. These deletions can impair the presence and/or function of TaHST1 and surrounding genes, thus rendering global wheat germplasm vulnerable to HS or other environmental adversities. Therefore, conscientious and urgent efforts are needed in global wheat breeding programmes to optimize the structure and function of 4AL distal terminus by ensuring the presence of TaHST1 and surrounding genes. The new information reported here will help to accelerate the ongoing global efforts in improving wheat HST.

Surface and Thermal Characterization of Cotton Fibers of Phenotypes Differing in Fiber Length

Cotton is one of the most important and widely grown crops in the world. Understanding the synthesis mechanism of cotton fiber elongation can provide valuable tools to the cotton industry for improving cotton fiber yield and quality at the molecular level. In this work, the surface and thermal characteristics of cotton fiber samples collected from a wild type (WT) and three mutant lines (Li1, Li2-short, Li2-long, Li2-mix, and liy) were comparatively investigated. Microimaging revealed a general similarity trend of WT ≥ Li2-long ≈ Li2-mix > Li1 > Li2 short ≈ liy with Ca detected on the surface of the last two. Attenuated total reflectance Fourier transform infrared (ATR FT-IR) spectroscopy and thermogravimetric measurements also showed that Li2-short and liy were more similar to each other, and Li2-long and Li2-mix closer to WT while Li1 was quite independent. FT-IR results further demonstrated that wax and amorphous cellulose were co-present in fiber structures during the fiber formation processes. The correlation analysis found that the FT-IR-based maturity parameter was well correlated (p ≤ 0.05) to the onset decomposition temperature and all three weight-loss parameters at onset, peak, and end decomposition stages, suggesting that the maturity degree is a better parameter than crystallinity index (CI) and other FT-IR parameters that reflect the thermal stability of the cotton fiber. In summary, this work demonstrated that genetic mutation altered the surface and thermal characteristics in the same way for Li2-short and liy, but with different mechanisms for the other three mutant cotton fiber samples.

Genetic mapping and transcriptomic characterization of a new fuzzless-tufted cottonseed mutant

Fiber mutants are unique and valuable resources for understanding the genetic and molecular mechanisms controlling initiation and development of cotton fibers that are extremely elongated single epidermal cells protruding from the seed coat of cottonseeds. In this study, we reported a new fuzzless-tufted cotton mutant (Gossypium hirsutum) and showed that fuzzless-tufted near-isogenic lines (NILs) had similar agronomic traits and a higher ginning efficiency compared to their recurrent parents with normal fuzzy seeds. Genetic analysis revealed that the mutant phenotype is determined by a single incomplete dominant locus, designated N5. The mutation was fine mapped to an approximately 250-kb interval containing 33 annotated genes using a combination of bulked segregant sequencing, SNP chip genotyping, and fine mapping. Comparative transcriptomic analysis using 0-6 days post-anthesis (dpa) ovules from NILs segregating for the phenotypes of fuzzless-tufted (mutant) and normal fuzzy cottonseeds (wild-type) uncovered candidate genes responsible for the mutant phenotype. It also revealed that the flanking region of the N5 locus is enriched with differentially expressed genes (DEGs) between the mutant and wild-type. Several of those DEGs are members of the gene families with demonstrated roles in cell initiation and elongation, such as calcium-dependent protein kinase and expansin. The transcriptome landscape of the mutant was significantly reprogrammed in the 6 dpa ovules and, to a less extent, in the 0 dpa ovules, but not in the 2 and 4 dpa ovules. At both 0 and 6 dpa, the reprogrammed mutant transcriptome was mainly associated with cell wall modifications and transmembrane transportation, while transcription factor activity was significantly altered in the 6 dpa mutant ovules. These results imply a similar molecular basis for initiation of lint and fuzz fibers despite certain differences.