Joint Linkage and Association Analysis Using GENEHUNTER-MODSCORE with an Application to Familial Pancreatic Cancer.

INTRODUCTION: Joint linkage and association (JLA) analysis combines two disease gene mapping strategies: linkage information contained in families and association information contained in populations. Such a JLA analysis can increase mapping power, especially when the evidence for both linkage and association is low to moderate. Similarly, an association analysis based on haplotypes instead of single markers can increase mapping power when the association pattern is complex. METHODS: In this paper, we present an extension to the GENEHUNTER-MODSCORE software package that enables a JLA analysis based on haplotypes and uses information from arbitrary pedigree types and unrelated individuals. Our new JLA method is an extension of the MOD score approach for linkage analysis, which allows the estimation of trait-model and linkage disequilibrium (LD) parameters, i.e., penetrance, disease-allele frequency, and haplotype frequencies. LD is modeled between alleles at a single diallelic disease locus and up to three diallelic test markers. Linkage information is contributed by additional multi-allelic flanking markers. We investigated the statistical properties of our JLA implementation using extensive simulations, and we compared our approach to another commonly used single-marker JLA test. To demonstrate the applicability of our new method in practice, we analyzed pedigree data from the German National Case Collection for Familial Pancreatic Cancer (FaPaCa). RESULTS: Based on the simulated data, we demonstrated the validity of our JLA-MOD score analysis implementation and identified scenarios in which haplotype-based tests outperformed the single-marker test. The estimated trait-model and LD parameters were in good accordance with the simulated values. Our method outperformed another commonly used JLA single-marker test when the LD pattern was complex. The exploratory analysis of the FaPaCa families led to the identification of a promising genetic region on chromosome 22q13.33, which can serve as a starting point for future mutation analysis and molecular research in pancreatic cancer. CONCLUSION: Our newly proposed JLA-MOD score method proves to be a valuable gene mapping and characterization tool, especially when either linkage or association information alone provide insufficient power to identify the disease-causing genetic variants.

Removal of Aqueous Uranyl and Arsenate Mixtures after Reaction with Limestone, PO43-, and Ca2.

The co-occurrence of uranyl and arsenate in contaminated water caused by natural processes and mining is a concern for impacted communities, including in Native American lands in the U.S. Southwest. We investigated the simultaneous removal of aqueous uranyl and arsenate after the reaction with limestone and precipitated hydroxyapatite (HAp, Ca10(PO4)6(OH)2). In benchtop experiments with an initial pH of 3.0 and initial concentrations of 1 mM U and As, uranyl and arsenate coprecipitated in the presence of 1 g L-1 limestone. However, related experiments initiated under circumneutral pH conditions showed that uranyl and arsenate remained soluble. Upon addition of 1 mM PO43- and 3 mM Ca2+ in solution (initial concentration of 0.05 mM U and As) resulted in the rapid removal of over 97% of U via Ca-U-P precipitation. In experiments with 2 mM PO43- and 10 mM Ca2+ at pH rising from 7.0 to 11.0, aqueous concentrations of As decreased (between 30 and 98%) circa pH 9. HAp precipitation in solids was confirmed by powder X-ray diffraction and scanning electron microscopy/energy dispersive X-ray. Electron microprobe analysis indicated U was coprecipitated with Ca and P, while As was mainly immobilized through HAp adsorption. The results indicate that natural materials, such as HAp and limestone, can effectively remove uranyl and arsenate mixtures.

Interfacial Interactions of Uranium and Arsenic with Microplastics: From Field Detection to Controlled Laboratory Tests.

