ProbABEL package for genome-wide association analysis of imputed data

Exploring Genetic Factors Associated with Moniezia spp. Tapeworm Resistance in Central Anatolian Merino Sheep via GWAS Approach

Gastrointestinal parasite (GIP) infections pose significant challenges in pasture-based sheep farming, leading to economic losses and welfare concerns. This study aimed to uncover the genetic basis of resistance to Moniezia spp. infections in Central Anatolian Merino (CAM) sheep. Genome-Wide Association Analysis (GWAS) was conducted between Moniezia spp. egg burden and genomic data from 226 CAM lambs. Thirteen significant Single-Nucleotide Polymorphisms (SNPs) were identified, with five surpassing the genome-wide threshold and eight exceeding the chromosome-wide threshold. Functional annotation revealed associations with genes involved in immune function, notably CD79A and MAP3K7. CD79A, integral to B-cell activation and antibody production, plays a key role in the immune response against parasitic infections. Its interaction with helminth-derived proteins modulates B-cell function, highlighting its potential as a therapeutic target. MAP3K7, a central regulator of immune signaling pathways, modulates host responses to helminth infections by influencing NF-κB activity. Additionally, it regulates macrophage function in bacterial infections, showcasing its versatility in mediating immune responses against diverse pathogens. From a practical perspective, the findings of the current research underscore the potential of integrating genomic information into breeding programs to bolster disease resilience in livestock populations for sustainable production purposes. However, further research is needed to elucidate the functional significance of identified SNPs and associated genes. This study underscores the potential of genomic approaches in combating parasitic diseases and promoting sustainable agriculture in sheep production systems.

Genome-Wide Interaction Analyses of Serum Calcium on Ventricular Repolarization Time in 125 393 Participants

BACKGROUND

Ventricular repolarization time (ECG QT and JT intervals) is associated with malignant arrhythmia. Genome-wide association studies have identified 230 independent loci for QT and JT; however, 50% of their heritability remains unexplained. Previous work supports a causal effect of lower serum calcium concentrations on longer ventricular repolarization time. We hypothesized calcium interactions with QT and JT variant associations could explain a proportion of the missing heritability.

METHODS AND RESULTS

We performed genome-wide calcium interaction analyses for QT and JT intervals. Participants were stratified by their calcium level relative to the study distribution (top or bottom 20%). We performed a 2-stage analysis (genome-wide discovery [N=62 532] and replication [N=59 861] of lead variants) and a single-stage genome-wide meta-analysis (N=122 393, [European ancestry N=117 581, African ancestry N=4812]). We also calculated 2-degrees of freedom joint main and interaction and 1-degree of freedom interaction P values. In 2-stage and single-stage analyses, 50 and 98 independent loci, respectively, were associated with either QT or JT intervals (2-degrees of freedom joint main and interaction P value <5×10-8). No lead variant had a significant interaction result after correcting for multiple testing and sensitivity analyses provided similar findings. Two loci in the single-stage meta-analysis were not reported previously (SPPL2B and RFX6).

CONCLUSIONS

We have found limited support for an interaction effect of serum calcium on QT and JT variant associations despite sample sizes with suitable power to detect relevant effects. Therefore, such effects are unlikely to explain a meaningful proportion of the heritability of QT and JT, and factors including rare variation and other environmental interactions need to be considered.

Genetic regulation of carnitine metabolism controls lipid damage repair and aging RBC hemolysis in vivo and in vitro

ABSTRACT

Recent large-scale multiomics studies suggest that genetic factors influence the chemical individuality of donated blood. To examine this concept, we performed metabolomics analyses of 643 blood units from volunteers who donated units of packed red blood cells (RBCs) on 2 separate occasions. These analyses identified carnitine metabolism as the most reproducible pathway across multiple donations from the same donor. We also measured l-carnitine and acyl-carnitines in 13 091 packed RBC units from donors in the Recipient Epidemiology and Donor Evaluation study. Genome-wide association studies against 879 000 polymorphisms identified critical genetic factors contributing to interdonor heterogeneity in end-of-storage carnitine levels, including common nonsynonymous polymorphisms in genes encoding carnitine transporters (SLC22A16, SLC22A5, and SLC16A9); carnitine synthesis (FLVCR1 and MTDH) and metabolism (CPT1A, CPT2, CRAT, and ACSS2), and carnitine-dependent repair of lipids oxidized by ALOX5. Significant associations between genetic polymorphisms on SLC22 transporters and carnitine pools in stored RBCs were validated in 525 Diversity Outbred mice. Donors carrying 2 alleles of the rs12210538 SLC22A16 single-nucleotide polymorphism exhibited the lowest l-carnitine levels, significant elevations of in vitro hemolysis, and the highest degree of vesiculation, accompanied by increases in lipid peroxidation markers. Separation of RBCs by age, via in vivo biotinylation in mice, and Percoll density gradients of human RBCs, showed age-dependent depletions of l-carnitine and acyl-carnitine pools, accompanied by progressive failure of the reacylation process after chemically induced membrane lipid damage. Supplementation of stored murine RBCs with l-carnitine boosted posttransfusion recovery, suggesting this could represent a viable strategy to improve RBC storage quality.

Biological and Genetic Determinants of Glycolysis: Phosphofructokinase Isoforms Boost Energy Status of Stored Red Blood Cells and Transfusion Outcomes

Mature red blood cells (RBCs) lack mitochondria, and thus exclusively rely on glycolysis to generate adenosine triphosphate (ATP) during aging in vivo or storage in the blood bank. Here we leveraged 13,029 volunteers from the Recipient Epidemiology and Donor Evaluation Study to identify an association between end-of-storage levels of glycolytic metabolites and donor age, sex, and ancestry-specific genetic polymorphisms in regions encoding phosphofructokinase 1, platelet (detected in mature RBCs), hexokinase 1, ADP-ribosyl cyclase 1 and 2 (CD38/BST1). Gene-metabolite associations were validated in fresh and stored RBCs from 525 Diversity Outbred mice, and via multi-omics characterization of 1,929 samples from 643 human RBC units during storage. ATP and hypoxanthine levels - and the genetic traits linked to them - were associated with hemolysis in vitro and in vivo, both in healthy autologous transfusion recipients and in 5,816 critically ill patients receiving heterologous transfusions, suggesting their potential as markers to improve transfusion outcomes.

eTOC and Highlights

Highlights

Blood donor age and sex affect glycolysis in stored RBCs from 13,029 volunteers;Ancestry, genetic polymorphisms in PFKP, HK1, CD38/BST1 influence RBC glycolysis;Modeled PFKP effects relate to preventing loss of the total AXP pool in stored RBCs;ATP and hypoxanthine are biomarkers of hemolysis in vitro and in vivo.

Ferroptosis regulates hemolysis in stored murine and human red blood cells

Red blood cell (RBC) metabolism regulates hemolysis during aging in vivo and in the blood bank. Here, we leveraged a diversity outbred mouse population to map the genetic drivers of fresh/stored RBC metabolism and extravascular hemolysis upon storage and transfusion in 350 mice. We identify the ferrireductase Steap3 as a critical regulator of a ferroptosis-like process of lipid peroxidation. Steap3 polymorphisms were associated with RBC iron content, in vitro hemolysis, and in vivo extravascular hemolysis both in mice and 13,091 blood donors from the Recipient Epidemiology and Donor evaluation Study. Using metabolite Quantitative Trait Loci analyses, we identified a network of gene products (FADS1/2, EPHX2 and LPCAT3) - enriched in donors of African descent - associated with oxylipin metabolism in stored human RBCs and related to Steap3 or its transcriptional regulator, the tumor protein TP53. Genetic variants were associated with lower in vivo hemolysis in thousands of single-unit transfusion recipients.

Highlights

Steap3 regulates lipid peroxidation and extravascular hemolysis in 350 diversity outbred miceSteap3 SNPs are linked to RBC iron, hemolysis, vesiculation in 13,091 blood donorsmQTL analyses of oxylipins identified ferroptosis-related gene products FADS1/2, EPHX2, LPCAT3Ferroptosis markers are linked to hemoglobin increments in transfusion recipients.

Graphical abstract

Genome-Wide Analyses of Vocabulary Size in Infancy and Toddlerhood: Associations With Attention-Deficit/Hyperactivity Disorder, Literacy, and Cognition-Related Traits

BACKGROUND

The number of words children produce (expressive vocabulary) and understand (receptive vocabulary) changes rapidly during early development, partially due to genetic factors. Here, we performed a meta-genome-wide association study of vocabulary acquisition and investigated polygenic overlap with literacy, cognition, developmental phenotypes, and neurodevelopmental conditions, including attention-deficit/hyperactivity disorder (ADHD).

METHODS

We studied 37,913 parent-reported vocabulary size measures (English, Dutch, Danish) for 17,298 children of European descent. Meta-analyses were performed for early-phase expressive (infancy, 15-18 months), late-phase expressive (toddlerhood, 24-38 months), and late-phase receptive (toddlerhood, 24-38 months) vocabulary. Subsequently, we estimated single nucleotide polymorphism-based heritability (SNP-h2) and genetic correlations (rg) and modeled underlying factor structures with multivariate models.

RESULTS

