Evaluation and Recommendations for Routine Genotyping Using Skim Whole Genome Re-sequencing in Canola.

Whole genome sequencing offers genome wide, unbiased markers, and inexpensive library preparation. With the cost of sequencing decreasing rapidly, many plant genomes of modest size are amenable to skim whole genome resequencing (skim WGR). The use of skim WGR in diverse sample sets without the use of imputation was evaluated in silico in 149 canola samples representative of global diversity. Fastq files with an average of 10x coverage of the reference genome were used to generate skim samples representing 0.25x, 0.5x, 1x, 2x, 3x, 4x, and 5x sequencing coverage. Applying a pre-defined list of SNPs versus de novo SNP discovery was evaluated. As skim WGR is expected to result in some degree of insufficient allele sampling, all skim coverage levels were filtered at a range of minimum read depths from a relaxed minimum read depth of 2 to a stringent read depth of 5, resulting in 28 list-based SNP sets. As a broad recommendation, genotyping pre-defined SNPs between 1x and 2x coverage with relatively stringent depth filtering is appropriate for a diverse sample set of canola due to a balance between marker number, sufficient accuracy, and sequencing cost, but depends on the intended application. This was experimentally examined in two sample sets with different genetic backgrounds: 1x coverage of 1,590 individuals from 84 Australian spring type four-parent crosses aimed at maximizing diversity as well as one commercial F1 hybrid, and 2x coverage of 379 doubled haploids (DHs) derived from a subset of the four-parent crosses. To determine optimal coverage in a simpler genetic background, the DH sample sequence coverage was further down sampled in silico. The flexible and cost-effective nature of the protocol makes it highly applicable across a range of species and purposes.

Exploitation of data from breeding programs supports rapid implementation of genomic selection for key agronomic traits in perennial ryegrass.

KEY MESSAGE: Exploitation of data from a ryegrass breeding program has enabled rapid development and implementation of genomic selection for sward-based biomass yield with a twofold-to-threefold increase in genetic gain. Genomic selection, which uses genome-wide sequence polymorphism data and quantitative genetics techniques to predict plant performance, has large potential for the improvement in pasture plants. Major factors influencing the accuracy of genomic selection include the size of reference populations, trait heritability values and the genetic diversity of breeding populations. Global diversity of the important forage species perennial ryegrass is high and so would require a large reference population in order to achieve moderate accuracies of genomic selection. However, diversity of germplasm within a breeding program is likely to be lower. In addition, de novo construction and characterisation of reference populations are a logistically complex process. Consequently, historical phenotypic records for seasonal biomass yield and heading date over a 18-year period within a commercial perennial ryegrass breeding program have been accessed, and target populations have been characterised with a high-density transcriptome-based genotyping-by-sequencing assay. Ability to predict observed phenotypic performance in each successive year was assessed by using all synthetic populations from previous years as a reference population. Moderate and high accuracies were achieved for the two traits, respectively, consistent with broad-sense heritability values. The present study represents the first demonstration and validation of genomic selection for seasonal biomass yield within a diverse commercial breeding program across multiple years. These results, supported by previous simulation studies, demonstrate the ability to predict sward-based phenotypic performance early in the process of individual plant selection, so shortening the breeding cycle, increasing the rate of genetic gain and allowing rapid adoption in ryegrass improvement programs.

Genotyping-by-sequencing through transcriptomics: implementation in a range of crop species with varying reproductive habits and ploidy levels.

