Exome sequencing resolves apparent incidental findings and reveals further complexity of SH3TC2 variant alleles causing Charcot-Marie-Tooth neuropathy

Molecular findings in patients for whole exome sequencing and mitochondrial genome assessment

OBJECTIVE

Whole exome sequencing (WES) is becoming more widely used as a diagnostic tool in the field of medicine. In this article, we reported the diagnostic yield of WES and mitochondrial genome assessment in 2226 consecutive cases in a single clinical laboratory.

MATERIALS AND METHODS

We retrospectively analyzed consecutive WES reports from 2226 patients with various genetic disorders. WES-process was focused exclusively on the probands and aimed at a higher diagnostic capacity. We determined the diagnostic rate of WES overall and by phenotypic category, mode of inheritance, mitochondrial genome variant, and copy number variants (CNVs).

RESULTS

Among the 2226 patients who had diagnostic WES proband-only, the overall diagnostic yield of WES was 34.59% (770/2226). The highest diagnostic yield was observed in autosomal dominant disorders, at 45.58% (351/770), followed by autosomal recessive at 31.95%(246/770), X-linked disorder at 9.61%(74/770), and mitochondrial diseases at a notably lower 0.65%(5/770). The 12.21% (94/770) diagnoses were based on a total of 94 copy number variants reported from WES data. CNVs in children accounted for 67.02% of the total CNVs. While majority of the molecular diagnoses were related to nuclear genes, the inclusion of mitochondrial genome sequencing in the WES test contributed to five diagnoses. all mitochondrial diseases were identified in adults.

CONCLUSIONS

The proband-only WES provided a definitive molecular diagnosis for 34.59% of a large cohort of patients while analysis of WES simultaneously analyzed the SNVs, exons, mitochondrial genome, and CNVs, thereby improving the diagnostic yield significantly compared to the single-detection WES method; and facilitating the identification of novel candidate genes.

β-Actin G342D as a Cause of NK Cell Deficiency Impairing Lytic Synapse Termination

NK cell deficiency (NKD) occurs when an individual's major clinical immunodeficiency derives from abnormal NK cells and is associated with several genetic etiologies. Three categories of β-actin-related diseases with over 60 ACTB (β-actin) variants have previously been identified, none with a distinct NK cell phenotype. An individual with mild developmental delay, macrothrombocytopenia, and susceptibility to infections, molluscum contagiosum virus, and EBV-associated lymphoma had functional NKD for over a decade. A de novo ACTB variant encoding G342D β-actin was identified and was consistent with the individual's developmental and platelet phenotype. This novel variant also was found to have direct impact in NK cells because its expression in the human NK cell line YTS (YTS-NKD) caused increased cell spreading in lytic immune synapses created on activating surfaces. YTS-NKD cells were able to degranulate and perform cytotoxicity, but they demonstrated defective serial killing because of prolonged conjugation to the killed target cell and thus were effectively unable to terminate lytic synapses. G342D β-actin results in a novel, to our knowledge, mechanism of functional NKD via increased synaptic spreading and defective lytic synapse termination with resulting impaired serial killing, leading to overall reductions in NK cell cytotoxicity.

Exome Sequencing Implicates DGKZ , ESRRA , and GXYLT1 for Modulating Granuloma Formation in Crohn Disease

Non-caseating granulomas may indicate a more aggressive phenotype of Crohn disease (CD). Genetic associations of granulomatous CD (GCD) may help elucidate disease pathogenesis. Whole-exome sequencing was performed on peripheral blood-derived DNA from 17 pediatric patients with GCD and 19 with non-GCD (NGCD), and from an independent validation cohort of 44 GCD and 19 NGCD cases. PLINK (a tool set for whole-genome association and population-based linkage analyses) analysis was used to identify single nucleotide polymorphisms (SNPs) differentiating between groups, and subgroup allele frequencies were also compared to a public genomic database (gnomAD). The Combined Annotation Dependent Depletion scoring tool was used to predict deleteriousness of SNPs. Human leukocyte antigen (HLA) haplotype findings were compared to a control group (n = 8496). PLINK-based analysis between GCD and NGCD groups did not find consistently significant hits. gnomAD control comparisons, however, showed consistent subgroup associations with DGKZ , ESRRA , and GXYLT1 , genes that have been implicated in mammalian granulomatous inflammation. Our findings may guide future research and precision medicine.

Use of whole genome sequencing to determine the genetic basis of visceral myopathies including Prune Belly syndrome

Objectives/aims

The visceral myopathies (VM) are a group of disorders characterised by poorly contractile or acontractile smooth muscle. They manifest in both the GI and GU tracts, ranging from megacystis to Prune Belly syndrome. We aimed to apply a bespoke virtual genetic panel and describe novel variants associated with this condition using whole genome sequencing data within the Genomics England 100,000 Genomes Project.

Methods

We screened the Genomics England 100,000 Genomes Project rare diseases database for patients with VM-related phenotypes. These patients were screened for sequence variants and copy number variants (CNV) in ACTG2, ACTA2, MYH11, MYLK, LMOD1, CHRM3, MYL9, FLNA and KNCMA1 by analysing whole genome sequencing data. The identified variants were analysed using variant effect predictor online tool, and any possible segregation in other family members and novel missense mutations was modelled using in silico tools. The VM cohort was also used to perform a genome-wide variant burden test in order to identify confirm gene associations in this cohort.

Results

We identified 76 patients with phenotypes consistent with a diagnosis of VM. The range of presentations included megacystis/microcolon hypoperistalsis syndrome, Prune Belly syndrome and chronic intestinal pseudo-obstruction. Of the patients in whom we identified heterozygous ACTG2 variants, 7 had likely pathogenic variants including 1 novel likely pathogenic allele. There were 4 patients in whom we identified a heterozygous MYH11 variant of uncertain significance which leads to a frameshift and a predicted protein elongation. We identified one family in whom we found a heterozygous variant of uncertain significance in KCNMA1 which in silico models predicted to be disease causing and may explain the VM phenotype seen. We did not find any CNV changes in known genes leading to VM-related disease phenotypes. In this phenotype selected cohort, ACTG2 is the largest monogenic cause of VM-related disease accounting for 9% of the cohort, supported by a variant burden test approach, which identified ACTG2 variants as the largest contributor to VM-related phenotypes.

Conclusions

VM are a group of disorders that are not easily classified and may be given different diagnostic labels depending on their phenotype. Molecular genetic analysis of these patients is valuable as it allows precise diagnosis and aids understanding of the underlying disease manifestations. We identified ACTG2 as the most frequent genetic cause of VM. We recommend a nomenclature change to 'autosomal dominant ACTG2 visceral myopathy' for patients with pathogenic variants in ACTG2 and associated VM phenotypes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s44162-023-00012-z.

An ELF4 hypomorphic variant results in NK cell deficiency

NK cell deficiencies (NKD) are a type of primary immune deficiency in which the major immunologic abnormality affects NK cell number, maturity, or function. Since NK cells contribute to immune defense against virally infected cells, patients with NKD experience higher susceptibility to chronic, recurrent, and fatal viral infections. An individual with recurrent viral infections and mild hypogammaglobulinemia was identified to have an X-linked damaging variant in the transcription factor gene ELF4. The variant does not decrease expression but disrupts ELF4 protein interactions and DNA binding, reducing transcriptional activation of target genes and selectively impairing ELF4 function. Corroborating previous murine models of ELF4 deficiency (Elf4-/-) and using a knockdown human NK cell line, we determined that ELF4 is necessary for normal NK cell development, terminal maturation, and function. Through characterization of the NK cells of the proband, expression of the proband's variant in Elf4-/- mouse hematopoietic precursor cells, and a human in vitro NK cell maturation model, we established this ELF4 variant as a potentially novel cause of NKD.

Biology in balance: human diploid genome integrity, gene dosage, and genomic medicine

The path to completion of the functional annotation of the haploid human genome reference build, exploration of the clan genomics hypothesis, understanding human gene and genome functional biology, and gene genome and organismal evolution, is in reach.

A case of Bainbridge-Ropers syndrome with breath holding spells and intractable epilepsy: challenges in diagnosis and management

Background

