Dan Milbourne

Construction of Reference Chromosome-Scale Pseudomolecules for Potato: Integrating the Potato Genome with Genetic and Physical Maps

The genome of potato, a major global food crop, was recently sequenced. The work presented here details the integration of the potato reference genome (DM) with a new sequence-tagged site marker−based linkage map and other physical and genetic maps of potato and the closely related species tomato. Primary anchoring of the DM genome assembly was accomplished by the use of a diploid segregating population, which was genotyped with several types of molecular genetic markers to construct a new ~936 cM linkage map comprising 2469 marker loci. In silico anchoring approaches used genetic and physical maps from the diploid potato genotype RH89-039-16 (RH) and tomato. This combined approach has allowed 951 superscaffolds to be ordered into pseudomolecules corresponding to the 12 potato chromosomes. These pseudomolecules represent 674 Mb (~93%) of the 723 Mb genome assembly and 37,482 (~96%) of the 39,031 predicted genes. The superscaffold order and orientation within the pseudomolecules are closely collinear with independently constructed high density linkage maps. Comparisons between marker distribution and physical location reveal regions of greater and lesser recombination, as well as regions exhibiting significant segregation distortion. The work presented here has led to a greatly improved ordering of the potato reference genome superscaffolds into chromosomal “pseudomolecules”.

Sequencing the Potato Genome: Outline and First Results to Come from the Elucidation of the Sequence of the World's Third Most Important Food Crop

Potato is a member of the Solanaceae, a plant family that includes several other economically important species, such as tomato, eggplant, petunia, tobacco and pepper. The Potato Genome Sequencing Consortium (PGSC) aims to elucidate the complete genome sequence of potato, the third most important food crop in the world. The PGSC is a collaboration between 13 research groups from China, India, Poland, Russia, the Netherlands, Ireland, Argentina, Brazil, Chile, Peru, USA, New Zealand and the UK. The potato genome consists of 12 chromosomes and has a (haploid) length of approximately 840 million base pairs, making it a medium-sized plant genome. The sequencing project builds on a diploid potato genomic bacterial artificial chromosome (BAC) clone library of 78000 clones, which has been fingerprinted and aligned into ~7000 physical map contigs. In addition, the BAC-ends have been sequenced and are publicly available. Approximately 30000 BACs are anchored to the Ultra High Density genetic map of potato, composed of 10000 unique AFLPTM markers. From this integrated genetic-physical map, between 50 to 150 seed BACs have currently been identified for every chromosome. Fluorescent in situ hybridization experiments on selected BAC clones confirm these anchor points. The seed clones provide the starting point for a BAC-by-BAC sequencing strategy. This strategy is being complemented by whole genome shotgun sequencing approaches using both 454 GS FLX and Illumina GA2 instruments. Assembly and annotation of the sequence data will be performed using publicly available and tailor-made tools. The availability of the annotated data will help to characterize germplasm collections based on allelic variance and to assist potato breeders to more fully exploit the genetic potential of potato.ResumenLa papa es un miembro de las Solanaceae, una familia de plantas que incluye varias otras especies económicamente importantes, tales como tomate, berenjena, petunia, tabaco y ají o chili. El consorcio de secuenciación del genoma de la papa (PGSC) tiene por objeto dilucidar la secuencia completa del genoma de la papa, el tercer cultivo alimentario más importante del mundo. El PGSC es una colaboración entre 13 grupos de investigación procedentes de China, India, Polonia, Rusia, los Países Bajos, Irlanda, Argentina, Brasil, Chile, Perú, EE.UU., Nueva Zelanda y el Reino Unido. El genoma de la papa consiste de 12 cromosomas y tiene una longitud (haploide) de aproximadamente 840 millones de pares de bases, por lo que es una planta con un genoma de tamaño mediano. El proyecto de secuenciación se basa en una biblioteca de 78000 clones de cromosoma artificial bacteriano genomico de papa diploide (BAC), del que se ha obtenido la huella genética y alineado en 7000 ~ contigs de mapa físico. Además, los extremos terminales BAC se han secuenciado y están a disposición del público. Aproximadamente 30000 BACS están anclados al mapa genético de ultra alta densidad de la papa, compuesto de 10000 marcadores AFLPTM únicos. De esta mapa genético-físico integrado, entre 50 a 150 semillas BACs han sido identificadas para cada cromosoma. Experimentos de hibridación in situ fluorescente en clones BAC selectos confirman estos puntos de anclaje. La clones semilla proveen el punto de partida para la estrategia de secuenciación de BAC a BAC. Esta estrategia se complementa con los enfoques de secuenciación escopeta del genoma completo usando los instrumentos 454 GS FLX e Illumina GA2. El ensamblaje y anotación de los datos de la secuencia será realizados mediante herramientas publicas disponibles y hechas a la medida. La disponibilidad de los datos anotados ayudarán a caracterizar las colecciones de germoplasma basándose en variación alélica y ayudará a los fitomejoradores de papa a explotar más plenamente el potencial genético de la papa.