Early-life vocabulary size was modestly heritable (SNP-h2 = 0.08-0.24). Genetic overlap between infant expressive and toddler receptive vocabulary was negligible (rg = 0.07), although each measure was moderately related to toddler expressive vocabulary (rg = 0.69 and rg = 0.67, respectively), suggesting a multifactorial genetic architecture. Both infant and toddler expressive vocabulary were genetically linked to literacy (e.g., spelling: rg = 0.58 and rg = 0.79, respectively), underlining genetic similarity. However, a genetic association of early-life vocabulary with educational attainment and intelligence emerged only during toddlerhood (e.g., receptive vocabulary and intelligence: rg = 0.36). Increased ADHD risk was genetically associated with larger infant expressive vocabulary (rg = 0.23). Multivariate genetic models in the ALSPAC (Avon Longitudinal Study of Parents and Children) cohort confirmed this finding for ADHD symptoms (e.g., at age 13; rg = 0.54) but showed that the association effect reversed for toddler receptive vocabulary (rg = -0.74), highlighting developmental heterogeneity.

CONCLUSIONS

The genetic architecture of early-life vocabulary changes during development, shaping polygenic association patterns with later-life ADHD, literacy, and cognition-related traits.

Regulation of kynurenine metabolism by blood donor genetics and biology impacts red cell hemolysis in vitro and in vivo

ABSTRACT

In the field of transfusion medicine, the clinical relevance of the metabolic markers of the red blood cell (RBC) storage lesion is incompletely understood. Here, we performed metabolomics of RBC units from 643 donors enrolled in the Recipient Epidemiology and Donor Evaluation Study, REDS RBC Omics. These units were tested on storage days 10, 23, and 42 for a total of 1929 samples and also characterized for end-of-storage hemolytic propensity after oxidative and osmotic insults. Our results indicate that the metabolic markers of the storage lesion poorly correlated with hemolytic propensity. In contrast, kynurenine was not affected by storage duration and was identified as the top predictor of osmotic fragility. RBC kynurenine levels were affected by donor age and body mass index and were reproducible within the same donor across multiple donations from 2 to 12 months apart. To delve into the genetic underpinnings of kynurenine levels in stored RBCs, we thus tested kynurenine levels in stored RBCs on day 42 from 13 091 donors from the REDS RBC Omics study, a population that was also genotyped for 879 000 single nucleotide polymorphisms. Through a metabolite quantitative trait loci analysis, we identified polymorphisms in SLC7A5, ATXN2, and a series of rate-limiting enzymes (eg, kynurenine monooxygenase, indoleamine 2,3-dioxygenase, and tryptophan dioxygenase) in the kynurenine pathway as critical factors affecting RBC kynurenine levels. By interrogating a donor-recipient linkage vein-to-vein database, we then report that SLC7A5 polymorphisms are also associated with changes in hemoglobin and bilirubin levels, suggestive of in vivo hemolysis in 4470 individuals who were critically ill and receiving single-unit transfusions.

Genetic factors affecting survival in Japanese patients with sporadic amyotrophic lateral sclerosis: a genome-wide association study and verification in iPSC-derived motor neurons from patients

BACKGROUND

Several genetic factors are associated with the pathogenesis of sporadic amyotrophic lateral sclerosis (ALS) and its phenotypes, such as disease progression. Here, in this study, we aimed to identify the genes that affect the survival of patients with sporadic ALS.

METHODS

We enrolled 1076 Japanese patients with sporadic ALS with imputed genotype data of 7 908 526 variants. We used Cox proportional hazards regression analysis with an additive model adjusted for sex, age at onset and the first two principal components calculated from genotyped data to conduct a genome-wide association study. We further analysed messenger RNA (mRNA) and phenotype expression in motor neurons derived from induced pluripotent stem cells (iPSC-MNs) of patients with ALS.

RESULTS

Three novel loci were significantly associated with the survival of patients with sporadic ALS-FGF1 at 5q31.3 (rs11738209, HR=2.36 (95% CI, 1.77 to 3.15), p=4.85×10-9), THSD7A at 7p21.3 (rs2354952, 1.38 (95% CI, 1.24 to 1.55), p=1.61×10-8) and LRP1 at 12q13.3 (rs60565245, 2.18 (95% CI, 1.66 to 2.86), p=2.35×10-8). FGF1 and THSD7A variants were associated with decreased mRNA expression of each gene in iPSC-MNs and reduced in vitro survival of iPSC-MNs obtained from patients with ALS. The iPSC-MN in vitro survival was reduced when the expression of FGF1 and THSD7A was partially disrupted. The rs60565245 was not associated with LRP1 mRNA expression.

CONCLUSIONS

We identified three loci associated with the survival of patients with sporadic ALS, decreased mRNA expression of FGF1 and THSD7A and the viability of iPSC-MNs from patients. The iPSC-MN model reflects the association between patient prognosis and genotype and can contribute to target screening and validation for therapeutic intervention.

Mechanosensitive ion channel gene survey suggests potential roles in primary open angle glaucoma

Although glaucoma is a disease modulated by eye pressure, the mechanisms of pressure sensing in the eye are not well understood. Here, we investigated associations between mechanosensitive ion channel gene variants and primary open-angle glaucoma (POAG). Common (minor allele frequency > 5%) single nucleotide polymorphisms located within the genomic regions of 20 mechanosensitive ion channel genes in the K2P, TMEM63, PIEZO and TRP channel families were assessed using genotype data from the NEIGHBORHOOD consortium of 3853 cases and 33,480 controls. Rare (minor allele frequency < 1%) coding variants were assessed using exome array genotyping data for 2606 cases and 2606 controls. Association with POAG was analyzed using logistic regression adjusting for age and sex. Two rare PIEZO1 coding variants with protective effects were identified in the NEIGHBOR dataset: R1527H, (OR 0.17, P = 0.0018) and a variant that alters a canonical splice donor site, g.16-88737727-C-G Hg38 (OR 0.38, P = 0.02). Both variants showed similar effects in the UK Biobank and the R1527H also in the FinnGen database. Several common variants also reached study-specific thresholds for association in the NEIGHBORHOOD dataset. These results identify novel variants in several mechanosensitive channel genes that show associations with POAG, suggesting that these channels may be potential therapeutic targets.

Genetic insights into resting heart rate and its role in cardiovascular disease

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Genetic regulators of sputum mucin concentration and their associations with COPD phenotypes

Hyper-secretion and/or hyper-concentration of mucus is a defining feature of multiple obstructive lung diseases, including chronic obstructive pulmonary disease (COPD). Mucus itself is composed of a mixture of water, ions, salt and proteins, of which the gel-forming mucins, MUC5AC and MUC5B, are the most abundant. Recent studies have linked the concentrations of these proteins in sputum to COPD phenotypes, including chronic bronchitis (CB) and acute exacerbations (AE). We sought to determine whether common genetic variants influence sputum mucin concentrations and whether these variants are also associated with COPD phenotypes, specifically CB and AE. We performed a GWAS to identify quantitative trait loci for sputum mucin protein concentration (pQTL) in the Sub-Populations and InteRmediate Outcome Measures in COPD Study (SPIROMICS, n = 708 for total mucin, n = 215 for MUC5AC, MUC5B). Subsequently, we tested for associations of mucin pQTL with CB and AE using regression modeling (n = 822-1300). Replication analysis was conducted using data from COPDGene (n = 5740) and by examining results from the UK Biobank. We identified one genome-wide significant pQTL for MUC5AC (rs75401036) and two for MUC5B (rs140324259, rs10001928). The strongest association for MUC5B, with rs140324259 on chromosome 11, explained 14% of variation in sputum MUC5B. Despite being associated with lower MUC5B, the C allele of rs140324259 conferred increased risk of CB (odds ratio (OR) = 1.42; 95% confidence interval (CI): 1.10-1.80) as well as AE ascertained over three years of follow up (OR = 1.41; 95% CI: 1.02-1.94). Associations between rs140324259 and CB or AE did not replicate in COPDGene. However, in the UK Biobank, rs140324259 was associated with phenotypes that define CB, namely chronic mucus production and cough, again with the C allele conferring increased risk. We conclude that sputum MUC5AC and MUC5B concentrations are associated with common genetic variants, and the top locus for MUC5B may influence COPD phenotypes, in particular CB.

Genome-wide metabolite quantitative trait loci analysis (mQTL) in red blood cells from volunteer blood donors

The red blood cell (RBC)-Omics study, part of the larger NHLBI-funded Recipient Epidemiology and Donor Evaluation Study (REDS-III), aims to understand the genetic contribution to blood donor RBC characteristics. Previous work identified donor demographic, behavioral, genetic, and metabolic underpinnings to blood donation, storage, and (to a lesser extent) transfusion outcomes, but none have yet linked the genetic and metabolic bodies of work. We performed a genome-wide association (GWA) analysis using RBC-Omics study participants with generated untargeted metabolomics data to identify metabolite quantitative trait loci in RBCs. We performed GWA analyses of 382 metabolites in 243 individuals imputed using the 1000 Genomes Project phase 3 all-ancestry reference panel. Analyses were conducted using ProbABEL and adjusted for sex, age, donation center, number of whole blood donations in the past 2 years, and first 10 principal components of ancestry. Our results identified 423 independent genetic loci associated with 132 metabolites (p < 5×10^-8). Potentially novel locus-metabolite associations were identified for the region encoding heme transporter FLVCR1 and choline and for lysophosphatidylcholine acetyltransferase LPCAT3 and lysophosphatidylserine 16.0, 18.0, 18.1, and 18.2; these associations are supported by published rare disease and mouse studies. We also confirmed previous metabolite GWA results for associations, including N(6)-methyl-L-lysine and protein PYROXD2 and various carnitines and transporter SLC22A16. Association between pyruvate levels and G6PD polymorphisms was validated in an independent cohort and novel murine models of G6PD deficiency (African and Mediterranean variants). We demonstrate that it is possible to perform metabolomics-scale GWA analyses with a modest, trans-ancestry sample size.

