Genetic mapping and comparative analysis of seven mutants related to seed fiber development in cotton

Two duplicated GhMML3 genes coordinately control development of lint and fuzz fibers in cotton

Cotton produces two types of fibers: fuzz and lint. Cotton yield is determined by the number of epidermal cells that develop into lint fibers. Despite numerous studies, the genetic and molecular mechanisms that control lint and fuzz fiber development remain unclear. Here, using the recessive naked-seed or fuzzless-linted mutant (n2NSM) in combination with gene editing and complementation, we found that the recessive fuzzless gene n2 encodes the MYBMIXTA-like (MML) transcription factor GhMML3_D12. Overexpression of GhMML3_D12 in n2NSM restored fuzz fiber development, whereas CRISPR-Cas9 knockout of GhMML3_D12 in wild-type cotton (J668) resulted in a fuzzless-linted phenotype. Interestingly, simultaneous edits to GhMML3_D12 and its duplicate GhMML3_A12 resulted in plants with a fiberless (fuzzless-lintless) phenotype. Detailed investigation of the seed fiber phenotypes of segregating progeny derived from a cross between J668 and a fiberless gene-edited mutant of GhMML3 (#mml3s) not only identified progeny that mimicked natural fuzzless and fiberless mutants but also revealed that the duplicated GhMML3_A12 and GhMML3_D12 regulate the development of fuzz and lint fibers in a dose-dependent manner. Comparative transcriptome analysis and single-cell RNA sequencing identified GhMML3 as the central hub of the gene network that regulates fiber initiation and early-stage elongation. The gene regulatory network revealed potential candidate genes and key regulators that may contribute to fiber initiation and development, and a model for the control of lint and fuzz fiber development by GhMML3 was proposed. We also found that the GhMML3_D12 protein can bind directly to the promoters of GhHD-1 and GhMYB25, two key genes involved in fiber initiation, thereby activating their expression. This study provides new insights into the fundamental mechanisms that underlie cotton fiber development.

Genome-Wide and Expression Pattern Analysis of the DVL Gene Family Reveals GhM_A05G1032 Is Involved in Fuzz Development in G. hirsutum

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

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.

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.

Inheritance, QTLs, and Candidate Genes of Lint Percentage in Upland Cotton

Cotton (Gossypium spp.) is an important natural fiber plant. Lint percentage (LP) is one of the most important determinants of cotton yield and is a typical quantitative trait with high variation and heritability. Many cotton LP genetic linkages and association maps have been reported. This work summarizes the inheritance, quantitative trait loci (QTLs), and candidate genes of LP to facilitate LP genetic study and molecular breeding. More than 1439 QTLs controlling LP have been reported. Excluding replicate QTLs, 417 unique QTLs have been identified on 26 chromosomes, including 243 QTLs identified at LOD >3. More than 60 are stable, major effective QTLs that can be used in marker-assisted selection (MAS). More than 90 candidate genes for LP have been reported. These genes encode MYB, HOX, NET, and other proteins, and most are preferentially expressed during fiber initiation and elongation. A putative molecular regulatory model of LP was constructed and provides the foundation for the genetic study and molecular breeding of LP.

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.

Large-fragment insertion activates gene GaFZ (Ga08G0121) and is associated with the fuzz and trichome reduction in cotton (Gossypium arboreum)

Cotton seeds are typically covered by lint and fuzz fibres. Natural 'fuzzless' mutants are an ideal model system for identifying genes that regulate cell initiation and elongation. Here, using a genome-wide association study (GWAS), we identified a ~ 6.2 kb insertion, larINDEL_FZ , located at the end of chromosome 8, composed of a ~ 5.0 kb repetitive sequence and a ~ 1.2 kb fragment translocated from chromosome 12 in fuzzless Gossypium arboreum. The presence of larINDEL_FZ was associated with a fuzzless seed and reduced trichome phenotypes in G. arboreum. This distant insertion was predicted to be an enhancer, located ~ 18 kb upstream of the dominant-repressor GaFZ (Ga08G0121). Ectopic overexpression of GaFZ in Arabidopsis thaliana and G. hirsutum suggested that GaFZ negatively modulates fuzz and trichome development. Co-expression and interaction analyses demonstrated that GaFZ might impact fuzz fibre/trichome development by repressing the expression of genes in the very-long-chain fatty acid elongation pathway. Thus, we identified a novel regulator of fibre/trichome development while providing insights into the importance of noncoding sequences in cotton.

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.

A Retrotransposon Insertion in GhMML3_D12 Is Likely Responsible for the Lintless Locus li3 of Tetraploid Cotton

Cotton (Gossypium) seed fibers can be divided into lint (long) or fuzz (very short). Using fiberless (fuzzless-lintless) mutants, the lint initiation gene Li₃ was identified by map-based cloning. The gene is an R2R3-MYB transcription factor located on chromosome D12 (GhMML3_D12). Sequence analysis revealed that li₃ is a loss-of-function allele containing a retrotransposon insertion in the second exon that completely blocks the gene’s expression. The genetic loci n₂ and n₃ underlying the recessive fuzzless phenotype in Gossypium hirsutum were also mapped. The genomic location of n₃ overlapped with that of the dominant fuzzless locus N₁, and n₃ appeared to be a loss-of-function allele caused by a single nucleotide polymorphism (SNP) mutation in the coding region of GhMML3_A12. The n₂ allele was found to be co-located with li₃ and originated from G. babardense. n₂ and li₃ are possibly the multiple alleles of the GhMML3_D12 gene. Genetic analysis showed that Li₃ and N₃ are a pair of homologs with additive effects for the initiation of fibers (fuzz or lint). In addition, the presence of another locus was speculated, and it appeared to show an inhibitory effect on the expression of GhMML3. These findings provide new information about the genetic factors affecting the initiation of fibers in cotton.

Improved reconstruction and comparative analysis of chromosome 12 to rectify Mis-assemblies in Gossypium arboreum

BACKGROUND

Genome sequencing technologies have been improved at an exponential pace but precise chromosome-scale genome assembly still remains a great challenge. The draft genome of cultivated G. arboreum was sequenced and assembled with shotgun sequencing approach, however, it contains several misassemblies. To address this issue, we generated an improved reassembly of G. arboreum chromosome 12 using genetic mapping and reference-assisted approaches and evaluated this reconstruction by comparing with homologous chromosomes of G. raimondii and G. hirsutum.

RESULTS

In this study, we generated a high quality assembly of the 94.64 Mb length of G. arboreum chromosome 12 (A_A12) which comprised of 144 scaffolds and contained 3361 protein coding genes. Evaluation of results using syntenic and collinear analysis of reconstructed G. arboreum chromosome A_A12 with its homologous chromosomes of G. raimondii (D_D08) and G. hirsutum (AD_A12 and AD_D12) confirmed the significant improved quality of current reassembly as compared to previous one. We found major misassemblies in previously assembled chromosome 12 (A_Ca9) of G. arboreum particularly in anchoring and orienting of scaffolds into a pseudo-chromosome. Further, homologous chromosomes 12 of G. raimondii (D_D08) and G. arboreum (A_A12) contained almost equal number of transcription factor (TF) related genes, and showed good collinear relationship with each other. As well, a higher rate of gene loss was found in corresponding homologous chromosomes of tetraploid (AD_A12 and AD_D12) than diploid (A_A12 and D_D08) cotton, signifying that gene loss is likely a continuing process in chromosomal evolution of tetraploid cotton.

CONCLUSION

This study offers a more accurate strategy to correct misassemblies in sequenced draft genomes of cotton which will provide further insights towards its genome organization.

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.

Genetic Identification and Transcriptome Analysis of Lintless and Fuzzless Traits in Gossypium arboreum L

Cotton fibres, as single cells arising from the seed coat, can be classified as lint and fuzz according to their final length. Gossypium arboreum is a cultivated diploid cotton species and a potential donor of the A subgenome of the more widely grown tetraploid cottons. In this study, we performed genetic studies on one lintless and seven fuzzless G. arboreum accessions. Through association and genetic linkage analyses, a recessive locus on Chr06 containing GaHD-1 was found to be the likely gene underlying the lintless trait. GaHD-1 carried a mutation at a splicing acceptor site that resulted in alternative splicing and a deletion of 247 amino acid from the protein. The regions containing GaGIR1 and GaMYB25-like were found to be associated with fuzz development in G. arboreum, with the former being the major contributor. Comparative transcriptome analyses using 0-5 days post-anthesis (dpa) ovules from lintless, fuzzless, and normal fuzzy seed G. arboreum accessions revealed gene modules and hub genes potentially important for lint and fuzz initiation and development. Three significant modules and 26 hub genes associated with lint fibre initiation were detected by weighted gene co-expression network analysis. Similar analyses identified three vital modules and 10 hub genes to be associated with fuzz development. The findings in this study contribute to understanding the complex molecular mechanism(s) regulating fibre initiation and development and indicate that G. arboreum may have fibre developmental pathways different from tetraploid cotton. It also provides candidate genes for further investigation into modifying fibre development in G. arboreum.

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.

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.

Fine mapping and identification of the fuzzless gene GaFzl in DPL972 (Gossypium arboreum)

KEY MESSAGE

The fuzzless gene GaFzl was fine mapped to a 70-kb region containing a GIR1 gene, Cotton_A_11941, responsible for the fuzzless trait in Gossypium arboreum DPL972. Cotton fiber is the most important natural textile resource. The fuzzless mutant DPL972 (Gossypium arboreum) provides a useful germplasm resource to explore the molecular mechanism underlying fiber and fuzz initiation and development. In our previous research, the fuzzless gene in DPL972 was identified as a single dominant gene and named GaFzl. In the present study, we fine mapped this gene using F2 and BC1 populations. By combining traditional map-based cloning and next-generation sequencing, we mapped GaFzl to a 70-kb region containing seven annotated genes. RNA-Sequencing and re-sequencing analysis narrowed these candidates to two differentially expressed genes, Cotton_A_11941 and Cotton_A_11942. Sequence alignment uncovered no variation in coding or promoter regions of Cotton_A_11942 between DPL971 and DPL972, whereas two single-base mutations in the promoter region and a TTG insertion in the coding region were detected in Cotton_A_11941 in DPL972. Cotton_A_11941 encoding a homologous gene of GIR1 (GLABRA2-interacting repressor) in Arabidopsis thaliana is thus the candidate gene most likely responsible for the fuzzless trait in DPL972. Our findings should lead to a better understanding of cotton fuzz formation, thereby accelerating marker-assisted selection during cotton breeding.

Genetic dissection of the fuzzless seed trait in Gossypium barbadense

Cotton fibres are single-celled trichomes arising from the epidermal cells of the seed coat and may be either long (lint) or very short (fuzz). The dominant fuzzless N1 of Gossypium hirsutum is a defective allele of the At-subgenome homoeolog of MYB25-like, but the genetic components underlying the recessive fuzzless trait from G. barbadense (Gb) are unknown. We have identified five genetic loci, including a major contributing locus containing MYB25-like_Dt, associated with Gb fuzzless seeds based on genotyping of fuzzy and fuzzless near isogenic lines (NILs) from an interspecies cross (G. barbadense × G. hirsutum). At 3 d post-anthesis when fuzz fibres are initiating, expression of MYB25-like_Dt was significantly lower in fuzzless NILs than in fuzzy seeded NILs, while higher MYB25-like_Dt expression was associated with more seed fuzz across different cotton genotypes. Phenotypic and genotypic analysis of MYB25-like homoeoalleles in cottons showing different fibre phenotypes and their crossing progeny indicated that both MYB25-like_At and MYB25-like_Dt are associated with lint development, and that fuzz development is mainly determined by the expression level of MYB25-like_Dt at ~3 d post-anthesis. Expression of Gb fuzzless seeds depends on genetic background and interactions amongst the multiple loci identified. MYB25-like_Dt is one of the best candidates for N2.

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.

Dissecting Genetic Network of Fruit Branch Traits in Upland Cotton by Association Mapping Using SSR Markers

Genetic architecture of branch traits has large influences on the morphological structure, photosynthetic capacity, planting density, and yield of Upland cotton (Gossypium hirsutum L.). This research aims to reveal the genetic effects of six branch traits, including bottom fruit branch node number (BFBNN), bottom fruit branch length (BFBL), middle fruit branch node number (MFBNN), middle fruit branch length (MFBL), upper fruit branch node number (UFBNN), and upper fruit branch length (UFBL). Association mapping was conducted for these traits of 39 lines and their 178 F₁ hybrids in three environments. There were 20 highly significant Quantitative Trait SSRs (QTSs) detected by mixed linear model approach analyzing a full genetic model with genetic effects of additive, dominance, epistasis and their environment interaction. The phenotypic variation explained by genetic effects ranged from 32.64 ~ 91.61%, suggesting these branch traits largely influenced by genetic factors.

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.

Integrative transcriptome, proteome, phosphoproteome and genetic mapping reveals new aspects in a fiberless mutant of cotton

To investigate the molecular mechanisms of fiber initiation in cotton (Gossypium spp.), an integrated approach combining transcriptome, iTRAQ-based proteome and genetic mapping was taken to compare the ovules of the Xuzhou 142 wild type (WT) with its fuzzless-lintless (fl) mutant at −3 and 0 day post-anthesis. A total of 1,953 mRNAs, 187 proteins, and 131 phosphoproteins were differentially expressed (DE) between WT and fl, and the levels of transcripts and their encoded proteins and phosphoproteins were highly congruent. A functional analysis suggested that the abundance of proteins were mainly involved in amino sugar, nucleotide sugar and fatty acid metabolism, one carbon pool for folate metabolism and flavonoid biosynthesis. qRT-PCR, Western blotting, and enzymatic assays were performed to confirm the regulation of these transcripts and proteins. A molecular mapping located the lintless gene li3 in the fl mutant on chromosome 26 for the first time. A further in-silico physical mapping of DE genes with sequence variations between fl and WT identified one and four candidate genes in the li3 and n2 regions, respectively. Taken together, the transcript abundance, phosphorylation status of proteins at the fiber initiation stage and candidate genes have provided insights into regulatory processes underlying cotton fiber initiation.

RNA-Seq-Mediated Transcriptome Analysis of a Fiberless Mutant Cotton and Its Possible Origin Based on SNP Markers

As the longest known single-celled trichomes, cotton (Gossypium L.) fibers constitute a classic model system to investigate cell initiation and elongation. In this study, we used a high-throughput transcriptome sequencing technology to identify fiber-initiation-related single nucleotide polymorphism (SNP) markers and differentially expressed genes (DEGs) between the wild-type (WT) Upland cotton (G. hirsutum) Xuzhou 142 and its natural fuzzless-lintless mutant Xuzhou 142 fl. Approximately 700 million high-quality cDNA reads representing over 58 Gb of sequences were obtained, resulting in the identification of 28,610 SNPs--of which 17,479 were novel--from 13,960 expressed genes. Of these SNPs, 50% of SNPs in fl were identical to those of G. barbadense, which suggests the likely origin of the fl mutant from an interspecific hybridization between Xuzhou 142 and an unknown G. barbadense genotype. Of all detected SNPs, 15,555, 12,750, and 305 were classified as non-synonymous, synonymous, and pre-terminated ones, respectively. Moreover, 1,352 insertion/deletion polymorphisms (InDels) were also detected. A total of 865 DEGs were identified between the WT and fl in ovules at -3 and 0 days post-anthesis, with 302 candidate SNPs selected from these DEGs for validation by a high-resolution melting analysis and Sanger sequencing in seven cotton genotypes. The number of genotypic pairwise polymorphisms varied from 43 to 302, indicating that the identified SNPs are reliable. These SNPs should serve as good resources for breeding and genetic studies in cotton.

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.

Construction of BAC contig maps of homoeologous chromosomes A12 and D12 of Gossypium hirsutum L. acc. TM-1

Background

The Gossypium hirsutum homoeologous chromosome 12 encodes important genes that contribute to fiber fuzz, lethality, gland development and male sterility. In this study a physical map of the cotton TM-1 chromosome 12 was constructed. A number of large-insert cotton genome libraries are available, and genome-wide physical mapping using large insert segments combined with bacterial cloning is a thriving area of genome research. However, sequencing of the cotton genome is difficult due to sequence repeats and homoeologous regions. In order to effectively distinguish the homologous segments, a new method for adjusting the parameters of the FPC software was applied for contig map construction.

Results

All available markers on chromosomes A12 and D12 were used to screen the TM-1 BAC library by PCR. A total of 775 clones (387 for A12, 388 for D12) were obtained using Hind III fingerprinting and used for construction of the contig map. Seven pairs of SSR markers located on A12 and D12 were chosen for contig analysis. Following optimization of the tolerance (10) and cutoff (1e-12) parameters, combining all clones from A12 and D12 produced two separate contigs.

Conclusions

The BAC contig map of chromosomes A12 and D12 was constructed and FPC software parameters were optimized for analysis. The resulting approach is a powerful platform for genome-wide and evolutionary research on cotton.

The Hairless Stem Phenotype of Cotton (Gossypium barbadense) Is Linked to a Copia-Like Retrotransposon Insertion in a Homeodomain-Leucine Zipper Gene (HD1)

Cotton (Gossypium) stem trichomes are mostly single cells that arise from stem epidermal cells. In this study, a homeodomain-leucine zipper gene (HD1) was found to cosegregate with the dominant trichome locus previously designated as T1 and mapped to chromosome 6. Characterization of HD1 orthologs revealed that the absence of stem trichomes in modern Gossypium barbadense varieties is linked to a large retrotransposon insertion in the ninth exon, 2565 bp downstream from the initial codon in the At subgenome HD1 gene (At-GbHD1). In both the At and Dt subgenomes, reduced transcription of GbHD1 genes is caused by this insertion. The disruption of At-HD1 further affects the expression of downstream GbMYB25 and GbHOX3 genes. Analyses of primitive cultivated accessions identified another retrotransposon insertion event in the sixth exon of At-GbHD1 that might predate the previously identified retrotransposon in modern varieties. Although both retrotransposon insertions results in similar phenotypic changes, the timing of these two retrotransposon insertion events fits well with our current understanding of the history of cotton speciation and dispersal. Taken together, the results of genetics mapping, gene expression and association analyses suggest that GbHD1 is an important component that controls stem trichome development and is a promising candidate gene for the T1 locus. The interspecific phenotypic difference in stem trichome traits also may be attributable to HD1 inactivation associated with retrotransposon insertion.

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.

Construction of a high-density genetic map and lint percentage and cottonseed nutrient trait QTL identification in upland cotton (Gossypium hirsutum L.)

Upland cotton plays a critical role not only in the textile industry, but also in the production of important secondary metabolites, such as oil and proteins. Construction of a high-density linkage map and identifying yield and seed trait quantitative trail loci (QTL) are prerequisites for molecular marker-assisted selective breeding projects. Here, we update a high-density upland cotton genetic map from recombinant inbred lines. A total of 25,313 SSR primer pairs were screened for polymorphism between Yumian 1 and T586, and 1712 SSR primer pairs were used to genotype the mapping population and construct a map. An additional 1166 loci have been added to our previously published map with 509 SSR markers. The updated genetic map spans a total recombinant length of 3338.2 cM and contains 1675 SSR loci and nine morphological markers, with an average interval of 1.98 cM between adjacent markers. Green lint (Lg) mapped on chromosome 15 in a previous report is mapped in an interval of 2.6 cM on chromosome 21. Based on the map and phenotypic data from multiple environments, 79 lint percentage and seed nutrient trait QTL are detected. These include 8 lint percentage, 13 crude protein, 15 crude oil, 8 linoleic, 10 oleic, 13 palmitic, and 12 stearic acid content QTL. They explain 3.5-62.7 % of the phenotypic variation observed. Four morphological markers identified have a major impact on lint percentage and cottonseed nutrients traits. In this study, our genetic map provides new sights into the tetraploid cotton genome. Furthermore, the stable QTL and morphological markers could be used for fine-mapping and map-based cloning.

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.

Comparative proteomic and biochemical analyses reveal different molecular events occurring in the process of fiber initiation between wild-type allotetraploid cotton and its fuzzless-lintless mutant

To explore lint fiber initiation-related proteins in allotetraploid cotton (Gossypium hirsutum L.), a comparative proteomic analysis was performed between wild-type cotton (Xu-142) and its fuzzless-lintless mutant (Xu-142-fl) at five developmental time points for lint fiber initiation from -3 to +3 days post-anthesis (dpa). Using two-dimensional gel electrophoresis (2-DE) combined with mass spectrometry (MS) analyses, 91 differentially accumulated protein (DAP) species that are related to fiber initiation were successfully identified, of which 58 preferentially accumulated in the wild-type and 33 species in the fl mutant. These DAPs are involved in various cellular and metabolic processes, mainly including important energy/carbohydrate metabolism, redox homeostasis, amino acid and fatty acid biosynthesis, protein quality control, cytoskeleton dynamics, and anthocyanidin metabolism. Further physiological and biochemical experiments revealed dynamic changes in the carbohydrate flux and H₂O₂ levels in the cotton fiber initiation process. Compared with those in the fl mutant, the contents of glucose and fructose in wild-type ovules sharply increased after anthesis with a relatively higher rate of amino acid biosynthesis. The relative sugar starvation and lower rate of amino acid biosynthesis in the fl mutant ovules may impede the carbohydrate/energy supply and cell wall synthesis, which is consistent with the proteomic results. However, the H₂O₂ burst was only observed in the wild-type ovules on the day of anthesis. Cotton boll injection experiments in combination with electron microscope observation collectively indicated that H₂O₂ burst, which is negatively regulated by ascorbate peroxidases (APx), plays an important role in the fiber initiation process. Taken together, our study demonstrates a putative network of DAP species related to fiber initiation in cotton ovules and provides a foundation for future studies on the specific functions of these proteins in fiber development.

Next generation genetic mapping of the Ligon-lintless-2 (Li₂) locus in upland cotton (Gossypium hirsutum L.)

Mapping-by-sequencing and novel subgenome-specific SNP markers were used to fine map the Ligon-lintless 2 ( Li 2 ) short-fiber gene in tetraploid cotton. These methodologies will accelerate gene identification in polyploid species. Next generation sequencing offers new ways to identify the genetic mechanisms that underlie mutant phenotypes. The release of a reference diploid Gossypium raimondii (D5) genome and bioinformatics tools to sort tetraploid reads into subgenomes has brought cotton genetic mapping into the genomics era. We used multiple high-throughput sequencing approaches to identify the relevant region of reference sequence and identify single nucleotide polymorphisms (SNPs) near the short-fiber mutant Ligon-lintless 2 (Li 2) gene locus. First, we performed RNAseq on 8-day post-anthesis (DPA) fiber cells from the Li 2 mutant and its wild type near isogenic line (NIL) Gossypium hirsutum cv. DP5690. We aligned sequence reads to the D5 genome, sorted the reads into A and D subgenomes with PolyCat and called SNPs with InterSNP. We then identified SNPs that would result in non-synonymous substitutions to amino acid sequences of annotated genes. This step allowed us to identify a 1-Mb region with 24 non-synonymous SNPs, representing the introgressed region that differentiates Li 2 from its NIL. Next, we sequenced total DNA from pools of F2 plants, using a super bulked segregant analysis sequencing (sBSAseq) approach. The sBSAseq predicted 82 non-synonymous SNPs among 3,494 SNPs in a 3-Mb region that includes the region identified by RNAseq. We designed subgenome-specific SNP markers and tested them in an F2 population of 1,733 individuals to construct a genetic map. Our resulting genetic interval contains only one gene, an aquaporin, which is highly expressed in wild-type fibers and is significantly under-expressed in elongating Li 2 fiber cells.

Identification of centromeric regions on the linkage map of cotton using centromere-related repeats

Centromere usually contains high-copy-number retrotransposons and satellite repeats, which are difficult to map, clone and sequence. Currently, very little is known about the centromere in cotton. Here, we sequenced a bacterial artificial chromosome (BAC) mapping to the centromeric region and predicted four long-terminal-repeat (LTR) retrotransposons. They were located in the heterochromatic centromeric regions of all 52 pachytene chromosomes in Gossypium hirsutum. Fiber-FISH mapping revealed that these retrotransposons span an area of at least 1.8Mb in the centromeric region. Comparative analysis showed that these retrotransposons generated similar, strong fluorescent signals in the D progenitor Gossypium raimondii but not in the A progenitor Gossypium herbaceum, suggesting that the centromere sequence of tetraploid cotton might be derived from the D progenitor. Centromeric regions were anchored on 13 chromosomes of D-genome sequence. Characterization of these centromere-related repeats and regions will enhance cotton centromere mapping, sequencing and evolutionary studies.

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.

The Li2 mutation results in reduced subgenome expression bias in elongating fibers of allotetraploid cotton (Gossypium hirsutum L.)

Next generation sequencing (RNA-seq) technology was used to evaluate the effects of the Ligon lintless-2 (Li₂) short fiber mutation on transcriptomes of both subgenomes of allotetraploid cotton (Gossypium hirsutum L.) as compared to its near-isogenic wild type. Sequencing was performed on 4 libraries from developing fibers of Li₂ mutant and wild type near-isogenic lines at the peak of elongation followed by mapping and PolyCat categorization of RNA-seq data to the reference D₅ genome (G. raimondii) for homeologous gene expression analysis. The majority of homeologous genes, 83.6% according to the reference genome, were expressed during fiber elongation. Our results revealed: 1) approximately two times more genes were induced in the A_T subgenome comparing to the D_T subgenome in wild type and mutant fiber; 2) the subgenome expression bias was significantly reduced in the Li₂ fiber transcriptome; 3) Li₂ had a significantly greater effect on the D_T than on the A_T subgenome. Transcriptional regulators and cell wall homeologous genes significantly affected by the Li₂ mutation were reviewed in detail. This is the first report to explore the effects of a single mutation on homeologous gene expression in allotetraploid cotton. These results provide deeper insights into the evolution of allotetraploid cotton gene expression and cotton fiber development.

