Genotype imputation

PTPRS is a novel marker for early Tau pathology and synaptic integrity in Alzheimer's disease

We examined the role of protein tyrosine phosphatase receptor sigma (PTPRS) in the context of Alzheimer's disease and synaptic integrity. Publicly available datasets (BRAINEAC, ROSMAP, ADC1) and a cohort of asymptomatic but "at risk" individuals (PREVENT-AD) were used to explore the relationship between PTPRS and various Alzheimer's disease biomarkers. We identified that PTPRS rs10415488 variant C shows features of neuroprotection against early Tau pathology and synaptic degeneration in Alzheimer's disease. This single nucleotide polymorphism correlated with higher PTPRS transcript abundance and lower p(181)Tau and GAP-43 levels in the CSF. In the brain, PTPRS protein abundance was significantly correlated with the quantity of two markers of synaptic integrity: SNAP25 and SYT-1. We also found the presence of sexual dimorphism for PTPRS, with higher CSF concentrations in males than females. Male carriers for variant C were found to have a 10-month delay in the onset of AD. We thus conclude that PTPRS acts as a neuroprotective receptor in Alzheimer's disease. Its protective effect is most important in males, in whom it postpones the age of onset of the disease.

Disparities in ABO blood type determination across diverse ancestries: a systematic review and validation in the All of Us Research Program

OBJECTIVES

ABO blood types have widespread clinical use and robust associations with disease. The purpose of this study is to evaluate the portability and suitability of tag single-nucleotide polymorphisms (tSNPs) used to determine ABO alleles and blood types across diverse populations in published literature.

MATERIALS AND METHODS

Bibliographic databases were searched for studies using tSNPs to determine ABO alleles. We calculated linkage between tSNPs and functional variants across inferred continental ancestry groups from 1000 Genomes. We compared r2 across ancestry and assessed real-world consequences by comparing tSNP-derived blood types to serology in a diverse population from the All of Us Research Program.

RESULTS

Linkage between functional variants and O allele tSNPs was significantly lower in African (median r2 = 0.443) compared to East Asian (r2 = 0.946, P = 1.1 × 10-5) and European (r2 = 0.869, P = .023) populations. In All of Us, discordance between tSNP-derived blood types and serology was high across all SNPs in African ancestry individuals and linkage was strongly correlated with discordance across all ancestries (ρ = -0.90, P = 3.08 × 10-23).

DISCUSSION

Many studies determine ABO blood types using tSNPs. However, tSNPs with low linkage disequilibrium promote misinference of ABO blood types, particularly in diverse populations. We observe common use of inappropriate tSNPs to determine ABO blood type, particularly for O alleles and with some tSNPs mistyping up to 58% of individuals.

CONCLUSION

Our results highlight the lack of transferability of tSNPs across ancestries and potential exacerbation of disparities in genomic research for underrepresented populations. This is especially relevant as more diverse cohorts are made publicly available.

PTPRS is a novel marker for early tau pathology and synaptic integrity in Alzheimer's disease

We examined the role of protein tyrosine phosphatase receptor sigma (PTPRS) in the context of Alzheimer's disease and synaptic integrity. Publicly available datasets (BRAINEAC, ROSMAP, ADC1) and a cohort of asymptomatic but "at risk" individuals (PREVENT-AD) were used to explore the relationship between PTPRS and various Alzheimer's disease biomarkers. We identified that PTPRS rs10415488 variant C shows features of neuroprotection against early tau pathology and synaptic degeneration in Alzheimer's disease. This single nucleotide polymorphism correlated with higher PTPRS transcript abundance and lower P-tau181 and GAP-43 levels in the CSF. In the brain, PTPRS protein abundance was significantly correlated with the quantity of two markers of synaptic integrity: SNAP25 and SYT-1. We also found the presence of sexual dimorphism for PTPRS, with higher CSF concentrations in males than females. Male carriers for variant C were found to have a 10-month delay in the onset of AD. We thus conclude that PTPRS acts as a neuroprotective receptor in Alzheimer's disease. Its protective effect is most important in males, in whom it postpones the age of onset of the disease.

Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Imputation of ancient canid genomes reveals inbreeding history over the past 10,000 years

The multi-millenia long history between dogs and humans has justly placed them at the forefront of archeological and genomic research. Despite ongoing efforts including the analysis of ancient dog and wolf genomes, many questions remain regarding their geographic and temporal origins, and the microevolutionary processes that led to the huge diversity of breeds today. Although ancient genomes provide valuable information, their use is significantly hindered by low depth of coverage and post-mortem damage, which often inhibits confident genotype calling. In the present study, we assess how genotype imputation of ancient dog and wolf genomes, utilising a large reference panel, can improve the resolution afforded by ancient genomic datasets. Imputation accuracy was evaluated by down-sampling 10 high coverage ancient and modern dog and wolf genomes to 0.05-2x coverage and comparing concordance between imputed and high coverage genotypes. We also measured the impact of imputation on principal component analyses (PCA) and runs of homozygosity (ROH). Our findings show high (R2 > 0.9) imputation accuracy for dogs with coverage as low as 0.5x and for wolves as low as 1.0x. We then imputed a worldwide dataset of 81 published ancient dog and wolf genomes, in addition to nine newly sequenced medieval and early modern period European dogs, to assess changes in inbreeding during the last 10,000 years of dog evolution. Ancient dog and wolf populations generally exhibited lower inbreeding levels than present-day individuals, though with some exceptions occurring in ancient Arctic and European dogs. Interestingly, regions with low ROH density maintained across ancient and present-day samples were significantly associated with genes related to olfaction and immune response. Our study indicates that imputing ancient canine genomes is a viable strategy that allows for the use of analytical methods previously limited to high-quality genetic data.

Assessing the impact of post-mortem damage and contamination on imputation performance in ancient DNA

Low-coverage imputation is becoming ever more present in ancient DNA (aDNA) studies. Imputation pipelines commonly used for present-day genomes have been shown to yield accurate results when applied to ancient genomes. However, post-mortem damage (PMD), in the form of C-to-T substitutions at the reads termini, and contamination with DNA from closely related species can potentially affect imputation performance in aDNA. In this study, we evaluated imputation performance (i) when using a genotype caller designed for aDNA, ATLAS, compared to bcftools, and (ii) when contamination is present. We evaluated imputation performance with principal component analyses and by calculating imputation error rates. With a particular focus on differently imputed sites, we found that using ATLAS prior to imputation substantially improved imputed genotypes for a very damaged ancient genome (42% PMD). Trimming the ends of the sequencing reads led to similar improvements in imputation accuracy. For the remaining genomes, ATLAS brought limited gains. Finally, to examine the effect of contamination on imputation, we added various amounts of reads from two present-day genomes to a previously downsampled high-coverage ancient genome. We observed that imputation accuracy drastically decreased for contamination rates above 5%. In conclusion, we recommend (i) accounting for PMD by either trimming sequencing reads or using a genotype caller such as ATLAS before imputing highly damaged genomes and (ii) only imputing genomes containing up to 5% of contamination.

The Women4Health cohort: a unique cohort to study women-specific mechanisms of cardio-metabolic regulation

Aims

Epidemiological research has shown relevant differences between sexes in clinical manifestations, severity, and progression of cardiovascular and metabolic disorders. To date, the mechanisms underlying these differences remain unknown. Given the rising incidence of such diseases, gender-specific research on established and emerging risk factors, such as dysfunction of glycaemic and/or lipid metabolism, of sex hormones and of gut microbiome, is of paramount importance. The relationships between sex hormones, gut microbiome, and host glycaemic and/or lipid metabolism are largely unknown even in the homoeostasis status. Yet this knowledge gap would be pivotal to pinpoint to key mechanisms that are likely to be disrupted in disease context.

Methods and results

