How array design creates SNP ascertainment bias

Genetic diversity and selection signatures in sheep breeds

Natural and artificial selection in domesticated animals can cause specific changes in genomic regions known as selection signatures. Our study used the integrated haplotype score (iHS) and Tajima's D tests within non-overlapping windows of 100 kb to identify selection signatures, in addition to genetic diversity and linkage disequilibrium estimates in 9498 sheep from breeds in Ireland (Belclare, Charollais, Suffolk, Texel, and Vendeen). The mean observed and expected heterozygosity for all the sheep breeds were 0.353 and 0.355, respectively. Suffolk had the least genetic variation and, along with Texel, had slower linkage disequilibrium decay. iHS and Tajima's D detected selection signatures for all breeds, with some regions overlapping, thus forming longer segments of selection signatures. Common selection signatures were identified across iHS and Tajima's D methods for all breeds, with Belclare and Texel having several common regions under positive selection. Several genes were detected within the selection signature regions, including ITGA4, TLR3, and TGFB2 related to the immune system against endoparasites; DLG1, ROBO2, MXI1, MTMR2, CEP57, and FAM78B related to reproductive traits; WDR70 related to milk traits; SCHM1 and MYH15 related to meat traits; and TAS2R4, TAS2R39, and TAS2R40 related to adaptive traits. In conclusion, our results demonstrated moderate genetic diversity in the sheep breeds and detected and characterized selection signatures harboring genes associated with reproductive traits, milk production, meat production, and adaptive traits such as endoparasite resistance.

Genotyping Genebank Collections: Strategic Approaches and Considerations for Optimal Collection Management

The maintenance of plant germplasm and its genetic diversity is critical to preserving and making it available for food security, so this invaluable diversity is not permanently lost due to population growth and development, climate change, or changing needs from the growers and/or the marketplace. There are numerous genebanks worldwide that serve to preserve valuable plant germplasm for humankind's future and to serve as a resource for research, breeding, and training. The United States Department of Agriculture (USDA) National Plant Germplasm System (NPGS) and the Consultative Group for International Agricultural Research (CGIAR) both have a network of plant germplasm collections scattered across varying geographical locations preserving genetic resources for the future. Besides the USDA and CGIAR, there are germplasm collections established in many countries across the world that also aim to preserve crop and plant collections. Due to the advancement of technology, genotyping and sequencing whole genomes of plant germplasm collections is now feasible. Data from genotyping can help define genetic diversity within a collection, identify genetic gaps, reveal genetic redundancies and verify uniqueness, enable the comparison of collections of the same crop across genebanks (rationalization), and determine errors or mix-ups in genetic identity that may have occurred in a germplasm collection. Large-scale projects, such as genotyping germplasm collections, require strategic planning and the development of best practices. This article details strategies and best practices to consider when genotyping whole collections, considerations for the identity verification of germplasm and determining genetic replicates, quality management systems (QMS)/QC genotyping, and some use cases.

Development of a robust SNP marker set for genotyping diverse gene bank collections of polyploid roses

BACKGROUND

Due to genetic depletion in nature, gene banks play a critical role in the long-term conservation of plant genetic resources and the provision of a wide range of plant genetic diversity for research and breeding programs. Genetic information on accessions facilitates gene bank management and can help to conserve limited resources and to identify taxonomic misclassifications or mislabelling. Here, we developed SNP markers for genotyping 4,187 mostly polyploid rose accessions from large rose collections, including the German Genebank for Roses.

RESULTS

We filtered SNP marker information from the RhWag68k Axiom SNP array using call rates, uniformity of the four allelic dosage groups and chromosomal position to improve genotyping efficiency. After conversion to individual PACE® markers and further filtering, we selected markers with high discriminatory power. These markers were used to analyse 4,187 accessions with a mean call rate of 91.4%. By combining two evaluation methods, the mean call rate was increased to 95.2%. Additionally, the robustness against the genotypic groups used for calling was evaluated, resulting in a final set of 18 markers. Analyses of 94 pairs of assumed duplicate accessions included as controls revealed unexpected differences for eight pairs, which were confirmed using SSR markers. After removing the duplicates and filtering for accessions that were robustly called with all 18 markers, 141 out of the 1,957 accessions showed unexpected identical marker profiles with at least one other accession in our PACE® and SSR analysis. Given the attractiveness of NGS technologies, 13 SNPs from the marker set were also analysed using amplicon sequencing, with 76% agreement observed between PACE® and amplicon markers.

CONCLUSIONS

Although sampling error cannot be completely excluded, this is an indication that mislabelling occurs in rose collections and that molecular markers may be able to detect these cases. In future applications, our marker set could be used to develop a core reference set of representative accessions, and thus optimise the selection of gene bank accessions.

DArTseq genotyping facilitates identification of Aegilops biuncialis chromatin introgressed into bread wheat Mv9kr1

Wild wheat relative Aegilops biuncialis offers valuable traits for crop improvement through interspecific hybridization. However, gene transfer from Aegilops has been hampered by difficulties in detecting introgressed Ub- and Mb-genome chromatin in the wheat background at high resolution. The present study applied DArTseq technology to genotype two backcrossed populations (BC382, BC642) derived from crosses of wheat line Mv9kr1 with Ae. biuncialis accession, MvGB382 (early flowering and drought-tolerant) and MvGB642 (leaf rust-resistant). A total of 11,952 Aegilops-specific Silico-DArT markers and 8,998 wheat-specific markers were identified. Of these, 7,686 markers were assigned to Ub-genome chromosomes and 4,266 to Mb-genome chromosomes and were ordered using chromosome scale reference assemblies of hexaploid wheat and Ae. umbellulata. Ub-genome chromatin was detected in 5.7% of BC382 and 22.7% of BC642 lines, while 88.5% of BC382 and 84% of BC642 lines contained Mb-genome chromatin, predominantly the chromosomes 4Mb and 5Mb. The presence of alien chromatin was confirmed by microscopic analysis of mitotic metaphase cells using GISH and FISH, which allowed precise determination of the size and position of the introgression events. New Mv9kr1-Ae. biuncialis MvGB382 4Mb and 5Mb disomic addition lines together with a 5DS.5DL-5MbL recombination were identified. A possible effect of the 5MbL distal region on seed length has also been observed. Moreover, previously developed Mv9kr1-MvGB642 introgression lines were more precisely characterized. The newly developed cytogenetic stocks represent valuable genetic resources for wheat improvement, highlighting the importance of utilizing diverse genetic materials to enhance wheat breeding strategies.

Genetic characterization of cashmere goat (Capra hircus) populations in Mongolia

Objective

Characterization studies of the phenotypic and genetic diversity of Mongolian goats are limited, despite several goat breeds being registered in the country. This study aimed to evaluate the phenotypic and genetic diversity of 14 cashmere goat populations in Mongolia, consisting largely of identified goat breeds.

Methods

Body weight, cashmere quality, and coat color were the phenotypic traits considered in this study. A linear model was used to fit body weight and cashmere traits, and least squares means (LSMs) were estimated for the region and location classes. Genetic diversity and structure were assessed using a goat 50K SNP array.

