Bidirectional interlocus concerted evolution following allopolyploid speciation in cotton (Gossypium)

Polyploidy is a prominent process in plant evolution; yet few data address the question of whether homeologous sequences evolve independently subsequent to polyploidization. We report on ribosomal DNA (rDNA) evolution in five allopolyploid (AD genome) species of cotton (Gossypium) and species representing their diploid progenitors (A genome, D genome). Sequence data from the internal transcribed spacer regions (ITS1 and ITS2) and the 5.8S gene indicate that rDNA arrays are homogeneous, or nearly so, in all diploids and allopolyploids examined. Because these arrays occur at four chromosomal loci in allopolyploid cotton, two in each subgenome, repeats from different arrays must have become homogenized by interlocus concerted evolution. Southern hybridization analysis combined with copy-number estimation demonstrate that this process has gone to completion in the diploids and to completion or near-completion in all allopolyploid species and that it most likely involves the entire rDNA repeat. Phylogenetic analysis demonstrates that interlocus concerted evolution has been bidirectional in allopolyploid species--i.e., rDNA from four polyploid lineages has been homogenized to a D genome repeat type, whereas sequences from Gossypium mustelinum have concerted to an A genome repeat type. Although little is known regarding the functional significance of interlocus concerted evolution of homeologous sequences, this study demonstrates that the process occurs for tandemly repeated sequences in diploid and polyploid plants. That interlocus concerted evolution can occur bidirectionally subsequent to hybidization and polyploidization has significant implications for phylogeny reconstruction, especially when based on rDNA sequences.

Meta-analysis of polyploid cotton QTL shows unequal contributions of subgenomes to a complex network of genes and gene clusters implicated in lint fiber development

QTL mapping experiments yield heterogeneous results due to the use of different genotypes, environments, and sampling variation. Compilation of QTL mapping results yields a more complete picture of the genetic control of a trait and reveals patterns in organization of trait variation. A total of 432 QTL mapped in one diploid and 10 tetraploid interspecific cotton populations were aligned using a reference map and depicted in a CMap resource. Early demonstrations that genes from the non-fiber-producing diploid ancestor contribute to tetraploid lint fiber genetics gain further support from multiple populations and environments and advanced-generation studies detecting QTL of small phenotypic effect. Both tetraploid subgenomes contribute QTL at largely non-homeologous locations, suggesting divergent selection acting on many corresponding genes before and/or after polyploid formation. QTL correspondence across studies was only modest, suggesting that additional QTL for the target traits remain to be discovered. Crosses between closely-related genotypes differing by single-gene mutants yield profoundly different QTL landscapes, suggesting that fiber variation involves a complex network of interacting genes. Members of the lint fiber development network appear clustered, with cluster members showing heterogeneous phenotypic effects. Meta-analysis linked to synteny-based and expression-based information provides clues about specific genes and families involved in QTL networks.

Rapid DNA loss as a counterbalance to genome expansion through retrotransposon proliferation in plants

Transposable elements, particularly LTR-retrotransposons, comprise the primary vehicle for genome size expansion in plants, while DNA removal through illegitimate recombination and intrastrand homologous recombination serve as the most important counteracting forces to plant genomic obesity. Despite extensive research, the relative impact of these opposing forces and hence the directionality of genome size change remains unknown. In Gossypium (cotton), the 3-fold genome size variation among diploids is due largely to copy number variation of the gypsy-like retrotransposon Gorge3. Here we combine comparative sequence analysis with a modeling approach to study the directionality of genome size change in Gossypium. We demonstrate that the rate of DNA removal in the smaller genomes is sufficient to reverse genome expansion through Gorge3 proliferation. These data indicate that rates of DNA loss can be highly variable even within a single plant genus, and that the known mechanisms of DNA loss can indeed reverse the march toward genomic obesity.

Acute bronchoconstriction induced by cotton dust: dose-related responses to endotoxin and other dust factors

Fifty-four healthy humans, selected for their acute airway responsiveness to cotton dust, had spirometric tests immediately before and after 6 hours of exposure to card-generated cotton dust from seven different cottons (of several grades and growing regions). During exposures, we measured airborne concentrations of viable fungi and bacteria (total and gram negative), vertically elutriated gravimetric dust, and vertically elutriated endotoxin. Correlation between each of these five exposure indices and exposure-related acute changes in forced expiratory volume in 1 s showed a statistically significant relationship between all of the indices except concentration of viable fungi. Of the other four indices, endotoxin was the most highly correlated (r = -0.94; p less than 0.00001), and gravimetric dust was the least correlated (r = -0.34; p less than 0.05). These findings suggest that gram-negative endotoxin may play a major role in the acute pulmonary response to inhaled cotton dust.