Gene expression profile analysis of Ligon lintless-1 (Li1) mutant reveals important genes and pathways in cotton leaf and fiber development

Ligon lintless-1 (Li1) is a monogenic dominant mutant of Gossypium hirsutum (upland cotton) with a phenotype of impaired vegetative growth and short lint fibers. Despite years of research involving genetic mapping and gene expression profile analysis of Li1 mutant ovule tissues, the gene remains uncloned and the underlying pathway of cotton fiber elongation is still unclear. In this study, we report the whole genome-level deep-sequencing analysis of leaf tissues of the Li1 mutant. Differentially expressed genes in leaf tissues of mutant versus wild-type (WT) plants are identified, and the underlying pathways and potential genes that control leaf and fiber development are inferred. The results show that transcription factors AS2, YABBY5, and KANDI-like are significantly differentially expressed in mutant tissues compared with WT ones. Interestingly, several fiber development-related genes are found in the downregulated gene list of the mutant leaf transcriptome. These genes include heat shock protein family, cytoskeleton arrangement, cell wall synthesis, energy, H2O2 metabolism-related genes, and WRKY transcription factors. This finding suggests that the genes are involved in leaf morphology determination and fiber elongation. The expression data are also compared with the previously published microarray data of Li1 ovule tissues. Comparative analysis of the ovule transcriptomes of Li1 and WT reveals that a number of pathways important for fiber elongation are enriched in the downregulated gene list at different fiber development stages (0, 6, 9, 12, 15, 18dpa). Differentially expressed genes identified in both leaf and fiber samples are aligned with cotton whole genome sequences and combined with the genetic fine mapping results to identify a list of candidate genes for Li1.