Results

The studied populations exhibited greater phenotypic diversity at the regional level. A very small overall differentiation index (Fst: 0.017) was revealed by Wright's Fst and a very small overall inbreeding index (F ROH1 :0.019) was revealed based on runs of homozygosity. Genetic clustering of populations by principal components showed large variances for the two goat populations of the Russian admixture (Gobi Gurvan Saikhan and Uuliin Bor), and smaller but differentiated clusters for the remaining populations. Similar results were observed in the admixture analysis, which identified populations with the highest (Govi Gurvan Saikhan and Uuliin Bor) and lowest (Tsagaan Ovoo Khar) exotic admixtures. A genomewide association study (GWAS) of body weight and cashmere traits identified a few significant variants on chromosomes 2, 4, 5, 9, and 15, with the strongest variant for cashmere yield on chromosome 4. The GWAS on coat color yielded nine significant variants, with the strongest variants located on chromosomes 6, 13, and 18 and potential associations with KIT, ASIP, and MC1R genes. These signals were also found in other studies on coat color and patterns in goats.

Conclusion

Mongolian cashmere goats showed relatively low genetic differentiation and low inbreeding levels, possibly caused by the traditional pastoral livestock management system and the practice of trading breeding bucks across provinces, along with a recent increase in the goat population. Further investigation of cashmere traits using larger samples and alternative methods may help identify the genes or genomic regions underlying cashmere quality in goats.

Development and validation of a 66K SNP array for the hard clam (Mercenaria mercenaria)

BACKGROUND

The hard clam (Mercenaria mercenaria), a marine bivalve distributed along the U.S. eastern seaboard, supports a significant shellfish industry. Overharvest in the 1970s and 1980s led to a reduction in landings. While the transition of industry from wild harvest to aquaculture since that time has enhanced production, it has also exacerbated challenges such as disease outbreaks. In this study, we developed and validated a 66K SNP array designed to advance genetic studies and improve breeding programs in the hard clam, focusing particularly on the development of markers that could be useful in understanding disease resistance and environmental adaptability.

RESULTS

Whole-genome resequencing of 84 individual clam samples and 277 pooled clam libraries yielded over 305 million SNPs, which were filtered down to a set of 370,456 SNPs that were used as input for the design of a 66K SNP array. This medium-density array features 66,543 probes targeting coding and non-coding regions, including 70 mitochondrial SNPs, to capture the extensive genetic diversity within the species. The SNPs were distributed evenly throughout the clam genome, with an average interval of 25,641 bp between SNPs. The array incorporates markers for detecting the clam pathogen Mucochytrium quahogii (formerly QPX), enhancing its utility in disease management. Performance evaluation on 1,904 samples demonstrated a 72.7% pass rate with stringent quality control. Concordance testing affirmed the array's repeatability, with an average agreement of allele calls of 99.64% across multiple tissue types, highlighting its reliability. The tissue-specific analysis demonstrated that some tissue types yield better genotyping results than others. Importantly, the array, including its embedded mitochondrial markers, effectively elucidated complex genetic relationships across different clam groups, both wild populations and aquacultured stocks, showcasing its utility for detailed population genetics studies.

CONCLUSIONS

The 66K SNP array is a powerful and robust genotyping tool that offers unprecedented insights into the species' genomic architecture and population dynamics and that can greatly facilitate hard clam selective breeding. It represents an important resource that has the potential to transform clam aquaculture, thereby promoting industry sustainability and ecological and economic resilience.

Accurate genotype imputation from low-coverage whole-genome sequencing data of rainbow trout

With the rapid and significant cost reduction of next-generation sequencing, low-coverage whole-genome sequencing (lcWGS), followed by genotype imputation, is becoming a cost-effective alternative to single-nucleotide polymorphism (SNP)-array genotyping. The objectives of this study were 2-fold: (1) construct a haplotype reference panel for genotype imputation from lcWGS data in rainbow trout (Oncorhynchus mykiss); and (2) evaluate the concordance between imputed genotypes and SNP-array genotypes in 2 breeding populations. Medium-coverage (12×) whole-genome sequences were obtained from a total of 410 fish representing 5 breeding populations with various spawning dates. The short-read sequences were mapped to the rainbow trout reference genome, and genetic variants were identified using GATK. After data filtering, 20,434,612 biallelic SNPs were retained. The reference panel was phased with SHAPEIT5 and was used as a reference to impute genotypes from lcWGS data employing GLIMPSE2. A total of 90 fish from the Troutlodge November breeding population were sequenced with an average coverage of 1.3×, and these fish were also genotyped with the Axiom 57K rainbow trout SNP array. The concordance between array-based genotypes and imputed genotypes was 99.1%. After downsampling the coverage to 0.5×, 0.2×, and 0.1×, the concordance between array-based genotypes and imputed genotypes was 98.7, 97.8, and 96.7%, respectively. In the USDA odd-year breeding population, the concordance between array-based genotypes and imputed genotypes was 97.8% for 109 fish downsampled to 0.5× coverage. Therefore, the reference haplotype panel reported in this study can be used to accurately impute genotypes from lcWGS data in rainbow trout breeding populations.

Variant calling and genotyping accuracy of ddRAD-seq: Comparison with 20X WGS in layers

Whole Genome Sequencing (WGS) remains a costly or unsuitable method for routine genotyping of laying hens. Until now, breeding companies have been using or developing SNP chips. Nevertheless, alternatives methods based on sequencing have been developed. Among these, reduced representation sequencing approaches can offer sequencing quality and cost-effectiveness by reducing the genomic regions covered by sequencing. The aim of this study was to evaluate the ability of double digested Restriction site Associated DNA sequencing (ddRAD-seq) to identify and genotype SNPs in laying hens, by comparison with a presumed reliable WGS approach. Firstly, the sensitivity and precision of variant calling and the genotyping reliability of ddRADseq were determined. Next, the SNP Call Rate (CRSNP) and mean depth of sequencing per SNP (DPSNP) were compared between both methods. Finally, the effect of multiple combinations of thresholds for these parameters on genotyping reliability and amount of remaining SNPs in ddRAD-seq was studied. In raw form, the ddRAD-seq identified 349,497 SNPs evenly distributed on the genome with a CRSNP of 0.55, a DPSNP of 11X and a mean genotyping reliability rate per SNP of 80%. Considering genomic regions covered by expected enzymatic fragments (EFs), the sensitivity of the ddRAD-seq was estimated at 32.4% and its precision at 96.4%. The low CRSNP and DPSNP values were explained by the detection of SNPs outside the EFs theoretically generated by the ddRAD-seq protocol. Indeed, SNPs outside the EFs had significantly lower CRSNP (0.25) and DPSNP (1X) values than SNPs within the EFs (0.7 and 17X, resp.). The study demonstrated the relationship between CRSNP, DPSNP, genotyping reliability and the number of SNPs retained, to provide a decision-support tool for defining filtration thresholds. Severe quality control over ddRAD-seq data allowed to retain a minimum of 40% of the SNPs with a CcR of 98%. Then, ddRAD-seq was defined as a suitable method for variant calling and genotyping in layers.

