Interactive molecular causal networks of hypertension using a fast machine learning algorithm MRdualPC.

BACKGROUND: Understanding the complex interactions between genes and their causal effects on diseases is crucial for developing targeted treatments and gaining insight into biological mechanisms. However, the analysis of molecular networks, especially in the context of high-dimensional data, presents significant challenges. METHODS: This study introduces MRdualPC, a computationally tractable algorithm based on the MRPC approach, to infer large-scale causal molecular networks. We apply MRdualPC to investigate the upstream causal transcriptomics influencing hypertension using a comprehensive dataset of kidney genome and transcriptome data. RESULTS: Our algorithm proves to be 100 times faster than MRPC on average in identifying transcriptomics drivers of hypertension. Through clustering, we identify 63 modules with causal driver genes, including 17 modules with extensive causal networks. Notably, we find that genes within one of the causal networks are associated with the electron transport chain and oxidative phosphorylation, previously linked to hypertension. Moreover, the identified causal ancestor genes show an over-representation of blood pressure-related genes. CONCLUSIONS: MRdualPC has the potential for broader applications beyond gene expression data, including multi-omics integration. While there are limitations, such as the need for clustering in large gene expression datasets, our study represents a significant advancement in building causal molecular networks, offering researchers a valuable tool for analyzing big data and investigating complex diseases.

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Heritability and family-based GWAS analyses of the N-acyl ethanolamine and ceramide plasma lipidome.

Signalling lipids of the N-acyl ethanolamine (NAE) and ceramide (CER) classes have emerged as potential biomarkers of cardiovascular disease (CVD). We sought to establish the heritability of plasma NAEs (including the endocannabinoid anandamide) and CERs, to identify common DNA variants influencing the circulating concentrations of the heritable lipids, and assess causality of these lipids in CVD using 2-sample Mendelian randomization (2SMR). Nine NAEs and 16 CERs were analyzed in plasma samples from 999 members of 196 British Caucasian families, using targeted ultra-performance liquid chromatography with tandem mass spectrometry. All lipids were significantly heritable (h2 = 36-62%). A missense variant (rs324420) in the gene encoding the enzyme fatty acid amide hydrolase (FAAH), which degrades NAEs, associated at genome-wide association study (GWAS) significance (P < 5 × 10-8) with four NAEs (DHEA, PEA, LEA and VEA). For CERs, rs680379 in the SPTLC3 gene, which encodes a subunit of the rate-limiting enzyme in CER biosynthesis, associated with a range of species (e.g. CER[N(24)S(19)]; P = 4.82 × 10-27). We observed three novel associations between SNPs at the CD83, SGPP1 and DEGS1 loci, and plasma CER traits (P < 5 × 10-8). 2SMR in the CARDIoGRAMplusC4D cohorts (60 801 cases; 123 504 controls) and in the DIAGRAM cohort (26 488 cases; 83 964 controls), using the genetic instruments from our family-based GWAS, did not reveal association between genetically determined differences in CER levels and CVD or diabetes. Two of the novel GWAS loci, SGPP1 and DEGS1, suggested a casual association between CERs and a range of haematological phenotypes, through 2SMR in the UK Biobank, INTERVAL and UKBiLEVE cohorts (n = 110 000-350 000).

Rare GATA6 variants associated with risk of congenital heart disease phenotypes in 200,000 UK Biobank exomes.

Congenital heart disease (CHD) has a complex and largely uncharacterised genetic etiology. Using 200,000 UK Biobank (UKB) exomes, we assess the burden of ultra-rare, potentially pathogenic variants in the largest case/control cohort of predominantly mild CHD to date. We find an association with GATA6, a member of the GATA family of transcription factors that play an important role during heart development and has been linked with several CHD phenotypes previously. Several identified GATA6 variants are previously unreported and their roles in conferring risk to CHD warrants further study. We demonstrate that despite limitations regarding detailed familial phenotype information in large-scale biobank projects, through careful consideration of case and control cohorts it is possible to derive important associations.