Identification of coincident QTL for days to heading, spike length and spikelets per spike in Lolium perenne L.

Flowering time is a trait which has a major influence on the quality of forage. In addition, flowering and subsequent seed yields are important traits for seed production by grass breeders. In this study, we have identified quantitative trait loci (QTL) for flowering time and morphological traits of the flowering head in an F1 mapping population in Lolium perenne L (perennial ryegrass), a number of which have not previously been identified in L. perenne mapping studies. QTL for days to heading (DTH) were mapped in both outdoor and glasshouse experiments, revealing three and five QTL for DTH which explained 53% and 42% of the total phenotypic variation observed, respectively. Two QTL for DTH were detected in both environments, although they had contrasting relative magnitudes in each environment. One QTL for spike length and three QTL for spikelets per spike were also identified explaining, a total of 32 and 33% of the phenotypic variance, respectively. Furthermore, the QTL for spike length and spikelets per spike generally coincided with QTL for days to heading, implying co-ordinate regulation by underlying genes. Of particular interest was a region harbouring overlapping QTL for days to heading, spike length and spikelets per spike on the top of linkage group 4, containing the major QTL for spike length identified in this population.

Development of diagnostic markers for use in breeding potatoes resistant to Globodera pallida pathotype Pa2/3 using germplasm derived from Solanum tuberosum ssp. andigena CPC 2802

Quantitative resistance to Globodera pallida pathotype Pa2/3, originally derived from Solanum tuberosum ssp. andigena Commonwealth Potato Collection (CPC) accession 2802, is present in several potato cultivars and advanced breeding lines. One genetic component of this resistance, a large effect quantitative trait locus (QTL) on linkage group IV (which we have renamed GpaIVadgs) has previously been mapped in the tetraploid breeding line 12601ab1. In this study, we show that GpaIVadgs is also present in a breeding line called C1992/31 via genetic mapping in an F1 population produced by crossing C1992/31 with the G. pallida susceptible cultivar Record. C1992/31 is relatively divergent from 12601ab1, confirming that GpaIVadgs is an ideal target for marker-assisted selection in currently available germplasm. To generate markers exhibiting diagnostic potential for GpaIVadgs, three bacterial artificial chromosome clones were isolated from the QTL region, sequenced, and used to develop 15 primer sets generating single-copy amplicons, which were examined for polymorphisms exhibiting linkage to GpaIVadgs in C1992/31. Eight such polymorphisms were found. Subsequently, one insertion/deletion polymorphism, three single nucleotide polymorphisms and a specific allele of the microsatellite marker STM3016 were shown to exhibit diagnostic potential for the QTL in a panel of 37 potato genotypes, 12 with and 25 without accession CPC2082 in their pedigrees. STM3016 and one of the SNP polymorphisms, C237(119), were assayed in 178 potato genotypes, arising from crosses between C1992/31 and 16 G. pallida susceptible genotypes, undergoing selection in a commercial breeding programme. The results suggest that the diagnostic markers would most effectively be employed in MAS-based approaches to pyramid different resistance loci to develop cultivars exhibiting strong, durable resistance to G. pallida pathotype Pa2/3.