Promiscuous potato: elucidating genetic identity and the complex genetic relationships of a cultivated potato germplasm collection

A total of 3,860 accessions from the global in trust clonal potato germplasm collection w3ere genotyped with the Illumina Infinium SolCAP V2 12K potato SNP array to evaluate genetic diversity and population structure within the potato germplasm collection. Diploid, triploid, tetraploid, and pentaploid accessions were included representing the cultivated potato taxa. Heterozygosity ranged from 9.7% to 66.6% increasing with ploidy level with an average heterozygosity of 33.5%. Identity, relatedness, and ancestry were evaluated using hierarchal clustering and model-based Bayesian admixture analyses. Errors in genetic identity were revealed in a side-by-side comparison of in vitro clonal material with the original mother plants revealing mistakes putatively occurring during decades of processing and handling. A phylogeny was constructed to evaluate inter- and intraspecific relationships which together with a STRUCTURE analysis supported both commonly used treatments of potato taxonomy. Accessions generally clustered based on taxonomic and ploidy classifications with some exceptions but did not consistently cluster by geographic origin. STRUCTURE analysis identified putative hybrids and suggested six genetic clusters in the cultivated potato collection with extensive gene flow occurring among the potato populations, implying most populations readily shared alleles and that introgression is common in potato. Solanum tuberosum subsp. andigena (ADG) and S. curtilobum (CUR) displayed significant admixture. ADG likely has extensive admixture due to its broad geographic distribution. Solanum phureja (PHU), Solanum chaucha (CHA)/Solanum stenotomum subsp. stenotomum (STN), and Solanum tuberosum subsp. tuberosum (TBR) populations had less admixture from an accession/population perspective relative to the species evaluated. A core and mini core subset from the genebank material was also constructed. SNP genotyping was also carried out on 745 accessions from the Seed Savers potato collection which confirmed no genetic duplication between the two potato collections, suggesting that the collections hold very different genetic resources of potato. The Infinium SNP Potato Array is a powerful tool that can provide diversity assessments, fingerprint genebank accessions for quality management programs, use in research and breeding, and provide insights into the complex genetic structure and hybrid origin of the diversity present in potato genetic resource collections.

GWAS advancements to investigate disease associations and biological mechanisms

Genome-wide association studies (GWAS) have been instrumental in elucidating the genetic architecture of various traits and diseases. Despite the success of GWAS, inherent limitations such as identifying rare and ultra-rare variants, the potential for spurious associations, and in pinpointing causative agents can undermine diagnostic capabilities. This review provides an overview of GWAS and highlights recent advances in genetics that employ a range of methodologies, including Whole Genome Sequencing (WGS), Mendelian Randomization (MR), the Pangenome's high-quality T2T-CHM13 panel, and the Human BioMolecular Atlas Program (HuBMAP), as potential enablers of current and future GWAS research. State of the literature demonstrate the capabilities of these techniques in enhancing the statistical power of GWAS. WGS, with its comprehensive approach, captures the entire genome, surpassing the capabilities of the traditional GWAS technique focused on predefined Single Nucleotide Polymorphism (SNP) sites. The Pangenome's T2T-CHM13 panel, with its holistic approach, aids in the analysis of regions with high sequence identity, such as segmental duplications (SDs). Mendelian Randomization has advanced causative inference, improving clinical diagnostics and facilitating definitive conclusions. Furthermore, spatial biology techniques like HuBMAP, enable 3D molecular mapping of tissues at single-cell resolution, offering insights into pathology of complex traits. This study aims to elucidate and advocate for the increased application of these technologies, highlighting their potential to shape the future of GWAS research.

Anthropogenic events and responses to environmental stress are shaping the genomes of Ethiopian indigenous goats

Anthropological and biophysical processes have shaped livestock genomes over Millenia and can explain their current geographic distribution and genetic divergence. We analyzed 57 Ethiopian indigenous domestic goat genomes alongside 67 equivalents of east, west, and north-west African, European, South Asian, Middle East, and wild Bezoar goats. Cluster, ADMIXTURE (K = 4) and phylogenetic analysis revealed four genetic groups comprising African, European, South Asian, and wild Bezoar goats. The Middle Eastern goats had an admixed genome of these four genetic groups. At K = 5, the West African Dwarf and Moroccan goats were separated from East African goats demonstrating a likely historical legacy of goat arrival and dispersal into Africa via the coastal Mediterranean Sea and the Horn of Africa. FST, XP-EHH, and Hp analysis revealed signatures of selection in Ethiopian goats overlaying genes for thermo-sensitivity, oxidative stress response, high-altitude hypoxic adaptation, reproductive fitness, pathogen defence, immunity, pigmentation, DNA repair, modulation of renal function and integrated fluid and electrolyte homeostasis. Notable examples include TRPV1 (a nociception gene); PTPMT1 (a critical hypoxia survival gene); RETREG (a regulator of reticulophagy during starvation), and WNK4 (a molecular switch for osmoregulation). These results suggest that human-mediated translocations and adaptation to contrasting environments are shaping indigenous African goat genomes.

Unveiling the predominance of Saccharum spontaneum alleles for resistance to orange rust in sugarcane using genome-wide association

KEY MESSAGE

Six QTLs of resistance to sugarcane orange rust were identified in modern interspecific hybrids by GWAS. For five of them, the resistance alleles originated from S. spontaneum. Altogether, they efficiently predict disease resistance. Sugarcane orange rust (SOR) is a threatening emerging disease in many sugarcane industries worldwide. Improving the genetic resistance of commercial cultivars remains the most promising solution to control this disease. In this study, an association panel of 568 modern interspecific sugarcane hybrids (Saccharum officinarum x S. spontaneum) from Réunion's breeding program was evaluated for its resistance to SOR under natural conditions of infection. Two genome-wide association studies (GWAS) were conducted between disease reactions and 183,842 single nucleotide polymorphism (SNP) markers obtained by targeted genotyping-by-sequencing. Five resistance quantitative trait loci (QTLs), named Oru1, Oru2, Oru3, Oru4 and Oru5, were identified using a single-locus GWAS (SL-GWAS). These five QTLs all originated from the species S. spontaneum. A multi-locus GWAS (ML-GWAS) uncovered an additional but less significant resistance QTL named Oru6, which originated from S. officinarum. All six QTLs had a moderate to major phenotypic effect on disease resistance. Prediction accuracy estimated with linear regression models based on each of the five QTLs identified by SL-GWAS was between 0.16-0.41. Altogether, these five QTLs provided a relatively high prediction accuracy of 0.60. In comparison, accuracies obtained with six genome-wide prediction models (i.e., GBLUP, Bayes-A, Bayes-B, Bayes-C, Bayesian Lasso and RKHS) reached only 0.65. The good prediction accuracy of disease resistance provided by the QTLs and the predominant S. spontaneum origin of their resistance alleles pave the way for effective marker-assisted breeding strategies.

Genome-wide association study identifies genetic variants underlying footrot in Portuguese Merino sheep