Joint analysis of functionally related genes yields further candidates associated with Tetralogy of Fallot.

Although several genes involved in the development of Tetralogy of Fallot have been identified, no genetic diagnosis is available for the majority of patients. Low statistical power may have prevented the identification of further causative genes in gene-by-gene survey analyses. Thus, bigger samples and/or novel analytic approaches may be necessary. We studied if a joint analysis of groups of functionally related genes might be a useful alternative approach. Our reanalysis of whole-exome sequencing data identified 12 groups of genes that exceedingly contribute to the burden of Tetralogy of Fallot. Further analysis of those groups showed that genes with high-impact variants tend to interact with each other. Thus, our results strongly suggest that additional candidate genes may be found by studying the protein interaction network of known causative genes. Moreover, our results show that the joint analysis of functionally related genes can be a useful complementary approach to classical single-gene analyses.

Data-independent acquisition mass spectrometry in severe rheumatic heart disease (RHD) identifies a proteomic signature showing ongoing inflammation and effectively classifying RHD cases.

BACKGROUND: Rheumatic heart disease (RHD) remains a major source of morbidity and mortality in developing countries. A deeper insight into the pathogenetic mechanisms underlying RHD could provide opportunities for drug repurposing, guide recommendations for secondary penicillin prophylaxis, and/or inform development of near-patient diagnostics. METHODS: We performed quantitative proteomics using Sequential Windowed Acquisition of All Theoretical Fragment Ion Mass Spectrometry (SWATH-MS) to screen protein expression in 215 African patients with severe RHD, and 230 controls. We applied a machine learning (ML) approach to feature selection among the 366 proteins quantifiable in at least 40% of samples, using the Boruta wrapper algorithm. The case-control differences and contribution to Area Under the Receiver Operating Curve (AUC) for each of the 56 proteins identified by the Boruta algorithm were calculated by Logistic Regression adjusted for age, sex and BMI. Biological pathways and functions enriched for proteins were identified using ClueGo pathway analyses. RESULTS: Adiponectin, complement component C7 and fibulin-1, a component of heart valve matrix, were significantly higher in cases when compared with controls. Ficolin-3, a protein with calcium-independent lectin activity that activates the complement pathway, was lower in cases than controls. The top six biomarkers from the Boruta analyses conferred an AUC of 0.90 indicating excellent discriminatory capacity between RHD cases and controls. CONCLUSIONS: These results support the presence of an ongoing inflammatory response in RHD, at a time when severe valve disease has developed, and distant from previous episodes of acute rheumatic fever. This biomarker signature could have potential utility in recognizing different degrees of ongoing inflammation in RHD patients, which may, in turn, be related to prognostic severity.

Genetic analyses of circulating PUFA-derived mediators identifies heritable dihydroxyeicosatrienoic acid species.

Estimates of heritability are the first step in identifying a trait with substantial variation due to genetic factors. Large-scale genetic analyses can identify the DNA variants that influence the levels of circulating lipid species and the statistical technique Mendelian randomisation can use these DNA variants to address potential causality of these lipids in disease. We estimated the heritability of plasma eicosanoids, octadecanoids and docosanoids to identify those lipid species with substantial heritability. We analysed plasma lipid mediators in 31 White British families (196 participants) ascertained for high blood pressure and deeply clinically and biochemically phenotyped over a 25-year period. We found that the dihydroxyeicosatrienoic acid (DHET) species, 11,12-DHET and 14,15-DHET, products of arachidonic acid metabolism by cytochrome P450 (CYP) monooxygenase and soluble epoxide hydrolase (sEH), exhibited substantial heritability (h2 = 33%-37%; Padj<0.05). Identification of these two heritable bioactive lipid species allows for future large-scale, targeted, lipidomics-genomics analyses to address causality in cardiovascular and other diseases.

The use of missing values in proteomic data-independent acquisition mass spectrometry to enable disease activity discrimination.

