The cotton centromere contains a Ty3-gypsy-like LTR retroelement

The centromere is a repeat-rich structure essential for chromosome segregation; with the long-term aim of understanding centromere structure and function, we set out to identify cotton centromere sequences. To isolate centromere-associated sequences from cotton, (Gossypium hirsutum) we surveyed tandem and dispersed repetitive DNA in the genus. Centromere-associated elements in other plants include tandem repeats and, in some cases, centromere-specific retroelements. Examination of cotton genomic survey sequences for tandem repeats yielded sequences that did not localize to the centromere. However, among the repetitive sequences we also identified a gypsy-like LTR retrotransposon (Centromere Retroelement Gossypium, CRG) that localizes to the centromere region of all chromosomes in domestic upland cotton, Gossypium hirsutum, the major commercially grown cotton. The location of the functional centromere was confirmed by immunostaining with antiserum to the centromere-specific histone CENH3, which co-localizes with CRG hybridization on metaphase mitotic chromosomes. G. hirsutum is an allotetraploid composed of A and D genomes and CRG is also present in the centromere regions of other AD cotton species. Furthermore, FISH and genomic dot blot hybridization revealed that CRG is found in D-genome diploid cotton species, but not in A-genome diploid species, indicating that this retroelement may have invaded the A-genome centromeres during allopolyploid formation and amplified during evolutionary history. CRG is also found in other diploid Gossypium species, including B and E2 genome species, but not in the C, E1, F, and G genome species tested. Isolation of this centromere-specific retrotransposon from Gossypium provides a probe for further understanding of centromere structure, and a tool for future engineering of centromere mini-chromosomes in this important crop species.     

Analysis of flavanone 3-hydroxylase in Arabidopsis seedlings. Coordinate regulation with chalcone synthase and chalcone isomerase.

A genomic clone encoding flavanone 3-hydroxylase (F3H) was isolated from Arabidopsis thaliana. The deduced amino acid sequence is 72 to 94% identical to all previously reported F3H proteins. Low-stringency DNA blot analysis indicated that F3H is encoded by a single gene in Arabidopsis. The F3H locus was mapped to the bottom of chromosome 3 and therefore does not correspond to any of the 13 flavonoid-deficient transparent testa mutants for which a map position is known. Analysis of gene expression in etiolated seedlings exposed to white light and in two putative regulatory mutants, ttg and tt8, demonstrated that the Arabidopsis F3H gene is coordinately expressed with chalcone synthase and chalcone isomerases is seedlings, whereas dihydroflavonol reductase expression is controlled by distinct regulatory mechanisms. The F3H gene may represent a pivotal point in the regulation of flavonoid biosynthesis because its expression is coordinated with different subsets of genes in different plant species.     

Resistance gene analogue markers are mapped to homeologous chromosomes in cultivated tetraploid cotton

Degenerate primers designed from conserved motifs of known plant resistance gene products were used to amplify genomic DNA sequences from the root-knot nematode (Meloidogyne incognita) resistance genetic source, Upland cotton (Gossypium hirsutum) cultivar Auburn 634 RNR. A total of 165 clones were isolated, and sequence analysis revealed 57 of the clones to be novel nucleotide sequences, many containing the resistance (R)-protein nucleotide-binding site motif. A cluster analysis was performed with resistance gene analogue (RGA) nucleotide sequences isolated in this study, in addition to 99 cotton RGA nucleotide sequences already deposited in GenBank, to generate a phylogenetic tree of cotton R genes. The cotton RGA nucleotide sequences were arranged into 11 groups and 56 sub-groups, based on genetic distances. Multiple sequence alignments were performed on the RGA sequences of each sub-group, and either the consensus sequences or individual RGA sequences were used to design 61 RGA-sequence-tagged site primers. A recombinant inbred line (RIL) population of cultivated tetraploid cotton was genotyped using RGA-specific primers that amplified polymorphic fragments between the two RIL parents. Nine RGA markers were mapped to homeologous chromosomes 12 and 26, based on linkage to existing markers that are located on these chromosomes.     

Cot-based sampling of genomes for polymorphic low-copy DNA