Identification of the family of aquaporin genes and their expression in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Cotton (Gossypium spp.) is produced in over 30 countries and represents the most important natural fiber in the world. One of the primary factors affecting both the quantity and quality of cotton production is water. A major facilitator of water movement through cell membranes of cotton and other plants are the aquaporin proteins. Aquaporin proteins are present as diverse forms in plants, where they function as transport systems for water and other small molecules. The plant aquaporins belong to the large major intrinsic protein (MIP) family. In higher plants, they consist of five subfamilies including plasma membrane intrinsic proteins (PIP), tonoplast intrinsic proteins (TIP), NOD26-like intrinsic proteins (NIP), small basic intrinsic proteins (SIP), and the recently discovered X intrinsic proteins (XIP). Although a great deal is known about aquaporins in plants, very little is known in cotton.

RESULTS

From a molecular cloning effort, together with a bioinformatic homology search, 71 upland cotton (G. hirsutum) aquaporin genes were identified. The cotton aquaporins consist of 28 PIP and 23 TIP members with high sequence similarity. We also identified 12 NIP and 7 SIP members that showed more divergence. In addition, one XIP member was identified that formed a distinct 5th subfamily. To explore the physiological roles of these aquaporin genes in cotton, expression analyses were performed for a select set of aquaporin genes from each subfamily using semi-quantitative reverse transcription (RT)-PCR. Our results suggest that many cotton aquaporin genes have high sequence similarity and diverse roles as evidenced by analysis of sequences and their expression.

CONCLUSION

This study presents a comprehensive identification of 71 cotton aquaporin genes. Phylogenetic analysis of amino acid sequences divided the large and highly similar multi-gene family into the known 5 aquaporin subfamilies. Together with expression and bioinformatic analyses, our results support the idea that the genes identified in this study represent an important genetic resource providing potential targets to modify the water use properties of cotton.

Differential regulation of six novel MYB-domain genes defines two distinct expression patterns in allotetraploid cotton (Gossypium hirsutum L.)

A PCR-based strategy was employed to identify myb-related genes potentially involved in the differentiation and development of cotton seed trichomes. cDNA clones representing six newly identified cotton myb-domain genes (GhMYB) of the R2R3-MYB family were characterized in the allotetraploid species Gossypium hirsutum L. (2n = 4x = 52; AADD). Several interesting motifs and domains in the transregulatory region (TRR) were identified as potential candidates for modulating GhMYB activity. One such structural feature is a basic 40-amino acid stretch (TRR1) located immediately downstream of the DNA-binding domain (DBD) in five of the GhMYBs. Furthermore, the conserved motif GIDxxH identified in a subset of plant MYBs is also present in the same position in the TRR1 domains of GhMYB1 and GhMYB6, exactly 12 amino acid residues downstream of the last tryptophan in the R3 repeat of the DBD. At least two of the GhMYBs (GhMYB4 and GhMYB5) contain unidentified ORFS in the 5' leader sequence (5'-uORFs) that may serve to regulate the synthesis of these particular GhMYB proteins at the translational level. Multiple alignment of DBD sequences indicated that GhMYBs show structural similarity to plant R2R3-MYB factors implicated in phenylpropanoid biosynthesis. GhMYB5 is the most distantly related cotton R2R3-MYB and is found in an isolated cluster that includes the drought-inducible AtMYB2. Sequence comparisons of DBD and TRR domains from GhMYBs, MIXTA (AmMYBMx) and G11 (AtMYBG11) did not reveal any striking similarity beyond conserved motifs. However, based on earlier phylogenetic analysis, GhMYB2, GhMYB3, and GHMYB4 are members of a cluster that contains GLABROUS1, while GhMYB1 and GhMYB6 belong to a closely related cluster. Semi-quantitative RT-PCR analysis revealed two discrete patterns of GhMYB gene expression. Type I cotton MYB (GhMYB-1, -2, and -3) transcripts were found in all tissue-types examined and were relatively more abundant than those derived from type II GhMYB genes (GhMYB-4, -5, and -6), which showed distinct, tissue-specific expression patterns. The developmental regulation of GhMYBs is consistent with a role for these DNA-binding factors in the differentiation and expansion of cotton seed trichomes.

Comparative development of fiber in wild and cultivated cotton