MOTIVATION: Data-independent acquisition mass spectrometry allows for comprehensive peptide detection and relative quantification than standard data-dependent approaches. While less prone to missing values, these still exist. Current approaches for handling the so-called missingness have challenges. We hypothesized that non-random missingness is a useful biological measure and demonstrate the importance of analysing missingness for proteomic discovery within a longitudinal study of disease activity. RESULTS: The magnitude of missingness did not correlate with mean peptide concentration. The magnitude of missingness for each protein strongly correlated between collection time points (baseline, 3 months, 6 months; R = 0.95-0.97, confidence interval = 0.94-0.97) indicating little time-dependent effect. This allowed for the identification of proteins with outlier levels of missingness that differentiate between the patient groups characterized by different patterns of disease activity. The association of these proteins with disease activity was confirmed by machine learning techniques. Our novel approach complements analyses on complete observations and other missing value strategies in biomarker prediction of disease activity. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Rare variants in KDR, encoding VEGF Receptor 2, are associated with tetralogy of Fallot.

PURPOSE: Rare genetic variants in KDR, encoding the vascular endothelial growth factor receptor 2 (VEGFR2), have been reported in patients with tetralogy of Fallot (TOF). However, their role in disease causality and pathogenesis remains unclear. METHODS: We conducted exome sequencing in a familial case of TOF and large-scale genetic studies, including burden testing, in >1,500 patients with TOF. We studied gene-targeted mice and conducted cell-based assays to explore the role of KDR genetic variation in the etiology of TOF. RESULTS: Exome sequencing in a family with two siblings affected by TOF revealed biallelic missense variants in KDR. Studies in knock-in mice and in HEK 293T cells identified embryonic lethality for one variant when occurring in the homozygous state, and a significantly reduced VEGFR2 phosphorylation for both variants. Rare variant burden analysis conducted in a set of 1,569 patients of European descent with TOF identified a 46-fold enrichment of protein-truncating variants (PTVs) in TOF cases compared to controls (P = 7 × 10-11). CONCLUSION: Rare KDR variants, in particular PTVs, strongly associate with TOF, likely in the setting of different inheritance patterns. Supported by genetic and in vivo and in vitro functional analysis, we propose loss-of-function of VEGFR2 as one of the mechanisms involved in the pathogenesis of TOF.

Whole Exome Sequencing Reveals the Major Genetic Contributors to Nonsyndromic Tetralogy of Fallot.

RATIONALE: Familial recurrence studies provide strong evidence for a genetic component to the predisposition to sporadic, nonsyndromic Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease phenotype. Rare genetic variants have been identified as important contributors to the risk of congenital heart disease, but relatively small numbers of TOF cases have been studied to date. OBJECTIVE: We used whole exome sequencing to assess the prevalence of unique, deleterious variants in the largest cohort of nonsyndromic TOF patients reported to date. METHODS AND RESULTS: Eight hundred twenty-nine TOF patients underwent whole exome sequencing. The presence of unique, deleterious variants was determined; defined by their absence in the Genome Aggregation Database and a scaled combined annotation-dependent depletion score of ≥20. The clustering of variants in 2 genes, NOTCH1 and FLT4, surpassed thresholds for genome-wide significance (assigned as P<5×10-8) after correction for multiple comparisons. NOTCH1 was most frequently found to harbor unique, deleterious variants. Thirty-one changes were observed in 37 probands (4.5%; 95% CI, 3.2%-6.1%) and included 7 loss-of-function variants 22 missense variants and 2 in-frame indels. Sanger sequencing of the unaffected parents of 7 cases identified 5 de novo variants. Three NOTCH1 variants (p.G200R, p.C607Y, and p.N1875S) were subjected to functional evaluation, and 2 showed a reduction in Jagged1-induced NOTCH signaling. FLT4 variants were found in 2.4% (95% CI, 1.6%-3.8%) of TOF patients, with 21 patients harboring 22 unique, deleterious variants. The variants identified were distinct to those that cause the congenital lymphoedema syndrome Milroy disease. In addition to NOTCH1, FLT4 and the well-established TOF gene, TBX1, we identified potential association with variants in several other candidates, including RYR1, ZFPM1, CAMTA2, DLX6, and PCM1. CONCLUSIONS: The NOTCH1 locus is the most frequent site of genetic variants predisposing to nonsyndromic TOF, followed by FLT4. Together, variants in these genes are found in almost 7% of TOF patients.