HCLC-FC: A novel statistical method for phenome-wide association studies

The emergence of genetic data coupled to longitudinal electronic medical records (EMRs) offers the possibility of phenome-wide association studies (PheWAS). In PheWAS, the whole phenome can be divided into numerous phenotypic categories according to the genetic architecture across phenotypes. Currently, statistical analyses for PheWAS are mainly univariate analyses, which test the association between one genetic variant and one phenotype at a time. In this article, we derived a novel and powerful multivariate method for PheWAS. The proposed method involves three steps. In the first step, we apply the bottom-up hierarchical clustering method to partition a large number of phenotypes into disjoint clusters within each phenotypic category. In the second step, the clustering linear combination method is used to combine test statistics within each category based on the phenotypic clusters and obtain p-values from each phenotypic category. In the third step, we propose a new false discovery rate (FDR) control approach. We perform extensive simulation studies to compare the performance of our method with that of other existing methods. The results show that our proposed method controls FDR very well and outperforms other methods we compared with. We also apply the proposed approach to a set of EMR-based phenotypes across more than 300,000 samples from the UK Biobank. We find that the proposed approach not only can well-control FDR at a nominal level but also successfully identify 1,244 significant SNPs that are reported to be associated with some phenotypes in the GWAS catalog. Our open-access tools and instructions on how to implement HCLC-FC are available at https://github.com/XiaoyuLiang/HCLCFC.

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.

Genetic determinants of chromatin reveal prostate cancer risk mediated by context-dependent gene regulation

Many genetic variants affect disease risk by altering context-dependent gene regulation. Such variants are difficult to study mechanistically using current methods that link genetic variation to steady-state gene expression levels, such as expression quantitative trait loci (eQTLs). To address this challenge, we developed the cistrome-wide association study (CWAS), a framework for identifying genotypic and allele-specific effects on chromatin that are also associated with disease. In prostate cancer, CWAS identified regulatory elements and androgen receptor-binding sites that explained the association at 52 of 98 known prostate cancer risk loci and discovered 17 additional risk loci. CWAS implicated key developmental transcription factors in prostate cancer risk that are overlooked by eQTL-based approaches due to context-dependent gene regulation. We experimentally validated associations and demonstrated the extensibility of CWAS to additional epigenomic datasets and phenotypes, including response to prostate cancer treatment. CWAS is a powerful and biologically interpretable paradigm for studying variants that influence traits by affecting transcriptional regulation.

Deciphering the genetic basis of male fertility in Italian Brown Swiss dairy cattle

Improving reproductive performance remains a major goal in dairy cattle worldwide. Service sire has been recognized as an important factor affecting herd fertility. The main objective of this study was to reveal the genetic basis of male fertility in Italian Brown Swiss dairy cattle. Dataset included 1102 Italian Brown Swiss bulls with sire conception rate records genotyped with 454k single nucleotide polymorphisms. The analysis included whole-genome scans and gene-set analyses to identify genomic regions, individual genes and genetic mechanisms affecting Brown Swiss bull fertility. One genomic region on BTA1 showed significant additive effects. This region harbors gene RABL3 which is implicated cell proliferation and motility. Two genomic regions, located on BTA6 and BTA26, showed marked non-additive effects. These regions harbor genes, such as WDR19 and ADGRA1, that are directly involved in male fertility, including sperm motility, acrosome reaction, and embryonic development. The gene-set analysis revealed functional terms related to cell adhesion, cellular signaling, cellular transport, immune system, and embryonic development. Remarkably, a gene-set analysis also including Holstein and Jersey data, revealed significant processes that are common to the three dairy breeds, including cell migration, cell-cell interaction, GTPase activity, and the immune function. Overall, this comprehensive study contributes to a better understanding of the genetic basis of male fertility in cattle. In addition, our findings may guide the development of novel genomic strategies for improving service sire fertility in Brown Swiss cattle.

Genome-wide association analysis reveals insights into the genetic architecture of right ventricular structure and function

Right ventricular (RV) structure and function influence the morbidity and mortality from coronary artery disease (CAD), dilated cardiomyopathy (DCM), pulmonary hypertension and heart failure. Little is known about the genetic basis of RV measurements. Here we perform genome-wide association analyses of four clinically relevant RV phenotypes (RV end-diastolic volume, RV end-systolic volume, RV stroke volume, RV ejection fraction) from cardiovascular magnetic resonance images, using a state-of-the-art deep learning algorithm in 29,506 UK Biobank participants. We identify 25 unique loci associated with at least one RV phenotype at P < 2.27 ×10-8, 17 of which are validated in a combined meta-analysis (n = 41,830). Several candidate genes overlap with Mendelian cardiomyopathy genes and are involved in cardiac muscle contraction and cellular adhesion. The RV polygenic risk scores (PRSs) are associated with DCM and CAD. The findings substantially advance our understanding of the genetic underpinning of RV measurements.

Interaction between cigarette smoking and genetic polymorphisms on the associations with age of natural menopause and reproductive lifespan: the Singapore Chinese Health Study

STUDY QUESTION

Are there genetic variants that interact with smoking to reduce reproductive lifespan in East-Asian women?

SUMMARY ANSWER

Our study corroborates several recently identified genetic loci associated with reproductive lifespan and highlights specific genetic predispositions that may interact with smoking status to adversely affect reproductive lifespan in East-Asian women.

WHAT IS KNOWN ALREADY

Epidemiological data as well as evaluations on genetic predisposition to smoke indicate on the importance of smoking in adverse effects on reproductive lifespan in women. However, there are no previous smoking and gene interaction studies for reproductive traits in East-Asian women.

STUDY DESIGN, SIZE, DURATION

This population-based prospective cohort study comprised 11 643 East-Asian Chinese women with overlapping genome-wide genotyping and reproductive data.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We performed a genome-wide association study for reproductive lifespan in women (n = 11 643) from the Singapore Chinese Health Study (SCHS) and carried out a genome-wide interaction study to identify loci that interacted with smoking status to affect age of natural menopause and reproductive-time.

MAIN RESULTS AND THE ROLE OF CHANCE

Two known loci associated with menopause, rs113430717 (near HMCES, chromosome 3, Pmeta = 5.72 × 10-15) and rs3020136 (near RAD21, chromosome 8, Pmeta = 1.38 × 10-8) were observed beyond genome-wide levels of association with age at menopause in this study. For reproductive lifespan, the genome-wide association observed at rs79784106 (chromosome 3, Pmeta = 5.05 × 10-12) was in linkage disequilibrium with the menopause lead single-nucleotide polymorphism (SNP) (rs113430717). Four additional loci, first reported to be associated with menopause, were also associated with reproductive lifespan in our study (PAdj between 7.42 × 10-5 to 4.51 × 10-3). A significant interaction was observed between smoking and an East-Asian specific SNP, rs140146885, for reduced reproductive lifespan, per copy of the minor C allele (beta = -1.417 years, Pinteraction = 2.31 × 10-10). This interaction was successfully replicated in additional independent samples (beta = -1.389 years, Pinteraction = 6.78 × 10-3). Another known variant associated with menopause, rs11031006 (near FSHB), was also observed to interact with smoking status to reduce age at menopause in our dataset (beta = -0.450 years, Padj = 0.042).

LIMITATIONS, REASONS FOR CAUTION

The modest sample size of the replication datasets used likely affected the statistical power to firmly replicate all identified novel loci observed in our smoking interaction analyses.

WIDER IMPLICATIONS OF THE FINDINGS

Age of natural menopause and reproductive lifespan have clear genetic predispositions with distinct ethnic differences, and they may be adversely truncated by lifestyle factors such as smoking, which can pose a significant impact on the reproductive lifespan and future health outcomes in women.

STUDY FUNDING/COMPETING INTEREST(S)

The Singapore Chinese Health Study is funded by the National Medical Research Council, Singapore (NMRC/CIRG/1456/2016), National Institutes of Health (R01 CA144034 and UM1 CA182876) and National Research Foundation, Singapore (Project Number 370062002). W.-P.K. is supported by the National Medical Research Council, Singapore (MOH-CSASI19nov-0001). The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication. The authors do not report conflicts of interest.

TRIAL REGISTRATION NUMBER

N/A.

Genetic modifiers of Huntington disease differentially influence motor and cognitive domains

Genome-wide association studies (GWASs) of Huntington disease (HD) have identified six DNA maintenance gene loci (among others) as modifiers and implicated a two step-mechanism of pathogenesis: somatic instability of the causative HTT CAG repeat with subsequent triggering of neuronal damage. The largest studies have been limited to HD individuals with a rater-estimated age at motor onset. To capitalize on the wealth of phenotypic data in several large HD natural history studies, we have performed algorithmic prediction by using common motor and cognitive measures to predict age at other disease landmarks as additional phenotypes for GWASs. Combined with imputation with the Trans-Omics for Precision Medicine reference panel, predictions using integrated measures provided objective landmark phenotypes with greater power to detect most modifier loci. Importantly, substantial differences in the relative modifier signal across loci, highlighted by comparing common modifiers at MSH3 and FAN1, revealed that individual modifier effects can act preferentially in the motor or cognitive domains. Individual components of the DNA maintenance modifier mechanisms may therefore act differentially on the neuronal circuits underlying the corresponding clinical measures. In addition, we identified additional modifier effects at the PMS1 and PMS2 loci and implicated a potential second locus on chromosome 7. These findings indicate that broadened discovery and characterization of HD genetic modifiers based on additional quantitative or qualitative phenotypes offers not only the promise of in-human validated therapeutic targets but also a route to dissecting the mechanisms and cell types involved in both the somatic instability and toxicity components of HD pathogenesis.

