Search

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Naoumkina, Hinchliffe and Thyssen.)

Naturally-colored brown cotton (NBC) fiber is an environmentally friendly raw source of fiber for textile applications. The fiber of some NBC cultivars exhibits flame-retardant properties, which can be used in textiles that require flame resistance. Proanthocyanidins or their derivatives are responsible for the brown pigment in NBC; however, how flame retardancy is related to pigmentation in NBC is poorly understood. To gain insight into brown pigment biosynthesis, we conducted comparative transcripts and metabolites profiling analysis of developing cotton fibers between the brown (MC-BL) and white (MC-WL) cotton near-isogenic lines (NILs), genetically different only in the Lc1 locus. In this study, mass spectrometry was used to detect metabolites in BL and WL developing fibers at 8, 12, 16, 20, 24, 36, and 40 days post anthesis (DPA) and mature fibers. Transcripts analysis was performed at two critical fiber developmental points, 8 DPA (fiber elongation) and 20 DPA (secondary cell wall deposition). We found 5836 (ESI MS positive mode) and 4541 (ESI MS negative mode) metabolites significantly different accumulated between BL and WL. Among them, 142 were known non-redundant metabolites, including organic acids, amino acids, and derivatives of the phenylpropanoid pathway. Transcript analysis determined 1691 (8 DPA) and 5073 (20 DPA) differentially expressed genes (DEGs) between BL and WL, with the majority of DEGs down-regulated at 20 DPA. Organic acids of the citric acid cycle were induced, while most of the detected amino acids were reduced in the MC-BL line. Both cis - and trans -stereoisomers of flavan-3-ols were detected in developing MC-WL and MC-BL fibers; however, the gallocatechin and catechin accumulated multiple times higher. Gas chromatography-mass spectrometry (GC-MS) analysis of fatty acids determined that palmitic acid long-chain alcohols were the main constituents of waxes of mature fibers. Energy-dispersive X-ray spectrometry (EDS) analysis of mature fibers revealed that potassium accumulated three times greater in MC-BL than in MC-WL mature fibers. This study provides novel insights into the biosynthesis of pigments and its association with flame retardancy in NBC fibers., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Hinchliffe, Naoumkina, Thyssen, Nam, Chang, McCarty and Jenkins.)

Background: Cotton (Gossypium sp.) has been cultivated for centuries for its spinnable fibers, but its seed oil also possesses untapped economic potential if, improvements could be made to its oleic acid content., Results: Previous studies, including those from our laboratory, identified pima accessions containing approximately doubled levels of seed oil oleic acid, compared to standard upland cottonseed oil. Here, the molecular properties of a fatty acid desaturase encoded by a mutant allele identified by genome sequencing in an earlier analysis were analyzed. The mutant sequence is predicted to encode a C-terminally truncated protein lacking nine residues, including a predicted endoplasmic reticulum membrane retrieval motif. We determined that the mutation was caused by a relatively recent movement of a Ty1/copia type retrotransposon that is not found associated with this desaturase gene in other sequenced cotton genomes. The mutant desaturase, along with its repaired isozyme and the wild-type A-subgenome homoeologous protein were expressed in transgenic yeast and stably transformed Arabidopsis plants. All full-length enzymes efficiently converted oleic acid to linoleic acid. The mutant desaturase protein produced only trace amounts of linoleic acid, and only when strongly overexpressed in yeast cells, indicating that the missing C-terminal amino acid residues are not strictly required for enzyme activity, yet are necessary for proper subcellular targeting to the endoplasmic reticulum membrane., Conclusion: These results provide the biochemical underpinning that links a genetic lesion present in a limited group of South American pima cotton accessions and their rare seed oil oleic acid traits. Markers developed to the mutant desaturase allele are currently being used in breeding programs designed to introduce this trait into agronomic upland cotton varieties., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

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., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