We studied the co-occurrence of microplastics (MPs) and metals in field sites and further investigated their interfacial interaction in controlled laboratory conditions. First, we detected MPs in freshwater co-occurring with metals in rural and urban areas in New Mexico. Automated particle counting and fluorescence microscopy indicated that particles in field samples ranged from 7 to 149 particles/L. The urban location contained the highest count of confirmed MPs, including polyester, cellophane, and rayon, as indicated by Attenuated Total Reflectance-Fourier Transform Infrared (ATR-FTIR) spectroscopy analyses. Metal analyses using inductively coupled plasma (ICP) revealed that bodies of water in a rural site affected by mining legacy contained up to 332.8 µg/L of U, while all bodies of water contained As concentrations below 11.4 µg/L. These field findings motivated experiments in laboratory conditions, reacting MPs with 0.02-0.2 mM of As or U solutions at acidic and neutral pH with poly(methyl-methacrylate), polyethylene, and polystyrene MPs. In these experiments, As did not interact with any of the MPs tested at pH 3 and pH 7, nor U with any MPs at pH 3. Experiments supplied with U and MPs at pH 7 indicated that MPs served as substrate surface for the adsorption and nucleation of U precipitates. Chemical speciation modeling and microscopy analyses (i.e., Transmission Electron Microscopy [TEM]) suggest that U precipitates resemble sodium-compreignacite and schoepite. These findings have relevant implications to further understanding the occurrence and interfacial interaction of MPs and metals in freshwater.

Kinetics of Na- and K- uranyl arsenate dissolution.

We integrated aqueous chemistry analyses with geochemical modeling to determine the kinetics of the dissolution of Na and K uranyl arsenate solids (UAs(s)) at acidic pH. Improving our understanding of how UAs(s) dissolve is essential to predict transport of U and As, such as in acid mine drainage. At pH 2, Na0.48H0.52(UO2)(AsO4)(H2O)2.5(s) (NaUAs(s)) and K0.9H0.1(UO2)(AsO4)(H2O)2.5(s) (KUAs(s)) both dissolve with a rate constant of 3.2 × 10-7 mol m-2 s-1, which is faster than analogous uranyl phosphate solids. At pH 3, NaUAs(s) (6.3 × 10-8 mol m-2 s-1) and KUAs(s) (2.0 × 10-8 mol m-2 s-1) have smaller rate constants. Steady-state aqueous concentrations of U and As are similarly reached within the first several hours of reaction progress. This study provides dissolution rate constants for UAs(s), which may be integrated into reactive transport models for risk assessment and remediation of U and As contaminated waters.

Accurate tumor purity determination is critical for the analysis of homologous recombination deficiency (HRD).

Homologous recombination deficiency (HRD) is a predictive marker for response to poly (ADP-ribose) polymerase inhibitors (PARPi) in ovarian carcinoma. HRD scores have entered routine diagnostics, but the influence of algorithms, parameters and confounders has not been analyzed comprehensively. A series of 100 poorly differentiated ovarian carcinoma samples was analyzed using whole exome sequencing (WES) and genotyping. Tumor purity was determined using conventional pathology, digital pathology, and two bioinformatic methods. HRD scores were calculated from copy number profiles determined by Sequenza and by Sclust either with or without fixed tumor purity. Tumor purity determination by digital pathology combined with a tumory purity informed variant of Sequenza served as reference method for HRD scoring. Seven tumors had deleterious mutations in BRCA1/2, 12 tumors had deleterious mutations in other homologous recombination repair (HRR) genes, 18 tumors had variants of unknown significance (VUS) in BRCA1/2 or other HRR genes, while the remaining 63 tumors had no relevant alterations. Using the reference method for HRD scoring, 68 tumors were HRD-positive. HRDsum determined by WES correlated strongly with HRDsum determined by single nucleotide polymorphism (SNP) arrays (R = 0.85). Conventional pathology systematically overestimated tumor purity by 8% compared to digital pathology. All investigated methods agreed on classifying the deleterious BRCA1/2-mutated tumors as HRD-positive, but discrepancies were observed for some of the remaining tumors. Discordant HRD classification of 11% of the tumors was observed comparing the tumor purity uninformed default of Sequenza and the reference method. In conclusion, tumor purity is a critical factor for the determination of HRD scores. Assistance by digital pathology helps to improve accuracy and imprecision of its estimation.

EGFR and MMP-9 are associated with neointimal hyperplasia in systemic-to-pulmonary shunts in children with complex cyanotic heart disease.