Sex-specific genetic modifiers identified susceptibility of cold stored red blood cells to osmotic hemolysis

Background

Genetic variants have been found to influence red blood cell (RBC) susceptibility to hemolytic stress and affect transfusion outcomes and the severity of blood diseases. Males have a higher susceptibility to hemolysis than females, but little is known about the genetic mechanism contributing to the difference.

Results

To investigate the sex differences in RBC susceptibility to hemolysis, we conducted a sex-stratified genome-wide association study and a genome-wide gene-by-sex interaction scan in a multi-ethnic dataset with 12,231 blood donors who have in vitro osmotic hemolysis measurements during routine blood storage. The estimated SNP-based heritability for osmotic hemolysis was found to be significantly higher in males than in females (0.46 vs. 0.41). We identified SNPs associated with sex-specific susceptibility to osmotic hemolysis in five loci (SPTA1, KCNA6, SLC4A1, SUMO1P1, and PAX8) that impact RBC function and hemolysis.

Conclusion

Our study established a best practice to identify sex-specific genetic modifiers for sexually dimorphic traits in datasets with mixed ancestries, providing evidence of different genetic regulations of RBC susceptibility to hemolysis between sexes. These and other variants may help explain observed sex differences in the severity of hemolytic diseases, such as sickle cell and malaria, as well as the viability of red cell storage and recovery.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08461-4.

Genome-wide pharmacogenetics of anti-drug antibody response to bococizumab highlights key residues in HLA DRB1 and DQB1

In this largest to-date genetic analysis of anti-drug antibody (ADA) response to a therapeutic monoclonal antibody (MAb), genome-wide association was performed for five measures of ADA status among 8844 individuals randomized to bococizumab, which targets PCSK9 for LDL-C lowering and cardiovascular protection. Index associations prioritized specific amino acid substitutions at the DRB1 and DQB1 MHC class II genes rather than canonical haplotypes. Two clusters of missense variants at DRB1 were associated with general ADA measures (residues 9, 11, 13; and 96, 112, 120, 180) and a third cluster of missense variants in DQB1 was associated with ADA measures including neutralizing antibody (NAb) titers (residues 66, 67, 71, 74, 75). The structural disposition of the missense substitutions implicates peptide antigen binding and CD4 effector function, mechanisms that are potentially generalizable to other therapeutic mAbs.

Clinicaltrials.gov: NCT01968954, NCT01968967, NCT01968980, NCT01975376, NCT01975389, NCT02100514.

Genome-wide association and functional interrogation identified a variant at 3p26.1 modulating ovarian cancer survival among Chinese women

Ovarian cancer survival varies considerably among patients, to which germline variation may also contribute in addition to mutational signatures. To identify genetic markers modulating ovarian cancer outcome, we performed a genome-wide association study in 2130 Chinese ovarian cancer patients and found a hitherto unrecognized locus at 3p26.1 to be associated with the overall survival (P_combined = 8.90 × 10⁻¹⁰). Subsequent statistical fine-mapping, functional annotation, and eQTL mapping prioritized a likely casual SNP rs9311399 in the non-coding regulatory region. Mechanistically, rs9311399 altered its enhancer activity through an allele-specific transcription factor binding and a long-range interaction with the promoter of a lncRNA BHLHE40-AS1. Deletion of the rs9311399-associated enhancer resulted in expression changes in several oncogenic signaling pathway genes and a decrease in tumor growth. Thus, we have identified a novel genetic locus that is associated with ovarian cancer survival possibly through a long-range gene regulation of oncogenic pathways.

Genetic Studies of Metabolomics Change After a Liquid Meal Illuminate Novel Pathways for Glucose and Lipid Metabolism

Humans spend the greater part of the day in a postprandial state. However, the genetic basis of postprandial blood measures is relatively uncharted territory. We examined the genetics of variation in concentrations of postprandial metabolites (t = 150 min) in response to a liquid mixed meal through genome-wide association studies (GWAS) performed in the Netherlands Epidemiology of Obesity (NEO) study (n = 5,705). The metabolite response GWAS identified an association between glucose change and rs10830963:G in the melatonin receptor 1B (β [SE] -0.23 [0.03], P = 2.15 × 10^-19). In addition, the ANKRD55 locus led by rs458741:C showed strong associations with extremely large VLDL (XXLVLDL) particle response (XXLVLDL total cholesterol: β [SE] 0.17 [0.03], P = 5.76 × 10^-10; XXLVLDL cholesterol ester: β [SE] 0.17 [0.03], P = 9.74 × 10^-10), which also revealed strong associations with body composition and diabetes in the UK Biobank (P < 5 × 10^-8). Furthermore, the associations between XXLVLDL response and insulinogenic index, HOMA-β, Matsuda insulin sensitivity index, and HbA_1c in the NEO study implied the role of chylomicron synthesis in diabetes (with false discovery rate-corrected q <0.05). To conclude, genetic studies of metabolomics change after a liquid meal illuminate novel pathways for glucose and lipid metabolism. Further studies are warranted to corroborate biological pathways of the ANKRD55 locus underlying diabetes.

Phenotypic and Genotypic Associations Between Migraine and Lipoprotein Subfractions

BACKGROUND AND OBJECTIVE

To evaluate phenotypic and genetic relationships between migraine and lipoprotein subfractions.

METHODS

We evaluated phenotypic associations between migraine and 19 lipoprotein subfraction measures in the Women's Genome Health Study (n = 22,788). We then investigated genetic relationships between these traits using summary statistics from the International Headache Genetics Consortium for migraine (n_case = 54,552, n_control = 297,970) and combined summary data for lipoprotein subfractions (n up to 47,713).

RESULTS

There was a significant phenotypic association (odds ratio 1.27 [95% confidence interval 1.12-1.44]) and a significant genetic correlation at 0.18 (p = 0.001) between migraine and triglyceride-rich lipoproteins (TRLPs) concentration but not for low-density lipoprotein or high-density lipoprotein subfractions. Mendelian randomization (MR) estimates were largely null, implying that pleiotropy rather than causality underlies the genetic correlation between migraine and lipoprotein subfractions. Pleiotropy was further supported in cross-trait meta-analysis, revealing significant shared signals at 4 loci (chr2p21 harboring THADA, chr5q13.3 harboring HMGCR, chr6q22.31 harboring HEY2, and chr7q11.23 harboring MLXIPL) between migraine and lipoprotein subfractions. Three of these loci were replicated for migraine (p < 0.05) in a smaller sample from the UK Biobank. The shared signal at chr5q13.3 colocalized with expression of HMGCR, ANKDD1B, and COL4A3BP in multiple tissues.

CONCLUSIONS

The study supports the association between certain lipoprotein subfractions, especially for TRLP, and migraine in populations of European ancestry. The corresponding shared genetic components may help identify potential targets for future migraine therapeutics.

CLASSIFICATION OF EVIDENCE

This study provides Class I evidence that migraine is significantly associated with some lipoprotein subfractions.

Multi-ancestry genome-wide gene-sleep interactions identify novel loci for blood pressure

Long and short sleep duration are associated with elevated blood pressure (BP), possibly through effects on molecular pathways that influence neuroendocrine and vascular systems. To gain new insights into the genetic basis of sleep-related BP variation, we performed genome-wide gene by short or long sleep duration interaction analyses on four BP traits (systolic BP, diastolic BP, mean arterial pressure, and pulse pressure) across five ancestry groups in two stages using 2 degree of freedom (df) joint test followed by 1df test of interaction effects. Primary multi-ancestry analysis in 62,969 individuals in stage 1 identified three novel gene by sleep interactions that were replicated in an additional 59,296 individuals in stage 2 (stage 1 + 2 Pjoint < 5 × 10-8), including rs7955964 (FIGNL2/ANKRD33) that increases BP among long sleepers, and rs73493041 (SNORA26/C9orf170) and rs10406644 (KCTD15/LSM14A) that increase BP among short sleepers (Pint < 5 × 10-8). Secondary ancestry-specific analysis identified another novel gene by long sleep interaction at rs111887471 (TRPC3/KIAA1109) in individuals of African ancestry (Pint = 2 × 10-6). Combined stage 1 and 2 analyses additionally identified significant gene by long sleep interactions at 10 loci including MKLN1 and RGL3/ELAVL3 previously associated with BP, and significant gene by short sleep interactions at 10 loci including C2orf43 previously associated with BP (Pint < 10-3). 2df test also identified novel loci for BP after modeling sleep that has known functions in sleep-wake regulation, nervous and cardiometabolic systems. This study indicates that sleep and primary mechanisms regulating BP may interact to elevate BP level, suggesting novel insights into sleep-related BP regulation.