BACKGROUND

Ovine footrot caused by Dichelobacter nodosus (D. nodosus) is a contagious disease with serious economic and welfare impacts in sheep production systems worldwide. A better understanding of the host genetic architecture regarding footrot resistance/susceptibility is crucial to develop disease control strategies that efficiently reduce infection and its severity. A genome-wide association study was performed using a customized SNP array (47,779 SNPs in total) to identify genetic variants associated to footrot resistance/susceptibility in two Portuguese native breeds, i.e. Merino Branco and Merino Preto, and a population of crossbred animals. A cohort of 1375 sheep sampled across 17 flocks, located in the Alentejo region (southern Portugal), was included in the analyses.

RESULTS

Phenotypes were scored from 0 (healthy) to 5 (severe footrot) based on visual inspection of feet lesions, following the Modified Egerton System. Using a linear mixed model approach, three SNPs located on chromosome 24 reached genome-wide significance after a Bonferroni correction (p < 0.05). Additionally, six genome-wide suggestive SNPs were identified each on chromosomes 2, 4, 7, 8, 9 and 15. The annotation and KEGG pathway analyses showed that these SNPs are located within regions of candidate genes such as the nonsense mediated mRNA decay associated PI3K related kinase (SMG1) (chromosome 24) and the RALY RNA binding protein like (RALYL) (chromosome 9), both involved in immunity, and the heparan sulfate proteoglycan 2 (HSPG2) (chromosome 2) and the Thrombospodin 1 (THBS1) (chromosome 7) implicated in tissue repair and wound healing processes.

CONCLUSION

This is the first attempt to identify molecular markers associated with footrot in Portuguese Merino sheep. These findings provide relevant information on a likely genetic association underlying footrot resistance/susceptibility and the potential candidate genes affecting this trait. Genetic selection strategies assisted on the information obtained from this study could enhance Merino sheep-breeding programs, in combination with farm management strategies, for a more effective and sustainable long-term solution for footrot control.

Genome-wide genetic variation and population structure of native and cosmopolitan cattle breeds reared in Türkiye

This is the first comprehensive study to reveal genetic variation and population structure at genome level in six Anatolian (Anatolian Black, East Anatolian Red, South Anatolian Red, South Anatolian Yellow, Turkish Grey Steppe, and Zavot) and two cosmopolitan (Brown Swiss and Holstein Friesian) cattle breeds reared in Türkiye. Being 20 samples from each population, a total of 160 blood samples retrieved from representative herds were utilized to generate genomic libraries by ddRADseq method. Genomic libraries sequenced by Illumina HiSeq X Ten instrument revealed a total of 211,119 bi-allelic SNPs with high call rate. Compared to cosmopolitan cattle breeds, a higher genetic variation was observed in native Turkish cattle with an average of 0.380 observed heterozygosity. Genetic distances were comparatively low between native cattle breeds, whereas the highest genetic distance (0.064) was detected between South Anatolian Yellow and Brown Swiss. Population structure analyses showed that the native Turkish and cosmopolitan cattle breeds were clearly different from each other according to their phylogenetic origin. Besides, a high level of genetic admixture was detected among Anatolian Black, Turkish Grey Steppe, South Anatolian Red, and South Anatolian Yellow, whereas East Anatolian Red and Zavot were distinct from the other native and cosmopolitan cattle breeds. TreeMix algorithm under the assumption of one and two migration events revealed a migration route from Anatolian clade to Anatolian Black, while a second migration edge was drawn from Brown Swiss to East Anatolian Red. This study demonstrates the importance of national conservation studies in the native breeds whose population size has dramatically decreased. In addition, SNP arrays and next-generation sequencing platforms are recommended for future studies to reveal the genetic variation of other local Turkish livestock species to arrange effective conservation programs.

Landmarks in the history of selective sweeps

Half a century ago, a seminal article on the hitchhiking effect by Smith and Haigh inaugurated the concept of the selection signature. Selective sweeps are characterised by the rapid spread of an advantageous genetic variant through a population and hence play an important role in shaping evolution and research on genetic diversity. The process by which a beneficial allele arises and becomes fixed in a population, leading to a increase in the frequency of other linked alleles, is known as genetic hitchhiking or genetic draft. Kimura's neutral theory and hitchhiking theory are complementary, with Kimura's neutral evolution as the 'null model' and positive selection as the 'signal'. Both are widely accepted in evolution, especially with genomics enabling precise measurements. Significant advances in genomic technologies, such as next-generation sequencing, high-density SNP arrays and powerful bioinformatics tools, have made it possible to systematically investigate selection signatures in a variety of species. Although the history of selection signatures is relatively recent, progress has been made in the last two decades, owing to the increasing availability of large-scale genomic data and the development of computational methods. In this review, we embark on a journey through the history of research on selective sweeps, ranging from early theoretical work to recent empirical studies that utilise genomic data.

Modeling of African population history using f-statistics is biased when applying all previously proposed SNP ascertainment schemes

f-statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. Not only are they guaranteed to allow robust tests of the fits of proposed models of population history to data when analyzing full genome sequencing data-that is, all single nucleotide polymorphisms (SNPs) in the individuals being analyzed-but they are also guaranteed to allow robust tests of models for SNPs ascertained as polymorphic in a population that is an outgroup in a phylogenetic sense to all groups being analyzed. True "outgroup ascertainment" is in practice impossible in humans because our species has arisen from a substructured ancestral population that does not descend from a homogeneous ancestral population going back many hundreds of thousands of years into the past. However, initial studies suggested that non-outgroup-ascertainment schemes might produce robust enough results using f-statistics, and that motivated widespread fitting of models to data using non-outgroup-ascertained SNP panels such as the "Affymetrix Human Origins array" which has been genotyped on thousands of modern individuals from hundreds of populations, or the "1240k" in-solution enrichment reagent which has been the source of about 70% of published genome-wide data for ancient humans. In this study, we show that while analyses of population history using such panels work well for studies of relationships among non-African populations and one African outgroup, when co-modeling more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans), fitting of f-statistics to such SNP sets is expected to frequently lead to false rejection of true demographic histories, and failure to reject incorrect models. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, has limited statistical power and retains important biases. However, by carrying out simulations of diverse demographic histories, we show that bias in inferences based on f-statistics can be minimized by ascertaining on variants common in a union of diverse African groups; such ascertainment retains high statistical power while allowing co-analysis of archaic and modern groups.

From chip to SNP: Rapid development and evaluation of a targeted capture genotyping-by-sequencing approach to support research and management of a plaguing rodent