Here we present the Women4Health (W4H) cohort, a unique cohort comprising up to 300 healthy women followed up during a natural menstrual cycle, set up with the primary goal to investigate the combined role of sex hormones and gut microbiota variations in regulating host lipid and glucose metabolism during homoeostasis, using a multi-omics strategy. Additionally, the W4H cohort will take into consideration another ecosystem that is unique to women, the vaginal microbiome, investigating its interaction with gut microbiome and exploring-for the first time-its role in cardiometabolic disorders.

Conclusion

The W4H cohort study lays a foundation for improving current knowledge of women-specific mechanisms in cardiometabolic regulation. It aspires to transform insights on host-microbiota interactions into prevention and therapeutic approaches for personalized health care.

The predictive capacity of polygenic risk scores for disease risk is only moderately influenced by imputation panels tailored to the target population

MOTIVATION

Polygenic risk scores (PRSs) predict individuals' genetic risk of developing complex diseases. They summarize the effect of many variants discovered in genome-wide association studies (GWASs). However, to date, large GWASs exist primarily for the European population and the quality of PRS prediction declines when applied to other ethnicities. Genetic profiling of individuals in the discovery set (on which the GWAS was performed) and target set (on which the PRS is applied) is typically done by SNP arrays that genotype a fraction of common SNPs. Therefore, a key step in GWAS analysis and PRS calculation is imputing untyped SNPs using a panel of fully sequenced individuals. The imputation results depend on the ethnic composition of the imputation panel. Imputing genotypes with a panel of individuals of the same ethnicity as the genotyped individuals typically improves imputation accuracy. However, there has been no systematic investigation into the influence of the ethnic composition of imputation panels on the accuracy of PRS predictions when applied to ethnic groups that differ from the population used in the GWAS.

RESULTS

We estimated the effect of imputation of the target set on prediction accuracy of PRS when the discovery and the target sets come from different ethnic groups. We analyzed binary phenotypes on ethnically distinct sets from the UK Biobank and other resources. We generated ethnically homogenous panels, imputed the target sets, and generated PRSs. Then, we assessed the prediction accuracy obtained from each imputation panel. Our analysis indicates that using an imputation panel matched to the ethnicity of the target population yields only a marginal improvement and only under specific conditions.

AVAILABILITY AND IMPLEMENTATION

The source code used for executing the analyses is this paper is available at https://github.com/Shamir-Lab/PRS-imputation-panels.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E-11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

Construction of a strawberry breeding core collection to capture and exploit genetic variation

BACKGROUND

Genetic diversity is crucial for the success of plant breeding programs and core collections are important resources to capture this diversity. Many core collections have already been constructed by gene banks, whose main goal is to obtain a panel of a limited number of genotypes to simplify management practices and to improve shareability while retaining as much diversity as possible. However, as gene banks have a different composition and goal than plant breeding programs, constructing a core collection for a plant breeding program should consider different aspects.

RESULTS

In this study, we present a novel approach for constructing a core collection by integrating both genomic and pedigree information to maximize the representation of the breeding germplasm in a minimum subset of genotypes while accounting for future genetic variation within a strawberry breeding program. Our stepwise approach starts with selecting the most important crossing parents of advanced selections and genotypes included for specific traits, to represent also future genetic variation. We then use pedigree-genomic-based relationship coefficients combined with the 'accession to nearest entry' criterion to complement the core collection and maximize its representativeness of the current breeding program. Combined pedigree-genomic-based relationship coefficients allow for accurate relationship estimation without the need to genotype every individual in the breeding program.

CONCLUSIONS

This stepwise construction of a core collection in a strawberry breeding program can be applied in other plant breeding programs to construct core collections for various purposes.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than five-fold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here we sought to define the roles of PAX1 and newly-identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); P=7.07e-11, OR=1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc (IVD)-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wildtype. By genetic targeting we found that wildtype Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, this suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2, or tamoxifen treatment, significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a Pax1-Col11a1-Mmp3 signaling axis in spinal chondrocytes.

Review: Computational analysis of human skeletal remains in ancient DNA and forensic genetics

Degraded DNA is used to answer questions in the fields of ancient DNA (aDNA) and forensic genetics. While aDNA studies typically center around human evolution and past history, and forensic genetics is often more concerned with identifying a specific individual, scientists in both fields face similar challenges. The overlap in source material has prompted periodic discussions and studies on the advantages of collaboration between fields toward mutually beneficial methodological advancements. However, most have been centered around wet laboratory methods (sampling, DNA extraction, library preparation, etc.). In this review, we focus on the computational side of the analytical workflow. We discuss limitations and considerations to consider when working with degraded DNA. We hope this review provides a framework to researchers new to computational workflows for how to think about analyzing highly degraded DNA and prompts an increase of collaboration between the forensic genetics and aDNA fields.

Genetic QT Score and Sleep Apnea as Predictors of Sudden Cardiac Death in the UK Biobank

Introduction

The goal of this study was to evaluate the association between a polygenic risk score (PRS) for QT prolongation (QTc-PRS), QTc intervals and mortality in patients enrolled in the UK Biobank with and without sleep apnea.

Methods

The QTc-PRS was calculated using allele copy number and previously reported effect estimates for each single nuclear polymorphism SNP. Competing-risk regression models adjusting for age, sex, BMI, QT prolonging medication, race, and comorbid cardiovascular conditions were used for sudden cardiac death (SCD) analyses.

Results

500,584 participants were evaluated (56.5 ±8 years, 54% women, 1.4% diagnosed with sleep apnea). A higher QTc-PRS was independently associated with the increased QTc interval duration (p<0.0001). The mean QTc for the top QTc-PRS quintile was 15 msec longer than the bottom quintile (p<0.001). Sleep apnea was found to be an effect modifier in the relationship between QTc-PRS and SCD. The adjusted HR per 5-unit change in QTc-PRS for SCD was 1.64 (95% CI 1.16 - 2.31, p=0.005) among those with sleep apnea and 1.04 (95% CI 0.95 - 1.14, p=0.44) among those without sleep apnea (p for interaction =0.01). Black participants with sleep apnea had significantly elevated adjusted risk of SCD compared to White participants (HR=9.6, 95% CI 1.24 - 74, p=0.03).

Conclusion

In the UK Biobank population, the QTc-PRS was associated with SCD among participants with sleep apnea but not among those without sleep apnea, indicating that sleep apnea is a significant modifier of the genetic risk. Black participants with sleep apnea had a particularly high risk of SCD.

Additive and non-additive genetic variance in juvenile Sitka spruce (Picea sitchensis Bong. Carr)

Many quantitative genetic models assume that all genetic variation is additive because of a lack of data with sufficient structure and quality to determine the relative contribution of additive and non-additive variation. Here the fractions of additive (fa) and non-additive (fd) genetic variation were estimated in Sitka spruce for height, bud burst and pilodyn penetration depth. Approximately 1500 offspring were produced in each of three sib families and clonally replicated across three geographically diverse sites. Genotypes from 1525 offspring from all three families were obtained by RADseq, followed by imputation using 1630 loci segregating in all families and mapped using the newly developed linkage map of Sitka spruce. The analyses employed a new approach for estimating fa and fd, which combined all available genotypic and phenotypic data with spatial modelling for each trait and site. The consensus estimate for fa increased with age for height from 0.58 at 2 years to 0.75 at 11 years, with only small overlap in 95% support intervals (I95). The estimated fa for bud burst was 0.83 (I95=[0.78, 0.90]) and 0.84 (I95=[0.77, 0.92]) for pilodyn depth. Overall, there was no evidence of family heterogeneity for height or bud burst, or site heterogeneity for pilodyn depth, and no evidence of inbreeding depression associated with genomic homozygosity, expected if dominance variance was the major component of non-additive variance. The results offer no support for the development of sublines for crossing within the species. The models give new opportunities to assess more accurately the scale of non-additive variation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11295-023-01627-5.

