Molecular and computational characterization of ABCB11 and ABCG5 variants in Tunisian patients with neonatal/infantile low-GGT intrahepatic cholestasis: Genetic diagnosis and genotype-phenotype correlation assessment.

Many inherited conditions cause hepatocellular cholestasis in infancy, including progressive familial intrahepatic cholestasis (PFIC), a heterogeneous group of diseases with highly overlapping symptoms. In our study, six unrelated Tunisian infants with PFIC suspicion were the subject of a panel-target sequencing followed by an exhaustive bioinformatic and modeling investigations. Results revealed five disease-causative variants including known ones: (the p.Asp482Gly and p.Tyr354 * in the ABCB11 gene and the p.Arg446 * in the ABCC2 gene), a novel p.Ala98Cys variant in the ATP-binding cassette subfamily G member 5 (ABCG5) gene and a first homozygous description of the p.Gln312His in the ABCB11 gene. The p.Gln312His disrupts the interaction pattern of the bile salt export pump as well as the flexibility of the second intracellular loop domain harboring this residue. As for the p.Ala98Cys, it modulates both the interactions within the first nucleotide-binding domain of the bile transporter and its accessibility. Two additional potentially modifier variants in cholestasis-associated genes were retained based on their pathogenicity (p.Gly758Val in the ABCC2 gene) and functionality (p.Asp19His in the ABCG8 gene). Molecular findings allowed a PFIC2 diagnosis in five patients and an unexpected diagnosis of sisterolemia in one case. The absence of genotype/phenotype correlation suggests the implication of environmental and epigenetic factors as well as modifier variants involved directly or indirectly in the bile composition, which could explain the cholestasis phenotypic variability.

Characterization of a missense variant in COG5 in a Tunisian patient with COG5-CDG syndrome and insights into the effect of non-synonymous variants on COG5 protein.

The clinical diagnosis of patients with multisystem involvement including a pronounced neurologic damage is challenging. High-throughput sequencing methods remains crucial to provide an accurate diagnosis. In this study, we reported a Tunisian patient manifesting hypotonia and global developmental delay with visual and skin abnormalities. Exome sequencing was conducted followed by segregation analysis and, subsequently additional investigations. In silico analysis of non-synonymous variants (nsSNPs) described in COG5 in conserved positions was made. Results revealed a homozygous missense variant c.298 C > T (p.Leu100Phe) in the COG5 inherited from both parents. This variant altered both protein solubility and stability, in addition to a putative disruption of the COG5-COG7 interaction. This disruption has been confirmed using patient-derived cells in vitro in a COG5 co-immuno-precipitation, where interaction with binding partner COG7 was abrogated. Hence, we established the COG5-CDG diagnosis. Clinically, the patient shared common features with the already described cases with the report of the ichtyosis as a new manifestation. Conversely, the CADD scoring revealed 19 putatively pathogenic nsSNPs (Minor Allele Frequency MAF < 0.001, CADD > 30), 11 of which had a significant impact on the solubility and/or stability of COG5. These properties seem to be disrupted by six of the seven missense COG5-CDG variants. In conclusion, our study expands the genetic and phenotypic spectrum of COG5-CDG disease and highlight the utility of the next generation sequencing as a powerful tool in accurate diagnosis. Our results shed light on a likely molecular mechanism underlying the pathogenic effect of missense COG5 variants, which is the alteration of COG5 stability and solubility.

Molecular and in silico investigation of a novel ECHS1 gene mutation in a consanguine family with short-chain enoyl-CoA hydratase deficiency and Mt-DNA depletion: effect on trimer assembly and catalytic activity.

