Genome sequencing identifies a large non-coding region deletion of SNX10 causing autosomal recessive osteopetrosis.

Autosomal recessive osteopetrosis (ARO) is a rare genetic disorder caused by impaired osteoclast activity. In this study, we describe a 4-year-old boy with increased bone density due to osteopetrosis, autosomal recessive 8. Using genome sequencing, we identified a large deletion in the 5'-untranslated region (UTR) of SNX10 (sorting nexin 10), where the regulatory region of this gene is located. This large deletion resulted in the absence of the SNX10 transcript and led to abnormal osteoclast activity. SNX10 is one of the nine genes known to cause ARO, shown to interact with V-ATPase (vacuolar type H( + )-ATPase), as it plays an important role in bone resorption. Our study highlights the importance of regulatory regions in the 5'-UTR of SNX10 for its expression while also demonstrating the importance of genome sequencing for detecting large deletion of the regulatory region of SNX10.

A homozygous hypomorphic BNIP1 variant causes an increase in autophagosomes and reduced autophagic flux and results in a spondylo-epiphyseal dysplasia.

BNIP1 (BCL2 interacting protein 1) is a soluble N-ethylmaleimide-sensitive factor-attachment protein receptor involved in ER membrane fusion. We identified the homozygous BNIP1 intronic variant c.84+3A>T in the apparently unrelated patients 1 and 2 with disproportionate short stature. Radiographs showed abnormalities affecting both the axial and appendicular skeleton and spondylo-epiphyseal dysplasia. We detected ~80% aberrantly spliced BNIP1 pre-mRNAs, reduced BNIP1 mRNA level to ~80%, and BNIP1 protein level reduction by ~50% in patient 1 compared to control fibroblasts. The BNIP1 ortholog in Drosophila, Sec20, regulates autophagy and lysosomal degradation. We assessed lysosome positioning and identified a decrease in lysosomes in the perinuclear region and an increase in the cell periphery in patient 1 cells. Immunofluorescence microscopy and immunoblotting demonstrated an increase in LC3B-positive structures and LC3B-II levels, respectively, in patient 1 fibroblasts under steady-state condition. Treatment of serum-starved fibroblasts with or without bafilomycin A1 identified significantly decreased autophagic flux in patient 1 cells. Our data suggest a block at the terminal stage of autolysosome formation and/or clearance in patient fibroblasts. BNIP1 together with RAB33B and VPS16, disease genes for Smith-McCort dysplasia 2 and a multisystem disorder with short stature, respectively, highlight the importance of autophagy in skeletal development.

Biallelic loss-of-function variants in EXOC6B are associated with impaired primary ciliogenesis and cause spondylo-epi-metaphyseal dysplasia with joint laxity type 3.

Spondylo-epi-metaphyseal dysplasias with joint laxity, type 3 (SEMDJL3) is a genetic skeletal disorder characterized by multiple joint dislocations, caused by biallelic pathogenic variants in the EXOC6B gene. Only four individuals from two families have been reported to have this condition to date. The molecular pathogenesis related to primary ciliogenesis has not been enumerated in subjects with SEMDJL3. In this study, we report two additional affected individuals from unrelated families with biallelic pathogenic variants, c.2122+15447_2197-59588del and c.401T>G in EXOC6B identified by exome sequencing. One of the affected individuals had an intellectual disability and central nervous system anomalies, including hydrocephalus, hypoplastic mesencephalon, and thin corpus callosum. Using the fibroblast cell lines, we demonstrate the primary evidence for the abrogation of exocytosis in an individual with SEMDLJ3 leading to impaired primary ciliogenesis. Osteogenesis differentiation and pathways related to the extracellular matrix were also found to be reduced. Additionally, we provide a review of the clinical and molecular profile of all the mutation-proven patients reported hitherto, thereby further characterizing SEMDJL3. SEMDJL3 with biallelic pathogenic variants in EXOC6B might represent yet another ciliopathy with central nervous system involvement and joint dislocations.