Record-matching of STR profiles with fragmentary genomic SNP data

In many forensic settings, identity of a DNA sample is sought from poor-quality DNA, for which the typical STR loci tabulated in forensic databases are not possible to reliably genotype. Genome-wide SNPs, however, can potentially be genotyped from such samples via next-generation sequencing, so that queries can in principle compare SNP genotypes from DNA samples of interest to STR genotype profiles that represent proposed matches. We use genetic record-matching to evaluate the possibility of testing SNP profiles obtained from poor-quality DNA samples to identify exact and relatedness matches to STR profiles. Using simulations based on whole-genome sequences, we show that in some settings, similar match accuracies to those seen with full coverage of the genome are obtained by genetic record-matching for SNP data that represent 5-10% genomic coverage. Thus, if even a fraction of random genomic SNPs can be genotyped by next-generation sequencing, then the potential may exist to test the resulting genotype profiles for matches to profiles consisting exclusively of nonoverlapping STR loci. The result has implications in relation to criminal justice, mass disasters, missing-person cases, studies of ancient DNA, and genomic privacy.

Animal board invited review: Practical applications of genomic information in livestock

Access to high-dimensional genomic information in many livestock species is accelerating. This has been greatly aided not only by continual reductions in genotyping costs but also an expansion in the services available that leverage genomic information to create a greater return-on-investment. Genomic information on individual animals has many uses including (1) parentage verification and discovery, (2) traceability, (3) karyotyping, (4) sex determination, (5) reporting and monitoring of mutations conferring major effects or congenital defects, (6) better estimating inbreeding of individuals and coancestry among individuals, (7) mating advice, (8) determining breed composition, (9) enabling precision management, and (10) genomic evaluations; genomic evaluations exploit genome-wide genotype information to improve the accuracy of predicting an animal's (and by extension its progeny's) genetic merit. Genomic data also provide a huge resource for research, albeit the outcome from this research, if successful, should eventually be realised through one of the ten applications already mentioned. The process for generating a genotype all the way from sample procurement to identifying erroneous genotypes is described, as are the steps that should be considered when developing a bespoke genotyping panel for practical application.

Association of PDGFRA polymorphisms with the risk of corneal astigmatism in a Japanese population

Corneal astigmatism is reportedly associated with polymorphisms of the platelet-derived growth factor receptor alpha (PDGFRA) gene region in Asian populations of Chinese, Malay, and Indian ancestry and populations of European ancestry. In this study, we investigated whether these PDGFRA polymorphisms are associated with corneal astigmatism in a Japanese population. We recruited 1,535 cases with corneal astigmatism (mean corneal cylinder power across both eyes: ≤  - 0.75 diopters [D]) and 842 controls (> - 0.75 D) to genotype 13 single-nucleotide polymorphisms (SNPs) in the PDGFRA gene region. We also performed imputation analysis in the region, with 179 imputed SNPs included in the statistical analyses. The PDGFRA SNPs were not significantly associated with the cases with corneal astigmatism ≤  - 0.75 D. However, the odds ratios (ORs) of the minor alleles of SNPs in the upstream region of PDGFRA, including rs7673984, rs4864857, and rs11133315, tended to increase according to the degree of corneal astigmatism, and these SNPs were significantly associated with the cases with corneal astigmatism ≤  - 1.25 D or ≤  - 1.50 D (Pc < 0.05, OR = 1.34-1.39). These results suggest that PDGFRA SNPs play a potential role in the development of greater corneal astigmatism.

Methods for Assessing Population Relationships and History Using Genomic Data

Genetic data contain a record of our evolutionary history. The availability of large-scale datasets of human populations from various geographic areas and timescales, coupled with advances in the computational methods to analyze these data, has transformed our ability to use genetic data to learn about our evolutionary past. Here, we review some of the widely used statistical methods to explore and characterize population relationships and history using genomic data. We describe the intuition behind commonly used approaches, their interpretation, and important limitations. For illustration, we apply some of these techniques to genome-wide autosomal data from 929 individuals representing 53 worldwide populations that are part of the Human Genome Diversity Project. Finally, we discuss the new frontiers in genomic methods to learn about population history. In sum, this review highlights the power (and limitations) of DNA to infer features of human evolutionary history, complementing the knowledge gleaned from other disciplines, such as archaeology, anthropology, and linguistics.

A whole-genome reference panel of 14,393 individuals for East Asian populations accelerates discovery of rare functional variants

Underrepresentation of non-European (EUR) populations hinders growth of global precision medicine. Resources such as imputation reference panels that match the study population are necessary to find low-frequency variants with substantial effects. We created a reference panel consisting of 14,393 whole-genome sequences including more than 11,000 Asian individuals. Genome-wide association studies were conducted using the reference panel and a population-specific genotype array of 72,298 subjects for eight phenotypes. This panel yields improved imputation accuracy of rare and low-frequency variants within East Asian populations compared with the largest reference panel. Thirty-nine previously unidentified associations were found, and more than half of the variants were East Asian specific. We discovered genes with rare protein-altering variants, including LTBP1 for height and GPR75 for body mass index, as well as putative regulatory mechanisms for rare noncoding variants with cell type-specific effects. We suggest that this dataset will add to the potential value of Asian precision medicine.

Assessing the performance of European-derived cardiometabolic polygenic risk scores in South-Asians and their interplay with family history

BACKGROUND & AIMS

We aimed to assess the performance of European-derived polygenic risk scores (PRSs) for common metabolic diseases such as coronary artery disease (CAD), obesity, and type 2 diabetes (T2D) in the South Asian (SAS) individuals in the UK Biobank. Additionally, we studied the interaction between PRS and family history (FH) in the same population.

METHODS

To calculate the PRS, we used a previously published model derived from the EUR population and applied it to the individuals of SAS ancestry from the UKB study. Each PRS was adjusted according to an individual's genotype location in the principal components (PC) space to derive an ancestry adjusted PRS (aPRS). We calculated the percentiles based on aPRS and stratified individuals into three aPRS categories: low, intermediate, and high. Considering the intermediate-aPRS percentile as a reference, we compared the low and high aPRS categories and generated the odds ratio (OR) estimates. Further, we measured the combined role of aPRS and first-degree family history (FH) in the SAS population.

RESULTS

The risk of developing severe obesity for SAS individuals was almost twofold higher for individuals with high aPRS than for those with intermediate aPRS, with an OR of 1.95 (95% CI = 1.71-2.23, P < 0.01). At the same time, the risk of severe obesity was lower in the low-aPRS group (OR = 0.60, CI = 0.53-0.67, P < 0.01). Results in the same direction were found in the EUR data, where the low-PRS group had an OR of 0.53 (95% CI = 0.51-0.56, P < 0.01) and the high-PRS group had an OR of 2.06 (95% CI = 2.00-2.12, P < 0.01). We observed similar results for CAD and T2D. Further, we show that SAS individuals with a familial history of CAD and T2D with high-aPRS are associated with a higher risk of these diseases, implying a greater genetic predisposition.

CONCLUSION

Our findings suggest that CAD, obesity, and T2D GWAS summary statistics generated predominantly from the EUR population can be potentially used to derive aPRS in SAS individuals for risk stratification. With future GWAS recruiting more SAS participants and tailoring the PRSs towards SAS ancestry, the predictive power of PRS is likely to improve further.

Evidence of causality of low body mass index on risk of adolescent idiopathic scoliosis: a Mendelian randomization study

Introduction

Adolescent idiopathic scoliosis (AIS) is a disorder with a three-dimensional spinal deformity and is a common disease affecting 1-5% of adolescents. AIS is also known as a complex disease involved in environmental and genetic factors. A relation between AIS and body mass index (BMI) has been epidemiologically and genetically suggested. However, the causal relationship between AIS and BMI remains to be elucidated.

Material and methods