Systemic-to-pulmonary shunt malfunction contributes to morbidity in children with complex congenital heart disease after palliative procedure. Neointimal hyperplasia might play a role in the pathogenesis increasing risk for shunt obstruction. The aim was to evaluate the role of epidermal growth factor receptor (EGFR) and matrix-metalloproteinase 9 (MMP-9) in the formation of neointimal within shunts. Immunohistochemistry was performed with anti-EGFR and anti-MMP-9 on shunts removed at follow-up palliative or corrective procedure. Whole-genome single-nucleotide polymorphisms genotyping was performed on DNA extracted from patients´ blood samples and allele frequencies were compared between the group of patients with shunts displaying severe stenosis (≥ 40% of lumen) and the remaining group. Immunohistochemistry detected EGFR and MMP-9 in 24 of 31 shunts, located mainly in the luminal area. Cross-sectional area of EGFR and MMP-9 measured in median 0.19 mm2 (IQR 0.1-0.3 mm2) and 0.04 mm2 (IQR 0.03-0.09 mm2), respectively, and correlated positively with the area of neointimal measured on histology (r = 0.729, p < 0.001 and r = 0.0479, p = 0.018, respectively). There was a trend of inverse correlation between the dose of acetylsalicylic acid and the degree of EGFR, but not MMP-9, expression within neointima. Certain alleles in epidermal growth factor (EGF) and tissue inhibitor of metalloproteinases 1 (TIMP-1) were associated with increased stenosis and neointimal hyperplasia within shunts. EGFR and MMP-9 contribute to neointimal proliferation in SP shunts of children with complex cyanotic heart disease. SP shunts from patients carrying certain risk alleles in the genes encoding for EGF and TIMP-1 displayed increased neointima.

Investigation of Shared Genetic Risk Factors Between Parkinson's Disease and Cancers.

BACKGROUND: Epidemiological studies that examined the association between Parkinson's disease (PD) and cancers led to inconsistent results, but they face a number of methodological difficulties. OBJECTIVE: We used results from genome-wide association studies (GWASs) to study the genetic correlation between PD and different cancers to identify common genetic risk factors. METHODS: We used individual data for participants of European ancestry from the Courage-PD (Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease; PD, N = 16,519) and EPITHYR (differentiated thyroid cancer, N = 3527) consortia and summary statistics of GWASs from iPDGC (International Parkinson Disease Genomics Consortium; PD, N = 482,730), Melanoma Meta-Analysis Consortium (MMAC), Breast Cancer Association Consortium (breast cancer), the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (prostate cancer), International Lung Cancer Consortium (lung cancer), and Ovarian Cancer Association Consortium (ovarian cancer) (N comprised between 36,017 and 228,951 for cancer GWASs). We estimated the genetic correlation between PD and cancers using linkage disequilibrium score regression. We studied the association between PD and polymorphisms associated with cancers, and vice versa, using cross-phenotypes polygenic risk score (PRS) analyses. RESULTS: We confirmed a previously reported positive genetic correlation of PD with melanoma (Gcorr = 0.16 [0.04; 0.28]) and reported an additional significant positive correlation of PD with prostate cancer (Gcorr = 0.11 [0.03; 0.19]). There was a significant inverse association between the PRS for ovarian cancer and PD (odds ratio [OR] = 0.89 [0.84; 0.94]). Conversely, the PRS of PD was positively associated with breast cancer (OR = 1.08 [1.06; 1.10]) and inversely associated with ovarian cancer (OR = 0.95 [0.91; 0.99]). The association between PD and ovarian cancer was mostly driven by rs183211 located in an intron of the NSF gene (17q21.31). CONCLUSIONS: We show evidence in favor of a contribution of pleiotropic genes to the association between PD and specific cancers. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Solubility and Thermodynamic Investigation of Meta-Autunite Group Uranyl Arsenate Solids with Monovalent Cations Na and K.

We investigated the aqueous solubility and thermodynamic properties of two meta-autunite group uranyl arsenate solids (UAs). The measured solubility products (log Ksp) obtained in dissolution and precipitation experiments at equilibrium pH 2 and 3 for NaUAs and KUAs ranged from -23.50 to -22.96 and -23.87 to -23.38, respectively. The secondary phases (UO2)(H2AsO4)2(H2O)(s) and trögerite, (UO2)3(AsO4)2·12H2O(s), were identified by powder X-ray diffraction in the reacted solids of KUA precipitation experiments (pH 2) and NaUAs dissolution and precipitation experiments (pH 3), respectively. The identification of these secondary phases in reacted solids suggest that H3O+ co-occurring with Na or K in the interlayer region can influence the solubilities of uranyl arsenate solids. The standard-state enthalpy of formation from the elements (ΔHf-el) of NaUAs is -3025 ± 22 kJ mol-1 and for KUAs is -3000 ± 28 kJ mol-1 derived from measurements by drop solution calorimetry, consistent with values reported in other studies for uranyl phosphate solids. This work provides novel thermodynamic information for reactive transport models to interpret and predict the influence of uranyl arsenate solids on soluble concentrations of U and As in contaminated waters affected by mining legacy and other anthropogenic activities.