The chemical and physical properties of vegetable oils are largely dictated by the ratios of 4-6 common fatty acids contained within each oil. However, examples of plant species that accumulate from trace amounts to >90% of certain unusual fatty acids in seed triacylglycerols have been reported. Many of the general enzymatic reactions that drive both common and unusual fatty acid biosynthesis and accumulation in stored lipids are known, but which isozymes have evolved to specifically fill this role and how they coordinate in vivo is still poorly understood. Cotton (Gossypium sp.) is the very rare example of a commodity oilseed that produces biologically relevant amounts of unusual fatty acids in its seeds and other organs. In this case, unusual cyclopropyl fatty acids (named after the cyclopropane and cyclopropene moieties within the fatty acids) are found in membrane and storage glycerolipids (e.g. seed oils). Such fatty acids are useful in the synthesis of lubricants, coatings, and other types of valuable industrial feedstocks. To characterize the role of cotton acyltransferases in cyclopropyl fatty acid accumulation for bioengineering applications, we cloned and characterized type-1 and type-2 diacylglycerol acyltransferases from cotton and compared their biochemical properties to that of litchi (Litchi chinensis), another cyclopropyl fatty acid-producing plant. The results presented from transgenic microbes and plants indicate both cotton DGAT1 and DGAT2 isozymes efficiently utilize cyclopropyl fatty acid-containing substrates, which helps to alleviate biosynthetic bottlenecks and enhances total cyclopropyl fatty acid accumulation in the seed oil., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Published by Elsevier Masson SAS.)

Textiles made from cotton fibers are flammable and thus often include flame retardant additives for consumer safety. Transgressive segregation in multi-parent populations facilitates new combinations of alleles of genes and can result in traits that are superior to those of any of the parents. A screen of 257 recombinant inbred lines from a multi-parent advanced generation intercross (MAGIC) population for naturally enhance flame retardance (FR) was conducted. All eleven parents, like all conventional white fiber cotton cultivars produce flammable fabric. MAGIC recombinant inbred lines (RILs) that produced fibers with significantly lower heat release capacities (HRC) as measured by microscale combustion calorimetry (MCC) were identified and the stability of the phenotypes of the outliers were confirmed when the RILs were grown at an additional location. Of the textiles fabricated from the five superior RILs, four exhibited the novel characteristic of inherent flame resistance. When exposed to open flame by standard 45° incline flammability testing, these four fabrics self-extinguished. To determine the genetic architecture of this novel trait, linkage, epistatic and multi-locus genome wide association studies (GWAS) were conducted with 473k SNPs identified by whole genome sequencing (WGS). Transcriptomes of developing fiber cells from select RILs were sequenced (RNAseq). Together, these data provide insight into the genetic mechanism of the unexpected emergence of flame-resistant cotton by transgressive segregation in a breeding program. The incorporation of this trait into global cotton germplasm by breeding has the potential to greatly reduce the costs and impacts of flame-retardant chemicals., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