Mendelian randomization (MR) analysis was performed using summary statistics from genome-wide association studies (GWASs) of AIS (Japanese cohort, 5,327 cases, 73,884 controls; US cohort: 1,468 cases, 20,158 controls) and BMI (Biobank Japan: 173430 individual; meta-analysis of genetic investigation of anthropometric traits and UK Biobank: 806334 individuals; European Children cohort: 39620 individuals; Population Architecture using Genomics and Epidemiology: 49335 individuals). In MR analyses evaluating the effect of BMI on AIS, the association between BMI and AIS summary statistics was evaluated using the inverse-variance weighted (IVW) method, weighted median method, and Egger regression (MR-Egger) methods in Japanese.

Results

Significant causality of genetically decreased BMI on risk of AIS was estimated: IVW method (Estimate (beta) [SE] = -0.56 [0.16], p = 1.8 × 10-3), weighted median method (beta = -0.56 [0.18], p = 8.5 × 10-3) and MR-Egger method (beta = -1.50 [0.43], p = 4.7 × 10-3), respectively. Consistent results were also observed when using the US AIS summary statistic in three MR methods; however, no significant causality was observed when evaluating the effect of AIS on BMI.

Conclusions

Our Mendelian randomization analysis using large studies of AIS and GWAS for BMI summary statistics revealed that genetic variants contributing to low BMI have a causal effect on the onset of AIS. This result was consistent with those of epidemiological studies and would contribute to the early detection of AIS.

Genome-wide DNA methylation association study of recent and cumulative marijuana use in middle aged adults

Marijuana is a widely used psychoactive substance in the US and medical and recreational legalization has risen over the past decade. Despite the growing number of individuals using marijuana, studies investigating the association between epigenetic factors and recent and cumulative marijuana use remain limited. We therefore investigated the association between recent and cumulative marijuana use and DNA methylation levels. Participants from the Coronary Artery Risk Development in Young Adults Study with whole blood collected at examination years (Y) 15 and Y20 were randomly selected to undergo DNA methylation profiling at both timepoints using the Illumina MethylationEPIC BeadChip. Recent use of marijuana was queried at each examination and used to estimate cumulative marijuana use from Y0 to Y15 and Y20. At Y15 (n = 1023), we observed 22 and 31 methylation markers associated (FDR P ≤ 0.05) with recent and cumulative marijuana use and 132 and 16 methylation markers at Y20 (n = 883), respectively. We replicated 8 previously reported methylation markers associated with marijuana use. We further identified 640 cis-meQTLs and 198 DMRs associated with recent and cumulative use at Y15 and Y20. Differentially methylated genes were statistically overrepresented in pathways relating to cellular proliferation, hormone signaling, and infections as well as schizophrenia, bipolar disorder, and substance-related disorders. We identified numerous methylation markers, pathways, and diseases associated with recent and cumulative marijuana use in middle-aged adults, providing additional insight into the association between marijuana use and the epigenome. These results provide novel insights into the role marijuana has on the epigenome and related health conditions.

The Bayesian lens and Bayesian blinkers

I discuss the benefits of looking through the 'Bayesian lens' (seeking a Bayesian interpretation of ostensibly non-Bayesian methods), and the dangers of wearing 'Bayesian blinkers' (eschewing non-Bayesian methods as a matter of philosophical principle). I hope that the ideas may be useful to scientists trying to understand widely used statistical methods (including confidence intervals and [Formula: see text]-values), as well as teachers of statistics and practitioners who wish to avoid the mistake of overemphasizing philosophy at the expense of practical matters. This article is part of the theme issue 'Bayesian inference: challenges, perspectives, and prospects'.

Genome-wide association analyses identified novel susceptibility loci for pulmonary embolism among Han Chinese population

BACKGROUND

A large proportion of pulmonary embolism (PE) heritability remains unexplained, particularly among the East Asian (EAS) population. Our study aims to expand the genetic architecture of PE and reveal more genetic determinants in Han Chinese.

METHODS

We conducted the first genome-wide association study (GWAS) of PE in Han Chinese, then performed the GWAS meta-analysis based on the discovery and replication stages. To validate the effect of the risk allele, qPCR and Western blotting experiments were used to investigate possible changes in gene expression. Mendelian randomization (MR) analysis was employed to implicate pathogenic mechanisms, and a polygenic risk score (PRS) for PE risk prediction was generated.

RESULTS

After meta-analysis of the discovery dataset (622 cases, 8853 controls) and replication dataset (646 cases, 8810 controls), GWAS identified 3 independent loci associated with PE, including the reported loci FGG rs2066865 (p-value = 3.81 × 10^-14), ABO rs582094 (p-value = 1.16 × 10^-10) and newly reported locus FABP2 rs1799883 (p-value = 7.59 × 10^-17). Previously reported 10 variants were successfully replicated in our cohort. Functional experiments confirmed that FABP2-A163G(rs1799883) promoted the transcription and protein expression of FABP2. Meanwhile, MR analysis revealed that high LDL-C and TC levels were associated with an increased risk of PE. Individuals with the top 10% of PRS had over a fivefold increased risk for PE compared to the general population.

CONCLUSIONS

We identified FABP2, related to the transport of long-chain fatty acids, contributing to the risk of PE and provided more evidence for the essential role of metabolic pathways in PE development.

Accuracy of haplotype estimation and whole genome imputation affects complex trait analyses in complex biobanks

Sample recruitment for research consortia, biobanks, and personal genomics companies span years, necessitating genotyping in batches, using different technologies. As marker content on genotyping arrays varies, integrating such datasets is non-trivial and its impact on haplotype estimation (phasing) and whole genome imputation, necessary steps for complex trait analysis, remains under-evaluated. Using the iPSYCH dataset, comprising 130,438 individuals, genotyped in two stages, on different arrays, we evaluated phasing and imputation performance across multiple phasing methods and data integration protocols. While phasing accuracy varied by choice of method and data integration protocol, imputation accuracy varied mostly between data integration protocols. We demonstrate an attenuation in imputation accuracy within samples of non-European origin, highlighting challenges to studying complex traits in diverse populations. Finally, imputation errors can bias association tests, reduce predictive utility of polygenic scores. Carefully optimized data integration strategies enhance accuracy and replicability of complex trait analyses in complex biobanks.

The association between a genetic variant in the SULF2 gene, metabolic parameters and vascular disease in patients at high cardiovascular risk

Clearance of triglyceride-rich lipoproteins (TRLs) is mediated by several receptors, including heparan sulfate proteoglycans (HSPGs). Sulfate glucosamine-6-O-endosulfatase-2 is a gene related to the regulation of HSPG. A variant in this gene, rs2281279, has been shown to be associated with triglycerides and insulin resistance.

Objective

To determine the relationship between rs2281279, metabolic parameters and vascular events, and type 2 diabetes mellitus (T2DM) in patients at high cardiovascular risk and whether APOE genotype modifies this relationship.

Methods

Patients (n = 4386) at high cardiovascular risk from the Utrecht Cardiovascular Cohort-Second Manifestations of Arterial Disease study were stratified according to their imputed rs2281279 genotype: AA (n = 2438), AG (n = 1642) and GG (n = 306). Effects of rs2281279 on metabolic parameters, vascular events and T2DM were analyzed with linear regression and Cox models.

Results

There was no relationship between imputed rs2281279 genotype and triglycerides, non-high-density lipoprotein (HDL)-cholesterol, insulin and quantitative insulin sensitivity check index. During a median follow-up of 11.8 (IQR, 9.3-15.5) years, 1026 cardiovascular events and 320 limb events occurred. The presence of the G allele in rs2281279 did not affect the risk of vascular events [hazard ratio (HR), 1.03; 95% confidence interval (CI), 0.94-1.14] or limb events (HR, 0.92; 95% CI, 0.77-1.10). The presence of the G allele in rs2281279 did not affect the risk of T2DM (HR, 1.09; 95% CI, 0.94-1.27). The presence of the minor G allele of rs2281279 was associated with a beneficial risk profile in ε2ε2 patients, but not in ε3ε3 patients.

