Don C. Jones

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.

Genome-wide identification and expression analysis of TCP transcription factors in Gossypium raimondii

Plant-specific TEOSINTE-BRANCHED1/CYCLOIDEA/PCF (TCP) transcription factors play versatile functions in multiple aspects of plant growth and development. However, no systematical study has been performed in cotton. In this study, we performed for the first time the genome-wide identification and expression analysis of the TCP transcription factor family in Gossypium raimondii. A total of 38 non-redundant cotton TCP encoding genes were identified. The TCP transcription factors were divided into eleven subgroups based on phylogenetic analysis. Most TCP genes within the same subfamily demonstrated similar exon and intron organization and the motif structures were highly conserved among the subfamilies. Additionally, the chromosomal distribution pattern revealed that TCP genes were unevenly distributed across 11 out of the 13 chromosomes; segmental duplication is a predominant duplication event for TCP genes and the major contributor to the expansion of TCP gene family in G. raimondii. Moreover, the expression profiles of TCP genes shed light on their functional divergence.

Genome-Wide Identification of R2R3-MYB Genes and Expression Analyses During Abiotic Stress in Gossypium raimondii.

The R2R3-MYB is one of the largest families of transcription factors, which have been implicated in multiple biological processes. There is great diversity in the number of R2R3-MYB genes in different plants. However, there is no report on genome-wide characterization of this gene family in cotton. In the present study, a total of 205 putative R2R3-MYB genes were identified in cotton D genome (Gossypium raimondii), that are much larger than that found in other cash crops with fully sequenced genomes. These GrMYBs were classified into 13 groups with the R2R3-MYB genes from Arabidopsis and rice. The amino acid motifs and phylogenetic tree were predicted and analyzed. The sequences of GrMYBs were distributed across 13 chromosomes at various densities. The results showed that the expansion of the G. Raimondii R2R3-MYB family was mainly attributable to whole genome duplication and segmental duplication. Moreover, the expression pattern of 52 selected GrMYBs and 46 GaMYBs were tested in roots and leaves under different abiotic stress conditions. The results revealed that the MYB genes in cotton were differentially expressed under salt and drought stress treatment. Our results will be useful for determining the precise role of the MYB genes during stress responses with crop improvement.

RNA-Seq transcriptome profiling of upland cotton (Gossypium hirsutum L.) root tissue under water-deficit stress.

An RNA-Seq experiment was performed using field grown well-watered and naturally rain fed cotton plants to identify differentially expressed transcripts under water-deficit stress. Our work constitutes the first application of the newly published diploid D5 Gossypium raimondii sequence in the study of tetraploid AD1 upland cotton RNA-seq transcriptome analysis. A total of 1,530 transcripts were differentially expressed between well-watered and water-deficit stressed root tissues, in patterns that confirm the accuracy of this technique for future studies in cotton genomics. Additionally, putative sequence based genome localization of differentially expressed transcripts detected A2 genome specific gene expression under water-deficit stress. These data will facilitate efforts to understand the complex responses governing transcriptomic regulatory mechanisms and to identify candidate genes that may benefit applied plant breeding programs.

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.

Development and bin mapping of gene-associated interspecific SNPs for cotton ( Gossypium hirsutum L.) introgression breeding efforts

BackgroundCotton (Gossypium spp.) is the largest producer of natural fibers for textile and is an important crop worldwide. Crop production is comprised primarily of G. hirsutum L., an allotetraploid. However, elite cultivars express very small amounts of variation due to the species monophyletic origin, domestication and further bottlenecks due to selection. Conversely, wild cotton species harbor extensive genetic diversity of prospective utility to improve many beneficial agronomic traits, fiber characteristics, and resistance to disease and drought. Introgression of traits from wild species can provide a natural way to incorporate advantageous traits through breeding to generate higher-producing cotton cultivars and more sustainable production systems. Interspecific introgression efforts by conventional methods are very time-consuming and costly, but can be expedited using marker-assisted selection.ResultsUsing transcriptome sequencing we have developed the first gene-associated single nucleotide polymorphism (SNP) markers for wild cotton species G. tomentosum, G. mustelinum, G. armourianum and G. longicalyx. Markers were also developed for a secondary cultivated species G. barbadense cv. 3–79. A total of 62,832 non-redundant SNP markers were developed from the five wild species which can be utilized for interspecific germplasm introgression into cultivated G. hirsutum and are directly associated with genes. Over 500 of the G. barbadense markers have been validated by whole-genome radiation hybrid mapping. Overall 1,060 SNPs from the five different species have been screened and shown to produce acceptable genotyping assays.ConclusionsThis large set of 62,832 SNPs relative to cultivated G. hirsutum will allow for the first high-density mapping of genes from five wild species that affect traits of interest, including beneficial agronomic and fiber characteristics. Upon mapping, the markers can be utilized for marker-assisted introgression of new germplasm into cultivated cotton and in subsequent breeding of agronomically adapted types, including cultivar development.