The application of genomics in crops has the ability to significantly improve genetic gain for agriculture. Many marker-dense tools have been developed, but few have seen broad adoption in plant genomics due to issues of significant variations of genome size, levels of ploidy, single nucleotide polymorphism (SNP) frequency and reproductive habit. When combined with limited breeding activities, small research communities and scant sequence resources, the suitability of popular systems is often suboptimal and routinely fails to effectively balance cost-effectiveness and sample throughput. Genotyping-by-sequencing (GBS) encompasses a range of protocols including resequencing of the transcriptome. This study describes a skim GBS-transcriptomics (GBS-t) approach developed to be broadly applicable, cost-effective and high-throughput while still assaying a significant number of SNP loci. A range of crop species with differing levels of ploidy and degree of inbreeding/outbreeding were chosen, including perennial ryegrass, a diploid outbreeding forage grass; phalaris, a putative segmental allotetraploid outbreeding forage grass; lentil, a diploid inbreeding grain legume; and canola, an allotetraploid partially outbreeding oilseed. GBS-t was validated as a simple and largely automated, cost-effective method which generates sufficient SNPs (from 89 738 to 231 977) with acceptable levels of missing data and even genome coverage from c. 3 million sequence reads per sample. GBS-t is therefore a broadly applicable system suitable for many crops, offering advantages over other systems. The correct choice of subsequent sequence analysis software is important, and the bioinformatics process should be iterative and tailored to the specific challenges posed by ploidy variation and extent of heterozygosity.

Targeted genotyping-by-sequencing permits cost-effective identification and discrimination of pasture grass species and cultivars.

KEY MESSAGE: A targeted amplicon-based genotyping-by-sequencing approach has permitted cost-effective and accurate discrimination between ryegrass species (perennial, Italian and inter-species hybrid), and identification of cultivars based on bulked samples. Perennial ryegrass and Italian ryegrass are the most important temperate forage species for global agriculture, and are represented in the commercial pasture seed market by numerous cultivars each composed of multiple highly heterozygous individuals. Previous studies have identified difficulties in the use of morphophysiological criteria to discriminate between these two closely related taxa. Recently, a highly multiplexed single nucleotide polymorphism (SNP)-based genotyping assay has been developed that permits accurate differentiation between both species and cultivars of ryegrasses at the genetic level. This assay has since been further developed into an amplicon-based genotyping-by-sequencing (GBS) approach implemented on a second-generation sequencing platform, allowing accelerated throughput and ca. sixfold reduction in cost. Using the GBS approach, 63 cultivars of perennial, Italian and interspecific hybrid ryegrasses, as well as intergeneric Festulolium hybrids, were genotyped. The genetic relationships between cultivars were interpreted in terms of known breeding histories and indistinct species boundaries within the Lolium genus, as well as suitability of current cultivar registration methodologies. An example of applicability to quality assurance and control (QA/QC) of seed purity is also described. Rapid, low-cost genotypic assays provide new opportunities for breeders to more fully explore genetic diversity within breeding programs, allowing the combination of novel unique genetic backgrounds. Such tools also offer the potential to more accurately define cultivar identities, allowing protection of varieties in the commercial market and supporting processes of cultivar accreditation and quality assurance.

Identification of genetic factors influencing salt stress tolerance in white clover (Trifolium repens L.) by QTL analysis.

Allotetraploid (2n = 4x = 32) white clover (Trifolium repens L.) is the most commonly cultivated legume component of temperate pastures, sown in swards with a companion grass species. Genetic control of growth performance of white clover on saline land is highly important for dairy industries, due to increasing soil salinity problems. The objective of this study was to identify quantitative trait loci (QTLs) for salinity tolerance in terms of vegetative growth under stress. Two parental genetic maps consisting of 213 and 159 marker loci and spanning 1,973.0 and 1,837.6 cM, respectively, were constructed using simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers from a two-way pseudo-test cross F(1) population derived from pair-crossing of the Haifa(2) and LCL(2) genotypes. A total of 8 unique genomic regions on 8 linkage groups (LGs) of the Haifa(2) parental map and 6 unique regions on 5 LGs in the LCL(2) parental map were associated with plant growth under salt stress and relative growth under stress, as compared to control conditions. The results of this study indicate that salt tolerance in white clover is controlled by multiple QTLs, some at common locations, but each of limited magnitude. Location of these QTLs provides the genetic basis and potential for pyramiding of salt tolerance genes in breeding improvement.

