Combined exome and whole-genome sequencing identifies mutations in ARMC4 as a cause of primary ciliary dyskinesia with defects in the outer dynein arm

Primary ciliary dyskinesia (PCD) is a rare, genetically heterogeneous ciliopathy disorder affecting cilia and sperm motility. A range of ultrastructural defects of the axoneme underlie the disease, which is characterised by chronic respiratory symptoms and obstructive lung disease, infertility and body axis laterality defects. We applied a next-generation sequencing approach to identify the gene responsible for this phenotype in two consanguineous families

DNAAF1 links heart laterality with the AAA+ ATPase RUVBL1 and ciliary intraflagellar transport

DNAAF1 (LRRC50) is a cytoplasmic protein required for dynein heavy chain assembly and cilia motility, and DNAAF1 mutations cause primary ciliary dyskinesia (PCD; MIM 613193). We describe four families with DNAAF1 mutations and complex congenital heart disease (CHD). In three families, all affected individuals have typical PCD phenotypes. However, an additional family demonstrates isolated CHD (heterotaxy) in two affected siblings, but no clinical evidence of PCD. We identified a homozygous DNAAF1 missense mutation, p.Leu191Phe, as causative for heterotaxy in this family. Genetic complementation in dnaaf1-null zebrafish embryos demonstrated the rescue of normal heart looping with wild-type human DNAAF1, but not the p.Leu191Phe variant, supporting the conserved pathogenicity of this DNAAF1 missense mutation. This observation points to a phenotypic continuum between CHD and PCD, providing new insights into the pathogenesis of isolated CHD. In further investigations of the function of DNAAF1 in dynein arm assembly, we identified interactions with members of a putative dynein arm assembly complex. These include the ciliary intraflagellar transport protein IFT88 and the AAA+ (ATPases Associated with various cellular Activities) family proteins RUVBL1 (Pontin) and RUVBL2 (Reptin). Co-localization studies support these findings, with the loss of RUVBL1 perturbing the co-localization of DNAAF1 with IFT88. We show that RUVBL1 orthologues have an asymmetric left-sided distribution at both the mouse embryonic node and the Kupffer’s vesicle in zebrafish embryos, with the latter asymmetry dependent on DNAAF1. These results suggest that DNAAF1-RUVBL1 biochemical and genetic interactions have a novel functional role in symmetry breaking and cardiac development

Phenotypic Characterization of EIF2AK4 Mutation Carriers in a Large Cohort of Patients Diagnosed Clinically With Pulmonary Arterial Hypertension.

BACKGROUND: Pulmonary arterial hypertension (PAH) is a rare disease with an emerging genetic basis. Heterozygous mutations in the gene encoding the bone morphogenetic protein receptor type 2 (BMPR2) are the commonest genetic cause of PAH, whereas biallelic mutations in the eukaryotic translation initiation factor 2 alpha kinase 4 gene (EIF2AK4) are described in pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis. Here, we determine the frequency of these mutations and define the genotype-phenotype characteristics in a large cohort of patients diagnosed clinically with PAH. METHODS: Whole-genome sequencing was performed on DNA from patients with idiopathic and heritable PAH and with pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis recruited to the National Institute of Health Research BioResource-Rare Diseases study. Heterozygous variants in BMPR2 and biallelic EIF2AK4 variants with a minor allele frequency of <1:10 000 in control data sets and predicted to be deleterious (by combined annotation-dependent depletion, PolyPhen-2, and sorting intolerant from tolerant predictions) were identified as potentially causal. Phenotype data from the time of diagnosis were also captured. RESULTS: Eight hundred sixty-four patients with idiopathic or heritable PAH and 16 with pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis were recruited. Mutations in BMPR2 were identified in 130 patients (14.8%). Biallelic mutations in EIF2AK4 were identified in 5 patients with a clinical diagnosis of pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis. Furthermore, 9 patients with a clinical diagnosis of PAH carried biallelic EIF2AK4 mutations. These patients had a reduced transfer coefficient for carbon monoxide (Kco; 33% [interquartile range, 30%-35%] predicted) and younger age at diagnosis (29 years; interquartile range, 23-38 years) and more interlobular septal thickening and mediastinal lymphadenopathy on computed tomography of the chest compared with patients with PAH without EIF2AK4 mutations. However, radiological assessment alone could not accurately identify biallelic EIF2AK4 mutation carriers. Patients with PAH with biallelic EIF2AK4 mutations had a shorter survival. CONCLUSIONS: Biallelic EIF2AK4 mutations are found in patients classified clinically as having idiopathic and heritable PAH. These patients cannot be identified reliably by computed tomography, but a low Kco and a young age at diagnosis suggests the underlying molecular diagnosis. Genetic testing can identify these misclassified patients, allowing appropriate management and early referral for lung transplantation

