Next Article in Journal
Time Series RNA-seq in Pigeonpea Revealed the Core Genes in Metabolic Pathways under Aluminum Stress
Next Article in Special Issue
Altered Expression of GABAergic Markers in the Forebrain of Young and Adult Engrailed-2 Knockout Mice
Previous Article in Journal
Increasing the Genetic Diagnosis Yield in Inherited Retinal Dystrophies: Assigning Pathogenicity to Novel Non-canonical Splice Site Variants
Previous Article in Special Issue
Small Molecules Targeting H3K9 Methylation Prevent Silencing of Reactivated FMR1 Alleles in Fragile X Syndrome Patient Derived Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compound Phenotype Due to Recessive Variants in LARP7 and OTOG Genes Disclosed by an Integrated Approach of SNP-Array and Whole Exome Sequencing

1
Division of Medical Genetics, Fondazione IRCCS-Casa Sollievo della Sofferenza, 71013 San Giovanni Rotondo (Foggia), Italy
2
Bioinformatics Unit, Fondazione IRCCS-Casa Sollievo della Sofferenza, 71013 San Giovanni Rotondo (Foggia), Italy
*
Author to whom correspondence should be addressed.
Submission received: 24 February 2020 / Revised: 25 March 2020 / Accepted: 29 March 2020 / Published: 31 March 2020

Abstract

:
Neurodevelopmental disorders are a challenge in medical genetics due to genetic heterogeneity and complex genotype-phenotype correlations. For this reason, the resolution of single cases not belonging to well-defined syndromes often requires an integrated approach of multiple whole-genome technologies. Such an approach has also unexpectedly revealed a complex molecular basis in an increasing number of patients, for whom the original suspect of a pleiotropic syndrome has been resolved as the summation effect of multiple genes. We describe a 10-year-old boy, the third son of first-cousin parents, with global developmental delay, facial dysmorphism, and bilateral deafness. SNP-array analysis revealed regions of homozygosity (ROHs) in multiple chromosome regions. Whole-exome sequencing prioritized on gene-mapping into the ROHs showed homozygosity for the likely pathogenic c.1097_1098delAG p. (Arg366Thrfs*2) frameshift substitution in LARP7 and the likely pathogenic c.5743C>T p.(Arg1915*) nonsense variant in OTOG. Recessive variants in LARP7 cause Alazami syndrome, while variants in OTOG cause an extremely rare autosomal recessive form of neurosensorial deafness. Previously unreported features were acrocyanosis and palmoplantar hyperhidrosis. This case highlights the utility of encouraging technological updates in medical genetics laboratories involved in the study of neurodevelopmental disorders and integrating laboratory outputs with the competencies of next-generation clinicians.

1. Introduction

Neurodevelopmental disorders (NDDs) are the most frequent cause of disability in children and, currently, the main reason for referral in clinical genetic services. Nevertheless, they remain a major challenge due to high genetic-heterogeneity and weak genotype-phenotype correlations in most cases. The introduction of next-generation technologies has been demonstrated as effective in reducing the costs of the diagnostic trajectory of intellectual disability and in clinical decision-making for focal epilepsy [1,2]. In addition, whole-exome sequencing (WES)-based diagnostics is predicted to increase the number of individuals in whom the suspected pleiotropic syndrome is resolved as compound phenotypes, due to the simultaneous involvement of two or more disease-genes with separated biological functions (i.e., compound phenotypes). In this context, consanguinity facilitates the concurrence of multiple autosomal recessive disorders in the same individual [3].
Alazami syndrome (MIM 615071) is an ultra-rare autosomal recessive syndromal form of NDD, described for the first time in 2012 in an inbred Saudi Arabian family with a novel form of primordial dwarfism, due to a homozygous variant in LARP7 (OMIM 612026) gene [4]. Twenty-four patients have been reported to date [5], often presenting subtle but recurrent facial features, such as deep-set eyes, broad nose, wide mouth and teeth anomalies in addition to short stature, intellectual disabilities/global developmental delay and poor speech. LARP7 encodes for the La-related protein 7, which acts as a chaperone of the noncoding RNA 7SK [4]. Alazami syndrome is usually caused by frameshift variants in LARP7 [5]. Neurosensorial deafness is not a feature of Alazami syndrome.
OTOG gene (OMIM 604487) encodes for the Otogelin protein, a noncollagenous component of the acellular structures that cover the sensory epithelia of the inner ear [6]. Recessive variants in OTOG are a very rare cause of hereditary hearing loss, known as DFNB18B (MIM 614945) [7]. To date, only four causative variants in OTOG have been found in seven patients with DFNB18B: two nonsense, one frameshift, and one missense variant [8].
Here, we report a 10 year-old boy, born to consanguineous parents, with global developmental delay, absent speech, facial dysmorphisms, and bilateral congenital hearing loss. Multiple regions of homozygosity (ROH) identified by SNP-array-supported whole-exome sequencing analysis, allowing us to identify two homozygous deleterious variants in LARP7 and OTOG. This study reinforces the utility of combining high-resolution SNP-array technologies with WES in resolving highly complex phenotypes, especially in the presence of parental consanguinity.