The management of invasive species has been greatly enhanced by population genetic analyses of multilocus single-nucleotide polymorphism (SNP) datasets that provide critical information regarding pest population structure, invasion pathways, and reproductive biology. For many applications there is a need for protocols that offer rapid, robust and efficient genotyping on the order of hundreds to thousands of SNPs, that can be tailored to specific study populations and that are scalable for long-term monitoring schemes. Despite its status as a model laboratory species, there are few existing resources for studying wild populations of house mice (Mus musculus spp.) that strike this balance between data density and laboratory efficiency. Here we evaluate the utility of a custom targeted capture genotyping-by-sequencing approach to support research on plaguing house mouse populations in Australia. This approach utilizes 3,651 hybridization capture probes targeting genome-wide SNPs identified from a sample of mice collected in grain-producing regions of southeastern Australia genotyped using a commercially available microarray platform. To assess performance of the custom panel, we genotyped wild caught mice (N = 320) from two adjoining farms and demonstrate the ability to correctly assign individuals to source populations with high confidence (mean >95%), as well as robust kinship inference within sites. We discuss these results in the context of proposed applications for future genetic monitoring of house mice in Australia.

A Multibreed Genome-Wide Association Study for Cattle Leukocyte Telomere Length

Telomeres are terminal DNA regions of chromosomes that prevent chromosomal fusion and degradation during cell division. In cattle, leukocyte telomere length (LTL) is associated with longevity, productive lifespan, and disease susceptibility. However, the genetic basis of LTL in this species is less studied than in humans. In this study, we utilized the whole-genome resequencing data of 239 animals from 17 cattle breeds for computational leukocyte telomere length estimation and subsequent genome-wide association study of LTL. As a result, we identified 42 significant SNPs, of which eight were found in seven genes (EXOC6B, PTPRD, RPS6KC1, NSL1, AGBL1, ENSBTAG00000052188, and GPC1) when using covariates for two major breed groups (Turano-Mongolian and European). Association analysis with covariates for breed effect detected 63 SNPs, including 13 in five genes (EXOC6B, PTPRD, RPS6KC1, ENSBTAG00000040318, and NELL1). The PTPRD gene, demonstrating the top signal in analysis with breed effect, was previously associated with leukocyte telomere length in cattle and likely is involved in the mechanism of alternative lengthening of telomeres. The single nucleotide variants found could be tested for marker-assisted selection to improve telomere-length-associated traits.

Using species richness calculations to model the global profile of unsampled pathogenic variants: Examples from BRCA1 and BRCA2

There have been many surveys of genetic variation in BRCA1 and BRCA2 to identify variant prevalence and catalogue population specific variants, yet none have evaluated the magnitude of unobserved variation. We applied species richness estimation methods from ecology to estimate "variant richness" and determine how many germline pathogenic BRCA1/2 variants have yet to be identified and the frequency of these missing variants in different populations. We also estimated the prevalence of germline pathogenic BRCA1/2 variants and identified those expected to be most common. Data was obtained from a literature search including studies conducted globally that tested the entirety of BRCA1/2 for pathogenic variation. Across countries, 45% to 88% of variants were estimated to be missing, i.e., present in the population but not observed in study data. Estimated variant frequencies in each country showed a higher proportion of rare variants compared to recurrent variants. The median prevalence estimate of BRCA1/2 pathogenic variant carriers was 0.64%. BRCA1 c.68_69del is likely the most recurrent BRCA1/2 variant globally due to its estimated prevalence in India. Modeling variant richness using ecology methods may assist in evaluating clinical targeted assays by providing a picture of what is observed with estimates of what is still unknown.

Modeling of African population history using f -statistics can be highly biased and is not addressed by previously suggested SNP ascertainment schemes

f -statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. These statistics can provide strong evidence for either admixture or cladality, which can be robust to substantial rates of errors or missing data. f -statistics are guaranteed to be unbiased under "SNP ascertainment" (analyzing non-randomly chosen subsets of single nucleotide polymorphisms) only if it relies on a population that is an outgroup for all groups analyzed. However, ascertainment on a true outgroup that is not co-analyzed with other populations is often impractical and uncommon in the literature. In this study focused on practical rather than theoretical aspects of SNP ascertainment, we show that many non-outgroup ascertainment schemes lead to false rejection of true demographic histories, as well as to failure to reject incorrect models. But the bias introduced by common ascertainments such as the 1240K panel is mostly limited to situations when more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans) or non-human outgroups are co-modelled, for example, f 4 -statistics involving one non-African group, two African groups, and one archaic group. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, cannot fix all these problems since for some classes of f -statistics it is not a clean outgroup ascertainment, and in other cases it demonstrates relatively low power to reject incorrect demographic models since it provides a relatively small number of variants common in anatomically modern humans. And due to the paucity of high-coverage archaic genomes, archaic individuals used for ascertainment often act as sole representatives of the respective groups in an analysis, and we show that this approach is highly problematic. By carrying out large numbers of simulations of diverse demographic histories, we find that bias in inferences based on f -statistics introduced by non-outgroup ascertainment can be minimized if the derived allele frequency spectrum in the population used for ascertainment approaches the spectrum that existed at the root of all groups being co-analyzed. Ascertaining on sites with variants common in a diverse group of African individuals provides a good approximation to such a set of SNPs, addressing the great majority of biases and also retaining high statistical power for studying population history. Such a "pan-African" ascertainment, although not completely problem-free, allows unbiased exploration of demographic models for the widest set of archaic and modern human populations, as compared to the other ascertainment schemes we explored.

Identification of major quantitative trait loci and candidate genes for seed weight in soybean

Four major quantitative trait loci for 100-seed weight were identified in a soybean RIL population under five environments, and the most likely candidate genes underlying these loci were identified. Seed weight is an important target of soybean breeding. However, the genes underlying the major quantitative trait loci (QTL) controlling seed weight remain largely unknown. In this study, a soybean population of 300 recombinant inbred lines (RILs) derived from a cross between PI595843 (PI) and WH was used to map the QTL and identify candidate genes for seed weight. The RIL population was genotyped through whole genome resequencing, and phenotyped for 100-seed weight under five environments. A total of 38 QTL were detected, and four major QTL, each explained at least 10% of the variation in 100-seed weight, were identified. Six candidate genes within these four major QTL regions were identified by analyses of their tissue expression patterns, gene annotations, and differential gene expression levels in soybean seeds during four developmental stages between two parental lines. Further sequence variation analyses revealed a C to T substitution in the first exon of the Glyma.19G143300, resulting in an amino acid change between PI and WH, and thus leading to a different predicted kinase domain, which might affect its protein function. Glyma.19G143300 is highly expressed in soybean seeds and encodes a leucine-rich repeat receptor-like protein kinase (LRR-RLK). Its predicted protein has typical domains of LRR-RLK family, and phylogenetic analyses reveled its similarity with the known LRR-RLK protein XIAO (LOC_Os04g48760), which is involved in controlling seed size. The major QTL and candidate genes identified in this study provide useful information for molecular breeding of new soybean cultivars with desirable seed weight.

Testing the generalizability of ancestry-specific polygenic risk scores to predict prostate cancer in sub-Saharan Africa

BACKGROUND

Genome-wide association studies do not always replicate well across populations, limiting the generalizability of polygenic risk scores (PRS). Despite higher incidence and mortality rates of prostate cancer in men of African descent, much of what is known about cancer genetics comes from populations of European descent. To understand how well genetic predictions perform in different populations, we evaluated test characteristics of PRS from three previous studies using data from the UK Biobank and a novel dataset of 1298 prostate cancer cases and 1333 controls from Ghana, Nigeria, Senegal, and South Africa.