One of the most striking examples of plant hairs is the single-celled epidermal seed trichome of cultivated cotton. The developmental morphology of these commercial "fibers" has been well-characterized in Gossypium hirsutum, but little is known about the pattern and tempo of fiber development in wild Gossypium species, all of which have short, agronomically inferior fiber. To identify developmental differences that account for variation in fiber length, and to place these differences in a phylogenetic context, we conducted SEM studies of ovules at and near the time of flowering, and generated growth curves for cultivated and wild diploid and tetraploid species. Trichome initiation was found to be similar in all taxa, with few notable differences in trichome density or early growth. Developmental profiles of the fibers of most wild species are similar, with fiber elongation terminating at about two weeks post-anthesis. In contrast, growth is extended to three weeks in the A- and F-genome diploids. This prolonged elongation period is diagnosed as a key evolutionary event in the origin of long fiber. A second evolutionary innovation is that absolute growth rate is higher in species with long fibers. Domestication of species is associated with a further prolongation of elongation at both the diploid and allopolyploid levels, suggesting the effects of parallel artificial selection. Comparative analysis of fiber growth curves lends developmental support to previous quantitative genetic suggestions that genes for fiber "improvement" in tetraploid cotton were contributed by the agronomically inferior D-genome diploid parent.

The complete chloroplast genome sequence of Gossypium hirsutum: organization and phylogenetic relationships to other angiosperms

BACKGROUND

Cotton (Gossypium hirsutum) is the most important fiber crop grown in 90 countries. In 2004-2005, US farmers planted 79% of the 5.7-million hectares of nuclear transgenic cotton. Unfortunately, genetically modified cotton has the potential to hybridize with other cultivated and wild relatives, resulting in geographical restrictions to cultivation. However, chloroplast genetic engineering offers the possibility of containment because of maternal inheritance of transgenes. The complete chloroplast genome of cotton provides essential information required for genetic engineering. In addition, the sequence data were used to assess phylogenetic relationships among the major clades of rosids using cotton and 25 other completely sequenced angiosperm chloroplast genomes.

RESULTS

The complete cotton chloroplast genome is 160,301 bp in length, with 112 unique genes and 19 duplicated genes within the IR, containing a total of 131 genes. There are four ribosomal RNAs, 30 distinct tRNA genes and 17 intron-containing genes. The gene order in cotton is identical to that of tobacco but lacks rpl22 and infA. There are 30 direct and 24 inverted repeats 30 bp or longer with a sequence identity > or = 90%. Most of the direct repeats are within intergenic spacer regions, introns and a 72 bp-long direct repeat is within the psaA and psaB genes. Comparison of protein coding sequences with expressed sequence tags (ESTs) revealed nucleotide substitutions resulting in amino acid changes in ndhC, rpl23, rpl20, rps3 and clpP. Phylogenetic analysis of a data set including 61 protein-coding genes using both maximum likelihood and maximum parsimony were performed for 28 taxa, including cotton and five other angiosperm chloroplast genomes that were not included in any previous phylogenies.

CONCLUSION

Cotton chloroplast genome lacks rpl22 and infA and contains a number of dispersed direct and inverted repeats. RNA editing resulted in amino acid changes with significant impact on their hydropathy. Phylogenetic analysis provides strong support for the position of cotton in the Malvales in the eurosids II clade sister to Arabidopsis in the Brassicales. Furthermore, there is strong support for the placement of the Myrtales sister to the eurosid I clade, although expanded taxon sampling is needed to further test this relationship.

Genome-wide analysis of codon usage bias in four sequenced cotton species

Codon usage bias (CUB) is an important evolutionary feature in a genome which provides important information for studying organism evolution, gene function and exogenous gene expression. The CUB and its shaping factors in the nuclear genomes of four sequenced cotton species, G. arboreum (A₂), G. raimondii (D₅), G. hirsutum (AD₁) and G. barbadense (AD₂) were analyzed in the present study. The effective number of codons (ENC) analysis showed the CUB was weak in these four species and the four subgenomes of the two tetraploids. Codon composition analysis revealed these four species preferred to use pyrimidine-rich codons more frequently than purine-rich codons. Correlation analysis indicated that the base content at the third position of codons affect the degree of codon preference. PR2-bias plot and ENC-plot analyses revealed that the CUB patterns in these genomes and subgenomes were influenced by combined effects of translational selection, directional mutation and other factors. The translational selection (P2) analysis results, together with the non-significant correlation between GC12 and GC3, further revealed that translational selection played the dominant role over mutation pressure in the codon usage bias. Through relative synonymous codon usage (RSCU) analysis, we detected 25 high frequency codons preferred to end with T or A, and 31 low frequency codons inclined to end with C or G in these four species and four subgenomes. Finally, 19 to 26 optimal codons with 19 common ones were determined for each species and subgenomes, which preferred to end with A or T. We concluded that the codon usage bias was weak and the translation selection was the main shaping factor in nuclear genes of these four cotton genomes and four subgenomes.