Genome-wide Association and Meta-analysis of Age at Onset in Parkinson Disease: Evidence From the COURAGE-PD Consortium.

BACKGROUND AND OBJECTIVES: Considerable heterogeneity exists in the literature concerning genetic determinants of the age at onset (AAO) of Parkinson disease (PD), which could be attributed to a lack of well-powered replication cohorts. The previous largest genome-wide association studies (GWAS) identified SNCA and TMEM175 loci on chromosome (Chr) 4 with a significant influence on the AAO of PD; these have not been independently replicated. This study aims to conduct a meta-analysis of GWAS of PD AAO and validate previously observed findings in worldwide populations. METHODS: A meta-analysis was performed on PD AAO GWAS of 30 populations of predominantly European ancestry from the Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in Parkinson's Disease (COURAGE-PD) Consortium. This was followed by combining our study with the largest publicly available European ancestry dataset compiled by the International Parkinson Disease Genomics Consortium (IPDGC). RESULTS: The COURAGE-PD Consortium included a cohort of 8,535 patients with PD (91.9%: Europeans and 9.1%: East Asians). The average AAO in the COURAGE-PD dataset was 58.9 years (SD = 11.6), with an underrepresentation of females (40.2%). The heritability estimate for AAO in COURAGE-PD was 0.083 (SE = 0.057). None of the loci reached genome-wide significance (p < 5 × 10-8). Nevertheless, the COURAGE-PD dataset confirmed the role of the previously published TMEM175 variant as a genetic determinant of the AAO of PD with Bonferroni-corrected nominal levels of significance (p < 0.025): (rs34311866: ß(SE)COURAGE = 0.477(0.203), p COURAGE = 0.0185). The subsequent meta-analysis of COURAGE-PD and IPDGC datasets (Ntotal = 25,950) led to the identification of 2 genome-wide significant association signals on Chr 4, including the previously reported SNCA locus (rs983361: ß(SE)COURAGE+IPDGC = 0.720(0.122), p COURAGE+IPDGC = 3.13 × 10-9) and a novel BST1 locus (rs4698412: ß(SE)COURAGE+IPDGC = -0.526(0.096), p COURAGE+IPDGC = 4.41 × 10-8). DISCUSSION: Our study further refines the genetic architecture of Chr 4 underlying the AAO of the PD phenotype through the identification of BST1 as a novel AAO PD locus. These findings open a new direction for the development of treatments to delay the onset of PD.

The Interaction between HLA-DRB1 and Smoking in Parkinson's Disease Revisited.

BACKGROUND: Two studies that examined the interaction between HLA-DRB1 and smoking in Parkinson's disease (PD) yielded findings in opposite directions. OBJECTIVE: To perform a large-scale independent replication of the HLA-DRB1 × smoking interaction. METHODS: We genotyped 182 single nucleotide polymorphism (SNPs) associated with smoking initiation in 12 424 cases and 9480 controls to perform a Mendelian randomization (MR) analysis in strata defined by HLA-DRB1. RESULTS: At the amino acid level, a valine at position 11 (V11) in HLA-DRB1 displayed the strongest association with PD. MR showed an inverse association between genetically predicted smoking initiation and PD only in absence of V11 (odds ratio, 0.74, 95% confidence interval, 0.59-0.93, PInteraction  = 0.028). In silico predictions of the influence of V11 and smoking-induced modifications of α-synuclein on binding affinity showed findings consistent with this interaction pattern. CONCLUSIONS: Despite being one of the most robust findings in PD research, the mechanisms underlying the inverse association between smoking and PD remain unknown. Our findings may help better understand this association. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Human skin-resident host T cells can persist long term after allogeneic stem cell transplantation and maintain recirculation potential.