Adverse Cardiovascular Outcomes and Antihypertensive Treatment: A Genome-Wide Interaction Meta-Analysis in the International Consortium for Antihypertensive Pharmacogenomics Studies

We sought to identify genome-wide variants influencing antihypertensive drug response and adverse cardiovascular outcomes, utilizing data from four randomized controlled trials in the International Consortium for Antihypertensive Pharmacogenomics Studies (ICAPS). Genome-wide antihypertensive drug-single nucleotide polymorphism (SNP) interaction tests for four drug classes (β-blockers, n = 9,195; calcium channel blockers (CCBs), n = 10,511; thiazide/thiazide-like diuretics, n = 3,516; ACE-inhibitors/ARBs, n = 2,559) and cardiovascular outcomes (incident myocardial infarction, stroke, or death) were analyzed among patients with hypertension of European ancestry. Top SNPs from the meta-analyses were tested for replication of cardiovascular outcomes in an independent Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) study (n = 21,267), blood pressure (BP) response in independent ICAPS studies (n = 1,552), and ethnic validation in African Americans from the Genetics of Hypertension Associated Treatment study (GenHAT; n = 5,115). One signal reached genome-wide significance in the β-blocker-SNP interaction analysis (rs139945292, Interaction P = 1.56 × 10-8 ). rs139945292 was validated through BP response to β-blockers, with the T-allele associated with less BP reduction (systolic BP response P = 6 × 10-4 , Beta = 3.09, diastolic BP response P = 5 × 10-3 , Beta = 1.53). The T-allele was also associated with increased adverse cardiovascular risk within the β-blocker treated patients' subgroup (P = 2.35 × 10-4 , odds ratio = 1.57, 95% confidence interval = 1.23-1.99). The locus showed nominal replication in CHARGE, and consistent directional trends in β-blocker treated African Americans. rs139945292 is an expression quantitative trait locus for the 50 kb upstream gene NTM (neurotrimin). No SNPs attained genome-wide significance for any other drugs classes. Top SNPs were located near CALB1 (CCB), FLJ367777 (ACE-inhibitor), and CES5AP1 (thiazide). The NTM region is associated with increased risk for adverse cardiovascular outcomes and less BP reduction in β-blocker treated patients. Further investigation into this region is warranted.

Sex-specific genetic factors affect the risk of early-onset periodontitis in Europeans

AIMS

Various studies have reported that young European women are more likely to develop early-onset periodontitis compared to men. A potential explanation for the observed variations in sex and age of disease onset is the natural genetic variation within the autosomal genomes. We hypothesized that genotype-by-sex (G × S) interactions contribute to the increased prevalence and severity.

MATERIALS AND METHODS

Using the case-only design, we tested for differences in genetic effects between men and women in 896 North-West European early-onset cases, using imputed genotypes from the OmniExpress genotyping array. Population-representative 6823 controls were used to verify that the interacting variables G and S were uncorrelated in the general population.

RESULTS

In total, 20 loci indicated G × S associations (P < 0.0005), 3 of which were previously suggested as risk genes for periodontitis (ABLIM2, CDH13, and NELL1). We also found independent G × S interactions of the related gene paralogs MACROD1/FLRT1 (chr11) and MACROD2/FLRT3 (chr20). G × S-associated SNPs at CPEB4, CDH13, MACROD1, and MECOM were genome-wide-associated with heel bone mineral density (CPEB4, MECOM), waist-to-hip ratio (CPEB4, MACROD1), and blood pressure (CPEB4, CDH13).

CONCLUSIONS

Our results indicate that natural genetic variation affects the different heritability of periodontitis among sexes and suggest genes that contribute to inter-sex phenotypic variation in early-onset periodontitis.

Genetics of PlGF plasma levels highlights a role of its receptors and supports the link between angiogenesis and immunity

Placental growth factor (PlGF) is a member of the vascular endothelial growth factor family and is involved in bone marrow-derived cell activation, endothelial stimulation and pathological angiogenesis. High levels of PlGF have been observed in several pathological conditions especially in cancer, cardiovascular, autoimmune and inflammatory diseases. Little is known about the genetics of circulating PlGF levels. Indeed, although the heritability of circulating PlGF levels is around 40%, no studies have assessed the relation between PlGF plasma levels and genetic variants at a genome-wide level. In the current study, PlGF plasma levels were measured in a population-based sample of 2085 adult individuals from three isolated populations of South Italy. A GWAS was performed in a discovery cohort (N = 1600), followed by a de novo replication (N = 468) from the same populations. The meta-analysis of the discovery and replication samples revealed one signal significantly associated with PlGF circulating levels. This signal was mapped to the PlGF co-receptor coding gene NRP1, indicating its important role in modulating the PlGF plasma levels. Two additional signals, at the PlGF receptor coding gene FLT1 and RAPGEF5 gene, were identified at a suggestive level. Pathway and TWAS analyses highlighted genes known to be involved in angiogenesis and immune response, supporting the link between these processes and PlGF regulation. Overall, these data improve our understanding of the genetic variation underlying circulating PlGF levels. This in turn could lead to new preventive and therapeutic strategies for a wide variety of PlGF-related pathologies.

A variant near DHCR24 associates with microstructural properties of white matter and peripheral lipid metabolism in adolescents

Visceral adiposity has been associated with altered microstructural properties of white matter in adolescents. Previous evidence suggests that circulating phospholipid PC(16:0/2:0) may mediate this association. To investigate the underlying biology, we performed a genome-wide association study (GWAS) of the shared variance of visceral fat, PC(16:0/2:0), and white matter microstructure in 872 adolescents from the Saguenay Youth Study. We further studied the metabolomic profile of the GWAS-lead variant in 931 adolescents. Visceral fat and white matter microstructure were assessed with magnetic resonance imaging. Circulating metabolites were quantified with serum lipidomics and metabolomics. We identified a genome-wide significant association near DHCR24 (Seladin-1) encoding a cholesterol-synthesizing enzyme (rs588709, p = 3.6 × 10-8); rs588709 was also associated nominally with each of the three traits (white matter microstructure: p = 2.1 × 10-6, PC(16:0/2:0): p = 0.005, visceral fat: p = 0.010). We found that the metabolic profile associated with rs588709 resembled that of a TM6SF2 variant impacting very low-density lipoprotein (VLDL) secretion and was only partially similar to that of a HMGCR variant. This suggests that the effect of rs588709 on VLDL lipids may arise due to altered phospholipid rather than cholesterol metabolism. The rs588709 was also nominally associated with circulating concentrations of omega-3 fatty acids in interaction with visceral fat and PC(16:0/2:0), and these fatty acid measures showed robust associations with white matter microstructure. Overall, the present study provides evidence that the DHCR24 locus may link peripheral metabolism to brain microstructure, an association with implications for cognitive impairment.

Genome-wide gene-smoking interaction study identified novel susceptibility loci for non-small cell lung cancer in Chinese populations

Gene-smoking interactions play important roles in the development of non-small cell lung cancer (NSCLC). To identify single nucleotide polymorphisms (SNPs) that modify the association of smoking behavior with NSCLC risk, we conducted a genome-wide gene-smoking interaction study in Chinese populations. The genome-wide interaction analysis between SNPs and smoking status (ever- versus never-smokers) was carried out using genome-wide association studies (GWAS) of NSCLC, which included 13,327 cases and 13,328 controls. Stratified analysis by histological subtypes was also conducted. We used a genome-wide significance threshold of 5×10 -8 for identifying significant gene-smoking interactions and 1×10 -6 for identifying suggestive results. Functional annotation was performed to identify potential functional SNPs and target genes. We identified three novel loci with significant or suggestive gene-smoking interaction. For NSCLC, the interaction between rs2746087 (20q11.23) and smoking status reached genome-wide significance threshold (OR = 0.63, 95%CI: 0.54-0.74, P = 3.31×10 -8), and the interaction between rs11912498 (22q12.1) and smoking status reached suggestive significance threshold (OR = 0.72, 95%CI: 0.63-0.82, P = 8.10×10 -7). Stratified analysis by histological subtypes identified suggestive interactions between rs459724 (5q11.2) and smoking status (OR = 0.61, 95%CI: 0.51-0.73, P = 7.55×10 -8) in the risk of lung squamous cell carcinoma. Functional annotation indicated that both classic and novel biological processes, including nicotine addiction and airway clearance, may modulate the susceptibility to NSCLC. These novel loci provide new insights into the biological mechanisms underlying NSCLC risk. Independent replication in large-scale studies is needed and experimental studies are warranted to functionally validate these associations.

Multiple-ancestry genome-wide association study identifies 27 loci associated with measures of hemolysis following blood storage

BackgroundThe evolutionary pressure of endemic malaria and other erythrocytic pathogens has shaped variation in genes encoding erythrocyte structural and functional proteins, influencing responses to hemolytic stress during transfusion and disease.MethodsWe sought to identify such genetic variants in blood donors by conducting a genome-wide association study (GWAS) of 12,353 volunteer donors, including 1,406 African Americans, 1,306 Asians, and 945 Hispanics, whose stored erythrocytes were characterized by quantitative assays of in vitro osmotic, oxidative, and cold-storage hemolysis.ResultsGWAS revealed 27 significant loci (P < 5 × 10-8), many in candidate genes known to modulate erythrocyte structure, metabolism, and ion channels, including SPTA1, ALDH2, ANK1, HK1, MAPKAPK5, AQP1, PIEZO1, and SLC4A1/band 3. GWAS of oxidative hemolysis identified variants in genes encoding antioxidant enzymes, including GLRX, GPX4, G6PD, and SEC14L4 (Golgi-transport protein). Genome-wide significant loci were also tested for association with the severity of steady-state (baseline) in vivo hemolytic anemia in patients with sickle cell disease, with confirmation of identified SNPs in HBA2, G6PD, PIEZO1, AQP1, and SEC14L4.ConclusionsMany of the identified variants, such as those in G6PD, have previously been shown to impair erythrocyte recovery after transfusion, associate with anemia, or cause rare Mendelian human hemolytic diseases. Candidate SNPs in these genes, especially in polygenic combinations, may affect RBC recovery after transfusion and modulate disease severity in hemolytic diseases, such as sickle cell disease and malaria.