Molecular diversity and association mapping of fiber quality traits in exotic G. hirsutum L. germplasm

The narrow genetic base of cultivated cotton germplasm is hindering the cotton productivity worldwide. Although potential genetic diversity exists in Gossypium genus, it is largely 'underutilized' due to photoperiodism and the lack of innovative tools to overcome such challenges. The application of linkage disequilibrium (LD)-based association mapping is an alternative powerful molecular tool to dissect and exploit the natural genetic diversity conserved within cotton germplasm collections, greatly accelerating still 'lagging' cotton marker-assisted selection (MAS) programs. However, the extent of genome-wide linkage disequilibrium (LD) has not been determined in cotton. We report the extent of genome-wide LD and association mapping of fiber quality traits by using a 95 core set of microsatellite markers in a total of 285 exotic Gossypium hirsutum accessions, comprising of 208 landrace stocks and 77 photoperiodic variety accessions. We demonstrated the existence of useful genetic diversity within exotic cotton germplasm. In this germplasm set, 11-12% of SSR loci pairs revealed a significant LD. At the significance threshold (r(2)>/=0.1), a genome-wide average of LD declines within the genetic distance at <10 cM in the landrace stocks germplasm and >30 cM in variety germplasm. Genome wide LD at r(2)>/=0.2 was reduced on average to approximately 1-2 cM in the landrace stock germplasm and 6-8 cM in variety germplasm, providing evidence of the potential for association mapping of agronomically important traits in cotton. We observed significant population structure and relatedness in assayed germplasm. Consequently, the application of the mixed liner model (MLM), considering both kinship (K) and population structure (Q) detected between 6% and 13% of SSR markers associated with the main fiber quality traits in cotton. Our results highlight for the first time the feasibility and potential of association mapping, with consideration of the population structure and stratification existing in cotton germplasm resources. The number of SSR markers associated with fiber quality traits in diverse cotton germplasm, which broadly covered many historical meiotic events, should be useful to effectively exploit potentially new genetic variation by using MAS programs.

An L1 box binding protein, GbML1, interacts with GbMYB25 to control cotton fibre development

Transcription factors play key roles in plant development through their interaction with cis-elements and/or other transcription factors. A HD-Zip IV family transcription factor, Gossypium barbadense Meristem Layer 1 (GbML1) has been identified and characterized here. GbML1 specifically bound to the L1 box and the promoters of GbML1 and GbRDL1. GbML1 physically interacted with a key regulator of cotton fibre development, GbMYB25. Truncated and point mutation assays indicated the START-SAD domain was required for the binding to the C terminal domain (CTD) of GbMYB25. GbML1 overexpression in Arabidopsis increased the number of trichomes on stems and leaves and increased the accumulation of anthocyanin in leaves. Taken together, the L1 box binding protein, GbML1 was identified as the first partner for GbMYB25 and the role of START domain was discovered to be a protein binding domain in plants. Our findings will help the improvement of cotton fibre production and the understanding of the key role of HD-Zip family and MYB family in plants.

Construction of a high-density genetic map by specific locus amplified fragment sequencing (SLAF-seq) and its application to Quantitative Trait Loci (QTL) analysis for boll weight in upland cotton (Gossypium hirsutum.)

BACKGROUND

Upland Cotton (Gossypium hirsutum) is one of the most important worldwide crops it provides natural high-quality fiber for the industrial production and everyday use. Next-generation sequencing is a powerful method to identify single nucleotide polymorphism markers on a large scale for the construction of a high-density genetic map for quantitative trait loci mapping.

RESULTS

In this research, a recombinant inbred lines population developed from two upland cotton cultivars 0-153 and sGK9708 was used to construct a high-density genetic map through the specific locus amplified fragment sequencing method. The high-density genetic map harbored 5521 single nucleotide polymorphism markers which covered a total distance of 3259.37 cM with an average marker interval of 0.78 cM without gaps larger than 10 cM. In total 18 quantitative trait loci of boll weight were identified as stable quantitative trait loci and were detected in at least three out of 11 environments and explained 4.15-16.70 % of the observed phenotypic variation. In total, 344 candidate genes were identified within the confidence intervals of these stable quantitative trait loci based on the cotton genome sequence. These genes were categorized based on their function through gene ontology analysis, Kyoto Encyclopedia of Genes and Genomes analysis and eukaryotic orthologous groups analysis.

