Thea A. Wilkins

A 3347-Locus Genetic Recombination Map of Sequence-Tagged Sites Reveals Features of Genome Organization, Transmission and Evolution of Cotton (Gossypium)

We report genetic maps for diploid (D) and tetraploid (AtDt) Gossypium genomes composed of sequence-tagged sites (STS) that foster structural, functional, and evolutionary genomic studies. The maps include, respectively, 2584 loci at 1.72-cM (∼600 kb) intervals based on 2007 probes (AtDt) and 763 loci at 1.96-cM (∼500 kb) intervals detected by 662 probes (D). Both diploid and tetraploid cottons exhibit negative crossover interference; i.e., double recombinants are unexpectedly abundant. We found no major structural changes between Dt and D chromosomes, but confirmed two reciprocal translocations between At chromosomes and several inversions. Concentrations of probes in corresponding regions of the various genomes may represent centromeres, while genome-specific concentrations may represent heterochromatin. Locus duplication patterns reveal all 13 expected homeologous chromosome sets and lend new support to the possibility that a more ancient polyploidization event may have predated the A-D divergence of 6–11 million years ago. Identification of SSRs within 312 RFLP sequences plus direct mapping of 124 SSRs and exploration for CAPS and SNPs illustrate the “portability” of these STS loci across populations and detection systems useful for marker-assisted improvement of the worlds leading fiber crop. These data provide new insights into polyploid evolution and represent a foundation for assembly of a finished sequence of the cotton genome.

Functional genomics of cell elongation in developing cotton fibers.

Cotton fibers are single-celled seed trichomes of major economic importance. Factors that regulate the rate and duration of cell expansion control fiber morphology and important agronomic traits. For genetic characterization of rapid cell elongation in cotton fibers, ∼ 14,000 unique genes were assembled from 46,603 expressed sequence tags (ESTs) from developmentally staged fiber cDNAs of a cultivated diploid species (Gossypium arboreumL.). Conservatively, the fiber transcriptome represents 35–40% of the genes in the cotton genome. In silico expression analysis revealed that rapidly elongating fiber cells exhibit significant metabolic activity, with the bulk of gene transcripts, represented by three major functional groups – cell wall structure and biogenesis, the cytoskeleton and energy/carbohydrate metabolism. Oligonucleotide microarrays revealed dynamic changes in gene expression between primary and secondary cell wall biogenesis showing that fiber genes in the dbEST are highly stage-specific for cell expansion – a conclusion supported by the absence of known secondary cell wall-specific genes from our fiber dbEST. During the developmental switch from primary to secondary cell wall syntheses, 2553 “expansion-associated” fiber genes are significantly down regulated. Genes (81) significantly up-regulated during secondary cell wall synthesis are involved in cell wall biogenesis and energy/carbohydrate metabolism, which is consistent with the stage of cellulose synthesis during secondary cell wall modification in developing fibers. This work provides the first in-depth view of the genetic complexity of the transcriptome of an expanding cell, and lays the groundwork for studying fundamental biological processes in plant biology with applications in agricultural biotechnology.

Toward Sequencing Cotton (Gossypium) Genomes

Despite rapidly decreasing costs and innovative technologies, sequencing of angiosperm genomes is not yet undertaken lightly. Generating larger amounts of sequence data more quickly does not address the difficulties of sequencing and assembling complex genomes de novo. The cotton ( Gossypium spp.)

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.

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

Abstract 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 Gl1 (AtMYBGl1) 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.

Genetic mapping and comparative analysis of seven mutants related to seed fiber development in cotton

Mapping of genes that play major roles in cotton fiber development is an important step toward their cloning and manipulation, and provides a test of their relationships (if any) to agriculturally-important QTLs. Seven previously identified fiber mutants, four dominant (Li1, Li2, N1 and Fbl) and three recessive (n2, sma-4(ha), and sma-4(fz)), were genetically mapped in six F2 populations comprising 124 or more plants each. For those mutants previously assigned to chromosomes by using aneuploids or by linkage to other morphological markers, all map locations were concordant except n2, which mapped to the homoeolog of the chromosome previously reported. Three mutations with primary effects on fuzz fibers (N1, Fbl, n2) mapped near the likelihood peaks for QTLs that affected lint fiber productivity in the same populations, perhaps suggesting pleiotropic effects on both fiber types. However, only Li1 mapped within the likelihood interval for 191 previously detected lint fiber QTLs discovered in non-mutant crosses, suggesting that these mutations may occur in genes that played early roles in cotton fiber evolution, and for which new allelic variants are quickly eliminated from improved germplasm. A close positional association between sma-4(ha), two leaf and stem-borne trichome mutants (t1, t2), and a gene previously implicated in fiber development, sucrose synthase, raises questions about the possibility that these genes may be functionally related. Increasing knowledge of the correspondence of the cotton and Arabidopsis genomes provides several avenues by which genetic dissection of cotton fiber development may be accelerated.

Sampling nucleotide diversity in cotton

BackgroundCultivated cotton is an annual fiber crop derived mainly from two perennial species, Gossypium hirsutum L. or upland cotton, and G. barbadense L., extra long-staple fiber Pima or Egyptian cotton. These two cultivated species are among five allotetraploid species presumably derived monophyletically between G. arboreum and G. raimondii. Genomic-based approaches have been hindered by the limited variation within species. Yet, population-based methods are being used for genome-wide introgression of novel alleles from G. mustelinum and G. tomentosum into G. hirsutum using combinations of backcrossing, selfing, and inter-mating. Recombinant inbred line populations between genetics standards TM-1, (G. hirsutum) × 3-79 (G. barbadense) have been developed to allow high-density genetic mapping of traits.ResultsThis paper describes a strategy to efficiently characterize genomic variation (SNPs and indels) within and among cotton species. Over 1000 SNPs from 270 loci and 279 indels from 92 loci segregating in G. hirsutum and G. barbadense were genotyped across a standard panel of 24 lines, 16 of which are elite cotton breeding lines and 8 mapping parents of populations from six cotton species. Over 200 loci were genetically mapped in a core mapping population derived from TM-1 and 3-79 and in G. hirsutum breeding germplasm.ConclusionIn this research, SNP and indel diversity is characterized for 270 single-copy polymorphic loci in cotton. A strategy for SNP discovery is defined to pre-screen loci for copy number and polymorphism. Our data indicate that the A and D genomes in both diploid and tetraploid cotton remain distinct from each such that paralogs can be distinguished. This research provides mapped DNA markers for intra-specific crosses and introgression of exotic germplasm in cotton.

Expression of Two Related Vacuolar H+-ATPase 16-Kilodalton Proteolipid Genes Is Differentially Regulated in a Tissue-Specific Manner

The 16-kD proteolipid subunit is the principal integral membrane protein of the vaculor H+-ATPase (V-ATPase) complex that forms the proton channel responsible for translocating protons across lipid bilayers. Two degenerate synthetic oligonucleotides, COT11 and COT12, corresponding to highly conserved transmembrane domains in all 16-kD subunits sequenced so far, were used to amplify a partial cDNA of the V-ATPase proteolipid subunit from cotton (Gossypium hirsutum L.) by polymerase chain reaction (PCR). These PCR products were used to isolate two full-length cDNAs from a -3 d postanthesis cotton ovule library. Both clones, CVA16.2 and CVA16.4, consisting of 816 and 895 bp, respectively, encode the 16-kD proteolipid subunit of the V-ATPase. At the nucleotide level, the complete sequences of the two clones show 73.5% identity, but share about 95% identity within the coding region, although the two polypeptides differ by only one amino acid. Comparison of deduced amino acid sequences of the proteolipid subunits revealed that the four transmembrane domains and the two cytosolic extramembrane domains are highly conserved in all eukaryotes. Southern blot analysis of cotton genomic DNA showed that these clones belong to small gene families in related diploid and allotetraploid species. Northern blot analysis suggested that the three major V-ATPase subunits (69, 60, and 16 kD) are coordinately regulated, in part, at the transcriptional level. RNA analysis and reverse-transcription PCR established that 16-kD proteolipid transcripts differentially accumulate in different tissues and increase dramatically in tissues undergoing rapid expansion, particularly in anthers, ovules, and petals. The CVA16.4 proteolipid transcript is the most prevalent of the two proteolipid messages in expanding ovules harvested 10 d postanthesis. In contrast, the two proteolipid mRNAs accumulate to similar levels in developing petals.

Identification of salt responsive genes using comparative microarray analysis in Upland cotton (Gossypium hirsutum L.)

Salinity negatively impacts plant growth and productivity, and little is known about salt responsive genes in cotton. In this study, an intra-specific backcross population of cotton (Gossypium hirsutum L.) was treated with 200 mM NaCl after which differentially expressed genes were identified by comparison between salt tolerant and susceptible segregant bulks using comparative microarray analysis. Microarray analysis identified 720 salt-responsive genes, of which 695 were down-regulated and only 25 were up-regulated in the salt tolerant bulk. Gene ontology of annotated genes revealed that at least some of the identified salt responsive transcripts belong to pathways known to be associated with salt stress including osmolyte and lipid metabolism, cell wall structure, and membrane synthesis. About 48% of all salt-responsive genes were functionally unknown. Quantitative RT-PCR was used to validate 17 selected salt responsive genes. This work represents the first study in employing microarray to investigate the possible mechanisms of the salt response in cotton. Further analysis of salt-responsive genes associated with salt tolerance in cotton will assist in laying a foundation for molecular manipulation in development of new cultivars with improved salt tolerance.

Gene-rich islands for fiber development in the cotton genome

Cotton fiber is an economically important seed trichome and the worlds leading natural fiber used in the manufacture of textiles. As a step toward elucidating the genomic organization and distribution of gene networks responsible for cotton fiber development, we investigated the distribution of fiber genes in the cotton genome. Results revealed the presence of gene-rich islands for fiber genes with a biased distribution in the tetraploid cotton (Gossypium hirsutum L.) genome that was also linked to discrete fiber developmental stages based on expression profiles. There were 3 fiber gene-rich islands associated with fiber initiation on chromosome 5, 3 islands for the early to middle elongation stage on chromosome 10, 3 islands for the middle to late elongation stage on chromosome 14, and 1 island on chromosome 15 for secondary cell wall deposition, for a total of 10 fiber gene-rich islands. Clustering of functionally related gene clusters in the cotton genome displaying similar transcriptional regulation indicates an organizational hierarchy with significant implications for the genetic enhancement of particular fiber quality traits. The relationship between gene-island distribution and functional expression profiling suggests for the first time the existence of functional coupling gene clusters in the cotton genome.