Comprehensive Rare Variant Analysis via Whole-Genome Sequencing to Determine the Molecular Pathology of Inherited Retinal Disease

Inherited retinal disease is a common cause of visual impairment and represents a highly heterogeneous group of conditions. Here, we present findings from a cohort of 722 individuals with inherited retinal disease, who have had whole-genome sequencing (n = 605), whole-exome sequencing (n = 72), or both (n = 45) performed, as part of the NIHR-BioResource Rare Diseases research study. We identified pathogenic variants (single-nucleotide variants, indels, or structural variants) for 404/722 (56%) individuals. Whole-genome sequencing gives unprecedented power to detect three categories of pathogenic variants in particular: structural variants, variants in GC-rich regions, which have significantly improved coverage compared to whole-exome sequencing, and variants in non-coding regulatory regions. In addition to previously reported pathogenic regulatory variants, we have identified a previously unreported pathogenic intronic variant in $\textit{CHM}$ in two males with choroideremia. We have also identified 19 genes not previously known to be associated with inherited retinal disease, which harbor biallelic predicted protein-truncating variants in unsolved cases. Whole-genome sequencing is an increasingly important comprehensive method with which to investigate the genetic causes of inherited retinal disease.This work was supported by The National Institute for Health Research England (NIHR) for the NIHR BioResource – Rare Diseases project (grant number RG65966). The Moorfields Eye Hospital cohort of patients and clinical and imaging data were ascertained and collected with the support of grants from the National Institute for Health Research Biomedical Research Centre at Moorfields Eye Hospital, National Health Service Foundation Trust, and UCL Institute of Ophthalmology, Moorfields Eye Hospital Special Trustees, Moorfields Eye Charity, the Foundation Fighting Blindness (USA), and Retinitis Pigmentosa Fighting Blindness. M.M. is a recipient of an FFB Career Development Award. E.M. is supported by UCLH/UCL NIHR Biomedical Research Centre. F.L.R. and D.G. are supported by Cambridge NIHR Biomedical Research Centre

ZSWIM7 Is associated with human female meiosis and familial primary ovarian insufficiency

Background Primary ovarian insufficiency (POI) affects 1% of women and is associated with significant medical consequences. A genetic cause for POI can be found in up to 30% of women, elucidating key roles for these genes in human ovary development. Objective We aimed to identify the genetic mechanism underlying early-onset POI in 2 sisters from a consanguineous pedigree. Methods Genome sequencing and variant filtering using an autosomal recessive model was performed in the 2 affected sisters and their unaffected family members. Quantitative reverse transcriptase PCR (qRT-PCR) and RNA sequencing were used to study the expression of key genes at critical stages of human fetal gonad development (Carnegie Stage 22/23, 9 weeks post conception (wpc), 11 wpc, 15/16 wpc, 19/20 wpc) and in adult tissue. Results Only 1 homozygous variant cosegregating with the POI phenotype was found: a single nucleotide substitution in zinc finger SWIM-type containing 7 (ZSWIM7), NM_001042697.2: c.173C > G; resulting in predicted loss-of-function p.(Ser58*). qRT-PCR demonstrated higher expression of ZSWIM7 in the 15/16 wpc ovary compared with testis, corresponding to peak meiosis in the fetal ovary. RNA sequencing of fetal gonad samples showed that ZSWIM7 has a similar temporal expression profile in the developing ovary to other homologous recombination genes. Main conclusions Disruption of ZSWIM7 is associated with POI in humans. ZSWIM7 is likely to be important for human homologous recombination; these findings expand the range of genes associated with POI in women

STAG3 truncating variant as the cause of primary ovarian insufficiency

Primary ovarian insufficiency (POI) is a distressing cause of infertility in young women. POI is heterogeneous with only a few causative genes having been discovered so far. Our objective was to determine the genetic cause of POI in a consanguineous Lebanese family with two affected sisters presenting with primary amenorrhoea and an absence of any pubertal development. Multipoint parametric linkage analysis was performed. Whole-exome sequencing was done on the proband. Linkage analysis identified a locus on chromosome 7 where exome sequencing successfully identified a homozygous two base pair duplication (c.1947_48dupCT), leading to a truncated protein p.(Y650Sfs*22) in the STAG3 gene, confirming it as the cause of POI in this family. Exome sequencing combined with linkage analyses offers a powerful tool to efficiently find novel genetic causes of rare, heterogeneous disorders, even in small single families. This is only the second report of a STAG3 variant; the first STAG3 variant was recently described in a phenotypically similar family with extreme POI. Identification of an additional family highlights the importance of STAG3 in POI pathogenesis and suggests it should be evaluated in families affected with POI

Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