CONCLUSIONS

This research reported the first high-density genetic map for Upland Cotton (Gossypium hirsutum) with a recombinant inbred line population using single nucleotide polymorphism markers developed by specific locus amplified fragment sequencing. We also identified quantitative trait loci of boll weight across 11 environments and identified candidate genes within the quantitative trait loci confidence intervals. The results of this research would provide useful information for the next-step work including fine mapping, gene functional analysis, pyramiding breeding of functional genes as well as marker-assisted selection.

Incongruent patterns of local and global genome size evolution in cotton

Genome sizes in plants vary over several orders of magnitude, reflecting a combination of differentially acting local and global forces such as biases in indel accumulation and transposable element proliferation or removal. To gain insight into the relative role of these and other forces, approximately 105 kb of contiguous sequence surrounding the cellulose synthase gene CesA1 was compared for the two coresident genomes (AT and DT) of the allopolyploid cotton species, Gossypium hirsutum. These two genomes differ approximately twofold in size, having diverged from a common ancestor approximately 5-10 million years ago (Mya) and been reunited in the same nucleus at the time of polyploid formation, approximately 1-2 Mya. Gene content, order, and spacing are largely conserved between the two genomes, although a few transposable elements and a single cpDNA fragment distinguish the two homoeologs. Sequence conservation is high in both intergenic and genic regions, with 14 conserved genes detected in both genomes yielding a density of 1 gene every 7.5 kb. In contrast to the twofold overall difference in DNA content, no disparity in size was observed for this 105-kb region, and 555 indels were detected that distinguish the two homoeologous BACs, approximately equally distributed between AT and DT in number and aggregate size. The data demonstrate that genome size evolution at this phylogenetic scale is not primarily caused by mechanisms that operate uniformly across different genomic regions and components; instead, the twofold overall difference in DNA content must reflect locally operating forces between gene islands or in largely gene-free regions.

Analysis of complete nucleotide sequences of 12 Gossypium chloroplast genomes: origin and evolution of allotetraploids

BACKGROUND

Cotton (Gossypium spp.) is a model system for the analysis of polyploidization. Although ascertaining the donor species of allotetraploid cotton has been intensively studied, sequence comparison of Gossypium chloroplast genomes is still of interest to understand the mechanisms underlining the evolution of Gossypium allotetraploids, while it is generally accepted that the parents were A- and D-genome containing species. Here we performed a comparative analysis of 13 Gossypium chloroplast genomes, twelve of which are presented here for the first time.

METHODOLOGY/PRINCIPAL FINDINGS

The size of 12 chloroplast genomes under study varied from 159,959 bp to 160,433 bp. The chromosomes were highly similar having >98% sequence identity. They encoded the same set of 112 unique genes which occurred in a uniform order with only slightly different boundary junctions. Divergence due to indels as well as substitutions was examined separately for genome, coding and noncoding sequences. The genome divergence was estimated as 0.374% to 0.583% between allotetraploid species and A-genome, and 0.159% to 0.454% within allotetraploids. Forty protein-coding genes were completely identical at the protein level, and 20 intergenic sequences were completely conserved. The 9 allotetraploids shared 5 insertions and 9 deletions in whole genome, and 7-bp substitutions in protein-coding genes. The phylogenetic tree confirmed a close relationship between allotetraploids and the ancestor of A-genome, and the allotetraploids were divided into four separate groups. Progenitor allotetraploid cotton originated 0.43-0.68 million years ago (MYA).

CONCLUSION

Despite high degree of conservation between the Gossypium chloroplast genomes, sequence variations among species could still be detected. Gossypium chloroplast genomes preferred for 5-bp indels and 1-3-bp indels are mainly attributed to the SSR polymorphisms. This study supports that the common ancestor of diploid A-genome species in Gossypium is the maternal source of extant allotetraploid species and allotetraploids have a monophyletic origin. G. hirsutum AD1 lineages have experienced more sequence variations than other allotetraploids in intergenic regions. The available complete nucleotide sequences of 12 Gossypium chloroplast genomes should facilitate studies to uncover the molecular mechanisms of compartmental co-evolution and speciation of Gossypium allotetraploids.

Differential evolutionary dynamics of duplicated paralogous Adh loci in allotetraploid cotton (Gossypium)

