Integration of Infinium and Axiom SNP array data in the outcrossing species Malus × domestica and causes for seemingly incompatible calls

Development of a next generation SNP genotyping array for wheat

High-throughput genotyping arrays have provided a cost-effective, reliable and interoperable system for genotyping hexaploid wheat and its relatives. Existing, highly cited arrays including our 35K Wheat Breeder's array and the Illumina 90K array were designed based on a limited amount of varietal sequence diversity and with imperfect knowledge of SNP positions. Recent progress in wheat sequencing has given us access to a vast pool of SNP diversity, whilst technological improvements have allowed us to fit significantly more probes onto a 384-well format Axiom array than previously possible. Here we describe a novel Axiom genotyping array, the 'Triticum aestivum Next Generation' array (TaNG), largely derived from whole genome skim sequencing of 204 elite wheat lines and 111 wheat landraces taken from the Watkins 'Core Collection'. We used a novel haplotype optimization approach to select SNPs with the highest combined varietal discrimination and a design iteration step to test and replace SNPs which failed to convert to reliable markers. The final design with 43 372 SNPs contains a combination of haplotype-optimized novel SNPs and legacy cross-platform markers. We show that this design has an improved distribution of SNPs compared to previous arrays and can be used to generate genetic maps with a significantly higher number of distinct bins than our previous array. We also demonstrate the improved performance of TaNGv1.1 for Genome-wide association studies (GWAS) and its utility for Copy Number Variation (CNV) analysis. The array is commercially available with supporting marker annotations and initial genotyping results freely available.

Why we thrive beneath a northern sky - genomic signals of selection in apple for adaptation to northern Sweden

Good understanding of the genomic regions underlying adaptation of apple to boreal climates is needed to facilitate efficient breeding of locally adapted apple cultivars. Proper infrastructure for phenotyping and evaluation is essential for identification of traits responsible for adaptation, and dissection of their genetic composition. However, such infrastructure is costly and currently not available for the boreal zone of northern Sweden. Therefore, we used historical pomological data on climate adaptation of 59 apple cultivars and whole genome sequencing to identify genomic regions that have undergone historical selection among apple cultivars recommended for cultivation in northern Sweden. We found the apple collection to be composed of two ancestral groups that are largely concordant with the grouping into 'hardy' and 'not hardy' cultivars based on the pomological literature. Using a number of genome-wide scans for signals of selection, we obtained strong evidence of positive selection at a genomic region around 29 MbHFTH1 of chromosome 1 among apple cultivars in the 'hardy' group. Using phased genotypic data from the 20 K apple Infinium® SNP array, we identified haplotypes associated with the two cultivar groups and traced transmission of these haplotypes through the pedigrees of some apple cultivars. This demonstrates that historical data from pomological literature can be analyzed by population genomic approaches as a step towards revealing the genomic control of a key property for a horticultural niche market. Such knowledge is needed to facilitate efficient breeding strategies for development of locally adapted apple cultivars in the future. The current study illustrates the response to a very strong selective pressure imposed on tree crops by climatic factors, and the importance of genetic research on this topic and feasibility of breeding efforts in the light of the ongoing climate change.

Genomic prediction and genome-wide association study using combined genotypic data from different genotyping systems: application to apple fruit quality traits

With advances in next-generation sequencing technologies, various marker genotyping systems have been developed for genomics-based approaches such as genomic selection (GS) and genome-wide association study (GWAS). As new genotyping platforms are developed, data from different genotyping platforms must be combined. However, the potential use of combined data for GS and GWAS has not yet been clarified. In this study, the accuracy of genomic prediction (GP) and the detection power of GWAS increased for most fruit quality traits of apples when using combined data from different genotyping systems, Illumina Infinium single-nucleotide polymorphism array and genotyping by random amplicon sequencing-direct (GRAS-Di) systems. In addition, the GP model, which considered the inbreeding effect, further improved the accuracy of the seven fruit traits. Runs of homozygosity (ROH) islands overlapped with the significantly associated regions detected by the GWAS for several fruit traits. Breeders may have exploited these regions to select promising apples by breeders, increasing homozygosity. These results suggest that combining genotypic data from different genotyping platforms benefits the GS and GWAS of fruit quality traits in apples. Information on inbreeding could be beneficial for improving the accuracy of GS for fruit traits of apples; however, further analysis is required to elucidate the relationship between the fruit traits and inbreeding depression (e.g. decreased vigor).

Whole-Genome Sequencing Data Reveal New Loci Affecting Milk Production in German Black Pied Cattle (DSN)

