Functional classification, genomic organization, putatively cis-acting regulatory elements, and relationship to quantitative trait loci, of sorghum genes with rhizome-enriched expression.

Rhizomes are organs of fundamental importance to plant competitiveness and invasiveness. We have identified genes expressed at substantially higher levels in rhizomes than other plant parts, and explored their functional categorization, genomic organization, regulatory motifs, and association with quantitative trait loci (QTLs) conferring rhizomatousness. The finding that genes with rhizome-enriched expression are distributed across a wide range of functional categories suggests some degree of specialization of individual members of many gene families in rhizomatous plants. A disproportionate share of genes with rhizome-enriched expression was implicated in secondary and hormone metabolism, and abiotic stimuli and development. A high frequency of unknown-function genes reflects our still limited knowledge of this plant organ. A putative oligosaccharyl transferase showed the highest degree of rhizome-specific expression, with several transcriptional or regulatory protein complex factors also showing high (but lesser) degrees of specificity. Inferred by the upstream sequences of their putative rice (Oryza sativa) homologs, sorghum (Sorghum bicolor) genes that were relatively highly expressed in rhizome tip tissues were enriched for cis-element motifs, including the pyrimidine box, TATCCA box, and CAREs box, implicating the gibberellins in regulation of many rhizome-specific genes. From cDNA clones showing rhizome-enriched expression, expressed sequence tags forming 455 contigs were plotted on the rice genome and aligned to QTL likelihood intervals for ratooning and rhizomatous traits in rice and sorghum. Highly expressed rhizome genes were somewhat enriched in QTL likelihood intervals for rhizomatousness or ratooning, with specific candidates including some of the most rhizome-specific genes. Some rhizomatousness and ratooning QTLs were shown to be potentially related to one another as a result of ancient duplication, suggesting long-term functional conservation of the underlying genes. Insight into genes and pathways that influence rhizome growth set the stage for genetic and/or exogenous manipulation of rhizomatousness, and for further dissection of the molecular evolution of rhizomatousness.

A global assembly of cotton ESTs.

Approximately 185,000 Gossypium EST sequences comprising >94,800,000 nucleotides were amassed from 30 cDNA libraries constructed from a variety of tissues and organs under a range of conditions, including drought stress and pathogen challenges. These libraries were derived from allopolyploid cotton (Gossypium hirsutum; A(T) and D(T) genomes) as well as its two diploid progenitors, Gossypium arboreum (A genome) and Gossypium raimondii (D genome). ESTs were assembled using the Program for Assembling and Viewing ESTs (PAVE), resulting in 22,030 contigs and 29,077 singletons (51,107 unigenes). Further comparisons among the singletons and contigs led to recognition of 33,665 exemplar sequences that represent a nonredundant set of putative Gossypium genes containing partial or full-length coding regions and usually one or two UTRs. The assembly, along with their UniProt BLASTX hits, GO annotation, and Pfam analysis results, are freely accessible as a public resource for cotton genomics. Because ESTs from diploid and allotetraploid Gossypium were combined in a single assembly, we were in many cases able to bioinformatically distinguish duplicated genes in allotetraploid cotton and assign them to either the A or D genome. The assembly and associated information provide a framework for future investigation of cotton functional and evolutionary genomics.

Organization and evolution of resistance gene analogs in peanut.

The scarcity of genetic polymorphism in Arachis hypogaea (peanut), as in other monophyletic polyploid species, makes it especially vulnerable to nematode, bacterial, fungal, and viral pathogens. Although no disease resistance genes have been cloned from peanut itself, the conserved motifs in cloned resistance genes from other plant species provide a means to isolate and analyze similar genes from peanut. To survey the number, diversity, evolutionary history, and genomic organization of resistance gene-like sequences in peanut, we isolated 234 resistance gene analogs (RGAs) by using primers designed from conserved regions of different classes of resistance genes including NBS-LRR, and LRR-TM classes. Phylogenetic and sequence analyses were performed to explore evolutionary relationships both among peanut RGAs and with orthologous genes from other plant taxa. Fifty-six overgos designed from the RGA sequences on the basis of their phyletic association were applied to a peanut BAC library; 736 hybridizing BAC clones were fingerprinted and contigs were formed in order to gain insights into the genomic organization of these genes. All the fingerprinting gels were blotted and screened with the respective overgos in order to verify the authenticity of the hits from initial screens, and to explore the physical organization of these genes in terms of both copy number and distribution in the genome. As a result, we identified 250 putative resistance gene loci. A correlation was found between the phyletic positions of the sequences and their physical locations. The BACs isolated here will serve as a valuable resource for future applications, such as map-based cloning, and will help improve our understanding of the evolution and organization of these genes in the peanut genome.