Introgression of superior fiber traits from Pima cotton (Gossypium barbadense, GB) into high yield Upland cotton (G. hirsutum) has been a breeding objective for many years in a few breeding programs in the world. However, progress has been very slow due to introgression barriers resulting from whole genome hybridization between the two species. To minimize such barriers, chromosome substitution lines (CS-B) from Pima cotton 3-79 in an Upland cotton cultivar TM-1 were developed. A multiparent advanced generation inter-cross (MAGIC) population consisting of 180 recombinant inbred lines (RILs) was subsequently made using the 18 CS-B lines and three Upland cotton cultivars as parents. In this research, we sequenced the whole genomes of the 21 parents and 180 RILs to examine the G. barbadense introgression. Of the 18 CS-B lines, 11 contained the target GB chromosome or chromosome segment, two contained more than two GB chromosomes, and five did not have the expected introgression. Residual introgression in non-target chromosomes was prevalent in all CS-B lines. A clear structure existed in the MAGIC population and the 180 RILs were distributed into three groups, i.e., high, moderate, and low GB introgression. Large blocks of GB chromosome introgression were still present in some RILs after five cycles of random-mating, an indication of recombination suppression or other unknown reasons present in the population. Identity by descent analysis revealed that the MAGIC RILs contained less introgression than expected. This research presents an insight on understanding the complex problems of introgression between cotton species., (© 2022. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

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., (Published by Oxford University Press on behalf of American Society of Plant Biologists 2022.)

Key Message: A major QTL conferring resistance to Fusarium wilt race 4 in a narrow region of chromosome D02 was identified in a MAGIC population of 550 RILs of Upland cotton. Numerous studies have been conducted to investigate the genetic basis of Fusarium wilt (FW, caused by Fusarium oxysporum f. sp. vasinfectum, FOV) resistance using bi-parental and association mapping populations in cotton. In this study, a multi-parent advanced generation inter-cross (MAGIC) population of 550 recombinant inbred lines (RILs), together with their 11 Upland cotton (Gossypium hirsutum) parents, was used to identify QTLs for FOV race 4 (FOV4) resistance. Among the parents, Acala Ultima, M-240 RNR, and Stoneville 474 were the most resistant, while Deltapine Acala 90, Coker 315, and Stoneville 825 were the most susceptible. Twenty-two MAGIC lines were consistently resistant to FOV4. Through a genome-wide association study (GWAS) based on 473,516 polymorphic SNPs, a major FOV4 resistance QTL within a narrow region on chromosomes D02 was detected, allowing identification of 14 candidate genes. Additionally, a meta-analysis of 133 published FW resistance QTLs showed a D subgenome and individual chromosome bias and no correlation between homeologous chromosome pairs. This study represents the first GWAS study using a largest genetic population and the most comprehensive meta-analysis for FW resistance in cotton. The results illustrated that 550 lines were not enough for high resolution mapping to pinpoint a candidate gene, and experimental errors in phenotyping cotton for FW resistance further compromised the accuracy and precision in QTL localization and identification of candidate genes. This study identified FOV4-resistant parents and MAGIC lines, and the first major QTL for FOV4 resistance in Upland cotton, providing useful information for developing FOV4-resistant cultivars and further genomic studies towards identification of causal genes for FOV4 resistance in cotton., (© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Cotton (Gossypium hirsutum L.) fiber is the most important resource of natural and renewable fiber for the textile industry. However, the understanding of genetic components and their genome-wide interactions controlling fiber quality remains fragmentary. Here, we sequenced a multiple-parent advanced-generation inter-cross (MAGIC) population, consisting of 550 individuals created by inter-crossing 11 founders, and established a mosaic genome map through tracing the origin of haplotypes that share identity-by-descent (IBD). We performed two complementary GWAS methods-SNP-based GWAS (sGWAS) and IBD-based haplotype GWAS (hGWAS). A total of 25 sQTLs and 14 hQTLs related to cotton fiber quality were identified, of which 26 were novel QTLs. Two major QTLs detected by both GWAS methods were responsible for fiber strength and length. The gene Ghir_D11G020400 (GhZF14) encoding the MATE efflux family protein was identified as a novel candidate gene for fiber length. Beyond the additive QTLs, we detected prevalent epistatic interactions that contributed to the genetics of fiber quality, pinpointing another layer for trait variance. This study provides new targets for future molecular design breeding of superior fiber quality., (© 2022. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

Thickness of cotton fiber, referred to as fiber maturity, is a key determinant of fiber quality, lint yield, and textile performance. The cotton immature fiber (im) mutant has been used to study fiber maturity since its fiber is thinner than the wild type near isogeneic line (NIL), Texas Marker-1 (TM-1). The im phenotype is caused by a single recessive mutation of a pentatricopeptide repeat (PPR) gene that reduces the activity of mitochondrial complex I and up-regulates stress responsive genes. However, the mechanisms altering the stress responses in im mutant are not well understood. Thus, we characterized growth and gas exchange in im and TM-1 under no stress and also investigated their stress responses by comparing gas exchange and transcriptomic profiles under high temperature. Phenotypic differences were detected between the NILs in non-fiber tissues although less pronounced than the variation in fibers. At near optimum temperature (28±3°C), im maintained the same photosynthetic performance as TM-1 by means of greater stomatal conductance. In contrast, under high temperature stress (>34°C), im leaves reduced photosynthesis by decreasing the stomatal conductance disproportionately more than TM-1. Transcriptomic analyses showed that the genes involved in heat stress responses were differentially expressed between the NIL leaves. These results indicate that the im mutant previously reported to have low activity of mitochondrial complex I displays increased thermosensitivity by impacting stomatal conductance. They also support a notion that mitochondrial complex I activity is required for maintenance of optimal photosynthetic performance and acclimation of plants to high temperature stress. These findings may be useful in the future efforts to understand how physiological mechanisms play a role in determining cotton fiber maturity and may influence stress responses in other crops., Competing Interests: This research was partially supported by Cotton Incorporated (#19-858). The funder provided supports in the form of research materials, but did not have any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. This does not alter our adherence to PLOS ONE policies on sharing data and materials. There are no patents, products in development or marketed products associated with this research to declare.

Cotton fiber mutants are valuable resources for studying functions of altered genes and their roles in fiber development. The n 4 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 2 population of 602 progeny plants and sequenced the genomes of the parents and two DNA bulks from F 2 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 4 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 4 t mutant along with other fiber mutants will facilitate our understanding of the molecular mechanisms of cotton fiber cell growth and development., (© 2021. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

Most commercially produced cotton cultivars have two types of fibers on the seed coat, short fuzz and long lint. Lint fiber is used in the textile industry, while fuzz is considered an undesirable trait. Both types of fibers are believed to be controlled by the same regulators; however, their mechanisms of actions are still obscure. Cotton fiber mutants provide an excellent system to study the genes that regulate fiber development. Here we described four uncharacterized and three previously reported cotton mutants with fuzzless seed phenotypes. To evaluate whether or not the genes previously associated with fuzzless seed phenotypes have mutations we sequenced whole genomic DNA of seven mutants and wild type varieties. We identified multiple polymorphic changes among the tested genes. Non-synonymous SNPs in the coding region of the MML3-A gene was common in the six mutant lines tested in this study, showing both dominant and recessive fuzzless phenotypes. We have mapped the locus of the causative mutation for one of the uncharacterized fuzzless lines using an F 2 population that originated from a cross between the dominant fuzzless mutant and a wild type. Further, we have clarified the current knowledge about the causative n 2 mutations by analyzing the sequence data and previously reported mapping data. The key genes and possible mechanisms of fiber differentiation are discussed in this study.

Cotton is grown in arid and semi-arid regions where abiotic stresses such as drought and salt are prevalent. There is a lack of studies that simultaneously address the genetic and genomic basis of tolerance to drought and salt stress. In this study, a multi-parent advanced generation inter-cross (MAGIC) population of 550 recombinant inbred lines (RILs) together with their 11 Upland cotton parents with a total of 473,516 polymorphic SNP markers was used to identify quantitative trait loci (QTL) for drought tolerance (DT) and salt tolerance (ST) at the seedling stage based on two replicated greenhouse tests. Transgressive segregation occurred in the MAGIC-RILs, indicating that tolerant and sensitive alleles recombined for tolerance to the abiotic stress during the intermating process for the population development. A total of 20 QTL were detected for DT including 13 and 7 QTL based on plant height (PH) and dry shoot weight (DSW), respectively; and 23 QTL were detected for ST including 12 and 11 QTL for PH and DSW, respectively. There were several chromosomes with QTL clusters for abiotic stress tolerance including four QTL on chromosome A13 and three QTL on A01 for DT, and four QTL on D08 and three QTL on A11 for ST. Nine QTL (21% of the 43 QTL) detected were in common between DT and ST, indicating a common genetic basis for DT and ST. The narrow chromosomal regions for most of the QTL detected in this study allowed identification of 53 candidate genes associated with responses to salt and drought stress and abiotic stimulus. The QTL identified for both DT and ST have significantly augmented the repertoire of QTL for abiotic stress tolerance that can be used for marker-assisted selection to develop cultivars with resilience to drought and/or salt and further genomic studies towards the identification of drought and salt tolerance genes in cotton.

The cotton fibre serves as a valuable experimental system to study cell wall synthesis in plants, but our understanding of the genetic regulation of this process during fibre development remains limited. We performed a genome-wide association study (GWAS) and identified 28 genetic loci associated with fibre quality in allotetraploid cotton. To investigate the regulatory roles of these loci, we sequenced fibre transcriptomes of 251 cotton accessions and identified 15 330 expression quantitative trait loci (eQTL). Analysis of local eQTL and GWAS data prioritised 13 likely causal genes for differential fibre quality in a transcriptome-wide association study (TWAS). Characterisation of distal eQTL revealed unequal genetic regulation patterns between two subgenomes, highlighted by an eQTL hotspot (Hot216) that established a genome-wide genetic network regulating the expression of 962 genes. The primary regulatory role of Hot216, and specifically the gene encoding a KIP-related protein, was found to be the transcriptional regulation of genes responsible for cell wall synthesis, which contributes to fibre length by modulating the developmental transition from rapid cell elongation to secondary cell wall synthesis. This study uncovered the genetic regulation of fibre-cell development and revealed the molecular basis of the temporal modulation of secondary cell wall synthesis during plant cell elongation., (© 2020 The Authors. New Phytologist © 2020 New Phytologist Trust.)

Key Message: 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 2 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 2 populations and F 3 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 2 and F 3 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.

Key Message: MAGIC population sequencing and virus-induced gene silencing identify Gh_D02G0276 as a novel root-knot nematode resistance gene on chromosome 14 in Upland cotton. The southern root-knot nematode [RKN; Meloidogyne incognita (Kofoid & White)] remains the primary yield-limiting biotic stress to Upland cotton (Gossypium hirsutum L.) throughout the southeastern USA. While useful genetic markers have been developed for two major RKN resistance loci on chromosomes 11 (A11) and 14 (D02), these markers are not completely effective because the causative genes have not been identified. Here, we sequenced 550 recombinant inbred lines (RILs) from a multi-parent advanced generation intercross (MAGIC) population to identify five RILs that had informative recombinations near the D02-RKN resistance locus. The RKN resistance phenotypes of these five RILs narrowed the D02-RKN locus to a 30-kb region with four candidate genes. We conducted virus-induced gene silencing (VIGS) on each of these genes and found that Gh_D02G0276 was required for suppression of RKN egg production conferred by the Chr. D02 resistance gene. The resistant lines all possessed an allele of Gh_D02G0276 that showed non-synonymous mutations and was prematurely truncated. Furthermore, a Gh_D02G0276-specific marker for the resistance allele variant was able to identify RKN-resistant germplasm from a collection of 367 cotton accessions. The Gh_D02G0276 peptide shares similarity with domesticated hAT-like transposases with additional novel N- and C-terminal domains that resemble the target of known RKN effector molecules and a prokaryotic motif, respectively. The truncation in the resistance allele results in a loss of a plant nuclear gene-specific C-terminal motif, potentially rendering this domain antigenic due to its high homology with bacterial proteins. The conclusive identification of this RKN resistance gene opens new avenues for understanding plant resistance mechanisms to RKN as well as opportunities to develop more efficient marker-assisted selection in cotton breeding programs.

Key Message: Significant associations between candidate genes and six major cotton fiber quality traits were identified in a MAGIC population using GWAS and whole genome sequencing. Upland cotton (Gossypium hirsutum L.) is the world's major renewable source of fibers for textiles. To identify causative genetic variants that influence the major agronomic measures of cotton fiber quality, which are used to set discount or premium prices on each bale of cotton in the USA, we measured six fiber phenotypes from twelve environments, across three locations and 7 years. Our 550 recombinant inbred lines were derived from a multi-parent advanced generation intercross population and were whole-genome-sequenced at 3× coverage, along with the eleven parental cultivars at 20× coverage. The segregation of 473,517 single nucleotide polymorphisms (SNPs) in this population, including 7506 non-synonymous mutations, was combined with phenotypic data to identify seven highly significant fiber quality loci. At these loci, we found fourteen genes with non-synonymous SNPs. Among these loci, some had simple additive effects, while others were only important in a subset of the population. We observed additive effects for elongation and micronaire, when the three most significant loci for each trait were examined. In an informative subset where the major multi-trait locus on chromosome A07:72-Mb was fixed, we unmasked the identity of another significant fiber strength locus in gene Gh_D13G1792 on chromosome D13. The micronaire phenotype only revealed one highly significant genetic locus at one environmental location, demonstrating a significant genetic by environment component. These loci and candidate causative variant alleles will be useful to cotton breeders for marker-assisted selection with minimal linkage drag and potential biotechnological applications.

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.

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

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 2 ) is a monogenic dominant cotton fiber mutation that causes extreme reduction in lint fiber length with no pleiotropic effects on vegetative growth. Li 2 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 2 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 2 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 2 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 2 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 2 ., (Copyright © 2017. Published by Elsevier B.V.)

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

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 1 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 1 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 1 fiber cells accumulated disoriented transverse cables. The polymerization-defective Gly65Val allele in Li 1 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 1 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., (Published 2017. This article is a U.S. Government work and is in the public domain in the USA.)

Background: Cotton supplies a great majority of natural fiber for the global textile industry. The negative correlation between yield and fiber quality has hindered breeders' ability to improve these traits simultaneously. A multi-parent advanced generation inter-cross (MAGIC) population developed through random-mating of multiple diverse parents has the ability to break this negative correlation. Genotyping-by-sequencing (GBS) is a method that can rapidly identify and genotype a large number of single nucleotide polymorphisms (SNP). Genotyping a MAGIC population using GBS technologies will enable us to identify marker-trait associations with high resolution., Results: An Upland cotton MAGIC population was developed through random-mating of 11 diverse cultivars for five generations. In this study, fiber quality data obtained from four environments and 6071 SNP markers generated via GBS and 223 microsatellite markers of 547 recombinant inbred lines (RILs) of the MAGIC population were used to conduct a genome wide association study (GWAS). By employing a mixed linear model, GWAS enabled us to identify markers significantly associated with fiber quantitative trait loci (QTL). We identified and validated one QTL cluster associated with four fiber quality traits [short fiber content (SFC), strength (STR), length (UHM) and uniformity (UI)] on chromosome A07. We further identified candidate genes related to fiber quality attributes in this region. Gene expression and amino acid substitution analysis suggested that a regeneration of bulb biogenesis 1 (GhRBB1_A07) gene is a candidate for superior fiber quality in Upland cotton. The DNA marker CFBid0004 designed from an 18 bp deletion in the coding sequence of GhRBB1_A07 in Acala Ultima is associated with the improved fiber quality in the MAGIC RILs and 105 additional commercial Upland cotton cultivars., Conclusion: Using GBS and a MAGIC population enabled more precise fiber QTL mapping in Upland cotton. The fiber QTL and associated markers identified in this study can be used to improve fiber quality through marker assisted selection or genomic selection in a cotton breeding program. Target manipulation of the GhRBB1_A07 gene through biotechnology or gene editing may potentially improve cotton fiber quality.

Key Message: Three QTL regions controlling three fiber quality traits were validated and further fine-mapped with 27 new single nucleotide polymorphism (SNP) markers. Transcriptome analysis suggests that receptor-like kinases found within the validated QTLs are potential candidate genes responsible for superior fiber strength in cotton line MD52ne. Fiber strength, length, maturity and fineness determine the market value of cotton fibers and the quality of spun yarn. Cotton fiber strength has been recognized as a critical quality attribute in the modern textile industry. Fine mapping along with quantitative trait loci (QTL) validation and candidate gene prediction can uncover the genetic and molecular basis of fiber quality traits. Four previously-identified QTLs (qFBS-c3, qSFI-c14, qUHML-c14 and qUHML-c24) related to fiber bundle strength, short fiber index and fiber length, respectively, were validated using an F3 population that originated from a cross of MD90ne × MD52ne. A group of 27 new SNP markers generated from mapping-by-sequencing (MBS) were placed in QTL regions to improve and validate earlier maps. Our refined QTL regions spanned 4.4, 1.8 and 3.7 Mb of physical distance in the Gossypium raimondii reference genome. We performed RNA sequencing (RNA-seq) of 15 and 20 days post-anthesis fiber cells from MD52ne and MD90ne and aligned reads to the G. raimondii genome. The QTL regions contained 21 significantly differentially expressed genes (DEGs) between the two near-isogenic parental lines. SNPs that result in non-synonymous substitutions to amino acid sequences of annotated genes were identified within these DEGs, and mapped. Taken together, transcriptome and amino acid mutation analysis indicate that receptor-like kinase pathway genes are likely candidates for superior fiber strength and length in MD52ne. MBS along with RNA-seq demonstrated a powerful strategy to elucidate candidate genes for the QTLs that control complex traits in a complex genome like tetraploid upland cotton.

Cotton seed trichomes are the most important source of natural fibers globally. The major fiber thickness properties influence the price of the raw material, and the quality of the finished product. The recessive immature fiber (im) gene reduces the degree of fiber cell wall thickening by a process that was previously shown to involve mitochondrial function in allotetraploid Gossypium hirsutum Here, we present the fine genetic mapping of the im locus, gene expression analysis of annotated proteins near the locus, and association analysis of the linked markers. Mapping-by-sequencing identified a 22-bp deletion in a pentatricopeptide repeat (PPR) gene that is completely linked to the immature fiber phenotype in 2837 F2 plants, and is absent from all 163 cultivated varieties tested, although other closely linked marker polymorphisms are prevalent in the diversity panel. This frame-shift mutation results in a transcript with two long open reading frames: one containing the N-terminal transit peptide that targets mitochondria, the other containing only the RNA-binding PPR domains, suggesting that a functional PPR protein cannot be targeted to mitochondria in the im mutant. Taken together, these results suggest that PPR gene Gh_A03G0489 is involved in the cotton fiber wall thickening process, and is a promising candidate gene at the im locus. Our findings expand our understanding of the molecular mechanisms that modulate cotton fiber fineness and maturity, and may facilitate the development of cotton varieties with superior fiber attributes., (Copyright © 2016 Thyssen et al.)

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.

Background: Individual fiber strength is an important quality attribute that greatly influences the strength of the yarn spun from cotton fibers. Fiber strength is usually measured from bundles of fibers due to the difficulty of reliably measuring strength from individual cotton fibers. However, bundle fiber strength (BFS) is not always correlated with yarn strength since it is affected by multiple fiber properties involved in fiber-to-fiber interactions within a bundle in addition to the individual fiber strength. Molecular mechanisms responsible for regulating individual fiber strength remain unknown. Gossypium hirsutum near isogenic lines (NILs), MD52ne and MD90ne showing variations in BFS provide an opportunity for dissecting the regulatory mechanisms involved in individual fiber strength., Results: Comprehensive fiber property analyses of the NILs revealed that the superior bundle strength of MD52ne fibers resulted from high individual fiber strength with minor contributions from greater fiber length. Comparative transcriptome analyses of the NILs showed that the superior bundle strength of MD52ne fibers was potentially related to two signaling pathways: one is ethylene and the interconnected phytohormonal pathways that are involved in cotton fiber elongation, and the other is receptor-like kinases (RLKs) signaling pathways that are involved in maintaining cell wall integrity. Multiple RLKs were differentially expressed in MD52ne fibers and localized in genomic regions encompassing the strength quantitative trait loci (QTLs). Several candidate genes involved in crystalline cellulose assembly were also up-regulated in MD52ne fibers while the secondary cell wall was produced., Conclusion: Comparative phenotypic and transcriptomic analyses revealed differential expressions of the genes involved in crystalline cellulose assembly, ethylene and RLK signaling pathways between the MD52ne and MD90ne developing fibers. Ethylene and its phytohormonal network might promote the elongation of MD52ne fibers and indirectly contribute to the bundle strength by potentially improving fiber-to-fiber interactions. RLKs that were suggested to mediate a coordination of cell elongation and SCW biosynthesis in other plants might be candidate genes for regulating cotton fiber cell wall assembly and strength.

Key Message: 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.

The number of cotton (Gossypium sp.) ovule epidermal cells differentiating into fiber initials is an important factor affecting cotton yield and fiber quality. Despite extensive efforts in determining the molecular mechanisms regulating fiber initial differentiation, only a few genes responsible for fiber initial differentiation have been discovered. To identify putative genes directly involved in the fiber initiation process, we used a cotton ovule culture technique that controls the timing of fiber initial differentiation by exogenous phytohormone application in combination with comparative expression analyses between wild type and three fiberless mutants. The addition of exogenous auxin and gibberellins to pre-anthesis wild type ovules that did not have visible fiber initials increased the expression of genes affecting auxin, ethylene, ABA and jasmonic acid signaling pathways within 1 h after treatment. Most transcripts expressed differentially by the phytohormone treatment in vitro were also differentially expressed in the ovules of wild type and fiberless mutants that were grown in planta. In addition to MYB25-like, a gene that was previously shown to be associated with the differentiation of fiber initials, several other differentially expressed genes, including auxin/indole-3-acetic acid (AUX/IAA) involved in auxin signaling, ACC oxidase involved in ethylene biosynthesis, and abscisic acid (ABA) 8'-hydroxylase an enzyme that controls the rate of ABA catabolism, were co-regulated in the pre-anthesis ovules of both wild type and fiberless mutants. These results support the hypothesis that phytohormonal signaling networks regulate the temporal expression of genes responsible for differentiation of cotton fiber initials in vitro and in planta.

The presence of two closely related subgenomes in the allotetraploid Upland cotton, combined with a narrow genetic base of the cultivated varieties, has hindered the identification of polymorphic genetic markers and their use in improving this important crop. Genotyping-by-sequencing (GBS) is a rapid way to identify single nucleotide polymorphism (SNP) markers; however, these SNPs may be specific to the sequenced cotton lines. Our objective was to obtain a large set of polymorphic SNPs with broad applicability to the cultivated cotton germplasm. We selected 11 diverse cultivars and their random-mated recombinant inbred progeny for SNP marker development via GBS. Two different GBS methodologies were used by Data2Bio (D2B) and the Institute for Genome Diversity (IGD) to identify 4441 and 1176 polymorphic SNPs with minor allele frequency of ≥0.1, respectively. We further filtered the SNPs and aligned their sequences to the diploid Gossypium raimondii reference genome. We were able to use homeologous SNPs to assign 1071 SNP loci to the At subgenome and 1223 to the Dt subgenome. These filtered SNPs were located in genic regions about twice as frequently as expected by chance. We tested 111 of the SNPs in 154 diverse Upland cotton lines, which confirmed the utility of the SNP markers developed in such approach. Not only were the SNPs identified in the 11 cultivars present in the 154 cotton lines, no two cultivars had identical SNP genotypes. We conclude that GBS can be easily used to discover SNPs in Upland cotton, which can be converted to functional genotypic assays for use in breeding and genetic studies., (© 2015 The Authors.)

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.

Key Message: 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.

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., (Published by Elsevier B.V.)