A comprehensive meta QTL analysis for fiber quality, yield, yield related and morphological traits, drought tolerance, and disease resistance in tetraploid cotton

BACKGROUND

The study of quantitative trait loci (QTL) in cotton (Gossypium spp.) is focused on traits of agricultural significance. Previous studies have identified a plethora of QTL attributed to fiber quality, disease and pest resistance, branch number, seed quality and yield and yield related traits, drought tolerance, and morphological traits. However, results among these studies differed due to the use of different genetic populations, markers and marker densities, and testing environments. Since two previous meta-QTL analyses were performed on fiber traits, a number of papers on QTL mapping of fiber quality, yield traits, morphological traits, and disease resistance have been published. To obtain a better insight into the genome-wide distribution of QTL and to identify consistent QTL for marker assisted breeding in cotton, an updated comparative QTL analysis is needed.

RESULTS

In this study, a total of 1,223 QTL from 42 different QTL studies in Gossypium were surveyed and mapped using Biomercator V3 based on the Gossypium consensus map from the Cotton Marker Database. A meta-analysis was first performed using manual inference and confirmed by Biomercator V3 to identify possible QTL clusters and hotspots. QTL clusters are composed of QTL of various traits which are concentrated in a specific region on a chromosome, whereas hotspots are composed of only one trait type. QTL were not evenly distributed along the cotton genome and were concentrated in specific regions on each chromosome. QTL hotspots for fiber quality traits were found in the same regions as the clusters, indicating that clusters may also form hotspots.