Levels and patterns of nucleotide diversity vary widely among lineages. Because allopolyploid species contain duplicated (homoeologous) genes, studies of nucleotide diversity at homoeologous loci may facilitate insight into the evolutionary dynamics of duplicated loci. In this study, we describe patterns of sequence diversity from an alcohol dehydrogenase homoeologous locus pair (AdhC) in allotetraploid cotton (Gossypium, Malvaceae). These data are compared with equivalent information from another homoeologous alcohol dehydrogenase gene pair (AdhA, Small, Ryburn, and Wendel 1999. Mol. Biol. Evol. 16:491-501) which has an overall slower evolutionary rate than AdhC. As expected from the predicted correlation between nucleotide diversity and evolutionary rate, nucleotide diversity was higher for AdhC than for AdhA. In addition, nucleotide diversity is higher in the D-subgenome of allotetraploid cotton for AdhC, confirming earlier observations for AdhA. These observations indicate that for these two pairs of Adh loci, the null hypothesis of equivalent evolutionary dynamics for duplicated genes in allotetraploid cotton is rejected.

Ancient gene duplicates in Gossypium (cotton) exhibit near-complete expression divergence

Whole genome duplication (WGD) is widespread in flowering plants and is a driving force in angiosperm diversification. The redundancy introduced by WGD allows the evolution of novel gene interactions and functions, although the patterns and processes of diversification are poorly understood. We identified ∼ 2,000 pairs of paralogous genes in Gossypium raimondii (cotton) resulting from an approximately 60 My old 5- to 6-fold ploidy increase. Gene expression analyses revealed that, in G. raimondii, 99.4% of the gene pairs exhibit differential expression in at least one of the three tissues (petal, leaf, and seed), with 93% to 94% exhibiting differential expression on a per-tissue basis. For 1,666 (85%) pairs, differential expression was observed in all tissues. These observations were mirrored in a time series of G. raimondii seed, and separately in leaf, petal, and seed of G. arboreum, indicating expression level diversification before species divergence. A generalized linear model revealed 92.4% of the paralog pairs exhibited expression divergence, with most exhibiting significant gene and tissue interactions indicating complementary expression patterns in different tissues. These data indicate massive, near-complete expression level neo- and/or subfunctionalization among ancient gene duplicates, suggesting these processes are essential in their maintenance over ∼ 60 Ma.

Genetic diversity and population structure in the US Upland cotton (Gossypium hirsutum L.)

Genetic diversity and population structure in the US Upland cotton was established and core sets of allelic richness were identified for developing association mapping populations in cotton. Elite plant breeding programs could likely benefit from the unexploited standing genetic variation of obsolete cultivars without the yield drag typically associated with wild accessions. A set of 381 accessions comprising 378 Upland (Gossypium hirsutum L.) and 3 G. barbadense L. accessions of the United States cotton belt were genotyped using 120 genome-wide SSR markers to establish the genetic diversity and population structure in tetraploid cotton. These accessions represent more than 100 years of Upland cotton breeding in the United States. Genetic diversity analysis identified a total of 546 alleles across 141 marker loci. Twenty-two percent of the alleles in Upland accessions were unique, specific to a single accession. Population structure analysis revealed extensive admixture and identified five subgroups corresponding to Southeastern, Midsouth, Southwest, and Western zones of cotton growing areas in the United States, with the three accessions of G. barbadense forming a separate cluster. Phylogenetic analysis supported the subgroups identified by STRUCTURE. Average genetic distance between G. hirsutum accessions was 0.195 indicating low levels of genetic diversity in Upland cotton germplasm pool. The results from both population structure and phylogenetic analysis were in agreement with pedigree information, although there were a few exceptions. Further, core sets of different sizes representing different levels of allelic richness in Upland cotton were identified. Establishment of genetic diversity, population structure, and identification of core sets from this study could be useful for genetic and genomic analysis and systematic utilization of the standing genetic variation in Upland cotton.

Evolution and Stress Responses of Gossypium hirsutum SWEET Genes

The SWEET (sugars will eventually be exported transporters) proteins are sugar efflux transporters containing the MtN3_saliva domain, which affects plant development as well as responses to biotic and abiotic stresses. These proteins have not been functionally characterized in the tetraploid cotton, Gossypium hirsutum, which is a widely cultivated cotton species. In this study, we comprehensively analyzed the cotton SWEET gene family. A total of 55 putative G. hirsutumSWEET genes were identified. The GhSWEET genes were classified into four clades based on a phylogenetic analysis and on the examination of gene structural features. Moreover, chromosomal localization and an analysis of homologous genes in Gossypium arboreum, Gossypium raimondii, and G. hirsutum suggested that a whole-genome duplication, several tandem duplications, and a polyploidy event contributed to the expansion of the cotton SWEET gene family, especially in Clade III and IV. Analyses of cis-acting regulatory elements in the promoter regions, expression profiles, and artificial selection revealed that the GhSWEET genes were likely involved in cotton developmental processes and responses to diverse stresses. These findings may clarify the evolution of G. hirsutum SWEET gene family and may provide a foundation for future functional studies of SWEET proteins regarding cotton development and responses to abiotic stresses.