Short-chain enoyl-CoA hydratase deficiency (ECHS1D) is a rare congenital metabolic disorder that follows an autosomal recessive inheritance pattern. It is caused by mutations in the ECHS1 gene, which encodes a mitochondrial enzyme involved in the second step of mitochondrial ß-oxidation of fatty acids. The main characteristics of the disease are severe developmental delay, regression, seizures, neurodegeneration, high blood lactate, and a brain MRI pattern consistent with Leigh syndrome. Here, we report three patients belonging to a consanguineous family who presented with mitochondrial encephalomyopathy. Whole-exome sequencing revealed a new homozygous mutation c.619G > A (p.Gly207Ser) at the last nucleotide position in exon 5 of the ECHS1 gene. Experimental analysis showed that normal ECHS1 pre-mRNA splicing occurred in all patients compared to controls. Furthermore, three-dimensional models of wild-type and mutant echs1 proteins revealed changes in catalytic site interactions, conformational changes, and intramolecular interactions, potentially disrupting echs1 protein trimerization and affecting its function. Additionally, the quantification of mtDNA copy number variation in blood leukocytes showed severe mtDNA depletion in all probands.

miRNA implication in the pathogenesis and the outcome of Tunisian endemic pemphigus foliaceus.

Pemphigus foliaceus (PF) is a bullous autoimmune skin disease diagnosed through sera and skin analyses. PF severity is associated with maintained anti-Dsg1 sera levels and its prognosis is unpredictable. MicroRNA (miRNA), dynamic regulators of immune function, have been identified as potential biomarkers for some autoimmune diseases. This study aimed to assess the miRNA expression of miR-17-5p, miR-21-5p, miR-146a-5p, miR-155-5p and miR-338-3p using quantitative real-time PCR in peripheral blood mononuclear cells (PBMC) and lesional skin samples from untreated and treated PF patients (both remittent and chronic) over 3 months. Overall, miRNA expression was significantly higher in PBMC than in biopsy samples. Blood miR-21 expression was increased in untreated patients compared to controls and had a diagnostic value with an AUC of 0.78. After 6 weeks, it decreased significantly, similar to anti-Dsg1 antibodies and the PDAI score. In addition, a positive correlation was observed between cutaneous miR-21 expression and the disease activity score. Conversely, cutaneous expressions of miR-17, miR-146a and miR-155 were significantly higher in treated chronic patients compared to remittent ones. The cutaneous level of miR-155 positively correlated with pemphigus activity, making it a potential predictive marker for patients' clinical stratification with an AUC of 0.86.These findings suggest that blood miR-21 and cutaneous miR-155 can be used as supplemental markers for PF diagnosis and activity, respectively in addition to classical parameters.

Potential dysfunctional effects of synonymous variants: Insights from an exhaustive in silico analysis of the ABCB4 gene.

The multiple drug resistance 3 (MDR3) protein is a canalicular phospholipid translocator involved in the bile secretion and encoded by the ABCB4 gene. Its deficiency is related to a large spectrum of liver diseases. Taking into account the increased evidence about the involvement of synonymous variants in inherited diseases, this study aims to explore the putative effects of silent genetic variants on the ABCB4 expression. We performed an exhaustive computational approach using ESE finder, RegRNA 2.0, MFOLD, SNPfold, and %MinMax software added to the measurement of the Relative Synonymous Codon Usage. This analysis included 216 synonymous variants distributed throughout the ABCB4 gene. Results have shown that 11 synonymous coding SNPs decrease the ESE activity, while 8 of them change the codon frequency. Besides, the c.24C>T variation, located 21 nucleotides downstream the start A (Adenine) U (Uracil) G (Glutamine) AUG causes an increase in the local stability. Moreover, the computational analysis of the 3'UTR region showed that six of the eight variants located in this region affected the Wild Type (WT) pattern of the miRNA targets sites and/or their proper display. The 26 sSNPs retained as putatively functional possessed a very low allele frequency, supporting their pathogenicity. In conclusion, the obtained results suggest that some synonymous SNPs in the ABCB4 gene, considered up to now as neutral, may be involved in the MDR3 deficiency.

Clinical, Molecular, and Computational Analysis in two cases with mitochondrial encephalomyopathy associated with SUCLG1 mutation in a consanguineous family.

