# Targeted-Capture Next-Generation Sequencing in Diagnosis Approach of Pediatric Cholestasis

Marion Almes, Anne Spraul, Mathias Ruiz, Muriel Girard, Bertrand

Roquelaure, Nolwenn Laborde, Frederic Gottrand, Anne Turquet, Thierry

Lamireau, Alain Dabadie, et al.

► To cite this version:

Marion Almes, Anne Spraul, Mathias Ruiz, Muriel Girard, Bertrand Roquelaure, et al. Targeted-Capture Next-Generation Sequencing in Diagnosis Approach of Pediatric Cholestasis. Diagnostics, 2022, 12 (5), 10.3390/diagnostics12051169. [https://doi.org/10.3390/diagnostics12051169](https://doi.org/10.3390/diagnostics12051169)

HAL Id: hal-03697994

https://hal.science/hal-03697994

Submitted on 17 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Targeted-Capture Next-Generation Sequencing in Diagnosis Approach of Pediatric Cholestasis

Marion Almes<sup>1,2,3,*†</sup>, Anne Spraul<sup>2,3,4,†</sup>, Mathias Ruiz<sup>5</sup>, Muriel Girard<sup>6</sup>, Bertrand Roquelaure<sup>7</sup>, Nolwenn Laborde<sup>8</sup>, Frédéric Gottrand<sup>9</sup>, Anne Turquet<sup>10</sup>, Thierry Lamireau<sup>11</sup>, Alain Dabadie<sup>12</sup>, Marjorie Bonneton<sup>13</sup>, Alice Thebaut<sup>1,2</sup>, Barbara Rohmer<sup>5</sup>, Florence Lacaille<sup>6</sup>, Pierre Broué<sup>8</sup>, Alexandre Fabre<sup>7</sup>, Karine Mention-Mulliez<sup>9</sup>, Jérôme Bouligand<sup>3,14,15</sup>, Emmanuel Jacquemin<sup>1,2,3,†</sup> and Emmanuel Gonzales<sup>1,2,3,†</sup>

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Almes, M.; Spraul, A.; Ruiz, M.; Girard, M.; Roquelaure, B.; Laborde, N.; Gottrand, F.; Turquet, A.; Lamireau, T.; Dabadie, A.; et al. Targeted-Capture Next-Generation Sequencing in Diagnosis Approach of Pediatric Cholestasis. *Diagnostics* **2022**, *12*, 1169. https://doi.org/10.3390/diagnostics12051169

Academic Editor: Consolato M. Sergi

Received: 15 March 2022

Accepted: 29 April 2022

Published: 7 May 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Background: Cholestasis is a frequent and severe condition during childhood. Genetic cholestatic diseases represent up to 25% of pediatric cholestasis. Molecular analysis by targeted-capture next generation sequencing (NGS) has recently emerged as an efficient diagnostic tool. The objective of this study is to evaluate the use of NGS in children with cholestasis. Methods: Children presenting cholestasis were included between 2015 and 2020. Molecular sequencing was performed by targeted capture of a panel of 34 genes involved in cholestasis and jaundice. Patients were classified into three categories: certain diagnosis; suggested diagnosis (when genotype was consistent with phenotype for conditions without any available OMIM or ORPHANET-number); uncertain diagnosis (when clinical and para-clinical findings were not consistent enough with molecular findings). Results: A certain diagnosis was established in 169 patients among the 602 included (28.1%). Molecular studies led to a suggested diagnosis in 40 patients (6.6%) and to an uncertain diagnosis in 21 patients (3.5%). In 372 children (61.7%), no molecular defect was identified. Conclusions: NGS is a useful diagnostic tool in pediatric cholestasis, providing a certain diagnosis in 28.1% of the patients included in

<footnote n="1">¹ Pediatric Hepatology and Pediatric Liver Transplant Unit, Reference Center for Biliary Atresia and Genetic Cholestasis, FSMR FILFOIE, ERN RARE LIVER, Bicêtre Hospital, Assistance Publique-Hôpitaux de Paris (AP-HP), Paris-Saclay University, 94270 Le Kremlin-Bicêtre, France; alicethebaut@aphp.fr (A.T.); emmanuel.jacquemin@aphp.fr (E.J.); emmanuel.gonzales@aphp.fr (E.G.)</footnote>

<footnote n="2">² INSERM UMR-S 1193, Paris-Saclay University, Hépatinov, 91400 Orsay, France; anne.spraul@aphp.fr</footnote>

<footnote n="3">³ Plateforme d'Expertise Maladies Rares Paris-Saclay, AP-HP, 94270 Le Kremlin-Bicêtre, France; jerome.bouligand@aphp.fr</footnote>

<footnote n="4">⁴ Biochemistry Unit, DMU15, Bicêtre Hospital, AP-HP, Paris-Saclay University, 94270 Le Kremlin-Bicêtre, France</footnote>

<footnote n="5">5 Pediatric Gastroenterology, Hepatology and Nutrition Unit, Bron Hospital, 69677 Lyon, France; mathias.ruiz@chu-lyon.fr (M.R.); barbara.rohmer@chu-lyon.fr (B.R.)</footnote>

<footnote n="6">6 Pediatric Gastroenterology, Hepatology and Nutrition Unit, Necker Hospital, 75015 Paris, France; muriel.girard@aphp.fr (M.G.); florence.lacaille@aphp.fr (F.L.)</footnote>

<footnote n="7">7 Pediatric Gastroenterology, Hepatology and Nutrition Unit, Marseille University Hospital, 13288 Marseille, France; bertrand.roquelaure@ap-hm.fr (B.R.); alexandre.fabre@ap-hm.fr (A.F.)</footnote>

<footnote n="8">8 Pediatric Gastroenterology, Hepatology and Nutrition Unit, Toulouse Hospital, 31300 Toulouse, France; laborde.n@chu-toulouse.fr (N.L.); broue.p@chu-toulouse.fr (P.B.)</footnote>

<footnote n="9">9 Department of Pediatric Gastroenterology Hepatology and Nutrition, Univ. Lille, CHU Lille, INSERM U1286, 59000 Lille, France; frederic.gottrand@chru-lille.fr (F.G.); karine.mulliez@chru-lille.fr (K.M.-M.)</footnote>

<footnote n="10">10 Pediatric Gastroenterology, Hepatology and Nutrition Unit, Saint-Denis Hospital, 97405 La Réunion, France; anne.turquet@chu-reunion.fr</footnote>

<footnote n="11">11 Pediatric Gastroenterology, Hepatology and Nutrition Unit, Bordeaux Hospital, 33076 Bordeaux, France; thierry.lamireau@chu-bordeaux.fr</footnote>

<footnote n="12">12 Pediatric Gastroenterology, Hepatology and Nutrition Unit, Rennes Hospital, 35033 Rennes, France; alain.dabadie@chu-rennes.fr</footnote>

<footnote n="13">13 Pediatric Gastroenterology, Hepatology and Nutrition Unit, Nancy Hospital, 54035 Nancy, France; m.bonneton@chru-nancy.fr</footnote>

<footnote n="14">14 Molecular Genetics, Pharmacology and Hormonology Unit, Bicêtre Hospital, AP-HP, Paris-Saclay University, 94270 Le Kremlin-Bicêtre, France</footnote>

<footnote n="15">15 INSERM UMR 1185, Endocrine Physiology and Physiopathology, Paris-Saclay University, 94270 Le Kremlin-Bicêtre, France</footnote>

<footnote>* Correspondence: marionflorence.almes@aphp.fr</footnote>

<footnote>† These authors contributed equally to this work.</footnote>
this study. In the remaining patients, especially those with variants of uncertain significance, the imputability of the variants requires further investigations.

**Keywords:** genetic cholestasis; children; NGS; neonatal sclerosing cholangitis; PFIC; Alagille syndrome; transient neonatal cholestasis

# 1. Introduction

In children, cholestatic jaundice is a rare condition occurring with an estimated prevalence of 1/2500 term infants [1,2]. In older children, the prevalence is certainly lower but is not known. Several etiologies must be investigated in cholestatic children, depending on the patient's age. Genetic cholestatic diseases collectively represent up to 25% of these etiologies, and might prove difficult to distinguish one from another [2,3].

Many genes related to cholestatic diseases have been identified and continue to be reported, including those responsible for the different progressive familial intrahepatic cholestasis (PFIC), benign recurrent intrahepatic cholestasis (BRIC), Alagille syndrome (AGS), neonatal sclerosing cholangitis (NSC), or inborn errors of bile acids metabolism. Molecular analysis by targeted-capture next generation sequencing (NGS) has recently emerged as a powerful diagnostic tool, allowing sequencing of a panel of selected genes in a single exam [4,5] and reducing both the delay and the cost of the diagnosis. It can be used as the first-line molecular technique or when previous Sanger sequencing of one or more genes failed to identify a precise diagnosis [6].

The main objective of our study was to evaluate from a clinical point of view the performance of NGS using a gene panel dedicated to cholestasis in a large cohort of children with cholestasis. The secondary objective is to report the clinical and paraclinical presentation of patients presenting an unexpected phenotype with respect to genotype.

# 2. Patients and Methods

## 2.1. *Type of Study*

We conducted a retrospective, observational cohort study. Recruitment was multicentric (Besançon, Béziers, Bordeaux, Clermont-Ferrand, Créteil, Dijon, Lille, Lyon, Nancy, Nice, Orléans, Marseille, Mayotte, Montpellier, Paris-Bicêtre, Paris-Necker, Paris-Trousseau, La Réunion, Nice, Saint-Étienne, Rouen, Rennes, Reims, Strasbourg, Toulouse). All molecular studies were performed in Bicêtre hospital and analyzed by the same biologists (AS and JB).

## 2.2. *Patients*

We included all patients under 18 years old who underwent a NGS using a gene panel dedicated to cholestasis as part of the diagnostic workup of a cholestatic condition (chronic or transient), between January 2015 and October 2020. Patients whose primary indication for a NGS molecular study was not cholestasis (e.g., acute liver failure and vascular liver disorders) were excluded from the study, as were patients presenting with cholestasis due to biliary atresia.

Our population mostly included patients who were studied prospectively, as part of their initial cholestasis assessment. A minority of patients were included retrospectively after Sanger sequencing of one or more genes if no diagnosis was provided by this technique.

Cholestasis was defined by an increase of serum conjugated bilirubin above 17 µmol/L (1 mg/dL) as recommended [2], and/or an increase of serum bile acids above 15 µmol/L.

## 2.3. *Ethics*

This study complies with the General Data Protection Regulation (RGPD—Regulation (EU) No. 2016/679) and the French Data Protection Act. The research has been declared
and listed in the Data Processing Registry of the Assistance Publique-Hôpitaux de Paris under No. 20191112110112.

The child's referring physician obtained written informed consent of the patient's parents, in accordance with the local guidelines established by the hospital's ethics committee.

## 2.4. Data Collection

When available, the following data as part of the diagnosis work up and of the follow-up were collected from patients' paper or computerized medical records.

Demographic data included date of birth, sex, age at onset of cholestasis, and age at the time of sampling. Clinical data included prematurity, perinatal anoxo-ischemia, infections, or suspected early-onset neonatal bacterial infections requiring antibiotics, intrauterine growth restriction (IUGR), stool color, hepatomegaly, splenomegaly, and abnormal extra-hepatological features. Biological parameters collected included: liver tests (aspartate amino transferase (AST), alanine amino transferase (ALT), serum gamma-glutamyl transferase (GGT) activity, total and conjugated bilirubin, serum bile acid (sBA) levels), coagulation tests (clotting factor V, prothrombin time (PT) and alpha-fetoprotein.

Imaging studies include mainly abdominal ultrasound (US) and in some patients MRI-cholangiography or cholecystography. Histological analyses of liver biopsy (LB) specimens were recorded.

## 2.5. Clinical Presentation Categories

We defined four categories of clinical presentation, based on the latest clinical and paraclinical (biological, radiological, histological) data available: chronic cholestasis (CC), transient neonatal cholestasis (TNC), benign recurrent intrahepatic cholestasis (BRIC) and oestroprogrestative-induced cholestasis (CholOP). CC was considered if the patient presented chronic liver disease, based on abnormal liver function tests (LFT) or abnormal abdominal imaging persisting after at least one year of treatment with ursodeoxycholic acid (UDCA). TNC definition criteria were: (1) normalization of LFT, with or without UDCA treatment, before the age of one year, and (2) persisting normal liver tests and imaging assessed after the age of one year and at least (3 months after discontinuation of UDCA treatment. BRIC was defined as recurrent episodes of cholestasis, with normalization of LFT between episodes. CholOP was defined as oestroprogrestative-induced episodes of cholestasis, resolutive after treatment discontinuation.

## 2.6. Sequencing Technique

### 2.6.1. Panel Design

A set of genes known to be involved or likely involved in various pediatric liver diseases, including 34 genes involved in cholestasis and jaundice, was selected to form the gene panel. Table 1 lists the 34 genes of the panel known to be involved in genetic cholestasis or jaundice.

<table><caption>**Table 1.** Genes of the panel involved in genetic cholestasis or jaundice.</caption><tr><td>ABCB11 (NM_003742.2)</td><td>ABCB4 (NM_000443.3)</td><td>ABCC2 (NM_000392.3)</td><td>ACOX2 (NM_003500) *</td><td>AIRE (NM_000383)</td></tr><tr><td>AKR1D1 (NM_005989.3)</td><td>AMACR (NM_014324.5)</td><td>ATP7B (NM_000053)</td><td>ATP8B1 (NM_005603.4)</td><td>BAAT (NM_001701.3)</td></tr><tr><td>CFTR (NM_000492.3)</td><td>CIRH1A (NM_032830.2)</td><td>CLDN1 (NM_021101.4)</td><td>CYP27A1 (NM_000784.3)</td><td>CYP7B1 (NM_004820.3)</td></tr><tr><td>DCDC2 (NM_016356) *</td><td>GBE1 (NM_000158)</td><td>GPBAR1 (NM_001077191)</td><td>HSD3B7 (NM_025193.3)</td><td>JAG1 (NM_000214.2)</td></tr><tr><td>MYO5B (NM_001080467)</td><td>NOTCH2 (NM_024408.3)</td><td>NR1H4 (NM_00512)</td><td>SERPINA1 (NM_000295) *</td><td>SCYL1 (NM_020680) *</td></tr><tr><td>SLCO1B1/SLCO1B3 (NM_006446/NM_019844)</td><td>SLC25A13 (NM_014251.2)</td><td>SLC27A5 (NM_012254.2)</td><td>TJP2 (NM_004817.3)</td><td>TTC37 (NM_014639.3)</td></tr><tr><td>UGT1A1 (NM_000463.2)</td><td>UNC45A (NM_018671)</td><td>VIPAS39 (NM_022067.3)</td><td>VPS33B (NM_018668.3)</td><td></td></tr></table>

<footnote>* Only included in some versions of the panel (see detailed versions of each gene panel in Figure S1 of the Supplementary Materials).</footnote>
The design of this panel was realized via the Agilent design platform (Agilent, Santa Clara, CA, USA). We adjusted the gene panel 4 times during this work, based on the data available in the literature. In this work, we report the results of versions 1 to 4 of the panel (Version 1: from January 2015 to November 2016, Version 2: from December 2016 to February 2018, Version 3: from March to April 2019, Version 4: from May 2019 to October 2020). The different versions of the panel used during this work are presented in Figure S1.

### 2.6.2. Technical Validation of the Targeted Gene Panel NGS

Thirty patients with CC, in whom polymorphisms in one of the genes included in the panel (*ABCB4*, *ABCB11*, *ATP8B1*, *JAG1*, *NOTCH2*, *ATP7B*) had been previously identified by Sanger sequencing, were used as controls to validate the NGS technique. These patients were not included in the rest of the study.

### 2.6.3. DNA Extraction

Patients' DNA was extracted from circulating leukocytes previously obtained from peripheral blood, and purified using a QiaSymphony MidiKit® kit. It was then fragmented by sonication using a Covaris M220 kit (Covaris Inc., MS, Woburn, MA, USA).

### 2.6.4. Libraries Constitution

The Illumina® libraries of the entire genome were prepared using the NEBNext® DNA Library Prep preparation kit for Illumina® (NEB Inc., Ipswich, MA, USA) on a Biomek Span 8 workstation (Beckman, Villepinte, France). Enrichment was achieved by targeted capture of the coding exons of the panel genes and their intronic bases (+/-30 base pairs) by hybridization, using the SureselectXT kit (Agilent, Santa Clara, CA, USA) robotized on a Biomek 4000 substation (Beckman, Villepinte, France). Deep intronic regions were not covered by this analysis. Some poorly covered exons were verified in Sanger sequencing if the clinical suspicion was high.

### 2.6.5. Sequencing Device

Sequencing was performed on a MiSeq Illumina® medium flow sequencer.

## 2.7. Bioinformatics Analysis

### 2.7.1. Selection of Variants of Interest

Sequencing was considered for further analyses when a coverage superior to 20X was obtained in more than 95% of the targeted regions.

The first selection of non-single nucleotide variants (SNV) mutations was carried out automatically by the Galaxy® Bioinformatics Platform of Paris-Sud University. In a simplified way, the readings of the Fatsq end pairs were taken after alignment of the sequences with the reference human genome 19 (hg19) using the BWA-MEM 0.7.10 software. Variant selection was made using GATK 3.4–46 and SNV variants were annotated using Annovar software (version 2015) and Snpeff 4.0. This first selection allowed, for each gene, the extraction of variants whose allele frequency in the general population was less than 2% except for *CFTR* for which an allele frequency cut off of 5% was used. Copy number variations were analyzed using an in-house Python script that compared depth and coverage data generated by Picard Metrics, as previously reported [7].

A second phase of analysis was then conducted by a specialized biologist who took a specific look at the variants presented in the genes of interest, according to the clinical presentation of the patient, as described below.

### 2.7.2. Classification of Variants

All mutations are reported following Human Genome Variation Society nomenclature (http://varnomen.hgvs.org/, accessed on 14 March 2022).
The different variants were classified according to the American College of Medical Genetics (ACMG) classification: pathogenic (P, class 5), likely pathogenic (LP, class 4), undetermined significance (VOUS, class 3), likely benign (class 2), and benign (class 1) [8].

All variants, except CFTR variants, were studied through the VarSome® platform (, accessed on 14 March 2022), which provided a class of pathogenicity, according to the ACMG classification. The classification proposed by the platform was generally accepted, unless data supporting the pathogenicity of a variant was available in the literature but not included in the VarSome® database, or if the pathogenicity of a variant was supported by familial studies.

CFTR variants were classified according to the recommendations of the National Association of Practitioners of Molecular Genetics (ANPGM) and using the CFTR-France website database (https://cftr.iurc.montp.inserm.fr/cftr, accessed on 14 March 2022). They were sorted into 4 classes [9]: cystic fibrosis-causing (CF Causing), CFTR-related diseases causing (CFTR-RD causing), and variants of undetermined significance (VOUS) for CF or for CFTR-RD mutations.

## 2.8. Variants Reported

All P and LP variants were reported. VOUS were not reported unless they were relevant when confronted to phenotype and their pathogenicity was supported by literature or in silico studies (VOUS-likely pathogenic, VOUS-LP). Benign and likely benign variants were not reported.

When possible, allelic repartition was studied. Copy number variations detected by the NGS technique were confirmed using quantitative PCR, MLPA, or CGHarray.

## 2.9. Diagnostic Categories

Molecular data were compared with clinical and paraclinical findings, allowing to classify the patients into one of the following diagnostic categories.

Provided that molecular data were consistent with phenotypes, patients with biallelic P/LP gene variants for autosomal recessive diseases or P/LP variants for autosomal dominant diseases were classified as having a certain diagnosis for conditions with an available OMIM or ORPHANET-number, or as having a suggested diagnosis for conditions without any available OMIM or ORPHANET-number, namely transient neonatal cholestasis, oestroprogestative-induced cholestasis and chronic or recurrent cholestasis due to monoallelic ABCB4 deficiency.

Patients were classified as having an uncertain diagnosis when phenotype was not in accordance with the one expected given the molecular findings, or when the molecular findings did not allow the patient to be classified as having a certain or suggested diagnosis (because allelic repartition was not studied and/or because the variant identified was a VOUS).

When the NGS study did not identify any P or LP variants, patients were classified as having no diagnosis.

# 3. Results

## 3.1. Patients

Among the 706 children who were studied through the NGS panel between January 2015 and October 2020, 104 were excluded because they presented biliary atresia or primary non-cholestatic liver disease.

Among the 602 patients included, 309 (51.3%) were under six months old; 357 (59.3%) were male. Two hundred patients were recruited directly by our center (33.2%), while the rest of them were sampled in their local reference center and referred only for genetic testing.

Among the 602 patients included, 65 had been previously studied by Sanger Sequencing.
## 3.2. Technical Validation Phase of Gene Panel

All the polymorphisms identified by Sanger sequencing in the 30 control patients were also identified by targeted NGS.

## 3.3. Depth and Coverage of Sequenced Regions

The average reading depth of the targeted regions was 136X (SD 20.6). Coverage greater than 20X was obtained for more than 98% of the targeted regions.

## 3.4. Patients with a Certain Diagnosis

Targeted gene panel NGS established a certain diagnosis of OMIM/ORPHANET-identified diseases in 169 out of the 602 patients (28.1%), including 24 patients who had previously been studied by Sanger Sequencing (Table 2). All of them had an initial clinical presentation, a biological and clinical evolution consistent with the results of NGS sequencing.

<table><caption>**Table 2.** Certain molecular diagnosis of cholestasis or hyperbilirubinemia as determined by the NGS panel (n = 170).</caption><thead><tr><th>Gene</th><th>Transmission Mode</th><th>Disease (OMIM/Orphanet Number)</th><th>Number of Patients</th></tr></thead><tbody><tr><td colspan="4" style="text-align: center; font-weight: bold;">Cholestasis with High GGT Activity</td></tr><tr><td>JAG1</td><td>AD</td><td>Alagille syndrome (118450/261619)</td><td>92</td></tr><tr><td>NOTCH2</td><td>AD</td><td>Alagille syndrome (610205/261629)</td><td>48</td></tr><tr><td>ABCB4</td><td>AR</td><td>PFIC3 (602347/79305)</td><td>7</td></tr><tr><td>DCDC2</td><td>AR</td><td>Neonatal sclerosing cholangitis (617394/480556)</td><td>17</td></tr><tr><td>CFTR</td><td>AR</td><td>Cystic fibrosis (219700/586)</td><td>9</td></tr><tr><td>ATP7B</td><td>AR</td><td>Wilson disease (277900/905)</td><td>3</td></tr><tr><td>SERPINA1</td><td>AR</td><td>Alpha-1-antitrypsine deficiency (613490/60)</td><td>5</td></tr><tr><td colspan="4" style="text-align: center; font-weight: bold;">Cholestasis with Normal or LOW GGT Activity</td></tr><tr><td>ATP8B1</td><td>AR</td><td>PFIC1 (211600/79306)</td><td>1</td></tr><tr><td>ATP8B1</td><td>AR</td><td>BRIC1 (243300/99960)</td><td>1</td></tr><tr><td>ABCB11</td><td>AR</td><td>PFIC2 (601847/79304)</td><td>19</td></tr><tr><td>ABCB11</td><td>AR</td><td>BRIC2 (605479/99961)</td><td>3</td></tr><tr><td>TJP2</td><td>AR</td><td>PFIC4 (615878/480483)</td><td>8</td></tr><tr><td>NR1H4</td><td>AR</td><td>PFIC5 (617049/480476)</td><td>1</td></tr><tr><td>MYO5B</td><td>AR</td><td>Myosin 5b deficiency related cholestasis (in absence of MVID) (ND/480491)</td><td>10</td></tr><tr><td>VPS33B</td><td>AR</td><td>ARC syndrome (208085/2697)</td><td>5</td></tr><tr><td>AKR1D1</td><td>AR</td><td>Primary bile acid synthesis defect (235555/79303)</td><td>1</td></tr><tr><td>ACOX2</td><td>AR</td><td>Primary bile acid synthesis defect (617308/ND)</td><td>2</td></tr><tr><td>BAAT</td><td>AR</td><td>Primary bile acid synthesis defect (conjugation defect) (619232/238475)</td><td>1</td></tr><tr><td>CYP27A1</td><td>AR</td><td>Cerebrotendinous xanthomatosis (213700/909)</td><td>1</td></tr><tr><td>SCYL1</td><td>AR</td><td>Spinocerebellar ataxia (616719/466794)</td><td>3</td></tr><tr><td>TTC37</td><td>AR</td><td>Tricho-hepato-enteric syndrome (614602/84064)</td><td>1</td></tr><tr><td colspan="4" style="text-align: center; font-weight: bold;">Genetic Hyperbilirubinemia</td></tr><tr><td>UGT1A1</td><td>AR</td><td>Gilbert syndrome (143500/357)</td><td>1</td></tr><tr><td>UGT1A1</td><td>AR</td><td>Crigler-Najjar syndrome (218800/79235)</td><td>1</td></tr><tr><td>ABCC2</td><td>AR</td><td>Dubin-Johnson syndrome (237500/234)</td><td>17</td></tr><tr><td>SLCO1B1/B3</td><td>AR</td><td>Rotor syndrome (237450/3111)</td><td>1</td></tr></tbody></table>

AD: autosomal dominant; AR: autosomal recessive, ARC: Arthrogryposis—renal failure—cholestasis, BRIC: benign recurrent intrahepatic cholestasis, microvillous inclusion disease, ND: not defined; PFIC: progressive familial intrahepatic cholestasis (PFIC).

Among the 92 patients who presented with high GGT activity cholestasis, 55 were diagnosed with Alagille syndrome, bearing either a *JAG1* mutation (n = 48) or a *NOTCH2* mutation (n = 42).
mutation (n = 7), 17 patients were diagnosed with PFIC3, nine with neonatal sclerosing cholangitis related to DCDC2 deficiency, three with cystic fibrosis, three with Wilson disease and five with alpha one anti-trypsin (A1AT) deficiency.

Among the 57 patients who presented with low or GGT normal activity two were diagnosed with FIC1 deficiency (PFIC1, n = 1; BRIC, n = 1), 22 were diagnosed with BSEP deficiency (PFIC2, n = 19; BRIC2, n = 3), eight patients were diagnosed with TJP2 deficiency (PFIC4), 10 patients were diagnosed with isolated cholestasis related to MYO5B deficiency, five patients with an arthrogryposis—renal failure—cholestasis syndrome (ARC syndrome, related to VPS33B deficiency), primary bile acid synthesis defects due to acyl-CoA oxidase deficiency was diagnosed in two patients (ACOX2), and spinocerebellar ataxia was diagnosed in three patients (SCYL1). The following diseases were all diagnosed in a single patient: NR1H4 deficiency (PFIC5), primary bile acid synthesis deficiency due to delta 4-3-oxoid dehydrogenase deficiency (AKR1D1), primary bile acid synthesis deficiency due to bile acid conjugation defects (BAAT), cerebrotendinous xanthomatosis (CYP27A1) and tricho-hepato-enteric syndrome (TTC37 gene).

In addition, 20 patients were diagnosed with a genetic hyperbilirubinemia syndrome: Dubin-Johnson syndrome (ABCC2, n = 17), Rotor syndrome (SLCO1B1/B3, n = 1), Crigler-Najjar syndrome (UGT1A1, n = 1) and Gilbert syndrome (UGTA1, n = 1). Of note, all patients diagnosed with Dubin-Johnson syndrome presented with TNC.

Table S1 in Supplementary Materials reports all causal mutations identified in these patients, as well as other pathogenic variants found in other genes of interest.

### 3.5. *Patients with a Suggested Diagnosis*

A suggested diagnosis was proposed for 40 patients (6.6%), for whom molecular studies were consistent with known cholestatic conditions without any available OMIM/ORPHANET-number.

A TNC related to pathogenic or likely pathogenic monoallelic mutations in the main PFIC genes was identified in 22 patients (ATP8B1, n = 2; ABCB11, n = 12; ABCB4, n = 8) (Table 3). Importantly, usual risk factors for TNC were disclosed in only 5 of these 22 patients (neonatal CMV infection (n = 4) and mild perinatal anoxo-ischemia (n = 1)). Patients presenting TNC but in whom no mutation was found were not reported.

Fourteen patients presented with chronic or recurrent high GGT level cholestasis and one PV or LPV in ABCB4 (Table 3b). Of note, the c.1769G > A variant of ABCB4 was reclassified as a class 4 variant based on previous studies showing its pathogenic nature [10,11]. This variant was also identified in one compound heterozygous patient of our cohort allowing a certain diagnosis of PFIC3 (patient 187, Table S1 of the Supplementary Materials).

Four patients with oestroprogestative-induced cholestasis associated with one PV or LPV in ABCB11 are reported in Table 3c (n = 4).

### 3.6. *Patients with an Uncertain Diagnosis*

For 21 patients (3.5%), molecular studies led to an uncertain diagnosis (Table 4).

In ten patients, the phenotype was not consistent with the genotype (patients 41 to 50, Table 4). One patient presented with a PFIC3 phenotype, but the variants identified were not fully supportive of this diagnosis and the allelic repartition was not studied (patient 46). One patient presented one PV in HNF1B but his clinical presentation was not concordant with data previously reported (patient 44) [12]. Two patients presented a CC with heterozygous PV or LPV in two different genes (patients 47–48, Table 4). In patient 45 (Table 4), NGS revealed a homozygous VOUS variant in ATP7B, previously reported to cause an attenuated phenotype of Wilson disease [13]. He presented with CC, cirrhosis, and normal cupric balance but elevated cupric levels in the LB specimen. Patient 49 (Table 4) presented with a PFIC2 phenotype, but only one PV in ABCB11 was identified.
<table><caption>**Table 3. (a).** Suggested diagnosis of cholestatic conditions without OMIM/ORPHANET-number: patients with transient neonatal cholestasis due to pathogenic or likely pathogenic heterozygous mutations in one of the main PFIC genes (n = 22). (b). Suggested diagnosis of cholestatic conditions without OMIM/ORPHANET number: patients with chronic or recurrent cholestasis due to pathogenic or likely pathogenic heterozygous variants in ABCB4 (n = 14). (c). Suggested diagnosis of cholestatic conditions without any available OMIM/ORPHANET number: patients with oestroprogestative-induced cholestasis (n = 4).</caption><thead><tr><th rowspan="2">Phenotype</th><th rowspan="2">Gene</th><th colspan="3">a</th></tr><tr><th>Patient</th><th>Mutation</th><th>Notes</th></tr></thead><tbody><tr><td rowspan="8">ATP8B1</td><td rowspan="8"></td><td>1</td><td>c.3040C>T:p.Arg1014Ter</td><td>TJP2: p.Gly538Ala (He)—Pathogenic<br>MYO5B: p.Glu144Gly (He) -VOUS</td></tr><tr><td>2</td><td>Complete gene deletion</td><td></td></tr><tr><td>3</td><td>c.3691C>T:p.Arg1231Trp</td><td></td></tr><tr><td>4</td><td>c.1445A>G:p.Asp482Gly</td><td></td></tr><tr><td>5</td><td>c.2873_2874insCG:p.<br>Phe959GlyfsTer49</td><td></td></tr><tr><td>6</td><td>c.3329C>T:Ala1110Val</td><td></td></tr><tr><td>7</td><td>c.1243C>T:p.Arg415Ter</td><td></td></tr><tr><td>8</td><td>c.3181_3184del:p.Ile1061<br>ValfsTer35</td><td></td></tr><tr><td rowspan="8">Transient neonatal<br>cholestasis</td><td rowspan="8">ABCB11</td><td>9</td><td>c.2061C>G:p.Tyr687Ter</td><td>Neonatal cholestasis with pale stool and hypoglycemia,<br>normal cholangiogram, advanced fibrosis on LB (at the age<br>of 2 months), normalization of LFT at 1 year old with<br>UDCA treatment. Persistent normal LFT after UDCA<br>discontinuation.</td></tr><tr><td>10</td><td>c.3752C>T:p.Thr1251Ile</td><td></td></tr><tr><td>11</td><td>c.896-<br>897GA>TT:p.Arg299Ile</td><td></td></tr><tr><td>12</td><td>c.890A>G:p.Glu297Gly</td><td>Inherited from his mother who presented ICP—Neonatal<br>CMV infection.</td></tr><tr><td>13</td><td>c.3148C>T:p.Arg1050Cys</td><td></td></tr><tr><td>14</td><td>c.1396C>A:p.Gln466Lys</td><td>Neonatal CMV infection.</td></tr><tr><td>15</td><td>c.140G>A:p.Arg47Gln</td><td></td></tr><tr><td>16</td><td>c.2800G>A:p.Ala934Thr</td><td>Neonatal CMV infection.</td></tr><tr><td rowspan="7">ABCB4</td><td rowspan="7"></td><td>17</td><td>c.2800G>A:p.Ala934Thr</td><td></td></tr><tr><td>18</td><td>c.2800G>A:p.Ala934Thr</td><td></td></tr><tr><td>19</td><td>c.101C>T:p.Thr34Met</td><td></td></tr><tr><td>20</td><td>c.1769G>A:p.Arg590Gln *<br>[10,11]</td><td>Mild perinatal anoxo-ischemia.</td></tr><tr><td>21</td><td>c.3280-2A>G:p(?)</td><td></td></tr><tr><td>22</td><td>c.3676T>G:p.Cys1226Gly</td><td></td></tr></tbody></table>
<table><caption>**Table 3. Cont.**</caption><thead><tr><th rowspan="2">Phenotype</th><th rowspan="2">Patient</th><th colspan="3">b</th></tr><tr><th>Mutation</th><th>Status</th><th>Notes</th></tr></thead><tbody><tr><td rowspan="10">Chronic or recurrent high GGT level cholestasis</td><td>23</td><td>ABCB4: c.1436C>T:p.Pro479Leu</td><td>He: Pathogenic</td><td>Recurrent episodes of cholestasis and pruritus. Normalization of LFT with UDCA treatment.</td></tr><tr><td>24</td><td>ABCB4: c.1769G>A:p.Arg590Gln*</td><td>He: Likely Pathogenic [10,11]</td><td>Anicteric cholestasis associated with abdominal pain. Normal abdominal US images, especially no lithiasis. Normalization of LFT with UDCA treatment.</td></tr><tr><td>25</td><td>ABCB4:c.1714C>T:p.Gln572Ter (inherited from the father)</td><td>He: Pathogenic</td><td>Chronic cholestasis and pruritus. Hydrocholecystis with possible compression of extrahepatic bile ducts. Family history of liver disease (father: LPAC syndrome).</td></tr><tr><td>26</td><td>ABCB4:c.1769G>A:p.Arg590Gln*</td><td>He: Likely Pathogenic [10,11]</td><td>Trisomy 21, neonatal cholestasis with intra vesicular sludge, cirrhosis and ascitis. Complete atrioventricular canal. Duodenal atresia and necrotizing enterocolitis requiring parenteral nutrition. Possible ischemic cholangiopathy due to complicated surgeries during neonatal period. Normalization of LFT at 1 year old, no discontinuation.</td></tr><tr><td>27</td><td>ABCB4:c.2800G>A:p.Ala934Thr</td><td>He: Likely pathogenic</td><td>Chronic cholestasis, cirrhosis.</td></tr><tr><td>28</td><td>ABCB4: c.574G>T:p.Val192Phe (inherited from the mother) ABCB4: c.3655-6C>G:p.(?) (de novo)</td><td>He: Likely pathogenic He: VOUS</td><td>Neonatal cholestasis, pale stool, normal cholangiogram, family history of liver disease (mother and maternal aunt: ICP). Partial improvement of LFT with UDCA treatment.</td></tr><tr><td>29</td><td>ABCB4: c.662T>C:p.Met221Thr (inherited from the father)</td><td>He: Likely pathogenic</td><td>Fallot's tetralogy, nephrocalcinosis, persistent abnormal LFT with UDCA treatment, normal cholangiogram, cirrhosis. No information on neonatal period. No family history of liver disease.</td></tr><tr><td>30</td><td>ABCB4:c.2800G>A:p.Ala934Thr (inherited from the mother) ATP8B1:c.2097+1G>A: p.(?) (inherited from the father)</td><td>He: Pathogenic He: Likely Pathogenic</td><td>Chronic cholestasis, cirrhosis, portal hypertension, aortic arch abnormalities. Family history of liver disease (mother: ICP).</td></tr><tr><td>31</td><td>ABCB4:c.2144C>T:p.Thr715Ile</td><td>He: Likely Pathogenic</td><td>Neonatal cholestasis, deceased during neonatal period (multiorgan failure).</td></tr><tr><td>32</td><td>ABCB4:c.1769G>A:p.Arg590Gln* [10,11]</td><td>He: Likely Pathogenic</td><td>Mental retardation, butterfly vertebrae. Lost of follow-up.</td></tr><tr><td>33</td><td>ABCB4:c.1769G>A:p.Arg590Gln* [10,11]</td><td>He: Likely Pathogenic</td><td>Chronic cholestasis, cirrhosis, liver histology compatible with neonatal sclerosing cholangitis, LT at 2 years old.</td></tr></tbody></table>
<table><caption>**Table 3. Cont.**</caption><tbody><tr><td>34</td><td>ABCB4:c.1769G&gt;<br/>A:p.Arg590Gln *<br/>[10,11]<br/>ABCB11:c.470A&gt;<br/>G:p.Tyr157Cys<br/>(inherited from<br/>the mother)</td><td>He: Likely Pathogenic<br/>He: VOUS-LP</td><td>Extreme prematurity with severe anoxo-ischemia and long term parenteral nutrition. Neonatal cholestasis. Deceased during neonatal period (multiorgan failure).<br/>Post mortem LB: panlobular dissecting fibrosis, ductopenia, hepatocellular cholestasis, disseminated intravascular coagulopation. Possible consequences of a submassive destruction of liver parenchyma.<br/>Family history of liver disease (mother: ICP; twin brother: TNC (carrying only ABCB11 variant); healthy older sister (same genotype): asymptomatic vesicular lithiasis).</td></tr><tr><td>35</td><td>ABCB4:c.959C&gt;<br/>T:p.Ser320Phe</td><td>He: Likely Pathogenic</td><td>Complex cardiac malformation, chronic cholestasis, exsudative enteropathy.</td></tr><tr><td>36</td><td>ABCB4:c.140G&gt;<br/>A:p.Arg47Gln</td><td>He: Pathogenic</td><td>Patient also diagnosed with a glycogenose storage disease type III and deafness.</td></tr><tr><td colspan="4" style="text-align: center;">c</td></tr><tr><th>Phenotype</th><th>Gene</th><th>Patient</th><th>Mutation</th><th>Status</th></tr><tr><td rowspan="4">Oestro-<br/>progestative<br/>induced<br/>cholestasis</td><td rowspan="4">ABCB11</td><td>37</td><td>c.890A&gt;G:p.Glu297Gly<br/>(He)<br/>(inherited from her<br/>mother)<br/>c.403G&gt;A:p.Glu135Lys<br/>(He)<br/>(de novo or inherited from<br/>her father)</td><td>CHe</td></tr><tr><td>38</td><td>c.150+3A&gt;C:p(?)</td><td>He</td></tr><tr><td>39</td><td><b>c.3130delG:p.Ala1044<br/>LeufsTer53</b></td><td>He</td></tr><tr><td>40</td><td><b>c.477+6T&gt;G:p(?)</b></td><td>Ho</td></tr></tbody></table>

ICP: intrahepatic cholestasis of pregnancy, LFT: liver function tests, UDCA: ursodeoxycholic acid. The references cited were used to reclassify the corresponding variant from class 3 to class 4. * The pathogenicity of this variant is supported by previous publications [10,11]. Previously unreported variants according to ClinVar® and Varsome® are indicated in bold. BA: bile acids, ICP: intrahepatic cholestasis of pregnancy, He: heterozygous, LB: liver biopsy, LFT: liver function tests, UDCA: ursodeoxycholic acid, TNC: transient neonatal cholestasis, US: ultrasound. The references cited were used to reclassify the corresponding variant from class 3 to class 4. * The pathogenicity of this variant is supported by previous publications [10,11]. Previously unreported variants according to ClinVar and Varsome are indicated in bold. CHe: compound heterozygous, He: heterozygous, Ho: homozygous. Previously unreported variants according to ClinVar® and Varsome® are indicated in bold.

Eleven patients were carriers for CFTR RD variants, including one patient presenting a homozygous status (Table 4). Among them, four had a transient neonatal cholestasis and seven presented with CC. Of note, no variant (P, LP, VOUS) in the other genes included in the panel was identified in these patients. In addition, no usual risk factors for TNC were disclosed in the four patients with TNC.

### 3.7. *Patients with No Diagnosis*

372 patients (61.8%) remained undiagnosed after performing NGS sequencing.
<table><caption>**Table 4. (a).** Uncertain diagnosis of cholestatic conditions (n = 10). (b). Uncertain molecular diagnosis: carriers for CFTR-related disease (RD) or CF-causing mutations (n = 11).</caption><thead><tr><th colspan="4">a</th></tr><tr><th>Patient</th><th>Genotype</th><th>ACMG Classification</th><th>Phenotype</th></tr></thead><tbody><tr><td>41</td><td>NOTCH2:c.1396C>T:p.Gln466Ter (He)</td><td>Pathogenic</td><td>Chronic cholestasis, no extra-hepatic AGS feature, neonatal occlusion, neonatal CMV infection, no bile duct paucity in a LB containing only 6 portal spaces.</td></tr><tr><td>42</td><td>NOTCH2:c.5002+2T>C:p.(?) (He)</td><td>Pathogenic</td><td>Transient neonatal cholestasis, No extra-hepatic AGS feature. No LB, no UDCA treatment.</td></tr><tr><td>43</td><td>NOTCH2:c.6101C>A:p.Ala2034Asp (He)</td><td>VOUS-LP</td><td>Chronic cholestasis, No extra-hepatic AGS feature, no bile duct paucity on LB.</td></tr><tr><td>44</td><td>HNF1B:1.32Mb deletion including HNF1B (He)</td><td>Pathogenic</td><td>Acute hepatitis with anicteric cholestasis, spontaneous resolution, no relapse (3 years of follow-up). No diabetes. Normal renal function. Slight perisinusoidal fibrosis, moderate hyperplasia of stellate cells and normal bile ducts in LB. Normal liver and renal ultrasound.</td></tr><tr><td>45</td><td>ATP7B:c.4135C>T:p.Pro1379Ser* (Ho) [13]</td><td>VOUS</td><td>Chronic cholestasis, cirrhosis, normal cupric balance but elevated cupric levels in LB.</td></tr><tr><td>46</td><td>ABCB4:c.2581T>G:p.Leu861Val (He)<br>ABCB4:c.1584G>C:p.Glu528Asp* (Ho) [14]</td><td>VOUS-LP<br>VOUS</td><td>Chronic cholestasis, cirrhosis on LB (MDR3 staining: not available).</td></tr><tr><td>47</td><td>NR1H4: c.920-2A>G:p.(?) (He)<br>DGUOK: c.353G>A:p.Arg118His (He)</td><td>Pathogenic<br>Likely Pathogenic</td><td>Neonatal cholestasis with hypoglycemia; steatosis and hepatocellular cholestasis on LB. Treated with UDCA: normal LFT and abdominal US at age 2 years. No discontinuation of UDCA.</td></tr><tr><td>48</td><td>MYO5B: c.3190C>T:p.Arg1064Ter (He)<br>(inherited from the father)<br>CFTR: c.3485G>T p.Arg1162Leu (He)<br>(inherited from the mother)</td><td>Pathogenic<br>RD CAUSING</td><td>Pruritus, chronic diarrhea, chronic bronchial congestion, normal cholangiogram.</td></tr><tr><td>49</td><td>ABCB11:c.1408C>T:p.Arg470Ter (He)</td><td>Pathogenic</td><td>Chronic cholestasis persistent with UDCA treatment, cirrhosis.</td></tr><tr><td>50</td><td>ABCB11:c.890A>G:p.Glu297Gly (He)</td><td>Likely Pathogenic</td><td>Neonatal cholestasis. Lost of follow-up.</td></tr><tr><th colspan="4">b</th></tr><tr><th>Patient</th><th>Genotype</th><th>CFTR-France Classification</th><th>Phenotype</th></tr><tr><td>51</td><td>CFTR:[c.2002C>T:p.Arg668Cys<br>c.1727G>C:p.Gly576Ala] (He)</td><td>RD CAUSING<br>RD CAUSING</td><td rowspan="3">Transient neonatal cholestasis</td></tr><tr><td>52</td><td>CFTR:c.2991G>C:p.Leu997Phe (He)</td><td>RD CAUSING</td></tr><tr><td>53</td><td>CFTR:c.3909C>G:p.Asn1303Lys (He)</td><td>CF-causing</td></tr><tr><td>54</td><td>CFTR:[c.220C>T:p.Arg74Trp;<br>c.3808G>A:p.Asp1270Asn] (Ho)</td><td>RD-VOUS4<br>RD-VOUS4</td><td rowspan="3">Chronic cholestasis</td></tr><tr><td>55</td><td>CFTR:[c.220C>T:p.Arg74Trp;<br>c.3808G>A:p.Asp1270Asn; c.601G>A:<br>p.Val201Met] (He)</td><td>RD-VOUS4<br>RD-VOUS4<br>RD-VOUS4</td></tr><tr><td>56</td><td>CFTR: c.1523T>G:p. Phe508Cys (He)<br>[c.1727G>C:p.Gly576Ala;<br>c.2002C>T:p.Arg668Cys] (He)</td><td>RD CAUSING<br>RD CAUSING<br>RD CAUSING</td></tr><tr><td>57</td><td>CFTR:c.2991G>C:p.Leu997Phe (He)</td><td>RD CAUSING</td><td></td></tr><tr><td>58</td><td>CFTR:c.2991G>C:p.Leu997Phe (He)</td><td>RD CAUSING</td><td></td></tr><tr><td>59</td><td>CFTR:c.3485G>T:p.Arg1162Leu (He)</td><td>RD CAUSING</td><td></td></tr><tr><td>60</td><td>CFTR:c.2173G>A:p.Glu725Lys (He)</td><td>RD CAUSING</td><td></td></tr><tr><td>61</td><td>CFTR: c.1516A>G:p.Ile506Val (He)</td><td>RD-VOUS4</td><td></td></tr></tbody></table>

AGS: Alagille syndrome, He: heterozygous, Ho: homozygous, VOUS: variant of uncertain significance, VOUS-LP: variant of uncertain significance-likely pathogenic, RD CAUSING: CFTR-related disease causing, LB: liver biopsy, LFT: liver function tests, UDCA: ursodeoxycholic acid. *, the pathogenicity of this variant is supported by a previous publication [13], ° the pathogenicity of this variant is supported by a previous publication [14]. Previously unreported variants according to ClinVar® and Varsome® are indicated in bold. He: heterozygous, VOUS4: VOUS-likely pathogenic (CFTR specific denomination).
## 4. Discussion

Genetic cholestasis encompasses an increasing number of rare and severe diseases. Clinical findings and biological presentation are often similar. An early diagnosis in these patients may offer more specific therapeutics, delay invasive diagnostic procedures (biliary tree opacification, LB), and help clinicians with genetic counseling. Thus, NGS seems to be a first choice tool in the diagnostic process of these diseases [3,15].

NGS sequencing of genes involved—or potentially involved—in pediatric cholestatic diseases in children allowed us to establish a certain diagnosis in 169 patients (28.1%). This diagnostic performance rate was similar to those previously reported in literature (Table 5) [15–20]. Two studies reported significantly higher diagnostic performance rates of 61% and 68% [21,22]. However, the first study was performed in a population with a high inbreeding rate [21] while the second one reported on a population mainly made of retrospective patients that had not been previously studied by Sanger sequencing of the classical PFIC genes [22]. Among the 602 patients included in our study, 65 had been previously studied using the Sanger technique for one or more genes, and had remained undiagnosed. NGS sequencing established a diagnosis for 24 of them (36.9%) (ABCB4, n = 2; TJP2, n = 6; JAG1, n = 2; DCDC2, n = 2; ABCC2, n = 1; AKR1D1, n = 1; ATP7B, n = 1; CYP27A1, n = 1; MYO5B, n = 5; NR1H4, n = 1; VPS33B, n = 2). If we only consider the 537 patients in whom NGS was used as first-line diagnosis approach, the diagnostic performance rate decreased slightly to 27%.

<table><caption>**Table 5.** Results of NGS sequencing in pediatric patients in diagnosis approach of cholestasis: data already published.</caption><thead><tr><td>Reference</td><td>Number of Genes</td><td>Number of Patients</td><td>Age (Years)</td><td>NGS Indication</td><td>Diagnostic Rate (%)</td></tr></thead><tbody><tr><td>Matte U et al. JPGN. 2010</td><td>5</td><td>51</td><td>0–17</td><td>Chronic cholestasis</td><td>27%</td></tr><tr><td>Wang NL et al. PLoS One. 2016</td><td>61</td><td>141</td><td>0–17</td><td>Chronic cholestasis</td><td>22%</td></tr><tr><td>Togawa T et al. J Pediatr. 2016</td><td>18</td><td>109</td><td>&lt;1</td><td>Chronic cholestasis</td><td>26%</td></tr><tr><td>Stalke A et al. Clin Genet. 2018</td><td>21</td><td>135</td><td>0–20</td><td>Chronic cholestasis</td><td>17%</td></tr><tr><td>Chen HL et al. J. Of Ped. 2018</td><td>52</td><td>102</td><td>0–18</td><td>Cholestasis</td><td>32,4%</td></tr><tr><td>Shagrani et al. Clin Genet. 2018</td><td>189</td><td>98</td><td>0–17</td><td>Severe cholestasis</td><td>61%</td></tr><tr><td>Lipinski et al. Front. Pediatr. 2020</td><td>53</td><td>22</td><td>0–18</td><td>Chronic cholestasis</td><td>68%</td></tr><tr><td>Karpen et al. JPGN. 2021</td><td>66</td><td>2171</td><td>0–18</td><td>Cholestasis</td><td>12%</td></tr><tr><td><strong>Our study</strong></td><td><strong>34</strong></td><td><strong>603</strong></td><td><strong>0–17</strong></td><td><strong>Cholestasis</strong></td><td><strong>28%</strong></td></tr></tbody></table>

Among the 169 patients in whom a certain diagnosis was provided by NGS, we reported patients with an atypical phenotype with respect to genotype. One patient was diagnosed with PFIC3 although he was considered as having a NSC since the first year of life, on the basis of a MRI-cholangiography showing an abnormal intrahepatic biliary tree (patient 188, Tables S1 and S2, Figure S2a). Of note, no variant in the genes involved in NSC was identified in this patient. Interestingly, in our experience, abnormal biliary imaging (transhepatic cholecystography) was observed in one additional PFIC3 patient (not included in this study but detailed in Table S2 and Figure S2b. One hypothesis to explain this abnormal aspect of the biliary tree mimicking NSC, is that MDR3 deficiency reduces phospholipid levels in patient bile, resulting in toxic bile that leads to biliary injuries as has been observed in MDR2-/– mice [23,24]. Such secondary sclerosing cholangitis has not been reported so far in the PFIC3 series [25,26] but has been described in LPAC patients with MDR3 deficiency [27,28].

One patient was diagnosed as a teenager with ARC syndrome (patient 124, Table S1), while previously considered as having an idiopathic CC, associated with tubulopathy and benign osseous abnormalities (talus foot, congenital luxation of both hips). He underwent a LT at four years old before the diagnosis of ARC syndrome was made. Renal dysfunc-
tion worsened after LT. Additional explorations, performed after the molecular diagnosis was made, showed perception deafness, frontal syndrome, and thrombopathy (agranular platelets). He died at age 30 of chronic respiratory failure, probably due to long-term medication toxicity (mycophenolate mofetil). Such hypomorphic forms of ARC syndrome have been previously reported [29].

Interestingly, in one patient (patient 82, Table S1) presenting with a clinical diagnosis of AGS, targeted gene panel NGS identified not only a PV in *JAG1* but also one VOUS-LP and one PV in *MYO5B*. While the allelic repartition was not studied, MYO5B immunostaining showed an abnormal pattern consistent with a diagnosis of MYO5B deficiency associated with a *JAG1*-related AGS (unpublished personal data, E. Gonzales). In another patient (patient 218, Table S1), NGS revealed a compound heterozygous status for pathogenic variants both in *MYO5B* and *ATP8B1*. Myo5B deficiency and PFIC1 are both responsible for low GGT chronic cholestasis and possible gastrointestinal manifestations, it is difficult to determine the contribution of each molecular defect to the global phenotype of this patient. No Myo5B immunostaining had been performed in this patient.

Forty patients (6.6%) of the study were classified as having a suggested diagnosis based on a consistent genotype with respect to a phenotype that is not OMIM/ORHANET labeled. In adults, it has been reported that monoallelic PV or LPV variant in *ABCB4* can be associated with CC and cirrhosis consistent with an attenuated phenotype of PFIC3 [10,25]. Interestingly, in this study, we report on 14 children presenting with chronic or recurrent cholestasis with such genotype. Phenotypes ranged from recurrent episodes of high GGT cholestasis to high GGT chronic cholestasis leading to cirrhosis. These observations are in line with those of a recent study and suggest that monoallelic P or LP *ABCB4* variants can also be responsible for chronic cholestatic disease in children [30,31]. We also reported on four female teenagers who presented a first episode of normal GGT cholestasis when started with an oestroprogestative drug and in whom a mono-allelic P or LP variant of *ABCB11* was identified. Such oestroprogestative induced cholestasis has been reported in female patients carrying a monoallelic variant in *ABCB4* [25] and *ABCB11* [32]. These reports highlight the need for specialized gynecological management in collaboration with hepatologists not only in female patients diagnosed with a PFIC but also in female patients known to carry a P or an LP variant in the genes involved in genetic cholestasis especially *ABCB11* and *ABCB4*. Reciprocally, these data suggest that molecular analysis of the genes involved in cholestasis should be considered in female patients presenting a bout of cholestasis after having been started with an oestroprogestative drug. While most of the published cohorts included patients with severe and/or persistent cholestasis [16,17,19,21], our study also included patients presenting with a TNC. Many non-genetic risk factors of TNC have been reported [33,34], and few studies have suggested that monoallelic pathogenic variants in genes involved in genetic cholestasis could contribute to TNC physiopathology [35–39]. Here we report 22 patients presenting with a TNC in whom one monoallelic pathogenic variant of *ABCB11* (*n* = 12), *ABCB4* (*n* = 8), *ATP8B1* (*n* = 2) was identified. Of note, 17 out of these 22 patients did not present any identifiable non-genetic risk factors for CNT. These findings further suggest the role of genetic variants in the pathophysiology of CNT. Whether or not the TNC features of these patients are different from those without genetic risk factors of TNC remains to be studied. In any case, such patients should be discussed with tertiary care center as they may benefit from specific treatment (e.g., prolonged treatment with UDCA) and follow-up (e.g., during pregnancy or when considering oestroprogestative therapy, the possible occurrence of biliary lithiasis and even slow progression fibrosis, as discussed above).

Among the 21 patients (3.5%) classified as having an uncertain diagnosis, we identified a combination of a monoallelic pathogenic variant in two genes in four patients with a CC. One patient (patient 30) presented with a PV in *ABCB4* and a LPV variant in *ATP8B1*. The two proteins encoded by these genes, MDR3 and FIC1 respectively, had been shown to work synergistically to maintain the lipid composition of the canalicular membranes of hepatocytes. Another patient (patient 48) presented with a PV in *MYO5B*, involved
in the intracellular traffic of proteins, and a CFTR-RD variant. We cannot exclude that the combination of these variants could explain CC in these patients. NGS sequencing revealed cystic fibrosis in three patients and detected heterozygous CFTR variants in 11 patients (ten CFTR-RD variants, one CF causing variant). Four of them had a CC phenotype and seven a TNC phenotype. For these variants, we referred to data available in CFTR specialized databases. Our work raises the hypothesis of a specific cholestasis entity associated with mutations in CFTR, different from chronic liver disease associated with CF. This clinical entity could be explained through the dysfunction of the CFTR protein (encoding a chloride channel), leading to disturbances in the hydro-electrolytic balance of the bile and, consequently, to an alteration of bile flow. This hypothesis has already been mentioned in animal models (ICP, in mice [40,41]). Interestingly, patients presenting neonatal cholestasis and a CFTR variant have been reported by two other teams [15,38]. Both genetic and phenotypic features of this subtype of cholestasis related to CFTR-RD variants appear to be different from the classical chronic liver diseases associated with CF. Regarding clinical features, data in the literature indicate that CF-associated liver disease usually becomes clinically patent in school-age children and manifests as cirrhosis and portal hypertension rather than as CC. We suggest that cholestasis associated with CFTR-RD variants described in this work should better be integrated into "CFTR spectrum pathologies not associated with cystic fibrosis", such as isolated vas deferens agenesis and diffuse bronchiectasis [9]. This hypothesis is supported by the fact that 6 of the above-mentioned variants are CFTR-RD variants (Table 4b). Further clinical and fundamental studies are needed to strengthen our hypothesis.

Our work also presents limitations that deserve to be detailed. First, regarding the NGS technique itself, detection of mutations involving large variations in CNV (large insertion or deletion mutations) requires a very careful analysis of the databases generated by sequencing, especially in poorly covered regions. qPCR, MLPA, or CGHarray may be required when CNV is suspected in front of abnormal deep coverage ratio (<0.5 or >1.6) in targeted regions. The biologist's viewpoint in selecting variants of interest in the second phase of bio analysis is very important. This selection is operator dependent: we considered it reproducible in our work because the same person had carried it out for all patients. In addition, mutations in non-coding sequences (deep intronic, or in promoter regions) can also influence the phenotype without being detected by the NGS technique. Moreover, since the number of genes involved in neonatal and infantile cholestasis is constantly evolving, the gene panel must be frequently updated with the addition of new genes. Among the 372 patients remaining undiagnosed in this study (61.8%), those with chronic cholestasis and progressive liver disease are candidates for further molecular analyses, such as last update versions of the gene panel, exome analysis (Whole Exome Sequencing—WES), and genome analysis (Whole Genome Sequencing—WGS). In the setting of a French initiative for genomic medicine aiming to change methods for diagnosis, prevention, and treatment of patients ("Plan France médecine génomique 2025"), these patients with severe liver diseases of pediatric onset can benefit from a WGS (trio analysis). However, it is to be reminded that these studies generate important bioinformatics analysis, and are more expensive and time-consuming than the analysis of data provided by a panel of targeted genes. Finally, as in all genetic diseases, environmental or epigenetic factors may be involved and must be investigated [42].

**5. Conclusions**

Herein we report the results of more than five years of experience with NGS using a targeted gene panel in the diagnostic work up of children presenting with cholestasis. Used as a first-line diagnostic tool, it provides an etiological diagnosis in 28.1% of patients and completes the routine etiological assessment. The hypothesis of multigenic contribution and imputability of VOUS in patients' phenotypes require further investigations, both in clinical data collection and functional aspects. Cellular and animal models could contribute to a better understanding of these data. Thus, for some patients, the interpretation of
molecular findings should be considered an ongoing process, which deserves to be updated
taking into account newly published data and the clinical evolution of the patients.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/diagnostics12051169/s1, Figure S1: Details and dates of use of the 4 genes panels, Table S1:Genotype of patients with certain diagnosis (ranged in alphabetic order, regarding the disease, Table S2: Main characteristics of 2 patients with PFIC3 and cholangiopathy mimicking neonatal sclerosing cholangitis, Figure S2. MRI-cholangiography in patient 188 with MDR3 deficiency (PFIC3) and transhepatic cholecystography in an additional PFIC3 patient not included in the study.

**Author Contributions:** Conceptualization, M.A., A.S., J.B., E.J. and E.G.; methodology, M.A., A.S., J.B., E.J. and E.G.; software, A.S. and J.B.; validation, M.A., A.S., E.J. and E.G.; formal analysis, M.A., A.S. and E.G.; investigation, M.A., A.S. and E.G.; writing—original draft preparation, M.A.; writing—review and editing, M.A., A.S. and E.G. All other authors furnished clinical data for the majority of the patients included in the study. All doctors cited in the acknowledgments section furnished clinical data for the rest of the patients. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** This study complies with the General Data Protection Regulation (RGPD—Regulation (EU) No. 2016/679) and the French Data Protection Act. The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Bicêtre Hospital (declared and listed in the Data Processing Registry of the Assistance Publique-Hôpitaux de Paris, under protocol number: 20191112110112; date of approval: 12 November 2019" The child's referring physician, obtained written, informed consent of the patients' parents, in accordance with the local guidelines established by the hospital's ethics committee.

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Acknowledgments:** The authors thank all referring physicians who provided clinical data of their patients: M. Al Khouri (Montpellier), J.B. Arnoux (Necker), D. Astruc (Strasbourg), G. Barjonet (Montélimar), A. Béranger (Necker), M. Bonneton (Nancy), L. Bridoux-Henno (Rennes), M. Butori (Nice), P. Callamand (Béziers), A. Comte (Besançon), C. Cremilleux (Saint Etienne), G. Dassieu (Créteil), D. Debray (Paris), Z. Demir (Paris), S. Destombe (Saint Etienne), B. Digeon (Reims), G. Dimitrov (Orléans), C. Dubuisson (Paris), C. Dumant Forest (Rouen), I. Etienne (Mayotte), L. Restier (Lyon), R. Maudinas (Dijon), K. Mention (Lille), A. Turquet (La Réunion), C. Rivet (Lyon), B. Rohmer (Lyon), J. Rebeuh, (Strasbourg), C. Wicker (Strasbourg). They also thank A. Remil and J. Honey for their help in the redaction of this manuscript. We thank AMFE (Association Maladie Foie Enfants, Malakof, France), MLD (Monaco Liver Disorder, Monaco), Association "Pour Louis 1000 Foie Merci" (Fournet-Luisans, France), and Fondation Rumsey-Cartier (Genève, Switzerland) for their support.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Harpavat, S.; Garcia-Prats, J.A.; Anaya, C.; Brandt, M.L.; Lupo, P.J.; Finegold, M.J.; Obuobi, A.; ElHennawy, A.A.; Jarriel, W.S.; Shneider, B.L. Diagnostic Yield of Newborn Screening for Biliary Atresia Using Direct or Conjugated Bilirubin Measurements. *JAMA* **2020**, *323*, 1141–1150. [CrossRef]

2. Fawaz, R.; Baumann, U.; Ekong, U.; Fischler, B.; Hadzic, N.; Mack, C.L.; McLin, V.A.; Molleston, J.P.; Neimark, E.; Ng, V.L.; et al. Guideline for the Evaluation of Cholestatic Jaundice in Infants: Joint Recommendations of the North American Society for Pediatric Gastroenterology, Hepatology, and Nutrition and the European Society for Pediatric Gastroenterology, Hepatology, and Nutrition. *J. Pediatr. Gastroenterol. Nutr.* **2017**, *64*, 154–168.

3. Nicastro, E.; Di Giorgio, A.; Marchetti, D.; Barboni, C.; Cereda, A.; Iascone, M.; D'Antiga, L. Diagnostic Yield of an Algorithm for Neonatal and Infantile Cholestasis Integrating Next-Generation Sequencing. *J. Pediatr.* **2019**, *211*, 54–62.e4. [CrossRef]

4. Bahassi, E.M.; Stambrook, P.J. Next-generation sequencing technologies: Breaking the sound barrier of human genetics. *Mutagenesis* **2014**, *29*, 303–310. [CrossRef]

5. Herbst, S.M.; Schirmer, S.; Posovszky, C.; Jochum, F.; Rödl, T.; Schroeder, J.A.; Barth, T.F.; Hehr, U.; Melter, M.; Vermehren, J. Taking the next step forward-Diagnosing inherited infantile cholestatic disorders with next generation sequencing. *Mol. Cell. Probes* **2015**, *29*, 291–298. [CrossRef]

6. Nicastro, E.; D'Antiga, L. Next generation sequencing in pediatric hepatology and liver transplantation. *Liver Transplant.* **2018**, *24*, 282–293. [CrossRef]
7. Chappell, K.; Francou, B.; Habib, C.; Huby, T.; Leoni, M.; Cottin, A.; Nadal, F.; Adnet, E.; Paoli, E.; Oliveira, C.; et al. Galaxy Is a Suitable Bioinformatics Platform for the Molecular Diagnosis of Human Genetic Disorders Using High-Throughput Sequencing Data Analysis: Five Years of Experience in a Clinical Laboratory. *Clin. Chem.* **2022**, *68*, 313–321. [CrossRef]

8. Richards, S.; Aziz, N.; Bale, S.; Bick, D.; Das, S.; Gastier-Foster, J.; Grody, W.W.; Hegde, M.; Lyon, E.; Spector, E.; et al. 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. [CrossRef]

9. Castellani, C.; Cuppens, H.; Macek, M.; Cassiman, J.J.; Kerem, E.; Durie, P.; Tullis, E.; Assael, B.M.; Bombieri, C.; Brown, A.; et al. Consensus on the use and interpretation of cystic fibrosis mutation analysis in clinical practice. *J. Cyst. Fibros.* **2008**, *7*, 179–196. [CrossRef]

10. Ziol, M.; Barbu, V.; Rosmorduc, O.; Frassati-Biaggi, A.; Barget, N.; Hermelin, B.; Scheffer, G.L.; Bennouna, S.; Trinchet, J.-C.; Beaugrand, M.; et al. ABCB4 heterozygous gene mutations associated with fibrosing cholestatic liver disease in adults. *Gastroenterology* **2008**, *135*, 131–141. [CrossRef]

11. Dong, C.; Condat, B.; Picon-Coste, M.; Chrétien, Y.; Potier, P.; Noblinski, B.; Arrivé, L.; Hauuy, M.-P.; Barbu, V.; Maftouh, A.; et al. Low-phospholipid-associated cholelithiasis syndrome: Prevalence, clinical features, and comorbidities. *JHEP Rep. Innov. Hepatol.* **2021**, *3*, 100201. [CrossRef]

12. Roelandt, P.; Antoniou, A.; Libbrecht, L.; Van Steenbergen, W.; Laleman, W.; Verslype, C.; Van der Merwe, S.; Nevens, F.; De Vos, R.; Fischer, E.; et al. HNF1B deficiency causes ciliary defects in human cholangiocytes. *Hepatology* **2012**, *56*, 1178–1181. [CrossRef]

13. Yi, F.; Poskanzer, S.A.; Myers, C.T.; Thies, J.; Collins, C.J.; Dayuha, R.; Duong, P.; Houwen, R.; Hahn, S.H. p.P1379S, a benign variant with reduced ATP7B protein level in Wilson Disease. *JIMD Rep.* **2020**, *54*, 32–36. [CrossRef]

14. Khabou, B.; Durand-Schneider, A.-M.; Delaunay, J.-L.; Aït-Slimane, T.; Barbu, V.; Fakhfakh, F.; Housset, C.; Maurice, M. Comparison of in silico prediction and experimental assessment of ABCB4 variants identified in patients with biliary diseases. *Int. J. Biochem. Cell Biol.* **2017**, *89*, 101–109. [CrossRef]

15. Karpen, S.J.; Kamath, B.M.; Alexander, J.J.; Ichetovkin, I.; Rosenthal, P.; Sokol, R.J.; Dunn, S.; Thompson, R.J.; Heubi, J.E. Use of a Comprehensive 66-Gene Cholestasis Sequencing Panel in 2171 Cholestatic Infants, Children, and Young Adults. *J. Pediatr. Gastroenterol. Nutr.* **2021**, *72*, 654–660. [CrossRef]

16. Matte, U.; Mourya, R.; Miethke, A.; Liu, C.; Kauffmann, G.; Moyer, K.; Zhang, K.; Bezerra, J.A. Analysis of gene mutations in children with cholestasis of undefined etiology. *J. Pediatr. Gastroenterol. Nutr.* **2010**, *51*, 488–493. [CrossRef]

17. Wang, N.-L.; Lu, Y.-L.; Zhang, P.; Zhang, M.-H.; Gong, J.-Y.; Lu, Y.; Xie, X.-B.; Qiu, Y.-L.; Yan, Y.-Y.; Wu, B.-B.; et al. A Specially Designed Multi-Gene Panel Facilitates Genetic Diagnosis in Children with Intrahepatic Cholestasis: Simultaneous Test of Known Large Insertions/Deletions. *PLoS ONE* **2016**, *11*, e0164058. [CrossRef]

18. Togawa, T.; Sugiura, T.; Ito, K.; Endo, T.; Aoyama, K.; Ohashi, K.; Negishi, Y.; Kudo, T.; Ito, R.; Kikuchi, A.; et al. Molecular Genetic Dissection and Neonatal/Infantile Intrahepatic Cholestasis Using Targeted Next-Generation Sequencing. *J. Pediatr.* **2016**, *171*, 171–177.e4. [CrossRef]

19. Stalke, A.; Skawran, B.; Auber, B.; Illig, T.; Schlegelberger, B.; Junge, N.; Goldschmidt, I.; Leiskau, C.; von Neuhoff, N.; Baumann, U.; et al. Diagnosis of monogenic liver diseases in childhood by next-generation sequencing. *Clin. Genet.* **2018**, *93*, 665–670. [CrossRef]

20. Chen, H.-L.; Li, H.-Y.; Wu, J.-F.; Wu, S.-H.; Chen, H.-L.; Yang, Y.-H.; Hsu, Y.-H.; Liou, B.-Y.; Chang, M.-H.; Ni, Y.-H. Panel-Based Next-Generation Sequencing for the Diagnosis of Cholestatic Genetic Liver Diseases: Clinical Utility and Challenges. *J. Pediatr.* **2019**, *205*, 153–159.e6. [CrossRef]

21. Shagrani, M.; Burkholder, J.; Broering, D.; Abouelhoda, M.; Faquih, T.; El-Kalioby, M.; Subhani, S.N.; Goljan, E.; Albar, R.; Monies, D.; et al. Genetic profiling of children with advanced cholestatic liver disease. *Clin. Genet.* **2017**, *92*, 52–61. [CrossRef]

22. Lipiński, P.; Ciara, E.; Jurkiewicz, D.; Pollak, A.; Wypchło, M.; Płoski, R.; Cielecka-Kuszyk, J.; Socha, P.; Pawłowska, J.; Jankowska, I. Targeted Next-Generation Sequencing in Diagnostic Approach to Monogenic Cholestatic Liver Disorders-Single-Center Experience. *Front. Pediatr.* **2020**, *8*, 414. [CrossRef]

23. Fickert, P.; Fuchsbichler, A.; Wagner, M.; Zollner, G.; Kaser, A.; Tilg, H.; Krause, R.; Lammert, F.; Langner, C.; Zatloukal, K.; et al. Regurgitation of bile acids from leaky bile ducts causes sclerosing cholangitis in Mdr2 (Abcb4) knockout mice. *Gastroenterology* **2004**, *127*, 261–274. [CrossRef]

24. Fickert, P.; Zollner, G.; Fuchsbichler, A.; Stumptner, C.; Weiglein, A.H.; Lammert, F.; Marschall, H.-U.; Tsybrovskyy, O.; Zatloukal, K.; Denk, H.; et al. Ursodeoxycholic acid aggravates bile infarcts in bile duct-ligated and Mdr2 knockout mice via disruption of cholangioles. *Gastroenterology* **2002**, *123*, 1238–1251. [CrossRef]

25. Jacquemin, E.; De Vree, J.M.; Cresteil, D.; Sokal, E.M.; Sturm, E.; Dumont, M.; Scheffer, G.L.; Paul, M.; Burdelski, M.; Bosma, P.J.; et al. The wide spectrum of multidrug resistance 3 deficiency: From neonatal cholestasis to cirrhosis of adulthood. *Gastroenterology* **2001**, *120*, 1448–1458. [CrossRef]

26. Schatz, S.B.; Jüngst, C.; Keitel-Anselmo, V.; Kubitz, R.; Becker, C.; Gerner, P.; Pfister, E.-D.; Goldschmidt, I.; Junge, N.; Wenning, D.; et al. Phenotypic spectrum and diagnostic pitfalls of ABCB4 deficiency depending on age of onset. *Hepatol. Commun.* **2018**, *2*, 504–514. [CrossRef]

27. Wendum, D.; Barbu, V.; Rosmorduc, O.; Arrivé, L.; Fléjou, J.-F.; Poupon, R. Aspects of liver pathology in adult patients with MDR3/ABCB4 gene mutations. *Virchows Arch.* **2012**, *460*, 291–298. [CrossRef]
28. Benzimra, J.; Derhy, S.; Rosmorduc, O.; Menu, Y.; Poupon, R.; Arrivé, L. Hepatobiliary anomalies associated with ABCB4/MDR3 deficiency in adults: A pictorial essay. *Insights Imaging* **2013**, *4*, 331–338. [CrossRef]

29. Qiu, Y.-L.; Liu, T.; Abuduxikuer, K.; Hao, C.-Z.; Gong, J.-Y.; Zhang, M.-H.; Li, L.-T.; Yan, Y.-Y.; Li, J.-Q.; Wang, J.-S. Novel missense mutation in VPS33B is associated with isolated low gamma-glutamyltransferase cholestasis: Attenuated, incomplete phenotype of arthrogryposis, renal dysfunction, and cholestasis syndrome. *Hum. Mutat.* **2019**, *40*, 2247–2257. [CrossRef]

30. Aamann, L.; Ørntoft, N.; Vogel, I.; Grønbaek, H.; Becher, N.; Vilstrup, H.; Ott, P.; Lildballe, D.L. Unexplained cholestasis in adults and adolescents: Diagnostic benefit of genetic examination. *Scand. J. Gastroenterol.* **2018**, *53*, 305–311. [CrossRef]

31. Hertel, P.M.; Bull, L.N.; Thompson, R.J.; Goodrich, N.P.; Ye, W.; Magee, J.C.; Squires, R.H.; Bass, L.M.; Heubi, J.E.; Kim, G.E.; et al. Mutation Analysis and Disease Features at Presentation in a Multi-Center Cohort of Children With Monogenic Cholestasis. *J. Pediatr. Gastroenterol. Nutr.* **2021**, *73*, 169–177. [CrossRef]

32. Zhou, J.; Yao, Y.; Huang, X.; Chen, D.; Zhang, T.; Zhang, Y. Association between bile salt export pump polymorphisms and intrahepatic cholestasis of pregnancy susceptibility: A meta-analysis of case-control studies. *Gynecol. Endocrinol. Nutr.* **2019**, *35*, 179–183. [CrossRef]

33. Jacquemin, E.; Lykavieris, P.; Chaoui, N.; Hadchouel, M.; Bernard, O. Transient neonatal cholestasis: Origin and outcome. *J. Pediatr.* **1998**, *133*, 563–567. [CrossRef]

34. Champion, V.; Carbajal, R.; Lozar, J.; Girard, I.; Mitanchez, D. Risk factors for developing transient neonatal cholestasis. *J. Pediatr. Gastroenterol. Nutr.* **2012**, *55*, 592–598. [CrossRef]

35. Liu, L.-Y.; Wang, X.-H.; Lu, Y.; Zhu, Q.-R.; Wang, J.-S. Association of variants of ABCB11 with transient neonatal cholestasis. *Pediatr. Int.* **2013**, *55*, 138–144. [CrossRef]

36. Hermeziu, B.; Sanlaville, D.; Girard, M.; Léonard, C.; Lyonnet, S.; Jacquemin, E. Heterozygous bile salt export pump deficiency: A possible genetic predisposition to transient neonatal cholestasis. *J. Pediatr. Gastroenterol. Nutr.* **2006**, *42*, 114–116. [CrossRef]

37. Jacquemin, E.; Malan, V.; Rio, M.; Davit-Spraul, A.; Cohen, J.; Landrieu, P.; Bernard, O. Heterozygous FIC1 deficiency: A new genetic predisposition to transient neonatal cholestasis. *J. Pediatr. Gastroenterol. Nutr.* **2010**, *50*, 447–449. [CrossRef]

38. Santos Silva, E.; Almeida, A.; Frutuoso, S.; Martins, E.; Valente, M.J.; Santos-Silva, A.; Lopes, A.I. Neonatal Cholestasis Over Time: Changes in Epidemiology and Outcome in a Cohort of 154 Patients From a Portuguese Tertiary Center. *Front. Pediatr.* **2020**, *8*, 351. [CrossRef]

39. Jung, C.; Driancourt, C.; Baussan, C.; Zater, M.; Hadchouel, M.; Meunier-Rotival, M.; Guiochon-Mantel, A.; Jacquemin, E. Prenatal molecular diagnosis of inherited cholestatic diseases. *J. Pediatr. Gastroenterol. Nutr.* **2007**, *44*, 453–458. [CrossRef]

40. Aleksandrova, M.I.; Kushnareva, N.S.; Smirnova, O.V. Manifestation of cystic fibrosis transmembrane regulator (CFTR) in hepatic ductal structures and renal tubules of female rats with experimental cholestasis of pregnancy. *Bull. Exp. Biol. Med.* **2014**, *156*, 654–657. [CrossRef]

41. Miller, A.C.; Comellas, A.P.; Hornick, D.B.; Stoltz, D.A.; Cavanaugh, J.E.; Gerke, A.K.; Welsh, M.J.; Zabner, J.; Polgreen, P.M. Cystic fibrosis carriers are at increased risk for a wide range of cystic fibrosis-related conditions. *Proc. Natl. Acad. Sci. USA* **2020**, *117*, 1621–1627. [CrossRef]

42. Mann, D.A. Epigenetics in liver disease. *Hepatol. Baltim. Md* **2014**, *60*, 1418–1425. [CrossRef]