A hybrid next generation transcript sequencing-based approach to identify allelic and homeolog-specific single nucleotide polymorphisms in allotetraploid white clover

BackgroundWhite clover (Trifolium repens L.) is an allotetraploid species possessing two highly collinear ancestral sub-genomes. The apparent existence of highly similar homeolog copies for the majority of genes in white clover is problematic for the development of genome-based resources in the species. This is especially true for the development of genetic markers based on single nucleotide polymorphisms (SNPs), since it is difficult to distinguish between homeolog-specific and allelic variants. Robust methods for categorising single nucleotide variants as allelic or homeolog-specific in large transcript datasets are required. We illustrate one potential approach in this study.ResultsWe used 454-pyrosequencing sequencing to generate ~760,000 transcript sequences from an 8th generation white clover inbred line. These were assembled and partially annotated to yield a reference transcript set comprising 71,545 sequences. We subsequently performed Illumina sequencing on three further white clover samples, generating 14 million transcript reads from a mixed sample comprising 24 divergent white clover genotypes, and 50 million reads on two further eighth generation white clover inbred lines. Mapping these reads to the reference transcript set allowed us to develop a significant SNP resource for white clover, and to partition the SNPs from the inbred lines into categories reflecting allelic or homeolog-specific variation. The potential for using haplotype reconstruction and progenitor genome comparison to assign haplotypes to specific ancestral sub-genomes of white clover is demonstrated for sequences corresponding to genes encoding dehydration responsive element binding protein and acyl-coA oxidase.ConclusionsIn total, 208,854 independent SNPs in 31,715 reference sequences were discovered, approximately three quarters of which were categorised as representing allelic or homeolog-specific variation using two inbred lines. This represents a significant resource for white clover genomics and genetics studies. We discuss the potential to extend the analysis to identify a “core set” of ancestrally derived homeolog specific variants in white clover.

The genetic location of the self-incompatibility locus in white clover (Trifolium repens L.)

White clover (Trifolium repens L.) is a forage legume of considerable economic importance in temperate agricultural systems. It has a strong self-incompatibility system. The molecular basis of self-incompatibility in T. repens is unknown, but it is under the control of a single locus, which is expressed gametophytically. To locate the self-incompatibility locus (S locus) in T. repens, we carried out cross-pollination experiments in an F1 mapping population and constructed a genetic linkage map using amplified fragment length polymorphism and simple sequence repeat markers. As the first step in a map-based cloning strategy, we locate for the first time the S locus in T. repens on a genetic linkage map, on the homoeologous linkage group pair 1 (E), which is broadly syntenic to Medicago truncatula L. chromosome 1. On the basis of this syntenic relationship, the possibility that the S locus may or may not possess an S-RNase gene is discussed.

Breeding strategies for sustainable forage and turf grass improvement

The paper analyses actual trends in (European) ruminant agriculture and grassland based production systems. Consequences of reduced and/or zero grazing for grass breeding and grassland management are discussed. The impacts on ecoefficiency, recycling of minerals and ecosystem services are highlighted as well as the role of ley-arable farming. Special emphasis is on the potential use of tall fescue as a component of mixtures or as an interspecific cross. In grazed grassland, the role of white clover, the disease resistance and the nitrogen use efficiency of the grasses and the significance of biodiversity are considered. Based on an article published by Parsons et al. (2011) some reflections on the way ahead in grass and forage breeding are presented.

An ultra-high density genetic linkage map of perennial ryegrass (Lolium perenne) using genotyping by sequencing (GBS) based on a reference shotgun genome assembly

