Naturally colored cotton for wearable applications

Naturally colored cotton (NCC) offers an environmentally friendly fiber for textile applications. Processing white cotton fiber into textiles requires extensive energy, water, and chemicals, whereas processing of NCC skips the most polluting activity, scouring-bleaching and dyeing; therefore, NCC provides an avenue to minimize the harmful impacts of textile production. NCC varieties are suitable for organic agriculture since they are naturally insect and disease-resistant, salt and drought-tolerant. Various fiber shades, ranging from light green to tan and brown, are available in the cultivated NCC (Gossypium hirsutum L.) species. The pigments responsible for the color of brown cotton fiber are proanthocyanidins or their derivatives synthesized by the flavonoid pathway. Due to pigments, the NCC has excellent ultraviolet protection properties. Some brown cotton varieties exhibited superior thermal resistance of fiber that can be used to make fabrics with enhanced flame retardancy. Here, we review molecular mechanisms involved in the pigment production of brown cotton and challenges in breeding NCC varieties with a wide range of colors but without penalty in fiber quality. Also, we discuss opportunities for NCC with flame-retarding properties in textile applications.

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

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

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

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

Mapping-by-sequencing the locus of EMS-induced mutation responsible for tufted-fuzzless seed phenotype in cotton

Cotton fiber mutants are valuable resources for studying functions of altered genes and their roles in fiber development. The n₄^t is a recessive tufted-fuzzless seed mutant created through chemical mutagenesis with ethyl methanesulfonate. Genetic analysis indicated that the tufted-fuzzless phenotype is controlled by a single recessive locus. In this study, we developed an F₂ population of 602 progeny plants and sequenced the genomes of the parents and two DNA bulks from F₂ progenies showing the mutant phenotype. We identified DNA sequence variants between the tufted-fuzzless mutant and wild type by aligning the sequence reads to the reference TM-1 genome and designed subgenome-specific SNP markers. We mapped the n₄^t locus on chromosome D04 within a genomic interval of about 411 kb. In this region, seven genes showed significant differential expression between the tufted-fuzzless mutant and wild type. Possible candidate genes are discussed in this study. The utilization of the n₄^t mutant along with other fiber mutants will facilitate our understanding of the molecular mechanisms of cotton fiber cell growth and development.

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.

Genetic and transcriptomic dissection of the fiber length trait from a cotton (Gossypium hirsutum L.) MAGIC population

BACKGROUND

Improving cotton fiber length without reducing yield is one of the major goals of cotton breeding. However, genetic improvement of cotton fiber length by breeding has been a challenge due to the narrow genetic diversity of modern cotton cultivars and negative correlations between fiber quality and yield traits. A multi-parent advanced generation inter-cross (MAGIC) population developed through random mating provides an excellent genetic resource that allows quantitative trait loci (QTL) and causal genes to be identified.

RESULTS

An Upland cotton MAGIC population, consisting of 550 recombinant inbred lines (RILs) derived from eleven different cultivars, was used to identify fiber length QTLs and potential genes that contribute to longer fibers. A genome wide association study (GWAS) identified a cluster of single nucleotide polymorphisms (SNPs) on chromosome (Chr.) D11 that is significantly associated with fiber length. Further evaluation of the Chr. D11 genomic region among lines of the MAGIC population detected that 90% of RILs have a D11 haplotype similar to the reference TM-1 genome (D11-ref), whereas 10% of RILs inherited an alternative haplotype from one of the parents (D11-alt). The average length of fibers of D11-alt RILs was significantly shorter compared to D11-ref RILs, suggesting that alleles in the D11-alt haplotype contributed to the inferior fiber quality. RNAseq analysis of the longest and shortest fiber length RILs from D11-ref and D11-alt populations identified 949 significantly differentially expressed genes (DEGs). Gene set enrichment analysis revealed that different functional categories of genes were over-represented during fiber elongation between the four selected RILs. We found 12 genes possessing non-synonymous SNPs (nsSNPs) significantly associated with the fiber length, and three that were highly significant and were clustered at D11:24-Mb, including D11G1928, D11G1929 and D11G1931.