Tissue-resident memory T cells (TRM) have recently emerged as crucial cellular players for host defense in a wide variety of tissues and barrier sites. Insights into the maintenance and regulatory checkpoints of human TRM cells remain scarce, especially due to the difficulties associated with tracking T cells through time and space in humans. We therefore sought to identify and characterize skin-resident T cells in humans defined by their long-term in situ lodgment. Allogeneic hematopoietic stem cell transplantation (allo-HSCT) preceded by myeloablative chemotherapy unmasked long-term sequestration of host T cell subsets in human skin despite complete donor T cell chimerism in the blood. Single-cell chimerism analysis paired with single-cell transcriptional profiling comprehensively characterized these bona fide long-term skin-resident T cells and revealed differential tissue maintenance for distinct T cell subsets, specific TRM cell markers such as galectin-3, but also tissue exit potential with retention of the transcriptomic TRM cell identity. Analysis of 26 allo-HSCT patients revealed profound interindividual variation in the tissue maintenance of host skin T cells. The long-term persistence of host skin T cells in a subset of these patients did not correlate with the development of chronic GvHD. Our data exemplify the power of exploiting a clinical situation as a proof of concept for the existence of bona fide human skin TRM cells and reveal long-term persistence of host T cells in a peripheral tissue but not in the circulation or bone marrow in a subset of allo-HSCT patients.

Dairy Intake and Parkinson's Disease: A Mendelian Randomization Study.

BACKGROUND: Previous prospective studies highlighted dairy intake as a risk factor for Parkinson's disease (PD), particularly in men. It is unclear whether this association is causal or explained by reverse causation or confounding. OBJECTIVE: The aim is to examine the association between genetically predicted dairy intake and PD using two-sample Mendelian randomization (MR). METHODS: We genotyped a well-established instrumental variable for dairy intake located in the lactase gene (rs4988235) within the Courage-PD consortium (23 studies; 9823 patients and 8376 controls of European ancestry). RESULTS: Based on a dominant model, there was an association between genetic predisposition toward higher dairy intake and PD (odds ratio [OR] per one serving per day = 1.70, 95% confidence interval = 1.12-2.60, P = 0.013) that was restricted to men (OR = 2.50 [1.37-4.56], P = 0.003; P-difference with women = 0.029). CONCLUSIONS: Using MR, our findings provide further support for a causal relationship between dairy intake and higher PD risk, not biased by confounding or reverse causation. Further studies are needed to elucidate the underlying mechanisms. © 2022 International Parkinson and Movement Disorder Society.

Mendelian Randomisation Study of Smoking, Alcohol, and Coffee Drinking in Relation to Parkinson's Disease.

BACKGROUND: Previous studies showed that lifestyle behaviors (cigarette smoking, alcohol, coffee) are inversely associated with Parkinson's disease (PD). The prodromal phase of PD raises the possibility that these associations may be explained by reverse causation. OBJECTIVE: To examine associations of lifestyle behaviors with PD using two-sample Mendelian randomisation (MR) and the potential for survival and incidence-prevalence biases. METHODS: We used summary statistics from publicly available studies to estimate the association of genetic polymorphisms with lifestyle behaviors, and from Courage-PD (7,369 cases, 7,018 controls; European ancestry) to estimate the association of these variants with PD. We used the inverse-variance weighted method to compute odds ratios (ORIVW) of PD and 95%confidence intervals (CI). Significance was determined using a Bonferroni-corrected significance threshold (p = 0.017). RESULTS: We found a significant inverse association between smoking initiation and PD (ORIVW per 1-SD increase in the prevalence of ever smoking = 0.74, 95%CI = 0.60-0.93, p = 0.009) without significant directional pleiotropy. Associations in participants ≤67 years old and cases with disease duration ≤7 years were of a similar size. No significant associations were observed for alcohol and coffee drinking. In reverse MR, genetic liability toward PD was not associated with smoking or coffee drinking but was positively associated with alcohol drinking. CONCLUSION: Our findings are in favor of an inverse association between smoking and PD that is not explained by reverse causation, confounding, and survival or incidence-prevalence biases. Genetic liability toward PD was positively associated with alcohol drinking. Conclusions on the association of alcohol and coffee drinking with PD are hampered by insufficient statistical power.

First mitochondrial genome-wide association study with metabolomics.

