Patrick E. McGuire

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies

BackgroundThe size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination.ResultsWe develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome.ConclusionsIn addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.

Genome comparisons reveal a dominant mechanism of chromosome number reduction in grasses and accelerated genome evolution in Triticeae

Single-nucleotide polymorphism was used in the construction of an expressed sequence tag map of Aegilops tauschii, the diploid source of the wheat D genome. Comparisons of the map with the rice and sorghum genome sequences revealed 50 inversions and translocations; 2, 8, and 40 were assigned respectively to the rice, sorghum, and Ae. tauschii lineages, showing greatly accelerated genome evolution in the large Triticeae genomes. The reduction of the basic chromosome number from 12 to 7 in the Triticeae has taken place by a process during which an entire chromosome is inserted by its telomeres into a break in the centromeric region of another chromosome. The original centromere–telomere polarity of the chromosome arms is maintained in the new chromosome. An intrachromosomal telomere–telomere fusion resulting in a pericentric translocation of a chromosome segment or an entire arm accompanied or preceded the chromosome insertion in some instances. Insertional dysploidy has been recorded in three grass subfamilies and appears to be the dominant mechanism of basic chromosome number reduction in grasses. A total of 64% and 66% of Ae. tauschii genes were syntenic with sorghum and rice genes, respectively. Synteny was reduced in the vicinity of the termini of modern Ae. tauschii chromosomes but not in the vicinity of the ancient termini embedded in the Ae. tauschii chromosomes, suggesting that the dependence of synteny erosion on gene location along the centromere–telomere axis either evolved recently in the Triticeae phylogenetic lineage or its evolution was recently accelerated.

A 4-gigabase physical map unlocks the structure and evolution of the complex genome of Aegilops tauschii, the wheat D-genome progenitor

The current limitations in genome sequencing technology require the construction of physical maps for high-quality draft sequences of large plant genomes, such as that of Aegilops tauschii, the wheat D-genome progenitor. To construct a physical map of the Ae. tauschii genome, we fingerprinted 461,706 bacterial artificial chromosome clones, assembled contigs, designed a 10K Ae. tauschii Infinium SNP array, constructed a 7,185-marker genetic map, and anchored on the map contigs totaling 4.03 Gb. Using whole genome shotgun reads, we extended the SNP marker sequences and found 17,093 genes and gene fragments. We showed that collinearity of the Ae. tauschii genes with Brachypodium distachyon, rice, and sorghum decreased with phylogenetic distance and that structural genome evolution rates have been high across all investigated lineages in subfamily Pooideae, including that of Brachypodieae. We obtained additional information about the evolution of the seven Triticeae chromosomes from 12 ancestral chromosomes and uncovered a pattern of centromere inactivation accompanying nested chromosome insertions in grasses. We showed that the density of noncollinear genes along the Ae. tauschii chromosomes positively correlates with recombination rates, suggested a cause, and showed that new genes, exemplified by disease resistance genes, are preferentially located in high-recombination chromosome regions.

Annotation-based genome-wide SNP discovery in the large and complex Aegilops tauschii genome using next-generation sequencing without a reference genome sequence