CONCLUSION

The results of this study provide insights into molecular aspects of genetic variation in fiber length and suggests candidate genes for genetic manipulation for cotton improvement.

Genetic and transcriptomic dissection of the fiber length trait from a cotton (Gossypium hirsutum L.) MAGIC population

Background

Improving cotton fiber length without reducing yield is one of the major goals of cotton breeding. However, genetic improvement of cotton fiber length by breeding has been a challenge due to the narrow genetic diversity of modern cotton cultivars and negative correlations between fiber quality and yield traits. A multi-parent advanced generation inter-cross (MAGIC) population developed through random mating provides an excellent genetic resource that allows quantitative trait loci (QTL) and causal genes to be identified.

Results

An Upland cotton MAGIC population, consisting of 550 recombinant inbred lines (RILs) derived from eleven different cultivars, was used to identify fiber length QTLs and potential genes that contribute to longer fibers. A genome wide association study (GWAS) identified a cluster of single nucleotide polymorphisms (SNPs) on chromosome (Chr.) D11 that is significantly associated with fiber length. Further evaluation of the Chr. D11 genomic region among lines of the MAGIC population detected that 90% of RILs have a D11 haplotype similar to the reference TM-1 genome (D11-ref), whereas 10% of RILs inherited an alternative haplotype from one of the parents (D11-alt). The average length of fibers of D11-alt RILs was significantly shorter compared to D11-ref RILs, suggesting that alleles in the D11-alt haplotype contributed to the inferior fiber quality. RNAseq analysis of the longest and shortest fiber length RILs from D11-ref and D11-alt populations identified 949 significantly differentially expressed genes (DEGs). Gene set enrichment analysis revealed that different functional categories of genes were over-represented during fiber elongation between the four selected RILs. We found 12 genes possessing non-synonymous SNPs (nsSNPs) significantly associated with the fiber length, and three that were highly significant and were clustered at D11:24-Mb, including D11G1928, D11G1929 and D11G1931.

Conclusion

The results of this study provide insights into molecular aspects of genetic variation in fiber length and suggests candidate genes for genetic manipulation for cotton improvement.

Electronic supplementary material

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

The P450 gene CYP749A16 is required for tolerance to the sulfonylurea herbicide trifloxysulfuron sodium in cotton (Gossypium hirsutum L.)

BACKGROUND

Weed management is critical to global crop production and is complicated by rapidly evolving herbicide resistance in weeds. New sources of herbicide resistance are needed for crop plants so that applied herbicides can be rotated or combined to thwart the evolution of resistant weeds. The diverse family of cytochrome P450 proteins has been suggested to be a source of detoxifying herbicide metabolism in both weed and crop plants, and greater understanding of these genes will offer avenues for crop improvement and novel weed management practices.

RESULTS

Here, we report the identification of CYP749A16 (Gh_D10G1401) which is responsible for the natural tolerance exhibited by most cotton, Gossypium hirsutum L., cultivars to the herbicide trifloxysulfuron sodium (TFS, CGA 362622, commercial formulation Envoke). A 1-bp frameshift insertion in the third exon of CYP749A16 results in the loss of tolerance to TFS. The DNA marker designed from this insertion perfectly co-segregated with the phenotype in 2145 F₂ progeny of a cross between the sensitive cultivar Paymaster HS26 and tolerant cultivar Stoneville 474, and in 550 recombinant inbred lines of a multi-parent advanced generation inter-cross population. Marker analysis of 382 additional cotton cultivars identified twelve cultivars containing the 1-bp frameshift insertion. The marker genotypes matched perfectly with phenotypes in 188 plants from the selected twelve cultivars. Virus-induced gene silencing of CYP749A16 generated sensitivity in the tolerant cotton cultivar Stoneville 474.