Deficiency of the mitochondrial enzyme succinyl COA ligase (SUCL) is associated with encephalomyopathic mtDNA depletion syndrome and methylmalonic aciduria. This disorder is caused by mutations in both SUCL subunits genes: SUCLG1 (α subnit) and SUCLA2 (ß subnit). We report here, two Tunisian patients belonging to a consanguineous family with mitochondrial encephalomyopathy, hearing loss, lactic acidosis, hypotonia, psychomotor retardation and methylmalonic aciduria. Mutational analysis of SUCLG1 gene showed, for the first time, the presence of c.41T > C in the exon 1 at homozygous state. In-silico analysis revealed that this mutation substitutes a conserved methionine residue to a threonine at position 14 (p.M14T) located at the SUCLG1 protein mitochondrial targeting sequence. Moreover, these analysis predicted that this mutation alter stability structure and mitochondrial translocation of the protein. In Addition, a decrease in mtDNA copy number was revealed by real time PCR in the peripheral blood leukocytes in the two patients compared with controls.

Exploring the role of vitamin D-VDR pathway in pemphigus foliaceous: a novel perspective on disease pathogenesis.

Several auto-immune diseases have been linked to vitamin D deficiency as a contributing environmental factor. Its pleiotropic effects on the immune system, especially its essential role in maintaining immune tolerance, make the vitamin D pathway of great interest. In this study, we focused on Pemphigus foliaceous (PF) in Tunisian population. we aimed to quantify the Serum 25[OH]D levels using chemiluminescence assay and to analyze the differential expression of the VDR, CYP27B1 and CYP24A1 genes in the circulating blood cells and lesional skin tissue of PF patients using Q-PCR. A genetic explanation was then sought to explore any direct relationship between tag polymorphisms and the inherited features of PF. Results confirmed a vitamin D hypovitaminosis in Tunisian PF patients. Interestingly, a differential gene expression correlated to the disease stratification was noted. Indeed, at the systemic level, an upregulation of VDR and CYP27B1 genes was observed in healthy controls compared to PF patients. Notably, in lesional skin tissue, the clinical and serological remission phase was correlated with high transcriptional levels of the VDR gene and conversely a drop in expression of the CYP24A1 gene. Genetic analysis indicated the involvement of the most appealing polymorphisms, rs2228570 and poly (A) microsatellite, in PF etiopathogenesis. Indeed, CAC13 haplotype was associated with a higher risk of PF development. Our findings suggest that alterations in the vitamin D-VDR pathway may influence PF physiopathology, making this pathway a potential target for pharmacological modulation, especially for cortico-resistant PF patients.

Characterization of a novel ABCC2 mutation in infantile Dubin Johnson syndrome.

BACKGROUND AND AIMS: The Dubin Johnson Syndrome (DJS) occurs mostly in young adults but an early-onset of the disease has been reported in less common forms (Neonatal DJS and Infantile DJS). In this case, the clinical findings are of limit for the DJS diagnosis. Hence, the genetic testing remains the method of choice to provide an accurate diagnosis. In our study, we aimed to perform a genetic analysis for two siblings presented with an intrahepatic cholestasis before the age of 1 year to provide a molecular explanation for the developed phenotype. PATIENTS & METHODS: A Tunisian family, having two siblings, manifesting signs of a hepatopathy, was enrolled in our study. A molecular analysis was performed, using a panel-based next generation sequencing, supplying results that were the subject of computational analysis. Then, a clinical follow-up was carried out to assess the evolution of the disease. RESULTS: The genetic analysis revealed the presence of a novel missense c.4179G > T, (p.M1393I) mutation in ABCC2 gene associated with a substitution c.2789G > A (R930Q) in ATP8B1 gene. Predictive results consolidated the pathogenic effect of both variants. These results confirmed the DJS diagnosis in the studied patients. The clinical course of both patients fit well with the benign nature of DJS. CONCLUSION: We described here a novel ABCC2 mutation associated with a putative ATP8B1 modifier variant. This finding constituted the first report of a complex genotype in DJS. Hence, genetic analysis by a panel-based next generation sequencing permits an accurate diagnosis and the identification of putative variants that could influence the developed phenotype.