Comparative genomics of Gossypium and Arabidopsis: unraveling the consequences of both ancient and recent polyploidy.

Both ancient and recent polyploidy, together with post-polyploidization loss of many duplicated gene copies, complicates angiosperm comparative genomics. To explore an approach by which these challenges might be mitigated, genetic maps of extant diploid and tetraploid cottons (Gossypium spp.) were used to infer the approximate order of 3016 loci along the chromosomes of their hypothetical common ancestor. The inferred Gossypium gene order corresponded more closely than the original maps did to a similarly inferred ancestral gene order predating an independent paleopolyploidization (alpha) in Arabidopsis. At least 59% of the cotton map and 53% of the Arabidopsis transcriptome showed correspondence in multilocus gene arrangements based on one or both of two software packages (CrimeStatII, FISH). Genomic regions in which chromosome structural rearrangement has been rapid (obscuring gene order correspondence) have also been subject to greater divergence of individual gene sequences. About 26%-44% of corresponding regions involved multiple Arabidopsis or cotton chromosomes, in some cases consistent with known, more ancient, duplications. The genomic distributions of multiple-locus probes provided early insight into the consequences for chromosome structure of an ancient large-scale duplication in cotton. Inferences that mitigate the consequences of ancient duplications improve leveraging of genomic information for model organisms in the study of more complex genomes.

The repetitive landscape of the chicken genome.

Cot-based cloning and sequencing (CBCS) is a powerful tool for isolating and characterizing the various repetitive components of any genome, combining the established principles of DNA reassociation kinetics with high-throughput sequencing. CBCS was used to generate sequence libraries representing the high, middle, and low-copy fractions of the chicken genome. Sequencing high-copy DNA of chicken to about 2.7 x coverage of its estimated sequence complexity led to the initial identification of several new repeat families, which were then used for a survey of the newly released first draft of the complete chicken genome. The analysis provided insight into the diversity and biology of known repeat structures such as CR1 and CNM, for which only limited sequence data had previously been available. Cot sequence data also resulted in the identification of four novel repeats (Birddawg, Hitchcock, Kronos, and Soprano), two new subfamilies of CR1 repeats, and many elements absent from the chicken genome assembly. Multiple autonomous elements were found for a novel Mariner-like transposon, Galluhop, in addition to nonautonomous deletion derivatives. Phylogenetic analysis of the high-copy repeats CR1, Galluhop, and Birddawg provided insight into two distinct genome dispersion strategies. This study also exemplifies the power of the CBCS method to create representative databases for the repetitive fractions of genomes for which only limited sequence data is available.

An SNP resource for rice genetics and breeding based on subspecies indica and japonica genome alignments.

Dense coverage of the rice genome with polymorphic DNA markers is an invaluable tool for DNA marker-assisted breeding, positional cloning, and a wide range of evolutionary studies. We have aligned drafts of two rice subspecies, indica and japonica, and analyzed levels and patterns of genetic diversity. After filtering multiple copy and low quality sequence, 408,898 candidate DNA polymorphisms (SNPs/INDELs) were discerned between the two subspecies. These filters have the consequence that our data set includes only a subset of the available SNPs (in particular excluding large numbers of SNPs that may occur between repetitive DNA alleles) but increase the likelihood that this subset is useful: Direct sequencing suggests that 79.8% +/- 7.5% of the in silico SNPs are real. The SNP sample in our database is not randomly distributed across the genome. In fact, 566 rice genomic regions had unusually high (328 contigs/48.6 Mb/13.6% of genome) or low (237 contigs/64.7 Mb/18.1% of genome) polymorphism rates. Many SNP-poor regions were substantially longer than most SNP-rich regions, covering up to 4 Mb, and possibly reflecting introgression between the respective gene pools that may have occurred hundreds of years ago. Although 46.2% +/- 8.3% of the SNPs differentiate other pairs of japonica and indica genotypes, SNP rates in rice were not predictive of evolutionary rates for corresponding genes in another grass species, sorghum. The data set is freely available at http://www.plantgenome.uga.edu/snp.

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.

A 3347-locus genetic recombination map of sequence-tagged sites reveals features of genome organization, transmission and evolution of cotton (Gossypium).