German Black Pied (DSN) is considered an ancestral population of the Holstein breed. The goal of the current study was to fine-map genomic loci for milk production traits and to provide sequence variants for selection. We studied genome-wide associations for milk-production traits in 2160 DSN cows. Using 11.7 million variants from whole-genome sequencing of 304 representative DSN cattle, we identified 1980 associated variants (-log10(p) ≥ 7.1) in 13 genomic loci on 9 chromosomes. The highest significance was found for the MGST1 region affecting milk fat content (-log10(p) = 11.93, MAF = 0.23, substitution effect of the minor allele (ßMA) = -0.151%). Different from Holstein, DGAT1 was fixed (0.97) for the alanine protein variant for high milk and protein yield. A key gene affecting protein content was CSN1S1 (-log10(p) = 8.47, MAF = 049, ßMA = -0.055%) and the GNG2 region (-log10(p) = 10.48, MAF = 0.34, ßMA = 0.054%). Additionally, we suggest the importance of FGF12 for protein and fat yield, HTR3C for milk yield, TLE4 for milk and protein yield, and TNKS for milk and fat yield. Selection for favored alleles can improve milk yield and composition. With respect to maintaining the dual-purpose type of DSN, unfavored linkage to genes affecting muscularity has to be investigated carefully, before the milk-associated variants can be applied for selection in the small population.

DNA profiling with the 20K apple SNP array reveals Malus domestica hybridization and admixture in M. sieversii, M. orientalis, and M. sylvestris genebank accessions

The USDA-ARS National Plant Germplasm System (NPGS) apple collection in Geneva, NY, USA maintains accessions of the primary Malus domestica (Suckow) Borkh. progenitor species M. sieversii (Ledeb.) M. Roem., M. orientalis Uglitzk., and M. sylvestris (L.) Mill. Many of these accessions originated from seeds that were collected from wild populations in the species' centers of diversity. Some of these accessions have fruit phenotypes that suggest recent M. domestica hybridization, which if true would represent crop contamination of wild species populations and mislabeled species status of NPGS accessions. Pedigree connections and admixture between M. domestica and its progenitor species can be readily identified with apple SNP array data, despite such arrays not being designed for these purposes. To investigate species purity, most (463 accessions) of the NPGS accessions labeled as these three progenitor species were genotyped using the 20K apple SNP array. DNA profiles obtained were compared with a dataset of more than 5000 unique M. domestica apple cultivars. Only 212 accessions (151 M. sieversii, 26 M. orientalis, and 35 M. sylvestris) were identified as "pure" species representatives because their DNA profiles did not exhibit genotypic signatures of recent hybridization with M. domestica. Twenty-one accessions (17 M. sieversii, 1 M. orientalis, and 3 M. sylvestris) previously labeled as wild species were instead fully M. domestica. Previously unrealized hybridization and admixture between wild species and M. domestica was identified in 230 accessions (215 M. sieversii, 9 M. orientalis, and 6 M. sylvestris). Among these species-mislabeled accessions, 'Alexander', 'Gold Reinette', 'Charlamoff', 'Rosmarina Bianca', and 'King of the Pippins' were the most frequently detected M. domestica parents or grandparents. These results have implications for collection management, including germplasm distribution, and might affect conclusions of previous research focused on these three progenitor species in the NPGS apple collection. Specifically, accessions received from the NPGS for breeding and genomics, genetics, and evolutionary biology research might not be truly representative of their previously assigned species.

Multi-omics analyses reveal MdMYB10 hypermethylation being responsible for a bud sport of apple fruit color

Apple bud sports offer a rich resource for clonal selection of numerous elite cultivars. The accumulation of somatic mutations as plants develop may potentially impact the emergence of bud sports. Previous studies focused on somatic mutation in the essential genes associated with bud sports. However, the rate and function of genome-wide somatic mutations that accumulate when a bud sport arises remain unclear. In this study, we identified a branch from a 10-year-old tree of the apple cultivar 'Oregon Spur II' as a bud sport. The mutant branch showed reduced red coloration on fruit skin. Using this plant material, we assembled a high-quality haplotype reference genome consisting of 649.61 Mb sequences with a contig N50 value of 2.04 Mb. We then estimated the somatic mutation rate of the apple tree to be 4.56 × 10 -8 per base per year, and further identified 253 somatic single-nucleotide polymorphisms (SNPs), including five non-synonymous SNPs, between the original type and mutant samples. Transcriptome analyses showed that 69 differentially expressed genes between the original type and mutant fruit skin were highly correlated with anthocyanin content. DNA methylation in the promoter of five anthocyanin-associated genes was increased in the mutant compared with the original type as determined using DNA methylation profiling. Among the genetic and epigenetic factors that directly and indirectly influence anthocyanin content in the mutant apple fruit skin, the hypermethylated promoter of MdMYB10 is important. This study indicated that numerous somatic mutations accumulated at the emergence of a bud sport from a genome-wide perspective, some of which contribute to the low coloration of the bud sport.

A new method to reconstruct the direction of parent-offspring duo relationships using SNP array data and its demonstration on ancient and modern cultivars in the outcrossing species malus × domestica