CONCLUSIONS

Putative QTL clusters were identified via meta-analysis and will be useful for breeding programs and future studies involving Gossypium QTL. The presence of QTL clusters and hotspots indicates consensus regions across cultivated tetraploid Gossypium species, environments, and populations which contain large numbers of QTL, and in some cases multiple QTL associated with the same trait termed a hotspot. This study combines two previous meta-analysis studies and adds all other currently available QTL studies, making it the most comprehensive meta-analysis study in cotton to date.

A whole-genome DNA marker map for cotton based on the D-genome sequence of Gossypium raimondii L

We constructed a very-high-density, whole-genome marker map (WGMM) for cotton by using 18,597 DNA markers corresponding to 48,958 loci that were aligned to both a consensus genetic map and a reference genome sequence. The WGMM has a density of one locus per 15.6 kb, or an average of 1.3 loci per gene. The WGMM was anchored by the use of colinear markers to a detailed genetic map, providing recombinational information. Mapped markers occurred at relatively greater physical densities in distal chromosomal regions and lower physical densities in the central regions, with all 1 Mb bins having at least nine markers. Hotspots for quantitative trait loci and resistance gene analog clusters were aligned to the map and DNA markers identified for targeting of these regions of high practical importance. Based on the cotton D genome reference sequence, the locations of chromosome structural rearrangements plotted on the map facilitate its translation to other Gossypium genome types. The WGMM is a versatile genetic map for marker assisted breeding, fine mapping and cloning of genes and quantitative trait loci, developing new genetic markers and maps, genome-wide association mapping, and genome evolution studies.

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