CONCLUSIONS

CYP749A16 located on chromosome D10 is required for TFS herbicide tolerance in cotton. This finding should add to the repertoire of tools available to farmers and breeders for the advancement of agricultural productivity.

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₂.

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.

The GhTT2_A07 gene is linked to the brown colour and natural flame retardancy phenotypes of Lc1 cotton (Gossypium hirsutum L.) fibres

Some naturally coloured brown cotton fibres from accessions of Gossypium hirsutum L. can be used to make textiles with enhanced flame retardancy (FR). Several independent brown fibre loci have been identified and mapped to chromosomes, but the underlying genes have not yet been identified, and the mechanism of lint fibre FR is not yet fully understood. In this study, we show that both the brown colour and enhanced FR of the Lc1 lint colour locus are linked to a 1.4Mb inversion on chromosome A07 that is immediately upstream of a gene with similarity to Arabidopsis TRANSPARENT TESTA 2 (TT2). As a result of the alternative upstream sequence, the transcription factor GhTT2_A07 is highly up-regulated in developing fibres. In turn, genes in the phenylpropanoid metabolic pathway are activated, leading to biosynthesis of proanthocyanidins and accumulation of inorganic elements. We show that enhanced FR and anthocyanin precursors appear in developing brown fibres well before the brown colour is detectible, demonstrating for the first time that the polymerized proanthocyanidins that constitute the brown colour are not the source of enhanced FR. Identifying the particular colourless metabolite that provides Lc1 cotton with enhanced FR could help minimize the use of synthetic chemical flame retardant additives in textiles.

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.

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.

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.

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.

Independent replication of mitochondrial genes supports the transcriptional program in developing fiber cells of cotton (Gossypium hirsutum L.)

The mitochondrial genomes of flowering plants exist both as a "master circle" chromosome and as numerous subgenomic sublimons that are generated by intramolecular recombination. Differential stability or replication of these sublimons allows individual mitochondrial gene copy numbers to vary independently between different cell types and developmental stages. Our objective was to determine the relationship between mitochondrial gene copy number and transcript abundance in the elongating fiber cells of Upland cotton (Gossypium hirsutum L.). We compared RNA and DNA from cotton fiber cells at five developmental time points from early elongation through secondary cell wall thickening from the Ligon-lintless 2 (Li2) short fiber mutant and its wild type near isogenic line (NIL) DP5690. Mitochondrial gene copy number decreased from 3 to 8-DPA in the developing cotton fiber cells while transcript levels remained low. As secondary cell wall biosynthesis began in developing fibers, the expression levels and copy numbers of mitochondrial genes involved in energy production and respiration were up-regulated in wild type cotton DP5690. However, the short fiber mutant Li2, failed to increase expression of these genes, which include three subunits of ATP synthase, atp1, atp8 and atp9 and two cytochrome genes cox1 and cob. At the same time, Li2 failed to increase the copy numbers of these highly expressed genes. Surprisingly, we found that when mitochondrial genes were highly transcribed, they also had very high copy numbers. This observation suggests that in developing cotton fibers, increased mitochondrial sublimon replication may support increases in gene transcription.

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.

A transcript profiling approach reveals an abscisic acid-specific glycosyltransferase (UGT73C14) induced in developing fiber of Ligon lintless-2 mutant of cotton (Gossypium hirsutum L.)

