Transcript profiling by microarray and marker analysis of the short cotton (Gossypium hirsutum L.) fiber mutant Ligon lintless-1 (Li1)

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.

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.

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.

Transcriptional Landscape of Cotton Fiber Development and Its Alliance With Fiber-Associated Traits

Cotton fiber development is still an intriguing question to understand fiber commitment and development. At different fiber developmental stages, many genes change their expression pattern and have a pivotal role in fiber quality and yield. Recently, numerous studies have been conducted for transcriptional regulation of fiber, and raw data were deposited to the public repository for comprehensive integrative analysis. Here, we remapped > 380 cotton RNAseq data with uniform mapping strategies that span ∼400 fold coverage to the genome. We identified stage-specific features related to fiber cell commitment, initiation, elongation, and Secondary Cell Wall (SCW) synthesis and their putative cis-regulatory elements for the specific regulation in fiber development. We also mined Exclusively Expressed Transcripts (EETs) that were positively selected during cotton fiber evolution and domestication. Furthermore, the expression of EETs was validated in 100 cotton genotypes through the nCounter assay and correlated with different fiber-related traits. Thus, our data mining study reveals several important features related to cotton fiber development and improvement, which were consolidated in the "CottonExpress-omics" database.

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.

A Modified Actin (Gly65Val Substitution) Expressed in Cotton Disrupts Polymerization of Actin Filaments Leading to the Phenotype of Ligon Lintless-1 (Li1) Mutant

Dynamic remodeling of the actin cytoskeleton plays a central role in the elongation of cotton fibers, which are the most important natural fibers in the global textile industry. Here, a high-resolution mapping approach combined with comparative sequencing and a transgenic method revealed that a G65V substitution in the cotton actin Gh_D04G0865 (GhACT17D in the wild-type) is responsible for the Gossypium hirsutum Ligon lintless-1 (Li₁) mutant (GhACT17DM). In the mutant GhACT17DM from Li₁ plant, Gly65 is substituted with valine on the lip of the nucleotide-binding domain of GhACT17D, which probably affects the polymerization of F-actin. Over-expression of GhACT17DM, but not GhACT17D, driven by either a CaMV35 promoter or a fiber-specific promoter in cotton produced a Li₁-like phenotype. Compared with the wild-type control, actin filaments in Li₁ fibers showed higher growth and shrinkage rates, decreased filament skewness and parallelness, and increased filament density. Taken together, our results indicate that the incorporation of GhACT17DM during actin polymerization disrupts the establishment and dynamics of the actin cytoskeleton, resulting in defective fiber elongation and the overall dwarf and twisted phenotype of the Li₁ mutant.

GaHD1, a candidate gene for the Gossypium arboreum SMA-4 mutant, promotes trichome and fiber initiation by cellular H2O2 and Ca2+ signals

