Organellar transcripts dominate the cellular mRNA pool across plants of varying ploidy levels.

Mitochondrial and plastid functions depend on coordinated expression of proteins encoded by genomic compartments that have radical differences in copy number of organellar and nuclear genomes. In polyploids, doubling of the nuclear genome may add challenges to maintaining balanced expression of proteins involved in cytonuclear interactions. Here, we use ribo-depleted RNA sequencing (RNA-seq) to analyze transcript abundance for nuclear and organellar genomes in leaf tissue from four different polyploid angiosperms and their close diploid relatives. We find that even though plastid genomes contain <1% of the number of genes in the nuclear genome, they generate the majority (69.9 to 82.3%) of messenger RNA (mRNA) transcripts in the cell. Mitochondrial genes are responsible for a much smaller percentage (1.3 to 3.7%) of the leaf mRNA pool but still produce much higher transcript abundances per gene compared to nuclear genome. Nuclear genes encoding proteins that functionally interact with mitochondrial or plastid gene products exhibit mRNA expression levels that are consistently more than 10-fold lower than their organellar counterparts, indicating an extreme cytonuclear imbalance at the RNA level despite the predominance of equimolar interactions at the protein level. Nevertheless, interacting nuclear and organellar genes show strongly correlated transcript abundances across functional categories, suggesting that the observed mRNA stoichiometric imbalance does not preclude coordination of cytonuclear expression. Finally, we show that nuclear genome doubling does not alter the cytonuclear expression ratios observed in diploid relatives in consistent or systematic ways, indicating that successful polyploid plants are able to compensate for cytonuclear perturbations associated with nuclear genome doubling.

Gene disruption by structural mutations drives selection in US rice breeding over the last century.

The genetic basis of general plant vigor is of major interest to food producers, yet the trait is recalcitrant to genetic mapping because of the number of loci involved, their small effects, and linkage. Observations of heterosis in many crops suggests that recessive, malfunctioning versions of genes are a major cause of poor performance, yet we have little information on the mutational spectrum underlying these disruptions. To address this question, we generated a long-read assembly of a tropical japonica rice (Oryza sativa) variety, Carolina Gold, which allowed us to identify structural mutations (>50 bp) and orient them with respect to their ancestral state using the outgroup, Oryza glaberrima. Supporting prior work, we find substantial genome expansion in the sativa branch. While transposable elements (TEs) account for the largest share of size variation, the majority of events are not directly TE-mediated. Tandem duplications are the most common source of insertions and are highly enriched among 50-200bp mutations. To explore the relative impact of various mutational classes on crop fitness, we then track these structural events over the last century of US rice improvement using 101 resequenced varieties. Within this material, a pattern of temporary hybridization between medium and long-grain varieties was followed by recent divergence. During this long-term selection, structural mutations that impact gene exons have been removed at a greater rate than intronic indels and single-nucleotide mutations. These results support the use of ab initio estimates of mutational burden, based on structural data, as an orthogonal predictor in genomic selection.

Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids.

Polyploidy is a widespread phenomenon throughout eukaryotes. Due to the coexistence of duplicated genomes, polyploids offer unique challenges for estimating gene expression levels, which is essential for understanding the massive and various forms of transcriptomic responses accompanying polyploidy. Although previous studies have explored the bioinformatics of polyploid transcriptomic profiling, the causes and consequences of inaccurate quantification of transcripts from duplicated gene copies have not been addressed. Using transcriptomic data from the cotton genus (Gossypium) as an example, we present an analytical workflow to evaluate a variety of bioinformatic method choices at different stages of RNA-seq analysis, from homoeolog expression quantification to downstream analysis used to infer key phenomena of polyploid expression evolution. In general, EAGLE-RC and GSNAP-PolyCat outperform other quantification pipelines tested, and their derived expression dataset best represents the expected homoeolog expression and co-expression divergence. The performance of co-expression network analysis was less affected by homoeolog quantification than by network construction methods, where weighted networks outperformed binary networks. By examining the extent and consequences of homoeolog read ambiguity, we illuminate the potential artifacts that may affect our understanding of duplicate gene expression, including an overestimation of homoeolog co-regulation and the incorrect inference of subgenome asymmetry in network topology. Taken together, our work points to a set of reasonable practices that we hope are broadly applicable to the evolutionary exploration of polyploids.

Innovations in double digest restriction-site associated DNA sequencing (ddRAD-Seq) method for more efficient SNP identification.

