Second-generation PLINK: rising to the challenge of larger and richer datasets

Genome-Wide Association Studies of Growth Trait Heterosis in Crossbred Meat Rabbits

The application of heterosis can not only effectively improve the disease resistance and meat quality of livestock, but also significantly enhance the reproduction and growth of livestock and poultry. We conducted genome-wide association studies using data from F2 crossbred meat rabbits to screen out candidate genes with significant dominant effects associated with economic trait variation. High-throughput sequencing technology was used to obtain SNPs covering the whole genome to evaluate the homozygosity of the population genome, and analyze the number, length, frequency, and distribution of ROHs in the population. Candidate genes related to economic traits of meat rabbits were searched based on high-frequency ROH regions. After quality control filtering of genotype data, 380 F2 crossbred rabbits were identified with 78,579 SNPs and 42,018 ROHs on the autosomes. The fitting of the Logistic growth curve model showed that 49-day-old rabbits were a growth inflection point. Then, through genome-wide association studies, 10 SNP loci and seven growth trait candidate genes were found to be significantly related to body weight in meat rabbits at 84 days of age. In addition, we revealed the functional roles and locations of 20 candidate genes in the high-frequency ROH region associated with economic traits in meat rabbits. This study identified potential genes associated with growth and development in the high-frequency ROH region of meat rabbits. In this study, the identified candidate genes can be used as molecular markers for assisted selection in meat rabbits. At the same time, the inbreeding situation based on ROH assessment can provide reference for breeding and breeding preservation of meat rabbits.

Chromosome X-Wide Common Variant Association Study (XWAS) in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) displays a notable male bias in prevalence. Research into rare (<0.1) genetic variants on the X chromosome has implicated over 20 genes in ASD pathogenesis, such as MECP2, DDX3X, and DMD. The "female protective effect" in ASD suggests that females may require a higher genetic burden to manifest similar symptoms as males, yet the mechanisms remain unclear. Despite technological advances in genomics, the complexity of the biological nature of sex chromosomes leave them underrepresented in genome-wide studies. Here, we conducted an X chromosome-wide association study (XWAS) using whole-genome sequencing data from 6,873 individuals with ASD (82% males) across Autism Speaks MSSNG, Simons Simplex Cohort SSC, and Simons Foundation Powering Autism Research SPARK, alongside 8,981 population controls (43% males). We analyzed 418,652 X-chromosome variants, identifying 59 associated with ASD (p-values 7.9×10-6 to 1.51×10-5), surpassing Bonferroni-corrected thresholds. Key findings include significant regions on chrXp22.2 (lead SNP=rs12687599, p=3.57×10-7) harboring ASB9/ASB11, and another encompassing DDX53/PTCHD1-AS long non-coding RNA (lead SNP=rs5926125, p=9.47×10-6). When mapping genes within 10kb of the 59 most significantly associated SNPs, 91 genes were found, 17 of which yielded association with ASD (GRPR, AP1S2, DDX53, HDAC8, PCDH19, PTCHD1, PCDH11X, PTCHD1-AS, DMD, SYAP1, CNKSR2, GLRA2, OFD1, CDKL5, GPRASP2, NXF5, SH3KBP1). FGF13 emerged as a novel X-linked ASD candidate gene, highlighted by sex-specific differences in minor allele frequencies. These results reveal significant new insights into X chromosome biology in ASD, confirming and nominating genes and pathways for further investigation.

Landscape of genomic structural variations in Indian population-based cohorts: Deeper insights into their prevalence and clinical relevance

Structural variations (SV) are large (>50 base pairs) genomic rearrangements comprising deletions, duplications, insertions, inversions, and translocations. Studying SVs is important because they play active and critical roles in regulating gene expression, determining disease predispositions, and identifying population-specific differences among individuals of diverse ancestries. However, SV discoveries in the Indian population using whole-genome sequencing (WGS) have been limited. In this study, using short-read WGS having an average 42X depth of coverage, we identify and characterize 36,210 SVs from 529 individuals enrolled in population-based cohorts in India. These SVs include 24,574 deletions, 2,913 duplications, 8,710 insertions, and 13 inversions; 1.26% (456 out of 36,210) of the identified SVs can potentially impact the coding regions of genes. Furthermore, 56 of these SVs are highly intolerant to loss-of-function changes to the mapped genes, and five SVs impacting ADAMTS17, CCDC40, and RHCE are common in our study individuals. Seven rare SVs significantly impact dosage sensitivity of genes known to be associated with various clinical phenotypes. Most of the SVs in our study are rare and heterozygous. This fine-scale SV discovery in the underrepresented Indian population provides valuable insights that extend beyond Eurocentric human genetic studies.

Epigenetic associations of GPNMB rs199347 variant with alcohol consumption in Parkinson's disease

Introduction

Alcohol consumption can induce a neuroinflammatory response and contribute to the progression of neurodegeneration. However, its association with Parkinson's disease (PD), the second most common neurodegenerative disorder, remains undetermined. Recent studies suggest that the glycoprotein non-metastatic melanoma protein B (GPNMB) is a potential biomarker for PD. We evaluated the association of rs199347, a variant of the GPNMB gene, with alcohol consumption and methylation upstream of GPNMB.

Methods

We retrieved genetic and DNA methylation data obtained from participants enrolled in the Taiwan Biobank (TWB) between 2008 and 2016. After excluding individuals with incomplete or missing information about potential PD risk factors, we included 1,357 participants in our final analyses. We used multiple linear regression to assess the association of GPNMB rs199347 and chronic alcohol consumption (and other potential risk factors) with GPNMB cg17274742 methylation.

Results

There was no difference between the distribution of GPNMB rs199347 genotypes between chronic alcohol consumers and the other study participants. A significant interaction was observed between the GPNMB rs199347 variant and alcohol consumption (p = 0.0102) concerning cg17274742 methylation. Compared to non-chronic alcohol consumers with the AA genotype, alcohol drinkers with the rs199347 GG genotype had significantly lower levels (hypomethylation) of cg17274742 (p = 0.0187).

Conclusion

Alcohol consumption among individuals with the rs199347 GG genotype was associated with lower levels of cg17274742 methylation, which could increase expression of the GPNMB gene, an important neuroinflammatory-related risk gene for PD.

Distinct positions of genetic and oral histories: Perspectives from India

Over the past decade, genomic data have contributed to several insights on global human population histories. These studies have been met both with interest and critically, particularly by populations with oral histories that are records of their past and often reference their origins. While several studies have reported concordance between oral and genetic histories, there is potential for tension that may stem from genetic histories being prioritized or used to confirm community-based knowledge and ethnography, especially if they differ. To investigate the interplay between oral and genetic histories, we focused on the southwestern region of India and analyzed whole-genome sequence data from 156 individuals identifying as Bunt, Kodava, Nair, and Kapla. We supplemented limited anthropological records on these populations with oral history accounts from community members and historical literature, focusing on references to non-local origins such as the ancient Scythians in the case of Bunt, Kodava, and Nair, members of Alexander the Great's army for the Kodava, and an African-related source for Kapla. We found these populations to be genetically most similar to other Indian populations, with the Kapla more similar to South Indian tribal populations that maximize a genetic ancestry related to Ancient Ancestral South Indians. We did not find evidence of additional genetic sources in the study populations than those known to have contributed to many other present-day South Asian populations. Our results demonstrate that oral and genetic histories may not always provide consistent accounts of population origins and motivate further community-engaged, multi-disciplinary investigations of non-local origin stories in these communities.

A robust cis-Mendelian randomization method with application to drug target discovery

Mendelian randomization (MR) uses genetic variants as instrumental variables (IVs) to investigate causal relationships between traits. Unlike conventional MR, cis-MR focuses on a single genomic region using only cis-SNPs. For example, using cis-pQTLs for a protein as exposure for a disease opens a cost-effective path for drug target discovery. However, few methods effectively handle pleiotropy and linkage disequilibrium (LD) of cis-SNPs. Here, we propose cisMR-cML, a method based on constrained maximum likelihood, robust to IV assumption violations with strong theoretical support. We further clarify the severe but largely neglected consequences of the current practice of modeling marginal, instead of conditional genetic effects, and only using exposure-associated SNPs in cis-MR analysis. Numerical studies demonstrated our method's superiority over other existing methods. In a drug-target analysis for coronary artery disease (CAD), including a proteome-wide application, we identified three potential drug targets, PCSK9, COLEC11 and FGFR1 for CAD.