A homozygous ABCB4 mutation causing an LPAC syndrome evolves into cholangiocarcinoma.

Low phospholipid-associated cholelithiasis (LPAC) is characterized by the association of ABCB4 mutations and low biliary phospholipid concentration with symptomatic and recurring cholelithiasis. In the present study, we reported a case of a 63-year-old woman, who presented a biliary pain beginning at the age of 30, recurrent after cholecystectomy, along with "comet-tail shadows" revealed by ultrasonography thus, fulfilling the diagnosis of LPAC. This disease evolved into a cholangiocarcinoma. To understand the molecular basis of this phenotype, we performed the ABCB4 gene sequencing, followed by in silico analysis and Q-RT-PCR assay. The results displayed a homozygous missense sequence variation (c.140G > A, p.Arg47Gln), predicted as pathogenic according to MutPred. Accordingly, this gave rise to a decreased hepatic ABCB4 mRNA level and structural alterations of the mutated protein. Eventually, we reported, here, the first description of an ABCB4 missense mutation (p.Arg47Gln) at homozygous state in a Tunisian LPAC syndrome. An elucidation of its functional consequences was performed. Besides, this case suggests that the delayed diagnosis of LPAC syndrome and the lack of UDCA treatment may contribute in the development of complications, such as cholangiocarcinoma.

Phenotypic variability in Tunisian PFIC3 patients harboring a complex genotype with a differential clinical outcome of UDCA treatment.

INTRODUCTION: Progressive familial intrahepatic cholestasis type 3 (PFIC3) is a chronic autosomal recessive disorder characterized by a wide spectrum of clinical severity generally related to the degree of pathogenicity of the causal sequence variation in ABCB4 gene. PATIENTS AND METHODS: The present study reports the molecular investigation by Next Generation Sequencing (NGS) of five related patients with PFIC3 disease followed by bioinformatic analysis. A biochemical follow-up is also performed to assess the response of the ursodeoxycholic acid treatment. RESULTS: The molecular results revealed complex genotype in homozygous state in all patients including a pathogenic c.1436C > T (P479L) variation in the ABCB4 gene and two well-known polymorphisms, the V444A in ABCB11 gene and the D19H in the ABCG8 gene. Although the presence of the same genetic background, all patients present the disease at different ages and clinical signs with a variable degree of clinical severity at diagnosis. Additionally, a differential outcome to the treatment has been pointed out. CONCLUSION: Our results provide evidence regarding the putative intervention of modifier factors in the phenotypic variability reported for the first time in the PFIC3 disease and highlight the importance of an early administration of the UDCA as a good solution to ovoid the disease progression.

Clinical, Molecular, and Computational Analysis Showed a Novel Homozygous Mutation Among the Substrate-Binding Site of ARSA Protein in Consanguineous Family with Late-Infantile MLD.

Metachromatic leukodystrophy (MLD) is a neurodegenerative disorder characterized by progressive demyelination resulting from impaired degradation and thus the accumulation of cerebroside-3-sulfate (sulfatide). It is caused by the deficiency of arylsulfatase A (ARSA) enzyme which is encoded by the ARSA gene. The present study reports the clinical, molecular, and bioinformatic investigation of three patients belonging to a consanguineous family with late-infantile MLD disorder. The results revealed a novel homozygous missense mutation c.699C>A (p.His231Gln) in exon 4 of ARSA gene in the three patients inherited from their heterozygous parents. Interestingly, this novel mutation is the second mutation identified in the substrate-binding site of ARSA protein and it was classified as damaging and deleterious by several bioinformatics tools. The c.699C>A (p.His231Gln) leads to changes in the pre-mRNA secondary structure and in the ARSA protein 3D structure with a significant root mean square deviation value which could probably affect its stability and function.

Novel mutations in the CDKL5 gene in complex genotypes associated with West syndrome with variable phenotype: First description of somatic mosaic state.