We present an improved ddRAD-Seq protocol for identifying single nucleotide polymorphisms (SNPs). It utilizes selected restriction enzyme digestion fragments, quick acting ligases that are neutral with the restriction enzyme buffer eliminating buffer exchange steps, and adapters designed to be compatible with Illumina index primers. Library amplification and barcoding are completed in one PCR step, and magnetic beads are used to purify the genomic fragments from the ligation and library generation steps. Our protocol increases the efficiency and decreases the time to complete a ddRAD-Seq experiment. To demonstrate its utility, we compared SNPs from our protocol with those from whole genome resequencing data from Gossypium herbaceum and Gossypium arboreum. Principal component analysis demonstrated that the variability of the combined data was explained by the genotype (PC1) and methodology applied (PC2). Phylogenetic analysis showed that the SNPs from our method clustered with SNPs from the resequencing data of the corresponding genotype. Sequence alignments illustrated that for homozygous loci, more than 90% of the SNPs from the resequencing data were discovered by our method. Our analyses suggest that our ddRAD-Seq method is reliable in identifying SNPs suitable for phylogenetic and association genetic studies while reducing cost and time over known methods.

Follicular-fluid proteomics during equine follicle development.

The complex composition of the follicular fluid (FF), the intimate proximity to the oocyte, and the continual changes in their composition have a major effect on folliculogenesis and oogenesis. To date, the profiling of FF proteomes during follicle selection, development, and ovulation has not been comprehensively investigated. Therefore, a shotgun proteomics approach and bioinformatics analyses were used to profile the proteomes of equine FF harvested in vivo from follicles at the following development stages: predeviation (18-20 mm), deviation (22-25 mm), postdeviation (26-29 mm), preovulatory (30-35 mm), and impending ovulation. A total of 294 proteins were detected in FF (FDR <1%), corresponding to 65 common proteins and 124, 142, 167, 132, and 142 proteins in the predeviation, deviation, postdeviation, preovulatory, and impending ovulation groups, respectively. The higher expression of properdin and several other proteins belonging to the complement system during the deviation time and ovulation suggested their contribution in the selection of the future dominant follicle and ovulation. Apolipoprotein A-1 and antithrombin-III appeared to be important throughout folliculogenesis. The "complement and coagulation cascades" was the major KEGG pathway across all stages of follicle development. The significant expression of several proteins belonging to the serine-type endopeptidase indicated their likely contribution to follicle and oocyte development. Our data provide an extensive description and functional analyses of the equine FF proteome during follicle selection, development, and ovulation. This information will help improve understanding of the ovarian function and ovulatory dysfunctions and might serve as a reference for future biomarker discovery for oocyte quality assessment.

Dysregulated FXR-FGF19 signaling and choline metabolism are associated with gut dysbiosis and hyperplasia in a novel pig model of pediatric NASH.

To investigate the role of bile acids (BAs) in the pathogenesis of diet-induced nonalcoholic steatohepatitis (NASH), we fed a "Western-style diet" [high fructose, high fat (HFF)] enriched with fructose, cholesterol, and saturated fat for 10 wk to juvenile Iberian pigs. We also supplemented probiotics with in vitro BA deconjugating activity to evaluate their potential therapeutic effect in NASH. Liver lipid and function, cytokines, and hormones were analyzed using commercially available kits. Metabolites, BAs, and fatty acids were measured by liquid chromatography-mass spectrometry. Histology and gene and protein expression analyses were performed using standard protocols. HFF-fed pigs developed NASH, cholestasis, and impaired enterohepatic Farnesoid-X receptor (FXR)-fibroblast growth factor 19 (FGF19) signaling in the absence of obesity and insulin resistance. Choline depletion in HFF livers was associated with decreased lipoprotein and cholesterol in serum and an increase of choline-containing phospholipids in colon contents and trimethylamine-N-oxide in the liver. Additionally, gut dysbiosis and hyperplasia increased with the severity of NASH, and were correlated with increased colonic levels of choline metabolites and secondary BAs. Supplementation of probiotics in the HFF diet enhanced NASH, inhibited hepatic autophagy, increased excretion of taurine and choline, and decreased gut microbial diversity. In conclusion, dysregulation of BA homeostasis was associated with injury and choline depletion in the liver, as well as increased biliary secretion, gut metabolism and excretion of choline-based phospholipids. Choline depletion limited lipoprotein synthesis, resulting in hepatic steatosis, whereas secondary BAs and choline-containing phospholipids in colon may have promoted dysbiosis, hyperplasia, and trimethylamine synthesis, causing further damage to the liver.NEW & NOTEWORTHY Impaired Farnesoid-X receptor (FXR)-fibroblast growth factor 19 (FGF19) signaling and cholestasis has been described in nonalcoholic fatty liver disease (NAFLD) patients. However, therapeutic interventions with FXR agonists have produced contradictory results. In a swine model of pediatric nonalcoholic steatohepatitis (NASH), we show that the uncoupling of intestinal FXR-FGF19 signaling and a decrease in FGF19 levels are associated with a choline-deficient phenotype of NASH and increased choline excretion in the gut, with the subsequent dysbiosis, colonic hyperplasia, and accumulation of trimethylamine-N-oxide in the liver.