CLUE: Exact maximal reduction of kinetic models by constrained lumping of differential equations

MOTIVATION

Detailed mechanistic models of biological processes can pose significant challenges for analysis and parameter estimations due to the large number of equations used to track the dynamics of all distinct configurations in which each involved biochemical species can be found. Model reduction can help tame such complexity by providing a lower-dimensional model in which each macro-variable can be directly related to the original variables.

RESULTS

We present CLUE, an algorithm for exact model reduction of systems of polynomial differential equations by constrained linear lumping. It computes the smallest dimensional reduction as a linear mapping of the state space such that the reduced model preserves the dynamics of user-specified linear combinations of the original variables. Even though CLUE works with nonlinear differential equations, it is based on linear algebra tools, which makes it applicable to high-dimensional models. Using case studies from the literature, we show how CLUE can substantially lower model dimensionality and help extract biologically intelligible insights from the reduction.

AVAILABILITY

An implementation of the algorithm and relevant resources to replicate the experiments herein reported are freely available for download at https://github.com/pogudingleb/CLUE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Association between MTNR1B polymorphisms and obesity in African American: findings from the Jackson Heart Study

BACKGROUND

Melatonin is a hormone that is secreted at night by the pineal gland. It exerts its function by binding to the MT1 and MT2 receptors, which are encoded by the MTNR1A and MTNR1B genes, respectively. Previous studies reveal that MTNR1B variants are associated with insulin secretion impairments and an increased body mass index (BMI) in individuals of European and Asian ancestries. Obesity is highly prevalent in the US and disproportionately affects African Americans. Here, we hypothesized that common single nucleotide polymorphisms (SNPs) imputed in 1000 Genomes in the MTNR1B gene are associated with adiposity in African American adult men and women and that the association is modified by insomnia.

METHODS

We used an additive genetic model to describe the association between the adiposity traits (BMI and waist circumference) and selected MTNR1B variants in 3,029 Jackson Heart Study participants, with an average age of 55.13 ± 12.84 years, and 62% were women. We regressed the adiposity measures on the estimated allelic or genotypic dosage at every selected SNP and adjusted for age, sex, population stratification, and insomnia. Thirty common SNPs, spanning the MTNR1B gene, with a minor allele frequency ≥ 5%, a call rate ≥ 90%, a Hardy-Weinberg equilibrium p value > 10-6, were available for the analysis.

RESULTS

The allele T of rs76371840 was associated with adiposity (OR = 1.47 [1.13-1.82]; PFDR-adjusted = 0.0499), and the allele A of rs8192552 showed a significant association with waist circumference (β = 0.023 ± 0.007; PFDR-adjusted = 0.0077) after correcting for multiple testing. When insomnia was included in the adiposity analysis model, the following four variants became significantly associated with adiposity: rs6483208; rs4388843; rs4601728; and rs12804291.

CONCLUSIONS

Our data indicate that polymorphisms in the MTNR1B gene are associated with obesity traits in African Americans. To the best of our knowledge, this is the first study to explore the effect of insomnia on the association between the circadian MTNR1B genetic variants and metabolic traits in an African American sample population. We observed that insomnia affected the association between the MTNR1B variants and adiposity.

Lifestyle Risk Score: handling missingness of individual lifestyle components in meta-analysis of gene-by-lifestyle interactions

Recent studies consider lifestyle risk score (LRS), an aggregation of multiple lifestyle exposures, in identifying association of gene-lifestyle interaction with disease traits. However, not all cohorts have data on all lifestyle factors, leading to increased heterogeneity in the environmental exposure in collaborative meta-analyses. We compared and evaluated four approaches (Naïve, Safe, Complete and Moderator Approaches) to handle the missingness in LRS-stratified meta-analyses under various scenarios. Compared to "benchmark" results with all lifestyle factors available for all cohorts, the Complete Approach, which included only cohorts with all lifestyle components, was underpowered due to lower sample size, and the Naïve Approach, which utilized all available data and ignored the missingness, was slightly inflated. The Safe Approach, which used all data in LRS-exposed group and only included cohorts with all lifestyle factors available in the LRS-unexposed group, and the Moderator Approach, which handled missingness via moderator meta-regression, were both slightly conservative and yielded almost identical p values. We also evaluated the performance of the Safe Approach under different scenarios. We observed that the larger the proportion of cohorts without missingness included, the more accurate the results compared to "benchmark" results. In conclusion, we generally recommend the Safe Approach, a straightforward and non-inflated approach, to handle heterogeneity among cohorts in the LRS based genome-wide interaction meta-analyses.

Effect of plasma polyunsaturated fatty acid levels on leukocyte telomere lengths in the Singaporean Chinese population

BACKGROUND

Shorter telomere length (TL) has been associated with poor health behaviors, increased risks of chronic diseases and early mortality. Excessive shortening of telomere is a marker of accelerated aging and can be influenced by oxidative stress and nutritional deficiency. Plasma n6:n3 polyunsaturated fatty acid (PUFA) ratio may impact cell aging. Increased dietary intake of marine n-3 PUFA is associated with reduced telomere attrition. However, the effect of plasma PUFA on leukocyte telomere length (LTL) and its interaction with genetic variants are not well established.

METHODS

A nested coronary artery disease (CAD) case-control study comprising 711 cases and 638 controls was conducted within the Singapore Chinese Health Study (SCHS). Samples genotyped with the Illumina ZhongHua-8 array. Plasma n-3 and n-6 PUFA were quantified using mass spectrometry (MS). LTL was measured with quantitative PCR method. Linear regression was used to test the association between PUFA and LTL. The interaction between plasma PUFAs and genetic variants was assessed by introducing an additional term (PUFA×genetic variant) in the regression model. Analysis was carried out in cases and controls separately and subsequently meta-analyzed using the inverse-variance weighted method. We further assessed the association of PUFA and LTL with CAD risk by Cox Proportional-Hazards model and whether the effect of PUFA on CAD was mediated through LTL by using structural equation modeling.

RESULTS

Higher n6:n3 ratio was significantly associated with shorter LTL (p = 0.018) and increased CAD risk (p = 0.005). These associations were mainly driven by elevated plasma total n-3 PUFAs, especially eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) (p < 0.05). There was a statistically significant interaction for an intergenic single nucleotide polymorphism (SNP) rs529143 with plasma total n-3 PUFA and DHA on LTL beyond the genome-wide threshold (p < 5 ×  10- 8). Mediation analysis showed that PUFA and LTL affected CAD risk independently.

CONCLUSIONS

Higher plasma n6:n3 PUFA ratio, and lower EPA and DHA n-3 PUFAs were associated with shorter LTL and increased CAD risk in this Chinese population. Furthermore, genetic variants may modify the effect of PUFAs on LTL. PUFA and LTL had independent effect on CAD risk in our study population.

Whole Genome Sequence Data Provides Novel Insights Into the Genetic Architecture of Meat Quality Traits in Beef

Tenderness is a major quality attribute for fresh beef steaks in the United States, and meat quality traits in general are suitable candidates for genomic research. The objectives of the present analysis were to (1) perform genome-wide association (GWA) analysis for marbling, Warner-Bratzler shear force (WBSF), tenderness, and connective tissue using whole-genome data in an Angus population, (2) identify enriched pathways in each GWA analysis; (3) construct a protein-protein interaction network using the associated genes and (4) perform a μ-calpain proteolysis assessment for associated structural proteins. An Angus-sired population of 2,285 individuals was assessed. Animals were transported to a commercial packing plant and harvested at an average age of 457 ± 46 days. After 48 h postmortem, marbling was recorded by graders' visual appraisal. Two 2.54-cm steaks were sampled from each muscle for recording of WBSF, and tenderness, and connective tissue by a sensory panel. The relevance of additive effects on marbling, WBSF, tenderness, and connective tissue was evaluated on a genome-wide scale using a two-step mixed model-based approach in single-trait analysis. A tissue-restricted gene enrichment was performed for each GWA where all polymorphisms with an association p-value lower than 1 × 10-3 were included. The genes identified as associated were included in a protein-protein interaction network and a candidate structural protein assessment of proteolysis analyses. A total of 1,867, 3,181, 3,926, and 3,678 polymorphisms were significantly associated with marbling, WBSF, tenderness, and connective tissue, respectively. The associate region on BTA29 (36,432,655-44,313,046 bp) harbors 13 highly significant markers for meat quality traits. Enrichment for the GO term GO:0005634 (Nucleus), which includes transcription factors, was evident. The final protein-protein network included 431 interations between 349 genes. The 42 most important genes based on significance that encode structural proteins were included in a proteolysis analysis, and 81% of these proteins were potential μ-Calpain substrates. Overall, this comprehensive study unraveled genetic variants, genes and mechanisms of action responsible for the variation in meat quality traits. Our findings can provide opportunities for improving meat quality in beef cattle via marker-assisted selection.