Bainbridge-Ropers syndrome is caused by monoallelic ASXL3 variants on chromosome 18. Clinical features include dysmorphic facies, developmental delay, intellectual disability, autistic traits, hypotonia, failure to thrive, seizures and hyperventilation. Breath-holding spells with choreathetoid movements have been previously described.

Case presentation

We describe an 11-year old boy who has daily intractable seizures reported since birth, developmental delay, autistic features and feeding difficulties. He was eventually found to have de novo, heterozygous pathogenic variant (c.1612G > T, p.E538*) in the ASXL3 gene. He has frequent episodes of breath-holding accompanied by dystonic posturing with right leg extension and head turning without ictal EEG correlate. The breath-holding spells have been refractory to several medication trials including iron supplementation, acetazolamide, and desipramine.

Conclusions

This case represents a more severe phenotype of Bainbridge-Ropers Syndrome than previously described with refractory breath-holding spells with dystonia, intractable epilepsy, and progressive cerebral/cerebellar atrophy. Breath-holding spells cause significant morbidity, are poorly understood, and have very limited treatment options.

Genetic errors of immunity distinguish pediatric nonmalignant lymphoproliferative disorders

BACKGROUND

Pediatric nonmalignant lymphoproliferative disorders (PLPDs) are clinically and genetically heterogeneous. Long-standing immune dysregulation and lymphoproliferation in children may be life-threatening, and a paucity of data exists to guide evaluation and treatment of children with PLPD.

OBJECTIVE

The primary objective of this study was to ascertain the spectrum of genomic immunologic defects in PLPD. Secondary objectives included characterization of clinical outcomes and associations between genetic diagnoses and those outcomes.

METHODS

PLPD was defined by persistent lymphadenopathy, lymph organ involvement, or lymphocytic infiltration for more than 3 months, with or without chronic or significant Epstein-Barr virus (EBV) infection. Fifty-one subjects from 47 different families with PLPD were analyzed using whole exome sequencing.

RESULTS

Whole exome sequencing identified likely genetic errors of immunity in 51% to 62% of families (53% to 65% of affected children). Presence of a genetic etiology was associated with younger age and hemophagocytic lymphohistiocytosis. Ten-year survival for the cohort was 72.4%, and patients with viable genetic diagnoses had a higher survival rate (82%) compared to children without a genetic explanation (48%, P = .03). Survival outcomes for individuals with EBV-associated disease and no genetic explanation were particularly worse than outcomes for subjects with EBV-associated disease and a genetic explanation (17% vs 90%; P = .002). Ascertainment of a molecular diagnosis provided targetable treatment options for up to 18 individuals and led to active management changes for 12 patients.

CONCLUSIONS

PLPD defines children at high risk for mortality, and whole exome sequencing informs clinical risks and therapeutic opportunities for this diagnosis.

PRINCESS: comprehensive detection of haplotype resolved SNVs, SVs, and methylation

Long-read sequencing has been shown to have advantages in structural variation (SV) detection and methylation calling. Many studies focus either on SV, methylation, or phasing of SNV; however, only the combination of variants provides a comprehensive insight into the sample and thus enables novel findings in biology or medicine. PRINCESS is a structured workflow that takes raw sequence reads and generates a fully phased SNV, SV, and methylation call set within a few hours. PRINCESS achieves high accuracy and long phasing even on low coverage datasets and can resolve repetitive, complex medical relevant genes that often escape detection. PRINCESS is publicly available at https://github.com/MeHelmy/princess under the MIT license.

IFIH1 loss-of-function variants contribute to very early-onset inflammatory bowel disease

Genetic defects of innate immunity impairing intestinal bacterial sensing are linked to the development of Inflammatory Bowel Disease (IBD). Although much evidence supports a role of the intestinal virome in gut homeostasis, most studies focus on intestinal viral composition rather than on host intestinal viral sensitivity. To demonstrate the association between the development of Very Early Onset IBD (VEOIBD) and variants in the IFIH1 gene which encodes MDA5, a key cytosolic sensor for viral nucleic acids. Whole exome sequencing (WES) was performed in two independent cohorts of children with VEOIBD enrolled in Italy (n = 18) and USA (n = 24). Luciferase reporter assays were employed to assess MDA5 activity. An enrichment analysis was performed on IFIH1 comparing 42 VEOIBD probands with 1527 unrelated individuals without gastrointestinal or immunological issues. We identified rare, likely loss-of-function (LoF), IFIH1 variants in eight patients with VEOIBD from a combined cohort of 42 children. One subject, carrying a homozygous truncating variant resulting in complete LoF, experienced neonatal-onset, pan-gastrointestinal, IBD-like enteropathy plus multiple infectious episodes. The remaining seven subjects, affected by VEOIBD without immunodeficiency, were carriers of one LoF variant in IFIH1. Among these, two patients also carried a second hypomorphic variant, with partial function apparent when MDA5 was weakly stimulated. Furthermore, IFIH1 variants were significantly enriched in children with VEOIBD as compared to controls (p = 0.007). Complete and partial MDA5 deficiency is associated with VEOIBD with variable penetrance and expressivity, suggesting a role for impaired intestinal viral sensing in IBD pathogenesis.

Dealing with Pseudogenes in Molecular Diagnostics in the Next Generation Sequencing Era

Presence of pseudogenes is a dreadful issue in next generation sequencing (NGS), because their contamination can interfere with the detection of variants in the genuine gene and generate false positive and false negative variants.In this chapter we focus on issues related to the application of NGS strategies for analysis of genes with pseudogenes in a clinical setting. The degree to which a pseudogene impacts the ability to accurately detect and map variants in its parent gene depends on the degree of similarity (homology) with the parent gene itself. Hereby, target enrichment and mapping strategies are crucial factors to avoid "contaminating" pseudogene sequences. For target enrichment, we describe advantages and disadvantages of PCR- and capture-based strategies. For mapping strategies, we discuss crucial parameters that need to be considered to accurately distinguish sequences of functional genes from pseudogenic sequences. Finally, we discuss some examples of genes associated with Mendelian disorders, for which interesting NGS approaches are described to avoid interference with pseudogene sequences.

Unmapped exome reads implicate a role for Anelloviridae in childhood HIV-1 long-term non-progression

Human immunodeficiency virus (HIV) infection remains a significant public health burden globally. The role of viral co-infection in the rate of progression of HIV infection has been suggested but not empirically tested, particularly among children. We extracted and classified 42 viral species from whole-exome sequencing (WES) data of 813 HIV-infected children in Botswana and Uganda categorised as either long-term non-progressors (LTNPs) or rapid progressors (RPs). The Ugandan participants had a higher viral community diversity index compared to Batswana (p = 4.6 × 10-13), and viral sequences were more frequently detected among LTNPs than RPs (24% vs 16%; p = 0.008; OR, 1.9; 95% CI, 1.6-2.3), with Anelloviridae showing strong association with LTNP status (p = 3 × 10-4; q = 0.004, OR, 3.99; 95% CI, 1.74-10.25). This trend was still evident when stratified by country, sex, and sequencing platform, and after a logistic regression analysis adjusting for age, sex, country, and the sequencing platform (p = 0.02; q = 0.03; OR, 7.3; 95% CI, 1.6-40.5). Torque teno virus (TTV), which made up 95% of the Anelloviridae reads, has been associated with reduced immune activation. We identify an association between viral co-infection and prolonged AIDs-free survival status that may have utility as a biomarker of LTNP and could provide mechanistic insights to HIV progression in children, demonstrating the added value of interrogating off-target WES reads in cohort studies.

Variants in FLRT3 and SLC35E2B identified using exome sequencing in seven high myopia families from Central Europe

PURPOSE

High myopia (HM) is an eye disorder with both environmental and genetic factors involved. Many genetic factors responsible for HM were recognized worldwide, but little is known about genetic variants underlying HM in Central Europe. Thus, the aim of this study was to identify rare sequence variants involved in HM in families from Central Europe to better understand the genetic basis of HM.

MATERIALS AND METHODS

We assessed 17 individuals from 7 unrelated Central European families with hereditary HM using exome sequencing (ES). Segregation of selected variants in other available family members was performed using Sanger sequencing.

