Search

Background: The molecular genetic basis of pulmonary arterial hypertension (PAH) is heterogeneous, with at least 26 genes displaying putative evidence for disease causality. Heterozygous variants in the ATP13A3 gene were recently identified as a new cause of adult-onset PAH. However, the contribution of ATP13A3 risk alleles to child-onset PAH remains largely unexplored. Methods and results: We report three families with a novel, autosomal recessive form of childhood-onset PAH due to biallelic ATP13A3 variants. Disease onset ranged from birth to 2.5 years and was characterised by high mortality. Using genome sequencing of parent-offspring trios, we identified a homozygous missense variant in one case, which was subsequently confirmed to cosegregate with disease in an affected sibling. Independently, compound heterozygous variants in ATP13A3 were identified in two affected siblings and in an unrelated third family. The variants included three loss of function variants (two frameshift, one nonsense) and two highly conserved missense substitutions located in the catalytic phosphorylation domain. The children were largely refractory to treatment and four died in early childhood. All parents were heterozygous for the variants and asymptomatic. Conclusion: Our findings support biallelic predicted deleterious ATP13A3 variants in autosomal recessive, childhood-onset PAH, indicating likely semidominant dose-dependent inheritance for this gene.

Expression quantitative trait loci (eQTL) can provide a link between disease susceptibility variants discovered by genetic association studies and biology. To date, eQTL mapping studies have been primarily conducted in healthy individuals from population-based cohorts. Genetic effects have been known to be context-specific and vary with changing environmental stimuli. We conducted a transcriptome- and genome-wide eQTL mapping study in a cohort of patients with idiopathic or heritable pulmonary arterial hypertension (PAH) using RNA sequencing (RNAseq) data from whole blood. We sought confirmation from three published population-based eQTL studies, including the GTEx Project, and followed up potentially novel eQTL not observed in the general population. In total, we identified 2314 eQTL of which 90% were cis-acting and 75% were confirmed by at least one of the published studies. While we observed a higher GWAS trait colocalization rate among confirmed eQTL, colocalisation rate of novel eQTL reported for lung-related phenotypes was twice as high as that of confirmed eQTL. Functional enrichment analysis of genes with novel eQTL in PAH highlighted immune-related processes, a suspected contributor to PAH. These potentially novel eQTL specific to or active in PAH could be useful in understanding genetic risk factors for other diseases that share common mechanisms with PAH.

RATIONALE: Idiopathic and hereditary pulmonary arterial hypertension (PAH) are rare but comprise a genetically heterogeneous patient group. RNA-sequencing linked to the underlying genetic architecture can be used to better understand the underlying pathology by identifying key signalling pathways and stratify patients more robustly according to clinical risk. OBJECTIVES: Using a three-stage design of RNA discovery, RNA validation/model construction and model validation to define a set of PAH-associated RNAs and a single summarising RNA model score. To define genes most likely to be involved in disease development, we performed Mendelian randomisation (MR) analysis. METHODS: RNA-sequencing was performed on whole blood samples from 359 patients with idiopathic, heritable and drug-induced PAH and 72 age- and sex-matched healthy volunteers. The score was evaluated against disease severity markers including survival analysis using all-cause mortality from diagnosis. MR used known eQTL and summary statistics from a PAH GWAS. MEASUREMENTS AND MAIN RESULTS: We identified 507 genes with differential RNA expression in PAH patients compared to controls. A model of 25 RNAs was able to distinguish PAH with 87% accuracy (AUC 95% CI: 0.791-0.945) in model validation. The RNA model score was associated with disease severity and long-term survival (p=4.66x10-6) in PAH. MR detected an association between SMAD5 levels and PAH disease susceptibility (OR:0.317, 95%CI:0.129-0.776, p=0.012). CONCLUSIONS: A whole blood RNA signature of PAH, which includes RNAs relevant to disease pathogenesis, associates with disease severity and identifies patients with poor clinical outcomes. Genetic variants associated with lower SMAD5 expression may increase susceptibility to PAH.