Identifying blood pressure loci whose effects are modulated by multiple lifestyle exposures

Although multiple lifestyle exposures simultaneously impact blood pressure (BP) and cardiovascular health, most analysis so far has considered each single lifestyle exposure (e.g., smoking) at a time. Here, we exploit gene-multiple lifestyle exposure interactions to find novel BP loci. For each of 6,254 Framingham Heart Study participants, we computed lifestyle risk score (LRS) value by aggregating the risk of four lifestyle exposures (smoking, alcohol, education, and physical activity) on BP. Using the LRS, we performed genome-wide gene-environment interaction analysis in systolic and diastolic BP using the joint 2 degree of freedom (DF) and 1 DF interaction tests. We identified one genome-wide significant (p < 5 × 10-8 ) and 11 suggestive (p < 1 × 10-6 ) loci. Gene-environment analysis using single lifestyle exposures identified only one of the 12 loci. Nine of the 12 BP loci detected were novel. Loci detected by the LRS were located within or nearby genes with biologically plausible roles in the pathophysiology of hypertension, including KALRN, VIPR2, SNX1, and DAPK2. Our results suggest that simultaneous consideration of multiple lifestyle exposures in gene-environment interaction analysis can identify additional loci missed by single lifestyle approaches.

High-throughput and efficient multilocus genome-wide association study on longitudinal outcomes

MOTIVATION

With the emerging of high-dimensional genomic data, genetic analysis such as genome-wide association studies (GWAS) have played an important role in identifying disease-related genetic variants and novel treatments. Complex longitudinal phenotypes are commonly collected in medical studies. However, since limited analytical approaches are available for longitudinal traits, these data are often underutilized. In this article, we develop a high-throughput machine learning approach for multilocus GWAS using longitudinal traits by coupling Empirical Bayesian Estimates from mixed-effects modeling with a novel ℓ0-norm algorithm.

RESULTS

Extensive simulations demonstrated that the proposed approach not only provided accurate selection of single nucleotide polymorphisms (SNPs) with comparable or higher power but also robust control of false positives. More importantly, this novel approach is highly scalable and could be approximately >1000 times faster than recently published approaches, making genome-wide multilocus analysis of longitudinal traits possible. In addition, our proposed approach can simultaneously analyze millions of SNPs if the computer memory allows, thereby potentially allowing a true multilocus analysis for high-dimensional genomic data. With application to the data from Alzheimer's Disease Neuroimaging Initiative, we confirmed that our approach can identify well-known SNPs associated with AD and were much faster than recently published approaches (≥6000 times).

AVAILABILITY AND IMPLEMENTATION

The source code and the testing datasets are available at https://github.com/Myuan2019/EBE_APML0.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide association study of Buruli ulcer in rural Benin highlights role of two LncRNAs and the autophagy pathway

Buruli ulcer, caused by Mycobacterium ulcerans and characterized by devastating necrotizing skin lesions, is the third mycobacterial disease worldwide. The role of host genetics in susceptibility to Buruli ulcer has long been suggested. We conduct the first genome-wide association study of Buruli ulcer on a sample of 1524 well characterized patients and controls from rural Benin. Two-stage analyses identify two variants located within LncRNA genes: rs9814705 in ENSG00000240095.1 (P = 2.85 × 10⁻⁷; odds ratio = 1.80 [1.43–2.27]), and rs76647377 in LINC01622 (P = 9.85 × 10⁻⁸; hazard ratio = 0.41 [0.28–0.60]). Furthermore, we replicate the protective effect of allele G of a missense variant located in ATG16L1, previously shown to decrease bacterial autophagy (rs2241880, P = 0.003; odds ratio = 0.31 [0.14–0.68]). Our results suggest LncRNAs and the autophagy pathway as critical factors in the development of Buruli ulcer.

Jeremy Manry, Quentin Vincent et al. report a genome-wide association study for susceptibility to Buruli ulcer in a rural population from the West African country of Benin. They identify two independently associated variants within LncRNA genes and confirm the protective effect of a missense variant in the bacterial autophagy gene ATG16L1.

Identifying potential causal effects of age at menarche: a Mendelian randomization phenome-wide association study

BACKGROUND

Age at menarche has been associated with various health outcomes. We aimed to identify potential causal effects of age at menarche on health-related traits in a hypothesis-free manner.

METHODS

We conducted a Mendelian randomization phenome-wide association study (MR-pheWAS) of age at menarche with 17,893 health-related traits in UK Biobank (n = 181,318) using PHESANT. The exposure of interest was the genetic risk score for age at menarche. We conducted a second MR-pheWAS after excluding SNPs associated with BMI from the genetic risk score, to examine whether results might be due to the genetic overlap between age at menarche and BMI. We followed up a subset of health-related traits to investigate MR assumptions and seek replication in independent study populations.

RESULTS

Of the 17,893 tests performed in our MR-pheWAS, we identified 619 associations with the genetic risk score for age at menarche at a 5% false discovery rate threshold, of which 295 were below a Bonferroni-corrected P value threshold. These included potential effects of younger age at menarche on lower lung function, higher heel bone-mineral density, greater burden of psychosocial/mental health problems, younger age at first birth, higher risk of childhood sexual abuse, poorer cardiometabolic health, and lower physical activity. After exclusion of variants associated with BMI, the genetic risk score for age at menarche was related to 37 traits at a 5% false discovery rate, of which 29 were below a Bonferroni-corrected P value threshold. We attempted to replicate findings for bone-mineral density, lung function, neuroticism, and childhood sexual abuse using 5 independent cohorts/consortia. While estimates for lung function, higher bone-mineral density, neuroticism, and childhood sexual abuse in replication cohorts were consistent with UK Biobank estimates, confidence intervals were wide and often included the null.

CONCLUSIONS

The genetic risk score for age at menarche was related to a broad range of health-related traits. Follow-up analyses indicated imprecise evidence of an effect of younger age at menarche on greater bone-mineral density, lower lung function, higher neuroticism score, and greater risk of childhood sexual abuse in the smaller replication samples available; hence, these findings need further exploration when larger independent samples become available.

Gene mapping and genomic prediction of bull fertility using sex chromosome markers

Service sire has been recognized as an important factor affecting dairy herd fertility. Our group has reported promising results on gene mapping and genomic prediction of dairy bull fertility using autosomal SNP markers. Little is known, however, about the genetic contribution of sex chromosomes, which are enriched in genes related to sexual development and reproduction. As such, the main goal of this study was to investigate the effect of SNP markers on X and Y chromosomes (BTAX and BTAY, respectively) on sire conception rate (SCR) in US Holstein bulls. The analysis included a total of 5,014 bulls with SCR records and genotypes for roughly 291k SNP located on the autosomes, 1.5k SNP located on the pseudoautosomal region (PAR), 13.7k BTAX-specific SNP, and 24 BTAY-specific SNP. We first performed genomic scans of the sex chromosomes, and then we evaluated the genomic prediction of SCR including BTAX SNP markers in the predictive models. Two markers located on PAR and 3 markers located on the X-specific region showed significant associations with sire fertility. Interestingly, these regions harbor genes, such as FAM9B, TBL1X, and PIH1D3, that are directly implicated in testosterone concentration, spermatogenesis, and sperm motility. On the other hand, BTAY showed very low genetic variability, and none of the segregating markers were associated with SCR. Notably, model predictive ability was largely improved by including BTAX markers. Indeed, the combination of autosomal with BTAX SNP delivered predictive correlations around 0.343, representing an increase in accuracy of about 7.5% compared with the standard whole autosomal genome approach. Overall, this study provides evidence of the importance of both PAR and X-specific regions in male fertility in dairy cattle. These findings may help to improve conception rates in dairy herds through accurate genome-guided decisions on bull fertility.

Radiogenomics Consortium Genome-Wide Association Study Meta-Analysis of Late Toxicity After Prostate Cancer Radiotherapy

BACKGROUND

A total of 10%-20% of patients develop long-term toxicity following radiotherapy for prostate cancer. Identification of common genetic variants associated with susceptibility to radiotoxicity might improve risk prediction and inform functional mechanistic studies.

METHODS

We conducted an individual patient data meta-analysis of six genome-wide association studies (n = 3871) in men of European ancestry who underwent radiotherapy for prostate cancer. Radiotoxicities (increased urinary frequency, decreased urinary stream, hematuria, rectal bleeding) were graded prospectively. We used grouped relative risk models to test associations with approximately 6 million genotyped or imputed variants (time to first grade 2 or higher toxicity event). Variants with two-sided Pmeta less than 5 × 10-8 were considered statistically significant. Bayesian false discovery probability provided an additional measure of confidence. Statistically significant variants were evaluated in three Japanese cohorts (n = 962). All statistical tests were two-sided.

RESULTS

Meta-analysis of the European ancestry cohorts identified three genomic signals: single nucleotide polymorphism rs17055178 with rectal bleeding (Pmeta = 6.2 × 10-10), rs10969913 with decreased urinary stream (Pmeta = 2.9 × 10-10), and rs11122573 with hematuria (Pmeta = 1.8 × 10-8). Fine-scale mapping of these three regions was used to identify another independent signal (rs147121532) associated with hematuria (Pconditional = 4.7 × 10-6). Credible causal variants at these four signals lie in gene-regulatory regions, some modulating expression of nearby genes. Previously identified variants showed consistent associations (rs17599026 with increased urinary frequency, rs7720298 with decreased urinary stream, rs1801516 with overall toxicity) in new cohorts. rs10969913 and rs17599026 had similar effects in the photon-treated Japanese cohorts.