DNA polymorphisms are powerful tools for many evolutionary and genomic studies in plants including molecular breeding. Single nucleotide polymorphisms (SNPs) are the most elemental DNA marker for genomic studies, but even with advances in DNA sequencing technology, SNP discovery remains costly and computationally demanding, especially in large genomes that are rich in repetitive DNA such as those of many plants. Here we report a method using DNA renaturation kinetics (Cot techniques), sequencing, and BLAST-based screening to identify low-copy, non-coding DNA sequences that were subsequently found to be relatively rich in polymorphisms. A total of of 63 such fragments isolated from a diploid D genome cotton species (Gossypium raimondii) revealed a higher frequency of polymorphisms than that observed for cotton expressed sequence tags or hypomethylated (PstI-susceptible) genomic DNA. While microsatellite-derived loci show still higher polymorphism rates, they often fall in repetitive elements and their sequence analysis is often complicated by alignment difficulties. The potential applications of Cot-filtered noncoding (CFNC) DNA in development of DNA markers are discussed.     

Genome-Wide Analysis of the Sus Gene Family in Cotton

Sucrose synthase (Sus) is a key enzyme in plant sucrose metabolism. In cotton, Sus (EC 2.4.1.13) is the main enzyme that degrades sucrose imported into cotton fibers from the phloem of the seed coat. This study demonstrated that the genomes of Gossypium arboreum L., G. raimondii Ulbr., and G. hirsutum L., contained 8, 8, and 15 Sus genes, respectively. Their structural organizations, phylogenetic relationships, and expression profiles were characterized. Comparisons of genomic and coding sequences identified multiple introns, the number and positions of which were highly conserved between diploid and allotetraploid cotton species. Most of the phylogenetic clades contained sequences from all three species, suggesting that the Sus genes of tetraploid G. hirsutum derived from those of its diploid ancestors. One Sus group (Sus I) underwent expansion during cotton evolution. Expression analyses indicated that most Sus genes were differentially expressed in various tissues and had development-dependent expression profiles in cotton fiber cells. Members of the same orthologous group had very similar expression patterns in all three species. These results provide new insights into the evolution of the cotton Sus gene family, and insight into its members' physiological functions during fiber growth and development.     

EST derived PCR-based markers for functional gene homologues in cotton.

We investigated the utility of the Gossypium arboreum EST sequences in the GenBank database for developing PCR-based markers targeting known-function genes in cultivated tetraploid cottons, G. hirsutum and G. barbadense. Four hundred sixty-five randomly selected ESTs from this library were subjected to BLASTn search against all GenBank databases, of which putative function was assigned to 93 ESTs based on high nucleotide homology to previously studied genes. PCR primers were synthesized for 89 of the known-function ESTs. A total of 57 primer pairs amplified G. arboreum genomic DNA, but only 39 amplified in G. hirsutum and G. barbadense, suggesting that sequence divergence may be a factor causing non-amplification for some sites. DNA sequence analysis showed that most primer pairs were targeting the expected homologous loci. While the amplified products that were of larger size than the corresponding EST sequences contain introns, the primer pairs with a smaller amplicon than predicted from the flanking EST sequences did not amplify the expected orthologous gene sequences. Among the 39 primer pairs that amplified tetraploid cotton DNA, 3 detected amplicon size polymorphisms and 10 detected polymorphisms after digestion with one of six restriction enzymes. Ten of the polymorphic loci were subsequently mapped to an anchor RFLP map. Digestion of PCR-amplified sequences offers one means by which cotton genes can be mapped to their chromosomal locations more quickly and economically than by RFLP analysis.     

Molecular evolution and phylogenetic analysis of genes related to cotton fibers development from wild and domesticated cotton species in Gossypium

