TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline.

Genotyping by sequencing (GBS) is a next generation sequencing based method that takes advantage of reduced representation to enable high throughput genotyping of large numbers of individuals at a large number of SNP markers. The relatively straightforward, robust, and cost-effective GBS protocol is currently being applied in numerous species by a large number of researchers. Herein we describe a bioinformatics pipeline, TASSEL-GBS, designed for the efficient processing of raw GBS sequence data into SNP genotypes. The TASSEL-GBS pipeline successfully fulfills the following key design criteria: (1) Ability to run on the modest computing resources that are typically available to small breeding or ecological research programs, including desktop or laptop machines with only 8-16 GB of RAM, (2) Scalability from small to extremely large studies, where hundreds of thousands or even millions of SNPs can be scored in up to 100,000 individuals (e.g., for large breeding programs or genetic surveys), and (3) Applicability in an accelerated breeding context, requiring rapid turnover from tissue collection to genotypes. Although a reference genome is required, the pipeline can also be run with an unfinished "pseudo-reference" consisting of numerous contigs. We describe the TASSEL-GBS pipeline in detail and benchmark it based upon a large scale, species wide analysis in maize (Zea mays), where the average error rate was reduced to 0.0042 through application of population genetic-based SNP filters. Overall, the GBS assay and the TASSEL-GBS pipeline provide robust tools for studying genomic diversity.

Grapevine powdery mildew resistance and susceptibility loci identified on a high-resolution SNP map.

Improved efficacy and durability of powdery mildew resistance can be enhanced via knowledge of the genetics of resistance and susceptibility coupled with the development of high-resolution maps to facilitate the stacking of multiple resistance genes and other desirable traits. We studied the inheritance of powdery mildew (Erysiphe necator) resistance and susceptibility of wild Vitis rupestris B38 and cultivated V. vinifera 'Chardonnay', finding evidence for quantitative variation. Molecular markers were identified using genotyping-by-sequencing, resulting in 16,833 single nucleotide polymorphisms (SNPs) based on alignment to the V. vinifera 'PN40024' reference genome sequence. With an average density of 36 SNPs/Mbp and uniform coverage of the genome, this 17K set was used to identify 11 SNPs on chromosome 7 associated with a resistance locus from V. rupestris B38 and ten SNPs on chromosome 9 associated with a locus for susceptibility from 'Chardonnay' using single marker association and linkage disequilibrium analysis. Linkage maps for V. rupestris B38 (1,146 SNPs) and 'Chardonnay' (1,215 SNPs) were constructed and used to corroborate the 'Chardonnay' locus named Sen1 (Susceptibility to Erysiphe necator 1), providing the first insight into the genetics of susceptibility to powdery mildew from V. vinifera. The identification of markers associated with a susceptibility locus in a V. vinifera background can be used for negative selection among breeding progenies. This work improves our understanding of the nature of powdery mildew resistance in V. rupestris B38 and 'Chardonnay', while applying next-generation sequencing tools to advance grapevine genomics and breeding.

Population genomic and genome-wide association studies of agroclimatic traits in sorghum.

Accelerating crop improvement in sorghum, a staple food for people in semiarid regions across the developing world, is key to ensuring global food security in the context of climate change. To facilitate gene discovery and molecular breeding in sorghum, we have characterized ~265,000 single nucleotide polymorphisms (SNPs) in 971 worldwide accessions that have adapted to diverse agroclimatic conditions. Using this genome-wide SNP map, we have characterized population structure with respect to geographic origin and morphological type and identified patterns of ancient crop diffusion to diverse agroclimatic regions across Africa and Asia. To better understand the genomic patterns of diversification in sorghum, we quantified variation in nucleotide diversity, linkage disequilibrium, and recombination rates across the genome. Analyzing nucleotide diversity in landraces, we find evidence of selective sweeps around starch metabolism genes, whereas in landrace-derived introgression lines, we find introgressions around known height and maturity loci. To identify additional loci underlying variation in major agroclimatic traits, we performed genome-wide association studies (GWAS) on plant height components and inflorescence architecture. GWAS maps several classical loci for plant height, candidate genes for inflorescence architecture. Finally, we trace the independent spread of multiple haplotypes carrying alleles for short stature or long inflorescence branches. This genome-wide map of SNP variation in sorghum provides a basis for crop improvement through marker-assisted breeding and genomic selection.

Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism.

Powdery mildews are phytopathogens whose growth and reproduction are entirely dependent on living plant cells. The molecular basis of this life-style, obligate biotrophy, remains unknown. We present the genome analysis of barley powdery mildew, Blumeria graminis f.sp. hordei (Blumeria), as well as a comparison with the analysis of two powdery mildews pathogenic on dicotyledonous plants. These genomes display massive retrotransposon proliferation, genome-size expansion, and gene losses. The missing genes encode enzymes of primary and secondary metabolism, carbohydrate-active enzymes, and transporters, probably reflecting their redundancy in an exclusively biotrophic life-style. Among the 248 candidate effectors of pathogenesis identified in the Blumeria genome, very few (less than 10) define a core set conserved in all three mildews, suggesting that most effectors represent species-specific adaptations.