Telomere length is a risk factor in disease and the dynamics of telomere length are crucial to our understanding of cell replication and vitality. The proliferation of whole genome sequencing represents an unprecedented opportunity to glean new insights into telomere biology on a previously unimaginable scale. To this end, a number of approaches for estimating telomere length from whole-genome sequencing data have been proposed. Here we present Telomerecat, a novel approach to the estimation of telomere length. Previous methods have been dependent on the number of telomeres present in a cell being known, which may be problematic when analysing aneuploid cancer data and non-human samples. Telomerecat is designed to be agnostic to the number of telomeres present, making it suited for the purpose of estimating telomere length in cancer studies. Telomerecat also accounts for interstitial telomeric reads and presents a novel approach to dealing with sequencing errors. We show that Telomerecat performs well at telomere length estimation when compared to leading experimental and computational methods. Furthermore, we show that it detects expected patterns in longitudinal data, repeated measurements, and cross-species comparisons. We also apply the method to a cancer cell data, uncovering an interesting relationship with the underlying telomerase genotype

Phenotypic Characterization of <i>EIF2AK4</i> Mutation Carriers in a Large Cohort of Patients Diagnosed Clinically With Pulmonary Arterial Hypertension

Background: Pulmonary arterial hypertension (PAH) is a rare disease with an emerging genetic basis. Heterozygous mutations in the gene encoding the bone morphogenetic protein receptor type 2 ( BMPR2 ) are the commonest genetic cause of PAH, whereas biallelic mutations in the eukaryotic translation initiation factor 2 alpha kinase 4 gene ( EIF2AK4 ) are described in pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis. Here, we determine the frequency of these mutations and define the genotype-phenotype characteristics in a large cohort of patients diagnosed clinically with PAH. Methods: Whole-genome sequencing was performed on DNA from patients with idiopathic and heritable PAH and with pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis recruited to the National Institute of Health Research BioResource–Rare Diseases study. Heterozygous variants in BMPR2 and biallelic EIF2AK4 variants with a minor allele frequency of &lt;1:10 000 in control data sets and predicted to be deleterious (by combined annotation-dependent depletion, PolyPhen-2, and sorting intolerant from tolerant predictions) were identified as potentially causal. Phenotype data from the time of diagnosis were also captured. Results: Eight hundred sixty-four patients with idiopathic or heritable PAH and 16 with pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis were recruited. Mutations in BMPR2 were identified in 130 patients (14.8%). Biallelic mutations in EIF2AK4 were identified in 5 patients with a clinical diagnosis of pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis. Furthermore, 9 patients with a clinical diagnosis of PAH carried biallelic EIF2AK4 mutations. These patients had a reduced transfer coefficient for carbon monoxide (K co ; 33% [interquartile range, 30%–35%] predicted) and younger age at diagnosis (29 years; interquartile range, 23–38 years) and more interlobular septal thickening and mediastinal lymphadenopathy on computed tomography of the chest compared with patients with PAH without EIF2AK4 mutations. However, radiological assessment alone could not accurately identify biallelic EIF2AK4 mutation carriers. Patients with PAH with biallelic EIF2AK4 mutations had a shorter survival. Conclusions: Biallelic EIF2AK4 mutations are found in patients classified clinically as having idiopathic and heritable PAH. These patients cannot be identified reliably by computed tomography, but a low K co and a young age at diagnosis suggests the underlying molecular diagnosis. Genetic testing can identify these misclassified patients, allowing appropriate management and early referral for lung transplantation. </jats:sec

COLEC10 is mutated in 3MC patients and regulates early craniofacial development

3MC syndrome is an autosomal recessive heterogeneous disorder with features linked to developmental abnormalities. The main features include facial dysmorphism, craniosynostosis and cleft lip/palate; skeletal structures derived from cranial neural crest cells (cNCC). We previously reported that lectin complement pathway genes COLEC11 and MASP1/3 are mutated in 3MC syndrome patients. Here we define a new gene, COLEC10, also mutated in 3MC families and present novel mutations in COLEC11 and MASP1/3 genes in a further five families. The protein products of COLEC11 and COLEC10, CL-K1 and CL-L1 respectively, form heteromeric complexes. We show COLEC10 is expressed in the base membrane of the palate during murine embryo development. We demonstrate how mutations in COLEC10 (c.25C>T; p.Arg9Ter, c.226delA; p.Gly77Glufs*66 and c.528C>G p.Cys176Trp) impair the expression and/or secretion of CL-L1 highlighting their pathogenicity. Together, these findings provide further evidence linking the lectin complement pathway and complement factors COLEC11 and COLEC10 to morphogenesis of craniofacial structures and 3MC etiology.New life Fundation for Disabled Childre

Bi-allelic loss-of-function CACNA1B mutations in progressive epilepsy-dyskinesia

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Cav2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca2+ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Cav2.2 in normal human neurodevelopment