Hypertension and renin-angiotensin system blockers are not associated with expression of angiotensin-converting enzyme 2 (ACE2) in the kidney.

AIMS: Angiotensin-converting enzyme 2 (ACE2) is the cellular entry point for severe acute respiratory syndrome coronavirus (SARS-CoV-2)-the cause of coronavirus disease 2019 (COVID-19). However, the effect of renin-angiotensin system (RAS)-inhibition on ACE2 expression in human tissues of key relevance to blood pressure regulation and COVID-19 infection has not previously been reported. METHODS AND RESULTS: We examined how hypertension, its major metabolic co-phenotypes, and antihypertensive medications relate to ACE2 renal expression using information from up to 436 patients whose kidney transcriptomes were characterized by RNA-sequencing. We further validated some of the key observations in other human tissues and/or a controlled experimental model. Our data reveal increasing expression of ACE2 with age in both human lungs and the kidney. We show no association between renal expression of ACE2 and either hypertension or common types of RAS inhibiting drugs. We demonstrate that renal abundance of ACE2 is positively associated with a biochemical index of kidney function and show a strong enrichment for genes responsible for kidney health and disease in ACE2 co-expression analysis. CONCLUSION: Our results indicate that neither hypertension nor antihypertensive treatment is likely to alter the expression of the key entry receptor for SARS-CoV-2 in the human kidney. Our data further suggest that in the absence of SARS-CoV-2 infection, kidney ACE2 is most likely nephro-protective but the age-related increase in its expression within lungs and kidneys may be relevant to the risk of SARS-CoV-2 infection.

Human Y Chromosome Exerts Pleiotropic Effects on Susceptibility to Atherosclerosis.

OBJECTIVE: The male-specific region of the Y chromosome (MSY) remains one of the most unexplored regions of the genome. We sought to examine how the genetic variants of the MSY influence male susceptibility to coronary artery disease (CAD) and atherosclerosis. Approach and Results: Analysis of 129 133 men from UK Biobank revealed that only one of 7 common MSY haplogroups (haplogroup I1) was associated with CAD-carriers of haplogroup I1 had ≈11% increase in risk of CAD when compared with all other haplogroups combined (odds ratio, 1.11; 95% CI, 1.04-1.18; P=6.8×10-4). Targeted MSY sequencing uncovered 235 variants exclusive to this haplogroup. The haplogroup I1-specific variants showed 2.45- and 1.56-fold respective enrichment for promoter and enhancer chromatin states, in cells/tissues relevant to atherosclerosis, when compared with other MSY variants. Gene set enrichment analysis in CAD-relevant tissues showed that haplogroup I1 was associated with changes in pathways responsible for early and late stages of atherosclerosis development including defence against pathogens, immunity, oxidative phosphorylation, mitochondrial respiration, lipids, coagulation, and extracellular matrix remodeling. UTY was the only Y chromosome gene whose blood expression was associated with haplogroup I1. Experimental reduction of UTY expression in macrophages led to changes in expression of 59 pathways (28 of which overlapped with those associated with haplogroup I1) and a significant reduction in the immune costimulatory signal. CONCLUSIONS: Haplogroup I1 is enriched for regulatory chromatin variants in numerous cells of relevance to CAD and increases cardiovascular risk through proatherosclerotic reprogramming of the transcriptome, partly through UTY.

Global prevalence of congenital heart disease in school-age children: a meta-analysis and systematic review.