Background

Cotton fiber length is very important to the quality of textiles. Understanding the genetics and physiology of cotton fiber elongation can provide valuable tools to the cotton industry by targeting genes or other molecules responsible for fiber elongation. Ligon Lintless-1 (Li₁) is a monogenic mutant in Upland cotton (Gossypium hirsutum) which exhibits an early cessation of fiber elongation resulting in very short fibers (< 6 mm) at maturity. This presents an excellent model system for studying the underlying molecular and cellular processes involved with cotton fiber elongation. Previous reports have characterized Li₁ at early cell wall elongation and during later secondary cell wall synthesis, however there has been very limited analysis of the transition period between these developmental time points.

Results

Physical and morphological measurements of the Li₁ mutant fibers were conducted, including measurement of the cellulose content during development. Affymetrix microarrays were used to analyze transcript profiles at the critical developmental time points of 3 days post anthesis (DPA), the late elongation stage of 12 DPA and the early secondary cell wall synthesis stage of 16 DPA. The results indicated severe disruption to key hormonal and other pathways related to fiber development, especially pertaining to the transition stage from elongation to secondary cell wall synthesis. Gene Ontology enrichment analysis identified several key pathways at the transition stage that exhibited altered regulation. Genes involved in ethylene biosynthesis and primary cell wall rearrangement were affected, and a primary cell wall-related cellulose synthase was transcriptionally repressed. Linkage mapping using a population of 2,553 F₂ individuals identified SSR markers associated with the Li₁ genetic locus on chromosome 22. Linkage mapping in combination with utilizing the diploid G. raimondii genome sequences permitted additional analysis of the region containing the Li₁ gene.