BackgroundMany plants have large and complex genomes with an abundance of repeated sequences. Many plants are also polyploid. Both of these attributes typify the genome architecture in the tribe Triticeae, whose members include economically important wheat, rye and barley. Large genome sizes, an abundance of repeated sequences, and polyploidy present challenges to genome-wide SNP discovery using next-generation sequencing (NGS) of total genomic DNA by making alignment and clustering of short reads generated by the NGS platforms difficult, particularly in the absence of a reference genome sequence.ResultsAn annotation-based, genome-wide SNP discovery pipeline is reported using NGS data for large and complex genomes without a reference genome sequence. Roche 454 shotgun reads with low genome coverage of one genotype are annotated in order to distinguish single-copy sequences and repeat junctions from repetitive sequences and sequences shared by paralogous genes. Multiple genome equivalents of shotgun reads of another genotype generated with SOLiD or Solexa are then mapped to the annotated Roche 454 reads to identify putative SNPs. A pipeline program package, AGSNP, was developed and used for genome-wide SNP discovery in Aegilops tauschii- the diploid source of the wheat D genome, and with a genome size of 4.02 Gb, of which 90% is repetitive sequences. Genomic DNA of Ae. tauschii accession AL8/78 was sequenced with the Roche 454 NGS platform. Genomic DNA and cDNA of Ae. tauschii accession AS75 was sequenced primarily with SOLiD, although some Solexa and Roche 454 genomic sequences were also generated. A total of 195,631 putative SNPs were discovered in gene sequences, 155,580 putative SNPs were discovered in uncharacterized single-copy regions, and another 145,907 putative SNPs were discovered in repeat junctions. These SNPs were dispersed across the entire Ae. tauschii genome. To assess the false positive SNP discovery rate, DNA containing putative SNPs was amplified by PCR from AL8/78 and AS75 and resequenced with the ABI 3730 xl. In a sample of 302 randomly selected putative SNPs, 84.0% in gene regions, 88.0% in repeat junctions, and 81.3% in uncharacterized regions were validated.ConclusionAn annotation-based genome-wide SNP discovery pipeline for NGS platforms was developed. The pipeline is suitable for SNP discovery in genomic libraries of complex genomes and does not require a reference genome sequence. The pipeline is applicable to all current NGS platforms, provided that at least one such platform generates relatively long reads. The pipeline package, AGSNP, and the discovered 497,118 Ae. tauschii SNPs can be accessed at (http://avena.pw.usda.gov/wheatD/agsnp.shtml).

Synteny perturbations between wheat homoeologous chromosomes caused by locus duplications and deletions correlate with recombination rates

Loci detected by Southern blot hybridization of 3,977 expressed sequence tag unigenes were mapped into 159 chromosome bins delineated by breakpoints of a series of overlapping deletions. These data were used to assess synteny levels along homoeologous chromosomes of the wheat A, B, and D genomes, in relation to both bin position on the centromere-telomere axis and the gradient of recombination rates along chromosome arms. Synteny level decreased with the distance of a chromosome region from the centromere. It also decreased with an increase in recombination rates along the average chromosome arm. There were twice as many unique loci in the B genome than in the A and D genomes, and synteny levels between the B genome chromosomes and the A and D genome homoeologues were lower than those between the A and D genome homoeologues. These differences among the wheat genomes were attributed to differences in the mating systems of wheat diploid ancestors. Synteny perturbations were characterized in 31 paralogous sets of loci with perturbed synteny. Both insertions and deletions of loci were detected and both preferentially occurred in high recombination regions of chromosomes.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes.

BackgroundA genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.ResultsNucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.ConclusionsIn a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Reconstruction of the Synthetic W7984 × Opata M85 wheat reference population

Reference populations are valuable resources in genetics studies for determining marker order, marker selection, trait mapping, construction of large-insert libraries, cross-referencing marker platforms, and genome sequencing. Reference populations can be propagated indefinitely, they are polymorphic and have normal segregation. Described are two new reference populations who share the same parents of the original wheat reference population Synthetic W7984 (Altar84/ Aegilops tauschii (219) CIGM86.940) x Opata M85, an F(1)-derived doubled haploid population (SynOpDH) of 215 inbred lines and a recombinant inbred population (SynOpRIL) of 2039 F(6) lines derived by single-plant self-pollinations. A linkage map was constructed for the SynOpDH population using 1446 markers. In addition, a core set of 42 SSR markers was genotyped on SynOpRIL. A new approach to identifying a core set of markers used a step-wise selection protocol based on polymorphism, uniform chromosome distribution, and reliability to create nested sets starting with one marker per chromosome, followed by two, four, and six. It is suggested that researchers use these markers as anchors for all future mapping projects to facilitate cross-referencing markers and chromosome locations. To enhance this public resource, researchers are strongly urged to validate line identities and deposit their data in GrainGenes so that others can benefit from the accumulated information.

Biodiversity in Agriculture: Domestication, Evolution, and Sustainability

List of contributors Foreword Bruce D. Smith Acknowledgments Introduction. The domestication of plants and animals: ten unanswered questions Paul Gepts, Robert Bettinger, Stephen Brush, Ardeshir Damania, Thomas Famula, Patrick McGuire and Calvin Qualset 1. The local origins of domestication Jared Diamond Part I. Early Steps in Agricultural Domestication: 2. Evolution of agro-ecosystems: biodiversity, origins and differential development David R. Harris 3. From foraging to farming in western and eastern Asia Ofer Bar-Yosef 4. Predomestic cultivation during the late Pleistocene and early Holocene in the Northern Levant George Willcox 5. New archaeobotanical information on plant domestication from macro-remains: tracking the evolution of domestication syndrome traits Dorian Q. Fuller 6. New archaeobotanical information on early cultivation and plant domestication involving microplant remains Dolores R. Piperno 7. How and why did agriculture spread? Peter Bellwood 8. California Indian proto-agriculture: its characterization and legacy M. Kat Anderson and Eric Wohlgemuth Part II. Domestication of Animals and Impacts on Humans: 9. Pathways to animal domestication Melinda A. Zeder 10. Genetics of animal domestication Leif Andersson 11. Genome-wide approaches for the study of dog domestication Bridgett M. vonHoldt, Melissa M. Gray and Robert K. Wayne 12. Malaria and rickets represent selective forces for the convergent evolution of adult lactase persistence Loren Cordain, Matthew S. Hickey and Kami Kim Part III. Issues in Plant Domestication: 13. The dynamics of rice domestication: a balance between gene flow and genetic isolation Susan R. McCouch, Michael J. Kovach, Megan Sweeney, Hui Jiang and Mande Semon 14. Domestication of lima beans: a new look at an old problem M. I. Chacon S., J. R. Motta-Aldana, M. L. Serrano S. and D. G. Debouck 15. Genetic characterization of cassava (Manihot esculenta Crantz) and yam (Dioscorea trifida L.) landraces in swidden agriculture systems in Brazil Elizabeth A. Veasey, Eduardo A. Bressan, Marcos V. B. M. Siqueira, Aline Borges, Jurema R. Queiroz-Silva, Kayo J. C. Pereira, Gustavo H. Recchia and Lin Chau Ming 16. Pigeonpea - from an orphan to a leader in food legumes C. L. Laxmipathi Gowda, K. B. Saxena, R. K. Srivastava, H. D. Upadhyaya and S. N. Silim Part IV. Traditional Management of Biodiversity: 17. Ecological approaches to crop domestication D. B. McKey, M. Elias, B. Pujol and A. Duputie 18. Agrobiodiversity shifts on three continents since Vavilov and Harlan: assessing causes, processes and implications for food security Gary Paul Nabhan, Ken Wilson, Ogonazar Aknazarov, Karim-Aly Kassam, Laurie Monti, David Cavagnaro, Shawn Kelly, Tai Johnson and Ferrell Sekacucu 19. Indigenous peoples conserving, managing, and creating biodiversity Jan Salick 20. Land architecture in the Maya lowlands: implications for sustainability B. L. Turner II and Deborah Lawrence 21. Agrobiodiversity and water resources in agricultural landscape evolution (Andean Valley irrigation, Bolivia, 1986 to 2008) Karl S. Zimmerer Part V. Uses of Biodiversity and New and Future Domestications: 22. Participatory domestication of indigenous fruit and nut trees: new crops for sustainable agriculture in developing countries Roger R. B. Leakey 23. The introduction and dispersal of Vitis vinifera into California: a case study of the interaction of man, plants, economics, and environment James Lapsley 24. Genetic resources of yeast and other micro-organisms Charles W. Bamforth 25. Biodiversity of native bees and crop pollination with emphasis on California Robbin W. Thorp 26. Aquaculture, the next wave of domestication Dennis Hedgecock 27. Genetic sustainability and biodiversity: challenges to the California dairy industry Juan F. Medrano Index.

Genome sequence of the progenitor of the wheat D genome Aegilops tauschii

Aegilops tauschii is the diploid progenitor of the D genome of hexaploid wheat (Triticum aestivum, genomes AABBDD) and an important genetic resource for wheat. The large size and highly repetitive nature of the Ae. tauschii genome has until now precluded the development of a reference-quality genome sequence. Here we use an array of advanced technologies, including ordered-clone genome sequencing, whole-genome shotgun sequencing, and BioNano optical genome mapping, to generate a reference-quality genome sequence for Ae. tauschii ssp. strangulata accession AL8/78, which is closely related to the wheat D genome. We show that compared to other sequenced plant genomes, including a much larger conifer genome, the Ae. tauschii genome contains unprecedented amounts of very similar repeated sequences. Our genome comparisons reveal that the Ae. tauschii genome has a greater number of dispersed duplicated genes than other sequenced genomes and its chromosomes have been structurally evolving an order of magnitude faster than those of other grass genomes. The decay of colinearity with other grass genomes correlates with recombination rates along chromosomes. We propose that the vast amounts of very similar repeated sequences cause frequent errors in recombination and lead to gene duplications and structural chromosome changes that drive fast genome evolution.

Sequence of the Sugar Pine Megagenome.

Until very recently, complete characterization of the megagenomes of conifers has remained elusive. The diploid genome of sugar pine (Pinus lambertiana Dougl.) has a highly repetitive, 31 billion bp genome. It is the largest genome sequenced and assembled to date, and the first from the subgenus Strobus, or white pines, a group that is notable for having the largest genomes among the pines. The genome represents a unique opportunity to investigate genome “obesity” in conifers and white pines. Comparative analysis of P. lambertiana and P. taeda L. reveals new insights on the conservation, age, and diversity of the highly abundant transposable elements, the primary factor determining genome size. Like most North American white pines, the principal pathogen of P. lambertiana is white pine blister rust (Cronartium ribicola J.C. Fischer ex Raben.). Identification of candidate genes for resistance to this pathogen is of great ecological importance. The genome sequence afforded us the opportunity to make substantial progress on locating the major dominant gene for simple resistance hypersensitive response, Cr1. We describe new markers and gene annotation that are both tightly linked to Cr1 in a mapping population, and associated with Cr1 in unrelated sugar pine individuals sampled throughout the species’ range, creating a solid foundation for future mapping. This genomic variation and annotated candidate genes characterized in our study of the Cr1 region are resources for future marker-assisted breeding efforts as well as for investigations of fundamental mechanisms of invasive disease and evolutionary response.