Identification and characterization of microRNAs (miRNAs) and their transposable element origins in the saltwater crocodile, Crocodylus porosus.

MicroRNAs (miRNAs) are 18-24 nucleotide regulatory RNAs. They are involved in the regulation of genetic and biological pathways through post transcriptional gene silencing and/or translational repression. Data suggests a slow evolutionary rate for the saltwater crocodile (Crocodylus porosus) over the past several million years when compared to birds, the closest extant relatives of crocodilians. Understanding gene regulation in the saltwater crocodile in the context of relatively slow genomic change thus holds potential for the investigation of genomics, evolution, and adaptation. Utilizing eleven tissue types and sixteen small RNA libraries, we report 644 miRNAs in the saltwater crocodile with >78% of miRNAs being novel to crocodilians. We also identified potential targets for the miRNAs and analyzed the relationship of the miRNA repertoire to transposable elements (TEs). Results suggest an increased association of DNA transposons with miRNAs when compared to retrotransposons. This work reports the first comprehensive analysis of miRNAs in Crocodylus porosus and addresses the potential impacts of miRNAs in regulating the genome in the saltwater crocodile. In addition, the data suggests a supporting role of TEs as a source for miRNAs, adding to the increasing evidence that TEs play a significant role in the evolution of gene regulation.

Development of a qPCR Protocol to Detect the Cotton Bacterial Blight Pathogen, Xanthomonas citri pv. malvacearum, from Cotton Leaves and Seeds.

Bacterial blight, historically a seed-borne disease of cotton (Gossypium hirsutum) is caused by Xanthomonas citri pv. malvacearum, resulted in significant economic losses prior to development of resistant varieties and implementation of acid-delinting of planting seed. Periodic outbreaks have been associated with seed since the early twentieth century; of note, the disease has experienced a resurgence since 2011. Effective management of bacterial blight is dependent on accurate diagnosis and detection of the pathogen. Currently, detection of X. citri pv. malvacearum is performed by time-consuming microbiological methods. In this study, a novel and sensitive TaqMan-based qPCR protocol was developed to test for X. citri pv. malvacearum in cotton plant tissue. The primers developed are specific to five races of X. citri pv. malvacearum, but not to other Xanthomonas species or cotton-associated nonpathogenic bacteria. The efficiency of this assay was evaluated on artificially inoculated cotton leaves and seed, on naturally infected cotton leaves, and on bolls and seed originating from bacterial blight symptomatic bolls. The protocol's efficiency from artificially inoculated plant tissue was 102 copies g-1 and 37 copies from 1 g seed for leaves and seed, respectively. In addition, X. citri pv. malvacearum was detected from 94% of the seed samples originating from blight symptomatic bolls. The qPCR protocol provides a rapid and accurate method for diagnosis and detection of bacterial blight and offers a tool for monitoring X. citri pv. malvacearum and potentially reducing its spread in seed.

A genomic resource for the sedentary semi-endoparasitic reniform nematode, Rotylenchulus reniformis Linford & Oliveira.

The reniform nematode (Rotylenchulus reniformis) is a sedentary semi-endoparasitic species that is pathogenic on many row crops, fruits, and vegetables. Here, the authors present a draft genome assembly of R. reniformis using small- and large-insert libraries sequenced on the Illumina GAIIx and MiSeq platforms.The reniform nematode (Rotylenchulus reniformis) is a sedentary semi-endoparasitic species that is pathogenic on many row crops, fruits, and vegetables. Here, the authors present a draft genome assembly of R. reniformis using small- and large-insert libraries sequenced on the Illumina GAIIx and MiSeq platforms.