Comprehensive analysis of TCP transcription factors and their expression during cotton ( Gossypium arboreum ) fiber early development

TCP proteins are plant-specific transcription factors implicated to perform a variety of physiological functions during plant growth and development. In the current study, we performed for the first time the comprehensive analysis of TCP gene family in a diploid cotton species, Gossypium arboreum, including phylogenetic analysis, chromosome location, gene duplication status, gene structure and conserved motif analysis, as well as expression profiles in fiber at different developmental stages. Our results showed that G. arboreum contains 36 TCP genes, distributing across all of the thirteen chromosomes. GaTCPs within the same subclade of the phylogenetic tree shared similar exon/intron organization and motif composition. In addition, both segmental duplication and whole-genome duplication contributed significantly to the expansion of GaTCPs. Many these TCP transcription factor genes are specifically expressed in cotton fiber during different developmental stages, including cotton fiber initiation and early development. This suggests that TCP genes may play important roles in cotton fiber development.

Genome-wide identification of auxin response factor (ARF) genes and its tissue-specific prominent expression in Gossypium raimondii

Auxin response factors (ARFs) are recently discovered transcription factors that bind with auxin response elements (AuxRE, TGTCTC) to regulate the expression of early auxin-responsive genes. To our knowledge, the ARF gene family has never been characterized in cotton, the most important fiber crop in the world. In this study, a total of 35 ARF genes, named as GrARFs, were identified in a diploid cotton species Gossypium raimondii. The 35 ARF genes were located in 12 of the 13 cotton chromosomes; the intron/exon distribution of the GrARF genes was similar among sister pairs, whereas the divergence of some GrARF genes suggests the possibility of functional diversification. Our results show that the middle domains of nine GrARF proteins rich in glutamine (Q) are activators, while 26 other GrARF proteins rich in proline (P), serine (S), and threonine (T) are repressors. Our results also show that the expression of GrARF genes is diverse in different tissues. The expression of GrARF1 was significantly higher in leaves, whereas GrARF2a had higher expression level in shoots, which implicates different roles in the tested tissues. The GrARF11 has a higher expression level in buds than that in leaves, while GrARF19.2 shows contrasting expression patterns, having higher expression in leaves than that in buds. This suggests that they play different roles in leaves and buds. During long-term evolution of G. raimondii, some ARF genes were lost and some arose. The identification and characterization of the ARF genes in G. raimondii elucidate its important role in cotton that ARF genes regulate the development of flower buds, sepals, shoots, and leaves.

Modifications to a LATE MERISTEM IDENTITY1 gene are responsible for the major leaf shapes of Upland cotton (Gossypium hirsutum L.)

Significance Leaves are the primary source of photoassimilate in crop plants. A precise understanding of the genetic architecture underlying leaf morphology is critical to engineering climate-resilient crop varieties. An ideal cotton cultivar would produce a lower canopy of broad, normal leaves before transitioning to an upper canopy of highly lobed, okra leaves. Here we show that the major leaf shapes of cotton are controlled by the okra locus, which encodes an HD-Zip transcription factor Gossypium hirsutum LATE MERISTEM IDENTITY1-D1b (GhLMI1-D1b). Using gene silencing, we temporarily induced normal leaf formation in okra, thus validating the candidate gene and creating the leaf shape ideotype in cotton. This study, identifying a single locus responsible for cotton leaf shape, expands the genetic toolbox for breeders to produce superior cotton varieties. Leaf shape varies spectacularly among plants. Leaves are the primary source of photoassimilate in crop plants, and understanding the genetic basis of variation in leaf morphology is critical to improving agricultural productivity. Leaf shape played a unique role in cotton improvement, as breeders have selected for entire and lobed leaf morphs resulting from a single locus, okra (l-D1), which is responsible for the major leaf shapes in cotton. The l-D1 locus is not only of agricultural importance in cotton, but through pioneering chimeric and morphometric studies, it has contributed to fundamental knowledge about leaf development. Here we show that an HD-Zip transcription factor homologous to the LATE MERISTEM IDENTITY1 (LMI1) gene of Arabidopsis is the causal gene underlying the l-D1 locus. The classical okra leaf shape allele has a 133-bp tandem duplication in the promoter, correlated with elevated expression, whereas an 8-bp deletion in the third exon of the presumed wild-type normal allele causes a frame-shifted and truncated coding sequence. Our results indicate that subokra is the ancestral leaf shape of tetraploid cotton that gave rise to the okra allele and that normal is a derived mutant allele that came to predominate and define the leaf shape of cultivated cotton. Virus-induced gene silencing (VIGS) of the LMI1-like gene in an okra variety was sufficient to induce normal leaf formation. The developmental changes in leaves conferred by this gene are associated with a photosynthetic transcriptomic signature, substantiating its use by breeders to produce a superior cotton ideotype.