Abstract

Background. Neonatal intrahepatic cholestasis caused by citrin deficiency (NICCD) is a Mendelian disorder arising from biallelic SLC25A13 mutations, and SLC25A13 genetic analysis was indispensable for its definite diagnosis. However, conventional SLC25A13 analysis could not detect all mutations, especially obscure large insertions/deletions. This paper aimed to explore the obscure SLC25A13 mutation in an NICCD infant. Methods. Genomic DNA was extracted to screen for 4 high-frequency SLC25A13 mutations, and then all 18 exons and their flanking sequences were analyzed by Sanger sequencing. Subsequently, cDNA cloning, SNP analyses, and semiquantitative PCR were performed to identify the obscure mutation. Results. A maternally inherited mutation IVS16ins3kb was screened out, and then cDNA cloning unveiled paternally inherited alternative splicing variants (ASVs) featuring exon 5 skipping. Ultimately, a large deletion c.329-1687_c.468+3865del5692bp, which has never been described in any other references, was identified via intensive study on the genomic DNA around exon 5 of SLC25A13 gene. Conclusions. An NICCD patient was definitely diagnosed as a compound heterozygote of IVS16ins3kb and c.329-1687_c.468+3865del5692bp. The large deletion enriched the SLC25A13 mutation spectrum, and its identification supported the concept that cDNA cloning analysis, along with other molecular tools such as semiquantitative PCR, could provide valuable clues, facilitating the identification of obscure SLC25A13 deletions.

1. Introduction

Human citrin deficiency is an autosomal recessive disease due to dysfunction of citrin, the liver-type calcium-stimulated aspartate-glutamate carrier isoform 2 (AGC2) encoded by the SLC25A13 gene [1, 2]. Up to now, three age-dependent clinical phenotypes had been reported for this disorder, that is, adult-onset citrullinemia type 2 (CTLN2) in adolescents and adults, neonatal intrahepatic cholestasis caused by citrin deficiency (NICCD) in neonates or infants, and failure to thrive and dyslipidemia caused by citrin deficiency (FTTDCD) in older children [38]. Although most reported patients were Asian [914], citrin deficiency has been recognized as a worldwide panethnic disorder nowadays [1519].

Due to the lack of well-recognized clinical/biochemical diagnostic criteria for NICCD, SLC25A13 genetic analysis has been regarded as reliable tool for the definite diagnosis of such patients. However, routine SLC25A13 genetic analyses, such as polymerase chain reaction (PCR), long and accurate-PCR (LA-PCR), PCR-restriction fragment length polymorphism (RFLP), and Sanger sequencing, could not detect all SLC25A13 mutations, especially large deletions or insertions of an obscure nature. It has been estimated that approximately 15% of compound heterozygotes or homozygotes carrying SLC25A13 mutations in both alleles could not be definitely diagnosed just by the above conventional approaches [20]. In such cases, other molecular tools, although usually labor-intensive and cost-expensive, were needed to identify the obscure mutations.

In this study, a large SLC25A13 deletion was identified via sophisticated molecular analyses using peripheral blood lymphocytes (PBLs) in an NICCD patient, who responded well to a lactose-free and medium-chain triglycerides- (MCTs-) enriched formula. We herein reported the clinical and molecular findings.

2. Subjects and Methods

2.1. Subjects and Ethics

The research subjects in this study were a male patient suspected to have NICCD and his parents as well. With the informed consent from the parents and the ethical approval by the medical ethical committee of our hospital, we performed intensive clinical and genetic study on this family.

2.2. Conventional DNA Analysis

The DNA was extracted from the peripheral venous blood following the genomic DNA extraction kit (Omega, USA) instructions. Four high-frequency SLC25A13 mutations, c.851_854del, c.1638_1660dup, IVS6+5G>A, and IVS16ins3kb, were screened by PCR, LA-PCR, and PCR-RFLP procedures, and then Sanger sequencing of all the 18 exons and their flanking sequences was undertaken, using direct sequencing of DNA fragments amplified by genomic DNA-PCR to identify novel mutation in the gene SLC25A13 [21].

2.3. Reverse Transcriptional PCR (RT-PCR) and Nested PCR

PBLs were separated using lymphocyte separation medium (LSM, MP) from 2-3 milliliters of ethylene diamine tetraacetic acid- (EDTA-) anticoagulant peripheral venous blood [20]. PBLs were cracked with ribonucleic acid (RNA) trizol reagent (The Life Technologies) and total RNA was extracted using a variety of organic solvents, as described previously [12, 2225]. The cDNAs were synthesized from total RNA by Moloney murine leukemia virus (MMLV) reverse transcriptase (TaKaRa). Nested PCR was used to amplify SLC25A13 open reading frame (ORF), and the primers and PCR temperature profile were described in detail as in previous references [21, 22]. The two primer pairs in Nested PCR were RAS2 and RACEA1 in the first PCR and RAS3 and Ex18R in the second one, with expected PCR products of 3107 bp and 2191 bp in size, respectively. Both of them contained the citrin-coding sequence (CDS).