Conclusions

The early termination of fiber elongation in the Li₁ mutant is likely controlled by an early upstream regulatory factor resulting in the altered regulation of hundreds of downstream genes. Several elongation-related genes that exhibited altered expression profiles in the Li₁ mutant were identified. Molecular markers closely associated with the Li₁ locus were developed. Results presented here will lay the foundation for further investigation of the genetic and molecular mechanisms of fiber elongation.

Discovery and identification of a novel Ligon lintless-like mutant (Lix) similar to the Ligon lintless (Li1) in allotetraploid cotton

Mutants are a powerful resource for studying gene structure, function, and evolution. In this present study, a novel Ligon lintless-like mutant (Lix), that has short fibers and deformed leaves and stems, was isolated from the progeny of transgenic cottons. The Lix mutant is similar in morphology to the Ligon lintless (Li1) mutant. Genetic analysis and molecular mapping were performed for the Lix and Li1 mutants. These two mutants are monogenic dominant mutants with distorted growth of vegetative and reproductive structures. Seedlings of the dominant homozygote Li 1 Li 1 genotype are lethal, while LixLix plants are viable but show no reproductive growth. Molecular tagging showed that the Lix gene is located on Chr. 04 in a 30.9-cM region spanned by NAU8376 and NAU3469. In a previous study, the Li 1 gene was mapped to Chr. 22, and Chr. 04 and Chr. 22 are homoelogous chromosomes in tetraploid cotton. So, we propose that Lix and Li1 mutants have similar mutated morphology, and Lix is mapped to a homoelogous chromosome carrying Li 1 . The identification and genetic mapping of Lix/Li 1 genes using mutants provides a foundation for isolating these genes. In turn, this will permit studies to elucidate the functional and evolutionary roles for these genes in cotton growth and development.