Transcriptomic changes in Echinochloa colona in response to treatment with the herbicide imazamox.

MAIN CONCLUSION: Presented here is the first Echinochloa colona leaf transcriptome. Analysis of gene expression before and after herbicide treatment reveals that E. colona mounts a stress response upon exposure to herbicide. Herbicides are the most frequently used means of controlling weeds. For many herbicides, the target site is known; however, it is considerably less clear how plant gene expression changes in response to herbicide exposure. In this study, changes in gene expression in response to herbicide exposure in imazamox-sensitive (S) and- resistant (R) junglerice (Echinochloa colona L.) biotypes was examined. As no reference genome is available for this weed, a reference leaf transcriptome was generated. Messenger RNA was isolated from imazamox-treated- and untreated R and S plants and the resulting cDNA libraries were sequenced on an Illumina HiSeq2000. The transcriptome was assembled, annotated, and differential gene expression analysis was performed to identify transcripts that were upregulated or downregulated in response to herbicide exposure for both biotypes. Differentially expressed transcripts included transcription factors, protein-modifying enzymes, and enzymes involved in metabolism and signaling. A literature search revealed that members of the families represented in this analysis were known to be involved in abiotic stress response in other plants, suggesting that imazamox exposure induced a stress response. A time course study examining a subset of transcripts showed that expression peaked within 4-12 h and then returned to untreated levels within 48 h of exposure. Testing of plants from two additional biotypes showed a similar change in gene expression 4 h after herbicide exposure compared to the resistant and sensitive biotypes. This study shows that within 48 h junglerice mounts a stress response to imazamox exposure.

Quack: A quality assurance tool for high throughput sequence data.

The quality of data generated by high-throughput DNA sequencing tools must be rapidly assessed in order to determine how useful the data may be in making biological discoveries; higher quality data leads to more confident results and conclusions. Due to the ever-increasing size of data sets and the importance of rapid quality assessment, tools that analyze sequencing data should quickly produce easily interpretable graphics. Quack addresses these issues by generating information-dense visualizations from FASTQ files at a speed far surpassing other publicly available quality assurance tools in a manner independent of sequencing technology.

Five life stage-specific transcriptome assemblies for the reniform nematode, Rotylenchulus reniformis Linford & Oliveira.

The reniform nematode (Rotylenchulus reniformis Linford and Oliveira) is a semi-endoparasitic nematode that is a pathogen of numerous major crops such as cotton and soybean. Here, the authors present transcriptome assemblies of the egg, second-stage juvenile (J2), J3, vermiform adult, and sedentary female life stages of this important plant pathogen.The reniform nematode (Rotylenchulus reniformis Linford and Oliveira) is a semi-endoparasitic nematode that is a pathogen of numerous major crops such as cotton and soybean. Here, the authors present transcriptome assemblies of the egg, second-stage juvenile (J2), J3, vermiform adult, and sedentary female life stages of this important plant pathogen.

Bacterial Diversity Patterns Differ in Soils Developing in Sub-tropical and Cool-Temperate Ecosystems.

Microbial diversity patterns have been surveyed in many different soils and ecosystems, but we are unaware of studies comparing similar soils developing from similar parent materials in contrasting climates. In 2008, developmental chronosequences with ages ranging from 105 to 500,000 years across Georgia (GA) and Michigan (MI) were studied to investigate how bacterial community composition and diversity change as a result of local environmental gradients that develop during pedogenesis. Geographic factors were studied between and within locations spanning two scales: (1) regionally between 0.1 and 50 and (2) ∼1700 km apart. The diversity was surveyed using high-throughput pyrosequencing, and variance partitioning was used to describe the effects of spatial, environmental, and spatio-environmental factors on bacterial community composition. At the local scale, variation in bacterial communities was most closely related to environmental factors (rM = 0.59, p = 0.0001). There were differences in bacterial communities between the two locations, indicating spatial biogeography. Estimates of bacterial diversity were much greater in MI (numbers of OTU, ACE, and Chao1) and remained 2-3× greater in MI than GA after removing the effect of soil properties. The large differences in diversity between geographically separated bacterial communities in different climates need further investigation. It is not known if the rare members of the community, which contributed to greater bacterial diversity in GA relative to MI, play an important role in ecosystem function but has been hypothesized to play a role in ecosystem resiliency, resistance, and stability. Further research on the link between bacterial diversity and spatial variability related to climate needs further investigation.