BACKGROUND: Congenital heart disease (CHD) is the commonest birth defect. Studies estimating the prevalence of CHD in school-age children could therefore contribute to quantifying unmet health needs for diagnosis and treatment, particularly in lower-income countries. Data at school age are considerably sparser, and individual studies have generally been of small size. We conducted a literature-based meta-analysis to investigate global trends over a 40-year period. METHODS AND RESULTS: Studies reporting on CHD prevalence in school-age children (4-18 years old) from 1970 to 2017 were identified from PubMed, EMBASE, Web of Science and Google Scholar. According to the inclusion criteria, 42 studies including 2,638,475 children, reporting the prevalence of unrepaired CHDs (both pre-school diagnoses and first-time school-age diagnoses), and nine studies including 395,571 children, specifically reporting the prevalence of CHD first diagnosed at school ages, were included. Data were combined using random-effects models. The prevalence of unrepaired CHD in school children during the entire period of study was 3.809 (95% confidence intervals 3.075-4.621)/1000. A lower proportion of male than female school children had unrepaired CHD (OR = 0.84 [95% CI 0.74-0.95]; p = 0.001). Between 1970-1974 and 1995-1999, there was no significant change in the prevalence of unrepaired CHD at school age; subsequently there was an approximately 2.5-fold increase from 1.985 (95% CI 1.074-3.173)/1000 in 1995-1999 to 4.832 (95% CI 3.425-6.480)/1000 in 2010-2014, (p = 0.009). Among all CHD conditions, atrial septal defects and ventricular septal defects chiefly accounted for this increasing trend. The summarised prevalence (1970-2017) of CHD diagnoses first made in childhood was 1.384 (0.955, 1.891)/1000; during this time there was a fall from 2.050 [1.362, 2.877]/1000 pre-1995 to 0.848 [0.626, 1.104]/1000 in 1995-2014 (p = 0.04). CONCLUSIONS: Globally, these data show an increased prevalence of CHD (mainly mild CHD conditions) recognised at birth/infancy or early childhood, but remaining unrepaired at school-age. In parallel there has been a decrease of first-time CHD diagnoses in school-age children. These together imply a favourable shift of CHD recognition time to earlier in the life course. Despite this, substantial inequalities between higher and lower income countries remain. Increased healthcare resources for people born with CHD, particularly in poorer countries, are required.

Non-muscle myosin IIB (Myh10) is required for epicardial function and coronary vessel formation during mammalian development.

The coronary vasculature is an essential vessel network providing the blood supply to the heart. Disruptions in coronary blood flow contribute to cardiac disease, a major cause of premature death worldwide. The generation of treatments for cardiovascular disease will be aided by a deeper understanding of the developmental processes that underpin coronary vessel formation. From an ENU mutagenesis screen, we have isolated a mouse mutant displaying embryonic hydrocephalus and cardiac defects (EHC). Positional cloning and candidate gene analysis revealed that the EHC phenotype results from a point mutation in a splice donor site of the Myh10 gene, which encodes NMHC IIB. Complementation testing confirmed that the Myh10 mutation causes the EHC phenotype. Characterisation of the EHC cardiac defects revealed abnormalities in myocardial development, consistent with observations from previously generated NMHC IIB null mouse lines. Analysis of the EHC mutant hearts also identified defects in the formation of the coronary vasculature. We attribute the coronary vessel abnormalities to defective epicardial cell function, as the EHC epicardium displays an abnormal cell morphology, reduced capacity to undergo epithelial-mesenchymal transition (EMT), and impaired migration of epicardial-derived cells (EPDCs) into the myocardium. Our studies on the EHC mutant demonstrate a requirement for NMHC IIB in epicardial function and coronary vessel formation, highlighting the importance of this protein in cardiac development and ultimately, embryonic survival.

Induced pluripotent stem cell modelling of HLHS underlines the contribution of dysfunctional NOTCH signalling to impaired cardiogenesis.

