The diagnostic odyssey in children and adolescents with X-linked hypophosphataemia: population-based, case-control study.

CONTEXT: X-linked hypophosphataemia (XLH) is a rare genetic disorder causing renal phosphate wasting, which predicates musculoskeletal manifestations such as rickets. Diagnosis is often delayed. OBJECTIVE: To explore the recording of clinical features, and the diagnostic odyssey of children and adolescents with XLH in primary care electronic healthcare records (EHR) in the United Kingdom. METHODS: Using the Optimum Patient Care Research Database, individuals aged 20 years or younger after 1st Jan 2000 at date of recorded XLH diagnosis were identified using SNOMED/Read codes, and age-matched to 100 controls. Recording of XLH-related clinical features was summarised, then compared between cases and controls using Chi-squared or Fisher's exact tests. RESULTS: 261 XLH cases were identified; 99 met inclusion criteria. Of these, 84/99 had at least 1 XLH-related clinical feature recorded in their primary care EHR. Clinical codes for rickets, genu varum and low phosphate were recorded prior to XLH diagnosis in under 20% of cases (median of 1, 1, and 3 years prior, respectively). Rickets, genu varum, low phosphate, nephrocalcinosis and growth delay were significantly more likely to be recorded in cases. CONCLUSION: This characterisation of the EHR phenotypes of children and adolescents with XLH may inform future case-finding approaches to expedite diagnosis in primary care.

Whole genome sequences discriminate hereditary hemorrhagic telangiectasia phenotypes by non-HHT deleterious DNA variation.

The abnormal vascular structures of hereditary hemorrhagic telangiectasia (HHT) often cause severe anemia due to recurrent hemorrhage, but HHT causal genes do not predict the severity of hematological complications. We tested for chance inheritance and clinical associations of rare deleterious variants in which loss-of-function causes bleeding or hemolytic disorders in the general population. In double-blinded analyses, all 104 patients with HHT from a single reference center recruited to the 100 000 Genomes Project were categorized on new MALO (more/as-expected/less/opposite) sub-phenotype severity scales, and whole genome sequencing data were tested for high impact variants in 75 HHT-independent genes encoding coagulation factors, or platelet, hemoglobin, erythrocyte enzyme, and erythrocyte membrane constituents. Rare variants (all gnomAD allele frequencies <0.003) were identified in 56 (75%) of these 75 HHT-unrelated genes. Deleteriousness assignments by Combined Annotation Dependent Depletion (CADD) scores >15 were supported by gene-level mutation significance cutoff scores. CADD >15 variants were identified in 38/104 (36.5%) patients with HHT, found for 1 in 10 patients within platelet genes; 1 in 8 within coagulation genes; and 1 in 4 within erythrocyte hemolytic genes. In blinded analyses, patients with greater hemorrhagic severity that had been attributed solely to HHT vessels had more CADD-deleterious variants in platelet (Spearman ρ = 0.25; P = .008) and coagulation (Spearman ρ = 0.21; P = .024) genes. However, the HHT cohort had 60% fewer deleterious variants in platelet and coagulation genes than expected (Mann-Whitney test P = .021). In conclusion, patients with HHT commonly have rare variants in genes of relevance to their phenotype, offering new therapeutic targets and opportunities for informed, personalized medicine strategies.

Identification and validation of a novel pathogenic variant in GDF2 (BMP9) responsible for hereditary hemorrhagic telangiectasia and pulmonary arteriovenous malformations.

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant multisystemic vascular dysplasia, characterized by arteriovenous malformations (AVMs), mucocutaneous telangiectasia and nosebleeds. HHT is caused by a heterozygous null allele in ACVRL1, ENG, or SMAD4, which encode proteins mediating bone morphogenetic protein (BMP) signaling. Several missense and stop-gain variants identified in GDF2 (encoding BMP9) have been reported to cause a vascular anomaly syndrome similar to HHT, however none of these patients met diagnostic criteria for HHT. HHT families from UK NHS Genomic Medicine Centres were recruited to the Genomics England 100,000 Genomes Project. Whole genome sequencing and tiering protocols identified a novel, heterozygous GDF2 sequence variant in all three affected members of one HHT family who had previously screened negative for ACVRL1, ENG, and SMAD4. All three had nosebleeds and typical HHT telangiectasia, and the proband also had severe pulmonary AVMs from childhood. In vitro studies showed the mutant construct expressed the proprotein but lacked active mature BMP9 dimer, suggesting the mutation disrupts correct cleavage of the protein. Plasma BMP9 levels in the patients were significantly lower than controls. In conclusion, we propose that this heterozygous GDF2 variant is a rare cause of HHT associated with pulmonary AVMs.

Evaluating the performance of a clinical genome sequencing program for diagnosis of rare genetic disease, seen through the lens of craniosynostosis.