The domestication of both diploid and tetraploid cotton species was carried out for fiber utilization. To understand the origin and domestication of fibers, 18 genes related to fiber development were individually cloned and sequenced from 22 different cotton species. Their structures, phylogenetic relationship and molecular evolution were further studied. In the orthologous and homeologous loci of the 18 genes, the sequence and structure of 72.22% were conserved and 27.78% were diverse. Tree topologies constructed based on the combined sequences showed that all 13 D-genome species were congruent with Fryxell's subsection taxonomy, the A- and D-subgenomes independently evolved in the allopolyploid after polyploid formation, and Gossypium raimondii had the closest relationship with all allotetraploids of D-subgenomes. The molecular evolutionary rates revealed approximately equivalent rates among different D-genome species, and purifying selection acted on all genes in the wild D-genome species. Among orthologs and homeologs, the D-subgenomes had higher evolutionary rates than the A-subgenomes in tetraploid cotton species, and the cultivars had higher evolutionary rates than either the semi-domesticated or wild species. Our study revealed that human domestication altered the molecular evolutionary pattern of genes related to fiber development, and Gossypium hirsutum endured greater selective pressures than Gossypium barbadense during the domestication process.     

Sequence Divergence of Microsatellites and Phylogeny Analysis in Tetraploid Cotton Species and Their Putative Diploid Ancestors

To determine the level of microsatellite sequence differences and to use the information to construct a phylogenetic relationship for cultivated tetraploid cotton （Gossypium spp.） species and their putative diploid ancestors, 10 genome-derived microsatellite primer pairs were used to amplify eight species, including two tetraploid and six diploid species, in Gossypium. A total of 92 unique amplicons were resolved using polyacrylamide gel electrophoresis. Each amplicon was cloned, sequenced, and analyzed using standard phylogenetic software. Allelic diversities were caused mostly by changes in the number of simple sequence repeat （SSR） motif repeats and only a small proportion resulted from interruption of the SSR motif within the locus for the same genome. The frequency of base substitutions was 0.5%-1.0% in different genomes, with only few indels found. Based on the combined 10 SSR flanking sequence data, the homology of A-genome diploid species averaged 98.9%, even though most of the amplicons were of the same size, and the sequence homology between G. gossypioides （Ulbr.） Standl. and three other D-genome species （G. raimondii Ulbr., G. davidsonii Kell., and G. thurberi Tod.） was 98.5%, 98.6%, and 98.5%, respectively. Phylogenetic trees of the two allotetraploid species and their putative diploid progenitors showed that homoelogous sequences from the A- and D-subgenome were still present in the polyploid subgenomes and they evolved independently. Meanwhile, homoelogous sequence interaction that duplicated loci in the polyploid subgenomes became phylogenetic sisters was also found in the evolutionary history of tetraploid cotton species. The results of the present study suggest that evaluation of SSR variation at the sequence level can be effective in exploring the evolutionary relationships among Gossypuim species.     

Cloning and Expression Studies of Novel Small RNAs in Tetraploid Cotton

Small RNAs are a group of non-coding RNAs that downregulate gene expression in a sequence-specific manner to control plant growth and development. The objective of the present study was to clone and characterize several small RNAs in cotton. To identify small RNAs that are involved in the development of cotton bolls and fibers, we generated cDNA libraries from cotton bolls at 13?days post-anthesis from two cotton cultivars, Pima Phy 76 (Gossypium bardadense) and Acala 1517?C99 (Gossypium hirsutum). Screening of these libraries identified eight small RNAs, seven of which have not been reported in other plant species and appear to be absent in the known sequences of other plant species. Their predicted target genes are known to be involved in cotton fiber development. The cloned small RNAs displayed lower and differential expression in the examined boll developmental stages using RT-PCR and quantitative RT-PCR. The genetic polymorphism of the small RNAs at the DNA level was evaluated by miRNA-amplified fragment length polymorphism (AFLP) analysis using primers designed from the small miRNA genes in combination with AFLP primers. Homologous small RNA gene sequences were further isolated using this homology-based genotyping approach, and potential hairpin structures were identified. The results represent a novel method to isolate small including miRNA genes at the RNA and DNA levels in many plant species where genome sequences are not available or expressed sequence tags are limited.     

Development of chromosome-specific markers with high polymorphism for allotetraploid cotton based on genome-wide characterization of simple sequence repeats in diploid cottons (Gossypium arboreum L. and Gossypium raimondii Ulbrich)

Background

