Background & aims
The etiology of cholestasis remains unknown in many children. We surveyed the genome of children with chronic cholestasis for variants in genes not previously associated with liver disease and validated their biological relevance in zebrafish and murine models.

Method
Whole-exome (n = 4) and candidate gene sequencing (n = 89) was completed on 93 children with cholestasis and normal serum γ-glutamyl transferase (GGT) levels without pathogenic variants in genes known to cause low GGT cholestasis such as ABCB11 or ATP8B1. CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 genome editing was used to induce frameshift pathogenic variants in the candidate gene in zebrafish and mice.

Results
In a 1-year-old female patient with normal GGT cholestasis and bile duct paucity, we identified a homozygous truncating pathogenic variant (c.198delA, p.Gly67Alafs∗6) in the ABCC12 gene (NM_033226). Five additional rare ABCC12 variants, including a pathogenic one, were detected in our cohort. ABCC12 encodes multidrug resistance-associated protein 9 (MRP9) that belongs to the adenosine 5′-triphosphate-binding cassette transporter C family with unknown function and no previous implication in liver disease. Immunohistochemistry and Western blotting revealed conserved MRP9 protein expression in the bile ducts in human, mouse, and zebrafish. Zebrafish abcc12-null mutants were prone to cholangiocyte apoptosis, which caused progressive bile duct loss during the juvenile stage. MRP9-deficient mice had fewer well-formed interlobular bile ducts and higher serum alkaline phosphatase levels compared with wild-type mice. They exhibited aggravated cholangiocyte apoptosis, hyperbilirubinemia, and liver fibrosis upon cholic acid challenge.