Biallelic deep intronic variant c.5457+81T>A in TRIP11 causes loss of function and results in achondrogenesis 1A.

Biallelic loss of function variants in TRIP11 encoding for the Golgi microtubule-associated protein 210 (GMAP-210) causes the lethal chondrodysplasia achondrogenesis type 1A (ACG1A). Loss of TRIP11 activity has been shown to impair Golgi structure, vesicular transport, and results in loss of IFT20 anchorage to the Golgi that is vital for ciliary trafficking and ciliogenesis. Here, we report four fetuses, two each from two families, who were ascertained antenatally with ACG1A. Affected fetuses in both families are homozygous for the deep intronic TRIP11 variant, c.5457+81T>A, which was found in a shared region of homozygosity. This variant was found to cause aberrant transcript splicing and the retention of 77 base pairs of intron 18. The TRIP11 messenger RNA and protein levels were drastically reduced in fibroblast cells derived from one of the affected fetuses. Using immunofluorescence we also detected highly compacted Golgi apparatus in affected fibroblasts. Further, we observed a significant reduction in the frequency of ciliated cells and in the length of primary cilia in subject-derived cell lines, not reported so far in patient cells with TRIP11 null or hypomorphic variants. Our findings illustrate how pathogenic variants in intronic regions of TRIP11 can impact transcript splicing, expression, and activity, resulting in ACG1A.

A data set of variants derived from 1455 clinical and research exomes is efficient in variant prioritization for early-onset monogenic disorders in Indians.

Given the genomic uniqueness, a local data set is most desired for Indians, who are underrepresented in existing public databases. We hypothesize patients with rare monogenic disorders and their family members can provide a reliable source of common variants in the population. Exome sequencing (ES) data from families with rare Mendelian disorders was aggregated from five centers in India. The dataset was refined by excluding related individuals and removing the disease-causing variants (refined cohort). The efficiency of these data sets was assessed in a new set of 50 exomes against gnomAD and GenomeAsia. Our original cohort comprised 1455 individuals from 1203 families. The refined cohort had 836 unrelated individuals that retained 1,251,064 variants with 181,125 population-specific and 489,618 common variants. The allele frequencies from our cohort helped to define 97,609 rare variants in gnomAD and 44,520 rare variants in GenomeAsia as common variants in our population. Our variant dataset provided an additional 1.7% and 0.1% efficiency for prioritizing heterozygous and homozygous variants respectively for rare monogenic disorders. We observed additional 19 genes/human knockouts. We list carrier frequency for 142 recessive disorders. This is a large and useful resource of exonic variants for Indians. Despite limitations, datasets from patients are efficient tools for variant prioritization in a resource-limited setting.

Clinical and genetic spectrum of 104 Indian families with central nervous system white matter abnormalities.

Genetic disorders with predominant central nervous system white matter abnormalities (CNS WMAs), also called leukodystrophies, are heterogeneous entities. We ascertained 117 individuals with CNS WMAs from 104 unrelated families. Targeted genetic testing was carried out in 16 families and 13 of them received a diagnosis. Chromosomal microarray (CMA) was performed for three families and one received a diagnosis. Mendeliome sequencing was used for testing 11 families and all received a diagnosis. Whole exome sequencing (WES) was performed in 80 families and was diagnostic in 52 (65%). Singleton WES was diagnostic for 50/75 (66.67%) families. Overall, genetic diagnoses were obtained in 77 families (74.03%). Twenty-two of 47 distinct disorders observed in this cohort have not been reported in Indian individuals previously. Notably, disorders of nuclear mitochondrial pathology were most frequent (9 disorders in 20 families). Thirty-seven of 75 (49.33%) disease-causing variants are novel. To sum up, the present cohort describes the phenotypic and genotypic spectrum of genetic disorders with CNS WMAs in our population. It demonstrates WES, especially singleton WES, as an efficient tool in the diagnosis of these heterogeneous entities. It also highlights possible founder events and recurrent disease-causing variants in our population and their implications on the testing strategy.

An emerging ribosomopathy affecting the skeleton due to biallelic variations in NEPRO.