CONCLUSIONS

This study increases the understanding of the architecture of common genetic variants affecting radiotoxicity, points to novel radio-pathogenic mechanisms, and develops risk models for testing in clinical studies. Further multinational radiogenomics studies in larger cohorts are worthwhile.

Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure

Heart failure (HF) is a leading cause of morbidity and mortality worldwide. A small proportion of HF cases are attributable to monogenic cardiomyopathies and existing genome-wide association studies (GWAS) have yielded only limited insights, leaving the observed heritability of HF largely unexplained. We report results from a GWAS meta-analysis of HF comprising 47,309 cases and 930,014 controls. Twelve independent variants at 11 genomic loci are associated with HF, all of which demonstrate one or more associations with coronary artery disease (CAD), atrial fibrillation, or reduced left ventricular function, suggesting shared genetic aetiology. Functional analysis of non-CAD-associated loci implicate genes involved in cardiac development (MYOZ1, SYNPO2L), protein homoeostasis (BAG3), and cellular senescence (CDKN1A). Mendelian randomisation analysis supports causal roles for several HF risk factors, and demonstrates CAD-independent effects for atrial fibrillation, body mass index, and hypertension. These findings extend our knowledge of the pathways underlying HF and may inform new therapeutic strategies.

In vivo single-cell lineage tracing in zebrafish using high-resolution infrared laser-mediated gene induction microscopy

Heterogeneity broadly exists in various cell types both during development and at homeostasis. Investigating heterogeneity is crucial for comprehensively understanding the complexity of ontogeny, dynamics, and function of specific cell types. Traditional bulk-labeling techniques are incompetent to dissect heterogeneity within cell population, while the new single-cell lineage tracing methodologies invented in the last decade can hardly achieve high-fidelity single-cell labeling and long-term in-vivo observation simultaneously. In this work, we developed a high-precision infrared laser-evoked gene operator heat-shock system, which uses laser-induced CreERT2 combined with loxP-DsRedx-loxP-GFP reporter to achieve precise single-cell labeling and tracing. In vivo study indicated that this system can precisely label single cell in brain, muscle and hematopoietic system in zebrafish embryo. Using this system, we traced the hematopoietic potential of hemogenic endothelium (HE) in the posterior blood island (PBI) of zebrafish embryo and found that HEs in the PBI are heterogeneous, which contains at least myeloid unipotent and myeloid-lymphoid bipotent subtypes.

The impact of disregarding family structure on genome-wide association analysis of complex diseases in cohorts with simple pedigrees

The generalized linear mixed models (GLMMs) methodology is the standard framework for genome-wide association studies (GWAS) of complex diseases in family-based cohorts. Fitting GLMMs in very large cohorts, however, can be computationally demanding. Also, the modified versions of GLMM using faster algorithms may underperform, for instance when a single nucleotide polymorphism (SNP) is correlated with fixed-effects covariates. We investigated the extent to which disregarding family structure may compromise GWAS in cohorts with simple pedigrees by contrasting logistic regression models (i.e., with no family structure) to three LMMs-based ones. Our analyses showed that the logistic regression models in general resulted in smaller P values compared with the LMMs-based models; however, the differences in P values were mostly minor. Disregarding family structure had little impact on determining disease-associated SNPs at genome-wide level of significance (i.e., P < 5E-08) as the four P values resulted from the tested methods for any SNP were all below or all above 5E-08. Nevertheless, larger discrepancies were detected between logistic regression and LMMs-based models at suggestive level of significance (i.e., of 5E-08 ≤ P < 5E-06). The SNP effects estimated by the logistic regression models were not statistically different from those estimated by GLMMs that implemented Wald's test. However, several SNP effects were significantly different from their counterparts in LMMs analyses. We suggest that fitting GLMMs with Wald's test on a pre-selected subset of SNPs obtained from logistic regression models can ensure the balance between the speed of analyses and the accuracy of parameters.

Multi-ancestry sleep-by-SNP interaction analysis in 126,926 individuals reveals lipid loci stratified by sleep duration

Both short and long sleep are associated with an adverse lipid profile, likely through different biological pathways. To elucidate the biology of sleep-associated adverse lipid profile, we conduct multi-ancestry genome-wide sleep-SNP interaction analyses on three lipid traits (HDL-c, LDL-c and triglycerides). In the total study sample (discovery + replication) of 126,926 individuals from 5 different ancestry groups, when considering either long or short total sleep time interactions in joint analyses, we identify 49 previously unreported lipid loci, and 10 additional previously unreported lipid loci in a restricted sample of European-ancestry cohorts. In addition, we identify new gene-sleep interactions for known lipid loci such as LPL and PCSK9. The previously unreported lipid loci have a modest explained variance in lipid levels: most notable, gene-short-sleep interactions explain 4.25% of the variance in triglyceride level. Collectively, these findings contribute to our understanding of the biological mechanisms involved in sleep-associated adverse lipid profiles.

Role of Immune Response, Inflammation, and Tumor Immune Response-Related Cytokines/Chemokines in Melanoma Progression

To investigate the role of tumor cytokines/chemokines in melanoma immune response, we estimated the proportions of immune cell subsets in melanoma tumors from The Cancer Genome Atlas, followed by evaluation of the association between cytokine/chemokine expression and these subsets. We then investigated the association of immune cell subsets, chemokines, and cytokines with patient survival. Finally, we evaluated the immune cell tumor-infiltrating lymphocyte (TIL) score for correlation with melanoma patient outcome in a separate cohort. There was good agreement between RNA sequencing estimation of T-cell subset and pathologist-determined TIL score. Expression levels of cytokines IL-12A, IFNG, and IL-10, and chemokines CXCL9 and CXCL10 were positively correlated with PDCD1, CTLA-4, and CD8+ T-cell subset, but negatively correlated with tumor purity (Bonferroni-corrected P < 0.05). In multivariable analysis, higher expression levels of cytokines IFN-γ and TGFB1, but not chemokines, were associated with improved overall survival. A higher expression level of CD8+ T-cell subset was also associated with improved overall survival (hazard ratio [HR] = 0.06, 95% confidence interval [CI] = 0.01-0.35, P = 0.002). Finally, multivariable analysis showed that patients with a brisk TIL score had improved melanoma-specific survival than those with a nonbrisk score (HR = 0.51, 95% CI = 0.27-0.98, P = 0.0423). These results suggest that the expression of specific tumor cytokines represents important biomarkers of melanoma immune response.

Novel Genetic Locus Influencing Retinal Venular Tortuosity Is Also Associated With Risk of Coronary Artery Disease

Supplemental Digital Content is available in the text.

The retina may provide readily accessible imaging biomarkers of global cardiovascular health. Increasing evidence suggests variation in retinal vascular traits is highly heritable. This study aimed to identify the genetic determinants of retinal vascular traits.

Target genes, variants, tissues and transcriptional pathways influencing human serum urate levels

Elevated serum urate levels cause gout and correlate with cardiometabolic diseases via poorly understood mechanisms. We performed a trans-ancestry genome-wide association study of serum urate in 457,690 individuals, identifying 183 loci (147 previously unknown) that improve the prediction of gout in an independent cohort of 334,880 individuals. Serum urate showed significant genetic correlations with many cardiometabolic traits, with genetic causality analyses supporting a substantial role for pleiotropy. Enrichment analysis, fine-mapping of urate-associated loci and colocalization with gene expression in 47 tissues implicated the kidney and liver as the main target organs and prioritized potentially causal genes and variants, including the transcriptional master regulators in the liver and kidney, HNF1A and HNF4A. Experimental validation showed that HNF4A transactivated the promoter of ABCG2, encoding a major urate transporter, in kidney cells, and that HNF4A p.Thr139Ile is a functional variant. Transcriptional coregulation within and across organs may be a general mechanism underlying the observed pleiotropy between urate and cardiometabolic traits.

Evidence of a Causal Association Between Cancer and Alzheimer's Disease: a Mendelian Randomization Analysis

While limited observational evidence suggests that cancer survivors have a decreased risk of developing Alzheimer's disease (AD), and vice versa, it is not clear whether this relationship is causal. Using a Mendelian randomization approach that provides evidence of causality, we found that genetically predicted lung cancer (OR 0.91, 95% CI 0.84-0.99, p = 0.019), leukemia (OR 0.98, 95% CI 0.96-0.995, p = 0.012), and breast cancer (OR 0.94, 95% CI 0.89-0.99, p = 0.028) were associated with 9.0%, 2.4%, and 5.9% lower odds of AD, respectively, per 1-unit higher log odds of cancer. When genetic predictors of all cancers were pooled, cancer was associated with 2.5% lower odds of AD (OR 0.98, 95% CI 0.96-0.988, p = 0.00027) per 1-unit higher log odds of cancer. Finally, genetically predicted smoking-related cancers showed a more robust inverse association with AD than non-smoking related cancers (OR 0.95, 95% CI 0.92-0.98, p = 0.0026, vs. OR 0.98, 95% CI 0.97-0.995, p = 0.0091).