Fine mapping and identification of candidate genes for a QTL affecting Meloidogyne incognita reproduction in Upland cotton.

BACKGROUND: The southern root-knot nematode (Meloidogyne incognita; RKN) is one of the most important economic pests of Upland cotton (Gossypium hirsutum L.). Host plant resistance, the ability of a plant to suppress nematode reproduction, is the most economical, practical, and environmentally sound method to provide protection against this subterranean pest. The resistant line Auburn 623RNR and a number of elite breeding lines derived from it remain the most important source of root-knot nematode (RKN) resistance. Prior genetic analysis has identified two epistatically interacting RKN resistance QTLs, qMi-C11 and qMi-C14, affecting gall formation and RKN reproduction, respectively. RESULTS: We developed a genetic population segregating only for the qMi-C14 locus and evaluated the genetic effects of this QTL on RKN resistance in the absence of the qMi-C11 locus. The qMi-C14 locus had a LOD score of 12 and accounted for 24.5 % of total phenotypic variation for egg production. In addition to not being significantly associated with gall formation, this locus had a lower main effect on RKN reproduction than found in our previous study, which lends further support to evidence of epistasis with qMi-C11 in imparting RKN resistance in the Auburn 623RNR source. The locus qMi-C14 was fine-mapped with the addition of 16 newly developed markers. By using the reference genome sequence of G. raimondii, we identified 20 candidate genes encoding disease resistance protein homologs in the newly defined 2.3 Mb region flanked by two SSR markers. Resequencing of an RKN resistant and susceptible G. hirsutum germplasm revealed non-synonymous mutations in only four of the coding regions of candidate genes, and these four genes are consequently of high interest. CONCLUSIONS: Our mapping results validated the effects of the qMi-C14 resistance locus, delimiting the QTL to a smaller region, and identified tightly linked SSR markers to improve the efficiency of marker-assisted selection. The candidate genes identified warrant functional studies that will help in identifying and characterizing the actual qMi-C14 defense gene(s) against root-knot nematodes.

Epigenomics and bolting tolerance in sugar beet genotypes.

In sugar beet (Beta vulgaris altissima), bolting tolerance is an essential agronomic trait reflecting the bolting response of genotypes after vernalization. Genes involved in induction of sugar beet bolting have now been identified, and evidence suggests that epigenetic factors are involved in their control. Indeed, the time course and amplitude of DNA methylation variations in the shoot apical meristem have been shown to be critical in inducing sugar beet bolting, and a few functional targets of DNA methylation during vernalization have been identified. However, molecular mechanisms controlling bolting tolerance levels among genotypes are still poorly understood. Here, gene expression and DNA methylation profiles were compared in shoot apical meristems of three bolting-resistant and three bolting-sensitive genotypes after vernalization. Using Cot fractionation followed by 454 sequencing of the isolated low-copy DNA, 6231 contigs were obtained that were used along with public sugar beet DNA sequences to design custom Agilent microarrays for expression (56k) and methylation (244k) analyses. A total of 169 differentially expressed genes and 111 differentially methylated regions were identified between resistant and sensitive vernalized genotypes. Fourteen sequences were both differentially expressed and differentially methylated, with a negative correlation between their methylation and expression levels. Genes involved in cold perception, phytohormone signalling, and flowering induction were over-represented and collectively represent an integrative gene network from environmental perception to bolting induction. Altogether, the data suggest that the genotype-dependent control of DNA methylation and expression of an integrative gene network participate in bolting tolerance in sugar beet, opening up perspectives for crop improvement.

The Sorghum bicolor genome and the diversification of grasses.

Sorghum, an African grass related to sugar cane and maize, is grown for food, feed, fibre and fuel. We present an initial analysis of the approximately 730-megabase Sorghum bicolor (L.) Moench genome, placing approximately 98% of genes in their chromosomal context using whole-genome shotgun sequence validated by genetic, physical and syntenic information. Genetic recombination is largely confined to about one-third of the sorghum genome with gene order and density similar to those of rice. Retrotransposon accumulation in recombinationally recalcitrant heterochromatin explains the approximately 75% larger genome size of sorghum compared with rice. Although gene and repetitive DNA distributions have been preserved since palaeopolyploidization approximately 70 million years ago, most duplicated gene sets lost one member before the sorghum-rice divergence. Concerted evolution makes one duplicated chromosomal segment appear to be only a few million years old. About 24% of genes are grass-specific and 7% are sorghum-specific. Recent gene and microRNA duplications may contribute to sorghum's drought tolerance.