Conclusions

Imputed rs2281279 genotype is not associated with metabolic parameters and does not increase the risk of vascular events or T2DM in patients at high risk for cardiovascular disease.

A Pipeline for Phasing and Genotype Imputation on Mixed Human Data (Parents-Offspring Trios and Unrelated Subjects) by Reviewing Current Methods and Software

Genotype imputation has become an essential prerequisite when performing association analysis. It is a computational technique that allows us to infer genetic markers that have not been directly genotyped, thereby increasing statistical power in subsequent association studies, which consequently has a crucial impact on the identification of causal variants. Many features need to be considered when choosing the proper algorithm for imputation, including the target sample on which it is performed, i.e., related individuals, unrelated individuals, or both. Problems could arise when dealing with a target sample made up of mixed data, composed of both related and unrelated individuals, especially since the scientific literature on this topic is not sufficiently clear. To shed light on this issue, we examined existing algorithms and software for performing phasing and imputation on mixed human data from SNP arrays, specifically when related subjects belong to trios. By discussing the advantages and limitations of the current algorithms, we identified LD-based methods as being the most suitable for reconstruction of haplotypes in this specific context, and we proposed a feasible pipeline that can be used for imputing genotypes in both phased and unphased human data.

Transfer learning for genotype-phenotype prediction using deep learning models

BACKGROUND

For some understudied populations, genotype data is minimal for genotype-phenotype prediction. However, we can use the data of some other large populations to learn about the disease-causing SNPs and use that knowledge for the genotype-phenotype prediction of small populations. This manuscript illustrated that transfer learning is applicable for genotype data and genotype-phenotype prediction.

RESULTS

Using HAPGEN2 and PhenotypeSimulator, we generated eight phenotypes for 500 cases/500 controls (CEU, large population) and 100 cases/100 controls (YRI, small populations). We considered 5 (4 phenotypes) and 10 (4 phenotypes) different risk SNPs for each phenotype to evaluate the proposed method. The improved accuracy with transfer learning for eight different phenotypes was between 2 and 14.2 percent. The two-tailed p-value between the classification accuracies for all phenotypes without transfer learning and with transfer learning was 0.0306 for five risk SNPs phenotypes and 0.0478 for ten risk SNPs phenotypes.

CONCLUSION

The proposed pipeline is used to transfer knowledge for the case/control classification of the small population. In addition, we argue that this method can also be used in the realm of endangered species and personalized medicine. If the large population data is extensive compared to small population data, expect transfer learning results to improve significantly. We show that Transfer learning is capable to create powerful models for genotype-phenotype predictions in large, well-studied populations and fine-tune these models to populations were data is sparse.

Achieving improved accuracy for imputation of ancient DNA

MOTIVATION

Genotype imputation has the potential to increase the amount of information that can be gained from the often limited biological material available in ancient samples. As many widely used tools have been developed with modern data in mind, their design is not necessarily reflective of the requirements in studies of ancient DNA. Here, we investigate if an imputation method based on the full probabilistic Li and Stephens model of haplotype frequencies might be beneficial for the particular challenges posed by ancient data.

RESULTS

We present an implementation called prophaser and compare imputation performance to two alternative pipelines that have been used in the ancient DNA community based on the Beagle software. Considering empirical ancient data downsampled to lower coverages as well as present-day samples with artificially thinned genotypes, we show that the proposed method is advantageous at lower coverages, where it yields improved accuracy and ability to capture rare variation. The software prophaser is optimized for running in a massively parallel manner and achieved reasonable runtimes on the experiments performed when executed on a GPU.

AVAILABILITY AND IMPLEMENTATION

The C++ code for prophaser is available in the GitHub repository https://github.com/scicompuu/prophaser.

SUPPLEMENTARY INFORMATION

Supplementary information is available at Bioinformatics online.

MagicalRsq: Machine-learning-based genotype imputation quality calibration

Whole-genome sequencing (WGS) is the gold standard for fully characterizing genetic variation but is still prohibitively expensive for large samples. To reduce costs, many studies sequence only a subset of individuals or genomic regions, and genotype imputation is used to infer genotypes for the remaining individuals or regions without sequencing data. However, not all variants can be well imputed, and the current state-of-the-art imputation quality metric, denoted as standard Rsq, is poorly calibrated for lower-frequency variants. Here, we propose MagicalRsq, a machine-learning-based method that integrates variant-level imputation and population genetics statistics, to provide a better calibrated imputation quality metric. Leveraging WGS data from the Cystic Fibrosis Genome Project (CFGP), and whole-exome sequence data from UK BioBank (UKB), we performed comprehensive experiments to evaluate the performance of MagicalRsq compared to standard Rsq for partially sequenced studies. We found that MagicalRsq aligns better with true R2 than standard Rsq in almost every situation evaluated, for both European and African ancestry samples. For example, when applying models trained from 1,992 CFGP sequenced samples to an independent 3,103 samples with no sequencing but TOPMed imputation from array genotypes, MagicalRsq, compared to standard Rsq, achieved net gains of 1.4 million rare, 117k low-frequency, and 18k common variants, where net gains were gained numbers of correctly distinguished variants by MagicalRsq over standard Rsq. MagicalRsq can serve as an improved post-imputation quality metric and will benefit downstream analysis by better distinguishing well-imputed variants from those poorly imputed. MagicalRsq is freely available on GitHub.

Quality Control Procedures for Genome-Wide Association Studies

Genome-wide association studies (GWAS) are being conducted at an unprecedented rate in population-based cohorts and have increased our understanding of the pathophysiology of many complex diseases. Regardless of the context, the practical utility of this information ultimately depends upon the quality of the data used for statistical analyses. Quality control (QC) procedures for GWAS are constantly evolving. Here, we enumerate some of the challenges in QC of genotyped GWAS data and describe the approaches involving genotype imputation of a sample dataset along with post-imputation quality assurance, thereby minimizing potential bias and error in GWAS results. We discuss common issues associated with QC of the GWAS data (genotyped and imputed), including data file formats, software packages for data manipulation and analysis, sex chromosome anomalies, sample identity, sample relatedness, population substructure, batch effects, and marker quality. We provide detailed guidelines along with a sample dataset to suggest current best practices and discuss areas of ongoing and future research. © 2022 Wiley Periodicals LLC.

Genome-Wide Integration of Genetic and Genomic Studies of Atopic Dermatitis: Insights into Genetic Architecture and Pathogenesis

Atopic dermatitis (AD) is a common heterogeneous, chronic, itching, and inflammatory skin disease. Genetic studies have identified multiple AD susceptibility genes. However, the genetic architecture of AD has not been elucidated. In this study, we conducted a large-scale meta-analysis of AD (35,647 cases and 1,013,885 controls) to characterize the genetic basis of AD. The heritability of AD in different datasets varied from 0.6 to 7.1%. We identified 31 previously unreported genes by integrating multiomics data. Among the 31 genes, MCL1 was identified as a potential treatment target for AD by mediating gene‒drug interactions. Tissue enrichment analyses and phenome-wide association study provided strong support for the role of the hemic and immune systems in AD. Across 1,207 complex traits and diseases, genetic correlations indicated that AD shared links with multiple respiratory phenotypes. The phenome-wide Mendelian randomization analysis (Mendelian randomization‒phenome-wide association study) revealed that the age of onset of diabetes exhibited a positive causal effect on AD (inverse-variance weighted β = 0.39, SEM = 0.09, P = 2.77 × 10^-5). Overall, these results provide important insights into the genetic architecture of AD and will lead to a more thorough and complete understanding of the molecular mechanisms underlying AD.

GWAS of Reproductive Traits in Large White Pigs on Chip and Imputed Whole-Genome Sequencing Data