Cotton fibers are initiated from the epidermal cells of the ovule before or on the day of anthesis. Gossypium arboreum SMA-4 mutant contains recessive mutation (sma-4(ha)) and has the phenotypes of fibreless seeds and glabrous stems. In this study, fine mapping and alternative splicing analysis indicated a nucleotide substitution (AG → AC) at splicing site in a homeodomain-leucine zipper IV family gene (GaHD1) might cause gene A3S (Alternative 3' splicing) mistake, suggested that GaHD1 was the candidate gene of sma-4(ha). Many genes related to the fiber initiation are identified to be differentially expressed in the mutant which could result in the blocked fiber initiation signals such as H₂O_2, or Ca in the mutant. Further comparative physiological analysis of H₂O₂ production and Ca²⁺ flux in the SMA-4 and wide type cotton confirmed that H₂O₂ and Ca were important fiber initiation signals and regulated by GaHD1. The in vitro ovule culture of the mutant with hormones recovered the fibered phenotype coupled with the restoration of these signals. Overexpressing of GaHD1 in Arabidopsis increased trichome densities on the sepal, leaf, and stem tissues while transient silencing of the GaHD1 gene in G. arboreum reduced the trichome densities. These phenotypes indicated that GaHD1 is the candidate gene of SMA-4 with a crucial role in acting upstream molecular switch of signal transductions for cotton trichome and fiber initiations.

An EMS-induced mutation in a tetratricopeptide repeat-like superfamily protein gene (Ghir_A12G008870) on chromosome A12 is responsible for the liy short fiber phenotype in cotton

The EMS-induced threonine/isoleucine substitution in a tetratricopeptide repeat-like superfamily protein encoded by gene Ghir_A12G008870 is responsible for the Ligon-lintless-y (li_y) short fiber phenotype in cotton. A short fiber mutant Ligon-lintless-y was created through treating the seeds of the cotton line MD15 with ethyl methanesulfonate. Genetic analysis indicated that the short fiber phenotype is controlled by a single recessive locus designated li_y. From F₂ populations derived from crosses between the mutant and its wild type (WT), we selected 132 short fiber progeny (li_y/li_y) and made two DNA bulks. We sequenced these DNA bulks along with the two parents of the population. The li_y locus was located on chromosome A12. Using multiple F₂ populations and F₃ progeny plants, we mapped the li_y locus within a genomic region of 1.18 Mb. In this region, there is only one gene, i.e., Ghir_A12G008870 encoding a tetratricopeptide repeat-like superfamily protein that has a non-synonymous mutation between the li_y mutant and its WT. Analysis of a SNP marker representing this gene in the F₂ and F₃ progeny plants demonstrated its complete linkage with the li_y short fiber phenotype. We further analyzed this SNP marker in a panel of 384 cotton varieties. The mutant allele is absent in all varieties analyzed. RNAseq and RT-qPCR analysis of the gene Ghir_A12G008870 during fiber development showed a significant expression difference between the li_y mutant and its WT in developing fiber cells beginning at 12 days post-anthesis. Virus-induced gene silencing of the gene Ghir_A12G008870 significantly reduced the fiber length of the WT cotton line MD15. Taken together, our results suggest that the gene Ghir_A12G008870 is involved in the cotton fiber cell elongation process and is a promising candidate gene responsible for the li_y short fiber phenotype.

G65V Substitution in Actin Disturbs Polymerization Leading to Inhibited Cell Elongation in Cotton

The importance of the actin cytoskeleton for proper cell development has been well established in a variety of organisms. Actin protein sequences are highly conserved, and each amino acid residue may be essential for its function. In this study, we report the isolation and characterization of GhLi ₁ from an upland cotton mutant Ligon lintless-1 (Li₁), which harbors the G65V substitution in its encoded actin protein. Li₁ mutants exhibit pleiotropic malformed phenotypes, including dwarf plants, distorted organs, and extremely shortened fibers. Cytological analysis showed that the actin cytoskeleton was disorganized and the abundance of F-actin was decreased in the Li₁ cells. Vesicles were aggregated into patches, and excessive cellulose synthase complexes were inserted into the plasma membrane during the secondary cell wall biosynthesis stage, which dramatically affected the morphology of the Li₁ cells. Molecular model prediction suggested that the G65V substitution may affect the three-bodied G-actin interaction during F-actin assembly. Biochemical assays demonstrated that the recombinant GhLi₁ protein disturbs actin dynamics by inhibiting the nucleation and elongation processes. Therefore, our findings demonstrate that the G65V substitution in actin had dominant-negative effects on cell elongation, by disturbing actin polymerization and actin cytoskeleton-based biological processes such as intracellular transportation.

Genome-wide analysis of cotton C2H2-zinc finger transcription factor family and their expression analysis during fiber development

BACKGROUND

C2H2-zinc finger protein family is commonly found in the plant, and it is known as the key actors in the regulation of transcription and vital component of chromatin structure. A large number of the C2H2-zinc finger gene members have not been well characterized based on their functions and structure in cotton. However, in other plants, only a few C2H2-zinc finger genes have been studied.

RESULTS

In this work, we performed a comprehensive analysis and identified 386, 196 and 195 C2H2-zinc finger genes in Gossypium hirsutum (upland cotton), Gossypium arboreum and Gossypium raimondii, respectively. Phylogenetic tree analysis of the C2H2-zinc finger proteins encoding the C2H2-zinc finger genes were classified into seven (7) subgroups. Moreover, the C2H2-zinc finger gene members were distributed in all cotton chromosomes though with asymmetrical distribution patterns. All the orthologous genes were detected between tetraploid and the diploid cotton, with 154 orthologous genes pair detected between upland cotton and Gossypium arboreum while 165 orthologous genes were found between upland cotton and Gossypium raimondii. Synonymous (Ks) and non-synonymous (Ka) nucleotide substitution rates (Ka/Ks) analysis indicated that the cotton C2H2-zinc finger genes were highly influenced mainly by negative selection, which maintained their protein levels after the duplication events. RNA-seq data and RT-qPCR validation of the RNA seq result revealed differential expression pattern of some the C2H2-zinc finger genes at different stages of cotton fiber development, an indication that the C2H2-zinc finger genes play an important role in initiating and regulating fiber development in cotton.

CONCLUSIONS

This study provides a strong foundation for future practical genome research on C2H2-zinc finger genes in upland cotton. The expression levels of C2H2-zinc finger genes family is a pointer of their involvement in various biochemical and physiological functions which are directly related to cotton fiber development during initiation and elongation stages. This work not only provides a basis for determining the nominal role of the C2H2-zinc finger genes in fiber development but also provide valuable information for characterization of potential candidate genes involved in regulation of cotton fiber development.

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.

Protecting cotton crops under elevated CO2 from waterlogging by managing ethylene

Soil waterlogging and subsequent ethylene release from cotton (Gossypium hirsutum L.) tissues has been linked with abscission of developing cotton fruits. This glasshouse study investigates the effect of a 9-day waterlogging event and CO2 enrichment (eCO2, 700 parts per million (ppm)) on a fully linted cultivar 'Empire' and a lintless cotton mutant (5B). We hypothesised that cotton performance in extreme environments such as waterlogging can be improved through mitigating ethylene action. Plants were grown at 28:20°C day:night temperature, 50-70% relative humidity and a 14:10 light:dark photoperiod under natural light and were exposed to waterlogging and eCO2 at early reproductive growth. Ethylene synthesis was inhibited by spraying aminoethoxyvinylglycine (830ppm) 1 day before waterlogging. Waterlogging significantly increased ethylene release from both cotton genotypes, although fruit production was significantly inhibited only in Empire. Aminoethoxyvinylglycine consistently reduced waterlogging-induced abscission of fruits, mainly in Empire. Limited damage to fruits in 5B, despite increased ethylene production during waterlogging, suggested that fruit abscission in 5B was inhibited by disrupting ethylene metabolism genetically. Elevated CO2 promoted fruit production in both genotypes and was more effective in 5B than in Empire plants. Hence 5B produced more fruits than Empire, providing additional sinks (existing and new fruit) that enhanced the response to CO2 enrichment.

Endogenous Ethylene Concentration Is Not a Major Determinant of Fruit Abscission in Heat-Stressed Cotton (Gossypium hirsutum L.)

We investigated the role of ethylene in the response of cotton to high temperature using cotton genotypes with genetically interrupted ethylene metabolism. In the first experiment, Sicot 71BRF and 5B (a lintless variant with compromised ethylene metabolism) were exposed to 45°C, either by instantaneous heat shock or by ramping temperatures by 3°C daily for 1 week. One day prior to the start of heat treatment, half the plants were sprayed with 0.8 mM of the ethylene synthesis inhibitor, aminoethoxyvinylglycine (AVG). In a subsequent experiment, Sicot 71BRF and a putatively heat-tolerant line, CIM 448, were exposed to 36 or 45°C for 1 week, and half the plants were sprayed with 20 μM of the ethylene precursor, 1-aminocyclopropane-1-carboxylic acid, (ACC). High temperature exposure of plants in both experiments was performed at the peak reproductive phase (65-68 days after sowing). Elevated temperature (heat shock or ramping to 45°C) significantly reduced production and retention of fruits in all cotton lines used in this study. At the termination of heat treatment, cotton plants exposed to 45°C had at least 50% fewer fruits than plants under optimum temperature in all three genotypes, while plants at 36°C remained unaffected. Heat-stressed plants continued producing new squares (fruiting buds) after termination of heat stress but these squares did not turn into cotton bolls due to pollen infertility. In vitro inhibition of pollen germination by high temperatures supported this observation. Leaf photosynthesis (Pn) of heat-stressed plants (45°C) measured at the end of heat treatments remained significantly inhibited, despite an increased leaf stomatal conductance (gs), suggesting that high temperature impairs Pn independently of stomatal behavior. Metabolic injury was supported by high relative cellular injury and low photosystem II yield of the heat-stressed plants, indicating that high temperature impaired photosynthetic electron transport. Both heat shock and ramping of heat significantly reduced ethylene release from cotton leaf tissues measured at the end of heat treatment but modulating ethylene production via AVG or ACC application had no significant effect on fruit production or retention in heat-stressed cotton plants. Instead, high temperature accelerated fruit abortion by impairing pollen development and/or restricting leaf photosynthesis.

Role of xyloglucan in cotton (Gossypium hirsutum L.) fiber elongation of the short fiber mutant Ligon lintless-2 (Li2)

Xyloglucan is a matrix polysaccharide found in the cell walls of all land plants. In growing cells, xyloglucan is thought to connect cellulose microfibrils and regulate their separation during wall extension. Ligon lintless-2 (Li₂) is a monogenic dominant cotton fiber mutation that causes extreme reduction in lint fiber length with no pleiotropic effects on vegetative growth. Li₂ represents an excellent model system to study fiber elongation. To understand the role of xyloglucan in cotton fiber elongation we used the short fiber mutant Li₂ and its near isogenic wild type for analysis of xyloglucan content and expression of xyloglucan-related genes in developing fibers. Accumulation of xyloglucan was significantly higher in Li₂ developing fibers than in wild type. Genes encoding enzymes for nine family members of xyloglucan biosynthesis were identified in the draft Gossypium hirsutum genome. RNAseq analysis revealed that most differentially expressed xyloglucan-related genes were down-regulated in Li₂ fiber cells. RT-qPCR analysis revealed that the peak of expression for the majority of xyloglucan-related genes in wild type developing fibers was 5-16days post anthesis (DPA) compared to 1-3 DPA in Li₂ fibers. Thus, our results suggest that early activation of xyloglucan-related genes and down regulation of xyloglucan degradation genes during the elongation phase lead to elevated accumulation of xyloglucan that restricts elongation of fiber cells in Li₂.

Genome-wide analysis of gene expression of EMS-induced short fiber mutant Ligon lintless-y (liy) in cotton (Gossypium hirsutum L.)

In this work we describe a chemically-induced short fiber mutant cotton line, Ligon-lintless-y (li_y), which is controlled by a single recessive locus and affects multiple traits, including height of the plant, and length and maturity of fiber. An RNAseq analysis was used to evaluate global transcriptional changes during cotton fiber development at 3, 8 and 16days post anthesis. We found that 613, 2629 and 3397 genes were significantly down-regulated, while 2700, 477 and 3260 were significantly up-regulated in li_y at 3, 8 and 16 DPA. Gene set enrichment analysis revealed that many metabolic pathways, including carbohydrate, cell wall, hormone metabolism and transport were substantially altered in li_y developing fibers. We discuss perturbed expression of genes involved in signal transduction and biosynthesis of phytohormones, such as auxin, abscisic acid, gibberellin and ethylene. The results of this study provide new insights into transcriptional regulation of cotton fiber development.

A Gly65Val substitution in an actin, GhACT_LI1, disrupts cell polarity and F-actin organization resulting in dwarf, lintless cotton plants

Actin polymerizes to form part of the cytoskeleton and organize polar growth in all eukaryotic cells. Species with numerous actin genes are especially useful for the dissection of actin molecular function due to redundancy and neofunctionalization. Here, we investigated the role of a cotton (Gossypium hirsutum) actin gene in the organization of actin filaments in lobed cotyledon pavement cells and the highly elongated single-celled trichomes that comprise cotton lint fibers. Using mapping-by-sequencing, virus-induced gene silencing, and molecular modeling, we identified the causative mutation of the dominant dwarf Ligon lintless Li₁ short fiber mutant as a single Gly65Val amino acid substitution in a polymerization domain of an actin gene, GhACT_LI1 (Gh_D04G0865). We observed altered cell morphology and disrupted organization of F-actin in Li₁ plant cells by confocal microscopy. Mutant leaf cells lacked interdigitation of lobes and F-actin did not uniformly decorate the nuclear envelope. While wild-type lint fiber trichome cells contained long longitudinal actin cables, the short Li₁ fiber cells accumulated disoriented transverse cables. The polymerization-defective Gly65Val allele in Li₁ plants likely disrupts processive elongation of F-actin, resulting in a disorganized cytoskeleton and reduced cell polarity, which likely accounts for the dominant gene action and diverse pleiotropic effects associated with the Li₁ mutation. Lastly, we propose a model to account for these effects, and underscore the roles of actin organization in determining plant cell polarity, shape and plant growth.

A comparative transcriptome analysis of two sets of backcross inbred lines differing in lint-yield derived from a Gossypium hirsutum × Gossypium barbadense population

Upland cotton (Gossypium hirsutum L.) is the most important fiber crop, and its lint-yield improvement is impeded due to its narrow genetic base and the lack of understanding of the genetic basis of yield. Backcross inbred lines (BILs) or near-isogenic lines (NILs) in the same genetic background differing in lint yield, developed through advanced backcrossing, provide an important genomic resource to study the molecular genetic basis of lint yield. In the present study, a high-yield (HY) group and a low-yield (LY) group each with three BILs were selected from a BIL population between G. hirsutum and G. barbadense. Using a microarray-based comparative transcriptome analysis on developing fibers at 10 days post-anthesis (DPA) between the two groups, 1486 differentially expressed genes (DEGs) were identified. A total of 212 DEGs were further mapped in the regions of 24 yield QTL and 11 yield trait QTL hotspots as reported previously, and 81 DEGs mapped with the 7 lint-yield QTL identified in the BIL population from which the two sets of BILs were selected. Gene Ontology annotations and Blast-Mapping-Annotation-KEGG analysis via Blast2GO revealed that more DEGs were associated with catalytic activity and binding, followed by transporters, nucleic acid binding transcription factors, structural molecules and molecular transducer activities. Six DEGs were chosen for a quantitative RT-PCR assay, and the results were consistent with the microarray analysis. The development of DEGs-based markers revealed that 7 single strand conformation polymorphism-based single nucleotide polymorphic (SSCP-SNP) markers were associated with yield traits, and 3 markers with lint yield. In the present study, we identified a number of yield and yield component QTL-co-localizing DEGs and developed several DEG-based SSCP-SNP markers for the traits, thereby providing a set of candidate genes for molecular breeding and genetic manipulation of lint yield in cotton.

Small RNA sequencing and degradome analysis of developing fibers of short fiber mutants Ligon-lintles-1 (Li 1 ) and -2 (Li 2 ) revealed a role for miRNAs and their targets in cotton fiber elongation

BACKGROUND

The length of cotton fiber is an important agronomic trait that directly affects the quality of yarn and fabric. Understanding the molecular basis of fiber elongation would provide a means for improvement of fiber length. Ligon-lintless-1 (Li 1 ) and -2 (Li 2 ) are monogenic and dominant mutations that result in an extreme reduction in the length of lint fiber on mature seeds. In a near-isogenic state with wild type cotton these two short fiber mutants provide an effective model system to study the mechanisms of fiber elongation. Plant miRNAs regulate many aspects of growth and development. However, the mechanism underlying the miRNA-mediated regulation of fiber development is largely unknown.

RESULTS

Small RNA libraries constructed from developing fiber cells of the short fiber mutants Li 1 and Li 2 and their near-isogenic wild type lines were sequenced. We identified 24 conservative and 147 novel miRNA families with targets that were detected through degradome sequencing. The distribution of the target genes into functional categories revealed the largest set of genes were transcription factors. Expression profiles of 20 miRNAs were examined across a fiber developmental time course in wild type and short fiber mutations. We conducted correlation analysis between miRNA transcript abundance and the length of fiber for 11 diverse Upland cotton lines. The expression patterns of 4 miRNAs revealed significant negative correlation with fiber lengths of 11 cotton lines.

CONCLUSIONS

Our results suggested that the mutations have changed the regulation of miRNAs expression during fiber development. Further investigations of differentially expressed miRNAs in the Li 1 and Li 2 mutants will contribute to better understanding of the regulatory mechanisms of cotton fiber development. Four miRNAs negatively correlated with fiber length are good candidates for further investigations of miRNA regulation of important genotype dependent fiber traits. Thus, our results will contribute to further studies on the role of miRNAs in cotton fiber development and will provide a tool for fiber improvement through molecular breeding.

Fine mapping and RNA-Seq unravels candidate genes for a major QTL controlling multiple fiber quality traits at the T1 region in upland cotton

BACKGROUND

Improving fiber quality is a major challenge in cotton breeding, since the molecular basis of fiber quality traits is poorly understood. Fine mapping and candidate gene prediction of quantitative trait loci (QTL) controlling cotton fiber quality traits can help to elucidate the molecular basis of fiber quality. In our previous studies, one major QTL controlling multiple fiber quality traits was identified near the T1 locus on chromosome 6 in Upland cotton.

RESULTS

To finely map this major QTL, the F2 population with 6975 individuals was established from a cross between Yumian 1 and a recombinant inbred line (RIL118) selected from a recombinant inbred line population (T586 × Yumian 1). The QTL was mapped to a 0.28-cM interval between markers HAU2119 and SWU2302. The QTL explained 54.7 % (LOD = 222.3), 40.5 % (LOD = 145.0), 50.0 % (LOD = 194.3) and 30.1 % (LOD = 100.4) of phenotypic variation with additive effects of 2.78, -0.43, 2.92 and 1.90 units for fiber length, micronaire, strength and uniformity, respectively. The QTL region corresponded to a 2.7-Mb interval on chromosome 10 in the G. raimondii genome sequence and a 5.3-Mb interval on chromosome A06 in G. hirsutum. The fiber of Yumian 1 was much longer than that of RIL118 from 3 DPA to 7 DPA. RNA-Seq of ovules at 0 DPA and fibers at 5 DPA from Yumian 1 and RIL118 showed four genes in the QTL region of the G. raimondii genome to be extremely differentially expressed. RT-PCR analysis showed three genes in the QTL region of the G. hirsutum genome to behave similarly.

CONCLUSIONS

This study mapped a major QTL influencing four fiber quality traits to a 0.28-cM interval and identified three candidate genes by RNA-Seq and RT-PCR analysis. Integration of fine mapping and RNA-Seq is a powerful strategy to uncover candidates for QTL in large genomes.

Transcriptome Analysis of Short Fiber Mutant Ligon lintless-1 (Li1) Reveals Critical Genes and Key Pathways in Cotton Fiber Elongation and Leaf Development

For efficient spinning and superior fabric production, long fiber length is a desired trait for cotton production. To unveil the molecular basis of the cotton fiber length regulation, a short fiber mutant, Ligon lintless-1 (Li₁), is selected to compare with its corresponding wild type (WT). Li₁ is a monogenic dominant cotton mutant causing extremely short fibers (<6mm) on mature seeds with visible pleiotropic effects on vegetative growth and development. In this research, we compared the transcriptome of fiber bearing ovules at 1 DPA, 3 DPA, 8 DPA and leaf between Li₁ mutant and WT. A total of 7,852 differentially expressed genes (DEGs) were detected in ovules and leaves, which mainly participated in sugar, secondary metabolite and lipid metabolism pathways based on KEGG analysis. The common DEGs at 1 DPA and 3 DPA were involved in the responses to endogenous stimulus, signal transduction and long-chain fatty acid biosynthesis. For 3 DPA, 8 DPA and leaf, the common DEGs were involved in the responses to auxin and receptor kinases related pathway. Further analysis showed that 37 genes involved in very-long-chain fatty acid biosynthesis were suppressed in Li₁ mutant during fiber fast elongation development. Most of the DEGs involved in cell wall metabolism, such cellulose synthase, expansin family, and glycosyl hydrolase were differentially expressed at 3 DPA and 8 DPA. Our results provide new insights into the mechanisms of fiber elongation, and offer novel genes as potential objects for fiber length improvement.

Mapping-by-sequencing of Ligon-lintless-1 (Li 1 ) reveals a cluster of neighboring genes with correlated expression in developing fibers of Upland cotton (Gossypium hirsutum L.)

Mapping-by-sequencing and SNP marker analysis were used to fine map the Ligon-lintless-1 ( Li 1 ) short fiber mutation in tetraploid cotton to a 255-kb region that contains 16 annotated proteins. The Ligon-lintless-1 (Li 1 ) mutant of cotton (Gossypium hirsutum L.) has been studied as a model for cotton fiber development since its identification in 1929; however, the causative mutation has not been identified yet. Here we report the fine genetic mapping of the mutation to a 255-kb region that contains only 16 annotated genes in the reference Gossypium raimondii genome. We took advantage of the incompletely dominant dwarf vegetative phenotype to identify 100 mutants (Li 1 /Li 1 ) and 100 wild-type (li 1 /li 1 ) homozygotes from a mapping population of 2567 F2 plants, which we bulked and deep sequenced. Since only homozygotes were sequenced, we were able to use a high stringency in SNP calling to rapidly narrow down the region harboring the Li 1 locus, and designed subgenome-specific SNP markers to test the population. We characterized the expression of all sixteen genes in the region by RNA sequencing of elongating fibers and by RT-qPCR at seven time points spanning fiber development. One of the most highly expressed genes found in this interval in wild-type fiber cells is 40-fold under-expressed at the day of anthesis (DOA) in the mutant fiber cells.  This gene is a major facilitator superfamily protein, part of the large family of proteins that includes auxin and sugar transporters. Interestingly, nearly all genes in this region were most highly expressed at DOA and showed a high degree of co-expression. Further characterization is required to determine if transport of hormones or carbohydrates is involved in both the dwarf and lintless phenotypes of Li 1 plants.

Genetic fine mapping and candidate gene analysis of the Gossypium hirsutum Ligon lintless-1 (Li1) mutant on chromosome 22(D)

Ligon lintless-1 (Li1) is a Gossypium hirsutum mutant that is controlled by a dominant gene that arrests the development of cotton fiber after anthesis. Two F2 mapping populations were developed from mutant (Li1 × H7124) F1 plants in 2012 and 2013; each was composed of 142 and 1024 plants, respectively. Using these populations, Li1 was mapped to a 0.3-cM region in which nine single-strand conformation polymorphism markers co-segregated with the Li1 locus. In the published G. raimondii genome, these markers were mapped to a region of about 1.2 Mb (the Li1 region) and were separated by markers that flanked the Li1 locus in the genetic map, dividing the Li1 region into three segments. Thirty-six genes were annotated by the gene prediction software FGENESH (Softberry) in the Li1 region. Twelve genes were candidates of Li1, while the remaining 24 genes were identified as transposable elements, DNA/RNA polymerase superfamily or unknown function genes. Among the 12 candidate genes, those encoding ribosomal protein s10, actin protein, ATP synthase, and beta-tubulin 5 were the most-promising candidates of the Li1 mutant because the function of these genes is closely related to fiber development. High-throughput RNA sequencing and quantitative PCR revealed that these candidate genes had obvious differential gene expression between mutant and wild-type plants at the fiber elongation stage, strengthening the inference that they could be the most likely candidate gene of the Li1 mutant phenotype.

RNA-seq analysis of short fiber mutants Ligon-lintless-1 (Li 1 ) and - 2 (Li 2 ) revealed important role of aquaporins in cotton (Gossypium hirsutum L.) fiber elongation

BACKGROUND

Cotton fiber length is a key determinant of fiber quality for the textile industry. Understanding the molecular basis of fiber elongation would provide a means for improvement of fiber length. Ligon lintless-1 (Li 1 ) and Ligon lintless-2 (Li 2 ) are monogenic and dominant mutations, that result in an extreme reduction in the length of lint fiber to approximately 6 mm on mature seeds. In a near-isogenic state with wild type (WT) cotton these two short fiber mutants provide an excellent model system to study mechanisms of fiber elongation.

RESULTS

We used next generation sequencing (RNA-seq) to identify common fiber elongation related genes in developing fibers of Li 1 and Li 2 mutants growing in the field and a greenhouse. We found a large number of differentially expressed genes common to both mutants, including 531 up-regulated genes and 652 down-regulated genes. Major intrinsic proteins or aquaporins were one of the most significantly over-represented gene families among common down-regulated genes in Li 1 and Li 2 fibers. The members of three subfamilies of aquaporins, including plasma membrane intrinsic proteins, tonoplast intrinsic proteins and NOD26-like intrinsic proteins were down-regulated in short fiber mutants. The osmotic concentration and the concentrations of soluble sugars were lower in fiber cells of both short fiber mutants than in WT, whereas the concentrations of K+ and malic acid were significantly higher in mutants during rapid cell elongation.

CONCLUSIONS

We found that the aquaporins were the most down-regulated gene family in both short fiber mutants. The osmolality and concentrations of soluble sugars were less in saps of Li 1 - Li 2 , whereas the concentrations of malic acid, K+ and other detected ions were significantly higher in saps of mutants than in WT. These results suggest that higher accumulation of ions in fiber cells, reduced osmotic pressure and low expression of aquaporins, may contribute to the cessation of fiber elongation in Li 1 and Li 2 short-fiber mutants. The research presented here provides new insights into osmoregulation of short fiber mutants and the role of aquaporins in cotton fiber elongation.

Understanding the relationship between cotton fiber properties and non-cellulosic cell wall polysaccharides

A detailed knowledge of cell wall heterogeneity and complexity is crucial for understanding plant growth and development. One key challenge is to establish links between polysaccharide-rich cell walls and their phenotypic characteristics. It is of particular interest for some plant material, like cotton fibers, which are of both biological and industrial importance. To this end, we attempted to study cotton fiber characteristics together with glycan arrays using regression based approaches. Taking advantage of the comprehensive microarray polymer profiling technique (CoMPP), 32 cotton lines from different cotton species were studied. The glycan array was generated by sequential extraction of cell wall polysaccharides from mature cotton fibers and screening samples against eleven extensively characterized cell wall probes. Also, phenotypic characteristics of cotton fibers such as length, strength, elongation and micronaire were measured. The relationship between the two datasets was established in an integrative manner using linear regression methods. In the conducted analysis, we demonstrated the usefulness of regression based approaches in establishing a relationship between glycan measurements and phenotypic traits. In addition, the analysis also identified specific polysaccharides which may play a major role during fiber development for the final fiber characteristics. Three different regression methods identified a negative correlation between micronaire and the xyloglucan and homogalacturonan probes. Moreover, homogalacturonan and callose were shown to be significant predictors for fiber length. The role of these polysaccharides was already pointed out in previous cell wall elongation studies. Additional relationships were predicted for fiber strength and elongation which will need further experimental validation.

Related to ABA-Insensitive3(ABI3)/Viviparous1 and AtABI5 transcription factor coexpression in cotton enhances drought stress adaptation

Drought tolerance is an important trait being pursued by the agbiotech industry. Abscisic acid (ABA) is a stress hormone that mediates a multitude of processes in growth and development, water use efficiency (WUE) and gene expression during seed development and in response to environmental stresses. Arabidopsis B3-domain transcription factor Related to ABA-Insensitive3 (ABI3)/Viviparous1 (namely AtRAV2) and basic leucine zipper (bZIPs) AtABI5 or AtABF3 transactivated ABA-inducible promoter:GUS reporter expression in a maize mesophyll protoplast transient assay and showed synergies in reporter transactivation when coexpressed. Transgenic cotton (Gossypium hirsutum) expressing AtRAV1/2 and/or AtABI5 showed resistance to imposed drought stress under field and greenhouse conditions and exhibited improved photosynthesis and WUEs associated with absorption through larger root system and greater leaf area. We observed synergy for root biomass accumulation in the greenhouse, intrinsic WUE in the field and drought tolerance in stacked AtRAV and AtABI5 double-transgenic cotton. We assessed AtABI5 and AtRAV1/2 involvement in drought stress adaptations through reactive oxygen species scavenging and osmotic adjustment by marker gene expression in cotton. Deficit irrigation-grown AtRAV1/2 and AtABI5 transgenics had 'less-stressed' molecular and physiological phenotypes under drought, likely due to improved photoassimilation and root and shoot sink strengths and enhanced expression of endogenous GhRAV and genes for antioxidant and osmolyte biosynthesis. Overexpression of bZIP and RAV TFs could impact sustainable cotton agriculture and potentially other crops under limited irrigation conditions.

Quantitative trait loci analysis of fiber quality traits using a random-mated recombinant inbred population in Upland cotton (Gossypium hirsutum L.)

BACKGROUND

Upland cotton (Gossypium hirsutum L.) accounts for about 95% of world cotton production. Improving Upland cotton cultivars has been the focus of world-wide cotton breeding programs. Negative correlation between yield and fiber quality is an obstacle for cotton improvement. Random-mating provides a potential methodology to break this correlation. The suite of fiber quality traits that affect the yarn quality includes the length, strength, maturity, fineness, elongation, uniformity and color. Identification of stable fiber quantitative trait loci (QTL) in Upland cotton is essential in order to improve cotton cultivars with superior quality using marker-assisted selection (MAS) strategy.

RESULTS

Using 11 diverse Upland cotton cultivars as parents, a random-mated recombinant inbred (RI) population consisting of 550 RI lines was developed after 6 cycles of random-mating and 6 generations of self-pollination. The 550 RILs were planted in triplicates for two years in Mississippi State, MS, USA to obtain fiber quality data. After screening 15538 simple sequence repeat (SSR) markers, 2132 were polymorphic among the 11 parents. One thousand five hundred eighty-two markers covering 83% of cotton genome were used to genotype 275 RILs (Set 1). The marker-trait associations were analyzed using the software program TASSEL. At p < 0.01, 131 fiber QTLs and 37 QTL clusters were identified. These QTLs were responsible for the combined phenotypic variance ranging from 62.3% for short fiber content to 82.8% for elongation. The other 275 RILs (Set 2) were analyzed using a subset of 270 SSR markers, and the QTLs were confirmed. Two major QTL clusters were observed on chromosomes 7 and 16. Comparison of these 131 QTLs with the previously published QTLs indicated that 77 were identified before, and 54 appeared novel.

CONCLUSIONS

The 11 parents used in this study represent a diverse genetic pool of the US cultivated cotton, and 10 of them were elite commercial cultivars. The fiber QTLs, especially QTL clusters reported herein can be readily implemented in a cotton breeding program to improve fiber quality via MAS strategy. The consensus QTL regions warrant further investigation to better understand the genetics and molecular mechanisms underlying fiber development.

Quantitative trait loci analysis of fiber quality traits using a random-mated recombinant inbred population in Upland cotton (Gossypium hirsutum L.)

BACKGROUND

Upland cotton (Gossypium hirsutum L.) accounts for about 95% of world cotton production. Improving Upland cotton cultivars has been the focus of world-wide cotton breeding programs. Negative correlation between yield and fiber quality is an obstacle for cotton improvement. Random-mating provides a potential methodology to break this correlation. The suite of fiber quality traits that affect the yarn quality includes the length, strength, maturity, fineness, elongation, uniformity and color. Identification of stable fiber quantitative trait loci (QTL) in Upland cotton is essential in order to improve cotton cultivars with superior quality using marker-assisted selection (MAS) strategy.

RESULTS

Using 11 diverse Upland cotton cultivars as parents, a random-mated recombinant inbred (RI) population consisting of 550 RI lines was developed after 6 cycles of random-mating and 6 generations of self-pollination. The 550 RILs were planted in triplicates for two years in Mississippi State, MS, USA to obtain fiber quality data. After screening 15538 simple sequence repeat (SSR) markers, 2132 were polymorphic among the 11 parents. One thousand five hundred eighty-two markers covering 83% of cotton genome were used to genotype 275 RILs (Set 1). The marker-trait associations were analyzed using the software program TASSEL. At p < 0.01, 131 fiber QTLs and 37 QTL clusters were identified. These QTLs were responsible for the combined phenotypic variance ranging from 62.3% for short fiber content to 82.8% for elongation. The other 275 RILs (Set 2) were analyzed using a subset of 270 SSR markers, and the QTLs were confirmed. Two major QTL clusters were observed on chromosomes 7 and 16. Comparison of these 131 QTLs with the previously published QTLs indicated that 77 were identified before, and 54 appeared novel.

CONCLUSIONS

The 11 parents used in this study represent a diverse genetic pool of the US cultivated cotton, and 10 of them were elite commercial cultivars. The fiber QTLs, especially QTL clusters reported herein can be readily implemented in a cotton breeding program to improve fiber quality via MAS strategy. The consensus QTL regions warrant further investigation to better understand the genetics and molecular mechanisms underlying fiber development.

Quantitative trait loci analysis of fiber quality traits using a random-mated recombinant inbred population in Upland cotton (Gossypium hirsutum L.)

BACKGROUND

Upland cotton (Gossypium hirsutum L.) accounts for about 95% of world cotton production. Improving Upland cotton cultivars has been the focus of world-wide cotton breeding programs. Negative correlation between yield and fiber quality is an obstacle for cotton improvement. Random-mating provides a potential methodology to break this correlation. The suite of fiber quality traits that affect the yarn quality includes the length, strength, maturity, fineness, elongation, uniformity and color. Identification of stable fiber quantitative trait loci (QTL) in Upland cotton is essential in order to improve cotton cultivars with superior quality using marker-assisted selection (MAS) strategy.

RESULTS

Using 11 diverse Upland cotton cultivars as parents, a random-mated recombinant inbred (RI) population consisting of 550 RI lines was developed after 6 cycles of random-mating and 6 generations of self-pollination. The 550 RILs were planted in triplicates for two years in Mississippi State, MS, USA to obtain fiber quality data. After screening 15538 simple sequence repeat (SSR) markers, 2132 were polymorphic among the 11 parents. One thousand five hundred eighty-two markers covering 83% of cotton genome were used to genotype 275 RILs (Set 1). The marker-trait associations were analyzed using the software program TASSEL. At p < 0.01, 131 fiber QTLs and 37 QTL clusters were identified. These QTLs were responsible for the combined phenotypic variance ranging from 62.3% for short fiber content to 82.8% for elongation. The other 275 RILs (Set 2) were analyzed using a subset of 270 SSR markers, and the QTLs were confirmed. Two major QTL clusters were observed on chromosomes 7 and 16. Comparison of these 131 QTLs with the previously published QTLs indicated that 77 were identified before, and 54 appeared novel.

CONCLUSIONS

The 11 parents used in this study represent a diverse genetic pool of the US cultivated cotton, and 10 of them were elite commercial cultivars. The fiber QTLs, especially QTL clusters reported herein can be readily implemented in a cotton breeding program to improve fiber quality via MAS strategy. The consensus QTL regions warrant further investigation to better understand the genetics and molecular mechanisms underlying fiber development.

Comparative transcriptome analysis of short fiber mutants Ligon-lintless 1 and 2 reveals common mechanisms pertinent to fiber elongation in cotton (Gossypium hirsutum L.)

Understanding the molecular processes affecting cotton (Gossypium hirsutum) fiber development is important for developing tools aimed at improving fiber quality. Short fiber cotton mutants Ligon-lintless 1 (Li₁) and Ligon-lintless 2 (Li₂) are naturally occurring, monogenic mutations residing on different chromosomes. Both mutations cause early cessation in fiber elongation. These two mutants serve as excellent model systems to elucidate molecular mechanisms relevant to fiber length development. Previous studies of these mutants using transcriptome analysis by our laboratory and others had been limited by the fact that very large numbers of genes showed altered expression patterns in the mutants, making a targeted analysis difficult or impossible. In this research, a comparative microarray analysis was conducted using these two short fiber mutants and their near isogenic wild type (WT) grown under both field and greenhouse environments in order to identify key genes or metabolic pathways common to fiber elongation. Analyses of three transcriptome profiles obtained from different growth conditions and mutant types showed that most differentially expressed genes (DEGs) were affected by growth conditions. Under field conditions, short fiber mutants commanded higher expression of genes related to energy production, manifested by the increasing of mitochondrial electron transport activity or responding to reactive oxygen species when compared to the WT. Eighty-eight DEGs were identified to have altered expression patterns common to both short fiber mutants regardless of growth conditions. Enrichment, pathway and expression analyses suggested that these 88 genes were likely involved in fiber elongation without being affected by growth conditions.

Elucidation of Nuclear and Organellar Genomes of Gossypium hirsutum: Furthering Studies of Species Evolution and Applications for Crop Improvement

Plant genomes are larger and more complex than other eukaryotic organisms, due to small and large duplication events, recombination and subsequent reorganization of the genetic material. Commercially important cotton is the result of a polyploidization event between Old and New World cottons that occurred over one million years ago. Allotetraploid cotton has properties that are dramatically different from its progenitors-most notably, the presence of long, spinnable fibers. Recently, the complete genome of a New World cotton ancestral species, Gossypium raimondii, was completed. Future genome sequencing efforts are focusing on an Old World progenitor, G. arboreum. This sequence information will enable us to gain insights into the evolution of the cotton genome that may be used to understand the evolution of other plant species. The chloroplast genomes of multiple cotton species and races have been determined. This information has also been used to gain insight into the evolutionary history of cotton. Analysis of the database of nuclear and organellar sequences will facilitate the identification of potential genes of interest and subsequent development of strategies for improving cotton.