RESULTS

Allele frequency differences cause predicted risks of prostate cancer to vary across populations. However, natural selection is not the primary driver of these differences. Comparing continental datasets, we find that polygenic predictions of case vs. control status are more effective for European individuals (AUC 0.608-0.707, OR 2.37-5.71) than for African individuals (AUC 0.502-0.585, OR 0.95-2.01). Furthermore, PRS that leverage information from African Americans yield modest AUC and odds ratio improvements for sub-Saharan African individuals. These improvements were larger for West Africans than for South Africans. Finally, we find that existing PRS are largely unable to predict whether African individuals develop aggressive forms of prostate cancer, as specified by higher tumor stages or Gleason scores.

CONCLUSIONS

Genetic predictions of prostate cancer perform poorly if the study sample does not match the ancestry of the original GWAS. PRS built from European GWAS may be inadequate for application in non-European populations and perpetuate existing health disparities.

Integrated Approach in Genomic Selection to Accelerate Genetic Gain in Sugarcane

Marker-assisted selection (MAS) has been widely used in the last few decades in plant breeding programs for the mapping and introgression of genes for economically important traits, which has enabled the development of a number of superior cultivars in different crops. In sugarcane, which is the most important source for sugar and bioethanol, marker development work was initiated long ago; however, marker-assisted breeding in sugarcane has been lagging, mainly due to its large complex genome, high levels of polyploidy and heterozygosity, varied number of chromosomes, and use of low/medium-density markers. Genomic selection (GS) is a proven technology in animal breeding and has recently been incorporated in plant breeding programs. GS is a potential tool for the rapid selection of superior genotypes and accelerating breeding cycle. However, its full potential could be realized by an integrated approach combining high-throughput phenotyping, genotyping, machine learning, and speed breeding with genomic selection. For better understanding of GS integration, we comprehensively discuss the concept of genetic gain through the breeder’s equation, GS methodology, prediction models, current status of GS in sugarcane, challenges of prediction accuracy, challenges of GS in sugarcane, integrated GS, high-throughput phenotyping (HTP), high-throughput genotyping (HTG), machine learning, and speed breeding followed by its prospective applications in sugarcane improvement.

Genomic analysis of the rare British Lop pig and identification of distinctive genomic markers

Concentration of production on a few commercial pig breeds has led to the marginalization of many native, numerically small breeds, increasing their risk of endangerment. In the UK, one such rare breed is the British Lop, a lop-eared breed, of similar origin to the Welsh breed. The objective of the present study was to address the genomic status of the British Lop and its relationship with other breeds and identify a small set of genomic markers that uniquely characterize and distinguish British Lop animals. Results have shown that the British Lop is a relatively distinct population with reduced genomic diversity and effective size consistent with its status as a rare breed. Furthermore, we demonstrated the genetic closeness of the British Lop to phenotypically similar breeds such as Landrace and Welsh as well Large White, Middle White and Pietrain. Finally, a set of 75 Single Nucleotide Polymorphisms distributed across multiple chromosomes were identified and validated as markers that can consistently distinguish British Lops from other closely related breeds. Results may inform breeding and management strategies aiming to enhance diversity as well as the development of a breed purity test.

Spatiotemporal monitoring of the rare northern dragonhead (Dracocephalum ruyschiana, Lamiaceae) - SNP genotyping and environmental niche modeling herbarium specimens

The species we have studied the spatiotemporal genetic change in the northern dragonhead, a plant species that has experienced a drastic population decline and habitat loss in Europe. We have added a temporal perspective to the monitoring of northern dragonhead in Norway by genotyping herbarium specimens up to 200 years old. We have also assessed whether northern dragonhead has achieved its potential distribution in Norway. To obtain the genotype data from 130 herbarium specimens collected from 1820 to 2008, mainly from Norway (83) but also beyond (47), we applied a microfluidic array consisting of 96 SNP markers. To assess temporal genetic change, we compared our new genotype data with existing data from modern samples. We used sample metadata and observational records to model the species' environmental niche and potential distribution in Norway. Our results show that the SNP array successfully genotyped all included herbarium specimens. Hence, with the appropriate design procedures, the SNP array technology appears highly promising for genotyping old herbarium specimens. The captured genetic diversity correlates negatively with distance from Norway. The historical-modern comparisons reveal similar genetic structure and diversity across space and limited genetic change through time in Norway, providing no signs of any regional bottleneck (i.e., spatiotemporal stasis). The regional areas in Norway have remained genetically divergent, however, both from each other and more so from populations outside of Norway, rendering continued protection of the species in Norway relevant. The ENM results suggest that northern dragonhead has not fully achieved its potential distribution in Norway and corroborate that the species is anchored in warmer and drier habitats.

Diversity and population structure of Nordic potato cultivars and breeding clones

BACKGROUND

The genetic diversity and population structure of breeding germplasm is central knowledge for crop improvement. To gain insight into the genetic potential of the germplasm used for potato breeding in a Nordic breeding program as well as all available accessions from the Nordic genebank (NordGen), 133 potato genotypes were genotyped using the Infinium Illumina 20 K SNP array. After SNP filtering, 11 610 polymorphic SNPs were included in the analysis. In addition, data from three important breeding traits - percent dry matter and uniformity of tuber shape and eye - were scored to measure the variation potato cultivars and breeding clones.

RESULTS

The genetic diversity among the genotypes was estimated using principal coordinate analysis based on the genetic distance between individuals, as well as by using the software STRUCTURE. Both methods suggest that the collected breeding material and the germplasm from the gene-bank are closely related, with a low degree of population structure between the groups. The phenotypic distribution among the genotypes revealed significant differences, especially between farmer's cultivars and released cultivars and breeding clones. The percent heterozygosity was similar between the groups, with a mean average of 58-60%. Overall, the breeding germplasm and the accessions from the Nordic genebank seems to be closely related with similar genetic background.

CONCLUSION

The genetic potential of available Nordic potato breeding germplasm is low, and for genetic hybridization purposes, genotypes from outside the Nordic region should be employed.

Demographic history differences between Hispanics and Brazilians imprint haplotype features

Admixture is known to greatly impact the genetic landscape of a population and, while genetic variation underlying human phenotypes has been shown to differ among populations, studies on admixed subjects are still scarce. Latin American populations are the result of complex demographic history, such as 2 or 3-way admixing events, bottlenecks and/or expansions, and adaptive events unique to the American continent. To explore the impact of these events on the genetic structure of Latino populations, we evaluated the following haplotype features: linkage disequilibrium, shared identity by descent segments, runs of homozygosity, and extended haplotype homozygosity (integrated haplotype score) in Latinos represented in the 1000 Genome Project along with array data from 171 Brazilians sampled in the South and Southeast regions of Brazil. We found that linkage disequilibrium decay relates to the amount of American and African ancestry. The extent of identity by descent sharing positively correlates with historical effective population sizes, which we found to be steady or growing, except for Puerto Ricans and Colombians. Long runs of homozygosity, a particular instance of autozygosity, was only enriched in Peruvians and Native Americans. We used simulations to account for random sampling and linkage disequilibrium to filter positive selection indexes and found 244 unique markers under selection, 26 of which are common to 2 or more populations. Some markers exhibiting positive selection signals had estimated time to the most recent common ancestor consistent with human adaptation to the American continent. In conclusion, Latino populations present highly divergent haplotype characteristics that impact genetic architecture and underlie complex phenotypes.