Transcriptomic dissection of the rice-Burkholderia glumae interaction.

BACKGROUND: Bacterial panicle blight caused by the bacterium Burkholderia glumae is an emerging disease of rice in the United States. Not much is known about this disease, the disease cycle or any source of disease resistance. To understand the interaction between rice and Burkholderia glumae, we used transcriptomics via next-generation sequencing (RNA-Seq) and bioinformatics to identify differentially expressed transcripts between resistant and susceptible interactions and formulate a model for rice resistance to the disease. RESULTS: Using inoculated young seedlings as sample tissues, we identified unique transcripts involved with resistance to bacterial panicle blight, including a PIF-like ORF1 and verified differential expression of some selected genes using qRT-PCR. These transcripts, which include resistance genes of the NBS-LRR type, kinases, transcription factors, transporters and expressed proteins with functions that are not known, have not been reported in other pathosystems including rice blast or bacterial blight. Further, functional annotation analysis reveals enrichment of defense response and programmed cell death (biological processes); ATP and protein binding (molecular functions); and mitochondrion-related (cell component) transcripts in the resistant interaction. CONCLUSION: Taken together, we formulated a model for rice resistance to bacterial panicle blight that involves an activation of previously unknown resistance genes and their activation partners upon challenge with B. glumae. Other interesting findings are that 1) though these resistance transcripts were up-regulated upon inoculation in the resistant interaction, some of them were already expressed in the water-inoculated control from the resistant genotype, but not in the water- and bacterium-inoculated samples from the susceptible genotype; 2) rice may have co-opted an ORF that was previously a part of a transposable element to aid in the resistance mechanism; and 3) resistance may have existed immediately prior to rice domestication.

Molecular Mechanisms Associated with the Development of the Metritis Complex in Dairy Cattle.

The metritis complex (MC), a group of post-partum uterine diseases, is associated with increased treatment costs and reduced milk yield and fertility. The goal of this study was to identify genetic variants, genes, or genomic regions that modulate MC disease. A genome-wide association study was performed using a single-locus mixed linear model of 1967 genotypes (624,460 SNPs) and metritis complex records. Then, in-silico functional analyses were performed to detect biological mechanisms and pathways associated with the development of MC. The ATP8A2, COX16, AMN, and TRAF3 genes, located on chromosomes 12, 10, and 21, were associated with MC at p ≤ 0.0001. These genes are involved in the regulation of cholesterol metabolism in the stromal tissue of the uterus, which can be directly associated with the mode of transmission for pathogens causing the metritis complex. The modulation of cholesterol abundance alters the efficiency of virulence factors and may affect the susceptibility of the host to infection. The SIPA1L1, DEPDC5, and RNF122 genes were also significantly associated with MC at p ≤ 0.0001 and are involved in the PI3k-Akt pathway, responsible for activating the autophagic processes. Thus, the dysregulation of these genes allows for unhindered bacterial invasion, replication, and survival within the endometrium.

Complete genome sequences of two avian pathogenic Escherichia coli strains isolated from broilers exhibiting colibacillosis in Mississippi.

We report the whole-genome sequences of Escherichia coli strains APEC-O2-MS1266 and APEC-O2-MS1657 isolated from the liver and heart of infected broilers in Mississippi State, US. The genomic information of these two causative strains may provide a valuable reference for comparative studies of avian pathogenic E. coli.

Origin and diversity of the wild cottons (Gossypium hirsutum) of Mound Key, Florida.

Elucidating genetic diversity within wild forms of modern crops is essential for understanding domestication and the possibilities of wild germplasm utilization. Gossypium hirsutum is a predominant source of natural plant fibers and the most widely cultivated cotton species. Wild forms of G. hirsutum are challenging to distinguish from feral derivatives, and truly wild populations are uncommon. Here we characterize a population from Mound Key Archaeological State Park, Florida using genome-wide SNPs extracted from 25 individuals over three sites. Our results reveal that this population is genetically dissimilar from other known wild, landrace, and domesticated cottons, and likely represents a pocket of previously unrecognized wild genetic diversity. The unexpected level of divergence between the Mound Key population and other wild cotton populations suggests that the species may harbor other remnant and genetically distinct populations that are geographically scattered in suitable habitats throughout the Caribbean. Our work thus has broader conservation genetic implications and suggests that further exploration of natural diversity in this species is warranted.