Discovery and identification of candidate genes from the chitinase gene family for Verticillium dahliae resistance in cotton

Verticillium dahliae, a destructive and soil-borne fungal pathogen, causes massive losses in cotton yields. However, the resistance mechanism to V. dahilae in cotton is still poorly understood. Accumulating evidence indicates that chitinases are crucial hydrolytic enzymes, which attack fungal pathogens by catalyzing the fungal cell wall degradation. As a large gene family, to date, the chitinase genes (Chis) have not been systematically analyzed and effectively utilized in cotton. Here, we identified 47, 49, 92, and 116 Chis from four sequenced cotton species, diploid Gossypium raimondii (D5), G. arboreum (A2), tetraploid G. hirsutum acc. TM-1 (AD1), and G. barbadense acc. 3-79 (AD2), respectively. The orthologous genes were not one-to-one correspondence in the diploid and tetraploid cotton species, implying changes in the number of Chis in different cotton species during the evolution of Gossypium. Phylogenetic classification indicated that these Chis could be classified into six groups, with distinguishable structural characteristics. The expression patterns of Chis indicated their various expressions in different organs and tissues, and in the V. dahliae response. Silencing of Chi23, Chi32, or Chi47 in cotton significantly impaired the resistance to V. dahliae, suggesting these genes might act as positive regulators in disease resistance to V. dahliae.

Identification of the group IIa WRKY subfamily and the functional analysis of GhWRKY17 in upland cotton (Gossypium hirsutum L.)

WRKY transcription factors play important roles in plant defense, stress response, leaf senescence, and plant growth and development. Previous studies have revealed the important roles of the group IIa GhWRKY genes in cotton. To comprehensively analyze the group IIa GhWRKY genes in upland cotton, we identified 15 candidate group IIa GhWRKY genes in the Gossypium hirsutum genome. The phylogenetic tree, intron-exon structure, motif prediction and Ka/Ks analyses indicated that most group IIa GhWRKY genes shared high similarity and conservation and underwent purifying selection during evolution. In addition, we detected the expression patterns of several group IIa GhWRKY genes in individual tissues as well as during leaf senescence using public RNA sequencing data and real-time quantitative PCR. To better understand the functions of group IIa GhWRKYs in cotton, GhWRKY17 (KF669857) was isolated from upland cotton, and its sequence alignment, promoter cis-acting elements and subcellular localization were characterized. Moreover, the over-expression of GhWRKY17 in Arabidopsis up-regulated the senescence-associated genes AtWRKY53, AtSAG12 and AtSAG13, enhancing the plant’s susceptibility to leaf senescence. These findings lay the foundation for further analysis and study of the functions of WRKY genes in cotton.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array

BACKGROUND

Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits.

RESULTS

The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content.

CONCLUSIONS

Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

Genome-wide identification and analysis of the evolution and expression patterns of the GATA transcription factors in three species of Gossypium genus

GATA transcription factors (TFs), which bind to DNA in regulatory regions, are involved in cell differentiation and possess a type-IV zinc finger and a DNA-binding domain. GATA genes have been characterized in plant species such as Arabidopsis thaliana, Oryza sativa, and Glycine max, and their functions have been elucidated in A. thaliana. Although many Gossypium quantitative trait loci for fiber quality harbor GATA TFs, GATA genes have not yet been characterized in cotton. In this study, we identified 179 GATA genes from the genomes of three Gossypium species. We analyzed the phylogenetic relationships, chromosomal distribution, gene structure, expression pattern, and predicted promoters of all 179 Gossypium GATA genes (46 in G. raimondii, 46 in G. arboreum, and 87 in G. hirsutum). Phylogenetic analysis grouped the 179 GATA genes into four subfamilies. Domain analysis revealed that GATA domains in subfamilies I, II, and III were located near the C-terminal, whereas those in subfamily IV were adjacent to the N-terminal. RNA-seq and (Real-time PCR) qRT-PCR revealed that 39.1% (34/87) of GATA genes were expressed in growing plant tissues in G. hirsutum, but only 12.6% (11/87) were expressed during fiber development. In addition, 45.7% (21/46) and 26.1% (12/46) of GATA genes were expressed in G. arboreum and G. raimondii, respectively. Our results may be useful for elucidating the evolution, expression patterns, and functional divergence of GATA genes in Gossypium.