Tetraploid cotton contains two sets of homologous chromosomes, the At- and Dt-subgenomes. Consequently, many markers in cotton were mapped to multiple positions during linkage genetic map construction, posing a challenge to anchoring linkage groups and mapping economically-important genes to particular chromosomes. Chromosome-specific markers could solve this problem. Recently, the genomes of two diploid species were sequenced whose progenitors were putative contributors of the At- and Dt-subgenomes to tetraploid cotton. These sequences provide a powerful tool for developing chromosome-specific markers given the high level of synteny among tetraploid and diploid cotton genomes. In this study, simple sequence repeats (SSRs) on each chromosome in the two diploid genomes were characterized. Chromosome-specific SSRs were developed by comparative analysis and proved to distinguish chromosomes.

Results

A total of 200,744 and 142,409 SSRs were detected on the 13 chromosomes of Gossypium arboreum L. and Gossypium raimondii Ulbrich, respectively. Chromosome-specific SSRs were obtained by comparing SSR flanking sequences from each chromosome with those from the other 25 chromosomes. The average was 7,996 per chromosome. To confirm their chromosome specificity, these SSRs were used to distinguish two homologous chromosomes in tetraploid cotton through linkage group construction. The chromosome-specific SSRs and previously-reported chromosome markers were grouped together, and no marker mapped to another homologous chromosome, proving that the chromosome-specific SSRs were unique and could distinguish homologous chromosomes in tetraploid cotton. Because longer dinucleotide AT-rich repeats were the most polymorphic in previous reports, the SSRs on each chromosome were sorted by motif type and repeat length for convenient selection. The primer sequences of all chromosome-specific SSRs were also made publicly available.

Conclusion

Chromosome-specific SSRs are efficient tools for chromosome identification by anchoring linkage groups to particular chromosomes during genetic mapping and are especially useful in mapping of qualitative-trait genes or quantitative trait loci with just a few markers. The SSRs reported here will facilitate a number of genetic and genomic studies in cotton, including construction of high-density genetic maps, positional gene cloning, fingerprinting, and genetic diversity and comparative evolutionary analyses among Gossypium species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1265-2) contains supplementary material, which is available to authorized users.     

Analysis of xyloglucan endotransglycosylase/hydrolase (XTH) genes from allotetraploid (Gossypium hirsutum) cotton and its diploid progenitors expressed during fiber elongation

Multiple cellular pathways have been shown to be involved during fiber initiation and elongation stages in the cultivated allotetraploid cotton (Gossypium hirsutum). The cell wall enzymes xyloglucan endotransglycosylase/hydrolases (XTH) have been reported to be associated with the biosynthesis of the cell wall and the growth of cotton fibers, probably regulating the plasticity of the primary cell wall. Among various cotton fiber cDNAs found to be preferentially expressed in cotton fibers, a xyloglucan endotransglycosylase (XTH) cDNA was significantly up-regulated during the elongation stage of cotton fiber development. In the present study, we isolated and characterized genomic clones encoding cotton XTH from cultivated cotton (Gossypium hirsutum) and its diploid progenitors (Gossypium arboreum and Gossypium raimondii), designated GhXTH1-1, GhXTH1-2, GaXTH1 and GrXTH, respectively. In addition, we isolated and characterized, by in silico methods, the putative promoter of XTH1 from Gossypium hirsutum. Sequence analysis revealed more than 50% homology to XTH's at the protein level. DNA gel blot hybridization indicated that at least two copies of GhXTH1 are present in Gossypium hirsutum whereas the diploid progenitor species Gossypium arboreum and Gossypium raimondii has only a single copy. Quantitative real-time PCR and high-resolution melting experiments indicated that in Gossypium hirsutum cultivars, in cotton fibers during early stages of fiber elongation specifically expressing only the GhXTH1-1 gene and expression levels of GhXTH1-1 in fibers varies among cultivars differing in fiber percentage and fiber length.     

Characterization of the nucleotide sequence of a polyubiquitin gene (PUBC1) from Arabian camel, Camelus dromedarius

Molecular amplification and sequencing of genomic DNA that encodes camel polyubiquitin (PUBC1) was performed by a polymerase chain reaction (PCR) using various sets of primers. The amplification generated a number of DNA fragments, which were sequenced and compared with the polyubiquitin coding sequences of various species. One DNA fragment that conformed to 325 bp was found to be 95 and 88% homologous to the sequences of human polyubiquitin B and C, respectively. The DNA translated into 108 amino acids that corresponded to two fused units of ubiquitin with no intervening sequence, which indicates that it is a polyubiquitin and contains at least two units of ubiquitin. Although, variations were found in the nucleotide sequence when compared to those of other species, the amino acid sequence was 100% homologous to the polyubiquitin sequences of humans, mice, and rats. This is the first report of the polyubiquitin DNA coding sequence and its corresponding amino acid sequence from camels, amplified using direct genomic DNA preparations.     