2.4. Molecular Cloning and Alternative Splicing Variants (ASVs) Analyses

The Nested PCR products of paternal origin were purified by using a gel extraction kit (Omega) and then connected with the sequence of the pMD 18-T vector (TaKaRa) and transformed into DH5α Escherichia coli competent cells, as in our previous publications [12, 22, 24, 25]. Following that, the transformed cells were cultured in the shaking table for 4 h and then coated on plates after centrifugation on 1000 rpm for 5 min. After cultivation in 37°C for 16 hours, the positive clones which had the targeted bands after PCR amplification were selected and sequenced. The ASV sequences were analyzed by using the ContigExpress and DNAMAN software, and their designation was according to the nomenclature guidelines [26, 27].

2.5. Further Location and Identification of the Obscure Mutation

According to the ASV structures, to further locate the DNA span around exon 5 that might contain the obscure mutation, semiquantitative PCR was performed in a total volume of 50 μL, containing 5 μL of 10x Buffer (Mg2+ plus) (TaKaRa), 4 μL of dNTP (10 mM), 37.75 μL of sterilized distilled water, 0.25 μL of Taq (TaKaRa), 2 μL of the forward and reverse primers together, and 1 μL of DNA template. Then, all tubes were placed into a thermal cycler and the parameters were set in 94°C for 5 minutes followed by 28 to 30 cycles of 94°C for 30 seconds, 60°C for 40 seconds, 72°C for 40 seconds, and a final extension step at 72°C for 10 minutes. The two primer pairs were IVS5S5 and IVS5A6 in Set 1 and IVS 1NF and IVS 2NB in Set 2, whose sequences and locations were displayed in Figure 1, respectively.

Based on the findings above, a PCR approach using the primers IVS4S3 and J5.6KbDelR1 (Set 3) was carried out to identify the obscure mutation. The reaction mixture was a total volume of 50 μL, containing 5 μL of 10x LA Buffer (Mg2+ plus) (TaKaRa), 8 μL of dNTP (10 mM), 30.3 μL of sterilized distilled water, 0.2 μL of Taq, 0.5 μL of LA Taq (TaKaRa), 4 μL of the forward and reverse primers together, and 2 μL of DNA template. The temperature profile was 95°C for 5 minutes followed by 40 cycles of 94°C for 40 seconds, 55°C for 40 seconds, 68°C for 10 minutes, and a final extension step at 68°C for 10 minutes.

2.6. Statistical Analysis

The proportions of the ASVs with exon 5 skipping in all ASVs from the patient and 8 healthy volunteers were statistically assessed using SPSS (Version 13.0) chi-square test, and a value of less than 0.05 was considered statistically significant.

3. Results

3.1. Clinical Findings

A male infant at the age of 2 months was referred to the local hospital with the chief complaint of jaundiced skin over 10 days. Physical examination at referral revealed a mildly pale face and hemorrhagic spots scattered on the skin of whole body. The lungs were clear on auscultation. No abnormal sound or murmur was heard on heart auscultation. An enlarged liver, 5.0 cm below the right costal margin, was palpated. No pathological reflexes could be found on nervous system examination. Slightly visible pitting edema could be found in both lower extremities.

Biochemical test demonstrated elevated aspartate aminotransferase (AST), gamma-glutamyl transpeptidase (GGT), total bilirubin (Tbil), direct bilirubin (Dbil), indirect bilirubin (Ibil), and total bile acids (TBA), which indicated intrahepatic cholestasis. The albumin level was decreased while that of alpha-feto protein (AFP) elevated markedly (Table 1). Ultrasonography revealed hepatomegaly, ascites, and abdominal bloating, along with a patent oval foramen of the heart. The urinary gas chromatography mass spectrometry (GC-MS) analysis discovered large quantity of hexanedioic acid and 4-hydroxyphenyl lactate (4HPL), and on tandem mass spectrometry (MS-MS) analysis, elevated levels of citrulline and arginine were detected.

NICCD was thus suspected based on the above clinical and biochemical findings, and according to the pediatrician’s advice, breast feeding was stopped and a lactose-free and MCT-enriched therapeutic formula was introduced immediately. In the subsequent follow-up, the clinical manifestations were alleviated gradually. At his age of 3.7 months, the infant was referred to our hospital to confirm the NICCD diagnosis by SLC25A13 gene analyses.