Total number born (TNB), number of stillborn (NSB), and gestation length (GL) are economically important traits in pig production, and disentangling the molecular mechanisms associated with traits can provide valuable insights into their genetic structure. Genotype imputation can be used as a practical tool to improve the marker density of single-nucleotide polymorphism (SNP) chips based on sequence data, thereby dramatically improving the power of genome-wide association studies (GWAS). In this study, we applied Beagle software to impute the 50 K chip data to the whole-genome sequencing (WGS) data with average imputation accuracy (R2) of 0.876. The target pigs, 2655 Large White pigs introduced from Canadian and French lines, were genotyped by a GeneSeek Porcine 50K chip. The 30 Large White reference pigs were the key ancestral individuals sequenced by whole-genome resequencing. To avoid population stratification, we identified genetic variants associated with reproductive traits by performing within-population GWAS and cross-population meta-analyses with data before and after imputation. Finally, several genes were detected and regarded as potential candidate genes for each of the traits: for the TNB trait: NOTCH2, KLF3, PLXDC2, NDUFV1, TLR10, CDC14A, EPC2, ORC4, ACVR2A, and GSC; for the NSB trait: NUB1, TGFBR3, ZDHHC14, FGF14, BAIAP2L1, EVI5, TAF1B, and BCAR3; for the GL trait: PPP2R2B, AMBP, MALRD1, HOXA11, and BICC1. In conclusion, expanding the size of the reference population and finding an optimal imputation strategy to ensure that more loci are obtained for GWAS under high imputation accuracy will contribute to the identification of causal mutations in pig breeding.

Circulating Soluble CD163, Associations With Cardiovascular Outcomes and Mortality, and Identification of Genetic Variants in Older Individuals: The Cardiovascular Health Study

Background Monocytes/macrophages participate in cardiovascular disease. CD163 (cluster of differentiation 163) is a monocyte/macrophage receptor, and the shed sCD163 (soluble CD163) reflects monocyte/macrophage activation. We examined the association of sCD163 with incident cardiovascular disease events and performed a genome-wide association study to identify sCD163-associated variants. Methods and Results We measured plasma sCD163 in 5214 adults (aged ≥65 years, 58.7% women, 16.2% Black) of the CHS (Cardiovascular Health Study). We used Cox regression models (associations of sCD163 with incident events and mortality); median follow-up was 26 years. Genome-wide association study analyses were stratified on race. Adjusted for age, sex, and race and ethnicity, sCD163 levels were associated with all-cause mortality (hazard ratio [HR], 1.08 [95% CI, 1.04-1.12] per SD increase), cardiovascular disease mortality (HR, 1.15 [95% CI, 1.09-1.21]), incident coronary heart disease (HR, 1.10 [95% CI, 1.04-1.16]), and incident heart failure (HR, 1.18 [95% CI, 1.12-1.25]). When further adjusted (eg, cardiovascular disease risk factors), only incident coronary heart disease lost significance. In European American individuals, genome-wide association studies identified 38 variants on chromosome 2 near MGAT5 (top result rs62165726, P=3.3×10^-18),19 variants near chromosome 17 gene ASGR1 (rs55714927, P=1.5×10^-14), and 18 variants near chromosome 11 gene ST3GAL4. These regions replicated in the European ancestry ADDITION-PRO cohort, a longitudinal cohort study nested in the Danish arm of the Anglo-Danish-Dutch study of Intensive Treatment Intensive Treatment In peOple with screeNdetcted Diabetes in Primary Care. In Black individuals, we identified 9 variants on chromosome 6 (rs3129781 P=7.1×10^-9) in the HLA region, and 3 variants (rs115391969 P=4.3×10^-8) near the chromosome 16 gene MYLK3_. Conclusions Monocyte function, as measured by sCD163, may be predictive of overall and cardiovascular-specific mortality and incident heart failure.

Genome-wide association study for systemic lupus erythematosus in an egyptian population

Systemic lupus erythematosus (SLE) susceptibility has a strong genetic component. Genome-wide association studies (GWAS) across trans-ancestral populations show both common and distinct genetic variants of susceptibility across European and Asian ancestries, while many other ethnic populations remain underexplored. We conducted the first SLE GWAS on Egyptians-an admixed North African/Middle Eastern population-using 537 patients and 883 controls. To identify novel susceptibility loci and replicate previously known loci, we performed imputation-based association analysis with 6,382,276 SNPs while accounting for individual admixture. We validated the association analysis using adaptive permutation tests (n = 109). We identified a novel genome-wide significant locus near IRS1/miR-5702 (Pcorrected = 1.98 × 10-8) and eight novel suggestive loci (Pcorrected < 1.0 × 10-5). We also replicated (Pperm < 0.01) 97 previously known loci with at least one associated nearby SNP, with ITGAM, DEF6-PPARD and IRF5 the top three replicated loci. SNPs correlated (r 2 > 0.8) with lead SNPs from four suggestive loci (ARMC9, DIAPH3, IFLDT1, and ENTPD3) were associated with differential gene expression (3.5 × 10-95 < p < 1.0 × 10-2) across diverse tissues. These loci are involved in cellular proliferation and invasion-pathways prominent in lupus and nephritis. Our study highlights the utility of GWAS in an admixed Egyptian population for delineating new genetic associations and for understanding SLE pathogenesis.

Associations between forensic loci and expression levels of neighboring genes may compromise medical privacy

A set of 20 short tandem repeats (STRs) is used by the US criminal justice system to identify suspects and to maintain a database of genetic profiles for individuals who have been previously convicted or arrested. Some of these STRs were identified in the 1990s, with a preference for markers in putative gene deserts to avoid forensic profiles revealing protected medical information. We revisit that assumption, investigating whether forensic genetic profiles reveal information about gene-expression variation or potential medical information. We find six significant correlations (false discovery rate = 0.23) between the forensic STRs and the expression levels of neighboring genes in lymphoblastoid cell lines. We explore possible mechanisms for these associations, showing evidence compatible with forensic STRs causing expression variation or being in linkage disequilibrium with a causal locus in three cases and weaker or potentially spurious associations in the other three cases. Together, these results suggest that forensic genetic loci may reveal expression levels and, perhaps, medical information.

Population structure discovery in meta-analyzed microbial communities and inflammatory bowel disease using MMUPHin

Microbiome studies of inflammatory bowel diseases (IBD) have achieved a scale for meta-analysis of dysbioses among populations. To enable microbial community meta-analyses generally, we develop MMUPHin for normalization, statistical meta-analysis, and population structure discovery using microbial taxonomic and functional profiles. Applying it to ten IBD cohorts, we identify consistent associations, including novel taxa such as Acinetobacter and Turicibacter, and additional exposure and interaction effects. A single gradient of dysbiosis severity is favored over discrete types to summarize IBD microbiome population structure. These results provide a benchmark for characterization of IBD and a framework for meta-analysis of any microbial communities.

Rapid, Reference-Free human genotype imputation with denoising autoencoders

Genotype imputation is a foundational tool for population genetics. Standard statistical imputation approaches rely on the co-location of large whole-genome sequencing-based reference panels, powerful computing environments, and potentially sensitive genetic study data. This results in computational resource and privacy-risk barriers to access to cutting-edge imputation techniques. Moreover, the accuracy of current statistical approaches is known to degrade in regions of low and complex linkage disequilibrium. Artificial neural network-based imputation approaches may overcome these limitations by encoding complex genotype relationships in easily portable inference models. Here, we demonstrate an autoencoder-based approach for genotype imputation, using a large, commonly used reference panel, and spanning the entirety of human chromosome 22. Our autoencoder-based genotype imputation strategy achieved superior imputation accuracy across the allele-frequency spectrum and across genomes of diverse ancestry, while delivering at least fourfold faster inference run time relative to standard imputation tools.