Ligon lintless-2, a monogenic dominant cotton (Gossypium hirsutum L.) fiber mutation, causing extreme reduction in lint fiber length with no pleiotropic effects on vegetative growth, represents an excellent model system to study fiber elongation. A UDP-glycosyltransferase that was highly expressed in developing fibers of the mutant Ligon lintless-2 was isolated. The predicted amino acid sequence showed ~53% similarity with Arabidopsis UGT73C sub-family members and the UDP-glycosyltransferase was designated as UGT73C14. When expressed in Escherichia coli as a recombinant protein with a maltose binding protein tag, UGT73C14 displayed enzymatic activity toward ABA and utilized UDP-glucose and UDP-galactose as the sugar donors. The recombinant UGT73C14 converted natural occurring isoform (+)-cis, trans-ABA better than (+)-trans, trans-ABA and (-)-cis, trans-ABA. Transgenic Arabidopsis plants constitutively overexpressing UGT73C14 did not show phenotypic changes under standard growth conditions. However, the increased glycosylation of ABA resulted in phenotypic changes in post-germinative growth and seedling establishment, confirming in vivo activity of UGT73C14 for ABA. This suggests that the expression level of UGT73C14 is regulated by the observed elevated levels of ABA in developing fibers of the Li₂ mutant line and may be involved in the regulation of ABA homeostasis.

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.

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.

Characterization of the mannan synthase promoter from guar (Cyamopsis tetragonoloba)

Guar seed gum, consisting primarily of a high molecular weight galactomannan, is the most cost effective natural thickener, having broad applications in the food, cosmetics, paper, pharmaceutical and petroleum industries. The properties of the polymer can potentially be enhanced by genetic modification. Development of suitable endosperm-specific promoters for use in guar is desirable for metabolic engineering of the seed gum. A ~1.6 kb guar mannan synthase (MS) promoter region has been isolated. The MS promoter sequence was fused with the GUS reporter gene and overexpressed in the heterologous species alfalfa (Medicago sativa). The potential strength and specificity of the MS promoter was compared with those of the constitutive 35S promoter and the seed specific β-phaseolin promoter. Quantitative GUS assays revealed that the MS promoter directs GUS expression specifically in endosperm in transgenic alfalfa. Thus, the guar MS promoter could prove generally useful for directing endosperm-specific expression of transgenes in legume species.

Integrated metabolite and transcript profiling identify a biosynthetic mechanism for hispidol in Medicago truncatula cell cultures

Metabolic profiling of elicited barrel medic (Medicago truncatula) cell cultures using high-performance liquid chromatography coupled to photodiode and mass spectrometry detection revealed the accumulation of the aurone hispidol (6-hydroxy-2-[(4-hydroxyphenyl)methylidene]-1-benzofuran-3-one) as a major response to yeast elicitor. Parallel, large-scale transcriptome profiling indicated that three peroxidases, MtPRX1, MtPRX2, and MtPRX3, were coordinately induced with the accumulation of hispidol. MtPRX1 and MtPRX2 exhibited aurone synthase activity based upon in vitro substrate specificity and product profiles of recombinant proteins expressed in Escherichia coli. Hispidol possessed significant antifungal activity relative to other M. truncatula phenylpropanoids tested but has not been reported in this species before and was not found in differentiated roots in which high levels of the peroxidase transcripts accumulated. We propose that hispidol is formed in cell cultures by metabolic spillover when the pool of its precursor, isoliquiritigenin, builds up as a result of an imbalance between the upstream and downstream segments of the phenylpropanoid pathway, reflecting the plasticity of plant secondary metabolism. The results illustrate that integration of metabolomics and transcriptomics in genetically reprogrammed plant cell cultures is a powerful approach for the discovery of novel bioactive secondary metabolites and the mechanisms underlying their generation.

Metabolic and genetic perturbations accompany the modification of galactomannan in seeds of Medicago truncatula expressing mannan synthase from guar (Cyamopsis tetragonoloba L.)