Pulmonary arterial hypertension (PAH) is a rare disease that leads to premature death from right heart failure. It is strongly associated with elevated red cell distribution width (RDW), a correlate of several iron status biomarkers. High RDW values can signal early stage iron deficiency or iron deficiency anaemia. This study investigated if elevated RDW is causally associated with PAH.A two-sample Mendelian randomisation (MR) approach was applied to investigate whether genetic predisposition to higher levels of RDW increases the odds of developing PAH. Primary and secondary MR analyses were performed using all available genome-wide significant RDW variants (n=179) and five genome-wide significant RDW variants that act via systemic iron status, respectively.We confirmed the observed association between RDW and PAH (OR=1.90, 95% CI=1.80-2.01) in a multi-centre case-control study (N cases=642, N disease controls=15 889). The primary MR analysis was adequately powered to detect a causal effect (OR) from between 1.25 and 1.52 or greater based on estimates reported in the RDW GWAS or from our own data. There was no evidence for a causal association between RDW and PAH in either the primary (ORcausal=1.07, 95% CI=0.92-1.24) or the secondary (ORcausal=1.09, 95% CI=0.77-1.54) MR analysis.The results suggest that at least some of the observed association of RDW with PAH is secondary to disease progression. Results of iron therapeutic trials in PAH should be interpreted with caution as any improvements observed may not be mechanistically linked to the development of PAH.

While widespread genome sequencing ushers in a new era of preventive medicine, the tools for predictive genomics are still lacking. Time and resource limitations mean that human diseases remain uncharacterized because of an inability to predict clinically relevant genetic variants. A strategy of targeting highly conserved protein regions is used commonly in functional studies. However, this benefit is lost for rare diseases where the attributable genes are mostly conserved. An immunological disorder exemplifying this challenge occurs through damaging mutations in RAG1 and RAG2 which presents at an early age with a distinct phenotype of life-threatening immunodeficiency or autoimmunity. Many tools exist for variant pathogenicity prediction, but these cannot account for the probability of variant occurrence. Here, we present a method that predicts the likelihood of mutation for every amino acid residue in the RAG1 and RAG2 proteins. Population genetics data from approximately 146,000 individuals was used for rare variant analysis. Forty-four known pathogenic variants reported in patients and recombination activity measurements from 110 RAG1/2 mutants were used to validate calculated scores. Probabilities were compared with 98 currently known human cases of disease. A genome sequence dataset of 558 patients who have primary immunodeficiency but that are negative for RAG deficiency were also used as validation controls. We compared the difference between mutation likelihood and pathogenicity prediction. Our method builds a map of most probable mutations allowing pre-emptive functional analysis. This method may be applied to other diseases with hopes of improving preparedness for clinical diagnosis.

OBJECTIVES: To determine the frequency of rare and pertinent disease-causing variants in small vessel disease (SVD)-associated genes (such as NOTCH3, HTRA1, COL4A1, COL4A2, FOXC1, TREX1, and GLA) in cerebral SVD, we performed targeted gene sequencing in 950 patients with younger-onset apparently sporadic SVD stroke using a targeted sequencing panel. METHODS: We designed a high-throughput sequencing panel to identify variants in 15 genes (7 known SVD genes, 8 SVD-related disorder genes). The panel was used to screen a population of 950 patients with younger-onset (≤70 years) MRI-confirmed SVD stroke, recruited from stroke centers across the United Kingdom. Variants were filtered according to their frequency in control databases, predicted effect, presence in curated variant lists, and combined annotation dependent depletion scores. Whole genome sequencing and genotyping were performed on a subset of patients to provide a direct comparison of techniques. The frequency of known disease-causing and pertinent variants of uncertain significance was calculated. RESULTS: We identified previously reported variants in 14 patients (8 cysteine-changing NOTCH3 variants in 11 patients, 2 HTRA1 variants in 2 patients, and 1 missense COL4A1 variant in 1 patient). In addition, we identified 29 variants of uncertain significance in 32 patients. CONCLUSION: Rare monogenic variants account for about 1.5% of younger onset lacunar stroke. Most are cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy variants, but the second most common gene affected is HTRA1. A high-throughput sequencing technology platform is an efficient, reliable method to screen for such mutations.

