David D. Fang

A High Density Consensus Genetic Map of Tetraploid Cotton That Integrates Multiple Component Maps through Molecular Marker Redundancy Check

A consensus genetic map of tetraploid cotton was constructed using six high-density maps and after the integration of a sequence-based marker redundancy check. Public cotton SSR libraries (17,343 markers) were curated for sequence redundancy using 90% as a similarity cutoff. As a result, 20% of the markers (3,410) could be considered as redundant with some other markers. The marker redundancy information had been a crucial part of the map integration process, in which the six most informative interspecific Gossypium hirsutum×G. barbadense genetic maps were used for assembling a high density consensus (HDC) map for tetraploid cotton. With redundant markers being removed, the HDC map could be constructed thanks to the sufficient number of collinear non-redundant markers in common between the component maps. The HDC map consists of 8,254 loci, originating from 6,669 markers, and spans 4,070 cM, with an average of 2 loci per cM. The HDC map presents a high rate of locus duplications, as 1,292 markers among the 6,669 were mapped in more than one locus. Two thirds of the duplications are bridging homoeologous AT and DT chromosomes constitutive of allopolyploid cotton genome, with an average of 64 duplications per AT/DT chromosome pair. Sequences of 4,744 mapped markers were used for a mutual blast alignment (BBMH) with the 13 major scaffolds of the recently released Gossypium raimondii genome indicating high level of homology between the diploid D genome and the tetraploid cotton genetic map, with only a few minor possible structural rearrangements. Overall, the HDC map will serve as a valuable resource for trait QTL comparative mapping, map-based cloning of important genes, and better understanding of the genome structure and evolution of tetraploid cotton.

A High-Density Simple Sequence Repeat and Single Nucleotide Polymorphism Genetic Map of the Tetraploid Cotton Genome

Genetic linkage maps play fundamental roles in understanding genome structure, explaining genome formation events during evolution, and discovering the genetic bases of important traits. A high-density cotton (Gossypium spp.) genetic map was developed using representative sets of simple sequence repeat (SSR) and the first public set of single nucleotide polymorphism (SNP) markers to genotype 186 recombinant inbred lines (RILs) derived from an interspecific cross between Gossypium hirsutum L. (TM-1) and G. barbadense L. (3-79). The genetic map comprised 2072 loci (1825 SSRs and 247 SNPs) and covered 3380 centiMorgan (cM) of the cotton genome (AD) with an average marker interval of 1.63 cM. The allotetraploid cotton genome produced equivalent recombination frequencies in its two subgenomes (At and Dt). Of the 2072 loci, 1138 (54.9%) were mapped to 13 At-subgenome chromosomes, covering 1726.8 cM (51.1%), and 934 (45.1%) mapped to 13 Dt-subgenome chromosomes, covering 1653.1 cM (48.9%). The genetically smallest homeologous chromosome pair was Chr. 04 (A04) and 22 (D04), and the largest was Chr. 05 (A05) and 19 (D05). Duplicate loci between and within homeologous chromosomes were identified that facilitate investigations of chromosome translocations. The map augments evidence of reciprocal rearrangement between ancestral forms of Chr. 02 and 03 versus segmental homeologs 14 and 17 as centromeric regions show homeologous between Chr. 02 (A02) and 17 (D02), as well as between Chr. 03 (A03) and 14 (D03). This research represents an important foundation for studies on polyploid cottons, including germplasm characterization, gene discovery, and genome sequence assembly.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp.

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.

Linkage Map Construction and Quantitative Trait Locus Analysis of Agronomic and Fiber Quality Traits in Cotton

The superior fiber properties of Gossypium barbadense L. serve as a source of novel variation for improving fiber quality in Upland cotton (G. hirsutum L.), but introgression from G. barbadense has been largely unsuccessful due to hybrid breakdown and a lack of genetic and genomic resources. In an effort to overcome these limitations, we constructed a linkage map and conducted a quantitative trait locus (QTL) analysis of 10 agronomic and fiber quality traits in a recombinant inbred mapping population derived from a cross between TM‐1, an Upland cotton line, and NM24016, an elite G. hirsutum line with stabilized introgression from G. barbadense. The linkage map consisted of 429 simple‐sequence repeat (SSR) and 412 genotyping‐by‐sequencing (GBS)‐based single‐nucleotide polymorphism (SNP) marker loci that covered half of the tetraploid cotton genome. Notably, the 841 marker loci were unevenly distributed among the 26 chromosomes of tetraploid cotton. The 10 traits evaluated on the TM‐1 × NM24016 population in a multienvironment trial were highly heritable, and most of the fiber traits showed considerable transgressive variation. Through the QTL analysis, we identified a total of 28 QTLs associated with the 10 traits. Our study provides a novel resource that can be used by breeders and geneticists for the genetic improvement of agronomic and fiber quality traits in Upland cotton.

A combined functional and structural genomics approach identified an EST-SSR marker with complete linkage to the Ligon lintless-2 genetic locus in cotton (Gossypium hirsutum L.)