Integrated metabolomics and genomics analysis provides new insights into the fiber elongation process in Ligon lintless-2 mutant cotton (Gossypium hirsutum L.)

BACKGROUND

The length of cotton fiber is an important agronomic trait characteristic that directly affects the quality of yarn and fabric. The cotton (Gossypium hirsutum L.) fiber mutation, Ligon lintless-2, is controlled by a single dominant gene (Li(2)) and results in extremely shortened lint fibers on mature seeds with no visible pleiotropic effects on vegetative growth and development. The Li(2) mutant phenotype provides an ideal model system to study fiber elongation. To understand metabolic processes involved in cotton fiber elongation, changes in metabolites and transcripts in the Li(2) mutant fibers were compared to wild-type fibers during development.

RESULTS

Principal component analysis of metabolites from GC-MS data separated Li(2) mutant fiber samples from WT fiber samples at the WT elongation stage, indicating that the Li(2) mutation altered the metabolome of the mutant fibers. The observed alterations in the Li(2) metabolome included significant reductions in the levels of detected free sugars, sugar alcohols, sugar acids, and sugar phosphates. Biological processes associated with carbohydrate biosynthesis, cell wall loosening, and cytoskeleton were also down-regulated in Li(2) fibers. Gamma-aminobutyric acid, known as a signaling factor in many organisms, was significantly elevated in mutant fibers. Higher accumulation of 2-ketoglutarate, succinate, and malate suggested higher nitrate assimilation in the Li(2) line. Transcriptional activation of genes involved in nitrogen compound metabolism along with changes in the levels of nitrogen transport amino acids suggested re-direction of carbon flow into nitrogen metabolism in Li(2) mutant fibers.

CONCLUSIONS

This report provides the first comprehensive analysis of metabolite and transcript changes in response to the Li(2) mutation in elongating fibers. A number of factors associated with cell elongation found in this study will facilitate further research in understanding metabolic processes of cotton fiber elongation.

Comparative proteomic analysis reveals differentially expressed proteins correlated with fuzz fiber initiation in diploid cotton (Gossypium arboreum L.)

In this study, a comparative proteomic analysis was employed to identify fuzz fiber initiation-related proteins in wild-type diploid cotton (Gossypium arboreum L.) and its fuzzless mutant. Temporal changes in global proteomes were examined using 2-DE at five developmental time points for fuzz fiber initiation, and 71 differentially expressed protein species were identified by MS, 45 of which were preferentially accumulated in the wild-type. These proteins were assigned to several functional categories, mainly in cell response/signal transduction, redox homeostasis, protein metabolism and energy/carbohydrate metabolism. It was remarkable that more than ten key proteins with high-abundance were involved in gibberellic acid (GA) signaling and ROS scavenging, and increasing concentrations of active GAs and H2O2 were also detected approximately 5dpa in wild type ovules. Furthermore, in vivo GA and H2O2 treatments of ovules inside young bolls showed that these compounds can synergistically promote fuzz fiber initiation. Our findings not only described a dynamic protein network supporting fuzz initiation in diploid cotton fiber ovules, but also deepened our understanding of the molecular basis of cotton fiber initiation.

BIOLOGICAL SIGNIFICANCE