Hypoplastic left heart syndrome (HLHS) is among the most severe forms of congenital heart disease. Although the consensus view is that reduced flow through the left heart during development is a key factor in the development of the condition, the molecular mechanisms leading to hypoplasia of left heart structures are unknown. We have generated induced pluripotent stem cells (iPSC) from five HLHS patients and two unaffected controls, differentiated these to cardiomyocytes and identified reproducible in vitro cellular and functional correlates of the HLHS phenotype. Our data indicate that HLHS-iPSC have a reduced ability to give rise to mesodermal, cardiac progenitors and mature cardiomyocytes and an enhanced ability to differentiate to smooth muscle cells. HLHS-iPSC-derived cardiomyocytes are characterised by a lower beating rate, disorganised sarcomeres and sarcoplasmic reticulum and a blunted response to isoprenaline. Whole exome sequencing of HLHS fibroblasts identified deleterious variants in NOTCH receptors and other genes involved in the NOTCH signalling pathway. Our data indicate that the expression of NOTCH receptors was significantly downregulated in HLHS-iPSC-derived cardiomyocytes alongside NOTCH target genes confirming downregulation of NOTCH signalling activity. Activation of NOTCH signalling via addition of Jagged peptide ligand during the differentiation of HLHS-iPSC restored their cardiomyocyte differentiation capacity and beating rate and suppressed the smooth muscle cell formation. Together, our data provide firm evidence for involvement of NOTCH signalling in HLHS pathogenesis, reveal novel genetic insights important for HLHS pathology and shed new insights into the role of this pathway during human cardiac development.

A genome-wide association study of congenital cardiovascular left-sided lesions shows association with a locus on chromosome 20.

Congenital heart defects involving left-sided lesions (LSLs) are relatively common birth defects with substantial morbidity and mortality. Previous studies have suggested a high heritability with a complex genetic architecture, such that only a few LSL loci have been identified. We performed a genome-wide case-control association study to address the role of common variants using a discovery cohort of 778 cases and 2756 controls. We identified a genome-wide significant association mapping to a 200 kb region on chromosome 20q11 [P= 1.72 × 10-8 for rs3746446; imputed Single Nucleotide Polymorphism (SNP) rs6088703 P= 3.01 × 10-9, odds ratio (OR)= 1.6 for both]. This result was supported by transmission disequilibrium analyses using a subset of 541 case families (lowest P in region= 4.51 × 10-5, OR= 1.5). Replication in a cohort of 367 LSL cases and 5159 controls showed nominal association (P= 0.03 for rs3746446) resulting in P= 9.49 × 10-9 for rs3746446 upon meta-analysis of the combined cohorts. In addition, a group of seven SNPs on chromosome 1q21.3 met threshold for suggestive association (lowest P= 9.35 × 10-7 for rs12045807). Both regions include genes involved in cardiac development-MYH7B/miR499A on chromosome 20 and CTSK, CTSS and ARNT on chromosome 1. Genome-wide heritability analysis using case-control genotyped SNPs suggested that the mean heritability of LSLs attributable to common variants is moderately high ([Formula: see text] range= 0.26-0.34) and consistent with previous assertions. These results provide evidence for the role of common variation in LSLs, proffer new genes as potential biological candidates, and give further insight to the complex genetic architecture of congenital heart disease.

Increased Power for Detection of Parent-of-Origin Effects via the Use of Haplotype Estimation.

Parent-of-origin (or imprinting) effects relate to the situation in which traits are influenced by the allele inherited from only one parent and the allele from the other parent has little or no effect. Given SNP genotype data from case-parent trios, the parent of origin of each allele in the offspring can often be deduced unambiguously; however, this is not true when all three individuals are heterozygous. Most existing methods for investigating parent-of-origin effects operate on a SNP-by-SNP basis and either perform some sort of averaging over the possible parental transmissions or else discard ambiguous trios. If the correct parent of origin at a SNP could be determined, this would provide extra information and increase the power for detecting the effects of imprinting. We propose making use of the surrounding SNP information, via haplotype estimation, to improve estimation of parent of origin at a test SNP for case-parent trios, case-mother duos, and case-father duos. This extra information is then used in a multinomial modeling approach for estimating parent-of-origin effects at the test SNP. We show through computer simulations that our approach has increased power over previous approaches, particularly when the data consist only of duos. We apply our method to two real datasets and find a decrease in significance of p values in genomic regions previously thought to possibly harbor imprinting effects, thus weakening the evidence that such effects actually exist in these regions, although some regions retain evidence of significant effects.