Replication of Known and Identification of Novel Associations in Biobank-Scale Datasets: A Survey Using UK Biobank and FinnGen

Over the last two decades, numerous genome-wide association studies (GWAS) have been performed to unveil the genetic architecture of human complex traits. Despite multiple efforts aimed at the trans-biobank integration of GWAS results, no systematic analysis of the variant-level properties affecting the replication of known associations (or identifying novel ones) in genome-wide meta-analysis has yet been performed using biobank-scale data. To address this issue, we performed a systematic comparison of GWAS summary statistics for 679 complex traits in the UK Biobank (UKB) and FinnGen (FG) cohorts. We identified 37,148 index variants with genome-wide associations with at least one trait in either cohort or in the meta-analysis, only 3528 (9.5%) of which were shared between UKB and FG. Nearly twice as many variants (6577) were replicated in another dataset at the significance level adjusted for the number of variants selected for replication. However, as many as 9230 loci failed to be replicated. Moreover, as many as 5813 loci were observed as significant associations only in meta-analysis results, highlighting the importance of trans-biobank meta-analysis efforts. We showed that variants that failed to replicate in UKB or FG tend to correspond to rare, less pleiotropic variants with lower effect sizes and lower LD score values. Genome-wide associations specific to meta-analysis were also enriched in low-effect variants; however, such variants tended to be more common and have more consistent frequencies between populations. Taken together, our results show a relatively high rate of non-replication of genome-wide associations in the studied cohorts and highlight both widely appreciated and less acknowledged properties of the associations affecting their identification and replication.

Common and rare genetic variants predisposing females to unexplained recurrent pregnancy loss

Recurrent pregnancy loss (RPL) is a major reproductive health issue with multifactorial causes, affecting 2.6% of all pregnancies worldwide. Nearly half of the RPL cases lack clinically identifiable causes (e.g., antiphospholipid syndrome, uterine anomalies, and parental chromosomal abnormalities), referred to as unexplained RPL (uRPL). Here, we perform a genome-wide association study focusing on uRPL in 1,728 cases and 24,315 female controls of Japanese ancestry. We detect significant associations in the major histocompatibility complex (MHC) region at 6p21 (lead variant=rs9263738; P = 1.4 × 10-10; odds ratio [OR] = 1.51 [95% CI: 1.33-1.72]; risk allele frequency = 0.871). The MHC associations are fine-mapped to the classical HLA alleles, HLA-C*12:02, HLA-B*52:01, and HLA-DRB1*15:02 (P = 1.1 × 10-10, 1.5 × 10-10, and 1.2 × 10-9, respectively), which constitute a population-specific common long-range haplotype with a protective effect (P = 2.8 × 10-10; OR = 0.65 [95% CI: 0.57-0.75]; haplotype frequency=0.108). Genome-wide copy-number variation (CNV) calling demonstrates rare predicted loss-of-function (pLoF) variants of the cadherin-11 gene (CDH11) conferring the risk of uRPL (P = 1.3 × 10-4; OR = 3.29 [95% CI: 1.78-5.76]). Our study highlights the importance of reproductive immunology and rare variants in the uRPL etiology.

Small-group originating model: Optimized individual-level GWAS simulation featured by SLiM and using open-access data

The development of analytical methods for Genome-wide Association Studies (GWAS) has outpaced the evolution of simulation techniques and pipelines. This disparity underscores the importance of innovative simulation methods that can keep pace with the rapidly increasing scale of GWAS. The median sample size of GWAS over the past ten years has exceeded 50,000 individuals, a trend that emphasizes the need for simulation tools capable of generating data on a similar or larger scale. This paper introduces a novel method, the small-group originating (SGO) model, utilizing the SLiM software for simulating individual-level GWAS data. Our standardized protocol facilitates the generation of tens of thousands of pseudo-individuals with millions of variants from small (30-90) open-access datasets. SGO stands out, especially when compared to the widely-used resampling method in HapGen, showcasing superior simulation efficiency for large sample sizes (> 13,000) of unrelated individuals. This capability is particularly relevant given the current trajectory towards larger GWAS, necessitating tools that can simulate datasets reflective of this growth. Additionally, SGO provides customization options and can model dynamic life cycles and mating across generations, positioning it as a highly promising alternative for GWAS simulations. In a case study, sensitivity analyses of chromosome-level principal component analysis and kinship coefficient estimation were conducted. The results highlighted the poor robustness of chromosome-level quality control (QC) indexes and the uneven distribution of population structure across chromosomes and ancestries, advocating for the caution against relying solely on chromosome-level QC statistics. With its flexible and efficient approach to generating pseudo GWAS data, our standardized SGO protocol emerges as a crucial asset for method development, power analysis, and benchmarking in GWAS research. It is especially vital in the context of accommodating the demands for large-scale simulations, aligning with the current and future scale of GWAS.

Large-scale whole-exome sequencing analyses identified protein-coding variants associated with immune-mediated diseases in 350,770 adults