Identification of homologous, homoeologous and paralogous sequence variants in an outbreeding allopolyploid species based on comparison with progenitor taxa.

The combination of homologous, homoeologous and paralogous classes of sequence variation presents major challenges for SNP discovery in outbreeding allopolyploid species. Previous in vitro gene-associated SNP discovery studies in the allotetraploid forage legume white clover (Trifolium repens L.) were vulnerable to such effects, leading to prohibitive levels of attrition during SNP validation. Identification of T. occidentale and T. pallescens as the putative diploid progenitors of white clover has permitted discrimination of the different sequence variant categories. Amplicons from selected abiotic stress tolerance-related genes were obtained using mapping family parents and individuals from each diploid species. Following cloning, progenitor comparison allowed tentative assignment of individual haplotypes to one or other sub-genome, as well as to gene copies within sub-genomes. A high degree of coincidence and identity between SNPs and HSVs was observed. Close similarity was observed between the genome of T. occidentale and one white clover sub-genome, but the affinity between T. pallescens and the other sub-genome was weaker, suggesting that a currently uncharacterised taxon may be the true second progenitor. Selected validated SNPs were attributed to individual sub-genomes by assignment to and naming of homoeologous linkage groups, providing the basis for improved genetic trait-dissection studies. The approach described in this study is broadly applicable to a range of allopolyploid taxa of equivocal ancestry.

SNP discovery, validation, haplotype structure and linkage disequilibrium in full-length herbage nutritive quality genes of perennial ryegrass (Lolium perenne L.).

Development of accurate high-throughput molecular marker systems such as SNPs permits evaluation and selection of favourable gene variants to accelerate elite varietal production. SNP discovery in perennial ryegrass has been based on PCR amplification and sequencing of multiple amplicons designed to scan all components of the transcriptional unit. Full-length genes (with complete intron-exon structure and promoter information) corresponding to well-defined biochemical functions such as lignin biosynthesis and oligosaccharide metabolism are ideal for complete SNP haplotype determination. Multiple SNPs at regular intervals across the transcriptional unit were detected within and between the heterozygous parents and validated in the progeny of the F (1)(NA(6) x AU(6)) genetic mapping family. Haplotype structures in the parental genotypes were defined and haplotypic abundance, structure and variation were assessed in diverse germplasm sources. Decay of LD to r (2) values of c. 0.2 typically occurs over 500-3,000 bp, comparable with gene length and with little apparent variation between diverse, ecotypic and varietal population sub-groups. Similar patterns were revealed as limited blocks of intragenic LD. The results are compatible with the reproductive biology of perennial ryegrass and the effects of large ancestral population size. This analysis provides crucial information to validate strategies for correlation of haplotypic diversity and phenotypic variation through association mapping.

Gene-associated single nucleotide polymorphism discovery in perennial ryegrass (Lolium perenne L.).

Molecular genetic marker development in perennial ryegrass has largely been dependent on anonymous sequence variation. The availability of a large-scale EST resource permits the development of functionally-associated genetic markers based on SNP variation in candidate genes. Genic SNP loci and associated haplotypes are suitable for implementation in molecular breeding of outbreeding forage species. Strategies for in vitro SNP discovery through amplicon cloning and sequencing have been designed and implemented. Putative SNPs were identified within and between the parents of the F(1)(NA(6) x AU(6)) genetic mapping family and were validated among progeny individuals. Proof-of-concept for the process was obtained using the drought tolerance-associated LpASRa2 gene. SNP haplotype structures were determined and correlated with predicted amino acid changes. Gene-length LD was evaluated across diverse germplasm collections. A survey of SNP variation across 100 candidate genes revealed a high frequency of SNP incidence (c. 1 per 54 bp), with similar proportions in exons and introns. A proportion (c. 50%) of the validated genic SNPs were assigned to the F(1)(NA(6) x AU(6)) genetic map, showing high levels of coincidence with previously mapped RFLP loci. The perennial ryegrass SNP resource will enable genetic map integration, detailed LD studies and selection of superior allele content during varietal development.