2. Materials and Methods

2.1. Genomic DNA Extraction and Quantification

This family provided signed informed consent to molecular testing and to the full content of this publication. This study was conducted in accordance with the 1984 Declaration of Helsinki and its subsequent revisions. Molecular testing carried out in this report is based on the routine clinical care of our institution. Peripheral blood samples were taken from both the patient and his parents, and genomic DNA was isolated by using Bio Robot EZ1 (Quiagen, Solna, Sweden). The quality of DNA was tested on 1% electrophorese agarose gel, and the concentration was quantified by Nanodrop 2000 C spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).

2.2. SNP Array Analysis

SNP array analysis of the proband and his parents was carried out with the CytoScan HD Array (Thermo Fisher Scientific) as previously described [9]. Data analysis was performed using Chromosome Analysis Suite Software version 4.0 (Thermo Fisher Scientific) following a standardized pipeline. Briefly: (i) the raw data file (.CEL) was normalized using the default options; (ii) an unpaired analysis was performed using 270 HapMap samples as a baseline in order to obtain copy numbers value and regions of homozygosity (ROHs) from .CEL files. The amplified and/or deleted regions were detected using a standard Hidden Markov Model (HMM) method. Size threshold for analysis was kept as 5 Kb for copy number variations (CNVs), and 1 Mb for ROHs. In order to identify clinical or functionally relevant genomic variants, we compared all chromosomal alterations identified to those collected in our internal database of ~4,000 patients studied by SNP Arrays since 2010 and public databases, including the Database of Genomic Variants (DGV; available online at: http://projects.tcag.ca/variation/), DECIPHER (available online at: https://decipher.sanger.ac.uk/) and ClinVar (available online at: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/clinvar/). Base pair positions, information about genomic regions and genes affected by CNVs and/or ROHs, and known associated diseases will be derived from the University of California Santa Cruz (UCSC) Genome Browser (available online at: http://genome.ucsc.edu/cgi-bin/hgGateway), build GRCh37 (hg19). The clinical significance of each rearrangement detected was assessed following the American College of Medical Genetics (ACMG) guidelines [10].

2.3. Whole-Exome Sequencing

Proband DNA was analyzed by WES by using SureSelect Human Clinical Research Exome V6 (Agilent Technologies, Santa Clara, CA, USA) following manufacturer instructions. This is a combined shearing-free transposase-based library prep and target-enrichment solution, which enables comprehensive coverage of the entire exome. This system enables a specific mapping of reads to target deep coverage of protein-coding regions from RefSeq, GENCODE, CCDS, and UCSC Known Genes, with excellent overall exonic coverage and increased coverage of HGMD, OMIM, ClinVar, and ACMG targets. Sequencing was performed on a NextSeq 500 System (Illumina, San Diego, CA, USA) by using the High Output flow cells (300 cycles), with a minimum expected coverage depth of 70x. All variants obtained from WES were called by means of the HaplotypeCaller tool of GATK ver. 3.58 [11] and were annotated based on frequency, impact on the encoded protein, conservation, and expression using distinct tools, as appropriate (ANNOVAR, dbSNP, 1000 Genomes, EVS, ExAC, ESP, KAVIAR, and ClinVar) [12,13,14,15,16], and retrieving pre-computed pathogenicity predictions with dbNSFP v 3.0 (PolyPhen-2, SIFT, MutationAssessor, FATHMM, LRT and CADD) [17] and evolutionary conservation measures.
Next, variants prioritization was performed starting from homozygous variants, involving the genes located in the ROHs identified. Firstly, variants described as benign and likely benign were excluded. Then, remaining variants were classified based on their clinical relevance as pathogenic, likely pathogenic, or variant of uncertain significance according to following criteria: (i) nonsense/frameshift variant in genes previously described as disease-causing by haploinsufficiency or loss-of-function; (ii) missense variant located in a critical or functional domain; (iii) variant affecting canonical splicing sites (i.e., ±1 or ±2 positions); (iv) variant absent in allele frequency population databases; (v) variant reported in allele frequency population databases, but with a minor allele frequency (MAF) significantly lower than expected for the disease (<0.002 for autosomal recessive disease and <0.00001 for autosomal dominant disease); (vi) variant predicted and/or annotated as pathogenic/deleterious in ClinVar and/or LOVD; (vii) variants in GJB2 gene setting up as cut off a MAF <= 0.1 and variants for other genes associated with recessive hearing impairment and reported in Hereditary Hearing Loss Database (https://hereditaryhearingloss.org/) setting up as cut off a MAF <= 0.5.
The resulting putative pathogenic variants were confirmed by Sanger sequencing in both the proband and the parents’ DNA. PCR products were sequenced by using BigDye Terminator v1.1 Sequencing Kit (Applied Biosystems, Foster City, CA, USA) and ABI Prism 3100 Genetic Analyzer (Thermo Fisher Scientific). The clinical significance of the identified putative variants was interpreted according to the American College of Medical Genetics and Genomics (ACMG) [18]. Variant analysis was carried out considering the ethnicity of the patient.

2.4. Variant Designation

Nucleotide variants nomenclature follows the format indicated in the Human Genome Variation Society (HGVS, http://www.hgvs.org) recommendations and reported in the Leiden Open Variation Databases (LOVD) (https://databases.lovd.nl/shared/individuals/00276129).

3. Results

3.1. Clinical Description

The patient is a 10 year-old male, son of first cousins, unaffected parents of Caucasian origin (Southern Italy). His family history was unremarkable. The patient was born at term (39 weeks) with a low birth weight (2330 g). His birth length, head circumference, and Apgar score were unavailable. The neonatal period ran physiologically well, except for patent foramen ovale, but the following early psychomotor development was delayed. The patent foramen ovale resolved spontaneously in the following months. The patient held his head at 7 months, sat alone at about 18 months, and walked at 24 months. He never attained fully autonomous walking, language skills, or sphincter control. Due to the severe speech impairment, the patient underwent multiple audiological exams which diagnosed a profound bilateral neurosensorial deafness. Although the patient received auditory prosthesis, their language skills did not improve. At 1 year, a brain MRI scan was performed and did not reveal abnormalities except for a mild myelination delay. At examination, his height was 98 cm (<3rd centile), weight 19.8 kg (<3rd centile) and his head circumference was 49.5 cm (<3rd centile). His facial features included sparse eyebrows, prominent supraorbital ridges, short nose, high palate, and full lips. His language was absent. His gait was unstable with a broad base and mild spasticity of the lower limbs. The patient also displayed acrocyanosis and marked palmoplantar hyperhidrosis. The parents said that these two features were present since his infancy.

3.2. Molecular Findings

We did not identify any clinically significant CNV in the patient, while multiple regions of ROH were detected on chromosomes 2, 3, 4, 7, 8, 11, 16 and 18, respectively. The regions of homozygosity identified in the patient are depicted in Figure 1 and listed in Table 1.
This finding was made in accordance with parental consanguinity and offered us the opportunity to prioritize variant filtering from the WES data for homozygosity in disease-gene-mapping within these regions. Accordingly, WES allowed us to detect two previously unpublished, homozygous variants: a frameshift variant c.1097_1098del, p.(Arg366Thrfs*2) in LARP7 (NM_015454.2) mapping in the ROH region on chromosome 4, and a nonsense variant c.5743C>T, p.(Arg1915*) in OTOG (NM_001277269.1) mapping to the ROH region on chromosome 11. No further clinically relevant variant was identified at the homozygous state in the other ROH regions. Variants mapped outside these regions were prioritized according to the criteria listed in the Methods sections. No variant was identified as reaching the likely pathogenic or pathogenic rank according to the American College of Medical Genetics and Genomics and compatible with the observed phenotype.
The identified variants in LARP7 and OTOG were confirmed by Sanger sequencing (Figure 2) using specific primers (LARP7, exon 8, Forward Primer: GAATCCCTAGCTCCCCGATC; LARP7, exon 8, Reverse Primer: GTGCAGTTCTTGGCTACAGG. OTOG, exon 35, Forward Primer: CACTTAGCCCAGTACTGCCT; OTOG, exon 35, Reverse Primer: CTCAGGGCATAGGATGTGGG).
The LARP7 variant c.1097_1098del, p.(Arg366Thrfs*2) had an MAF of 0.00002 and 0.00003 in gnomAD and ExAC, respectively. The OTOG variant c.5743C>T, p.(Arg1915*) had an MAF of 0.00003 and 0.0001 in gnomAD and ExAC, respectively. The LARP7 variant is predicted to cause a frameshift leading to the deletion of the last 217 amino acids of the protein. The OTOG variant is predicted to truncate the protein and to cause loss of the last 1011 amino acids. Both variants were classified as likely pathogenic according to ACMG guidelines [18]. Complete bioinformatic details of these variants are reported in Table 2. No further variant classified as pathogenic or likely pathogenic, according to ACMG guidelines in other genes and previously associated with phenotypes compatible with the clinical features reported by the patients, were identified by the bioinformatics analysis. In particular, we did not find any additional deleterious variants in other genes responsible for hereditary hearing loss and annotated in the Hereditary Hearing Loss Database (https://hereditaryhearingloss.org/), or the other satellite symptoms, including palmoplantar hyperhidrosis and acrocyanosis.

4. Discussion

Here, we report a boy with severe NDD, absent speech with profound neurosensorial deafness, growth restriction, and facial dysmorphisms, resulting from double homozygosity in LARP7, causing Alazami syndrome and OTOG, associated with a very rare form of autosomal recessive neurosensorial deafness. Although the parental consanguinity predicted multiple ROHs at the SNP array analysis, it did not predict that there were multiple genes involved in the disease. Given our experience, we suggest, in cases like the one we describe, that the possibility of multiple genes involved in the etiology of the phenotype should be taken into account.
LARP7 gene encodes a chaperone protein required for both stability and function of the RNA. In particular, it forms a complex with the nuclear 7SK RNA to regulate RNA polymerase II transcription. To date, the clinical effect of deleterious variants in LARP7 is restricted to Alazami syndrome. Currently, 24 patients from 12 families are described as having Alazami syndrome and recessive variants in LARP7 [5]. Of the 14 identified variants, 12 were frameshift, one was nonsense, and one was predicted to cause abnormal splicing. Therefore, all deleterious variants causing Alazami syndrome are presumably null alleles [5]. Our findings are in line with the current literature, as the variant detected in our patient is a frameshift predicted to cause the deletion of the last 217 amino acids and, therefore, the loss of the RNA recognition motif 2 (RRM2). At a functional level, this event would likely lead to the formation of a truncated protein [p.(Arg366Thrfs*2)] or nonsense-mediated mRNA decay. This previously unpublished variant is located near to the variants reported in the original report, as well as those described by Hollink and colleagues [19]. This brings the number of identified deleterious variants in LARP7 to 15 (Figure 3).
The La protein RNA-binding motif is critical for the physiological function of the encoded protein as all identified deleterious variants cause the loss of this domain. On a clinical perspective, the overall phenotype observed in this patient is fully in agreement with the Alazami syndrome clinical spectrum, which mainly includes global developmental delay with absent speech, growth retardation (of prenatal onset), and peculiar facial appearance. Interestingly, the described patient also presented with bilateral profound neurosensorial deafness, acrocyanosis, and palmoplantar hyperhidrosis. While hearing impairment may be explained by the combined recessive variants in OTOG (see below), acrocyanosis and palmoplantar hyperhidrosis remain unexplained. The concurrence of these features might be causal, as both are unspecific and are common in the general population. Nevertheless, at the moment, we cannot exclude that acrocyanosis and palmoplantar hyperhidrosis might represent unusual phenotypic manifestations of LARP7, OTOG, or a combination of both. Further patients are needed in order to explore whether they are casually associated with Alazami syndrome or, rather, represent its unusual manifestations.
OTOG gene encodes for the Otogelin, an N-glycosylated protein that is expressed in the acellular membranes covering the six sensory epithelial patches of the inner ear, including the cochlea, the tectorial membrane over the organ of Corti, the vestibule, the otoconial membranes over the utricular and saccular maculae, and the cupulae over the cristae ampullares of the three semicircular canals. These membranes are collectively involved in the mechanotransduction process of external auditory stimuli. The movement of these membranes, which is induced by sound in the cochlea or acceleration in the vestibule, results in the deflection of the stereocilia bundle at the apex of the sensory hair cells, which in turn opens the mechanotransduction channels located at the tip of the stereocilia [6]. Several studies aimed at evaluating the expression level of otogelin in human tissues found that this protein is ubiquitous, but the transcript levels were the highest in the inner ear, followed by the kidney, the lung, the spleen, the thymus, and the liver [7]. To date, deleterious variants in OTOG are selectively associated with a rare non-syndromal form of autosomal recessive neurosensorial deafness. The present patient and previously published deleterious variants bring the number of described patients to eight and the number of reported variants to five (Figure 4).
The c.5743C>T variant is predicted to introduce a premature stop codon, which would, presumably, lead to the formation of a truncated protein [p.(Arg1915*)] or nonsense-mediated mRNA decay. This is in line with the predominance of null alleles in OTOG associated with hereditary deafness. In our patient, the hearing phenotype is apparently severe and of congenital onset, as suggested by previous papers speculating on the prenatal onset of the disease in OTOG-related autosomal recessive hereditary deafness [7,8].

5. Conclusions

Our case highlights the importance of combined genomic approaches in solving complex phenotypes and the role of parental consanguinity in prioritizing candidate gene-mapping in ROHs. We also identified novel deleterious variants in LARP7 and OTOG, which expand their molecular spectrum and confirm the associated phenotypes. This patient also presented acrocyanosis and palmoplantar hyperhidrosis, which might represent unusual phenotypic manifestations of LARP7, OTOG, or a combination of both.

Author Contributions

P.P., O.P. and M.C. (Marco Castori) designed the study; P.P. and O.P. wrote the draft; O.P. and M.C. (Massimo Carella) performed SNP array analysis; M.P.L, E.d.M., and L.B performed WES, data analysis, and Sanger sequencing; S.C. and T.M. carried out the bioinformatics analysis; L.B., M.P.L., O.P., and P.P. interpreted molecular data; M.C. (Marco Castori) provided the clinical evaluation of the patient and reviewed the final draft. M.C. (Massimo Carella) supervised the study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ricerca Corrente Program from the Italian Ministry of Health to M.Car.

Acknowledgments

The authors thank the family for their kind availability.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study or the collection, analyses, and interpretation of the data. They also had no role in the writing of the manuscript, or in the decision to publish the results.

References

  1. Vrijenhoek, T.; Middelburg, E.M.; Monroe, G.R.; van Gassen, K.L.I.; Geenen, J.W.; Hövels, A.M.; Knoers, N.V.; van Amstel, H.K.P.; Frederix, G.W.J. Whole-exome sequencing in intellectual disability; cost before and after a diagnosis. Eur. J. Hum. Genet. 2018, 26, 1566–1571. [Google Scholar] [CrossRef] [Green Version]
  2. Perucca, P.; Scheffer, I.E.; Harvey, A.S.; James, P.A.; Lunke, S.; Thorne, N.; Gaff, C.; Regan, B.M.; Damiano, J.A.; Hildebrand, M.S.; et al. Real-world utility of whole exome sequencing with targeted gene analysis for focal epilepsy. Epilepsy Res. 2017, 131, 1–8. [Google Scholar] [CrossRef]
  3. Lenfant, C.; Baz, P.; Degavre, A.; Philippi, A.; Senée, V.; Vandiedonck, C.; Derbois, C.; Nicolino, M.; Zalloua, P.; Julier, C. Juvenile-Onset Diabetes and Congenital Cataract: “Double-Gene” Mutations Mimicking a Syndromic Diabetes Presentation. Genes 2017, 8, 11. [Google Scholar] [CrossRef] [Green Version]
  4. Alazami, A.M.; Al-Owain, M.; Alzahrani, F.; Shuaib, T.; Al-Shamrani, H.; Al-Falki, Y.H.; Al-Qahtani, S.M.; Alsheddi, T.; Colak, D.; Alkuraya, F.S. Loss of function mutation in LARP7, chaperone of 7SK ncRNA, causes a syndrome of facial dysmorphism, intellectual disability, and primordial dwarfism. Hum. Mut. 2012, 33, 1429–1434. [Google Scholar] [CrossRef]
  5. Ivanovski, I.; Caraffi, S.G.; Magnani, E.; Rosato, S.; Pollazzon, M.; Matalonga, L.; Piana, S.; Nicoli, D.; Baldo, C.; Bernasconi, S.; et al. Alazami syndrome: The first case of papillary thyroid carcinoma. J. Hum. Genet. 2020, 65, 133–141. [Google Scholar] [CrossRef]
  6. Simmler, M.C.; Cohen-Salmon, M.; El-Amraoui, A.; Guillaud, L.; Benichou, J.C.; Petit, C.; Panthier, J.J. Targeted disruption of otog results in deafness and severe imbalance. Nat. Genet. 2000, 24, 139–143. [Google Scholar] [CrossRef]
  7. Schraders, M.; Ruiz-Palmero, L.; Kalay, E.; Oostrik, J.; del Castillo, F.J.; Sezgin, O.; Beynon, A.J.; Strom, T.M.; Pennings, R.J.; Zazo Seco, C.; et al. Mutations of the gene encoding otogelin are a cause of autosomal-recessive nonsyndromic moderate hearing impairment. Am. J. Hum. Genet. 2012, 91, 883–889. [Google Scholar] [CrossRef] [Green Version]
  8. Yu, S.; Choi, H.J.; Lee, J.S.; Lee, H.J.; Rim, J.H.; Choi, J.Y.; Gee, H.Y.; Jung, J. A novel early truncation mutation in OTOG causes prelingual mild hearing loss without vestibular dysfunction. Eur. J. Med. Genet. 2019, 62, 81–84. [Google Scholar] [CrossRef]
  9. Palumbo, O.; Fichera, M.; Palumbo, P.; Rizzo, R.; Mazzolla, E.; Cocuzza, D.M.; Carella, M.; Mattina, T. TBR1 is the candidate gene for intellectual disability in patients with a 2q24.2 interstitial deletion. Am. J. Med. Genet. 2014, 163, 828–833. [Google Scholar] [CrossRef]
  10. Kearney, H.M.; Thorland, E.C.; Brown, K.K.; Quintero-Rivera, F.; South, S.T.; Working Group of the American College of Medical Genetics Laboratory Quality Assurance Committee. American College of Medical Genetics standards and guidelines for interpretation and reporting of postnatal constitutional copy number variants. Genet. Med. 2011, 13, 680–685. [Google Scholar] [CrossRef] [Green Version]
  11. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from next-generation sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef]
  13. Sherry, S.T.; Ward, M.H.; Kholodov, M.; Baker, J.; Phan, L.; Smigielski, E.M.; Sirotkin, K. dbSNP: The NCBI database of genetic variation. Nucleic Acids Res. 2001, 29, 308–311. [Google Scholar] [CrossRef] [Green Version]
  14. Lek, M.; Karczewski, K.J.; Minikel, E.V.; Samocha, K.E.; Banks, E.; Fennell, T.; O’Donnell-Luria, A.H.; Ware, J.S.; Hill, A.J.; Cummings, B.B.; et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016, 536, 285–291. [Google Scholar] [CrossRef] [Green Version]
  15. Glusman, G.; Caballero, J.; Mauldin, D.E.; Hood, L.; Roach, J.C. KAVIAR: An accessible system for testing SNV novelty. Bioinformatics 2011, 27, 3216–3217. [Google Scholar] [CrossRef] [Green Version]
  16. Landrum, M.J.; Lee, J.M.; Riley, G.R.; Jang, W.; Rubinstein, W.S.; Church, D.M.; Maglott, D.R. ClinVar: Public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014, 42, D980–D985. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, X.; Wu, C.; Li, C.; Boerwinkle, E. dbNSFP v3.0: A One-Stop Database of Functional Predictions and Annotations for Human Non-synonymous and Splice Site SNVs. Hum. Mutat. 2016, 37, 235–241. [Google Scholar] [CrossRef] [Green Version]
  18. Richards, S.; Aziz, N.; Bale, S.; Bick, D.; Das, S.; Gastier-Foster, J.; Grody, W.W.; Hegde, M.; Lyon, E.; Spector, E.; et al. ACMG Laboratory Quality Assurance Committee. Standards and guidelines for the interpretation of sequence variants: A joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 2015, 17, 405–424. [Google Scholar] [CrossRef]
  19. Hollink, I.H.; Alfadhel, M.; Al-Wakeel, A.S.; Ababneh, F.; Pfundt, R.; de Man, S.A.; Jamra, R.A.; Rolfs, A.; Bertoli-Avella, A.M.; van de Laar, I.M. Broadening the phenotypic spectrum of pathogenic LARP7 variants: Two cases with intellectual disability, variable growth retardation and distinct facial features. J. Hum. Genet. 2016, 61, 229–233. [Google Scholar] [CrossRef]
  20. Uchikawa, E.; Natchiar, K.S.; Han, X.; Proux, F.; Roblin, P.; Zhang, E.; Durand, A.; Klaholz, B.P.; Dock-Bregeon, A.C. Structural insight into the mechanism of stabilization of the 7SK small nuclear RNA by LARP7. Nucleic Acids Res. 2015, 43, 3373–3388. [Google Scholar] [CrossRef]
  21. Eichhorn, C.D.; Chug, R.; Feigon, J. hLARP7 C-terminal domain contains an xRRM that binds the 3’ hairpin of 7SK RNA. Nucleic Acids Res. 2016, 44, 9977–9989. [Google Scholar] [CrossRef] [Green Version]
  22. Najmabadi, H.; Hu, H.; Garshasbi, M.; Zemojtel, T.; Abedini, S.S.; Chen, W.; Hosseini, M.; Behjati, F.; Haas, S.; Jamali, P.; et al. Deep sequencing reveals 50 novel genes for recessive cognitive disorders. Nature 2011, 44, 57–63. [Google Scholar] [CrossRef]
  23. Ling, T.T.; Sorrentino, S. Compound heterozygous variants in the LARP7 gene as a cause of Alazami syndrome in a Caucasian female with significant failure to thrive, short stature, and developmental disability. Am. J. Med. Genet. 2016, 170, 217–219. [Google Scholar] [CrossRef]
  24. Dateki, S.; Kitajima, T.; Kihara, T.; Watanabe, S.; Yoshiura, K.I.; Moriuchi, H. Novel compound heterozygous variants in the LARP7 gene in a patient with Alazami syndrome. Hum Genome Var. 2018, 5, 18014. [Google Scholar] [CrossRef]
  25. Imbert-Bouteille, M.; Mau Them, F.T.; Thevenon, J.; Guignard, T.; Gatinois, V.; Riviere, J.B.; Boland, A.; Meyer, V.; Deleuze, J.F.; Sanchez, E.; et al. LARP7 variants and further delineation of the Alazami syndrome phenotypic spectrum among primordial dwarfisms: 2 sisters. Eur. J. Med. Genet. 2019, 62, 161–166. [Google Scholar] [CrossRef]
  26. Wojcik, M.H.; Linnea, K.; Stoler, J.M.; Rappaport, L. Updating the neurodevelopmental profile of Alazami syndrome: Illustrating the role of developmental assessment in rare genetic disorders. Am. J. Med. Genet. 2019, 179, 1565–1569. [Google Scholar] [CrossRef]
  27. Goodyear, R.J.; Richardson, G.P. Structure, Function, and Development of the Tectorial Membrane: An Extracellular Matrix Essential for Hearing. Curr. Top. Dev. Biol. 2018, 130, 217–244. [Google Scholar]
Figure 1. Karyoview of the patient. Regions of homozygosity (ROHs) identified in the proband genome are indicated by purple bars.
Figure 1. Karyoview of the patient. Regions of homozygosity (ROHs) identified in the proband genome are indicated by purple bars.
Genes 11 00379 g001
Figure 2. Family tree and molecular findings in the family. (a) Double lines in pedigree indicate consanguinity. Filled and unfilled circles/squares represent affected and unaffected individuals respectively. (b) Sanger sequence of a PCR product amplified with primers targeting exon 8 of LARP7 of individual II.1 and his unaffected parents. (c) Sanger sequence of a PCR product amplified with primers targeting exon 35 of OTOG of individual II.1 and his unaffected parents. Vertical red arrows indicate variant position.
Figure 2. Family tree and molecular findings in the family. (a) Double lines in pedigree indicate consanguinity. Filled and unfilled circles/squares represent affected and unaffected individuals respectively. (b) Sanger sequence of a PCR product amplified with primers targeting exon 8 of LARP7 of individual II.1 and his unaffected parents. (c) Sanger sequence of a PCR product amplified with primers targeting exon 35 of OTOG of individual II.1 and his unaffected parents. Vertical red arrows indicate variant position.
Genes 11 00379 g002
Figure 3. Schematic representation of the protein domains of La-related protein 7 [20,21], encoded by LARP7, with positions of all known deleterious variants published to date and including the present one [4,5,19,22,23,24,25,26]. The variant identified here is in red. (LaM, Lupus antigen Motif; RRM1, RNA Recognition Motif 1; RRM2, RNA Recognition Motif 2).
Figure 3. Schematic representation of the protein domains of La-related protein 7 [20,21], encoded by LARP7, with positions of all known deleterious variants published to date and including the present one [4,5,19,22,23,24,25,26]. The variant identified here is in red. (LaM, Lupus antigen Motif; RRM1, RNA Recognition Motif 1; RRM2, RNA Recognition Motif 2).
Genes 11 00379 g003
Figure 4. Schematic representation of protein domains of Otogelin [27], encoded by OTOG, with positions of all known deleterious variants published so far, and including the present one [7,8]. The variant identified here is in red. (vWFD, von Willebrand factor type D domain; C8, cystein-rich domain; TIL, trypsin inhibitor-like cysteine-rich domain; vWFC, von Willebrand factor type C domain; ZP, zona pellucida domain; C-terminal cysteine knot-like domain).
Figure 4. Schematic representation of protein domains of Otogelin [27], encoded by OTOG, with positions of all known deleterious variants published so far, and including the present one [7,8]. The variant identified here is in red. (vWFD, von Willebrand factor type D domain; C8, cystein-rich domain; TIL, trypsin inhibitor-like cysteine-rich domain; vWFC, von Willebrand factor type C domain; ZP, zona pellucida domain; C-terminal cysteine knot-like domain).
Genes 11 00379 g004
Table 1. List of regions of homozygosity (ROHs) identified in the proband genome according to the International System for Human Cytogenetic Nomenclature (ISCN 2016). Size and number of genes included in each ROH are indicated. All genomic coordinates were based on the GRCh37/hg19 build of the Human Genome.
Table 1. List of regions of homozygosity (ROHs) identified in the proband genome according to the International System for Human Cytogenetic Nomenclature (ISCN 2016). Size and number of genes included in each ROH are indicated. All genomic coordinates were based on the GRCh37/hg19 build of the Human Genome.
Regions of Homozigosity (ISCN 2016)Size (Mb)Number of Genes
arr[GRCh37] 2p25.2p25.1(5573938_10537316)x2 hmz4.940
arr[GRCh37] 2p21p11.2(47453937_89129064)x2 hmz41.6317
arr[GRCh37] 2q11.1q22.1(95341387_137166151)x2 hmz41.8340
arr[GRCh37] 3p13p11.1(71199337_90485635)x2 hmz19.249
arr[GRCh37] 3q11.1q13.33(93536053_120421935)x2 hmz26.8173
arr[GRCh37] 4p15.2p11(21674784_49089181)x2 hmz27.4117
arr[GRCh37] 4q11q31.1(52686799_141024195)x2 hmz88.3468
arr[GRCh37] 7p14.2p12.2(36814859_50418506)x2 hmz13.698
arr[GRCh37] 7p12.1p11.1(51736013_58019983)x2 hmz6.238
arr[GRCh37] 7q11.21q11.22(62461703_68908285)x2 hmz6.450
arr[GRCh37] 8p23.3p23.1(168483_7011075)x2 hmz6.846
arr[GRCh37] 11p15.5p13(1751363_33420180)x2 hmz31.6367
arr[GRCh37] 16p13.3p13.2(6643315_8047081)x2 hmz1.41
arr[GRCh37] 18q12.2q21.1(34083335_43959703)x2 hmz9.833
Mb, Megabases.
Table 2. Characteristics of the variants identified in the LARP7 and OTOG genes.
Table 2. Characteristics of the variants identified in the LARP7 and OTOG genes.
ChromosomeStartEndReference AlleleAlternative AlleleGenotypeGeneNucleotide ChangeAmino Acid ChangeGene ImpactdbSNP IDgnomAD_exome Allele CountTOPMED Allele CountExAC_ALL
Allele Count
4113568939113568941AAGAhomo_altLARP7
NM_015454.2
c.1097_1098delp.(Arg366Thrfs*2)frameshift substitutionrs5664642494/237240
MAF 0.00002
6/125568
MAF 0.00005
3/113350
0.00003
111763255417632554CThomo_altOTOG
NM_001277269.1
c.5743C>Tp.(Arg1915*)stopgainrs7612870445/146398
MAF 0.00003
2/125568
MAF 0.00002
1/15676
0.0001

Share and Cite

MDPI and ACS Style

Palumbo, P.; Palumbo, O.; Leone, M.P.; di Muro, E.; Castellana, S.; Bisceglia, L.; Mazza, T.; Carella, M.; Castori, M. Compound Phenotype Due to Recessive Variants in LARP7 and OTOG Genes Disclosed by an Integrated Approach of SNP-Array and Whole Exome Sequencing. Genes 2020, 11, 379. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11040379

AMA Style

Palumbo P, Palumbo O, Leone MP, di Muro E, Castellana S, Bisceglia L, Mazza T, Carella M, Castori M. Compound Phenotype Due to Recessive Variants in LARP7 and OTOG Genes Disclosed by an Integrated Approach of SNP-Array and Whole Exome Sequencing. Genes. 2020; 11(4):379. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11040379

Chicago/Turabian Style

Palumbo, Pietro, Orazio Palumbo, Maria Pia Leone, Ester di Muro, Stefano Castellana, Luigi Bisceglia, Tommaso Mazza, Massimo Carella, and Marco Castori. 2020. "Compound Phenotype Due to Recessive Variants in LARP7 and OTOG Genes Disclosed by an Integrated Approach of SNP-Array and Whole Exome Sequencing" Genes 11, no. 4: 379. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11040379

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop