Sporadic Genetic Connectivity among Small Insular Populations of the Rare Geoendemic Plant Caulanthus amplexicaulis var. barbarae (Santa Barbara Jewelflower).

Globally, a small number of plants have adapted to terrestrial outcroppings of serpentine geology, which are characterized by soils with low levels of essential mineral nutrients (N, P, K, Ca, Mo) and toxic levels of heavy metals (Ni, Cr, Co). Paradoxically, many of these plants are restricted to this harsh environment. Caulanthus ampexlicaulis var. barbarae (Brassicaceae) is a rare annual plant that is strictly endemic to a small set of isolated serpentine outcrops in the coastal mountains of central California. The goals of the work presented here were to 1) determine the patterns of genetic connectivity among all known populations of C. ampexlicaulis var. barbarae, and 2) estimate contemporary effective population sizes (Ne), to inform ongoing genomic analyses of the evolutionary history of this taxon, and to provide a foundation upon which to model its future evolutionary potential and long-term viability in a changing environment. Eleven populations of this taxon were sampled, and population-genetic parameters were estimated using 11 nuclear microsatellite markers. Contemporary effective population sizes were estimated using multiple methods and found to be strikingly small (typically Ne < 10). Further, our data showed that a substantial component of genetic connectivity of this taxon is not at equilibrium, and instead showed sporadic gene flow. Several lines of evidence indicate that gene flow between isolated populations is maintained through long-distance seed dispersal (e.g., >1 km), possibly via zoochory.

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.

Development, genetic mapping and QTL association of cotton PHYA, PHYB, and HY5-specific CAPS and dCAPS markers.

BACKGROUND: Among SNP markers that become increasingly valuable in molecular breeding of crop plants are the CAPS and dCAPS markers derived from the genes of interest. To date, the number of such gene-based markers is small in polyploid crop plants such as allotetraploid cotton that has A- and D-sub-genomes. The objective of this study was to develop and map new CAPS and dCAPS markers for cotton developmental-regulatory genes that are important in plant breeding programs. RESULTS: Gossypium hirsutum and G. barbadense, are the two cultivated allotetraploid cotton species. These have distinct fiber quality and other agronomic traits. Using comparative sequence analysis of characterized GSTs of the PHYA1, PHYB, and HY5 genes of G. hirsutum and G. barbadense one PHYA1-specific Mbo I/Dpn II CAPS, one PHYB-specific Alu I dCAPS, and one HY5-specific Hinf I dCAPS cotton markers were developed. These markers have successfully differentiated the two allotetraploid genomes (AD1 and AD2) when tested in parental genotypes of 'Texas Marker-1' ('TM-1'), 'Pima 3-79' and their F1 hybrids. The genetic mapping and chromosome substitution line-based deletion analyses revealed that PHYA1 gene is located in A-sub-genome chromosome 11, PHYB gene is in A-sub-genome chromosome 10, and HY5 gene is in D-sub-genome chromosome 24, on the reference 'TM-1' x 'Pima 3-79' RIL genetic map. Further, it was found that genetic linkage map regions containing phytochrome and HY5-specific markers were associated with major fiber quality and flowering time traits in previously published QTL mapping studies. CONCLUSION: This study detailed the genome mapping of three cotton phytochrome genes with newly developed CAPS and dCAPS markers. The proximity of these loci to fiber quality and other cotton QTL was demonstrated in two A-subgenome and one D-subgenome chromosomes. These candidate gene markers will be valuable for marker-assisted selection (MAS) programs to rapidly introgress G. barbadense phytochromes and/or HY5 gene (s) into G. hirsutum cotton genotypes or vice versa.

Exploring origins, invasion history and genetic diversity of Imperata cylindrica (L.) P. Beauv. (Cogongrass) in the United States using genotyping by sequencing.

Imperata cylindrica (Cogongrass, Speargrass) is a diploid C4 grass that is a noxious weed in 73 countries and constitutes a significant threat to global biodiversity and sustainable agriculture. We used a cost-effective genotyping-by-sequencing (GBS) approach to identify the reproductive system, genetic diversity and geographic origins of invasions in the south-eastern United States. In this work, we demonstrated the advantage of employing the closely related, fully sequenced crop species Sorghum bicolor (L.) Moench as a proxy reference genome to identify a set of 2320 informative single nucleotide and insertion-deletion polymorphisms. Genetic analyses identified four clonal lineages of cogongrass and one clonal lineage of Imperata brasiliensis Trin. in the United States. Each lineage was highly homogeneous, and we found no evidence of hybridization among the different lineages, despite geographical overlap. We found evidence that at least three of these lineages showed clonal reproduction prior to introduction to the United States. These results indicate that cogongrass has limited evolutionary potential to adapt to novel environments and further suggest that upon arrival to its invaded range, this species did not require local adaptation through hybridization/introgression or selection of favourable alleles from a broad genetic base. Thus, cogongrass presents a clear case of broad invasive success, across a diversity of environments, in a clonal organism with limited genetic diversity.

De novo assembly and characterization of the transcriptome of the toxic dinoflagellate Karenia brevis.

BACKGROUND: Karenia brevis is a harmful algal species that blooms in the Gulf of Mexico and produces brevetoxins that cause neurotoxic shellfish poisoning. Elevated brevetoxin levels in K. brevis cells have been measured during laboratory hypo-osmotic stress treatments. To investigate mechanisms underlying K. brevis osmoacclimation and osmoregulation and establish a valuable resource for gene discovery, we assembled reference transcriptomes for three clones: Wilson-CCFWC268, SP3, and SP1 (a low-toxin producing variant). K. brevis transcriptomes were annotated with gene ontology terms and searched for putative transmembrane proteins that may elucidate cellular responses to hypo-osmotic stress. An analysis of single nucleotide polymorphisms among clones was used to characterize genetic divergence. RESULTS: K. brevis reference transcriptomes were assembled with 58.5 (Wilson), 78.0 (SP1), and 51.4 million (SP3) paired reads. Transcriptomes contained 86,580 (Wilson), 93,668 (SP1), and 84,309 (SP3) predicted transcripts. Approximately 40% of the transcripts were homologous to proteins in the BLAST nr database with an E value ≤ 1.00E-6. Greater than 80% of the highly conserved CEGMA core eukaryotic genes were identified in each transcriptome, which supports assembly completeness. Seven putative voltage-gated Na+ or Ca2+ channels, two aquaporin-like proteins, and twelve putative VATPase subunits were discovered in all clones using multiple bioinformatics approaches. Furthermore, 45% (Wilson) and 43% (SP1 and SP3) of the K. brevis putative peptides > 100 amino acids long produced significant hits to a sequence in the NCBI nr protein database. Of these, 77% (Wilson and SP1) and 73% (SP3) were successfully assigned gene ontology functional terms. The predicted single nucleotide polymorphism (SNP) frequencies between clones were 0.0028 (Wilson to SP1), 0.0030 (Wilson to SP3), and 0.0028 (SP1 to SP3). CONCLUSIONS: The K. brevis transcriptomes assembled here provide a foundational resource for gene discovery and future RNA-seq experiments. The identification of ion channels, VATPases, and aquaporins in all three transcriptomes indicates that K. brevis regulates cellular ion and water concentrations via transmembrane proteins. Additionally, > 40,000 unannotated loci may include potentially novel K. brevis genes. Ultimately, the SNPs identified among the three ecologically diverse clones with different toxin profiles may help to elucidate variations in K. brevis brevetoxin production.

Mapping genomic loci for cotton plant architecture, yield components, and fiber properties in an interspecific (Gossypium hirsutum L. × G. barbadense L.) RIL population.

A quantitative trait locus (QTL) mapping was conducted to better understand the genetic control of plant architecture (PA), yield components (YC), and fiber properties (FP) in the two cultivated tetraploid species of cotton (Gossypium hirsutum L. and G. barbadense L.). One hundred and fifty-nine genomic regions were identified on a saturated genetic map of more than 2,500 SSR and SNP markers, constructed with an interspecific recombinant inbred line (RIL) population derived from the genetic standards of the respective cotton species (G. hirsutum acc. TM-1 × G. barbadense acc. 3-79). Using the single nonparametric and MQM QTL model mapping procedures, we detected 428 putative loci in the 159 genomic regions that confer 24 cotton traits in three diverse production environments [College Station F&B Road (FB), TX; Brazos Bottom (BB), TX; and Shafter (SH), CA]. These putative QTL loci included 25 loci for PA, 60 for YC, and 343 for FP, of which 3, 12, and 60, respectively, were strongly associated with the traits (LOD score ≥ 3.0). Approximately 17.7 % of the PA putative QTL, 32.9 % of the YC QTL, and 48.3 % of the FP QTL had trait associations under multiple environments. The At subgenome (chromosomes 1-13) contributed 72.7 % of loci for PA, 46.2 % for YC, and 50.4 % for FP while the Dt subgenome (chromosomes 14-26) contributed 27.3 % of loci for PA, 53.8 % for YC, and 49.6 % for FP. The data obtained from this study augment prior evidence of QTL clusters or gene islands for specific traits or biological functions existing in several non-homoeologous cotton chromosomes. DNA markers identified in the 159 genomic regions will facilitate further dissection of genetic factors underlying these important traits and marker-assisted selection in cotton.

Novel nuclear markers inform the systematics and the evolution of serpentine use in Streptanthus and allies (Thelypodieae, Brassicaceae).

Streptanthus is a genus of ca. 35 species in the tribe Thelypodieae (Brassicaceae) that has remarkable morphological and ecological diversity, a large number of species in the group being edaphic specialists endemic to unusual soils such as serpentine. While ecological research has shed some light on adaptation to serpentine in Streptanthus, there have been few insights on the origins and evolution of serpentine tolerance in this group, largely due to limited success in resolving the phylogenetic relationships among Streptanthus and allied genera of the Thelypodieae (Streptanthoid complex). We present a well-resolved phylogenetic hypothesis for the Streptanthoid complex, based on three newly identified and highly variable single copy nuclear regions (AT4G34700, AT1G61620, and AT1G56590, and three others that are widely used (ITS, phyA, and PEPC). We also include data for two chloroplast regions (trnL and trnH-psbA). Collectively, our new markers provide 75% of the nuclear parsimony informative characters in our data. Taxonomically, our sampling is the most inclusive of any study of the Streptanthoid Complex to date, including 46 out of the 53 species of Streptanthus and Caulanthus, as well as representatives of several closely allied genera in the Thelypodieae. Our results reveal that Streptanthus, Caulanthus, and Thelypodium are not reciprocally monophyletic as currently defined. The species of Streptanthus form two rather distantly related clades. One clade (SC-I) is comprised of species with bilateral flowers and urn-shaped calyces that occur mainly within the California Floristic Province (CFP) hotspot; the other clade (SC-II) is composed of species with extant ranges mainly outside the CFP. Our data indicate that serpentine tolerance has evolved between eight and ten times in this group, of which between four and five have resulted in endemism. While serpentine endemism has been rarely lost, large and diverse clades composed mainly of serpentine endemics indicate that serpentine endemics in this group are more than mere 'dead-ends'.

Genetic analyses of nickel tolerance in a North American serpentine endemic plant, Caulanthus amplexicaulis var. barbarae (Brassicaceae).

PREMISE OF THE STUDY: The evolution of metal tolerance in plants is an important model for studies of adaptation to environment, population genetics, and speciation. Here, we investigated nickel tolerance in the North American serpentine endemic Caulanthus amplexicaulis var. barbarae in comparison with its nonserpentine sister taxon C. amplexicaulis var. amplexicaulis. We hypothesized that the serpentine endemic would have a heritable growth advantage on nickel-containing substrates. METHODS: We employed an artificial growth assay to quantify biomass accumulation. Study plants were crossed to create an F(2:3) population that was used to determine the heritability of nickel tolerance and to map quantitative trait loci (QTL). Nickel accumulation in both laboratory populations and native specimens was examined using energy-dispersive x-ray fluorescence (EDXRF). KEY RESULTS: The serpentine endemic had a dramatic growth advantage at concentrations of nickel >30 µmol/L. Caulanthus amplexicaulis var. barbarae and its nonserpentine sister taxon both accumulated nickel to substantial levels. Nickel tolerance was highly heritable (h(2) = 0.59) and not associated with accumulation. The QTL analyses identified two major loci for nickel tolerance, on linkage group 2 (LG2) and linkage group 9 (LG9). CONCLUSIONS: In our study, nickel tolerance was determined by two major loci with large effects. At both loci, alleles from the serpentine parent conferred positive effects on nickel tolerance, suggesting that they are adaptive in the natural serpentine environment. The mechanism of nickel tolerance in the serpentine plant was not exclusion of nickel. Nickel tolerance may have an inducible component in C. amplexicaulis var. barbarae.

A comparative genomic map for Caulanthus amplexicaulis and related species (Brassicaceae).

Adaptation to environment is the cornerstone of ecological genetics. The subject of this study is a wild relative of the sequenced and annotated model plant species, Arabidopsis thaliana. Caulanthus amplexicaulis var. barbarae lives on serpentine soils, known for high concentrations of heavy metals and low concentrations of essential plant macronutrients, and provides a compelling example of an organism's adaptation to environment. We constructed an F(2) linkage map, using a cross to the nonserpentine sister taxon, C. amplexicaulis var. amplexicaulis. C. amplexicaulis is a member of a highly diverse set of taxa (within the tribe Thelypodieae), described here as the 'Streptanthoid Complex' that are adapted to a broad range of environments, yet share a common n = 14 chromosome number and likely arose by a recent radiation. The linkage map consists of 97 polymorphic microsatellite markers, and 40 exon-primed intron-crossing markers based on A. thaliana exon sequences and Brassica ESTs. The map covers 14 linkage groups and has a total length of 1513 cM. Both the patterns of marker segregation and the comparative map indicate that C. amplexicaulis is a diploid organism with a compact genome. All exon-primed intron-crossing markers, and an unexpectedly large number of microsatellite markers (83%), had significant similarity to the A. thaliana genome, facilitating the development of a comparative genome map. As a proof of principle, we used the comparative map to identify candidate genes underlying differences in sepal colour between the two parent taxa. We demonstrate that the genomic tools developed here will be portable throughout the Streptanthoid Complex.

Isolation and characterization of microsatellite markers for Bothriochloa ischaemum (Poaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed for Bothriochloa ischaemum to investigate the structure of invasive populations within Texas and determine the origin of introduction from within the native range. METHODS AND RESULTS: We adapted the biotinylated nucleotide method of developing an enriched genomic library to isolate and characterize 10 polymorphic microsatellite markers. The number of alleles per locus (A) ranged from five to 18 (mean A = 10.45), expected heterozygosity (H(E)) ranged from 0.620 to 0.895 (mean H(E) = 0.785), and expected heterozygosity corrected for sample size (H(EC)) ranged from 0.635 to 0.909 (mean H(EC) = 0.799). The primers were also tested for amplification in Schizachyrium scoparium var. scoparium, Andropogon gerardii, Bothriochloa saccharoides, and Dichanthium annulatum. CONCLUSIONS: The use of microsatellite markers may assist in understanding the pattern of spread, determining the source of invasive populations, and developing biological control agents for invasive populations of Bothriochloa ischaemum.

Duplication, divergence and persistence in the Phytochrome photoreceptor gene family of cottons (Gossypium spp.).

BACKGROUND: Phytochromes are a family of red/far-red photoreceptors that regulate a number of important developmental traits in cotton (Gossypium spp.), including plant architecture, fiber development, and photoperiodic flowering. Little is known about the composition and evolution of the phytochrome gene family in diploid (G. herbaceum, G. raimondii) or allotetraploid (G. hirsutum, G. barbadense) cotton species. The objective of this study was to obtain a preliminary inventory and molecular-evolutionary characterization of the phytochrome gene family in cotton. RESULTS: We used comparative sequence resources to design low-degeneracy PCR primers that amplify genomic sequence tags (GSTs) for members of the PHYA, PHYB/D, PHYC and PHYE gene sub-families from A- and D-genome diploid and AD-genome allotetraploid Gossypium species. We identified two paralogous PHYA genes (designated PHYA1 and PHYA2) in diploid cottons, the result of a Malvaceae-specific PHYA gene duplication that occurred approximately 14 million years ago (MYA), before the divergence of the A- and D-genome ancestors. We identified a single gene copy of PHYB, PHYC, and PHYE in diploid cottons. The allotetraploid genomes have largely retained the complete gene complements inherited from both of the diploid genome ancestors, with at least four PHYA genes and two genes encoding PHYB, PHYC and PHYE in the AD-genomes. We did not identify a PHYD gene in any cotton genomes examined. CONCLUSIONS: Detailed sequence analysis suggests that phytochrome genes retained after duplication by segmental duplication and allopolyploidy appear to be evolving independently under a birth-and-death-process with strong purifying selection. Our study provides a preliminary phytochrome gene inventory that is necessary and sufficient for further characterization of the biological functions of each of the cotton phytochrome genes, and for the development of 'candidate gene' markers that are potentially useful for cotton improvement via modern marker-assisted selection strategies.

Linkage disequilibrium based association mapping of fiber quality traits in G. hirsutum L. variety germplasm.

Cotton is the world's leading cash crop, but it lags behind other major crops for marker-assisted breeding due to limited polymorphisms and a genetic bottleneck through historic domestication. This underlies a need for characterization, tagging, and utilization of existing natural polymorphisms in cotton germplasm collections. Here we report genetic diversity, population characteristics, the extent of linkage disequilibrium (LD), and association mapping of fiber quality traits using 202 microsatellite marker primer pairs in 335 G. hirsutum germplasm grown in two diverse environments, Uzbekistan and Mexico. At the significance threshold (r (2) >or= 0.1), a genome-wide average of LD extended up to genetic distance of 25 cM in assayed cotton variety accessions. Genome wide LD at r (2) >or= 0.2 was reduced to approximately 5-6 cM, providing evidence of the potential for association mapping of agronomically important traits in cotton. Results suggest linkage, selection, inbreeding, population stratification, and genetic drift as the potential LD-generating factors in cotton. In two environments, an average of ~20 SSR markers was associated with each main fiber quality traits using a unified mixed liner model (MLM) incorporating population structure and kinship. These MLM-derived significant associations were confirmed in general linear model and structured association test, accounting for population structure and permutation-based multiple testing. Several common markers, showing the significant associations in both Uzbekistan and Mexican environments, were determined. Between 7 and 43% of the MLM-derived significant associations were supported by a minimum Bayes factor at 'moderate to strong' and 'strong to very strong' evidence levels, suggesting their usefulness for marker-assisted breeding programs and overall effectiveness of association mapping using cotton germplasm resources.

Development of 15 microsatellite loci in the endangered Streptanthus glandulosus subsp. niger (Brassicaceae).

PREMISE OF THE STUDY: The endangered Streptanthus glandulosus subsp. niger (Brassicaceae) is endemic to a single peninsula in California and threatened by fragmentation. We developed microsatellite markers to investigate genetic diversity in the two extant populations and the degree to which they have diverged from one another. METHODS AND RESULTS: We used Illumina HiSeq high-throughput sequencing to develop 15 microsatellite markers, 14 of which were polymorphic. These di- and trinucleotide repeats yielded one to 11 alleles per locus in 61 plants across the two populations. Levels of observed and expected heterozygosities ranged from 0.108 to 0.946 and 0.257 to 0.839, respectively. We demonstrated cross-amplification in a second rare subspecies, S. glandulosus subsp. secundus, and in the widespread congener S. tortuosus. CONCLUSIONS: These are the first microsatellites reported for this subspecies, and they will aid in the inclusion of genetic information in conservation planning. Cross-amplification was demonstrated in two related taxa, including one of conservation concern.

Transcriptome Signatures of Selection, Drift, Introgression, and Gene Duplication in the Evolution of an Extremophile Endemic Plant.

Plants on serpentine soils provide extreme examples of adaptation to environment, and thus offer excellent models for the study of evolution at the molecular and genomic level. Serpentine outcrops are derived from ultramafic rock and have extremely low levels of essential plant nutrients (e.g., N, P, K, and Ca), as well as toxic levels of heavy metals (e.g., Ni, Cr, and Co) and low moisture availability. These outcrops provide habitat to a number of endemic plant species, including the annual mustard Caulanthus amplexicaulis var. barbarae (Cab) (Brassicaceae). Its sister taxon, C. amplexicaulis var. amplexicaulis (Caa), is intolerant to serpentine soils. Here, we assembled and annotated comprehensive reference transcriptomes of both Caa and Cab for use in protein coding sequence comparisons. A set of 29,443 reciprocal best Blast hit (RBH) orthologs between Caa and Cab was compared with identify coding sequence variants, revealing a high genome-wide dN/dS ratio between the two taxa (mean = 0.346). We show that elevated dN/dS likely results from the composite effects of genetic drift, positive selection, and the relaxation of negative selection. Further, analysis of paralogs within each taxon revealed the signature of a period of elevated gene duplication (∼10 Ma) that is shared with other species of the tribe Thelypodieae, and may have played a role in the striking morphological and ecological diversity of this tribe. In addition, distribution of the synonymous substitution rate, dS, is strongly bimodal, indicating a history of reticulate evolution that may have contributed to serpentine adaptation.

QTL mapping for flowering-time and photoperiod insensitivity of cotton Gossypium darwinii Watt.

Most wild and semi-wild species of the genus Gossypium are exhibit photoperiod-sensitive flowering. The wild germplasm cotton is a valuable source of genes for genetic improvement of modern cotton cultivars. A bi-parental cotton population segregating for photoperiodic flowering was developed by crossing a photoperiod insensitive irradiation mutant line with its pre-mutagenesis photoperiodic wild-type G. darwinii Watt genotype. Individuals from the F2 and F3 generations were grown with their parental lines and F1 hybrid progeny in the long day and short night summer condition (natural day-length) of Uzbekistan to evaluate photoperiod sensitivity, i.e., flowering-time during the seasons 2008-2009. Through genotyping the individuals of this bi-parental population segregating for flowering-time, linkage maps were constructed using 212 simple-sequence repeat (SSR) and three cleaved amplified polymorphic sequence (CAPS) markers. Six QTLs directly associated with flowering-time and photoperiodic flowering were discovered in the F2 population, whereas eight QTLs were identified in the F3 population. Two QTLs controlling photoperiodic flowering and duration of flowering were common in both populations. In silico annotations of the flanking DNA sequences of mapped SSRs from sequenced cotton (G. hirsutum L.) genome database has identified several potential 'candidate' genes that are known to be associated with regulation of flowering characteristics of plants. The outcome of this research will expand our understanding of the genetic and molecular mechanisms of photoperiodic flowering. Identified markers should be useful for marker-assisted selection in cotton breeding to improve early flowering characteristics.

RNA Interference for Functional Genomics and Improvement of Cotton (Gossypium sp.).

RNA interference (RNAi), is a powerful new technology in the discovery of genetic sequence functions, and has become a valuable tool for functional genomics of cotton (Gossypium sp.). The rapid adoption of RNAi has replaced previous antisense technology. RNAi has aided in the discovery of function and biological roles of many key cotton genes involved in fiber development, fertility and somatic embryogenesis, resistance to important biotic and abiotic stresses, and oil and seed quality improvements as well as the key agronomic traits including yield and maturity. Here, we have comparatively reviewed seminal research efforts in previously used antisense approaches and currently applied breakthrough RNAi studies in cotton, analyzing developed RNAi methodologies, achievements, limitations, and future needs in functional characterizations of cotton genes. We also highlighted needed efforts in the development of RNAi-based cotton cultivars, and their safety and risk assessment, small and large-scale field trials, and commercialization.

SNP discovery in complex allotetraploid genomes (Gossypium spp., Malvaceae) using genotyping by sequencing.

PREMISE OF THE STUDY: Single-nucleotide polymorphism (SNP) marker discovery in plants with complex allotetraploid genomes is often confounded by the presence of homeologous loci (along with paralogous and orthologous loci). Here we present a strategy to filter for SNPs representing orthologous loci. METHODS AND RESULTS: Using Illumina next-generation sequencing, 54 million reads were collected from restriction enzyme-digested DNA libraries of a diversity of Gossypium taxa. Loci with one to three SNPs were discovered using the Stacks software package, yielding 25,529 new cotton SNP combinations, including those that are polymorphic at both interspecific and intraspecific levels. Frequencies of predicted dual-homozygous (aa/bb) marker polymorphisms ranged from 6.7-11.6% of total shared fragments in intraspecific comparisons and from 15.0-16.4% in interspecific comparisons. CONCLUSIONS: This resource provides dual-homozygous (aa/bb) marker polymorphisms. Both in silico and experimental validation efforts demonstrated that these markers are enriched for single orthologous loci that are homozygous for alternative alleles.

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.

Phytochrome RNAi enhances major fibre quality and agronomic traits of the cotton Gossypium hirsutum L.

Simultaneous improvement of fibre quality, early-flowering, early-maturity and productivity in Upland cotton (G. hirsutum) is a challenging task for conventional breeding. The influence of red/far-red light ratio on the fibre length prompted us to examine the phenotypic effects of RNA interference (RNAi) of the cotton PHYA1 gene. Here we show a suppression of up to ~70% for the PHYA1 transcript, and compensatory overexpression of up to ~20-fold in the remaining phytochromes in somatically regenerated PHYA1 RNAi cotton plants. Two independent transformants of three generations exhibited vigorous root and vegetative growth, early-flowering, significantly improved upper half mean fibre length and an improvement in other major fibre characteristics. Small decreases in lint traits were observed but seed cotton yield was increased an average 10-17% compared with controls. RNAi-associated phenotypes were heritable and transferable via sexual hybridization. These results should aid in the development of early-maturing and productive Upland cultivars with superior fibre quality.

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.