Exploring the diversity and genetic structure of the U.S. National Cultivated Strawberry Collection

The cultivated strawberry (Fragaria ×ananassa) arose through a hybridization of two wild American octoploid strawberry species in a French garden in the 1750s. Since then, breeders have developed improved cultivars adapted to different growing regions. Diverse germplasm is crucial to meet the challenges strawberry breeders will continue to address. The USDA-ARS National Clonal Germplasm Repository (NCGR) in Corvallis, Oregon maintains the U.S. strawberry collection. Recent developments in high-throughput genotyping for strawberry can provide new insights about the diversity and structure of the collection, germplasm management, and future breeding strategies. Genotyping was conducted on 539 F. ×ananassa accessions using either the iStraw35 or FanaSNP 50 K Axiom array. Data for markers shared by the two arrays were curated for call quality, missing data, and minor allele frequency resulting in 4033 markers for structure assessment, diversity analysis, pedigree confirmation, core collection development, and the identification of haplotypes associated with desirable traits. The F. ×ananassa collection was equally diverse across the different geographic regions represented. K-means clustering, sNMF, and UPGMA hierarchal clustering revealed seven to nine sub-populations associated with different geographic breeding centers. Two 100 accession core collections were created. Pedigree linkages within the collection were confirmed. Finally, accessions containing disease resistance-associated haplotypes for FaRCa1, FaRCg1, FaRMp1, and FaRPc2 were identified. These new core collections will allow breeders and researchers to more efficiently utilize the F. ×ananassa collection. The core collections and other accessions of interest can be requested for research from the USDA-ARS NCGR via the Germplasm Resources Information Network (https://www.ars-grin.gov/).

A 20-SNP Panel as a Tool for Genetic Authentication and Traceability of Pig Breeds

Simple Summary

Given the high economic and qualitative values of local-breed meat products, it is not uncommon that substitution or mislabeling (either fraudulent or accidental) occurs at the market level. Therefore, to protect the interests of both producers and consumers, a reliable traceability tool should be developed. Nowadays, traceability usually relies on physical labeling systems (e.g., ear tags, tattoos, or electronic transponders). These systems do not, however, have good performances when dealing with carcasses or processed meat products. Molecular markers (i.e., based on the DNA sequence) can be a solution, since DNA is easily extracted from a wide variety of animal products and parts, and is not degraded during processing, even at the high temperatures involved. The aim of this study was to identify a small number of DNA mutations for breed-traceability purposes, in particular of the Italian Nero Siciliano pig and its derived products. A small panel of 12 DNA mutations was enough to discriminate Nero Siciliano pigs from other pig breeds and from wild boars.

Abstract

Food authentication in local breeds has important implications from both an economic and a qualitative point of view. Meat products from autochthonous breeds are of premium value, but can easily incur fraudulent or accidental substitution or mislabeling. The aim of this study was to identify a small number of SNPs using the Illumina PorcineSNP60 BeadChip for breed traceability, in particular of the Italian Nero Siciliano pig and its derived products. A panel of 12 SNPs was sufficient to discriminate Nero Siciliano pig from cosmopolitan breeds and wild boars. After adding 8 SNPs, the final panel of 20 SNPs allowed us to discriminate all the breeds involved in the study, to correctly assign each individual to its breed, and, moreover, to discriminate Nero Siciliano from first-generation hybrids. Almost all livestock breeds are being genotyped with medium- or high-density SNP panels, providing a large amount of information for many applications. Here, we proposed a method to select a reduced SNP panel to be used for the traceability of pig breeds.

Assessment of linkage disequilibrium patterns between structural variants and single nucleotide polymorphisms in three commercial chicken populations

BACKGROUND

Structural variants (SV) are causative for some prominent phenotypic traits of livestock as different comb types in chickens or color patterns in pigs. Their effects on production traits are also increasingly studied. Nevertheless, accurately calling SV remains challenging. It is therefore of interest, whether close-by single nucleotide polymorphisms (SNPs) are in strong linkage disequilibrium (LD) with SVs and can serve as markers. Literature comes to different conclusions on whether SVs are in LD to SNPs on the same level as SNPs to other SNPs. The present study aimed to generate a precise SV callset from whole-genome short-read sequencing (WGS) data for three commercial chicken populations and to evaluate LD patterns between the called SVs and surrounding SNPs. It is thereby the first study that assessed LD between SVs and SNPs in chickens.

RESULTS

The final callset consisted of 12,294,329 bivariate SNPs, 4,301 deletions (DEL), 224 duplications (DUP), 218 inversions (INV) and 117 translocation breakpoints (BND). While average LD between DELs and SNPs was at the same level as between SNPs and SNPs, LD between other SVs and SNPs was strongly reduced (DUP: 40%, INV: 27%, BND: 19% of between-SNP LD). A main factor for the reduced LD was the presence of local minor allele frequency differences, which accounted for 50% of the difference between SNP - SNP and DUP - SNP LD. This was potentially accompanied by lower genotyping accuracies for DUP, INV and BND compared with SNPs and DELs. An evaluation of the presence of tag SNPs (SNP in highest LD to the variant of interest) further revealed DELs to be slightly less tagged by WGS SNPs than WGS SNPs by other SNPs. This difference, however, was no longer present when reducing the pool of potential tag SNPs to SNPs located on four different chicken genotyping arrays.

CONCLUSIONS

The results implied that genomic variance due to DELs in the chicken populations studied can be captured by different SNP marker sets as good as variance from WGS SNPs, whereas separate SV calling might be advisable for DUP, INV, and BND effects.

SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding

Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.

A framework for research into continental ancestry groups of the UK Biobank

BACKGROUND

The UK Biobank is a large prospective cohort, based in the UK, that has deep phenotypic and genomic data on roughly a half a million individuals. Included in this resource are data on approximately 78,000 individuals with "non-white British ancestry." While most epidemiology studies have focused predominantly on populations of European ancestry, there is an opportunity to contribute to the study of health and disease for a broader segment of the population by making use of the UK Biobank's "non-white British ancestry" samples. Here, we present an empirical description of the continental ancestry and population structure among the individuals in this UK Biobank subset.

RESULTS

Reference populations from the 1000 Genomes Project for Africa, Europe, East Asia, and South Asia were used to estimate ancestry for each individual. Those with at least 80% ancestry in one of these four continental ancestry groups were taken forward (N = 62,484). Principal component and K-means clustering analyses were used to identify and characterize population structure within each ancestry group. Of the approximately 78,000 individuals in the UK Biobank that are of "non-white British" ancestry, 50,685, 6653, 2782, and 2364 individuals were associated to the European, African, South Asian, and East Asian continental ancestry groups, respectively. Each continental ancestry group exhibits prominent population structure that is consistent with self-reported country of birth data and geography.

CONCLUSIONS

Methods outlined here provide an avenue to leverage UK Biobank's deeply phenotyped data allowing researchers to maximize its potential in the study of health and disease in individuals of non-white British ancestry.

Increasing calling accuracy, coverage, and read-depth in sequence data by the use of haplotype blocks

High-throughput genotyping of large numbers of lines remains a key challenge in plant genetics, requiring geneticists and breeders to find a balance between data quality and the number of genotyped lines under a variety of different existing genotyping technologies when resources are limited. In this work, we are proposing a new imputation pipeline (“HBimpute”) that can be used to generate high-quality genomic data from low read-depth whole-genome-sequence data. The key idea of the pipeline is the use of haplotype blocks from the software HaploBlocker to identify locally similar lines and subsequently use the reads of all locally similar lines in the variant calling for a specific line. The effectiveness of the pipeline is showcased on a dataset of 321 doubled haploid lines of a European maize landrace, which were sequenced at 0.5X read-depth. The overall imputing error rates are cut in half compared to state-of-the-art software like BEAGLE and STITCH, while the average read-depth is increased to 83X, thus enabling the calling of copy number variation. The usefulness of the obtained imputed data panel is further evaluated by comparing the performance of sequence data in common breeding applications to that of genomic data generated with a genotyping array. For both genome-wide association studies and genomic prediction, results are on par or even slightly better than results obtained with high-density array data (600k). In particular for genomic prediction, we observe slightly higher data quality for the sequence data compared to the 600k array in the form of higher prediction accuracies. This occurred specifically when reducing the data panel to the set of overlapping markers between sequence and array, indicating that sequencing data can benefit from the same marker ascertainment as used in the array process to increase the quality and usability of genomic data.

High-throughput genotyping of large numbers of lines remains a key challenge in plant genetics and breeding. Cost, precision, and throughput must be balanced to achieve optimal efficiency given available technologies and finite resources. Although genotyping arrays are still considered the gold standard in high-throughput quantitative genetics, recent advances in sequencing provide new opportunities. Both the quality and cost of genomic data generated based on sequencing are highly dependent on the used read-depth. In this work, we propose a new imputation pipeline (“HBimpute”) that uses haplotype blocks to detect individuals of the same genetic origin and subsequently uses all reads of those individuals in the variant calling. Thus, the obtained virtual read-depth is artificially increased, leading to higher calling accuracy, coverage, and the ability to call copy number variation based on low read-depth sequencing data. To conclude, our approach makes sequencing a cost-competitive alternative to genotyping arrays with the added benefit of allowing the calling of structural variation.

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.

Exploring the legacy of Central European historical winter wheat landraces

Historical wheat landraces are rich sources of genetic diversity offering untapped reservoirs for broadening the genetic base of modern varieties. Using a 20K SNP array, we investigated the accessible genetic diversity in a Central European bread wheat landrace collection with great drought, heat stress tolerance and higher tillering capacity. We discovered distinct differences in the number of average polymorphisms between landraces and modern wheat cultivars, and identified a set of novel rare alleles present at low frequencies in the landrace collection. The detected polymorphisms were unevenly distributed along the wheat genome, and polymorphic markers co-localized with genes of great agronomic importance. The geographical distribution of the inferred Bayesian clustering revealed six genetically homogenous ancestral groups among the collection, where the Central European core bared an admixed background originating from four ancestral groups. We evaluated the effective population sizes (Ne) of the Central European collection and assessed changes in diversity over time, which revealed a dramatic ~ 97% genetic erosion between 1955 and 2015.

Informing conservation strategies with museum genomics: Long-term effects of past anthropogenic persecution on the elusive European wildcat

Like many carnivore species, European wildcats (Felis silvestris) have suffered severe anthropogenic population declines in the past, resulting in a strong population bottleneck at the beginning of the 20th century. In Germany, the species has managed to survive its near extinction in small isolated areas and is currently recolonizing former habitats owing to legal protection and concerted conservation efforts. Here, we SNP-genotyped and mtDNA-sequenced 56 historical and 650 contemporary samples to assess the impact of massive persecution on genetic diversity, population structure, and hybridization dynamics of wildcats. Spatiotemporal analyses suggest that the presumed postglacial differentiation between two genetically distinct metapopulations in Germany is in fact the result of the anthropogenic bottleneck followed by re-expansion from few secluded refugia. We found that, despite the bottleneck, populations experienced no severe genetic erosion, nor suffered from elevated inbreeding or showed signs of increased hybridization with domestic cats. Our findings have significant implications for current wildcat conservation strategies, as the data analyses show that the two presently recognized wildcat population clusters should be treated as a single conservation unit. Although current populations appear under no imminent threat from genetic factors, fostering connectivity through the implementation of forest corridors will facilitate the preservation of genetic diversity and promote long-term viability. The present study documents how museum collections can be used as essential resource for assessing long-term anthropogenic effects on natural populations, for example, regarding population structure and the delineation of appropriate conservation units, potentially informing todays' species conservation.

How imputation can mitigate SNP ascertainment Bias

Background

Population genetic studies based on genotyped single nucleotide polymorphisms (SNPs) are influenced by a non-random selection of the SNPs included in the used genotyping arrays. The resulting bias in the estimation of allele frequency spectra and population genetics parameters like heterozygosity and genetic distances relative to whole genome sequencing (WGS) data is known as SNP ascertainment bias. Full correction for this bias requires detailed knowledge of the array design process, which is often not available in practice. This study suggests an alternative approach to mitigate ascertainment bias of a large set of genotyped individuals by using information of a small set of sequenced individuals via imputation without the need for prior knowledge on the array design.

Results

The strategy was first tested by simulating additional ascertainment bias with a set of 1566 chickens from 74 populations that were genotyped for the positions of the Affymetrix Axiom™ 580 k Genome-Wide Chicken Array. Imputation accuracy was shown to be consistently higher for populations used for SNP discovery during the simulated array design process. Reference sets of at least one individual per population in the study set led to a strong correction of ascertainment bias for estimates of expected and observed heterozygosity, Wright’s Fixation Index and Nei’s Standard Genetic Distance. In contrast, unbalanced reference sets (overrepresentation of populations compared to the study set) introduced a new bias towards the reference populations. Finally, the array genotypes were imputed to WGS by utilization of reference sets of 74 individuals (one per population) to 98 individuals (additional commercial chickens) and compared with a mixture of individually and pooled sequenced populations. The imputation reduced the slope between heterozygosity estimates of array data and WGS data from 1.94 to 1.26 when using the smaller balanced reference panel and to 1.44 when using the larger but unbalanced reference panel. This generally supported the results from simulation but was less favorable, advocating for a larger reference panel when imputing to WGS.

Conclusions

The results highlight the potential of using imputation for mitigation of SNP ascertainment bias but also underline the need for unbiased reference sets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07663-6.