Background Raregenetic variantscause pulmonary arterial hypertension, but the contribution of commongenetic variationto disease risk and natural history is poorly characterised. We tested for genome-wide association for pulmonary arterial hypertension in large international cohorts and assessed the contribution of associated regions to outcomes. Methods We did two separate genome-wide association studies (GWAS) and a meta-analysis of pulmonary arterial hypertension. These GWAS used data from four international case-control studies across 11 744 individuals with European ancestry (including 2085 patients). One GWAS usedgenotypesfrom 5895 whole-genome sequences and the other GWAS used genotyping array data from an additional 5849 individuals.Cross-validationof loci reaching genome-wide significance was sought by meta-analysis. Conditional analysis corrected for the most significant variants at each locus was used to resolve signals for multiple associations. We functionally annotated associated variants and tested associations with duration ofsurvival.All-cause mortalitywas the primary endpoint in survival analyses. Findings A locus nearSOX17(rs10103692, odds ratio 1·80 [95% CI 1·55–2·08], p=5·13 × 10–15) and a second locus inHLA-DPA1andHLA-DPB1(collectively referred to asHLA-DPA1/DPB1here; rs2856830, 1·56 [1·42–1·71], p=7·65 × 10–20) within the class II MHC region were associated with pulmonary arterial hypertension. TheSOX17locus had two independent signals associated with pulmonary arterial hypertension (rs13266183, 1·36 [1·25–1·48], p=1·69 × 10–12; and rs10103692). Functional andepigenomicdata indicate that the risk variants nearSOX17altergene regulationvia anenhanceractive inendothelial cells. Pulmonary arterial hypertension risk variants determined haplotype-specific enhancer activity, and CRISPR-mediated inhibition of the enhancer reducedSOX17expression. TheHLA-DPA1/DPB1rs2856830 genotype was strongly associated with survival. Median survival fromdiagnosisin patients with pulmonary arterial hypertension with the C/C homozygous genotype was double (13·50 years [95% CI 12·07 to >13·50]) that of those with the T/T genotype (6·97 years [6·02–8·05]), despite similar baselinedisease severity. Interpretation This is the first study to report that common genetic variation at loci in an enhancer nearSOX17and inHLA-DPA1/DPB1is associated with pulmonary arterial hypertension. Impairment of SOX17 function might be more common in pulmonary arterial hypertension than suggested by raremutationsinSOX17. Further studies are needed to confirm the association betweenHLA typingor rs2856830 genotyping and survival, and to determine whether HLA typing or rs2856830 genotyping improvesrisk stratificationin clinical practice or trials.

Purpose Increasing numbers of genes are being implicated in Mendelian disorders and incorporated into clinical test panels. However, lack of evidence supporting the gene-disease relationship can hinder interpretation. We explored the utility of testing 51 additional genes for hypertrophic cardiomyopathy (HCM), one of the most commonly tested Mendelian disorders. Methods Using genome sequencing data from 240 sarcomere gene negative HCM cases and 6229 controls, we undertook case-control and individual variant analyses to assess 51 genes that have been proposed for HCM testing. Results We found no evidence to suggest that rare variants in these genes are prevalent causes of HCM. One variant, in a single case, was categorized as likely to be pathogenic. Over 99% of variants were classified as a variant of uncertain significance (VUS) and 54% of cases had one or more VUS. Conclusion For almost all genes, the gene-disease relationship could not be validated and lack of evidence precluded variant interpretation. Thus, the incremental diagnostic yield of extending testing was negligible, and would, we propose, be outweighed by problems that arise with a high rate of uninterpretable findings. These findings highlight the need for rigorous, evidence-based selection of genes for clinical test panels.

Purpose: Mutations in ARL2BP, encoding ADP-ribosylation factor-like 2 binding protein, have recently been implicated as a cause of autosomal recessive retinitis pigmentosa (arRP), with three homozygous variants identified to date. In this study, we performed next-generation sequencing to reveal additional arRP cases associated with ARL2BP variants. Methods: Whole-genome sequencing (WGS) or whole-exome sequencing (WES) was performed in 1,051 unrelated individuals recruited for the UK Inherited Retinal Disease Consortium and NIHR-BioResource Rare Diseases research studies. Sanger sequencing was used to validate the next-generation sequencing data, and reverse transcriptase (RT)–PCR analysis was performed on RNA extracted from blood from affected individuals to test for altered splicing of ARL2BP. Detailed phenotyping was performed, including clinical evaluation, electroretinography, fundus photography, fundus autofluorescence imaging, and spectral-domain optical coherence tomography. Results: Homozygous variants in ARL2BP (NM_012106.3) were identified in two unrelated individuals with RP. The variants, c.207+1G>A and c.390+5G>A, at conserved splice donor sites for intron 3 and intron 5, respectively, were predicted to alter the pre-mRNA splicing of ARL2BP. RT–PCR spanning the affected introns revealed that both variants caused abnormal splicing of ARL2BP in samples from affected individuals. Conclusions: This study identified two homozygous variants in ARL2BP as a rare cause of arRP. Further studies are required to define the underlying disease mechanism causing retinal degeneration as a result of mutations in ARL2BP and any phenotype-genotype correlation associated with residual levels of the wild-type transcript.