The genetic contribution of protein-coding variants to immune-mediated diseases (IMDs) remains underexplored. Through whole exome sequencing of 40 IMDs in 350,770 UK Biobank participants, we identified 162 unique genes in 35 IMDs, among which 124 were novel genes. Several genes, including FLG which is associated with atopic dermatitis and asthma, showed converging evidence from both rare and common variants. 91 genes exerted significant effects on longitudinal outcomes (interquartile range of Hazard Ratio: 1.12-5.89). Mendelian randomization identified five causal genes, of which four were approved drug targets (CDSN, DDR1, LTA, and IL18BP). Proteomic analysis indicated that mutations associated with specific IMDs might also affect protein expression in other IMDs. For example, DXO (celiac disease-related gene) and PSMB9 (alopecia areata-related gene) could modulate CDSN (autoimmune hypothyroidism-, psoriasis-, asthma-, and Graves' disease-related gene) expression. Identified genes predominantly impact immune and biochemical processes, and can be clustered into pathways of immune-related, urate metabolism, and antigen processing. Our findings identified protein-coding variants which are the key to IMDs pathogenesis and provided new insights into tailored innovative therapies.

Genome-Wide Association Study Points to Novel Locus for Gilles de la Tourette Syndrome

BACKGROUND

Tourette syndrome (TS) is a childhood-onset neurodevelopmental disorder of complex genetic architecture and is characterized by multiple motor tics and at least one vocal tic persisting for more than 1 year.

METHODS

We performed a genome-wide meta-analysis integrating a novel TS cohort with previously published data, resulting in a sample size of 6133 individuals with TS and 13,565 ancestry-matched control participants.

RESULTS

We identified a genome-wide significant locus on chromosome 5q15. Integration of expression quantitative trait locus, Hi-C (high-throughput chromosome conformation capture), and genome-wide association study data implicated the NR2F1 gene and associated long noncoding RNAs within the 5q15 locus. Heritability partitioning identified statistically significant enrichment in brain tissue histone marks, while polygenic risk scoring of brain volume data identified statistically significant associations with right and left thalamus volumes and right putamen volume.

CONCLUSIONS

Our work presents novel insights into the neurobiology of TS, thereby opening up new directions for future studies.

Vitamin D affects the risk of disease activity in multiple sclerosis

BACKGROUND

Vitamin D (VitD) affects the risk of multiple sclerosis (MS), but the impact on disease activity is controversial. We assessed whether VitD is associated with the No-Evidence of Disease Activity-3 (NEDA-3) status at 2 years from disease-modifying treatment (DMT) start, and whether this association is causal or the result of confounding factors. Furthermore, we explored if a genetic predisposition to higher VitD levels affects the risk of disease activity.

METHODS

230 untreated relapsing-remitting MS patients underwent serum 25-OH-vitamin-D measurement, and the association between seasonally adjusted VitD and disease activity was tested. Modelling a Polygenic Risk Score from a Genome-Wide Association Study on ~400 000 individuals, we studied the impact of genetic predisposition to higher VitD on the NEDA-3 status in 1408 independent MS patients. Two-sample Mendelian randomisation (MR) was used to assess causality.

RESULTS

Lower baseline VitD was associated with decreased probability of NEDA-3 at 2 years (p=0.019). Particularly, VitD levels <20 ng/mL conferred an over twofold risk of disease activity (OR 2.36, 95% CI 1.30 to 3.88, p=0.0037). Genetic predisposition to higher VitD levels was associated with delayed age at MS onset (p=0.018) and with a higher probability of NEDA-3 status (p=0.034). MR confirmed causality between VitD and the risk of disease activity (p=0.041).

CONCLUSIONS

VitD levels before DMT start affect the risk of disease activity in MS. Genetic predisposition to higher VitD levels confers a lower risk of disease activity and is associated with delayed MS onset. Our work prompts future prospective studies regarding VitD supplementation and lifestyle interventions to hamper disease activity in MS.

Genetic Variations in TrkB.T1 Isoform and Their Association With Somatic and Psychological Symptoms in Individuals With IBS

Irritable bowel syndrome (IBS), a disorder of gut-brain interaction, is often comorbid with somatic pain and psychological disorders. Dysregulated signaling of brain-derived neurotrophic factor (BDNF) and its receptor, tropomyosin-related kinase B (TrkB), has been implicated in somatic-psychological symptoms in individuals with IBS. We investigated the association of 10 single-nucleotide polymorphisms (SNPs) in the regulatory 3' untranslated region of neurotrophic receptor tyrosine kinase-2 (NTRK2) kinase domain-deficient truncated isoform (TrkB.T1) and BDNF Val66Met SNP with somatic and psychological symptoms and quality-of-life (QoL) in a cohort from the United States (IBS, n = 464; healthy controls, n = 156). We found that the homozygous recessive genotype (G/G) of rs2013566 in individuals with IBS is associated with worsened somatic symptoms, including headache, back pain, joint pain, muscle pain, and somatization as well as diminished sleep quality, energy level, and overall QoL. Validation using United Kingdom BioBank data confirmed the association of rs2013566 with an increased likelihood of headache. Several SNPs (rs1627784, rs1624327, and rs1147198) showed significant associations with muscle pain in our U.S. cohort. These 4 SNPs are predominantly located in H3K4Me1-enriched regions, suggesting their enhancer and/or transcription regulation potential. Our findings suggest that genetic variation within the 3' untranslated region region of the TrkB.T1 isoform may contribute to comorbid conditions in individuals with IBS, resulting in a spectrum of somatic and psychological symptoms impacting their QoL. These findings advance our understanding of the genetic interaction between BDNF/TrkB pathways and somatic-psychological symptoms in IBS, highlighting the importance of further exploring this interaction for potential clinical applications. PERSPECTIVE: This study aims to understand the genetic effects on IBS-related symptoms across somatic, psychological, and quality-of-life (QoL) domains, validated by United Kingdom BioBank data. The rs2013566 homozygous recessive genotype correlates with worsened somatic symptoms and reduced QoL, emphasizing its clinical significance.

Ancient genomic analysis of a Chinese hereditary elite from the Northern and Southern Dynasties

China's Northern and Southern Dynasties period (3rd-6th centuries AD) marked a significant era of ethnic integration in northern China. However, previous ancient DNA studies have primarily focused on northern ethnic groups, with limited research on the genetic formation of the hereditary elite family, especially considering their abundant archaeological record and clear material identity. In this study, we obtained the ancient genome of a hereditary elite family, Gao Bin (, 503-572 AD), at 0.6473-fold coverage with 475132 single-nucleotide polymorphisms (SNPs) on the 1240k panel. His mitochondrial haplogroup belonged to Z4 and Y-haplogroup to O1a1a2b-F2444*. The genetic profile of Gao Bin was most similar to that of the northern Han Chinese. He could be modelled as deriving all his ancestry from Late Neolithic to Iron Age Yellow River farmers without influence from Northeast Asia, Korea, or the Mongolian Plateau. Our study sheds light on the genetic formation of hereditary elite families in the context of the Southern and Northern Dynasties ethnic integration.

Ancient genomes illuminate the demographic history of Shandong over the past two millennia

Shandong province, located in the Lower Yellow River, is one of the birthplaces of ancient Chinese civilization. However, the comprehensive genetic histories of this region have remained largely unknown until now due to a lack of ancient human genomes. Here, we present 21 ancient genomes from Shandong dating from the Warring States period to the Jin-Yuan Dynasties. Unlike the early Neolithic samples from Shandong, the historical samples are most closely related to post-Late Neolithic populations of the Middle Yellow River Basin, suggesting a population turnover in Shandong from the Neolithic Age to the Historical era. In addition, we detect a close genetic affinity between the historical samples in Shandong and present-day Han Chinese, showing long-term genetic stability in Han Chinese at least since the Warring States period.

Natural selection drives emergent genetic homogeneity in a century-scale experiment with barley

Direct observation is central to our understanding of adaptation, but evolution is rarely documented in a large, multicellular organism for more than a few generations. In this study, we observed evolution across a century-scale competition experiment, barley composite cross II (CCII). CCII was founded in 1929 in Davis, California, with thousands of genotypes, but we found that natural selection has massively reduced genetic diversity, leading to a single lineage constituting most of the population by generation 50. Selection favored alleles originating from climates similar to that of Davis and targeted loci contributing to reproductive development, including the barley diversification loci Vrs1, HvCEN, Ppd-H1, and Vrn-H2. Our findings point to selection as the predominant force shaping genomic variation in one of the world's oldest biological experiments.

Genetic and multi-omic resources for Alzheimer disease and related dementia from the Knight Alzheimer Disease Research Center

The Knight-Alzheimer Disease Research Center (Knight-ADRC) at Washington University in St. Louis has pioneered and led worldwide seminal studies that have expanded our clinical, social, pathological, and molecular understanding of Alzheimer Disease. Over more than 40 years, research volunteers have been recruited to participate in cognitive, neuropsychologic, imaging, fluid biomarkers, genomic and multi-omic studies. Tissue and longitudinal data collected to foster, facilitate, and support research on dementia and aging. The Genetics and high throughput -omics core (GHTO) have collected of more than 26,000 biological samples from 6,625 Knight-ADRC participants. Samples available include longitudinal DNA, RNA, non-fasted plasma, cerebrospinal fluid pellets, and peripheral blood mononuclear cells. The GHTO has performed deep molecular profiling (genomic, transcriptomic, epigenomic, proteomic, and metabolomic) from large number of brain (n = 2,117), CSF (n = 2,012) and blood/plasma (n = 8,265) samples with the goal of identifying novel risk and protective variants, identify novel molecular biomarkers and causal and druggable targets. Overall, the resources available at GHTO support the increase of our understanding of Alzheimer Disease.

Genomic variation of 363 diverse tea accessions unveils the genetic diversity, domestication, and structural variations associated with tea adaptation

Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

Evaluation of polygenic scoring methods in five biobanks shows larger variation between biobanks than methods and finds benefits of ensemble learning

Methods of estimating polygenic scores (PGSs) from genome-wide association studies are increasingly utilized. However, independent method evaluation is lacking, and method comparisons are often limited. Here, we evaluate polygenic scores derived via seven methods in five biobank studies (totaling about 1.2 million participants) across 16 diseases and quantitative traits, building on a reference-standardized framework. We conducted meta-analyses to quantify the effects of method choice, hyperparameter tuning, method ensembling, and the target biobank on PGS performance. We found that no single method consistently outperformed all others. PGS effect sizes were more variable between biobanks than between methods within biobanks when methods were well tuned. Differences between methods were largest for the two investigated autoimmune diseases, seropositive rheumatoid arthritis and type 1 diabetes. For most methods, cross-validation was more reliable for tuning hyperparameters than automatic tuning (without the use of target data). For a given target phenotype, elastic net models combining PGS across methods (ensemble PGS) tuned in the UK Biobank provided consistent, high, and cross-biobank transferable performance, increasing PGS effect sizes (β coefficients) by a median of 5.0% relative to LDpred2 and MegaPRS (the two best-performing single methods when tuned with cross-validation). Our interactively browsable online-results and open-source workflow prspipe provide a rich resource and reference for the analysis of polygenic scoring methods across biobanks.

A multi-ancestry GWAS meta-analysis of facial features and its application in predicting archaic human features

Facial morphology, a complex trait influenced by genetics, holds great significance in evolutionary research. However, due to limited fossil evidence, the facial characteristics of Neanderthals and Denisovans have remained largely unknown. In this study, we conducted a large-scale multi-ethnic meta-analysis of the genome-wide association study (GWAS), including 9674 East Asians and 10,115 Europeans, quantitatively assessing 78 facial traits using 3D facial images. We identified 71 genomic loci associated with facial features, including 21 novel loci. We developed a facial polygenic score (FPS) that enables the prediction of facial features based on genetic information. Interestingly, the distribution of FPSs among populations from diverse continental groups exhibited relevant correlations with observed facial features. Furthermore, we applied the FPS to predict the facial traits of seven Neanderthals and one Denisovan using ancient DNA and aligned predictions with the fossil records. Our results suggested that Neanderthals and Denisovans likely shared similar facial features, such as a wider but shorter nose and a wider endocanthion distance. The decreased mouth width was characterized specifically in Denisovans. The integration of genomic data and facial trait analysis provides valuable insights into the evolutionary history and adaptive changes in human facial morphology.

Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed

BACKGROUND

The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity.

RESULTS

To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders.

CONCLUSIONS

The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.

Patterns of reproduction and autozygosity distinguish the breeding from nonbreeding gray wolves of Yellowstone National Park

For species of management concern, accurate estimates of inbreeding and associated consequences on reproduction are crucial for predicting their future viability. However, few studies have partitioned this aspect of genetic viability with respect to reproduction in a group-living social mammal. We investigated the contributions of foundation stock lineages, putative fitness consequences of inbreeding, and genetic diversity of the breeding versus nonreproductive segment of the Yellowstone National Park gray wolf population. Our dataset spans 25 years and seven generations since reintroduction, encompassing 152 nuclear families and 329 litters. We found more than 87% of the pedigree foundation genomes persisted and report influxes of allelic diversity from two translocated wolves from a divergent source in Montana. As expected for group-living species, mean kinship significantly increased over time but with minimal loss of observed heterozygosity. Strikingly, the reproductive portion of the population carried a significantly lower genome-wide inbreeding coefficients, autozygosity, and more rapid decay for linkage disequilibrium relative to the nonbreeding population. Breeding wolves had significantly longer lifespans and lower inbreeding coefficients than nonbreeding wolves. Our model revealed that the number of litters was negatively significantly associated with heterozygosity (R = -0.11). Our findings highlight genetic contributions to fitness, and the importance of the reproductively active individuals in a population to counteract loss of genetic variation in a wild, free-ranging social carnivore. It is crucial for managers to mitigate factors that significantly reduce effective population size and genetic connectivity, which supports the dispersion of genetic variation that aids in rapid evolutionary responses to environmental challenges.

The value of hybrid genomes: Building two highly contiguous reference genome assemblies to advance Canis genomic studies

Previous studies of canid population and evolutionary genetics have relied on high-quality domestic dog reference genomes that have been produced primarily for biomedical and trait mapping studies in dog breeds. However, the absence of highly contiguous genomes from other Canis species like the gray wolf and coyote, that represent additional distinct demographic histories, may bias inferences regarding interspecific genetic diversity and phylogenetic relationships. Here, we present single haplotype de novo genome assemblies for the gray wolf and coyote, generated by applying the trio-binning approach to long sequence reads generated from the genome of a female first-generation hybrid produced from a gray wolf and coyote mating. The assemblies were highly contiguous, with contig N50 sizes of 44.6 and 42.0 Mb for the wolf and coyote, respectively. Genome scaffolding and alignments between the two Canis assemblies and published dog reference genomes showed near complete collinearity, with one exception: a coyote-specific chromosome fission of chromosome 13 and fusion of the proximal portion of that chromosome with chromosome 8, retaining the Canis-typical haploid chromosome number of 2n = 78. We evaluated mapping quality for previous RADseq data from 334 canids and found nearly identical mapping quality and patterns among canid species and regional populations regardless of the genome used for alignment (dog, coyote, or gray wolf). These novel wolf and coyote genome reference assemblies will be important resources for proper and accurate inference of Canis demography, taxonomic evaluation, and conservation genetics.

eQTLs identify regulatory networks and drivers of variation in the individual response to sepsis

Sepsis is a clinical syndrome of life-threatening organ dysfunction caused by a dysregulated response to infection, for which disease heterogeneity is a major obstacle to developing targeted treatments. We have previously identified gene-expression-based patient subgroups (sepsis response signatures [SRS]) informative for outcome and underlying pathophysiology. Here, we aimed to investigate the role of genetic variation in determining the host transcriptomic response and to delineate regulatory networks underlying SRS. Using genotyping and RNA-sequencing data on 638 adult sepsis patients, we report 16,049 independent expression (eQTLs) and 32 co-expression module (modQTLs) quantitative trait loci in this disease context. We identified significant interactions between SRS and genotype for 1,578 SNP-gene pairs and combined transcription factor (TF) binding site information (SNP2TFBS) and predicted regulon activity (DoRothEA) to identify candidate upstream regulators. Overall, these approaches identified putative mechanistic links between host genetic variation, cell subtypes, and the individual transcriptomic response to infection.

Quantitative trait locus mapping in placenta: A comparative study of chorionic villus and birth placenta

The placenta, a pivotal player in the prenatal environment, holds crucial insights into early developmental pathways and future health outcomes. In this study, we explored genetic molecular regulation in chorionic villus samples (CVS) from the first trimester and placenta tissue at birth. We assessed quantitative trait locus (QTL) mapping on DNA methylation and gene expression data in a Finnish cohort of 574 individuals. We found more QTLs in birth placenta than in first-trimester placenta. Nevertheless, a substantial amount of associations overlapped in their effects and showed consistent direction in both tissues, with increasing molecular genetic effects from early pregnancy to birth placenta. The identified QTLs in birth placenta were most enriched in genes with placenta-specific expression. Conducting a phenome-wide-association study (PheWAS) on the associated SNPs, we observed numerous overlaps with genome-wide association study (GWAS) hits (spanning 57 distinct traits and 23 SNPs), with notable enrichments for immunological, skeletal, and respiratory traits. The QTL-SNP rs1737028 (chr6:29737993) presented with the highest number of GWAS hits. This SNP was related to HLA-G expression via DNA methylation and was associated with various immune, respiratory, and psychiatric traits. Our findings implicate increasing genetic molecular regulation during the course of pregnancy and support the involvement of placenta gene regulation, particularly in immunological traits. This study presents a framework for understanding placenta-specific gene regulation during pregnancy and its connection to health-related traits.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

HLA and KIR genetic association and NK cells in anti-NMDAR encephalitis

Introduction

Genetic predisposition to autoimmune encephalitis with antibodies against N-methyl-D-aspartate receptor (NMDAR) is poorly understood. Given the diversity of associated environmental factors (tumors, infections), we hypothesized that human leukocyte antigen (HLA) and killer-cell immunoglobulin-like receptors (KIR), two extremely polymorphic gene complexes key to the immune system, might be relevant for the genetic predisposition to anti-NMDAR encephalitis. Notably, KIR are chiefly expressed by Natural Killer (NK) cells, recognize distinct HLA class I allotypes and play a major role in anti-tumor and anti-infection responses.

Methods

We conducted a Genome Wide Association Study (GWAS) with subsequent control-matching using Principal Component Analysis (PCA) and HLA imputation, in a multi-ethnic cohort of anti-NMDAR encephalitis (n=479); KIR and HLA were further sequenced in a large subsample (n=323). PCA-controlled logistic regression was then conducted for carrier frequencies (HLA and KIR) and copy number variation (KIR). HLA-KIR interaction associations were also modeled. Additionally, single cell sequencing was conducted in peripheral blood mononuclear cells from 16 cases and 16 controls, NK cells were sorted and phenotyped.

Results

Anti-NMDAR encephalitis showed a weak HLA association with DRB1*01:01~DQA1*01:01~DQB1*05:01 (OR=1.57, 1.51, 1.45; respectively), and DRB1*11:01 (OR=1.60); these effects were stronger in European descendants and in patients without an underlying ovarian teratoma. More interestingly, we found increased copy number variation of KIR2DL5B (OR=1.72), principally due to an overrepresentation of KIR2DL5B*00201. Further, we identified two allele associations in framework genes, KIR2DL4*00103 (25.4% vs. 12.5% in controls, OR=1.98) and KIR3DL3*00302 (5.3% vs. 1.3%, OR=4.44). Notably, the ligands of these KIR2DL4 and KIR3DL3, respectively, HLA-G and HHLA2, are known to act as immune checkpoint with immunosuppressive functions. However, we did not find differences in specific KIR-HLA ligand interactions or HLA-G polymorphisms between cases and controls. Similarly, gene expression of CD56dim or CD56bright NK cells did not differ between cases and controls.

Discussion

Our observations for the first time suggest that the HLA-KIR axis might be involved in anti-NMDAR encephalitis. While the genetic risk conferred by the identified polymorphisms appears small, a role of this axis in the pathophysiology of this disease appears highly plausible and should be analyzed in future studies.

Characterising the genetic architecture of changes in adiposity during adulthood using electronic health records

Obesity is a heritable disease, characterised by excess adiposity that is measured by body mass index (BMI). While over 1,000 genetic loci are associated with BMI, less is known about the genetic contribution to adiposity trajectories over adulthood. We derive adiposity-change phenotypes from 24.5 million primary-care health records in over 740,000 individuals in the UK Biobank, Million Veteran Program USA, and Estonian Biobank, to discover and validate the genetic architecture of adiposity trajectories. Using multiple BMI measurements over time increases power to identify genetic factors affecting baseline BMI by 14%. In the largest reported genome-wide study of adiposity-change in adulthood, we identify novel associations with BMI-change at six independent loci, including rs429358 (APOE missense variant). The SNP-based heritability of BMI-change (1.98%) is 9-fold lower than that of BMI. The modest genetic correlation between BMI-change and BMI (45.2%) indicates that genetic studies of longitudinal trajectories could uncover novel biology of quantitative traits in adulthood.

Exome-wide evidence of compound heterozygous effects across common phenotypes in the UK Biobank

The phenotypic impact of compound heterozygous (CH) variation has not been investigated at the population scale. We phased rare variants (MAF ∼0.001%) in the UK Biobank (UKBB) exome-sequencing data to characterize recessive effects in 175,587 individuals across 311 common diseases. A total of 6.5% of individuals carry putatively damaging CH variants, 90% of which are only identifiable upon phasing rare variants (MAF < 0.38%). We identify six recessive gene-trait associations (p < 1.68 × 10-7) after accounting for relatedness, polygenicity, nearby common variants, and rare variant burden. Of these, just one is discovered when considering homozygosity alone. Using longitudinal health records, we additionally identify and replicate a novel association between bi-allelic variation in ATP2C2 and an earlier age at onset of chronic obstructive pulmonary disease (COPD) (p < 3.58 × 10-8). Genetic phase contributes to disease risk for gene-trait pairs: ATP2C2-COPD (p = 0.000238), FLG-asthma (p = 0.00205), and USH2A-visual impairment (p = 0.0084). We demonstrate the power of phasing large-scale genetic cohorts to discover phenome-wide consequences of compound heterozygosity.

Ethnic-specific genetic susceptibility loci for endometriosis in Taiwanese-Han population: a genome-wide association study

Endometriosis is a common gynecological disorder affecting around 10% of reproductive-age women. Although many hypotheses were proposed, genetic alteration has been considered as one of the key factors promoting pathogenesis. Due to racial/ethnic disparities in the process of hormone regulation and nutrition metabolism, a genome-wide association study (GWAS) with 2794 cases and 27,940 controls was conducted in a Taiwanese-Han population. Our study identified five significant susceptibility loci for endometriosis, and three of them, WNT4 (on the 1p36.12), RMND1 (6q25.1), and CCDC170 (6q25.1), have been previously associated with endometriosis across different populations, including European and Japanese descent cohorts. Other two including C5orf66/C5orf66-AS2 (5q31.1) and STN1 (10q24.33) are newly identified ones. Functional network analysis of potent risk genes revealed the involvement of cancer susceptibility and neurodevelopmental disorders in endometriosis development. In addition, long non-coding RNAs (lncRNAs) C5orf66 and C5orf66-AS2 can interact with many RNA-binding proteins (RBPs) which can influence RNA metabolic process, mRNA stabilization, and mRNA splicing, leading to dysregulation in tumor-promoting gene expression. Those findings support clinical observations of differences in the presentation of endometriosis in Taiwanese-Han population with higher risks of developing deeply infiltrating/invasive lesions and the associated malignancies.

Tall women with breast cancer have poorer survival than short women

Background

Tall women are more likely to develop breast cancer (BC). High Mobility Group AT-Hook 1(HMGA1), an oncofetal protein, plays a role in the progression of breast cancer. Non-coding sequences proximal to HMGA1 contain variants associated with 4.83 cm taller height. In the current study, we used UK Biobank data to examine the relationship of HMGA1 to height, risk, and prognosis of women with breast cancer.

Methods

Our analysis included all subjects with invasive BC that occurred either before or after participant enrollment and were recorded in the UK Biobank database using self-reported data and the International Classification of Diseases (ICD10, ICD9). We divided the subjects into three previously described three height groups: Short (< 155 cm), Medium (155 cm to 175 cm), Tall (> 175 cm). We analyzed the HMGA1 SNP rs41269028, a single nucleotide intron variant, C > T, minor allele frequency 0.044. SNP rs41269028 was previously evaluated in subjects with diabetes.

Results

Height of 9583 women with BC homozygous for the HMGA1 SNP rs41269028 major allele was 162.29 cm ± 6.18. Height of 944 women with BC who were carriers or homozygotes (CT + TT) of the minor allele T was 162.88 cm ± 6.001. This difference was significant (p = 0.005). The effect of height group on survival was significant (p = 0.032, log rank test). Tall women had the poorest survival. The effect of HMGA1 SNP rs41269028 genotype on BC risk (p = 0.602) and survival (p = 0.439, log rank test) was insignificant.

Conclusion

We conclude that HMGA1 influences height, but we were unable to demonstrate that HMGA1 is related to increased incidence or poor prognosis of tall women with breast cancer. We did find that tall women with breast cancer have poorer survival than short women. Our finding that tall women have a worse prognosis is important because it could help the oncologist decide, along with other prognostic factors, whether adjuvant therapy is warranted.

Characterization and distribution of a 14-Mb chromosomal inversion in native populations of rainbow trout (Oncorhynchus mykiss)

Multiple studies in a range of taxa have found links between structural variants and the development of ecologically important traits. Such variants are becoming easier to find due, in large part, to the increase in the amount of genome-wide sequence data in nonmodel organisms. The salmonids (salmon, trout, and charr) are a taxonomic group with abundant genome-wide datasets due to their importance in aquaculture, fisheries, and variation in multiple ecologically important life-history traits. Previous research on rainbow trout (Oncorhynchus mykiss) has documented a large pericentric (∼55 Mb) chromosomal inversion (CI) on chromosome 5 (Omy05) and a second smaller (∼14 Mb) chromosome inversion on Omy20. While the Omy05 inversion appears to be associated with multiple adaptive traits, the inversion on Omy20 has received far less attention. In this study, we re-analyze RAD-seq and amplicon data from several populations of rainbow trout (O. mykiss) to better document the structure and geographic distribution of variation in the Omy20 CI. Moreover, we utilize phylogenomic techniques to characterize both the age- and the protein-coding gene content of the Omy20 CI. We find that the age of the Omy20 inversion dates to the early stages of O. mykiss speciation and predates the Omy05 inversion by ∼450,000 years. The 2 CIs differ further in terms of the frequency of the homokaryotypes. While both forms of the Omy05 CI are found across the eastern Pacific, the ancestral version of the Omy20 CI is restricted to the southern portion of the species range in California. Furthermore, the Omy20 inverted haplotype is comparable in genetic diversity to the ancestral form, whereas derived CIs typically show substantially reduced genetic diversity. These data contribute to our understanding of the age and distribution of a large CI in rainbow trout and provide a framework for researchers looking to document CIs in other nonmodel species.

Increased frequency of repeat expansion mutations across different populations

Repeat expansion disorders (REDs) are a devastating group of predominantly neurological diseases. Together they are common, affecting 1 in 3,000 people worldwide with population-specific differences. However, prevalence estimates of REDs are hampered by heterogeneous clinical presentation, variable geographic distributions, and technological limitations leading to under-ascertainment. Here, leveraging whole genome sequencing data from 82,176 individuals from different populations, we found an overall disease allele frequency of REDs of 1 in 283 individuals. Modelling disease prevalence using genetic data, age at onset and survival, we show that the expected number of people with REDs would be two to three times higher than currently reported figures, indicating under-diagnosis and/or incomplete penetrance. While some REDs are population-specific, e.g. Huntington disease-like 2 in Africans, most REDs are represented in all broad genetic ancestries (i.e. Europeans, Africans, Americans, East Asians, and South Asians), challenging the notion that some REDs are found only in specific populations. These results have worldwide implications for local and global health communities in the diagnosis and counselling of REDs.

How group structure impacts the numbers at risk for coronary artery disease: polygenic risk scores and nongenetic risk factors in the UK Biobank cohort

The UK Biobank (UKB) is a large cohort study that recruited over 500,000 British participants aged 40-69 in 2006-2010 at 22 assessment centers from across the United Kingdom. Self-reported health outcomes and hospital admission data are 2 types of records that include participants' disease status. Coronary artery disease (CAD) is the most common cause of death in the UKB cohort. After distinguishing between prevalence and incidence CAD events for all UKB participants, we identified geographical variations in age-standardized rates of CAD between assessment centers. Significant distributional differences were found between the pooled cohort equation scores of UKB participants from England and Scotland using the Mann-Whitney test. Polygenic risk scores of UKB participants from England and Scotland and from different assessment centers differed significantly using permutation tests. Our aim was to discriminate between assessment centers with different disease rates by collecting data on disease-related risk factors. However, relying solely on individual-level predictions and averaging them to obtain group-level predictions proved ineffective, particularly due to the presence of correlated covariates resulting from participation bias. By using the Mundlak model, which estimates a random effects regression by including the group means of the independent variables in the model, we effectively addressed these issues. In addition, we designed a simulation experiment to demonstrate the functionality of the Mundlak model. Our findings have applications in public health funding and strategy, as our approach can be used to predict case rates in the future, as both population structure and lifestyle changes are uncertain.

An assessment of the genomic structural variation landscape in Sub-Saharan African populations

Structural variants are responsible for a large part of genomic variation between individuals and play a role in both common and rare diseases. Databases cataloguing structural variants notably do not represent the full spectrum of global diversity, particularly missing information from most African populations. To address this representation gap, we analysed 1,091 high-coverage African genomes, 545 of which are public data sets, and 546 which have been analysed for structural variants for the first time. Variants were called using five different tools and datasets merged and jointly called using SURVIVOR. We identified 67,795 structural variants throughout the genome, with 10,421 genes having at least one variant. Using a conservative overlap in merged data, 6,414 of the structural variants (9.5%) are novel compared to the Database of Genomic Variants. This study contributes to knowledge of the landscape of structural variant diversity in Africa and presents a reliable dataset for potential applications in population genetics and health-related research.

Museum genomics approach to study the taxonomy and evolution of Woolly-necked storks using historic specimens

The accessibility of genomic tools in evolutionary biology has allowed for a thorough exploration of various evolutionary processes associated with adaptation and speciation. However, genomic studies in natural systems present numerous challenges, reflecting the inherent complexities of studying organisms in their native habitats. The utilization of museum specimens for genomics research has received increased attention in recent times, facilitated by advancements in ancient DNA techniques. In this study, we have utilized a museum genomics approach to analyze historic specimens of Woolly-necked storks (Ciconia spp.) and examine their genetic composition and taxonomic status and explore the evolutionary and adaptive trajectories of populations over the years. The Woolly-necked storks are distributed in Asia and Africa with a taxonomic classification that has been a matter of ambiguity. Asian and African Woollynecks were recently recognized as different species based on their morphological differences; however, their genomic validation was lacking. In this study, we have used ∼70-year-old museum samples for whole-genome population-scale sequencing. Our study has revealed that Asian and African Woollynecks are genetically distinct, consistent with the current taxonomic classification based on morphological features. However, we also found a high genetic divergence between the Asian subspecies Ciconia episcopus neglecta and Ciconia episcopus episcopus, suggesting this classification requires a detailed examination to explore processes of ongoing speciation. Because taxonomic classification directly impacts conservation efforts, and there is evidence of declining populations of Asian Woollynecks in Southeast Asia, our results highlight that population-scale studies are urgent to determine the genetic, ecological, and phylogenetic diversity of these birds.

A Genome-Wide Association Study Suggests New Susceptibility Loci for Primary Antiphospholipid Syndrome

OBJECTIVE

Primary antiphospholipid syndrome (PAPS) is a rare autoimmune disease characterized by the presence of antiphospholipid antibodies and the occurrence of thrombotic events and pregnancy complications. Our study aimed to identify novel genetic susceptibility loci associated with PAPS.

METHODS

We performed a genome-wide association study comprising 5,485 individuals (482 affected individuals) of European ancestry. Significant and suggestive independent variants from a meta-analysis of approximately 7 million variants were evaluated for functional and biological process enrichment. The genetic risk variability for PAPS in different populations was also assessed. Hierarchical clustering, Mahalanobis distance, and Dirichlet Process Mixtures with uncertainty clustering methods were used to assess genetic similarities between PAPS and other immune-mediated diseases.

RESULTS

We revealed genetic associations with PAPS in a regulatory locus within the HLA class II region near HLA-DRA and in STAT1-STAT4 with a genome-wide level of significance; 34 additional suggestive genetic susceptibility loci for PAPS were also identified. The disease risk allele near HLA-DRA is associated with overexpression of HLA-DRB6, HLA-DRB9, HLA-DQA2, and HLA-DQB2 in immune cells, vascular tissue, and nervous tissue. This association is independent of the association between PAPS and HLA-DRB1*1302. Functional analyses highlighted immune-related pathways in PAPS-associated loci. The comparison with other immune-mediated diseases revealed a close genetic relatedness to neuromyelitis optica, systemic sclerosis, and Sjögren syndrome, suggesting co-localized causal variations close to STAT1-STAT4, TNPO3, and BLK.

CONCLUSION

This study represents a comprehensive large-scale genetic analysis for PAPS and provides new insights into the genetic basis and pathophysiology of this rare disease.

Genome scans reveal signals of selection associated with pollution in fish populations of Basilichthys microlepidotus, an endemic species of Chile

The Maipo River catchment is one of Chile's most polluted basins. In recent decades, discharges of untreated sewage and organic matter have caused eutrophication and water quality degradation. We employed the indigenous silverfish species Basilichthys microlepidotus as a model organism to investigate the process of adaptation and selection on genes influenced by pollution. Using variation at single nucleotide polymorphisms (SNPs), we determined the temporal stability of the population structure patterns previously identified in this species by varying SNPs. We also examined local adaptation to pollution-selected genes. Using the genotypes of 7684 loci in 180 individuals, we identified 429 and 700 loci that may be undergoing selection. We detected these loci using the FSTHET and ARLEQUIN outlier detection software, respectively. Both software packages simultaneously identified a total of 250 loci. B. microlepidotus' population structure did not change over time at contaminated or unpolluted sites. In addition, our analysis found: (i) selection of genes associated with pollution, consistent with observations in other organisms; (ii) identification of candidate genes that are functionally linked to the same biological processes, molecular functions and/or cellular components that previously showed differential expression in the same populations; and (iii) a candidate gene with differential expression and a non-synonymous substitution.

Genetic Associations of Anhedonia: Insights into Overlap of Mental and Somatic Disorders

BACKGROUND

Anhedonia is characterized by a reduced ability to anticipate, experience, and/or learn about pleasure. This phenomenon has a transdiagnostic nature and is one of the key symptoms of mood disorders, schizophrenia, addictions, and somatic conditions.

AIM

To evaluate the genetic architecture of anhedonia and its overlap with other mental disorders and somatic conditions.

METHODS

We performed a genome-wide association study of anhedonia on a sample of 4,520 individuals from a Russian non-clinical population. Using the available summary statistics, we calculated polygenic risk scores (PRS) to investigate the genetic relationship between anhedonia and other psychiatric or somatic phenotypes.

RESULTS

No variants with a genome-wide significant association were identified. PRS for major depression, bipolar disorder, and schizophrenia were significantly associated with anhedonia. Conversely, no significant associations were found between PRS for anxiety and anhedonia, which aligns well with existing clinical evidence. None of the PRS for somatic phenotypes attained a significance level after correction for multiple comparisons. A nominal significance for the anhedonia association was determined for omega-3 fatty acids, type 2 diabetes mellitus, and Crohn's disease.

CONCLUSION

Anhedonia has a complex polygenic architecture, and its presence in somatic diseases or normal conditions may be due to a genetic predisposition to mood disorders or schizophrenia.

Drawing human pedigree charts with DrawPed

Creating pedigree charts is a recurring task in biomedical research, but there are few online tools for drawing complex human pedigrees available and even fewer are free. With DrawPed we aim to close this gap. DrawPed automatically draws pedigree charts from standard PED format pedigree files. Users can also create pedigrees from scratch and interactively edit existing pedigrees. The application can display conditions not captured in a PED file such as deceased persons or suspected consanguinity of parents. Pedigree charts are displayed as SVGs, which are scalable and hence publication-ready. Pedigrees can be exported as PED files for storage, exchange, or use in other applications. DrawPed is open source and freely available at https://www.genecascade.org/DrawPed/.

Sequence variants associated with resilient responses in growing pigs

The current work aimed to identify genomic regions and candidate genes associated with resilience in pigs. In previous work, we proposed the body weight deviation from the expected growth curve (ΔBW) and the increase of the positive acute-phase protein haptoglobin (ΔHP) after a vaccine challenge as resilience indicators which may be improved through selective breeding in pigs. Individuals with steady growth rate and minor activation of haptoglobin (high ΔBW and low ΔHP values) were considered resilient. In contrast, pigs with perturbed growth rate and high activation of haptoglobin (low ΔBW and high ΔHP values) were considered susceptible. Both ∆BW and ∆HP were simultaneously considered to select the most resilient (N = 40) and susceptible (N = 40) pigs. A genome-wide association study was carried out for the pigs' response classification to the challenge test using whole-genome sequence data (7,760,720 variants). Eleven associated genomic regions were identified, harbouring relevant candidate genes related to the immune response (such as pro- and anti-inflammatory responses) and growth pathways. These associated genomic regions harboured 41 potential functional mutations (frameshift, splice donor, splice acceptor, start loss and stop loss/gain) in candidate genes. Overall, this study advances our knowledge about the genetic determinism of resilience, highlighting its polygenic nature and strong relationship with immunity and growth.

Imputation Server PGS: an automated approach to calculate polygenic risk scores on imputation servers

Polygenic scores (PGS) enable the prediction of genetic predisposition for a wide range of traits and diseases by calculating the weighted sum of allele dosages for genetic variants associated with the trait or disease in question. Present approaches for calculating PGS from genotypes are often inefficient and labor-intensive, limiting transferability into clinical applications. Here, we present 'Imputation Server PGS', an extension of the Michigan Imputation Server designed to automate a standardized calculation of polygenic scores based on imputed genotypes. This extends the widely used Michigan Imputation Server with new functionality, bringing the simplicity and efficiency of modern imputation to the PGS field. The service currently supports over 4489 published polygenic scores from publicly available repositories and provides extensive quality control, including ancestry estimation to report population stratification. An interactive report empowers users to screen and compare thousands of scores in a fast and intuitive way. Imputation Server PGS provides a user-friendly web service, facilitating the application of polygenic scores to a wide range of genetic studies and is freely available at https://imputationserver.sph.umich.edu.

Genetic and genomic analysis of Belgian Blue's susceptibility for psoroptic mange

BACKGROUND

Psoroptic mange, caused by Psoroptes ovis mites, is affecting Belgian Blue cattle's welfare and production potential. The Belgian Blue cattle-known for its high degree of muscling, low feed conversion ratio and high beef quality-is highly susceptible for this disease.

RESULTS

In this study, we phenotyped 1975 Belgian Blue cattle from more than 100 different groups on commercial beef farms for their psoroptic mange susceptibility. Substantial individual differences were observed within these management groups, with lesion extent differences up to ± 15%. Animal models showed that estimated heritabilities were low for lesion extent and severe lesion extent (0.07 and 0.09, respectively) and 0.12 for the number of mites. A genome wide association study for mange susceptibility revealed signals on BTA6, BTA11, BTA15 and BTA24. In these regions, candidate genes GBA3, RAG2, and TRAF6 were identified.

CONCLUSIONS

Despite the challenges in phenotyping for psoroptic mange due to the timing of screening, the continuous evolution of lesions and different management conditions, we successfully conducted a study on the genetic susceptibility to psoroptic mange in Belgian Blue cattle. Our results clearly indicate that psoroptic mange is under polygenic control and the underlying candidate genes should be studied more thoroughly. This is the first study providing candidate genes for this complex disease. These results are already valuable for Belgian Blue breeding, however, further research is needed to unravel the architecture of this disease and to identify causal mutations.

Congenital syndromic Chiari-like malformation (CSCM) in Holstein cattle: towards unravelling of possible genetic causes

BACKGROUND

Chiari malformation type II (CMII) was originally reported in humans as a rare disorder characterized by the downward herniation of the hindbrain and towering cerebellum. The congenital brain malformation is usually accompanied by spina bifida, a congenital spinal anomaly resulting from incomplete closure of the dorsal aspect of the spinal neural tube, and occasionally by other lesions. A similar disorder has been reported in several animal species, including cattle, particularly as a congenital syndrome. A cause of congenital syndromic Chiari-like malformation (CSCM) in cattle has not been reported to date. We collected a series of 14 CSCM-affected Holstein calves (13 purebred, one Red Danish Dairy F1 cross) and performed whole-genome sequencing (WGS). WGS was performed on 33 cattle, including eight cases with parents (trio-based; group 1), three cases with one parent (group 2), and three single cases (solo-based; group 3).

RESULTS

Sequencing-based genome-wide association study of the 13 Holstein calves with CSCM and 166 controls revealed no significantly associated genome region. Assuming a single Holstein breed-specific recessive allele, no region of shared homozygosity was detected suggesting heterogeneity. Subsequent filtering for protein-changing variants that were only homozygous in the genomes of the individual cases allowed the identification of two missense variants affecting different genes, SHC4 in case 4 in group 1 and WDR45B in case 13 in group 3. Furthermore, these two variants were only observed in Holstein cattle when querying WGS data of > 5,100 animals. Alternatively, potential de novo mutational events were assessed in each case. Filtering for heterozygous private protein-changing variants identified one DYNC1H1 frameshift variant as a candidate causal dominant acting allele in case 12 in group 3. Finally, the presence of larger structural DNA variants and chromosomal abnormalities was investigated in all cases. Depth of coverage analysis revealed two different partial monosomies of chromosome 2 segments in cases 1 and 7 in group 1 and a trisomy of chromosome 12 in the WDR45B homozygous case 13 in group 3.

CONCLUSIONS

This study presents for the first time a detailed genomic evaluation of CSCM in Holstein cattle and suggests an unexpected genetic and allelic heterogeneity considering the mode of inheritance, as well as the type of variant. For the first time, we propose candidate causal variants that may explain bovine CSCM in a certain proportion of affected calves. We present cattle as a large animal model for human CMII and propose new genes and genomic variants as possible causes for related diseases in both animals and humans.

Androgen receptor-mediated pharmacogenomic expression quantitative trait loci: implications for breast cancer response to AR-targeting therapy

BACKGROUND

Endocrine therapy is the most important treatment modality of breast cancer patients whose tumors express the estrogen receptor α (ERα). The androgen receptor (AR) is also expressed in the vast majority (80-90%) of ERα-positive tumors. AR-targeting drugs are not used in clinical practice, but have been evaluated in multiple trials and preclinical studies.

METHODS

We performed a genome-wide study to identify hormone/drug-induced single nucleotide polymorphism (SNP) genotype - dependent gene-expression, known as PGx-eQTL, mediated by either an AR agonist (dihydrotestosterone) or a partial antagonist (enzalutamide), utilizing a previously well characterized lymphoblastic cell line panel. The association of the identified SNPs-gene pairs with breast cancer phenotypes were then examined using three genome-wide association (GWAS) studies that we have published and other studies from the GWAS catalog.

RESULTS

We identified 13 DHT-mediated PGx-eQTL loci and 23 Enz-mediated PGx-eQTL loci that were associated with breast cancer outcomes post ER antagonist or aromatase inhibitors (AI) treatment, or with pharmacodynamic (PD) effects of AIs. An additional 30 loci were found to be associated with cancer risk and sex-hormone binding globulin levels. The top loci involved the genes IDH2 and TMEM9, the expression of which were suppressed by DHT in a PGx-eQTL SNP genotype-dependent manner. Both of these genes were overexpressed in breast cancer and were associated with a poorer prognosis. Therefore, suppression of these genes by AR agonists may benefit patients with minor allele genotypes for these SNPs.

CONCLUSIONS

We identified AR-related PGx-eQTL SNP-gene pairs that were associated with risks, outcomes and PD effects of endocrine therapy that may provide potential biomarkers for individualized treatment of breast cancer.

Rapid Adaptation and Interspecific Introgression in the North American Crop Pest Helicoverpa zea

Insect crop pests threaten global food security. This threat is amplified through the spread of nonnative species and through adaptation of native pests to control measures. Adaptations such as pesticide resistance can result from selection on variation within a population, or through gene flow from another population. We investigate these processes in an economically important noctuid crop pest, Helicoverpa zea, which has evolved resistance to a wide range of pesticides. Its sister species Helicoverpa armigera, first detected as an invasive species in Brazil in 2013, introduced the pyrethroid-resistance gene CYP337B3 to South American H. zea via adaptive introgression. To understand whether this could contribute to pesticide resistance in North America, we sequenced 237 H. zea genomes across 10 sample sites. We report H. armigera introgression into the North American H. zea population. Two individuals sampled in Texas in 2019 carry H. armigera haplotypes in a 4 Mbp region containing CYP337B3. Next, we identify signatures of selection in the panmictic population of nonadmixed H. zea, identifying a selective sweep at a second cytochrome P450 gene: CYP333B3. We estimate that its derived allele conferred a ∼5% fitness advantage and show that this estimate explains independently observed rare nonsynonymous CYP333B3 mutations approaching fixation over a ∼20-year period. We also detect putative signatures of selection at a kinesin gene associated with Bt resistance. Overall, we document two mechanisms of rapid adaptation: the introduction of fitness-enhancing alleles through interspecific introgression, and selection on intraspecific variation.

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

The multi-millenia long history between dogs and humans has 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 diversity of breeds today. Although ancient genomes provide valuable information, their use is hindered by low depth of coverage and post-mortem damage, which 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 provided by ancient datasets. Imputation accuracy was evaluated by down-sampling high coverage dog and wolf genomes to 0.05-2x coverage and comparing concordance between imputed and high coverage genotypes. We measured the impact of imputation on principal component analyses and runs of homozygosity. 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 dataset of 90 ancient dog and wolf genomes, 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. 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.

GenoTools: An Open-Source Python Package for Efficient Genotype Data Quality Control and Analysis

GenoTools, a Python package, streamlines population genetics research by integrating ancestry estimation, quality control (QC), and genome-wide association studies (GWAS) capabilities into efficient pipelines. By tracking samples, variants, and quality-specific measures throughout fully customizable pipelines, users can easily manage genetics data for large and small studies. GenoTools' "Ancestry" module renders highly accurate predictions, allowing for high-quality ancestry-specific studies, and enables custom ancestry model training and serialization, specified to the user's genotyping or sequencing platform. As the genotype processing engine that powers several large initiatives, including the NIH's Center for Alzheimer's and Related Dementias (CARD) and the Global Parkinson's Genetics Program (GP2). GenoTools was used to process and analyze the UK Biobank and major Alzheimer's Disease (AD) and Parkinson's Disease (PD) datasets with over 400,000 genotypes from arrays and 5000 sequences and has led to novel discoveries in diverse populations. It has provided replicable ancestry predictions, implemented rigorous QC, and conducted genetic ancestry-specific GWAS to identify systematic errors or biases through a single command. GenoTools is a customizable tool that enables users to efficiently analyze and scale genotype data with reproducible and scalable ancestry, QC, and GWAS pipelines.

Ancestral Origins and Admixture History of Kazakhs

Kazakh people, like many other populations that settled in Central Asia, demonstrate an array of mixed anthropological features of East Eurasian (EEA) and West Eurasian (WEA) populations, indicating a possible scenario of biological admixture between already differentiated EEA and WEA populations. However, their complex biological origin, genomic makeup, and genetic interaction with surrounding populations are not well understood. To decipher their genetic structure and population history, we conducted, to our knowledge, the first whole-genome sequencing study of Kazakhs residing in Xinjiang (KZK). We demonstrated that KZK derived their ancestries from 4 ancestral source populations: East Asian (∼39.7%), West Asian (∼28.6%), Siberian (∼23.6%), and South Asian (∼8.1%). The recognizable interactions of EEA and WEA ancestries in Kazakhs were dated back to the 15th century BCE. Kazakhs were genetically distinctive from the Uyghurs in terms of their overall genomic makeup, although the 2 populations were closely related in genetics, and both showed a substantial admixture of western and eastern peoples. Notably, we identified a considerable sex-biased admixture, with an excess of western males and eastern females contributing to the KZK gene pool. We further identified a set of genes that showed remarkable differentiation in KZK from the surrounding populations, including those associated with skin color (SLC24A5, OCA2), essential hypertension (HLA-DQB1), hypertension (MTHFR, SLC35F3), and neuron development (CNTNAP2). These results advance our understanding of the complex history of contacts between Western and Eastern Eurasians, especially those living or along the old Silk Road.

Genome sequence analyses identify novel risk loci for multiple system atrophy

Multiple system atrophy (MSA) is an adult-onset, sporadic synucleinopathy characterized by parkinsonism, cerebellar ataxia, and dysautonomia. The genetic architecture of MSA is poorly understood, and treatments are limited to supportive measures. Here, we performed a comprehensive analysis of whole genome sequence data from 888 European-ancestry MSA cases and 7,128 controls to systematically investigate the genetic underpinnings of this understudied neurodegenerative disease. We identified four significantly associated risk loci using a genome-wide association study approach. Transcriptome-wide association analyses prioritized USP38-DT, KCTD7, and lnc-KCTD7-2 as novel susceptibility genes for MSA within these loci, and single-nucleus RNA sequence analysis found that the associated variants acted as cis-expression quantitative trait loci for multiple genes across neuronal and glial cell types. In conclusion, this study highlights the role of genetic determinants in the pathogenesis of MSA, and the publicly available data from this study represent a valuable resource for investigating synucleinopathies.