A Gossypium BAC clone contains key repeat components distinguishing sub-genome of allotetraploidy cottons

Cloning and preliminary verification of telomere-associated sequences in upland cotton

Telomeres are structures enriched in repetitive sequences at the end of chromosomes. In this study, using the telomere primer AA(CCCTAAA)₃CCC for the single primer PCR, two DNA sequences were obtained from Gossypium hirsutum (Linnaeus, 1753) accession (acc.) TM-1. Sequence analysis showed that the two obtained sequences were all rich in A/T base, which was consistent with the characteristic of the telomere-associated sequence (TAS). They were designated as GhTAS1 and GhTAS2 respectively. GhTAS1 is 489 bp long, with 57.6% of A/T, and GhTAS2 is 539 bp long, with 63.9% of A/T. Fluorescence in situ hybridization results showed that both of the cloned TASs were located at the ends of the partial chromosomes of G. hirsutum, with the strong signals, which further confirmed that GhTAS1 and GhTAS2 were telomere-associated sequences including highly tandemly repetitive sequences. Results of blast against the assembled genome of G. hirsutum showed that GhTAS sequences may be missed on some assembled chromosomes. The results provide important evidence for the evaluation of the integrity of assembled chromosome end sequences, and will also contribute to the further perfection of the draft genomes of cotton.

Identification of a genome-specific repetitive element in the Gossypium D genome

The activity of genome-specific repetitive sequences is the main cause of genome variation between Gossypium A and D genomes. Through comparative analysis of the two genomes, we retrieved a repetitive element termed ICRd motif, which appears frequently in the diploid Gossypium raimondii (D₅) genome but rarely in the diploid Gossypium arboreum (A₂) genome. We further explored the existence of the ICRd motif in chromosomes of G. raimondii, G. arboreum, and two tetraploid (AADD) cotton species, Gossypium hirsutum and Gossypium barbadense, by fluorescence in situ hybridization (FISH), and observed that the ICRd motif exists in the D₅ and D-subgenomes but not in the A₂ and A-subgenomes. The ICRd motif comprises two components, a variable tandem repeat (TR) region and a conservative sequence (CS). The two constituents each have hundreds of repeats that evenly distribute across 13 chromosomes of the D₅genome. The ICRd motif (and its repeats) was revealed as the common conservative region harbored by ancient Long Terminal Repeat Retrotransposons. Identification and investigation of the ICRd motif promotes the study of A and D genome differences, facilitates research on Gossypium genome evolution, and provides assistance to subgenome identification and genome assembling.

Identification of a genome-specific repetitive element in the Gossypium D genome

The activity of genome-specific repetitive sequences is the main cause of genome variation between Gossypium A and D genomes. Through comparative analysis of the two genomes, we retrieved a repetitive element termed ICRd motif, which appears frequently in the diploid Gossypium raimondii (D₅) genome but rarely in the diploid Gossypium arboreum (A₂) genome. We further explored the existence of the ICRd motif in chromosomes of G. raimondii, G. arboreum, and two tetraploid (AADD) cotton species, Gossypium hirsutum and Gossypium barbadense, by fluorescence in situ hybridization (FISH), and observed that the ICRd motif exists in the D₅ and D-subgenomes but not in the A₂ and A-subgenomes. The ICRd motif comprises two components, a variable tandem repeat (TR) region and a conservative sequence (CS). The two constituents each have hundreds of repeats that evenly distribute across 13 chromosomes of the D₅genome. The ICRd motif (and its repeats) was revealed as the common conservative region harbored by ancient Long Terminal Repeat Retrotransposons. Identification and investigation of the ICRd motif promotes the study of A and D genome differences, facilitates research on Gossypium genome evolution, and provides assistance to subgenome identification and genome assembling.

Chromosome Painting Based on Bulked Oligonucleotides in Cotton

Chromosome painting is one of the key technologies in cytogenetic research, which can accurately identify chromosomes or chromosome regions. Oligonucleotide (oligo) probes designed based on genome sequences have both flexibility and specificity, which would be ideal probes for fluorescence in situ hybridization (FISH) analysis of genome structure. In this study, the bulked oligos of the two arms of chromosome seven of cotton were developed based on the genome sequence of Gossypium raimondii (DD, 2n = 2× = 26), and each arm contains 12,544 oligos. Chromosome seven was easily identified in both D genome and AD genome cotton species using the bulked chromosome-specific painting probes. Together with 45S ribosomal DNA (rDNA) probe, the chromosome-specific painting probe was also successfully used to correct the chromosomal localization of 45S rDNA in G. raimondii. The study reveals that bulked oligos specific to a chromosome is a useful tool for chromosome painting in cotton.