BACKGROUND: The genetic etiology of primary immunodeficiency disease (PID) carries prognostic information. OBJECTIVE: We conducted a whole-genome sequencing study assessing a large proportion of the NIHR-BioResource - Rare Disease cohort. METHODS: In the predominantly European study population of principally sporadic unrelated PID cases (n=846), a novel Bayesian method identified NFKB1 as one most strongly associated with PID, and the association was explained by 16 novel heterozygous truncating, missense and gene deletion variants. This accounted for 4% of common variable immunodeficiency (CVID) cases (n=390) in the cohort. Amino-acid substitutions predicted to be pathogenic were assessed by analysis of structural protein data. Immunophenotyping, immunoblotting and ex vivo stimulation of lymphocytes determined the functional effects of these variants. Detailed clinical and pedigree information was collected for genotype-phenotype co-segregation analyses. RESULTS: Both sporadic and familial cases demonstrated evidence of the non-infective complications of CVID, including massive lymphadenopathy (24%), unexplained splenomegaly (48%) and autoimmune disease (48%), features prior studies correlate with worse clinical prognosis. Although partial penetrance of clinical symptoms was noted in certain pedigrees, all carriers have a deficiency in B lymphocyte differentiation. Detailed assessment of B lymphocyte numbers, phenotype and function identifies the presence of a raised CD21lowB cell population: combined with identification of the disease-causing variant, this distinguishes between healthy individuals, asymptomatic carriers and clinically affected cases. CONCLUSION: We show that heterozygous loss-of-function variants in NFKB1 are the most common known monogenic cause of CVID that results in a temporally progressive defect in the formation of immunoglobulin-producing B cells.

BACKGROUND: The genetic cause of primary immunodeficiency disease (PID) carries prognostic information. OBJECTIVE: We conducted a whole-genome sequencing study assessing a large proportion of the NIHR BioResource-Rare Diseases cohort. METHODS: In the predominantly European study population of principally sporadic unrelated PID cases (n = 846), a novel Bayesian method identified nuclear factor κB subunit 1 (NFKB1) as one of the genes most strongly associated with PID, and the association was explained by 16 novel heterozygous truncating, missense, and gene deletion variants. This accounted for 4% of common variable immunodeficiency (CVID) cases (n = 390) in the cohort. Amino acid substitutions predicted to be pathogenic were assessed by means of analysis of structural protein data. Immunophenotyping, immunoblotting, and ex vivo stimulation of lymphocytes determined the functional effects of these variants. Detailed clinical and pedigree information was collected for genotype-phenotype cosegregation analyses. RESULTS: Both sporadic and familial cases demonstrated evidence of the noninfective complications of CVID, including massive lymphadenopathy (24%), unexplained splenomegaly (48%), and autoimmune disease (48%), features prior studies correlated with worse clinical prognosis. Although partial penetrance of clinical symptoms was noted in certain pedigrees, all carriers have a deficiency in B-lymphocyte differentiation. Detailed assessment of B-lymphocyte numbers, phenotype, and function identifies the presence of an increased CD21low B-cell population. Combined with identification of the disease-causing variant, this distinguishes between healthy subjects, asymptomatic carriers, and clinically affected cases. CONCLUSION: We show that heterozygous loss-of-function variants in NFKB1 are the most common known monogenic cause of CVID, which results in a temporally progressive defect in the formation of immunoglobulin-producing B cells.