We report genetic maps for diploid (D) and tetraploid (AtDt) Gossypium genomes composed of sequence-tagged sites (STS) that foster structural, functional, and evolutionary genomic studies. The maps include, respectively, 2584 loci at 1.72-cM ( approximately 600 kb) intervals based on 2007 probes (AtDt) and 763 loci at 1.96-cM ( approximately 500 kb) intervals detected by 662 probes (D). Both diploid and tetraploid cottons exhibit negative crossover interference; i.e., double recombinants are unexpectedly abundant. We found no major structural changes between Dt and D chromosomes, but confirmed two reciprocal translocations between At chromosomes and several inversions. Concentrations of probes in corresponding regions of the various genomes may represent centromeres, while genome-specific concentrations may represent heterochromatin. Locus duplication patterns reveal all 13 expected homeologous chromosome sets and lend new support to the possibility that a more ancient polyploidization event may have predated the A-D divergence of 6-11 million years ago. Identification of SSRs within 312 RFLP sequences plus direct mapping of 124 SSRs and exploration for CAPS and SNPs illustrate the "portability" of these STS loci across populations and detection systems useful for marker-assisted improvement of the world's leading fiber crop. These data provide new insights into polyploid evolution and represent a foundation for assembly of a finished sequence of the cotton genome.

A comparative phylogenetic approach for dating whole genome duplication events.

MOTIVATION: Whole genome duplications have played a major role in determining the structure of eukaryotic genomes. Current evidence revealing large blocks of duplicated chromatin yields new insights into the evolutionary history of species, but also presents a major challenge for researchers attempting to utilize comparative genomics techniques. Understanding the timing of duplication events relative to divergence among taxa is critical to accurate and comprehensive cross-species comparisons. RESULTS: We describe a large-scale approach to estimate the timing of duplication events in a phylogenetic context. The methodology has been previously utilized for analysis of Arabidopsis and Saccharomyces duplication events. This new implementation provides a more flexible and reusable framework for these analyses. Scripts written in the Python programming language drive a number of freely available bioinformatics programs, creating a no-cost tool for researchers. The usefulness of the approach is demonstrated through genome-scale analysis of Arabidopsis and Oryza (rice) duplications. AVAILABILITY: Software and documentation are freely available from http://plantgenome.agtec.uga.edu/bioinformatics/dating/

A high-density genetic recombination map of sequence-tagged sites for sorghum, as a framework for comparative structural and evolutionary genomics of tropical grains and grasses.

We report a genetic recombination map for Sorghum of 2512 loci spaced at average 0.4 cM ( approximately 300 kb) intervals based on 2050 RFLP probes, including 865 heterologous probes that foster comparative genomics of Saccharum (sugarcane), Zea (maize), Oryza (rice), Pennisetum (millet, buffelgrass), the Triticeae (wheat, barley, oat, rye), and Arabidopsis. Mapped loci identify 61.5% of the recombination events in this progeny set and reveal strong positive crossover interference acting across intervals of

Integration of Cot analysis, DNA cloning, and high-throughput sequencing facilitates genome characterization and gene discovery.

Cot-based sequence discovery represents a powerful means by which both low-copy and repetitive sequences can be selectively and efficiently fractionated, cloned, and characterized. Based upon the results of a Cot analysis, hydroxyapatite chromatography was used to fractionate sorghum (Sorghum bicolor) genomic DNA into highly repetitive (HR), moderately repetitive (MR), and single/low-copy (SL) sequence components that were consequently cloned to produce HRCot, MRCot, and SLCot genomic libraries. Filter hybridization (blotting) and sequence analysis both show that the HRCot library is enriched in sequences traditionally found in high-copy number (e.g., retroelements, rDNA, centromeric repeats), the SLCot library is enriched in low-copy sequences (e.g., genes and "nonrepetitive ESTs"), and the MRCot library contains sequences of moderate redundancy. The Cot analysis suggests that the sorghum genome is approximately 700 Mb (in agreement with previous estimates) and that HR, MR, and SL components comprise 15%, 41%, and 24% of sorghum DNA, respectively. Unlike previously described techniques to sequence the low-copy components of genomes, sequencing of Cot components is independent of expression and methylation patterns that vary widely among DNA elements, developmental stages, and taxa. High-throughput sequencing of Cot clones may be a means of "capturing" the sequence complexity of eukaryotic genomes at unprecedented efficiency.