In the era of personalized medicine with more and more patient-specific targeted therapies being used, we need reliable, dynamic, faster and sensitive biomarkers both to track the causes of disease and to develop and evolve therapies during the course of treatment. Metabolomics recently has shown substantial evidence to support its emerging role in disease diagnosis and prognosis. Aside from biomarkers and development of therapies, it is also an important goal to understand the involvement of mitochondrial DNA (mtDNA) in metabolic regulation, aging and disease development. Somatic mutations of the mitochondrial genome are also heavily implicated in age-related disease and aging. The general hypothesis is that an alteration in the concentration of metabolite profiles (possibly conveyed by lifestyle and environmental factors) influences the increase of mutation rate in the mtDNA and thereby contributes to a range of pathophysiological alterations observed in complex diseases. We performed an inverted mitochondrial genome-wide association analysis between mitochondrial nucleotide variants (mtSNVs) and concentration of metabolites. We used 151 metabolites and the whole sequenced mitochondrial genome from 2718 individuals to identify the genetic variants associated with metabolite profiles. Because of the high coverage, next-generation sequencing-based analysis of the mitochondrial genome allows for an accurate detection of mitochondrial heteroplasmy and for the identification of variants associated with the metabolome. The strongest association was found for mt715G > A located in the MT-12SrRNA with the metabolite ratio of C2/C10:1 (P-value = 6.82*10-09, ß = 0.909). The second most significant mtSNV was found for mt3714A > G located in the MT-ND1 with the metabolite ratio of phosphatidylcholine (PC) ae C42:5/PC ae C44:5 (P-value = 1.02*10-08, ß = 3.631). A large number of significant metabolite ratios were observed involving PC aa C36:6 and the variant mt10689G > A, located in the MT-ND4L gene. These results show an important interconnection between mitochondria and metabolite concentrations. Considering that some of the significant metabolites found in this study have been previously related to complex diseases, such as neurological disorders and metabolic conditions, these associations found here might play a crucial role for further investigations of such complex diseases. Understanding the mechanisms that control human health and disease, in particular, the role of genetic predispositions and their interaction with environmental factors is a prerequisite for the development of safe and efficient therapies for complex disorders.

The rare and the common: An Austrian DRPLA family harboring the European haplotype.

Dentatorubral-pallidoluysian atrophy (DRPLA) is a CAG trinucleotide repeat expansion disorder with an autosomal-dominant mode of inheritance and very low prevalence in Europe. We herein report the clinical characteristics of the first Austrian DRPLA family. Genetic analysis revealed the presence of a common European haplotype, suggesting a founder mutation in Europe.

Simple and reliable detection of CRISPR-induced on-target effects by qgPCR and SNP genotyping.

The recent CRISPR revolution has provided researchers with powerful tools to perform genome editing in a variety of organisms. However, recent reports indicate widespread occurrence of unintended CRISPR-induced on-target effects (OnTEs) at the edited site in mice and human induced pluripotent stem cells (iPSCs) that escape standard quality controls. By altering gene expression of targeted or neighbouring genes, OnTEs can severely affect phenotypes of CRISPR-edited cells and organisms and thus lead to data misinterpretation, which can undermine the reliability of CRISPR-based studies. Here we describe a broadly applicable framework for detecting OnTEs in genome-edited cells and organisms after non-homologous end joining-mediated and homology-directed repair-mediated editing. Our protocol enables identification of OnTEs such as large deletions, large insertions, rearrangements or loss of heterozygosity (LOH). This is achieved by subjecting genomic DNA first to quantitative genotyping PCR (qgPCR), which determines the number of intact alleles at the target site using the same PCR amplicon that has been optimized for genotyping. This combination of genotyping and quantitation makes it possible to exclude clones with monoallelic OnTEs and hemizygous editing, which are often mischaracterized as correctly edited in standard Sanger sequencing. Second, occurrence of LOH around the edited locus is detected by genotyping neighbouring single-nucleotide polymorphisms (SNPs), using either a Sanger sequencing-based method or SNP microarrays. All steps are optimized to maximize simplicity and minimize cost to promote wide dissemination and applicability across the field. The entire protocol from genomic DNA extraction to OnTE exclusion can be performed in 6-9 d.

Shared genetic pathways contribute to risk of hypertrophic and dilated cardiomyopathies with opposite directions of effect.