Background Rare genetic variants cause pulmonary arterial hypertension, but the contribution of common genetic variation to disease risk and natural history is poorly characterised. We tested for genome-wide association for pulmonary arterial hypertension in large international cohorts and assessed the contribution of associated regions to outcomes. Methods We did two separate genome-wide association studies (GWAS) and a meta-analysis of pulmonary arterial hypertension. These GWAS used data from four international case-control studies across 11 744 individuals with European ancestry (including 2085 patients). One GWAS used genotypes from 5895 whole-genome sequences and the other GWAS used genotyping array data from an additional 5849 individuals. Cross-validation of loci reaching genome-wide significance was sought by meta-analysis. Conditional analysis corrected for the most significant variants at each locus was used to resolve signals for multiple associations. We functionally annotated associated variants and tested associations with duration of survival. All-cause mortality was the primary endpoint in survival analyses. Findings A locus near SOX17 (rs10103692, odds ratio 1·80 [95% CI 1·55–2·08], p=5·13 × 10–15) and a second locus in HLA-DPA1 and HLA-DPB1 (collectively referred to as HLA-DPA1/DPB1 here; rs2856830, 1·56 [1·42–1·71], p=7·65 × 10–20) within the class II MHC region were associated with pulmonary arterial hypertension. The SOX17 locus had two independent signals associated with pulmonary arterial hypertension (rs13266183, 1·36 [1·25–1·48], p=1·69 × 10–12; and rs10103692). Functional and epigenomic data indicate that the risk variants near SOX17 alter gene regulation via an enhancer active in endothelial cells. Pulmonary arterial hypertension risk variants determined haplotype-specific enhancer activity, and CRISPR-mediated inhibition of the enhancer reduced SOX17 expression. The HLA-DPA1/DPB1 rs2856830 genotype was strongly associated with survival. Median survival from diagnosis in patients with pulmonary arterial hypertension with the C/C homozygous genotype was double (13·50 years [95% CI 12·07 to >13·50]) that of those with the T/T genotype (6·97 years [6·02–8·05]), despite similar baseline disease severity. Interpretation This is the first study to report that common genetic variation at loci in an enhancer near SOX17 and in HLA-DPA1/DPB1 is associated with pulmonary arterial hypertension. Impairment of SOX17 function might be more common in pulmonary arterial hypertension than suggested by rare mutations in SOX17. Further studies are needed to confirm the association between HLA typing or rs2856830 genotyping and survival, and to determine whether HLA typing or rs2856830 genotyping improves risk stratification in clinical practice or trials.

Multiple primary tumors (MPTs) affect a substantial proportion of cancer survivors and can result from various causes, including inherited predisposition. Currently, germline genetic testing of MPT-affected individuals for variants in cancer-predisposition genes (CPGs) is mostly targeted by tumor type. We ascertained pre-assessed MPT individuals (with at least two primary tumors by age 60 years or at least three by 70 years) from genetics centers and performed whole-genome sequencing (WGS) on 460 individuals from 440 families. Despite previous negative genetic assessment and molecular investigations, pathogenic variants in moderate- and high-risk CPGs were detected in 67/440 (15.2%) probands. WGS detected variants that would not be (or were not) detected by targeted resequencing strategies, including low-frequency structural variants (6/440 [1.4%] probands). In most individuals with a germline variant assessed as pathogenic or likely pathogenic (P/LP), at least one of their tumor types was characteristic of variants in the relevant CPG. However, in 29 probands (42.2% of those with a P/LP variant), the tumor phenotype appeared discordant. The frequency of individuals with truncating or splice-site CPG variants and at least one discordant tumor type was significantly higher than in a control population (χ2 = 43.642; p ≤ 0.0001). 2/67 (3%) probands with P/LP variants had evidence of multiple inherited neoplasia allele syndrome (MINAS) with deleterious variants in two CPGs. Together with variant detection rates from a previous series of similarly ascertained MPT-affected individuals, the present results suggest that first-line comprehensive CPG analysis in an MPT cohort referred to clinical genetics services would detect a deleterious variant in about a third of individuals. ispartof: AMERICAN JOURNAL OF HUMAN GENETICS vol:103 issue:1 pages:3-18 ispartof: location:United States status: published

Purpose: Recessive mutations in CLN7/MFSD8 usually cause variant late-infantile onset neuronal ceroid lipofuscinosis (vLINCL), a poorly understood neurodegenerative condition, though mutations may also cause nonsyndromic maculopathy. A series of 12 patients with nonsyndromic retinopathy due to novel CLN7/MFSD8 mutation combinations were investigated in this study. Methods: Affected patients and their family members were recruited in ophthalmic clinics at each center where they were examined by retinal imaging and detailed electrophysiology. Whole exome or genome next generation sequencing was performed on genomic DNA from at least one affected family member. Immunofluorescence confocal microscopy of murine retina cross-sections were used to localize the protein. Results: Compound heterozygous alleles were identified in six cases, one of which was always p.Glu336Gln. Such combinations resulted in isolated macular disease. Six further cases were homozygous for the variant p.Met454Thr, identified as a founder mutation of South Asian origin. Those patients had widespread generalized retinal disease, characterized by electroretinography as a rod-cone dystrophy with severe macular involvement. In addition, the photopic single flash electroretinograms demonstrated a reduced b- to a-wave amplitude ratio, suggesting dysfunction occurring after phototransduction. Immunohistology identified MFSD8 in the outer plexiform layer of the retina, a site rich in photoreceptor synapses. Conclusions: This study highlights a hierarchy of MFSD8 variant severity, predicting three consequences of mutation: (1) nonsyndromic localized maculopathy, (2) nonsyndromic widespread retinopathy, or (3) syndromic neurological disease. The data also shed light on the underlying pathogenesis by implicating the photoreceptor synaptic terminals as the major site of retinal disease.