BACKGROUND AND AIMS High density genetic linkage maps that are extensively anchored to assembled genome sequences of the organism in question are extremely useful in gene discovery. To facilitate this process in perennial ryegrass (Lolium perenne L.), a high density single nucleotide polymorphism (SNP)- and presence/absence variant (PAV)-based genetic linkage map has been developed in an F2 mapping population that has been used as a reference population in numerous studies. To provide a reference sequence to which to align genotyping by sequencing (GBS) reads, a shotgun assembly of one of the grandparents of the population, a tenth-generation inbred line, was created using Illumina-based sequencing. METHODS The assembly was based on paired-end Illumina reads, scaffolded by mate pair and long jumping distance reads in the range of 3-40 kb, with >200-fold initial genome coverage. A total of 169 individuals from an F2 mapping population were used to construct PstI-based GBS libraries tagged with unique 4-9 nucleotide barcodes, resulting in 284 million reads, with approx. 1·6 million reads per individual. A bioinformatics pipeline was employed to identify both SNPs and PAVs. A core genetic map was generated using high confidence SNPs, to which lower confidence SNPs and PAVs were subsequently fitted in a straightforward binning approach. KEY RESULTS The assembly comprises 424 750 scaffolds, covering 1·11 Gbp of the 2·5 Gbp perennial ryegrass genome, with a scaffold N50 of 25 212 bp and a contig N50 of 3790 bp. It is available for download, and access to a genome browser has been provided. Comparison of the assembly with available transcript and gene model data sets for perennial ryegrass indicates that approx. 570 Mbp of the gene-rich portion of the genome has been captured. An ultra-high density genetic linkage map with 3092 SNPs and 7260 PAVs was developed, anchoring just over 200 Mb of the reference assembly. CONCLUSIONS The combined genetic map and assembly, combined with another recently released genome assembly, represent a significant resource for the perennial ryegrass genetics community.

Markers associated with heading and aftermath heading in perennial ryegrass full-sib families

BackgroundHeading and aftermath heading are important traits in perennial ryegrass because they impact forage quality. So far, genome-wide association analyses in this major forage species have only identified a small number of genetic variants associated with heading date that overall explained little of the variation. Some possible reasons include rare alleles with large phenotypic affects, allelic heterogeneity, or insufficient marker density. We established a genome-wide association panel with multiple genotypes from multiple full-sib families. This ensured alleles were present at the frequency needed to have sufficient statistical power to identify associations. We genotyped the panel via partial genome sequencing and performed genome-wide association analyses with multi-year phenotype data collected for heading date, and aftermath heading.ResultsGenome wide association using a mixed linear model failed to identify any variants significantly associated with heading date or aftermath heading. Our failure to identify associations for these traits is likely due to the extremely low linkage disequilibrium we observed in this population. However, using single marker analysis within each full-sib family we could identify markers and genomic regions associated with heading and aftermath heading. Using the ryegrass genome we identified putative orthologs of key heading genes, some of which were located in regions of marker-trait associations.ConclusionGiven the very low levels of LD, genome wide association studies in perennial ryegrass populations are going to require very high SNP densities. Single marker analysis within full-sibs enabled us to identify significant marker-trait associations. One of these markers anchored proximal to a putative ortholog of TFL1, homologues of which have been shown to play a key role in continuous heading of some members of the rose family, Rosaceae.

Cadmium uptake and partitioning in potato (Solanum tuberosum L.) cultivars with different tuber-Cd concentration

Potatoes grown in soil with high Cd concentrations can accumulate high levels of Cd in the tubers. Although there is significant environmental variation involved in the trait of crop uptake of Cd, there are also distinctive cultivar differences. In order to understand this differential Cd accumulation mechanism, two potato cultivars were chosen that accumulate high and low levels of Cd in tubers. The patterns of Cd concentration, Cd content and dry weight accumulation of the two cultivars were examined at different stages of plant growth. The data suggest that differences in total Cd uptake and in Cd partitioning among organs are the mechanisms governing differential Cd-tuber accumulation in the two cultivars. The low tuber-Cd accumulator exhibited lower root-to-shoot and shoot-to-tuber translocation driven by higher root and shoot biomass that retained more Cd in roots and shoots, respectively, reducing its movement to the tubers. Higher remobilization and more efficient tuber loading was observed in the high tuber-Cd accumulator, indicating that remobilization of Cd from leaves to tubers was a major factor, not only in tuber-Cd loading, but also in the establishment of differential tuber-Cd levels. Regardless of cultivar differences, the concentration of Cd in the tuber was very low compared to that in other organs suggesting that, despite its high phloem mobility, Cd tends to be sequestered in the shoots.