The heart muscle diseases hypertrophic (HCM) and dilated (DCM) cardiomyopathies are leading causes of sudden death and heart failure in young, otherwise healthy, individuals. We conducted genome-wide association studies and multi-trait analyses in HCM (1,733 cases), DCM (5,521 cases) and nine left ventricular (LV) traits (19,260 UK Biobank participants with structurally normal hearts). We identified 16 loci associated with HCM, 13 with DCM and 23 with LV traits. We show strong genetic correlations between LV traits and cardiomyopathies, with opposing effects in HCM and DCM. Two-sample Mendelian randomization supports a causal association linking increased LV contractility with HCM risk. A polygenic risk score explains a significant portion of phenotypic variability in carriers of HCM-causing rare variants. Our findings thus provide evidence that polygenic risk score may account for variability in Mendelian diseases. More broadly, we provide insights into how genetic pathways may lead to distinct disorders through opposing genetic effects.

Congenital heart disease risk loci identified by genome-wide association study in European patients.

Genetic factors undoubtedly affect the development of congenital heart disease (CHD) but still remain ill defined. We sought to identify genetic risk factors associated with CHD and to accomplish a functional analysis of SNP-carrying genes. We performed a genome-wide association study (GWAS) of 4034 White patients with CHD and 8486 healthy controls. One SNP on chromosome 5q22.2 reached genome-wide significance across all CHD phenotypes and was also indicative for septal defects. One region on chromosome 20p12.1 pointing to the MACROD2 locus identified 4 highly significant SNPs in patients with transposition of the great arteries (TGA). Three highly significant risk variants on chromosome 17q21.32 within the GOSR2 locus were detected in patients with anomalies of thoracic arteries and veins (ATAV). Genetic variants associated with ATAV are suggested to influence the expression of WNT3, and the variant rs870142 related to septal defects is proposed to influence the expression of MSX1. We analyzed the expression of all 4 genes during cardiac differentiation of human and murine induced pluripotent stem cells in vitro and by single-cell RNA-Seq analyses of developing murine and human hearts. Our data show that MACROD2, GOSR2, WNT3, and MSX1 play an essential functional role in heart development at the embryonic and newborn stages.

SCN5A Mutation Type and a Genetic Risk Score Associate Variably With Brugada Syndrome Phenotype in SCN5A Families.

BACKGROUND: Brugada syndrome (BrS) is characterized by the type 1 Brugada ECG pattern. Pathogenic rare variants in SCN5A (mutations) are identified in 20% of BrS families in whom incomplete penetrance and genotype-negative phenotype-positive individuals are observed. E1784K-SCN5A is the most common SCN5A mutation identified. We determined the association of a BrS genetic risk score (BrS-GRS) and SCN5A mutation type on BrS phenotype in BrS families with SCN5A mutations. METHODS: Subjects with a spontaneous type 1 pattern or positive/negative drug challenge from cohorts harboring SCN5A mutations were recruited from 16 centers (n=312). Single nucleotide polymorphisms previously associated with BrS at genome-wide significance were studied in both cohorts: rs11708996, rs10428132, and rs9388451. An additive linear genetic model for the BrS-GRS was assumed (6 single nucleotide polymorphism risk alleles). RESULTS: In the total population (n=312), BrS-GRS ≥4 risk alleles yielded an odds ratio of 4.15 for BrS phenotype ([95% CI, 1.45-11.85]; P=0.0078). Among SCN5A-positive individuals (n=258), BrS-GRS ≥4 risk alleles yielded an odds ratio of 2.35 ([95% CI, 0.89-6.22]; P=0.0846). In SCN5A-negative relatives (n=54), BrS-GRS ≥4 alleles yielded an odds ratio of 22.29 ([95% CI, 1.84-269.30]; P=0.0146). Among E1784K-SCN5A positive family members (n=79), hosting ≥4 risk alleles gave an odds ratio=5.12 ([95% CI, 1.93-13.62]; P=0.0011). CONCLUSIONS: Common genetic variation is associated with variable expressivity of BrS phenotype in SCN5A families, explaining in part incomplete penetrance and genotype-negative phenotype-positive individuals. SCN5A mutation genotype and a BrS-GRS associate with BrS phenotype, but the strength of association varies according to presence of a SCN5A mutation and severity of loss of function.