INTRODUCTION: West syndrome is a rare epileptic encephalopathy of early infancy, characterized by epileptic spasms, hypsarrhythmia, and psychomotor retardation beginning in the first year of life. METHODS: The present study reports the clinical, molecular and bioinformatic investigation in the three studied West patients. RESULTS: The results revealed a complex genotype with more than one mutation in each patient including the known mutations c.1910C>G (P2, P3); c.2372A>C in P3 and c.2395C>G in P1 and novel variants including c.616G>A, shared by the three patients P1, P2 and P3; c.1403G>C shared by P2 and P3 and c.2288A>G in patient P1. CONCLUSIONS: All the mutations were at somatic mosaic state and were de novo in the patients except ones (c.2372A>C). To our knowledge; the somatic mosaic state is described for the first time in patients with West syndrome. Five identified mutations were located in the C-terminal domain of the protein, while the novel mutation (c.616G>A) was in the catalytic domain. Bioinformatic tools predicted that this latter is the most pathogenic substitution affecting 3D protein structure and the secondary mRNA structure. Complex genotype composed of different combinations of mutations in each patient seems to be related to the phenotype variability.

Comparison of in silico prediction and experimental assessment of ABCB4 variants identified in patients with biliary diseases.

Genetic variations of the phosphatidylcholine transporter, ABCB4 cause several biliary diseases. The large number of reported variations makes it difficult to foresee a comprehensive study of each variation. To appreciate the reliability of in silico prediction programs, 1) we confronted them with the assessment in cell models of two ABCB4 variations (E528D and P1161S) identified in patients with low phospholipid-associated cholelithiasis (LPAC); 2) we extended the confrontation to 19 variations that we had previously characterized in cellulo. Four programs (Provean, Polyphen-2, PhD-SNP and MutPred) were used to predict the degree of pathogenicity. The E528D and P1161S variants were studied in transfected HEK293 and HepG2 cells by immunofluorescence, immunoblotting and measurement of phosphatidylcholine secretion. All prediction tools qualified the P1161S variation as deleterious, but provided conflicting results for E528D. In cell models, both mutants were expressed and localized as the wild type but their activity was significantly reduced, by 48% (P1161S) and 33% (E528D). These functional defects best correlated with MutPred predictions. MutPred program also proved the most accurate to predict the pathogenicity of the 19 ABCB4 variants that we previously characterized in cell models, and the most sensitive to predict the pathogenicity of 65 additional mutations of the Human Gene Mutation Database. These results confirm the pathogenicity of E528D and P1161S variations and suggest that even a moderate decrease (by less than 50%) of phosphatidylcholine secretion can cause LPAC syndrome. They highlight the reliability of in silico prediction tools, most notably MutPred, as a first approach to predict the pathogenicity of ABCB4 variants.

In silico investigation of the impact of synonymous variants in ABCB4 gene on mRNA stability/structure, splicing accuracy and codon usage: Potential contribution to PFIC3 disease.

Progressive Familial Intrahepatic Cholestasis type 3 (PFIC3) is an autosomal-recessive liver disease due to mutations in the ABCB4 gene encoding for the MDR3 protein. In the present study, we performed molecular and bioinformatic analyses in PFIC3 patients in order to understand the molecular basis of the disease. The three studied patients with PFIC3 were screened by PCR amplification followed by direct sequencing of the 27 coding exons of ABCB4. In silico analysis was performed by bioinformatic programs. We revealed three synonymous polymorphisms c.175C>T, c.504C>T, c.711A>T respectively in exon 4, 6, 8 and an intronic c.3487-16T>C variation in intron 26. The computational study of these polymorphic variants using Human Splicing Finder, ex-skip, Mfold and kineFold tools showed the putative impact on the composition of the cis-acting regulatory elements of splicing as well as on the mRNA structure and stability. Moreover, the protein level was affected by codon usage changes estimated by the calculation of ΔRSCU and ΔLog Ratio of codon frequencies interfering as consequence with the accurate folding of the MDR3 protein. As the first initiative of the mutational study of ABCB4 genes in Tunisia, our results are suggestive of a potential downstream molecular effect for the described polymorphisms on the expression pattern of the ABCB4 underlining the importance of synonymous variants.