RESULTS

Detected 73 rare variants were selected for verification. We observed 2 missense variants, c.938C>T in SLC35E2B - encoding solute carrier family 35 member E2B, and c.1642G>C in FLRT3 - encoding fibronectin leucine rich transmembrane protein, segregating with HM in one family.

CONCLUSIONS

FLRT3 ​and/or ​SLC35E2B ​could represent disease candidate genes and identified sequence variants might be responsible for HM in the studied family.

Rare deleterious germline variants and risk of lung cancer

Recent studies suggest that rare variants exhibit stronger effect sizes and might play a crucial role in the etiology of lung cancers (LC). Whole exome plus targeted sequencing of germline DNA was performed on 1045 LC cases and 885 controls in the discovery set. To unveil the inherited causal variants, we focused on rare and predicted deleterious variants and small indels enriched in cases or controls. Promising candidates were further validated in a series of 26,803 LCs and 555,107 controls. During discovery, we identified 25 rare deleterious variants associated with LC susceptibility, including 13 reported in ClinVar. Of the five validated candidates, we discovered two pathogenic variants in known LC susceptibility loci, ATM p.V2716A (Odds Ratio [OR] 19.55, 95%CI 5.04-75.6) and MPZL2 p.I24M frameshift deletion (OR 3.88, 95%CI 1.71-8.8); and three in novel LC susceptibility genes, POMC c.*28delT at 3' UTR (OR 4.33, 95%CI 2.03-9.24), STAU2 p.N364M frameshift deletion (OR 4.48, 95%CI 1.73-11.55), and MLNR p.Q334V frameshift deletion (OR 2.69, 95%CI 1.33-5.43). The potential cancer-promoting role of selected candidate genes and variants was further supported by endogenous DNA damage assays. Our analyses led to the identification of new rare deleterious variants with LC susceptibility. However, in-depth mechanistic studies are still needed to evaluate the pathogenic effects of these specific alleles.

Absent B cells, agammaglobulinemia, and hypertrophic cardiomyopathy in folliculin-interacting protein 1 deficiency

Agammaglobulinemia is the most profound primary antibody deficiency that can occur due to an early termination of B-cell development. We here investigated 3 novel patients, including the first known adult, from unrelated families with agammaglobulinemia, recurrent infections, and hypertrophic cardiomyopathy (HCM). Two of them also presented with intermittent or severe chronic neutropenia. We identified homozygous or compound-heterozygous variants in the gene for folliculin interacting protein 1 (FNIP1), leading to loss of the FNIP1 protein. B-cell metabolism, including mitochondrial numbers and activity and phosphatidylinositol 3-kinase/AKT pathway, was impaired. These defects recapitulated the Fnip1-/- animal model. Moreover, we identified either uniparental disomy or copy-number variants (CNVs) in 2 patients, expanding the variant spectrum of this novel inborn error of immunity. The results indicate that FNIP1 deficiency can be caused by complex genetic mechanisms and support the clinical utility of exome sequencing and CNV analysis in patients with broad phenotypes, including agammaglobulinemia and HCM. FNIP1 deficiency is a novel inborn error of immunity characterized by early and severe B-cell development defect, agammaglobulinemia, variable neutropenia, and HCM. Our findings elucidate a functional and relevant role of FNIP1 in B-cell development and metabolism and potentially neutrophil activity.

Identifying a Minor Histocompatibility Antigen in Mauritian Cynomolgus Macaques Encoded by APOBEC3C

Allogeneic hematopoietic stem cell transplants can lead to dramatic reductions in human immunodeficiency virus (HIV) reservoirs. This effect is partially mediated by donor T cells recognizing lymphocyte-expressed minor histocompatibility antigens (mHAgs). The potential to mark malignant and latently infected cells for destruction makes mHAgs attractive targets for cellular immunotherapies. However, testing such HIV reservoir reduction strategies will likely require preclinical studies in non-human primates (NHPs). In this study, we used a combination of alloimmunization, whole exome sequencing, and bioinformatics to identify an mHAg in Mauritian cynomolgus macaques (MCMs). We mapped the minimal optimal epitope to a 10-mer peptide (SW10) in apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3C (APOBEC3C) and determined the major histocompatibility complex class I restriction element as Mafa-A1∗063, which is expressed in almost 90% of MCMs. APOBEC3C SW10-specific CD8+ T cells recognized immortalized B cells but not fibroblasts from an mHAg-positive MCM. These results provide a framework for identifying mHAgs in a non-transplant setting and suggest that APOBEC3C SW10 could be used as a model antigen to test mHAg-targeted therapies in NHPs.

Human NK cell deficiency as a result of biallelic mutations in MCM10

Human natural killer cell deficiency (NKD) arises from inborn errors of immunity that lead to impaired NK cell development, function, or both. Through the understanding of the biological perturbations in individuals with NKD, requirements for the generation of terminally mature functional innate effector cells can be elucidated. Here, we report a cause of NKD resulting from compound heterozygous mutations in minichromosomal maintenance complex member 10 (MCM10) that impaired NK cell maturation in a child with fatal susceptibility to CMV. MCM10 has not been previously associated with monogenic disease and plays a critical role in the activation and function of the eukaryotic DNA replisome. Through evaluation of patient primary fibroblasts, modeling patient mutations in fibroblast cell lines, and MCM10 knockdown in human NK cell lines, we have shown that loss of MCM10 function leads to impaired cell cycle progression and induction of DNA damage-response pathways. By modeling MCM10 deficiency in primary NK cell precursors, including patient-derived induced pluripotent stem cells, we further demonstrated that MCM10 is required for NK cell terminal maturation and acquisition of immunological system function. Together, these data define MCM10 as an NKD gene and provide biological insight into the requirement for the DNA replisome in human NK cell maturation and function.

Identification of potential crucial genes associated with the pathogenesis and prognosis of pancreatic adenocarcinoma

Pancreatic adenocarcinoma (PAAD) is a type of malignant tumor with the highest mortality rate among all neoplasms worldwide, and its exact pathogenesis is still poorly understood. Timely diagnosis and treatment are of great importance in order to decrease the mortality rate of PAAD. Therefore, identifying new biomarkers for diagnosis and prognosis is essential to enable early detection of PAAD and to improve the overall survival (OS) rate. In order to screen and integrate differentially expressed genes (DEGs) between PAAD and normal tissues, a total of seven datasets were downloaded from the Gene Expression Omnibus database and the ‘limma’ and ‘robustrankggreg’ packages in R software were used. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis of the DEGs was performed using the Database for Annotation, Visualization and Integrated Discovery website, and the protein-protein interaction network analysis was performed using the Search Tool for the Retrieval of Interacting Genes/Proteins database. A gene prognostic signature was constructed using the Cox regression model. A total of 10 genes (CDK1, CCNB1, CDC20, ASPM, UBE2C, TPX2, TOP2A, NUSAP1, KIF20A and DLGAP5) that may be associated with pancreatic adenocarcinoma were identified. According to the differentially expressed genes in The Cancer Genome Atlas, the present study set up four prognostic signatures (matrix metalloproteinase 12, sodium voltage-gated channel α subunit 11, tetraspanin 1 and SH3 domain and tetratricopeptide repeats-containing 2), which effectively predicted OS. The hub genes that were highly associated with the occurrence, development and prognosis of PAAD were identified, which may be helpful to further understand the molecular basis of pancreatic cancer and guide the synthesis of drugs for PPAD.

Functional biology of the Steel syndrome founder allele and evidence for clan genomics derivation of COL27A1 pathogenic alleles worldwide