Histone lysine methylation, mediated by mixed-lineage leukemia (MLL) proteins, is now known to be critical in the regulation of gene expression, genomic stability, cell cycle and nuclear architecture. Despite MLL proteins being postulated as essential for normal development, little is known about the specific functions of the different MLL lysine methyltransferases. Here we report heterozygous variants in the gene KMT2B (also known as MLL4) in 27 unrelated individuals with a complex progressive childhood-onset dystonia, often associated with a typical facial appearance and characteristic brain magnetic resonance imaging findings. Over time, the majority of affected individuals developed prominent cervical, cranial and laryngeal dystonia. Marked clinical benefit, including the restoration of independent ambulation in some cases, was observed following deep brain stimulation (DBS). These findings highlight a clinically recognizable and potentially treatable form of genetic dystonia, demonstrating the crucial role of KMT2B in the physiological control of voluntary movement.

Heritable platelet function disorders (PFDs) are genetically heterogeneous and poorly characterized. Pathogenic variants in RASGRP2, which encodes calcium and diacylglycerol-regulated guanine exchange factor I (CalDAG-GEFI), have been reported previously in 3 pedigrees with bleeding and reduced platelet aggregation responses. To better define the phenotype associated with pathogenic RASGRP2 variants, we compared high-throughput sequencing and phenotype data from 2042 cases in pedigrees with unexplained bleeding or platelet disorders to data from 5422 controls. Eleven cases harbored 11 different, previously unreported RASGRP2 variants that were biallelic and likely pathogenic. The variants included 5 high-impact variants predicted to prevent CalDAG-GEFI expression and 6 missense variants affecting the CalDAG-GEFI CDC25 domain, which mediates Rap1 activation during platelet inside-out αIIbβ3 signaling. Cases with biallelic RASGRP2 variants had abnormal mucocutaneous, surgical, and dental bleeding from childhood, requiring ≥1 blood or platelet transfusion in 78% of cases. Platelets displayed reduced aggregation in response to adenosine 5'-diphosphate and epinephrine, but variable aggregation defects with other agonists. There were no other consistent clinical or laboratory features. These data enable definition of human CalDAG-GEFI deficiency as a nonsyndromic, recessive PFD associated with a moderate or severe bleeding phenotype and complex defects in platelet aggregation.

Expression quantitative trait loci (eQTL) can provide a link between disease susceptibility variants discovered by genetic association studies and biology. To date, eQTL mapping studies have been primarily conducted in healthy individuals from population-based cohorts. Genetic effects have been known to be context-specific and vary with changing environmental stimuli. We conducted a transcriptome- and genome-wide eQTL mapping study in a cohort of patients with idiopathic or heritable pulmonary arterial hypertension (PAH) using RNA sequencing (RNAseq) data from whole blood. We sought confirmation from three published population-based eQTL studies, including the GTEx Project, and followed up potentially novel eQTL not observed in the general population. In total, we identified 2314 eQTL of which 90% were cis-acting and 75% were confirmed by at least one of the published studies. While we observed a higher GWAS trait colocalization rate among confirmed eQTL, colocalisation rate of novel eQTL reported for lung-related phenotypes was twice as high as that of confirmed eQTL. Functional enrichment analysis of genes with novel eQTL in PAH highlighted immune-related processes, a suspected contributor to PAH. These potentially novel eQTL specific to or active in PAH could be useful in understanding genetic risk factors for other diseases that share common mechanisms with PAH.