Cartilage hair hypoplasia (CHH), anauxetic dysplasia 1, and anauxetic dysplasia 2 are rare metaphyseal dysplasias caused by biallelic pathogenic variants in RMRP and POP1, which encode the components of RNAse-MRP endoribonuclease complex (RMRP) in ribosomal biogenesis pathway. Nucleolus and neural progenitor protein (NEPRO), encoded by NEPRO (C3orf17), is known to interact with multiple protein subunits of RMRP. We ascertained a 6-year-old girl with skeletal dysplasia and some features of CHH. RMRP and POP1 did not harbor any causative variant in the proband. Parents-child trio exomes revealed a candidate biallelic variant, c.435G>C, p.(Leu145Phe) in NEPRO. Two families with four affected individuals with skeletal dysplasia and a homozygous missense variant, c.280C>T, p.(Arg94Cys) in NEPRO, were identified from literature and their published phenotype was compared in detail to the phenotype of the child we described. All the five affected individuals have severe short stature, brachydactyly, skin laxity, joint hypermobility, and joint dislocations. They also have short metacarpals, broad middle phalanges, and metaphyseal irregularities. Protein modeling and stability prediction showed that the mutant protein has decreased stability. Both the reported variants are in the same domain of the protein. Our report delineates the clinical and radiological characteristics of an emerging ribosomopathy caused by biallelic variants in NEPRO.

Multilocus disease-causing genomic variations for Mendelian disorders: role of systematic phenotyping and implications on genetic counselling.

Multilocus disease-causing genomic variations (MGVs) and multiple genetic diagnoses (MGDs) are increasingly being recognised in individuals and families with Mendelian disorders. This can be mainly attributed to the widespread use of genomic tests for the evaluation of these disorders. We conducted a retrospective study of families evaluated over the last 6 years at our centre to identify families with MGVs and MGDs. MGVs were observed in fourteen families. We observed five different consequences: (i) individuals with MGVs presenting as blended phenotypes (ii) individuals with MGVs presenting with distinct phenotypes (iii) individuals with MGVs with age-dependent penetrance (iv) individuals with MGVs with one phenotype obscured by another more predominant phenotype (v) two distinct phenotypes in different individuals in families with MGVs. Consanguinity was present in eight (8/14, 57.1%) of them. Thirteen families had two Mendelian disorders and one had three Mendelian disorders. The risk of recurrence of one or more conditions in these families ranged from 25% to 75%. Our findings underline the importance of the role of a clinical geneticist in systematic phenotyping, challenges in genetic counselling and risk estimation in families with MGVs and MGDs, especially in highly inbred populations.

Clinically relevant variants in a large cohort of Indian patients with Marfan syndrome and related disorders identified by next-generation sequencing.

Marfan syndrome and related disorders are a group of heritable connective tissue disorders and share many clinical features that involve cardiovascular, skeletal, craniofacial, ocular, and cutaneous abnormalities. The majority of affected individuals have aortopathies associated with early mortality and morbidity. Implementation of targeted gene panel next-generation sequencing in these individuals is a powerful tool to obtain a genetic diagnosis. Here, we report on clinical and genetic spectrum of 53 families from India with a total of 83 patients who had a clinical diagnosis suggestive of Marfan syndrome or related disorders. We obtained a molecular diagnosis in 45/53 (85%) index patients, in which 36/53 (68%) had rare variants in FBN1 (Marfan syndrome; 63 patients in total), seven (13.3%) in TGFBR1/TGFBR2 (Loeys-Dietz syndrome; nine patients in total) and two patients (3.7%) in SKI (Shprintzen-Goldberg syndrome). 21 of 41 rare variants (51.2%) were novel. We did not detect a disease-associated variant in 8 (15%) index patients, and none of them met the Ghent Marfan diagnostic criteria. We found the homozygous FBN1 variant p.(Arg954His) in a boy with typical features of Marfan syndrome. Our study is the first reporting on the spectrum of variants in FBN1, TGFBR1, TGFBR2, and SKI in Indian individuals.