Assessment of genotyping array performance for genome-wide association studies and imputation in African cattle

BACKGROUND

In cattle, genome-wide association studies (GWAS) have largely focused on European or Asian breeds, using genotyping arrays that were primarily designed for European cattle. Because there is growing interest in performing GWAS in African breeds, we have assessed the performance of 23 commercial bovine genotyping arrays for capturing the diversity across African breeds and performing imputation. We used 409 whole-genome sequences (WGS) spanning global cattle breeds, and a real cohort of 2481 individuals (including African breeds) that were genotyped with the Illumina high-density (HD) array and the GeneSeek bovine 50 k array.

RESULTS

We found that commercially available arrays were not effective in capturing variants that segregate among African indicine animals. Only 6% of these variants in high linkage disequilibrium (LD) (r² > 0.8) were on the best performing arrays, which contrasts with the 17% and 25% in African and European taurine cattle, respectively. However, imputation from available HD arrays can successfully capture most variants (accuracies up to 0.93), mainly when using a global, not continent-specific, reference panel, which partially reflects the unusually high levels of admixture on the continent. When considering functional variants, the GGPF250 array performed best for tagging WGS variants and imputation. Finally, we show that imputation from low-density arrays can perform almost as well as HD arrays, if a two-stage imputation approach is adopted, i.e. first imputing to HD and then to WGS, which can potentially reduce the costs of GWAS.

CONCLUSIONS

Our results show that the choice of an array should be based on a balance between the objective of the study and the breed/population considered, with the HD and BOS1 arrays being the best choice for both taurine and indicine breeds when performing GWAS, and the GGPF250 being preferable for fine-mapping studies. Moreover, our results suggest that there is no advantage to using the indicus-specific arrays for indicus breeds, regardless of the objective. Finally, we show that using a reference panel that better represents global bovine diversity improves imputation accuracy, particularly for non-European taurine populations.

Extent to which array genotyping and imputation with large reference panels approximate deep whole-genome sequencing

Understanding the genetic basis of human diseases and traits is dependent on the identification and accurate genotyping of genetic variants. Deep whole-genome sequencing (WGS), the gold standard technology for SNP and indel identification and genotyping, remains very expensive for most large studies. Here, we quantify the extent to which array genotyping followed by genotype imputation can approximate WGS in studies of individuals of African, Hispanic/Latino, and European ancestry in the US and of Finnish ancestry in Finland (a population isolate). For each study, we performed genotype imputation by using the genetic variants present on the Illumina Core, OmniExpress, MEGA, and Omni 2.5M arrays with the 1000G, HRC, and TOPMed imputation reference panels. Using the Omni 2.5M array and the TOPMed panel, ≥90% of bi-allelic single-nucleotide variants (SNVs) are well imputed (r² > 0.8) down to minor-allele frequencies (MAFs) of 0.14% in African, 0.11% in Hispanic/Latino, 0.35% in European, and 0.85% in Finnish ancestries. There was little difference in TOPMed-based imputation quality among the arrays with >700k variants. Individual-level imputation quality varied widely between and within the three US studies. Imputation quality also varied across genomic regions, producing regions where even common (MAF > 5%) variants were consistently not well imputed across ancestries. The extent to which array genotyping and imputation can approximate WGS therefore depends on reference panel, genotype array, sample ancestry, and genomic location. Imputation quality by variant or genomic region can be queried with our new tool, RsqBrowser, now deployed on the Michigan Imputation Server.

Evaluation of vicinity-based hidden Markov models for genotype imputation

BACKGROUND

The decreasing cost of DNA sequencing has led to a great increase in our knowledge about genetic variation. While population-scale projects bring important insight into genotype-phenotype relationships, the cost of performing whole-genome sequencing on large samples is still prohibitive. In-silico genotype imputation coupled with genotyping-by-arrays is a cost-effective and accurate alternative for genotyping of common and uncommon variants. Imputation methods compare the genotypes of the typed variants with the large population-specific reference panels and estimate the genotypes of untyped variants by making use of the linkage disequilibrium patterns. Most accurate imputation methods are based on the Li-Stephens hidden Markov model, HMM, that treats the sequence of each chromosome as a mosaic of the haplotypes from the reference panel.

RESULTS

Here we assess the accuracy of vicinity-based HMMs, where each untyped variant is imputed using the typed variants in a small window around itself (as small as 1 centimorgan). Locality-based imputation is used recently by machine learning-based genotype imputation approaches. We assess how the parameters of the vicinity-based HMMs impact the imputation accuracy in a comprehensive set of benchmarks and show that vicinity-based HMMs can accurately impute common and uncommon variants.

CONCLUSIONS

Our results indicate that locality-based imputation models can be effectively used for genotype imputation. The parameter settings that we identified can be used in future methods and vicinity-based HMMs can be used for re-structuring and parallelizing new imputation methods. The source code for the vicinity-based HMM implementations is publicly available at https://github.com/harmancilab/LoHaMMer .

Exploring the optimal strategy of imputation from SNP array to whole-genome sequencing data in farm animals

Genotype imputation from BeadChip to whole-genome sequencing (WGS) data is a cost-effective method of obtaining genotypes of WGS variants. Beagle, one of the most popular imputation software programs, has been widely used for genotype inference in humans and non-human species. A few studies have systematically and comprehensively compared the performance of beagle versions and parameter settings of farm animals. Here, we investigated the imputation performance of three representative versions of Beagle (Beagle 4.1, Beagle 5.0, and Beagle 5.4), and the effective population size (Ne) parameter setting for three species (cattle, pig, and chicken). Six scenarios were investigated to explore the impact of certain key factors on imputation performance. The results showed that the default Ne (1,000,000) is not suitable for livestock and poultry in small reference or low-density arrays of target panels, with 2.47%–10.45% drops in accuracy. Beagle 5 significantly reduced the computation time (4.66-fold–13.24-fold) without an accuracy loss. In addition, using a large combined-reference panel or high-density chip provides greater imputation accuracy, especially for low minor allele frequency (MAF) variants. Finally, a highly significant correlation in the measures of imputation accuracy can be obtained with an MAF equal to or greater than 0.05.

Mining alleles for tar spot complex resistance from CIMMYT's maize Germplasm Bank

The tar spot complex (TSC) is a devastating disease of maize (Zea mays L.), occurring in 17 countries throughout Central, South, and North America and the Caribbean, and can cause grain yield losses of up to 80%. As yield losses from the disease continue to intensify in Central America, Phyllachora maydis, one of the causal pathogens of TSC, was first detected in the United States in 2015, and in 2020 in Ontario, Canada. Both the distribution and yield losses due to TSC are increasing, and there is a critical need to identify the genetic resources for TSC resistance. The Seeds of Discovery Initiative at CIMMYT has sought to combine next-generation sequencing technologies and phenotypic characterization to identify valuable alleles held in the CIMMYT Germplasm Bank for use in germplasm improvement programs. Individual landrace accessions of the “Breeders' Core Collection” were crossed to CIMMYT hybrids to form 918 unique accessions topcrosses (F1 families) which were evaluated during 2011 and 2012 for TSC disease reaction. A total of 16 associated SNP variants were identified for TSC foliar leaf damage resistance and increased grain yield. These variants were confirmed by evaluating the TSC reaction of previously untested selections of the larger F1 testcross population (4,471 accessions) based on the presence of identified favorable SNPs. We demonstrated the usefulness of mining for donor alleles in Germplasm Bank accessions for newly emerging diseases using genomic variation in landraces.

Genetic variants in DDO and PEX5L in peroxisome-related pathways predict non-small cell lung cancer survival