3.2. Screening and Sequencing Results

Screening of the 4 high-frequency mutations proved that the patient harbored a maternally inherited mutation IVS16ins3kb, but another mutation was not detected. Even after Sanger sequencing of all the 18 exons and their flanking sequences, the SLC25A13 mutation of paternal origin remained obscure. Therefore, cDNA cloning analysis was undertaken to facilitate the identification of the obscure mutation.

3.3. Findings of cDNA Cloning Analysis

As the maternal mutation IVS16ins3kb led to production of an aberrant mRNA molecule with a new exon 17, losing the normal exons 17 and 18 [21], only the ASVs transcribed from the paternal SLC25A13 allele were amplified by the Nested PCR approach using the two reverse primers both within the normal exon 18 of SLC25A13 gene. In other words, the 26 ASVs in Table 2 all resulted from the paternally inherited SLC25A13 allele.

Similar to our previous findings [22], these ASVs demonstrated remarkable structural heterogeneity. However, all these ASVs featured exon 5 skipping (r.329_468del), as shown in Table 2. When compared with the ASVs in the 8 healthy volunteers which had been reported previously [24], exon 5 skipping was a unique structural feature of the ASVs in this patient ( and as in Table 3), strongly suggesting the existence of a large insertion or deletion around exon 5 at the DNA level.

3.4. Positioning and Sequencing Analysis of the Obscure Mutation

Based on the findings above, a diversity of primers was designed and PCR amplification conducted to analyze a variety of SNPs within the DNA fragment around exon 5 in all 3 family members. Although no clues suggestive of the obscure mutation could be detected by these analyses, semiquantitative PCR had positive findings—when amplified using the primer Set 1 near exon 5, the DNA samples of the patient (P1) and the father (F1) had less PCR products in comparison to that of the mother (M1), but this is not the case when using the primer Set 2 covering the entire exon 2 (Figure 1), indicating that the patient and the father might have a large deletion or insertion involving the positions of the primer Set 1. According to this result, LA-PCR amplification using the primer Set 3 (Figure 2) yielded an unexpected band of 632 bp in size in the patient and father, besides the expected 6324 bp product as in the mother, and subsequent Sanger sequencing of the unexpected product uncovered a large deletion c.329-1687_c.468+3865del5692bp.

4. Discussion

Although a maternally originated SLC25A13 mutation IVS16ins3kb had been discovered as a high-frequency mutation in this study, conventional genetic analyses, such as PCR/LA-PCR, PCR-RFLP, and Sanger sequencing, could not unveil the paternally inherited mutation in the infant with typical clinical and biochemical features of intrahepatic cholestasis, making the definite diagnosis of NICCD a challenge. It was by using cDNA cloning along with SNP analysis and semiquantitative PCR in all family members that c.329-1687_c.468+3865del5692bp, a large deletion, was finally identified, which has never been described in any other references. Although time- and cost-consuming, the identification of this novel deletion confirmed the feasibility of SLC25A13 cDNA cloning analysis using PBLs as a molecular tool facilitating the identification of large deletions and provided reliable evidence for the definite diagnosis of NICCD in the patient along with the mutation IVS16ins3kb of maternal origin. This was important not only for the proband himself but also for future genetic counseling and antenatal diagnosis in the family.

This 5692 bp deletion resulted in exon 5 skipping at mRNA level, as described in Table 2, theoretically yielding a truncated citrin protein of 127 amino acids—the exon 5 skipping gave rise to a frameshift at codon 110, continued to translate 17 amino acids, and then encountered a terminator codon at the position of 127, eventually producing an aberrant citrin molecule p.E110fs127X. Interestingly, this truncated citrin molecule was the same as the one arising from the large insertion IVS4ins6kb which had been reported very recently by our group [24]. The aberrant mRNA with exon 5 skipping could be a target for nonsense-mediated mRNA decay (NMD). Moreover, according to the molecular structural feature of mitochondrial AGCs (aspartate/glutamate carriers) [28], such a truncated citrin molecule predictively lost the entire mitochondrial carrier domain and had no AGC2 function to export aspartate and import glutamate together with a proton into the mitochondrial matrix, thus causing a series of metabolic disturbances and finally leading to the formation of the laboratory and clinical presentations in the patient.