Characterization of the Brassinosteroid Insensitive 1 Genes of Cotton

Suppression of brassinosteroid (BR) biosynthesis in cotton ovules by treatment with brassinazole inhibits fiber formation, indicating that BR plays an important role in cotton fiber development. Plant responses to brassinosteroids (BR) are mediated through a plasma membrane-bound leucine-rich repeat (LRR) receptor-like protein kinase known as BRI1. Mutations in the BRI1 genes of several species result in dwarfed plants with reduced sensitivity to BR. A single expressed sequence tag (EST) from cotton with strong sequence similarity to Arabidopsis BRI1 ( AtBRI1 ) was identified in a search of publicly available databases. With this EST as a starting point, full-length cDNAs and genomic coding sequences from upland cotton ( Gossypium hirsutum ) BRI1 ( GhBRI1 ) were obtained and characterized. Ectopic expression of this coding sequence in BR-insensitive Arabidopsis plants resulted in recovery of normal growth indicating that GhBRI1 is a functional homologue of AtBRI1. G. hirsutum is an allotetraploid (AADD) derived from diploid ancestors. Analysis of several GhBRI1 cDNAs showed two distinct sequences indicating the presence of two GhBRI1 genes, denoted GhBRI1-1 and GhBRI1-2. Sequence comparisons between these GhBRI1 coding sequences and those from related A and D genome diploid Gossypium species ( G. arboreum and G. thurberi ) indicated that GhBRI1-1 is likely to the A sub-genome orthologue while GhBRI1-2 is from the D sub-genome.     

Cotton genome mapping with new microsatellites from Acala ‘Maxxa’ BAC-ends

Fine mapping and positional cloning will eventually improve with the anchoring of additional markers derived from genomic clones such as BACs. From 2,603 new BAC-end genomic sequences from Gossypium hirsutum Acala ‘Maxxa’, 1,316 PCR primer pairs (designated as MUSB) were designed to flank microsatellite or simple sequence repeat motif sequences. Most (1164 or 88%) MUSB primer pairs successfully amplified DNA from three species of cotton with an average of three amplicons per marker and 365 markers (21%) were polymorphic between G. hirsutum and G. barbadense. An interspecific RIL population developed from the above two entries was used to map 433 marker loci and 46 linkage groups with a genetic distance of 2,126.3 cM covering approximately 45% of the cotton genome and an average distance between two loci of 4.9 cM. Based on genome-specific chromosomes identified in G. hirsutum tetraploid (A and D), 56.9% of the coverage was located on the A subgenome while 39.7% was assigned to the D subgenome in the genetic map, suggesting that the A subgenome may be more polymorphic and recombinationally active than originally thought. The linkage groups were assigned to 23 of the 26 chromosomes. This is the first genetic map in which the linkage groups A01 and A02/D03 have been assigned to specific chromosomes. In addition the MUSB-derived markers from BAC-end sequences markers allows fine genetic and QTL mapping of important traits and for the first time provides reconciliation of the genetic and physical maps. Limited QTL analyses suggested that loci on chromosomes 2, 3, 12, 15 and 18 may affect variation in fiber quality traits. The original BAC clones containing the newly mapped MUSB that tag the QTLs provide critical DNA regions for the discovery of gene sequences involved in biological processes such as fiber development and pest resistance in cotton. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Development of one set of chromosome-specific microsatellite-containing BACs and their physical mapping in Gossypium hirsutum L.