BackgroundCotton fiber length is an important quality attribute to the textile industry and longer fibers can be more efficiently spun into yarns to produce superior fabrics. There is typically a negative correlation between yield and fiber quality traits such as length. An understanding of the regulatory mechanisms controlling fiber length can potentially provide a valuable tool for cotton breeders to improve fiber length while maintaining high yields. The cotton (Gossypium hirsutum L.) fiber mutation Ligon lintless-2 is controlled by a single dominant gene (Li2) that results in significantly shorter fibers than a wild-type. In a near-isogenic state with a wild-type cotton line, Li2 is a model system with which to study fiber elongation.ResultsTwo near-isogenic lines of Ligon lintless-2 (Li2) cotton, one mutant and one wild-type, were developed through five generations of backcrosses (BC5). An F2 population was developed from a cross between the two Li2 near-isogenic lines and used to develop a linkage map of the Li2 locus on chromosome 18. Five simple sequence repeat (SSR) markers were closely mapped around the Li2 locus region with two of the markers flanking the Li2 locus at 0.87 and 0.52 centimorgan. No apparent differences in fiber initiation and early fiber elongation were observed between the mutant ovules and the wild-type ones. Gene expression profiling using microarrays suggested roles of reactive oxygen species (ROS) homeostasis and cytokinin regulation in the Li2 mutant phenotype. Microarray gene expression data led to successful identification of an EST-SSR marker (NAU3991) that displayed complete linkage to the Li2 locus.ConclusionsIn the field of cotton genomics, we report the first successful conversion of gene expression data into an SSR marker that is associated with a genomic region harboring a gene responsible for a fiber trait. The EST-derived SSR marker NAU3991 displayed complete linkage to the Li2 locus on chromosome 18 and resided in a gene with similarity to a putative plectin-related protein. The complete linkage suggests that this expressed sequence may be the Li2 gene.

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

BackgroundCotton 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 (Li1) 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 Li1 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.ResultsPhysical and morphological measurements of the Li1 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 F2 individuals identified SSR markers associated with the Li1 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 Li1 gene.ConclusionsThe early termination of fiber elongation in the Li1 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 Li1 mutant were identified. Molecular markers closely associated with the Li1 locus were developed. Results presented here will lay the foundation for further investigation of the genetic and molecular mechanisms of fiber elongation.

Quantitative trait locus analysis of Verticillium wilt resistance in an introgressed recombinant inbred population of Upland cotton

Verticillium wilt (VW) of Upland cotton (Gossypium hirsutum L.) is caused by the soil-borne fungal pathogen Verticillium dahlia Kleb. The availability of VW-resistant cultivars is vital for control of this economically important disease, but there is a paucity of Upland cotton breeding lines and cultivars with a high level of resistance to VW. In general, G. barbadense L. (source of Pima cotton) is more VW-resistant than Upland cotton. However, the transfer of VW resistance from G. barbadense to Upland cotton is challenging because of hybrid breakdown in the F2 and successive generations of interspecific populations. We conducted two replicated greenhouse studies (tests 1 and 2) to assess the heritability of VW resistance to a defoliating V. dahliae isolate and identify genetic markers associated with VW resistance in an Upland cotton recombinant inbred mapping population that has stable introgression from Pima cotton. Disease ratings at the seedling stage on several different days after the first inoculation (DAI) in test 1, as well as the percentages of infected and defoliated leaves at 2 DAI in test 2, were found to be low to moderately heritable, indicating the importance of a replicated progeny test in selection for VW resistance. With a newly constructed linkage map consisting of 882 simple sequence repeat, single nucleotide polymorphism, and resistance gene analog–amplified fragment length polymorphism marker loci, we identified a total of 21 quantitative trait loci (QTLs) on 11 chromosomes and two linkage groups associated with VW resistance at several different DAIs in greenhouse tests 1 and 2. The markers associated with the VW resistance QTLs will facilitate fine mapping and cloning of VW resistance genes and genomics-assisted breeding for VW-resistant cultivars.

Development of a core set of SSR markers for the characterization of Gossypium germplasm

Molecular markers such as simple sequence repeats (SSR) are a useful tool for characterizing genetic diversity of Gossypium germplasm. Genetic profiles by DNA fingerprinting of cotton accessions can only be compared among different collections if a common set of molecular markers are used by different laboratories and/or research projects. Herein, we propose and report a core set of 105 SSR markers with wide genome coverage of at least four evenly distributed markers per chromosome for the 26 tetraploid cotton chromosomes. The core marker set represents the efforts of ten research groups involved in marker development, and have been systematically evaluated for DNA polymorphism on the 12 genotypes belonging to six Gossypium species [known collectively as the cotton marker database (CMD) panel]. A total of 35 marker bins in triplex sets were arranged from the 105 markers that were each labeled with one of the three fluorescent dyes (FAM, HEX, and NED). Results from this study indicated that the core marker set was robust in revealing DNA polymorphism either between and within species. Average value of polymorphism information content (PIC) among the CMD panel was 0.65, and that within the cultivated cotton species Gossypium hirsutum was 0.29. Based on the similarity matrix and phylogenetic analysis of the CMD panel, the core marker set appeared to be sufficient in characterizing the diversity within G. hirsutum and other Gossypium species. The portability of this core marker set would facilitate the systematic characterization and the simultaneous comparison among various research efforts involved in genetic diversity analysis and germplasm resource preservation.

Detection, Validation, and Application of Genotyping-by-Sequencing Based Single Nucleotide Polymorphisms in Upland Cotton

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.

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

BackgroundThe 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 (Li2) and results in extremely shortened lint fibers on mature seeds with no visible pleiotropic effects on vegetative growth and development. The Li2 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 Li2 mutant fibers were compared to wild-type fibers during development.ResultsPrincipal component analysis of metabolites from GC-MS data separated Li2 mutant fiber samples from WT fiber samples at the WT elongation stage, indicating that the Li2 mutation altered the metabolome of the mutant fibers. The observed alterations in the Li2 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 Li2 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 Li2 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 Li2 mutant fibers.ConclusionsThis report provides the first comprehensive analysis of metabolite and transcript changes in response to the Li2 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.