Galactomannan gums are widely used in the food and oil industries, and there is considerable interest in applying biotechnological approaches to improve their physical properties. A mannan synthase from guar (Cyamopsis tetragonoloba) was expressed under the control of a bean beta-phaseolin promoter in transgenic Medicago truncatula. Although the expression of exogenous mannan synthase caused a slight decrease in galactomannan levels in Medicago, the molecular weight and viscosity of the polymer were significantly increased, although the mannose to galactose ratio and degree of polydispersity remained unchanged. At the same time, expression of about 2.8% of the genes was altered significantly in the seeds of transgenic Medicago lines analysed by Affymetrix genome chip, with a particularly striking induction of putative trehalose phosphate synthase genes. Mannan synthase expression also caused large alterations in the levels of a number of sugars and sugar alcohols, suggesting that over-expression of a processive glycosyltransferase perturbs the mechanisms of sugar sensing and/or homeostasis, possibly involving signalling via trehalose-6-phosphate.

Analysis of cDNA libraries from developing seeds of guar (Cyamopsis tetragonoloba (L.) Taub)

BACKGROUND

Guar, Cyamopsis tetragonoloba (L.) Taub, is a member of the Leguminosae (Fabaceae) family and is economically the most important of the four species in the genus. The endosperm of guar seed is a rich source of mucilage or gum, which forms a viscous gel in cold water, and is used as an emulsifier, thickener and stabilizer in a wide range of foods and industrial applications. Guar gum is a galactomannan, consisting of a linear (1-->4)-beta-linked D-mannan backbone with single-unit, (1-->6)-linked, alpha-D-galactopyranosyl side chains. To better understand regulation of guar seed development and galactomannan metabolism we created cDNA libraries and a resulting EST dataset from different developmental stages of guar seeds.

RESULTS

A database of 16,476 guar seed ESTs was constructed, with 8,163 and 8,313 ESTs derived from cDNA libraries I and II, respectively. Library I was constructed from seeds at an early developmental stage (15-25 days after flowering, DAF), and library II from seeds at 30-40 DAF. Quite different sets of genes were represented in these two libraries. Approximately 27% of the clones were not similar to known sequences, suggesting that these ESTs represent novel genes or may represent non-coding RNA. The high flux of energy into carbohydrate and storage protein synthesis in guar seeds was reflected by a high representation of genes annotated as involved in signal transduction, carbohydrate metabolism, chaperone and proteolytic processes, and translation and ribosome structure. Guar unigenes involved in galactomannan metabolism were identified. Among the seed storage proteins, the most abundant contig represented a conglutin accounting for 3.7% of the total ESTs from both libraries.

CONCLUSION

The present EST collection and its annotation provide a resource for understanding guar seed biology and galactomannan metabolism.

Methyl jasmonate and yeast elicitor induce differential transcriptional and metabolic re-programming in cell suspension cultures of the model legume Medicago truncatula

Exposure of cell suspension cultures of Medicago truncatula Gaerth. to methyl jasmonate (MeJA) resulted in up to 50-fold induction of transcripts encoding the key triterpene biosynthetic enzyme beta-amyrin synthase (betaAS; EC 5.4.99.-). Transcripts reached maximum levels at 24 h post-elicitation with 0.5 mM MeJA. The entry point enzymes into the phenylpropanoid and flavonoid pathways, L: -phenylalanine ammonia-lyase (PAL; EC 4.3.1.5) and chalcone synthase (CHS; EC 2.3.1.74), respectively, were not induced by MeJA. In contrast, exposure of cells to yeast elicitor (YE) resulted in up to 45- and 14-fold induction of PAL and CHS transcripts, respectively, at only 2 h post-elicitation. betaAS transcripts were weakly induced at 12 h after exposure to YE. Over 30 different triterpene saponins were identified in the cultures, many of which were strongly induced by MeJA, but not by YE. In contrast, cinnamic acids, benzoic acids and isoflavone-derived compounds accumulated following exposure of cultures to YE, but few changes in phenylpropanoid levels were observed in response to MeJA. DNA microarray analysis confirmed the strong differential transcriptional re-programming of the cell cultures for multiple genes in the phenylpropanoid and triterpene pathways in response to MeJA and YE, and indicated different responses of individual members of gene families. This work establishes Medicago cell cultures as an excellent model for future genomics approaches to understand the regulation of legume secondary metabolism.