Fluorescence in situ hybridization (FISH), using bacterial artificial chromosome (BAC) clone as probe, is a reliable cytological technique for chromosome identification. It has been used in many plants, especially in those containing numerous small chromosomes. We previously developed eight chromosome-specific BAC clones from tetraploid cotton, which were used as excellent cytological markers for chromosomes identification. Here, we isolated the other chromosome-specific BAC clones to make a complete set for the identification of all 26 chromosome-pairs by this technology in tetraploid cotton (Gossypium hirsutum L.). This set of BAC markers was demonstrated to be useful to assign each chromosome to a genetic linkage group unambiguously. In addition, these BAC clones also served as convenient and reliable landmarks for establishing physical linkage with unknown targeted sequences. Moreover, one BAC containing an EST, with high sequence similarity to a G. hirsutum ethylene-responsive element-binding factor was located physically on the long arm of chromosome A7 with the help of a chromosome-A7-specific BAC FISH marker. Comparative analysis of physical marker positions in the chromosomes by BAC-FISH and genetic linkage maps demonstrated that most of the 26 BAC clones were localized close to or at the ends of their respective chromosomes, and indicated that the recombination active regions of cotton chromosomes are primarily located in the distal regions. This technology also enables us to make associations between chromosomes and their genetic linkage groups and re-assign each chromosome according to the corresponding genetic linkage group. This BAC clones and BAC-FISH technology will be useful for us to evaluate grossly the degree to which a linkage map provides adequate coverage for developing a saturated genetic map, and provides a powerful resource for cotton genomic researches.     

QTL mapping in A-genome diploid Asiatic cotton and their congruence analysis with AD-genome tetraploid cotton in genus Gossypium

Asiatic cotton(Gossypium arboreum L.) is an Old World cultivated cotton species.The sinense race was planted extensively in China.Due to the advances in spinning technology during the last century,the species was replaced by the New World allotetraploid cotton G.hirsutum L.Gossypium arboreum is still grown in India and Pakistan and also used as an elite in current cotton breeding programs.In addition,G.arboreum serves as a model for genomic research in Gossypium.In the present study,we generated an A-genome diploid cotton intraspecific genetic map including 264 SSR loci with three morphological markers mapped to 1 3 linkage groups.The map spans 2,508.71 cM with an average distance of 9.4 cM between adjacent loci.A population containing 1 76 F2:3 families was used to perform quantitative trait loci(QTL)mapping for 17 phenotypes using Multiple QTL Model(MQM)of MapQTL ver 5.0.Overall,108 QTLs were detected on 13 chromosomes.Thirty-one QTLs for yield and its components were detected in the F2 population.Forty-one QTLs for yield and its components were detected in the F2:3 families with a total of 43 QTLs for fiber qualities.Two QTLs for seed cotton weight/plant and lint index and three QTLs for seed index were consistently detected both in F2 and F2:3.Most QTLs for fiber qualities and yields were located at the same interval or neighboring intervals.These results indicated that the negative correlation between fiber qualities and yield traits may result from either pleiotropic effect of one gene or linkage effects of multiple closely linked genes.     

Molecular analysis of resveratrol synthase. cDNA, genomic clones and relationship with chalcone synthase

Resveratrol synthase (RS), a key enzyme in biosynthesis of stilbene-type phytoalexins, catalyzes the formation of resveratrol from coumaroyl-CoA and malonyl-CoA. Two cDNA clones, pGSC1 and pGSC2, have been isolated from cDNA libraries established with poly(A)-rich RNA from peanut (Arachis hypogaea) cell cultures specifically induced for RS. These cDNAs were used to identify two genomic clones (pGSG10 and pGSG11). Sequence analysis shows that the two clones overlap in a large stretch of nearly identical sequences, and that pGSG10 contains the 5' and pGSG11 the 3' end of RS genes. The sequences reveal a single intron, and the size of the predicted protein is 42.7 kDa, in close agreement with that observed in polyacrylamide gels (43 kDa). Chalcone synthase (CHS), a key enzyme of flavonoid biosynthesis, utilizes the same substrates as RS, but the product is different (naringenin chalcone). Comparison of RS with CHS consensus sequences shows that the two genes are related. Homology extends throughout the coding region, and the intron in RS is at the same position as a conserved intron in CHS. However, RS reveals a substantial number of amino acid differences to CHS in positions highly conserved in all CHS enzymes. It is proposed that the two proteins possess a common scaffold necessary for binding of the substrates and the type of enzyme reaction, and that the differences are responsible for the formation of different products.