Unordered parent-offspring (PO) relationships are an outstanding issue in pedigree reconstruction studies. Resolution of the order of these relationships would expand the results, conclusions, and usefulness of such studies; however, no such PO order resolution (POR) tests currently exist. This study describes two such tests, demonstrated using SNP array data in the outcrossing species apple (Malus × domestica) on a PO relationship of known order ("Keepsake" as a parent of "Honeycrisp") and two PO relationships previously ordered only via provenance information. The first test, POR-1, tests whether some of the extended haplotypes deduced from homozygous SNP calls from one individual in an unordered PO duo are composed of recombinant haplotypes from accurately phased SNP genotypes from the second individual. If so, the first individual would be the offspring of the second individual, otherwise the opposite relationship would be present. The second test, POR-2, does not require phased SNP genotypes and uses similar logic as the POR-1 test, albeit in a different approach. The POR-1 and POR-2 tests determined the correct relationship between "Keepsake" and "Honeycrisp". The POR-2 test confirmed "Reinette Franche" as a parent of "Nonpareil" and "Brabant Bellefleur" as a parent of "Court Pendu Plat". The latter finding conflicted with the recorded provenance information, demonstrating the need for these tests. The successful demonstration of these tests suggests they can add insights to future pedigree reconstruction studies, though caveats, like extreme inbreeding or selfing, would need to be considered where relevant.

Design and performance of a bovine 200 k SNP chip developed for endangered German Black Pied cattle (DSN)

BACKGROUND

German Black Pied cattle (DSN) are an endangered dual-purpose breed which was largely replaced by Holstein cattle due to their lower milk yield. DSN cattle are kept as a genetic reserve with a current herd size of around 2500 animals. The ability to track sequence variants specific to DSN could help to support the conservation of DSN's genetic diversity and to provide avenues for genetic improvement.

RESULTS

Whole-genome sequencing data of 304 DSN cattle were used to design a customized DSN200k SNP chip harboring 182,154 variants (173,569 SNPs and 8585 indels) based on ten selection categories. We included variants of interest to DSN such as DSN unique variants and variants from previous association studies in DSN, but also variants of general interest such as variants with predicted consequences of high, moderate, or low impact on the transcripts and SNPs from the Illumina BovineSNP50 BeadChip. Further, the selection of variants based on haplotype blocks ensured that the whole-genome was uniformly covered with an average variant distance of 14.4 kb on autosomes. Using 300 DSN and 162 animals from other cattle breeds including Holstein, endangered local cattle populations, and also a Bos indicus breed, performance of the SNP chip was evaluated. Altogether, 171,978 (94.31%) of the variants were successfully called in at least one of the analyzed breeds. In DSN, the number of successfully called variants was 166,563 (91.44%) while 156,684 (86.02%) were segregating at a minor allele frequency > 1%. The concordance rate between technical replicates was 99.83 ± 0.19%.

CONCLUSION

The DSN200k SNP chip was proved useful for DSN and other Bos taurus as well as one Bos indicus breed. It is suitable for genetic diversity management and marker-assisted selection of DSN animals. Moreover, variants that were segregating in other breeds can be used for the design of breed-specific customized SNP chips. This will be of great value in the application of conservation programs for endangered local populations in the future.

The use of shared haplotype length information for pedigree reconstruction in asexually propagated outbreeding crops, demonstrated for apple and sweet cherry

Pedigree information is of fundamental importance in breeding programs and related genetics efforts. However, many individuals have unknown pedigrees. While methods to identify and confirm direct parent-offspring relationships are routine, those for other types of close relationships have yet to be effectively and widely implemented with plants, due to complications such as asexual propagation and extensive inbreeding. The objective of this study was to develop and demonstrate methods that support complex pedigree reconstruction via the total length of identical by state haplotypes (referred to in this study as "summed potential lengths of shared haplotypes", SPLoSH). A custom Python script, HapShared, was developed to generate SPLoSH data in apple and sweet cherry. HapShared was used to establish empirical distributions of SPLoSH data for known relationships in these crops. These distributions were then used to estimate previously unknown relationships. Case studies in each crop demonstrated various pedigree reconstruction scenarios using SPLoSH data. For cherry, a full-sib relationship was deduced for 'Emperor Francis, and 'Schmidt', a half-sib relationship for 'Van' and 'Windsor', and the paternal grandparents of 'Stella' were confirmed. For apple, 29 cultivars were found to share an unknown parent, the pedigree of the unknown parent of 'Cox's Pomona' was reconstructed, and 'Fameuse' was deduced to be a likely grandparent of 'McIntosh'. Key genetic resources that enabled this empirical study were large genome-wide SNP array datasets, integrated genetic maps, and previously identified pedigree relationships. Crops with similar resources are also expected to benefit from using HapShared for empowering pedigree reconstruction.