A narrow germplasm base and a complex allotetraploid genome have historically made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). We conducted a genome reduction experiment to identify SNPs from two accessions of G. hirsutum and two accessions of G. barbadense. Approximately 2 million sequence reads were assembled into contigs with an N50 of 508 bp and analyzed for SNPs. A total of 11,834 and 1,679 SNPs between the accessions G. hirsutum and G. barbadense, respectively, were identified with highly conservative parameters (a minimum read depth of 8x at each SNP and a 100% identity of all reads within an accession at the SNP). Additionally, 4,327 SNPs were identified between accessions of G. hirsutum in and assembly of Expressed Sequence Tags (ESTs). 320 and 252 KASPAR assays were designed for SNP mapping in non-genic and genic regions respectively. 187 markers in total (136 non-genic, 51 genic) were mapped using KBioscience KASPar genotyping assays in a segregating F2 population using the Fluidigm EP1 system. EST The target genome of EST markers was successfully predicted bioinformaticly diploid reference sequences. Examination of nucleotide substitutions and SNP frequencies further confirms validity of new markers. A genetic map was constructed using a large G. hirsutum segregating F2 population. Genetic maps generated by these newly identified markers will be used to locate quantitative, economically important regions within the cotton genome.

Polyploidization instantly doubles all genome content by combining two genomes that have markedly different methylation and gene expression levels. This process may be accompanied by genetic and epigenetic changes in each genome. Sequencing of the transcriptome (RNA-seq) and the methylome (bisulfite treated libraries whole genome libraries) were used to measure gene expression and methylation levels of genic regions of allopolyploid cotton petals and petals of their diploid relatives. Many differentially expressed genes detected by RNA-seq were consistent with expression levels previously detected by microarrays. RNA-seq results also reconfirmed the presence of general polyploid gene expression trends like expression level dominance and homoeologous expression biases in Gossypium polyploid species. Expression biases between A- and D-genome homoeologs and expression level dominance was characterized for thousands of genes in tetraploids and a diploid F1-hybrid. Unlike the results of microarray study previously done we found a slightly greater number of genes showing A-genome bias vs genes showing D-genome bias. More commonly the overall expression level from homoeologs of polyploid is heterotic i.e the expression level is greater than the average of the expression levels from the two parent genomes. In addition, genome methylation (CG, CHG, and CHH contexts) of each genome was assessed in the diploid and tetraploid samples. The A- and D-genomes had distinct levels of DNA methylation for each context. DNA methylation may be independently regulating homoeologous expression levels of a small number of genes.