PURPOSE: Genome sequencing (GS) for diagnosis of rare genetic disease is being introduced into the clinic, but the complexity of the data poses challenges for developing pipelines with high diagnostic sensitivity. We evaluated the performance of the Genomics England 100,000 Genomes Project (100kGP) panel-based pipelines, using craniosynostosis as a test disease. METHODS: GS data from 114 probands with craniosynostosis and their relatives (314 samples), negative on routine genetic testing, were scrutinized by a specialized research team, and diagnoses compared with those made by 100kGP. RESULTS: Sixteen likely pathogenic/pathogenic variants were identified by 100kGP. Eighteen additional likely pathogenic/pathogenic variants were identified by the research team, indicating that for craniosynostosis, 100kGP panels had a diagnostic sensitivity of only 47%. Measures that could have augmented diagnoses were improved calling of existing panel genes (+18% sensitivity), review of updated panels (+12%), comprehensive analysis of de novo small variants (+29%), and copy-number/structural variants (+9%). Recent NHS England recommendations that partially incorporate these measures should achieve 85% overall sensitivity (+38%). CONCLUSION: GS identified likely pathogenic/pathogenic variants in 29.8% of previously undiagnosed patients with craniosynostosis. This demonstrates the value of research analysis and the importance of continually improving algorithms to maximize the potential of clinical GS.

Correction: SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

PURPOSE: Enrichment of heterozygous missense and truncating SMAD6 variants was previously reported in nonsyndromic sagittal and metopic synostosis, and interaction of SMAD6 variants with a common polymorphism nearBMP2 (rs1884302) was proposed to contribute to inconsistent penetrance. We determined the occurrence of SMAD6 variants in all types of craniosynostosis, evaluated the impact of different missense variants on SMAD6 function, and tested independently whether rs1884302 genotype significantly modifies the phenotype. METHODS: We performed resequencing of SMAD6 in 795 unsolved patients with any type of craniosynostosis and genotyped rs1884302 in SMAD6-positive individuals and relatives. We examined the inhibitory activity and stability of SMAD6 missense variants. RESULTS: We found 18 (2.3%) different rare damaging SMAD6 variants, with the highest prevalence in metopic synostosis (5.8%) and an 18.3-fold enrichment of loss-of-function variants comparedwith gnomAD data (P < 10-7). Combined with eight additional variants, ≥20/26 were transmitted from an unaffected parent but rs1884302 genotype did not predict phenotype. CONCLUSION: Pathogenic SMAD6 variants substantially increase the risk of both nonsyndromic and syndromic presentations of craniosynostosis, especially metopic synostosis. Functional analysis is important to evaluate missense variants. Genotyping of rs1884302 is not clinically useful. Mechanisms to explain the remarkable diversity of phenotypes associated with SMAD6 variants remain obscure.

Functional Analysis of a Carotid Intima-Media Thickness Locus Implicates BCAR1 and Suggests a Causal Variant.

BACKGROUND: Carotid intima-media thickness (IMT) is a marker of subclinical atherosclerosis that can predict cardiovascular disease events over traditional risk factors. This study examined the BCAR1-CFDP1-TMEM170A locus on chromosome 16, associated with carotid IMT and coronary artery disease in the IMT and IMT-Progression as Predictors of Vascular Events (IMPROVE) cohort, to identify the functional variant. METHODS AND RESULTS: In analysis of the locus lead single nucleotide polymorphism (SNP; rs4888378, intronic in CFDP1) in Progressione della Lesione Intimale Carotidea (PLIC), the protective AA genotype was associated with slower IMT progression in women (P=0.04) but not in men. Meta-analysis of 5 cohort studies also supported a protective effect of the A allele on common carotid IMT in women only (women: ß=-0.0047, P=1.63 × 10(-4); men: ß=-0.0029, P=0.0678). Two hundred fourteen noncoding variants in strong linkage disequilibrium (r(2) ≥ 0.8) with rs4888378 were identified from 1000 Genome Project. ENCODE regulatory chromatin marks were used to create a shortlist of 6 possible regulatory variants. Electrophoretic mobility shift assays on the shortlist detected allele-specific protein binding to the lead SNP rs4888378; multiplexed competitor electrophoretic mobility shift assays implicated FOXA as the protein. Luciferase reporter assays on rs4888378 showed a significant 35% to 92% (P=0.0057; P=4.0 × 10(-22)) decrease in gene expression with the A allele. Expression quantitative trait loci analysis confirmed previously reported associations of rs4888378 with BCAR1 in vascular tissues. CONCLUSIONS: Molecular studies suggest the lead SNP as a potentially causal SNP at the BCAR1-CFDP1-TMEM170A locus, and expression quantitative trait loci studies implicate BCAR1 as the causal gene. This variant showed stronger effects on common carotid IMT in women, raising questions about the mechanism of the causal SNP on atherosclerosis.