Histone lysine methylation, mediated by mixed-lineage leukemia (MLL) proteins, is now known to be critical in the regulation of gene expression, genomic stability, cell cycle and nuclear architecture. Despite MLL proteins being postulated as essential for normal development, little is known about the specific functions of the different MLL lysine methyltransferases. Here we report heterozygous variants in the gene KMT2B (also known as MLL4) in 27 unrelated individuals with a complex progressive childhood-onset dystonia, often associated with a typical facial appearance and characteristic brain magnetic resonance imaging findings. Over time, the majority of affected individuals developed prominent cervical, cranial and laryngeal dystonia. Marked clinical benefit, including the restoration of independent ambulation in some cases, was observed following deep brain stimulation (DBS). These findings highlight a clinically recognizable and potentially treatable form of genetic dystonia, demonstrating the crucial role of KMT2B in the physiological control of voluntary movement.

While widespread genome sequencing ushers in a new era of preventive medicine, the tools for predictive genomics are still lacking. Time and resource limitations mean that human diseases remain uncharacterized because of an inability to predict clinically relevant genetic variants. A strategy of targeting highly conserved protein regions is used commonly in functional studies. However, this benefit is lost for rare diseases where the attributable genes are mostly conserved. An immunological disorder exemplifying this challenge occurs through damaging mutations in RAG1 and RAG2 which presents at an early age with a distinct phenotype of life-threatening immunodeficiency or autoimmunity. Many tools exist for variant pathogenicity prediction, but these cannot account for the probability of variant occurrence. Here, we present a method that predicts the likelihood of mutation for every amino acid residue in the RAG1 and RAG2 proteins. Population genetics data from approximately 146,000 individuals was used for rare variant analysis. Forty-four known pathogenic variants reported in patients and recombination activity measurements from 110 RAG1/2 mutants were used to validate calculated scores. Probabilities were compared with 98 currently known human cases of disease. A genome sequence dataset of 558 patients who have primary immunodeficiency but that are negative for RAG deficiency were also used as validation controls. We compared the difference between mutation likelihood and pathogenicity prediction. Our method builds a map of most probable mutations allowing pre-emptive functional analysis. This method may be applied to other diseases with hopes of improving preparedness for clinical diagnosis.

Expression quantitative trait loci (eQTL) can provide a link between disease susceptibility variants discovered by genetic association studies and biology. To date, eQTL mapping studies have been primarily conducted in healthy individuals from population-based cohorts. Genetic effects have been known to be context-specific and vary with changing environmental stimuli. We conducted a transcriptome- and genome-wide eQTL mapping study in a cohort of patients with idiopathic or heritable pulmonary arterial hypertension (PAH) using RNA sequencing (RNAseq) data from whole blood. We sought confirmation from three published population-based eQTL studies, including the GTEx Project, and followed up potentially novel eQTL not observed in the general population. In total, we identified 2314 eQTL of which 90% were cis-acting and 75% were confirmed by at least one of the published studies. While we observed a higher GWAS trait colocalization rate among confirmed eQTL, colocalisation rate of novel eQTL reported for lung-related phenotypes was twice as high as that of confirmed eQTL. Functional enrichment analysis of genes with novel eQTL in PAH highlighted immune-related processes, a suspected contributor to PAH. These potentially novel eQTL specific to or active in PAH could be useful in understanding genetic risk factors for other diseases that share common mechanisms with PAH.