The lactose-free and MCT-enriched formulas have been reported to be clinically and biochemically effective for NICCD patients in increasing clinical cases [12, 25, 2931]. The clinical process in this study further supported this concept. Overall, the increased cytosolic NADH/NAD+ ratio in hepatocytes has been well-recognized as a key biochemical alteration for citrin deficiency [6]. Breast milk or common formula contained high carbohydrate in the form of lactose, a well-known disaccharide that was easily digested by lactase on the intestinal mucosa into glucose and galactose, both of which were then rapidly absorbed into the blood flow. In hepatocytes, galactose was conversed via Leloir pathway into glucose [32], the metabolism of which produced NADH, and thus increased the cytosolic NADH/NAD+ ratio in the liver and caused liver damage in patients with citrin deficiency [6, 31]. Moreover, secondary galactosemia due to citrin deficiency might be involved in the NICCD pathogenesis—the increased cytosolic NADH/NAD+ ratio in citrin-deficient hepatocytes might inhibit the activity of uridine diphosphate- (UDP-) galactose-4-epimerase [33], leading to accumulation of a large quantity of galactitol and galactonate [34], and galactitol has been suggested as one of the substrates causing jaundice, hepatosplenomegaly, hepatocellular insufficiency, and cataracts [30, 35, 36].

Energy shortage in the liver caused by an impairment of glycolysis due to an increased NADH/NAD+ ratio has been proposed as an important pathophysiology in citrin deficiency [31, 37]. Since the galactose kinase reaction in Leloir pathway was energy-consuming [32], the galactose metabolism in hepatocytes inevitably exacerbated such a pathophysiology. On the other hand, the absorption of MCTs was not bile acid-dependent, which might reduce the burden of the liver to synthesize and excrete bile salt into the gut [30]. Of particular note, MCTs can be better absorbed and transported via the portal vein into the liver and were mainly metabolized to acetyl-CoA along with FADH2 and NADH, which could supply more such substrates to hepatic cells as energy sources [31]. Taking all these factors together, it was not surprising for the NICCD infant in this study to respond well to the lactose-free and MCT-enriched therapeutic formula.

In summary, by sophisticated molecular analysis using PBLs, this study definitely diagnosed an NICCD patient who was a compound heterozygote of the IVS16ins3kb mutation and a novel deletion c.329-1687_c.468+3865del5692bp in the SLC25A13 gene. The large deletion constituted a novel component in the SLC25A13 mutation spectrum, and its identification lent further support to the concept that SLC25A13 cDNA cloning analysis using PBLs, along with other molecular tools such as semiquantitative PCR, could provide valuable clues, facilitating the identification of unknown large deletions.

Abbreviations

NICCD:Neonatal intrahepatic cholestasis caused by citrin deficiency
SLC25A13:Solute carrier family 25, member 13
DNA:Deoxyribonucleic acid
cDNA:Complementary deoxyribonucleic acid
SNP:Single nucleotide polymorphism
PCR:Polymerase chain reaction
ASVs:Alternative splicing variants
CTLN2:Adult-onset citrullinemia type 2
FTTDCD:Failure to thrive and dyslipidemia caused by citrin deficiency
LA-PCR:Long and accurate-polymerase chain reaction
PCR-RFLP:Polymerase chain reaction-restriction fragment length polymorphism
PBLs:Peripheral blood lymphocytes
MCTs:Medium-chain triglycerides
RT-PCR:Reverse transcriptional polymerase chain reaction
LSM:Lymphocyte separation medium
EDTA:Ethylene diamine tetraacetic acid
RNA:Ribonucleic acid
MMLV:Moloney murine leukemia virus
ORF:Open reading frame
CDS:Coding sequence
AST:Aspartate aminotransferase
GGT:Gamma-glutamyl transpeptidase
Tbil:Total bilirubin
Dbil:Direct bilirubin
Ibil:Indirect bilirubin
TBA:Total bile acids
AFP:Alpha-feto protein
GC-MS:Gas chromatography mass spectrometry
4HPL:4-Hydroxyphenyl lactate
MS-MS:Tandem mass spectrometry
NMD:Nonsense-mediated mRNA decay
AGC:Aspartate/glutamate carrier
NADH:Reduced form of nicotinamide adenine dinucleotide
:Nicotinamide adenine dinucleotide
UDP-galactose-4-epimerase:Uridine diphosphate-galactose-4-epimerase.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors would like to express their deep appreciation to all research subjects for their kind cooperation. The research was financially supported by the Cultivation Foundation for Scientific Research (nos. 2014208 and A2014101) approved by the First Affiliated Hospital, Jinan University, the Medical Scientific Research Foundation of Guangdong Province (A2014385), the Fundamental Research Funds for the Central Universities (no. 21615470), and the projects supported by the National Natural Science Foundation of China (nos. 81270957 and 81570793).