Our study reported the identification of differentially expressed proteins in wild-type diploid cotton (G. arboreum L.) and its fuzzless mutant by comparative proteomic approach. In total, 71 protein species related to fuzz initiation were identified by MS. These proteins were assigned to several functional categories, mainly in energy/carbohydrate metabolism, protein metabolism, signal transduction, redox homeostasis etc. Importantly, a number of key proteins were found to be associated with GA signaling and ROS scavenging. In consistence with these findings, we detected the increase of GAs and H2O2 concentrations during fiber initiation, and our in vivo ovule experiments with GA and H2O2 injection and following microscopy observation of fuzz fiber initiation supported promoting effects of GA and H2O2 on cotton fiber initiation. These findings depicted a dynamic protein network supporting cotton fiber initiation in diploid cotton ovules. Our study is of major significance for understanding the molecular mechanisms controlling fuzz initiation and also provides a solid basis for further functional research of single nodes of this network in relation to cotton fiber initiation.

Enhancing sucrose synthase activity results in increased levels of starch and ADP-glucose in maize (Zea mays L.) seed endosperms

Sucrose synthase (SuSy) is a highly regulated cytosolic enzyme that catalyzes the conversion of sucrose and a nucleoside diphosphate into the corresponding nucleoside diphosphate glucose and fructose. In cereal endosperms, it is widely assumed that the stepwise reactions of SuSy, UDPglucose pyrophosphorylase and ADPglucose (ADPG) pyrophosphorylase (AGP) take place in the cytosol to convert sucrose into ADPG necessary for starch biosynthesis, although it has also been suggested that SuSy may participate in the direct conversion of sucrose into ADPG. In this study, the levels of the major primary carbon metabolites, and the activities of starch metabolism-related enzymes were assessed in endosperms of transgenic maize plants ectopically expressing StSUS4, which encodes a potato SuSy isoform. A total of 29 fertile lines transformed with StSUS4 were obtained, five of them containing a single copy of the transgene that was still functional after five generations. The number of seeds per ear of the five transgenic lines containing a single StSUS4 copy was comparable with that of wild-type (WT) control seeds. However, transgenic seeds accumulated 10-15% more starch at the mature stage, and contained a higher amylose/amylopectin balance than WT seeds. Endosperms of developing StSUS4-expressing seeds exhibited a significant increase in SuSy activity, and in starch and ADPG contents when compared with WT endosperms. No significant changes could be detected in the transgenic seeds in the content of soluble sugars, and in activities of starch metabolism-related enzymes when compared with WT seeds. A suggested metabolic model is presented wherein both AGP and SuSy are involved in the production of ADPG linked to starch biosynthesis in maize endosperm cells.

No evidence for the occurrence of substrate inhibition of Arabidopsis thaliana sucrose synthase-1 (AtSUS1) by fructose and UDP-glucose

Sucrose synthase (SuSy) catalyzes the reversible conversion of sucrose and NDP into the corresponding nucleotide-sugars and fructose. The Arabidopsis genome possesses six SUS genes (AtSUS1-6) that code for proteins with SuSy activity. As a first step to investigate optimum fructose and UDP-glucose (UDPG) concentrations necessary to measure maximum sucrose-producing SuSy activity in crude extracts of Arabidopsis, in this work we performed kinetic analyses of recombinant AtSUS1 in two steps: (1) SuSy reaction at pH 7.5, and (2) chromatographic measurement of sucrose produced in step 1. These analyses revealed a typical Michaelis-Menten behavior with respect to both UDPG and fructose, with Km values of 50 μM and 25 mM, respectively. Unlike earlier studies showing the occurrence of substrate inhibition of UDP-producing AtSUS1 by fructose and UDP-glucose, these analyses also revealed no substrate inhibition of AtSUS1 at any UDPG and fructose concentration. By including 200 mM fructose and 1 mM UDPG in the SuSy reaction assay mixture, we found that sucrose-producing SuSy activity in leaves and stems of Arabidopsis were exceedingly higher than previously reported activities. Furthermore, we found that SuSy activities in organs of the sus1/sus2/sus3/sus4 mutant were ca. 80-90% of those found in WT plants.

Analyses of the sucrose synthase gene family in cotton: structure, phylogeny and expression patterns

BACKGROUND

In plants, sucrose synthase (Sus) is widely considered as a key enzyme involved in sucrose metabolism. Several paralogous genes encoding different isozymes of Sus have been identified and characterized in multiple plant genomes, while limited information of Sus genes is available to date for cotton.

RESULTS

Here, we report the molecular cloning, structural organization, phylogenetic evolution and expression profiles of seven Sus genes (GaSus1 to 7) identified from diploid fiber cotton (Gossypium arboreum). Comparisons between cDNA and genomic sequences revealed that the cotton GaSus genes were interrupted by multiple introns. Comparative screening of introns in homologous genes demonstrated that the number and position of Sus introns are highly conserved among Sus genes in cotton and other more distantly related plant species. Phylogenetic analysis showed that GaSus1, GaSus2, GaSus3, GaSus4 and GaSus5 could be clustered together into a dicot Sus group, while GaSus6 and GaSus7 were separated evenly into other two groups, with members from both dicot and monocot species. Expression profiles analyses of the seven Sus genes indicated that except GaSus2, of which the transcripts was undetectable in all tissues examined, and GaSus7, which was only expressed in stem and petal, the other five paralogues were differentially expressed in a wide ranges of tissues, and showed development-dependent expression profiles in cotton fiber cells.

CONCLUSIONS

This is a comprehensive study of the Sus gene family in cotton plant. The results presented in this work provide new insights into the evolutionary conservation and sub-functional divergence of the cotton Sus gene family in response to cotton fiber growth and development.