Comprehensive characterization and gene expression patterns of LBD gene family in Gossypium

A comprehensive account of the LBD gene family of Gossypium was provided in this work. Expression analysis and functional characterization revealed that LBD genes might play different roles in G. hirsutum and G. barbadense. The Lateral Organ Boundaries Domain (LBD) proteins comprise a plant-specific transcription factor family, which plays crucial roles in physiological processes of plant growth, development, and stress tolerance. In the present work, a systematical analysis of LBD gene family from two allotetraploid cotton species, G. hirsutum and G. barbadense, together with their genomic donor species, G. arboreum and G. raimondii, was conducted. There were 131, 128, 62, and 68 LBDs identified in G. hirsutum, G. barbadense, G. arboreum and G. raimondii, respectively. The LBD proteins could be classified into two main classes, class I and class II, based on the structure of their lateral organ boundaries domain and traits of phylogenetic tree, and class I was further divided into five subgroups. The gene structure and motif composition analyses conducted in both G. hirsutum and G. barbadense revealed that LBD genes kept relatively conserved within the subfamilies. Synteny analysis suggested that segmental duplication acted as an important mechanism in expansion of the cotton LBD gene family. Cis-element analysis predicated the possible functions of LBD genes. Public RNA-seq data were investigated to analyze the expression patterns of cotton LBD genes in various tissues as well as gene expression under abiotic stress treatments. Furthermore, RT-qPCR results found that GhLBDs had various expression regulation under MeJA treatments. Expression analysis indicated the differential functions of cotton LBD genes in response to abiotic stress and hormones.

Molecular and Physio-Biochemical Characterization of Cotton Species for Assessing Drought Stress Tolerance

Drought stress significantly limits cotton growth and production due to the necessity of water at every stage of crop growth. Hence, it is essential to identify tolerant genetic resources and understand the mechanisms of drought tolerance in economically and socially important plants such as cotton. In this study, molecular and physio-biochemical investigations were conducted by analyzing different parameters by following standard protocols in three different cotton species, namely TM-1 (Gossypium hirsutum), Zhongmian-16 (Gossypium arboreum), and Pima4-S (Gossypium barbadense). Drought stress significantly decreased plant growth, chlorophyll content, net photosynthetic rate (Pn), stomatal conductance (Gs), maximum photochemical efficiency of PSII (Fv/Fm), and relative water content. TM-1 resulted in more tolerance than the other two species. The accumulation of proline, soluble proteins, soluble sugars, hydrogen peroxide (H₂O₂), and superoxide radicals (O₂•-) increased significantly in TM-1. In addition, TM-1 maintained the integrity of the chloroplast structure under drought conditions. The relative expression level of drought-responsive genes including coding for transcription factors and other regulatory proteins or enzymes controlling genes (ERF, ERFB, DREB, WRKY6, ZFP1, FeSOD, CuZnSOD, MAPKKK17, P5CR, and PRP5) were higher in TM-1 under drought, conferring a more tolerant status than in Zhongmian-16 and Pima4-S. The findings of this research could be utilized for predicting a tolerant cotton genotype as well as evaluating prospective cotton species in the variety development program.

Molecular cloning, expression profiling, and yeast complementation of 19 beta-tubulin cDNAs from developing cotton ovules

Microtubules are a major structural component of the cytoskeleton and participate in cell division, intracellular transport, and cell morphogenesis. In the present study, 795 cotton tubulin expressed sequence tags were analysed and 19 beta-tubulin genes (TUB) cloned from a cotton cDNA library. Among the group, 12 cotton TUBs (GhTUBs) are reported for the first time here. Transcription profiling revealed that nine GhTUBs were highly expressed in elongating fibre cells as compared with fuzzless-lintless mutant ovules. Treating cultured wild-type cotton ovules with exogenous phytohormones showed that individual genes can be induced by different agents. Gibberellin induced expression of GhTUB1 and GhTUB3, ethylene induced expression of GhTUB5, GhTUB9, and GhTUB12, brassinosteroids induced expression of GhTUB1, GhTUB3, GhTUB9, and GhTUB12, and lignoceric acid induced expression of GhTUB1, GhTUB3, and GhTUB12. When GhTUBs were transformed into the Saccharomyces cerevisiae inviable mutant, tub2, which is deficient in beta-tubulin, one ovule-specific and eight of nine fibre-preferential GhTUBs rescued this lethality. This study suggests that the proteins encoded by cotton GhTUBs are involved during cotton fibre development.