Peroxisomes play a role in lipid metabolism and regulation of reactive oxygen species, but its role in development and progression of non-small cell lung cancer (NSCLC) is not well understood. Here, we investigated the associations between 9708 single-nucleotide polymorphisms (SNPs) in 113 genes in the peroxisome-related pathways and survival of NSCLC patients from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial (PLCO) and the Harvard Lung Cancer Susceptibility (HLCS) study. In 1185 NSCLC patients from the PLCO trial, we found that 213 SNPs were significantly associated with NSCLC overall survival (OS) (p ≤ 0.05, Bayesian false discovery probability [BFDP] ≤ 0.80), of which eight SNPs were validated in the HLCS data set. In a multivariate Cox proportional hazards regression model, two independent SNPs (rs9384742 DDO and rs9825224 PEX5L) were significantly associated with NSCLC survival (hazards ratios [HR] of 1.17 with 95% CI [confidence interval] of 1.06-1.28 and 0.86 with 95% CI of 0.77-0.96, respectively). Patients with one or two protective genotypes had a significantly higher OS (HR: 0.787 [95% CI: 0.620-0.998] and 0.691 [95% CI: 0.543-0.879], respectively). Further expression quantitative trait loci analysis using whole blood and lung tissue showed that the minor allele of rs9384742 DDO was significantly associated with decreased messenger RNA (mRNA) expression levels and that DDO expression was also decreased in NSCLC tumor tissue. Additionally, high PEX5L expression levels were significantly associated with lower survival of NSCLC. Our data suggest that variants in these peroxisome-related genes may influence gene regulation and are potential predictors of NSCLC OS, once validated by additional studies.

Genotype calling and haplotype inference from low coverage sequence data in heterozygous plant genome using HetMap

This study developed a new genotyping method that can accurately infer heterozygous genotype information from the complex plant genome sequence data, which helped discover new alleles in the association studies. Many software packages and pipelines had been developed to handle the sequence data of the model species. However, Genotyping from complex heterozygous plant genome needs further improvement on the previous methods. Here we present a new pipeline available at https://github.com/Ncgrhg/HetMapv1 ) for variant calling and missing genotype imputation from low coverage sequence data for heterozygous plant genomes. To check the performance of the HetMap on the real sequence data, HetMap was applied to both the F₁ hybrid rice population, which consists of 1495 samples and the wild rice population with 446 samples. The high coverage sequence data of four hybrid rice accessions and two wild rice accessions, which were also included in low coverage sequence data, were used to validate the accuracy of genotype inference. The validation results showed that HetMap archieved significant improvement in heterozygous genotype inference accuracy (13.65% for hybrid rice, 26.05% for wild rice) and total accuracy compared with similar software packages. The application of the new genotype with the genome-wide association study also showed improvement of association power in wild rice awn length phenotype. It could archive high genotype inference accuracy in low sequence coverage in a small population with both the natural and constructed recombination population. HetMap provided a powerful tool for the heterozygous plant genome sequence data analysis, which may help to discover new phenotype regions for the plant species with the complex heterozygous genome.

Aquaculture Molecular Breeding Platform (AMBP): a comprehensive web server for genotype imputation and genetic analysis in aquaculture

It is of vital importance to understand the population structure, dissect the genetic bases of performance traits, and make proper strategies for selection in breeding programs. However, there is no single webserver covering the specific needs in aquaculture. We present Aquaculture Molecular Breeding Platform (AMBP), the first web server for genetic data analysis in aquatic species of farming interest. AMBP integrates the haplotype reference panels of 18 aquaculture species, which greatly improves the accuracy of genotype imputation. It also supports multiple tools to infer genetic structures, dissect the genetic architecture of performance traits, estimate breeding values, and predict optimum contribution. All the tools are coherently linked in a web-interface for users to generate interpretable results and evaluate statistical appropriateness. The webserver supports standard VCF and PLINK (PED, MAP) files, and implements automated pipelines for format transformation and visualization to simplify the process of analysis. As a demonstration, we applied the webserver to Pacific white shrimp and Atlantic salmon datasets. In summary, AMBP constitutes comprehensive resources and analytical tools for exploring genetic data and guiding practical breeding programs. AMBP is available at http://mgb.qnlm.ac.

Imputation of non-genotyped F1 dams to improve genetic gain in swine crossbreeding programs

This study investigated using imputed genotypes from non-genotyped animals which were not in the pedigree for the purpose of genetic selection and improving genetic gain for economically relevant traits. Simulations were used to mimic a 3-breed crossbreeding system that resembled a modern swine breeding scheme. The simulation consisted of three purebred (PB) breeds A, B, and C each with 25 and 425 mating males and females, respectively. Males from A and females from B were crossed to produce AB females (n = 1,000), which were crossed with males from C to produce crossbreds (CB; n = 10,000). The genome consisted of three chromosomes with 300 quantitative trait loci and ~9,000 markers. Lowly heritable reproductive traits were simulated for A, B, and AB (h2 = 0.2, 0.2, and 0.15, respectively), whereas a moderately heritable carcass trait was simulated for C (h2 = 0.4). Genetic correlations between reproductive traits in A, B, and AB were moderate (rg = 0.65). The goal trait of the breeding program was AB performance. Selection was practiced for four generations where AB and CB animals were first produced in generations 1 and 2, respectively. Non-genotyped AB dams were imputed using FImpute beginning in generation 2. Genotypes of PB and CB were used for imputation. Imputation strategies differed by three factors: 1) AB progeny genotyped per generation (2, 3, 4, or 6), 2) known or unknown mates of AB dams, and 3) genotyping rate of females from breeds A and B (0% or 100%). PB selection candidates from A and B were selected using estimated breeding values for AB performance, whereas candidates from C were selected by phenotype. Response to selection using imputed genotypes of non-genotyped animals was then compared to the scenarios where true AB genotypes (trueGeno) or no AB genotypes/phenotypes (noGeno) were used in genetic evaluations. The simulation was replicated 20 times. The average increase in genotype concordance between unknown and known sire imputation strategies was 0.22. Genotype concordance increased as the number of genotyped CB increased with little additional gain beyond 9 progeny. When mates of AB were known and more than 4 progeny were genotyped per generation, the phenotypic response in AB did not differ (P > 0.05) from trueGeno yet was greater (P < 0.05) than noGeno. Imputed genotypes of non-genotyped animals can be used to increase performance when 4 or more progeny are genotyped and sire pedigrees of CB animals are known.

Genotype imputation for soybean nested association mapping population to improve precision of QTL detection

KEY MESSAGE

Software for high imputation accuracy in soybean was identified. Imputed dataset could significantly reduce the interval of genomic regions controlling traits, thus greatly improve the efficiency of candidate gene identification. Genotype imputation is a strategy to increase marker density of existing datasets without additional genotyping. We compared imputation performance of software BEAGLE 5.0, IMPUTE 5 and AlphaPlantImpute and tested software parameters that may help to improve imputation accuracy in soybean populations. Several factors including marker density, extent of linkage disequilibrium (LD), minor allele frequency (MAF), etc., were examined for their effects on imputation accuracy across different software. Our results showed that AlphaPlantImpute had a higher imputation accuracy than BEAGLE 5.0 or IMPUTE 5 tested in each soybean family, especially if the study progeny were genotyped with an extremely low number of markers. LD extent, MAF and reference panel size were positively correlated with imputation accuracy, a minimum number of 50 markers per chromosome and MAF of SNPs > 0.2 in soybean line were required to avoid a significant loss of imputation accuracy. Using the software, we imputed 5176 soybean lines in the soybean nested mapping population (NAM) with high-density markers of the 40 parents. The dataset containing 423,419 markers for 5176 lines and 40 parents was deposited at the Soybase. The imputed NAM dataset was further examined for the improvement of mapping quantitative trait loci (QTL) controlling soybean seed protein content. Most of the QTL identified were at identical or at similar position based on initial and imputed datasets; however, QTL intervals were greatly narrowed. The resulting genotypic dataset of NAM population will facilitate QTL mapping of traits and downstream applications. The information will also help to improve genotyping imputation accuracy in self-pollinated crops.