BACKGROUND: Rare genetic variants cause pulmonary arterial hypertension, but the contribution of common genetic variation to disease risk and natural history is poorly characterised. We tested for genome-wide association for pulmonary arterial hypertension in large international cohorts and assessed the contribution of associated regions to outcomes. METHODS: We did two separate genome-wide association studies (GWAS) and a meta-analysis of pulmonary arterial hypertension. These GWAS used data from four international case-control studies across 11 744 individuals with European ancestry (including 2085 patients). One GWAS used genotypes from 5895 whole-genome sequences and the other GWAS used genotyping array data from an additional 5849 individuals. Cross-validation of loci reaching genome-wide significance was sought by meta-analysis. Conditional analysis corrected for the most significant variants at each locus was used to resolve signals for multiple associations. We functionally annotated associated variants and tested associations with duration of survival. All-cause mortality was the primary endpoint in survival analyses. FINDINGS: A locus near SOX17 (rs10103692, odds ratio 1·80 [95% CI 1·55-2·08], p=5·13 × 10-15) and a second locus in HLA-DPA1 and HLA-DPB1 (collectively referred to as HLA-DPA1/DPB1 here; rs2856830, 1·56 [1·42-1·71], p=7·65 × 10-20) within the class II MHC region were associated with pulmonary arterial hypertension. The SOX17 locus had two independent signals associated with pulmonary arterial hypertension (rs13266183, 1·36 [1·25-1·48], p=1·69 × 10-12; and rs10103692). Functional and epigenomic data indicate that the risk variants near SOX17 alter gene regulation via an enhancer active in endothelial cells. Pulmonary arterial hypertension risk variants determined haplotype-specific enhancer activity, and CRISPR-mediated inhibition of the enhancer reduced SOX17 expression. The HLA-DPA1/DPB1 rs2856830 genotype was strongly associated with survival. Median survival from diagnosis in patients with pulmonary arterial hypertension with the C/C homozygous genotype was double (13·50 years [95% CI 12·07 to >13·50]) that of those with the T/T genotype (6·97 years [6·02-8·05]), despite similar baseline disease severity. INTERPRETATION: This is the first study to report that common genetic variation at loci in an enhancer near SOX17 and in HLA-DPA1/DPB1 is associated with pulmonary arterial hypertension. Impairment of SOX17 function might be more common in pulmonary arterial hypertension than suggested by rare mutations in SOX17. Further studies are needed to confirm the association between HLA typing or rs2856830 genotyping and survival, and to determine whether HLA typing or rs2856830 genotyping improves risk stratification in clinical practice or trials. FUNDING: UK NIHR, BHF, UK MRC, Dinosaur Trust, NIH/NHLBI, ERS, EMBO, Wellcome Trust, EU, AHA, ACClinPharm, Netherlands CVRI, Dutch Heart Foundation, Dutch Federation of UMC, Netherlands OHRD and RNAS, German DFG, German BMBF, APH Paris, INSERM, Université Paris-Sud, and French ANR.

International audience; Heritable platelet function disorders (PFDs) are genetically heterogeneous and poorly characterised. Pathogenic variants in RASGRP2, which encodes calcium and diacylglycerol-regulated guanine exchange factor I (CalDAG-GEFI), have been reported previously in three pedigrees with bleeding and reduced platelet aggregation responses. To better define the phenotype associated with pathogenic RASGRP2 variants, we compared high-throughput sequencing and phenotype data from 2,042 cases in pedigrees with unexplained bleeding or platelet disorders to data from 5,422 controls. Eleven cases harboured 11 different, previously unreported RASGRP2 variants that were biallelic and likely pathogenic. The variants included five high-impact variants predicted to prevent CalDAG-GEFI expression and six missense variants affecting the CalDAG-GEFI CDC25 domain, which mediates Rap1 activation during platelet inside-out αIIbβ3 signalling. Cases with biallelic RASGRP2 variants had abnormal mucocutaneous, surgical and dental bleeding from childhood, requiring at least one blood or platelet transfusion in 78% of cases. Platelets displayed reduced aggregation in response to ADP and epinephrine, but variable aggregation defects with other agonists. There were no other consistent clinical or laboratory features. These data enable definition of human CalDAG-GEFI deficiency as a non-syndromic, recessive PFD associated with a moderate or severe bleeding phenotype and complex defects in platelet aggregation.

Mg2+ plays a vital role in platelet function, but despite implications for life-threatening conditions such as stroke or myocardial infarction, the mechanisms controlling [Mg2+]i in megakaryocytes (MKs) and platelets are largely unknown. Transient receptor potential melastatin-like 7 channel (TRPM7) is a ubiquitous, constitutively active cation channel with a cytosolic α-kinase domain that is critical for embryonic development and cell survival. Here we report that impaired channel function of TRPM7 in MKs causes macrothrombocytopenia in mice (Trpm7fl/fl-Pf4Cre) and likely in several members of a human pedigree that, in addition, suffer from atrial fibrillation. The defect in platelet biogenesis is mainly caused by cytoskeletal alterations resulting in impaired proplatelet formation by Trpm7fl/fl-Pf4Cre MKs, which is rescued by Mg2+ supplementation or chemical inhibition of non-muscle myosin IIA heavy chain activity. Collectively, our findings reveal that TRPM7 dysfunction may cause macrothrombocytopenia in humans and mice., Although Mg2+ is vital for platelet activation and aggregation, its regulation in these cells is still largely unknown. Here, the authors show that TRPM7, a cation channel and a protein kinase, regulates thrombopoiesis and platelet size by affecting the cytoskeleton of these cells in mice and humans.