Previously we reported the identification of a homozygous COL27A1 (c.2089G>C; p.Gly697Arg) missense variant and proposed it as a founder allele in Puerto Rico segregating with Steel syndrome (STLS, MIM #615155); a rare osteochondrodysplasia characterized by short stature, congenital bilateral hip dysplasia, carpal coalitions, and scoliosis. We now report segregation of this variant in five probands from the initial clinical report defining the syndrome and an additional family of Puerto Rican descent with multiple affected adult individuals. We modeled the orthologous variant in murine Col27a1 and found it recapitulates some of the major Steel syndrome associated skeletal features including reduced body length, scoliosis, and a more rounded skull shape. Characterization of the in vivo murine model shows abnormal collagen deposition in the extracellular matrix and disorganization of the proliferative zone of the growth plate. We report additional COL27A1 pathogenic variant alleles identified in unrelated consanguineous Turkish kindreds suggesting Clan Genomics and identity-by-descent homozygosity contributing to disease in this population. The hypothesis that carrier states for this autosomal recessive osteochondrodysplasia may contribute to common complex traits is further explored in a large clinical population cohort. Our findings augment our understanding of COL27A1 biology and its role in skeletal development; and expand the functional allelic architecture in this gene underlying both rare and common disease phenotypes.

Recurrent arginine substitutions in the ACTG2 gene are the primary driver of disease burden and severity in visceral myopathy

Visceral myopathy with abnormal intestinal and bladder peristalsis includes a clinical spectrum with megacystis-microcolon intestinal hypoperistalsis syndrome and chronic intestinal pseudo-obstruction. The vast majority of cases are caused by dominant variants in ACTG2; however, the overall genetic architecture of visceral myopathy has not been well-characterized. We ascertained 53 families, with visceral myopathy based on megacystis, functional bladder/gastrointestinal obstruction, or microcolon. A combination of targeted ACTG2 sequencing and exome sequencing was used. We report a molecular diagnostic rate of 64% (34/53), of which 97% (33/34) is attributed to ACTG2. Strikingly, missense mutations in five conserved arginine residues involving CpG dinucleotides accounted for 49% (26/53) of disease in the cohort. As a group, the ACTG2-negative cases had a more favorable clinical outcome and more restricted disease. Within the ACTG2-positive group, poor outcomes (characterized by total parenteral nutrition dependence, death, or transplantation) were invariably due to one of the arginine missense alleles. Analysis of specific residues suggests a severity spectrum of p.Arg178>p.Arg257>p.Arg40 along with other less-frequently reported sites p.Arg63 and p.Arg211. These results provide genotype-phenotype correlation for ACTG2-related disease and demonstrate the importance of arginine missense changes in visceral myopathy.

DYRK1A-related intellectual disability: a syndrome associated with congenital anomalies of the kidney and urinary tract

PURPOSE

Haploinsufficiency of DYRK1A causes a recognizable clinical syndrome. The goal of this paper is to investigate congenital anomalies of the kidney and urinary tract (CAKUT) and genital defects (GD) in patients with DYRK1A variants.

METHODS

A large database of clinical exome sequencing (ES) was queried for de novo DYRK1A variants and CAKUT/GD phenotypes were characterized. Xenopus laevis (frog) was chosen as a model organism to assess Dyrk1a's role in renal development.

RESULTS

Phenotypic details and variants of 19 patients were compiled after an initial observation that one patient with a de novo pathogenic variant in DYRK1A had GD. CAKUT/GD data were available from 15 patients, 11 of whom presented with CAKUT/GD. Studies in Xenopus embryos demonstrated that knockdown of Dyrk1a, which is expressed in forming nephrons, disrupts the development of segments of embryonic nephrons, which ultimately give rise to the entire genitourinary (GU) tract. These defects could be rescued by coinjecting wild-type human DYRK1A RNA, but not with DYRK1AR205* or DYRK1AL245R RNA.

CONCLUSION

Evidence supports routine GU screening of all individuals with de novo DYRK1A pathogenic variants to ensure optimized clinical management. Collectively, the reported clinical data and loss-of-function studies in Xenopus substantiate a novel role for DYRK1A in GU development.

Mutations in ANKLE2, a ZIKA Virus Target, Disrupt an Asymmetric Cell Division Pathway in Drosophila Neuroblasts to Cause Microcephaly

The apical Par complex, which contains atypical protein kinase C (aPKC), Bazooka (Par-3), and Par-6, is required for establishing polarity during asymmetric division of neuroblasts in Drosophila, and its activity depends on L(2)gl. We show that loss of Ankle2, a protein associated with microcephaly in humans and known to interact with Zika protein NS4A, reduces brain volume in flies and impacts the function of the Par complex. Reducing Ankle2 levels disrupts endoplasmic reticulum (ER) and nuclear envelope morphology, releasing the kinase Ballchen-VRK1 into the cytosol. These defects are associated with reduced phosphorylation of aPKC, disruption of Par-complex localization, and spindle alignment defects. Importantly, removal of one copy of ballchen or l(2)gl suppresses Ankle2 mutant phenotypes and restores viability and brain size. Human mutational studies implicate the above-mentioned genes in microcephaly and motor neuron disease. We suggest that NS4A, ANKLE2, VRK1, and LLGL1 define a pathway impinging on asymmetric determinants of neural stem cell division.

Genotype-Phenotype Correlations in Charcot-Marie-Tooth Disease Due to MTMR2 Mutations and Implications in Membrane Trafficking

Charcot-Marie-Tooth type 4 (CMT4) is an autosomal recessive severe form of neuropathy with genetic heterogeneity. CMT4B1 is caused by mutations in the myotubularin-related 2 (MTMR2) gene and as a member of the myotubularin family, the MTMR2 protein is crucial for the modulation of membrane trafficking. To enable future clinical trials, we performed a detailed review of the published cases with MTMR2 mutations and describe four novel cases identified through whole-exome sequencing (WES). The four unrelated families harbor novel homozygous mutations in MTMR2 (NM_016156, Family 1: c.1490dupC; p.Phe498IlefsTer2; Family 2: c.1479+1G>A; Family 3: c.1090C>T; p.Arg364Ter; Family 4: c.883C>T; p.Arg295Ter) and present with CMT4B1-related severe early-onset motor and sensory neuropathy, generalized muscle atrophy, facial and bulbar weakness, and pes cavus deformity. The clinical description of the new mutations reported here overlap with previously reported CMT4B1 phenotypes caused by mutations in the phosphatase domain of MTMR2, suggesting that nonsense MTMR2 mutations, which are predicted to result in loss or disruption of the phosphatase domain, are associated with a severe phenotype and loss of independent ambulation by the early twenties. Whereas the few reported missense mutations and also those truncating mutations occurring at the C-terminus after the phosphatase domain cause a rather mild phenotype and patients were still ambulatory above the age 30 years. Charcot-Marie-Tooth neuropathy and Centronuclear Myopathy causing mutations have been shown to occur in proteins involved in membrane remodeling and trafficking pathway mediated by phosphoinositides. Earlier studies have showing the rescue of MTM1 myopathy by MTMR2 overexpression, emphasize the importance of maintaining the phosphoinositides equilibrium and highlight a potential compensatory mechanism amongst members of this pathway. This proved that the regulation of expression of these proteins involved in the membrane remodeling pathway may compensate each other's loss- or gain-of-function mutations by restoring the phosphoinositides equilibrium. This provides a potential therapeutic strategy for neuromuscular diseases resulting from mutations in the membrane remodeling pathway.

A Patient with Berardinelli-Seip Syndrome, Novel AGPAT2 Splicesite Mutation and Concomitant Development of Non-diabetic Polyneuropathy

Primary polyneuropathy in the context of Seip-Berardinelli type 1 seipinopathy, or congenital generalized lipodystrophy type 1 (CGL1) has not been previously reported. We report the case history of a 27 year old female CGL1 patient presenting with an unusual additional development of non-diabetic peripheral neuropathy and learning disabilities in early adolescence. Whole exome sequencing (WES) of the patient genome identified a novel variant, homozygous for a 52 bp intronic deletion in the AGPAT2 locus, coding for 1-acylglycerol-3-phosphate O-acyltransferase 2, which is uniquely associated with CGL1 seipinopathies, with no molecular evidence for dual diagnosis. Functional studies using RNA isolated from patient peripheral blood leucocytes showed abnormal RNA splicing resulting in the loss of 25 amino acids from the patient AGPAT2 protein coding sequence. Stability and transcription levels for the misspliced AGPAT2 mRNA in our patient nonetheless remained normal. Any AGPAT2 protein produced in our patient is therefore likely to be dysfunctional. However, formal linkage of this deletion to the neuropathy observed remains to be shown. The classical clinical presentation of a patient with AGPAT2-associated lipodystrophy shows normal cognition and no development of polyneuropathy. Cognitive disabilities and polyneuropathy are features associated exclusively with clinical CGL type 2 arising from seipin (BSCL2) gene mutations. This case study suggests that in some genetic contexts, AGPAT2 mutations can also produce phenotypes with primary polyneuropathy.

Heterozygous CTNNB1 and TBX4 variants in a patient with abnormal lung growth, pulmonary hypertension, microcephaly, and spasticity

The canonical wingless (Wnt) and fibroblast growth factor (FGF) signaling pathways involving CTNNB1 and TBX4, respectively, are crucial for the regulation of human development. Perturbations of these pathways and disruptions from biological homeostasis have been associated with abnormal morphogenesis of multiple organs, including the lung. The aim of this study was to identify the underlying genetic cause of abnormal lung growth, pulmonary hypertension (PAH), severe microcephaly, and muscle spasticity in a full-term newborn, who died at 4 months of age due to progressively worsening PAH and respiratory failure. Family trio exome sequencing showed a de novo heterozygous nonsense c.1603C>T (p.Arg535*) variant in CTNNB1 and a paternally inherited heterozygous missense c.1198G>A (p.Glu400Lys) variant in TBX4, both predicted to be likely deleterious. We expand the phenotypic spectrum associated with CTNNB1 and TBX4 variants and indicate that they could act synergistically to produce a distinct more severe phenotype. Our findings further support a recently proposed complex compound inheritance model in lethal lung developmental diseases and the contention that dual molecular diagnoses can parsimoniously explain blended phenotypes.

Implication of the SH3TC2 gene in Charcot-Marie-Tooth disease associated with deafness and/or scoliosis: Illustration with four new pathogenic variants

The autosomal recessive demyelinating form of Charcot-Marie-Tooth can be due to SH3TC2 gene pathogenic variants (CMT4C, AR-CMTde-SH3TC2). We report on a series of 13 patients with AR-CMTde-SH3TC2 among a French cohort of 350 patients suffering from all type of inheritance peripheral neuropathy. The SH3TC2 gene appeared to be the most frequently mutated gene for demyelinating neuropathy in this series by NGS. Four new pathogenic variants have been identified: two nonsense variants (p.(Tyr970*), p.(Trp1199*)) and two missense variants (p.(Leu1126Pro), p.(Ala1206Asp)). The recurrent variant p.Arg954* was present in 62%, and seems to be a founder mutation. The phenotype is fairly homogeneous, as all these patients, except the youngest ones, presented scoliosis and/or hearing loss.

Comparative analysis of de novo assemblers for variation discovery in personal genomes

Current variant discovery approaches often rely on an initial read mapping to the reference sequence. Their effectiveness is limited by the presence of gaps, potential misassemblies, regions of duplicates with a high-sequence similarity and regions of high-sequence divergence in the reference. Also, mapping-based approaches are less sensitive to large INDELs and complex variations and provide little phase information in personal genomes. A few de novo assemblers have been developed to identify variants through direct variant calling from the assembly graph, micro-assembly and whole-genome assembly, but mainly for whole-genome sequencing (WGS) data. We developed SGVar, a de novo assembly workflow for haplotype-based variant discovery from whole-exome sequencing (WES) data. Using simulated human exome data, we compared SGVar with five variation-aware de novo assemblers and with BWA-MEM together with three haplotype- or local de novo assembly-based callers. SGVar outperforms the other assemblers in sensitivity and tolerance of sequencing errors. We recapitulated the findings on whole-genome and exome data from a Utah residents with Northern and Western European ancestry (CEU) trio, showing that SGVar had high sensitivity both in the highly divergent human leukocyte antigen (HLA) region and in non-HLA regions of chromosome 6. In particular, SGVar is robust to sequencing error, k-mer selection, divergence level and coverage depth. Unlike mapping-based approaches, SGVar is capable of resolving long-range phase and identifying large INDELs from WES, more prominently from WGS. We conclude that SGVar represents an ideal platform for WES-based variant discovery in highly divergent regions and across the whole genome.

Mutations in the mitochondrial ribosomal protein MRPS22 lead to primary ovarian insufficiency

Primary ovarian insufficiency (POI) is characterized by amenorrhea and loss or dysfunction of ovarian follicles prior to the age of 40. POI has been associated with autosomal recessive mutations in genes involving hormonal signaling and folliculogenesis, however, the genetic etiology of POI most often remains unknown. Here we report MRPS22 homozygous missense variants c.404G>A (p.R135Q) and c.605G>A (p.R202H) identified in four females from two independent consanguineous families as a novel genetic cause of POI in adolescents. Both missense mutations identified in MRPS22 are rare, occurred in highly evolutionarily conserved residues, and are predicted to be deleterious to protein function. In contrast to prior reports of mutations in MRPS22 associated with severe mitochondrial disease, the POI phenotype is far less severe. Consistent with this genotype-phenotype correlation, mitochondrial defects in oxidative phosphorylation or rRNA levels were not detected in fibroblasts derived from the POI patients, suggesting a non-bioenergetic or tissue-specific mitochondrial defect. Furthermore, we demonstrate in a Drosophila model that mRpS22 deficiency specifically in somatic cells of the ovary had no effect on fertility, whereas flies with mRpS22 deficiency specifically in germ cells were infertile and agametic, demonstrating a cell autonomous requirement for mRpS22 in germ cell development. These findings collectively identify that MRPS22, a component of the small mitochondrial ribosome subunit, is critical for ovarian development and may therefore provide insight into the pathophysiology and treatment of ovarian dysfunction.

Complex Compound Inheritance of Lethal Lung Developmental Disorders Due to Disruption of the TBX-FGF Pathway

Primary defects in lung branching morphogenesis, resulting in neonatal lethal pulmonary hypoplasias, are incompletely understood. To elucidate the pathogenetics of human lung development, we studied a unique collection of samples obtained from deceased individuals with clinically and histopathologically diagnosed interstitial neonatal lung disorders: acinar dysplasia (n = 14), congenital alveolar dysplasia (n = 2), and other lethal lung hypoplasias (n = 10). We identified rare heterozygous copy-number variant deletions or single-nucleotide variants (SNVs) involving TBX4 (n = 8 and n = 2, respectively) or FGF10 (n = 2 and n = 2, respectively) in 16/26 (61%) individuals. In addition to TBX4, the overlapping ∼2 Mb recurrent and nonrecurrent deletions at 17q23.1q23.2 identified in seven individuals with lung hypoplasia also remove a lung-specific enhancer region. Individuals with coding variants involving either TBX4 or FGF10 also harbored at least one non-coding SNV in the predicted lung-specific enhancer region, which was absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies. The occurrence of rare coding variants involving TBX4 or FGF10 with the putative hypomorphic non-coding SNVs implies a complex compound inheritance of these pulmonary hypoplasias. Moreover, they support the importance of TBX4-FGF10-FGFR2 epithelial-mesenchymal signaling in human lung organogenesis and help to explain the histopathological continuum observed in these rare lethal developmental disorders of the lung.

Rare Variants in Known Susceptibility Loci and Their Contribution to Risk of Lung Cancer

BACKGROUND

Genome-wide association studies are widely used to map genomic regions contributing to lung cancer (LC) susceptibility, but they typically do not identify the precise disease-causing genes/variants. To unveil the inherited genetic variants that cause LC, we performed focused exome-sequencing analyses on genes located in 121 genome-wide association study-identified loci previously implicated in the risk of LC, chronic obstructive pulmonary disease, pulmonary function level, and smoking behavior.

METHODS

Germline DNA from 260 case patients with LC and 318 controls were sequenced by utilizing VCRome 2.1 exome capture. Filtering was based on enrichment of rare and potential deleterious variants in cases (risk alleles) or controls (protective alleles). Allelic association analyses of single-variant and gene-based burden tests of multiple variants were performed. Promising candidates were tested in two independent validation studies with a total of 1773 case patients and 1123 controls.

RESULTS

We identified 48 rare variants with deleterious effects in the discovery analysis and validated 12 of the 43 candidates that were covered in the validation platforms. The top validated candidates included one well-established truncating variant, namely, BRCA2, DNA repair associated gene (BRCA2) K3326X (OR = 2.36, 95% confidence interval [CI]: 1.38-3.99), and three newly identified variations, namely, lymphotoxin beta gene (LTB) p.Leu87Phe (OR = 7.52, 95% CI: 1.01-16.56), prolyl 3-hydroxylase 2 gene (P3H2) p.Gln185His (OR = 5.39, 95% CI: 0.75-15.43), and dishevelled associated activator of morphogenesis 2 gene (DAAM2) p.Asp762Gly (OR = 0.25, 95% CI: 0.10-0.79). Burden tests revealed strong associations between zinc finger protein 93 gene (ZNF93), DAAM2, bromodomain containing 9 gene (BRD9), and the gene LTB and LC susceptibility.

CONCLUSION

Our results extend the catalogue of regions associated with LC and highlight the importance of germline rare coding variants in LC susceptibility.

A comprehensive clinical and genetic study in 127 patients with ID in Kinshasa, DR Congo

Pathogenic variants account for 4 to 41% of patients with intellectual disability (ID) or developmental delay (DD). In Sub-Saharan Africa, the prevalence of ID is thought to be higher, but data in Central Africa are limited to some case reports. In addition, clinical descriptions of some syndromes are not available for this population. This study aimed at providing an estimate for the fraction of ID/DD for which an underlying etiological genetic cause may be elucidated and provide insights into their clinical presentation in special institutions in a Central African country. A total of 127 patients (33 females and 94 males, mean age 10.03 ± 4.68 years), were recruited from six institutions across Kinshasa. A clinical diagnosis was achieved in 44 but molecular confirmation was achieved in 21 of the 22 patients with expected genetic defect (95% clinical sensitivity). Identified diseases included Down syndrome (15%), submicroscopic copy number variants (9%), aminoacylase deficiency (0.8%), Partington syndrome in one patient (0.8%) and his similarly affected brother, X-linked syndromic Mental Retardation type 33 (0.8%), and two conditions without clear underlying molecular genetic etiologies (Oculo-Auriculo-Vertebral and Amniotic Bands Sequence). We have shown that genetic etiologies, similar to those reported in Caucasian subjects, are a common etiologic cause of ID in African patients from Africa. We have confirmed the diagnostic utility of clinical characterization prior to genetic testing. Finally, our clinical descriptions provide insights into the presentation of these genetic diseases in African patients.

Mutation in the intracellular chloride channel CLCC1 associated with autosomal recessive retinitis pigmentosa

We identified a homozygous missense alteration (c.75C>A, p.D25E) in CLCC1, encoding a presumptive intracellular chloride channel highly expressed in the retina, associated with autosomal recessive retinitis pigmentosa (arRP) in eight consanguineous families of Pakistani descent. The p.D25E alteration decreased CLCC1 channel function accompanied by accumulation of mutant protein in granules within the ER lumen, while siRNA knockdown of CLCC1 mRNA induced apoptosis in cultured ARPE-19 cells. TALEN KO in zebrafish was lethal 11 days post fertilization. The depressed electroretinogram (ERG) cone response and cone spectral sensitivity of 5 dpf KO zebrafish and reduced eye size, retinal thickness, and expression of rod and cone opsins could be rescued by injection of wild type CLCC1 mRNA. Clcc1+/- KO mice showed decreased ERGs and photoreceptor number. Together these results strongly suggest that intracellular chloride transport by CLCC1 is a critical process in maintaining retinal integrity, and CLCC1 is crucial for survival and function of retinal cells.

Genetic and mechanistic diversity in pediatric hemophagocytic lymphohistiocytosis

The HLH-2004 criteria are used to diagnose hemophagocytic lymphohistiocytosis (HLH), yet concern exists for their misapplication, resulting in suboptimal treatment of some patients. We sought to define the genomic spectrum and associated outcomes of a diverse cohort of children who met the HLH-2004 criteria. Genetic testing was performed clinically or through research-based whole-exome sequencing. Clinical metrics were analyzed with respect to genomic results. Of 122 subjects enrolled over the course of 17 years, 101 subjects received genetic testing. Biallelic familial HLH (fHLH) gene defects were identified in only 19 (19%) and correlated with presentation at younger than 1 year of age (P < .0001). Digenic fHLH variants were observed but lacked statistical support for disease association. In 28 (58%) of 48 subjects, research whole-exome sequencing analyses successfully identified likely molecular explanations, including underlying primary immunodeficiency diseases, dysregulated immune activation and proliferation disorders, and potentially novel genetic conditions. Two-thirds of patients identified by the HLH-2004 criteria had underlying etiologies for HLH, including genetic defects, autoimmunity, and malignancy. Overall survival was 45%, and increased mortality correlated with HLH triggered by infection or malignancy (P < .05). Differences in survival did not correlate with genetic profile or extent of therapy. HLH should be conceptualized as a phenotype of critical illness characterized by toxic activation of immune cells from different underlying mechanisms. In most patients with HLH, targeted sequencing of fHLH genes remains insufficient for identifying pathogenic mechanisms. Whole-exome sequencing, however, may identify specific therapeutic opportunities and affect hematopoietic stem cell transplantation options for these patients.

Heterozygous Truncating Variants in POMP Escape Nonsense-Mediated Decay and Cause a Unique Immune Dysregulatory Syndrome

The proteasome processes proteins to facilitate immune recognition and host defense. When inherently defective, it can lead to aberrant immunity resulting in a dysregulated response that can cause autoimmunity and/or autoinflammation. Biallelic or digenic loss-of-function variants in some of the proteasome subunits have been described as causing a primary immunodeficiency disease that manifests as a severe dysregulatory syndrome: chronic atypical neutrophilic dermatosis with lipodystrophy and elevated temperature (CANDLE). Proteasome maturation protein (POMP) is a chaperone for proteasome assembly and is critical for the incorporation of catalytic subunits into the proteasome. Here, we characterize and describe POMP-related autoinflammation and immune dysregulation disease (PRAID) discovered in two unrelated individuals with a unique constellation of early-onset combined immunodeficiency, inflammatory neutrophilic dermatosis, and autoimmunity. We also begin to delineate a complex genetic mechanism whereby de novo heterozygous frameshift variants in the penultimate exon of POMP escape nonsense-mediated mRNA decay (NMD) and result in a truncated protein that perturbs proteasome assembly by a dominant-negative mechanism. To our knowledge, this mechanism has not been reported in any primary immunodeficiencies, autoinflammatory syndromes, or autoimmune diseases. Here, we define a unique hypo- and hyper-immune phenotype and report an immune dysregulation syndrome caused by frameshift mutations that escape NMD.

Heterozygous Truncating Variants in POMP Escape Nonsense-Mediated Decay and Cause a Unique Immune Dysregulatory Syndrome

The proteasome processes proteins to facilitate immune recognition and host defense. When inherently defective, it can lead to aberrant immunity resulting in a dysregulated response that can cause autoimmunity and/or autoinflammation. Biallelic or digenic loss-of-function variants in some of the proteasome subunits have been described as causing a primary immunodeficiency disease that manifests as a severe dysregulatory syndrome: chronic atypical neutrophilic dermatosis with lipodystrophy and elevated temperature (CANDLE). Proteasome maturation protein (POMP) is a chaperone for proteasome assembly and is critical for the incorporation of catalytic subunits into the proteasome. Here, we characterize and describe POMP-related autoinflammation and immune dysregulation disease (PRAID) discovered in two unrelated individuals with a unique constellation of early-onset combined immunodeficiency, inflammatory neutrophilic dermatosis, and autoimmunity. We also begin to delineate a complex genetic mechanism whereby de novo heterozygous frameshift variants in the penultimate exon of POMP escape nonsense-mediated mRNA decay (NMD) and result in a truncated protein that perturbs proteasome assembly by a dominant-negative mechanism. To our knowledge, this mechanism has not been reported in any primary immunodeficiencies, autoinflammatory syndromes, or autoimmune diseases. Here, we define a unique hypo- and hyper-immune phenotype and report an immune dysregulation syndrome caused by frameshift mutations that escape NMD.

Whole-exome sequencing of sickle cell disease patients with hyperhemolysis syndrome suggests a role for rare variation in disease predisposition

BACKGROUND

Hyperhemolysis syndrome (HHS) is an uncommon, but life-threatening, transfusion-related complication of red blood cell transfusion. HHS has predominantly been described in patients with sickle cell disease (SCD) and is difficult to diagnose and treat. The pathogenesis of HHS, including its occurrence in only a subset of apparently susceptible individuals, is poorly understood. We undertook whole-exome sequencing (WES) of 12 SCD-HHS patients to identify shared genetic variants that might be relevant to the development of HHS.

STUDY DESIGN AND METHODS

DNA from adults with SCD having at least one previous episode of HHS were subject to WES. High-quality variants were passed through a series of bioinformatics filters to identify variants that were uncommon among African populations represented in public databases. Recurrent, putative loss-of-function variants occurring in biologically plausible genes were prioritized and then genotyped in a larger, ancestry-matched cohort of non-HHS controls.

RESULTS

A rare, heterozygous stop-gain variant (p.Glu210Ter) in MBL2 was significantly enriched among HHS cases (p = 0.002). This variant is predicted to result in a premature termination codon that escapes nonsense-mediated mRNA decay, potentially leading to a novel phenotype. We also observed a complex insertion-deletion variant in the final exon of KLRC3 that was enriched among cases (p = 0.0019), although neither variant was found among seven pediatric SCD-HHS patients.

CONCLUSION

Our results suggest a potential role for rare genetic defects in the development of HHS among adult SCD patients. Such enriched variants may ultimately be useful for identifying high-risk individuals and informing therapeutic approaches in HHS.

Mutations in PI3K110δ cause impaired natural killer cell function partially rescued by rapamycin treatment

BACKGROUND

Heterozygous gain-of-function mutations in PI3K110δ lead to lymphadenopathy, lymphoid hyperplasia, EBV and cytomegalovirus viremia, and sinopulmonary infections.

OBJECTIVE

The known role of natural killer (NK) cell function in the control of EBV and cytomegalovirus prompted us to investigate the functional and phenotypic effects of PI3K110δ mutations on NK cell subsets and cytotoxic function.

METHODS

Mutations in patients were identified by using whole-exome or targeted sequencing. We performed NK cell phenotyping and functional analysis of patients' cells using flow cytometry, standard Cr⁵¹ cytotoxicity assays, and quantitative confocal microscopy.

RESULTS

PI3K110δ mutations led to an altered NK cell developmental phenotype and cytotoxic dysfunction. Impaired NK cell cytotoxicity was due to decreased conjugate formation with susceptible target cells and abrogated activation of cell machinery required for target cell killing. These defects were restored partially after initiation of treatment with rapamycin in 3 patients.

CONCLUSION

We describe novel NK cell functional deficiency caused by PI3K110δ mutation, which is a likely contributor to the severe viremia observed in these patients. Rapamycin treatment partially restores NK cell function, providing a further rationale for its use in patients with this disease.

SVachra: a tool to identify genomic structural variation in mate pair sequencing data containing inward and outward facing reads

Background

Characterization of genomic structural variation (SV) is essential to expanding the research and clinical applications of genome sequencing. Reliance upon short DNA fragment paired end sequencing has yielded a wealth of single nucleotide variants and internal sequencing read insertions-deletions, at the cost of limited SV detection. Multi-kilobase DNA fragment mate pair sequencing has supplemented the void in SV detection, but introduced new analytic challenges requiring SV detection tools specifically designed for mate pair sequencing data. Here, we introduce SVachra – Structural Variation Assessment of CHRomosomal Aberrations, a breakpoint calling program that identifies large insertions-deletions, inversions, inter- and intra-chromosomal translocations utilizing both inward and outward facing read types generated by mate pair sequencing.

Results

We demonstrate SVachra’s utility by executing the program on large-insert (Illumina Nextera) mate pair sequencing data from the personal genome of a single subject (HS1011). An additional data set of long-read (Pacific BioSciences RSII) was also generated to validate SV calls from SVachra and other comparison SV calling programs. SVachra exhibited the highest validation rate and reported the widest distribution of SV types and size ranges when compared to other SV callers.

Conclusions

SVachra is a highly specific breakpoint calling program that exhibits a more unbiased SV detection methodology than other callers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4021-y) contains supplementary material, which is available to authorized users.

Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort

We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor–Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5% (82% for single-exonic and 88% for multi-exonic calls) and precision of 78% (53% single-exonic and 96% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17–50% of pathogenic CNVs in different disease cohorts where 7.1–11% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

A Children's Oncology Group and TARGET initiative exploring the genetic landscape of Wilms tumor

We performed genome-wide sequencing and analyzed mRNA and miRNA expression, DNA copy number, and DNA methylation in 117 Wilms tumors, followed by targeted sequencing of 651 Wilms tumors. In addition to genes previously implicated in Wilms tumors (WT1, CTNNB1, AMER1, DROSHA, DGCR8, XPO5, DICER1, SIX1, SIX2, MLLT1, MYCN, and TP53), we identified mutations in genes not previously recognized as recurrently involved in Wilms tumors, the most frequent being BCOR, BCORL1, NONO, MAX, COL6A3, ASXL1, MAP3K4, and ARID1A. DNA copy number changes resulted in recurrent 1q gain, MYCN amplification, LIN28B gain, and MIRLET7A loss. Unexpected germline variants involved PALB2 and CHEK2. Integrated analyses support two major classes of genetic changes that preserve the progenitor state and/or interrupt normal development.

Unraveling genetic predisposition to familial or early onset gastric cancer using germline whole-exome sequencing

Recognition of individuals with a genetic predisposition to gastric cancer (GC) enables preventive measures. However, the underlying cause of genetic susceptibility to gastric cancer remains largely unexplained. We performed germline whole-exome sequencing on leukocyte DNA of 54 patients from 53 families with genetically unexplained diffuse-type and intestinal-type GC to identify novel GC-predisposing candidate genes. As young age at diagnosis and familial clustering are hallmarks of genetic tumor susceptibility, we selected patients that were diagnosed below the age of 35, patients from families with two cases of GC at or below age 60 and patients from families with three GC cases at or below age 70. All included individuals were tested negative for germline CDH1 mutations before or during the study. Variants that were possibly deleterious according to in silico predictions were filtered using several independent approaches that were based on gene function and gene mutation burden in controls. Despite a rigorous search, no obvious candidate GC predisposition genes were identified. This negative result stresses the importance of future research studies in large, homogeneous cohorts.

REST Final-Exon-Truncating Mutations Cause Hereditary Gingival Fibromatosis

Hereditary gingival fibromatosis (HGF) is the most common genetic form of gingival fibromatosis that develops as a slowly progressive, benign, localized or generalized enlargement of keratinized gingiva. HGF is a genetically heterogeneous disorder and can be transmitted either as an autosomal-dominant or autosomal-recessive trait or appear sporadically. To date, four loci (2p22.1, 2p23.3-p22.3, 5q13-q22, and 11p15) have been mapped to autosomes and one gene (SOS1) has been associated with the HGF trait observed to segregate in a dominant inheritance pattern. Here we report 11 individuals with HGF from three unrelated families. Whole-exome sequencing (WES) revealed three different truncating mutations including two frameshifts and one nonsense variant in RE1-silencing transcription factor (REST) in the probands from all families and further genetic and genomic analyses confirmed the WES-identified findings. REST is a transcriptional repressor that is expressed throughout the body; it has different roles in different cellular contexts, such as oncogenic and tumor-suppressor functions and hematopoietic and cardiac differentiation. Here we show the consequences of germline final-exon-truncating mutations in REST for organismal development and the association with the HGF phenotype.

Novel Combined Immune Deficiency and Radiation Sensitivity Blended Phenotype in an Adult with Biallelic Variations in ZAP70 and RNF168

With the advent of high-throughput genomic sequencing techniques, novel genetic etiologies are being uncovered for previously unexplained Mendelian phenotypes, and the underlying genetic architecture of disease is being unraveled. Although most of these “mendelizing” disease traits represent phenotypes caused by single-gene defects, a percentage of patients have blended phenotypes caused by pathogenic variants in multiple genes. We describe an adult patient with susceptibility to bacterial, herpesviral, and fungal infections. Immunologic defects included CD8⁺ T cell lymphopenia, decreased T cell proliferative responses to mitogens, hypogammaglobulinemia, and radiation sensitivity. Whole-exome sequencing revealed compound heterozygous variants in ZAP70. Biallelic mutations in ZAP70 are known to produce a spectrum of immune deficiency that includes the T cell abnormalities observed in this patient. Analyses for variants in genes associated with radiation sensitivity identified the presence of a homozygous RNF168 variant of unknown significance. RNF168 deficiency causes radiosensitivity, immunodeficiency, dysmorphic features, and learning difficulties syndrome and may account for the radiation sensitivity. Thus, the patient was found to have a novel blended phenotype associated with multilocus genomic variation: i.e., separate and distinct genetic defects. These findings further illustrate the clinical utility of applying genomic testing in patients with primary immunodeficiency diseases.

Improved full-length killer cell immunoglobulin-like receptor transcript discovery in Mauritian cynomolgus macaques

Killer cell immunoglobulin-like receptors (KIRs) modulate disease progression of pathogens including HIV, malaria, and hepatitis C. Cynomolgus and rhesus macaques are widely used as nonhuman primate models to study human pathogens, and so, considerable effort has been put into characterizing their KIR genetics. However, previous studies have relied on cDNA cloning and Sanger sequencing that lack the throughput of current sequencing platforms. In this study, we present a high throughput, full-length allele discovery method utilizing Pacific Biosciences circular consensus sequencing (CCS). We also describe a new approach to Macaque Exome Sequencing (MES) and the development of the Rhexome1.0, an adapted target capture reagent that includes macaque-specific capture probe sets. By using sequence reads generated by whole genome sequencing (WGS) and MES to inform primer design, we were able to increase the sensitivity of KIR allele discovery. We demonstrate this increased sensitivity by defining nine novel alleles within a cohort of Mauritian cynomolgus macaques (MCM), a geographically isolated population with restricted KIR genetics that was thought to be completely characterized. Finally, we describe an approach to genotyping KIRs directly from sequence reads generated using WGS/MES reads. The findings presented here expand our understanding of KIR genetics in MCM by associating new genes with all eight KIR haplotypes and demonstrating the existence of at least one KIR3DS gene associated with every haplotype.

Resolution of Disease Phenotypes Resulting from Multilocus Genomic Variation

BACKGROUND

Whole-exome sequencing can provide insight into the relationship between observed clinical phenotypes and underlying genotypes.

METHODS

We conducted a retrospective analysis of data from a series of 7374 consecutive unrelated patients who had been referred to a clinical diagnostic laboratory for whole-exome sequencing; our goal was to determine the frequency and clinical characteristics of patients for whom more than one molecular diagnosis was reported. The phenotypic similarity between molecularly diagnosed pairs of diseases was calculated with the use of terms from the Human Phenotype Ontology.

RESULTS

A molecular diagnosis was rendered for 2076 of 7374 patients (28.2%); among these patients, 101 (4.9%) had diagnoses that involved two or more disease loci. We also analyzed parental samples, when available, and found that de novo variants accounted for 67.8% (61 of 90) of pathogenic variants in autosomal dominant disease genes and 51.7% (15 of 29) of pathogenic variants in X-linked disease genes; both variants were de novo in 44.7% (17 of 38) of patients with two monoallelic variants. Causal copy-number variants were found in 12 patients (11.9%) with multiple diagnoses. Phenotypic similarity scores were significantly lower among patients in whom the phenotype resulted from two distinct mendelian disorders that affected different organ systems (50 patients) than among patients with disorders that had overlapping phenotypic features (30 patients) (median score, 0.21 vs. 0.36; P=1.77×10^-7).

CONCLUSIONS

In our study, we found multiple molecular diagnoses in 4.9% of cases in which whole-exome sequencing was informative. Our results show that structured clinical ontologies can be used to determine the degree of overlap between two mendelian diseases in the same patient; the diseases can be distinct or overlapping. Distinct disease phenotypes affect different organ systems, whereas overlapping disease phenotypes are more likely to be caused by two genes encoding proteins that interact within the same pathway. (Funded by the National Institutes of Health and the Ting Tsung and Wei Fong Chao Foundation.).

Primary immunodeficiency diseases: Genomic approaches delineate heterogeneous Mendelian disorders

BACKGROUND

Primary immunodeficiency diseases (PIDDs) are clinically and genetically heterogeneous disorders thus far associated with mutations in more than 300 genes. The clinical phenotypes derived from distinct genotypes can overlap. Genetic etiology can be a prognostic indicator of disease severity and can influence treatment decisions.

OBJECTIVE

We sought to investigate the ability of whole-exome screening methods to detect disease-causing variants in patients with PIDDs.

METHODS

Patients with PIDDs from 278 families from 22 countries were investigated by using whole-exome sequencing. Computational copy number variant (CNV) prediction pipelines and an exome-tiling chromosomal microarray were also applied to identify intragenic CNVs. Analytic approaches initially focused on 475 known or candidate PIDD genes but were nonexclusive and further tailored based on clinical data, family history, and immunophenotyping.

RESULTS

A likely molecular diagnosis was achieved in 110 (40%) unrelated probands. Clinical diagnosis was revised in about half (60/110) and management was directly altered in nearly a quarter (26/110) of families based on molecular findings. Twelve PIDD-causing CNVs were detected, including 7 smaller than 30 Kb that would not have been detected with conventional diagnostic CNV arrays.

CONCLUSION

This high-throughput genomic approach enabled detection of disease-related variants in unexpected genes; permitted detection of low-grade constitutional, somatic, and revertant mosaicism; and provided evidence of a mutational burden in mixed PIDD immunophenotypes.

Biallelic mutations in IRF8 impair human NK cell maturation and function

Human NK cell deficiencies are rare yet result in severe and often fatal disease, particularly as a result of viral susceptibility. NK cells develop from hematopoietic stem cells, and few monogenic errors that specifically interrupt NK cell development have been reported. Here we have described biallelic mutations in IRF8, which encodes an interferon regulatory factor, as a cause of familial NK cell deficiency that results in fatal and severe viral disease. Compound heterozygous or homozygous mutations in IRF8 in 3 unrelated families resulted in a paucity of mature CD56dim NK cells and an increase in the frequency of the immature CD56bright NK cells, and this impairment in terminal maturation was also observed in Irf8-/-, but not Irf8+/-, mice. We then determined that impaired maturation was NK cell intrinsic, and gene expression analysis of human NK cell developmental subsets showed that multiple genes were dysregulated by IRF8 mutation. The phenotype was accompanied by deficient NK cell function and was stable over time. Together, these data indicate that human NK cells require IRF8 for development and functional maturation and that dysregulation of this function results in severe human disease, thereby emphasizing a critical role for NK cells in human antiviral defense.

Clinical exome sequencing: results from 2819 samples reflecting 1000 families

We report our results of 1000 diagnostic WES cases based on 2819 sequenced samples from 54 countries with a wide phenotypic spectrum. Clinical information given by the requesting physicians was translated to HPO terms. WES processes were performed according to standardized settings. We identified the underlying pathogenic or likely pathogenic variants in 307 families (30.7%). In further 253 families (25.3%) a variant of unknown significance, possibly explaining the clinical symptoms of the index patient was identified. WES enabled timely diagnosing of genetic diseases, validation of causality of specific genetic disorders of PTPN23, KCTD3, SCN3A, PPOX, FRMPD4, and SCN1B, and setting dual diagnoses by detecting two causative variants in distinct genes in the same patient. We observed a better diagnostic yield in consanguineous families, in severe and in syndromic phenotypes. Our results suggest that WES has a better yield in patients that present with several symptoms, rather than an isolated abnormality. We also validate the clinical benefit of WES as an effective diagnostic tool, particularly in nonspecific or heterogeneous phenotypes. We recommend WES as a first-line diagnostic in all cases without a clear differential diagnosis, to facilitate personal medical care.