Discovery and annotation of a novel transposable element family in Gossypium

BACKGROUND

Fluorescence in situ hybridization (FISH) is an efficient cytogenetic technology to study chromosome structure. Transposable element (TE) is an important component in eukaryotic genomes and can provide insights in the structure and evolution of eukaryotic genomes.

RESULTS

A FISH probe derived from bacterial artificial chromosome (BAC) clone 299N22 generated striking signals on all 26 chromosomes of the cotton diploid A genome (AA, 2x=26) but very few on the diploid D genome (DD, 2x=26). All 26 chromosomes of the A sub genome (At) of tetraploid cotton (AADD, 2n=4x=52) also gave positive signals with this FISH probe, whereas very few signals were observed on the D sub genome (Dt). Sequencing and annotation of BAC clone 299N22, revealed a novel Ty3/gypsy transposon family, which was named as 'CICR'. This family is a significant contributor to size expansion in the A (sub) genome but not in the D (sub) genome. Further FISH analysis with the LTR of CICR as a probe revealed that CICR is lineage-specific, since massive repeats were found in A and B genomic groups, but not in C-G genomic groups within the Gossypium genus. Molecular evolutionary analysis of CICR suggested that tetraploid cottons evolved after silence of the transposon family 1-1.5 million years ago (Mya). Furthermore, A genomes are more homologous with B genomes, and the C, E, F, and G genomes likely diverged from a common ancestor prior to 3.5-4 Mya, the time when CICR appeared. The genomic variation caused by the insertion of CICR in the A (sub) genome may have played an important role in the speciation of organisms with A genomes.

CONCLUSIONS

The CICR family is highly repetitive in A and B genomes of Gossypium, but not amplified in the C-G genomes. The differential amount of CICR family in At and Dt will aid in partitioning sub genome sequences for chromosome assemblies during tetraploid genome sequencing and will act as a method for assessing the accuracy of tetraploid genomes by looking at the proportion of CICR elements in resulting pseudochromosome sequences. The timeline of the expansion of CICR family provides a new reference for cotton evolutionary analysis, while the impact on gene function caused by the insertion of CICR elements will be a target for further analysis of investigating phenotypic differences between A genome and D genome species.

Construction and characterization of a bacterial artificial chromosome library for Gossypium mustelinum

A bacterial artificial chromosome (BAC) library for G. mustelinum Miers ex G. Watt (AD₄) was constructed. Intact nuclei from G. mustelinum (AD₄) were used to isolate high molecular weight DNA, which was partially cleaved with Hind III and cloned into pSMART BAC (Hind III) vectors. The BAC library consisted of 208,182 clones arrayed in 542 384-microtiter plates, with an average insert size of 121.72 kb ranging from 100 to 150 kb. About 2% of the clones did not contain inserts. Based on an estimated genome size of 2372 Mb for G. mustelinum, the BAC library was estimated to have a total coverage of 10.50 × genome equivalents. The high capacity library of G. mustelinum will serve as a giant gene resource for map-based cloning of quantitative trait loci or genes associated with important agronomic traits or resistance to Verticillium wilt, physical mapping and comparative genome analysis.

Identification and Characterization of miRNA Transcriptome in Asiatic Cotton (Gossypium arboreum) Using High Throughput Sequencing

MicroRNAs (miRNAs) are small 20-24nt molecules that have been well studied over the past decade due to their important regulatory roles in different cellular processes. The mature sequences are more conserved across vast phylogenetic scales than their precursors and some are conserved within entire kingdoms, hence, their loci and function can be predicted by homology searches. Different studies have been performed to elucidate miRNAs using de novo prediction methods but due to complex regulatory mechanisms or false positive in silico predictions, not all of them express in reality and sometimes computationally predicted mature transcripts differ from the actual expressed ones. With the availability of a complete genome sequence of Gossypium arboreum, it is important to annotate the genome for both coding and non-coding regions using high confidence transcript evidence, for this cotton species that is highly resistant to various biotic and abiotic stresses. Here we have analyzed the small RNA transcriptome of G. arboreum leaves and provided genome annotation of miRNAs with evidence from miRNA/miRNA∗ transcripts. A total of 446 miRNAs clustered into 224 miRNA families were found, among which 48 families are conserved in other plants and 176 are novel. Four short RNA libraries were used to shortlist best predictions based on high reads per million. The size, origin, copy numbers and transcript depth of all miRNAs along with their isoforms and targets has been reported. The highest gene copy number was observed for gar-miR7504 followed by gar-miR166, gar-miR8771, gar-miR156, and gar-miR7484. Altogether, 1274 target genes were found in G. arboreum that are enriched for 216 KEGG pathways. The resultant genomic annotations are provided in UCSC, BED format.