Replication of a GWAS signal in a Caucasian population implicates ADD3 in susceptibility to biliary atresia

Genetics of biliary atresia: Approaches, pathological insights and challenges

Biliary atresia (BA) is a severe neonatal cholestatic disorder marked by fibro-obliteration of the extrahepatic and intrahepatic bile ducts. It is the most common cause of pediatric end-stage liver disease and the leading indication for liver transplantation in children. There is significant heterogeneity in the etiology, involving various genetic and environmental factors such as viral infection, immune dysregulation and genetic predisposition to defective hepatobiliary development. In this review, we discuss the strategies to uncover the genetic factors underlying BA and highlight their associated molecular and pathological mechanisms, as well as the challenges faced in this area of research.

A narrative review of genes associated with liver fibrosis in biliary atresia

Background and Objective

Biliary atresia (BA) is characterized by biliary inflammation and obstruction. In the later phase, liver fibrosis occurs. Although the etiology of BA is believed to be multi-factorial, genetic predisposition has been proposed to play a critical role in the pathogenesis. This review aimed to provide an updated summary of the genes that have been reported to be involved in BA-associated liver fibrosis.

Methods

The review was conducted via evaluation of MalaCards (BA disease: MalaCards-research articles, drugs, genes, clinical trials) which is a universally applied website including various human disease database. The database of genes that are involved in liver fibrosis were studied.

Key Content and Findings

Thirty-one genes that are associated with BA according to the disease relevance score were reviewed after further evaluations. Eleven genes (GPT, NR1H4, TGF-B1, MMP7, CCN2, TIMP1, SPP1, ADD3, KRT7, ADD3-AS1, SOX9) that are specific and with a potential association with liver fibrosis were selected for detailed description. Increased expression of GPT, TGF-B1, MMP7, CCN2, TIMP1, SPP1, ADD3, KRT7 and ADD3-AS1 maybe associated with the development of liver fibrosis in BA patients, while the expression of NR1H4 and SOX9 are more likely to suppress liver fibrosis.

Conclusions

Current scientific evidence using gene database has revealed a close association between genetic anomalies and the pathogenesis of liver fibrosis in BA. With a better understanding of these anomalies, therapy targeting these related genes may be a new therapeutic approach to alleviate liver fibrosis in BA.

Biliary atresia

Biliary atresia (BA) is a progressive inflammatory fibrosclerosing disease of the biliary system and a major cause of neonatal cholestasis. It affects 1:5,000-20,000 live births, with the highest incidence in Asia. The pathogenesis is still unknown, but emerging research suggests a role for ciliary dysfunction, redox stress and hypoxia. The study of the underlying mechanisms can be conceptualized along the likely prenatal timing of an initial insult and the distinction between the injury and prenatal and postnatal responses to injury. Although still speculative, these emerging concepts, new diagnostic tools and early diagnosis might enable neoadjuvant therapy (possibly aimed at oxidative stress) before a Kasai portoenterostomy (KPE). This is particularly important, as timely KPE restores bile flow in only 50-75% of patients of whom many subsequently develop cholangitis, portal hypertension and progressive fibrosis; 60-75% of patients require liver transplantation by the age of 18 years. Early diagnosis, multidisciplinary management, centralization of surgery and optimized interventions for complications after KPE lead to better survival. Postoperative corticosteroid use has shown benefits, whereas the role of other adjuvant therapies remains to be evaluated. Continued research to better understand disease mechanisms is necessary to develop innovative treatments, including adjuvant therapies targeting the immune response, regenerative medicine approaches and new clinical tests to improve patient outcomes.

Progress in Biomarkers Related to Biliary Atresia

Biliary atresia (BA) is a congenital cholestatic disease that can seriously damage children's liver function. It is one of the main reasons for liver transplantation in children. Early diagnosis of BA is crucial to the prognosis of patients, but there is still a lack of reliable non-invasive diagnostic methods. Additionally, as some children are in urgent need of liver transplantation, evaluating the stage of liver fibrosis and postoperative native liver survival in children with BA using a straightforward, efficient, and less traumatic method is a major focus of doctors. In recent years, an increasing number of BA-related biomarkers have been identified and have shown great potential in the following three aspects of clinical practice: diagnosis, evaluation of the stage of liver fibrosis, and prediction of native liver survival. This review focuses on the pathophysiological function and clinical application of three novel BA-related biomarkers, namely MMP-7, FGF-19, and M2BPGi. Furthermore, progress in well-known biomarkers of BA such as gamma-glutamyltransferase, circulating cytokines, and other potential biomarkers is discussed, aiming to provide a reference for clinical practice.

Biliary atresia is associated with polygenic susceptibility in ciliogenesis and planar polarity effector genes

BACKGROUND & AIMS

Biliary atresia (BA) is poorly understood and leads to liver transplantation (LT), with the requirement for and associated risks of lifelong immunosuppression, in most children. We performed a genome-wide association study (GWAS) to determine the genetic basis of BA.

METHODS

We performed a GWAS in 811 European BA cases treated with LT in US, Canadian and UK centers, and 4,654 genetically matched controls. Whole-genome sequencing of 100 cases evaluated synthetic association with rare variants. Functional studies included whole liver transcriptome analysis of 64 BA cases and perturbations in experimental models.

RESULTS

A GWAS of common single nucleotide polymorphisms (SNPs), i.e. allele frequencies >1%, identified intronic SNPs rs6446628 in AFAP1 with genome-wide significance (p = 3.93E-8) and rs34599046 in TUSC3 at sub-threshold genome-wide significance (p = 1.34E-7), both supported by credible peaks of neighboring SNPs. Like other previously reported BA-associated genes, AFAP1 and TUSC3 are ciliogenesis and planar polarity effectors (CPLANE). In gene-set-based GWAS, BA was associated with 6,005 SNPs in 102 CPLANE genes (p = 5.84E-15). Compared with non-CPLANE genes, more CPLANE genes harbored rare variants (allele frequency <1%) that were assigned Human Phenotype Ontology terms related to hepatobiliary anomalies by predictive algorithms, 87% vs. 40%, p <0.0001. Rare variants were present in multiple genes distinct from those with BA-associated common variants in most BA cases. AFAP1 and TUSC3 knockdown blocked ciliogenesis in mouse tracheal cells. Inhibition of ciliogenesis caused biliary dysgenesis in zebrafish. AFAP1 and TUSC3 were expressed in fetal liver organoids, as well as fetal and BA livers, but not in normal or disease-control livers. Integrative analysis of BA-associated variants and liver transcripts revealed abnormal vasculogenesis and epithelial tube formation, explaining portal vein anomalies that co-exist with BA.

CONCLUSIONS

BA is associated with polygenic susceptibility in CPLANE genes. Rare variants contribute to polygenic risk in vulnerable pathways via unique genes.

IMPACT AND IMPLICATIONS

Liver transplantation is needed to cure most children born with biliary atresia, a poorly understood rare disease. Transplant immunosuppression increases the likelihood of life-threatening infections and cancers. To improve care by preventing this disease and its progression to transplantation, we examined its genetic basis. We find that this disease is associated with both common and rare mutations in highly specialized genes which maintain normal communication and movement of cells, and their organization into bile ducts and blood vessels during early development of the human embryo. Because defects in these genes also cause other birth defects, our findings could lead to preventive strategies to lower the incidence of biliary atresia and potentially other birth defects.

Contribution of ADD3 and the HLA Genes to Biliary Atresia Risk in Chinese

Nonsyndromic biliary atresia (BA) is a rare polygenic disease, with autoimmunity, virus infection and inflammation thought to play roles in its pathogenesis. We conducted a genome-wide association study in 336 nonsyndromic BA infants and 8900 controls. Our results validated the association of rs17095355 in ADD3 with BA risk (odds ratio (OR) = 1.70, 95% confidence interval (95% CI) = 1.49-1.99; p = 4.07 × 10-11). An eQTL analysis revealed that the risk allele of rs17095355 was associated with increased expression of ADD3. Single-cell RNA-sequencing data and immunofluorescence analysis revealed that ADD3 was moderately expressed in cholangiocytes and weakly expressed in hepatocytes. Immuno-fluorescent staining showed abnormal deposition of ADD3 in the cytoplasm of BA hepatocytes. No ADD3 auto-antibody was observed in the plasma of BA infants. In the HLA gene region, no variants achieved genome-wide significance. HLA-DQB1 residue Ala57 is the most significant residue in the MHC region (OR = 1.44, 95% CI = 1.20-1.74; p = 1.23 × 10-4), and HLA-DQB1 was aberrantly expressed in the bile duct cells. GWAS stratified by cytomegalovirus (CMV) IgM status in 87 CMV IgM (+) BA cases versus 141 CMV IgM (-) BA cases did not yield genome-wide significant associations. These findings support the notion that common variants of ADD3 account for BA risk. The HLA genes might have a minimal role in the genetic predisposition of BA due to the weak association signal. CMV IgM (+) BA patients might not have different genetic risk factor profiles compared to CMV IgM (-) subtype.

Genetic Contributions to Biliary Atresia: A Developmental Cholangiopathy

Biliary atresia (BA) is the most prevalent serious liver disease of infancy and childhood, and the principal indication for liver transplantation in pediatrics. BA is best considered as an idiopathic panbiliary cholangiopathy characterized by obstruction of bile flow and consequent cholestasis presenting during fetal and perinatal periods. While several etiologies have been proposed, each has significant drawbacks that have limited understanding of disease progression and the development of effective treatments. Recently, modern genetic analyses have uncovered gene variants contributing to BA, thereby shifting the paradigm for explaining the BA phenotype from an acquired etiology (e.g., virus, toxin) to one that results from genetically altered cholangiocyte development and function. Herein we review recently reported genetic contributions to BA, highlighting the enhanced representation of variants in biological pathways involving ciliary function, cytoskeletal structure, and inflammation. Finally, we blend these findings as a new framework for understanding the resultant BA phenotype as a developmental cholangiopathy.

Genetic alterations and molecular mechanisms underlying hereditary intrahepatic cholestasis

Hereditary cholestatic liver disease caused by a class of autosomal gene mutations results in jaundice, which involves the abnormality of the synthesis, secretion, and other disorders of bile acids metabolism. Due to the existence of a variety of gene mutations, the clinical manifestations of children are also diverse. There is no unified standard for diagnosis and single detection method, which seriously hinders the development of clinical treatment. Therefore, the mutated genes of hereditary intrahepatic cholestasis were systematically described in this review.

The correlation between rs2501577 gene polymorphism and biliary atresia: a systematic review and meta-analysis

IMPORTANCE

Multiple studies indicate a possible correlation between ADD3 rs2501577 and biliary atresia susceptibility; however, a conclusive determination has yet to be made.

OBJECTIVE

Investigate the role of ADD3 rs2501577 in biliary atresia susceptibility across diverse populations.

DATA SOURCES

The study protocol has been registered on PROSPERO, an international platform for systematic review registration (PROSPERO ID: CRD42023384641). The following databases will be searched until February 1, 2023: PubMed, Embase, Cochrane, CBM, Web of Science, and CNKI.

STUDY SELECTION

Eight studies were selected from seven papers to assess the data. A total of 7651 participants were included, consisting of 1662 in the BA group and 5989 in the NC group.

DATA EXTRACTION AND SYNTHESIS

Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines were followed while conducting the systematic reviews and meta-analyses. Two authors independently assessed the quality of the included studies using the Newcastle-Ottawa Quality Assessment Scale. The significance of the pooled odds ratio (OR) was evaluated with a Z test, and statistical heterogeneity across studies was assessed using the I2 and Q statistics. Publication bias was assessed using Egger's and Begg's tests.

MAIN OUTCOME(S) AND MEASURE(S)

The primary study outcome was the development of biliary atresia. Subgroup analysis was performed based on race, region, and assessment of Hardy-Weinberg equilibrium (HWE).

RESULTS

The studies indicate that the ADD3 rs2501577 susceptibility locus increases the risk of developing biliary atresia, regardless of allelic, homozygote, dominant, and recessive gene inheritance models. Furthermore, ADD3 has been found to be associated with apoptosis, cell cycle, and cell damage repair based on functional analysis.

CONCLUSIONS AND RELEVANCE

The ADD3 rs2501577 polymorphic locus is associated with an increased risk of biliary atresia, particularly in Asian populations. This study recommends further investigation of the ADD3 rs2501577 locus in Asian populations to validate its role in the diagnosis of biliary atresia.

A novel pathogenesis concept of biliary atresia approached by combined molecular strategies

Cholestatic jaundice is one of the most common neonatal conditions. BA, a correctable cholangiopathy, presents with cholestatic jaundice within the first weeks of life. The inflammation of bile ducts leads to progressive fibrosclerosis involving biliary trees, followed by cirrhosis and liver failure. With the use of modern molecular studies, this research aimed to define a novel pathogenesis by exploring variations. We performed genetic discovery by using supervised and unsupervised approaches. Ultimately, a combination of genetic variations and survival data was analyzed to strengthen the novel concept in this study. In this study, coding regions were explored to identify rare deleterious variants within genes from the first analysis together with gene sets reported in PFIC, and diseases with hyperbilirubinemia. Our unsupervised prioritization was primarily designed to identify novel causal genes from nonsynonymous variants derived by three biostatistical algorithms: enrichment analysis, burden test, and trio study. Survival analysis was integratively evaluated with a combination of identified causal genes. The individuals with identified variants from the supervised approach were frequently related to the severity of cirrhosis and poor postoperative outcome. In the unsupervised approach, nonsynonymous variants were enriched. Cilium and muscle related pathways had a significant correlation. CCDC8 was statistically significant gene in which six cases carried mutations identified through burden analysis. Individuals who carried variants in corresponding genes and significant pathways had significantly lower native-liver survival than individuals in whom none of these variants were identified (log-rank p value 0.016). This study explored genetic variations by multiple strategies. Different pathways of cholestatic diseases have been found to be associated with BA. Therefore, BA may be characterized as a shared sequela of many cholestatic disorders. Susceptibility in those pathways suggested an association with BA and strengthened this proposed novel hypothesis. The results emphasized the consequences of many disruptive pathophysiologies.

Impact of EFEMP1 on the survival outcome of biliary atresia in Thai infants

Genome-wide association studies (GWASs) have identified a genetic associated between EFEMP1 and biliary atresia (BA). To examine the susceptibility of single nucleotide polymorphisms (SNPs) in EFEMP1 in Thai BA patients, we performed an analysis of the genetic associations and biological interactions with previously reported key SNPs in ADD3, a key gene associated with BA. The study also used high-throughput sequencing to detect novel variants in both genes. In addition, the clinical impact of EFEMP1 SNPs in terms of survival association was also evaluated. The genotypes of 60 BA patients and 179 controls were evaluated using a TaqMan genotyping assay for rs2501577 and rs17095355 in ADD3 and rs6761893 and rs727878 in EFEMP1. The genotype frequencies were analyzed together with the SNP-SNP interactions. Fine mapping by whole-exome sequencing was performed to identify deleterious variants within both genes, and the survival analysis results were analyzed with the EFEMP1 SNPs. The recessive genotypes of rs2501577, rs17095355 and rs6761893 showed significantly higher frequencies in the BA patients than the controls, and a logistic regression showed that minor alleles of those SNPs increased the BA risk by ORs of 1.86, 1.67, and 1.84, respectively. Moreover, the SNP-SNP interference suggested that a combination of recessive alleles from the 2 genes resulted in an additive risk to BA. In addition, rare missense variants in the gene coding sequences were identified in 7 cases. Immunohistochemical studies revealed a pattern of ADD3 downregulation and EFEMP1 overexpression in the bile ducts of BA patients. Patients with the AA genotype of rs6761893 had significantly lower 5-year native liver survival (34.0%) than those with AT/TT (75.0%), with a log-rank p value of 0.041. Variants in EFEMP1 are associated with the occurrence of BA in Thai patients. In addition, these variants have an additive influence on BA risk when combined with ADD3 variants. Moreover, rs6761893 in EFEMP1 was indicative of survival in Thai BA patients.

Current Understanding in the Clinical Characteristics and Molecular Mechanisms in Different Subtypes of Biliary Atresia

Biliary atresia is a severe obliterative cholangiopathy in early infancy that is by far the most common cause of surgical jaundice and the most common indicator for liver transplantation in children. With the advanced knowledge gained from different clinical trials and the development of research models, a more precise clinical classification of BA (i.e., isolated BA (IBA), cystic BA (CBA), syndromic BA (SBA), and cytomegalovirus-associated BA (CMVBA)) is proposed. Different BA subtypes have similar yet distinguishable clinical manifestations. The clinical and etiological heterogeneity leads to dramatically different prognoses; hence, treatment needs to be specific. In this study, we reviewed the clinical characteristics of different BA subtypes and revealed the molecular mechanisms of their developmental contributors. We aimed to highlight the differences among these various subtypes of BA which ultimately contribute to the development of a specific management protocol for each subtype.

Role of long non-coding RNA-adducin 3 antisense RNA1 in liver fibrosis of biliary atresia

Biliary atresia (BA) is a devastating liver disease in neonates. Liver fibrosis is regarded as a universal and prominent feature of BA. Studies have revealed that long non-coding RNAs (lncRNAs) regulate cellular processes during the development of liver fibrosis in BA. Long non-coding RNA-adducin 3 antisense RNA1 (lnc-ADD3-AS1) has been shown to increase susceptibility to BA. However, the role of lnc-ADD3-AS1 in liver fibrosis in BA remains unclear. Here, we investigated the role of lnc-ADD3-AS1 in the proliferation, migration, and apoptosis of the immortalized human hepatic stellate cell (HSC) line, LX-2. We successfully overexpressed and silenced lnc-ADD3-AS1 in LX-2 cells using adenovirus vectors and evaluated the proliferation of transfected cells using the Cell Counting Kit-8 (CCK8) assay. Cell apoptosis was detected using annexin V-fluorescein isothiocyanate (FITC)/propidium iodide (PI) double staining and flow cytometry. We then analyzed cell migration by performing wound-scratch and transwell migration assays. Our results show that lnc-ADD3-AS1 significantly promoted LX-2 cell proliferation and attenuated apoptosis. More importantly, lncRNA-ADD3-AS1 significantly accelerated the migration of LX-2 cells. Our data indicated that lncRNA-ADD3-AS1 plays a role in the pathogenesis of liver fibrosis in patients with BA and may serve as a potential diagnostic marker for monitoring liver fibrosis in BA or as a therapeutic target for the disease.

Genetic Factors and Their Role in the Pathogenesis of Biliary Atresia

Biliary Atresia, a common basis for neonatal cholestasis and primary indication for Liver Transplantation, accounts for 60% of pediatric Liver Transplantations. While the pathogenesis of Biliary Atresia remains obscure, abnormalities within bile ducts and the liver, inflammation, fibrosis and cilia defects are thought to comprise the pathological basis for this condition. The findings of genetic variants in Biliary Atresia, such as Copy Number Variations and Single Nucleotide Polymorphism, are considered as essential factors in the development of this condition. In this review, we summarize and analyze these Biliary Atresia variants from a perspective of their pathological characteristics. In conclusion, such analyses may offer novel insights into the pathogenesis of Biliary Atresia and provide a foundation for future studies directed toward a better understanding and treatment of Biliary Atresia.

A joint transcriptional regulatory network and protein activity inference analysis identifies clinically associated master regulators for biliary atresia

Biliary atresia (BA) is a devastating cholangiopathy in neonate. Transcription factors (TFs), a type of master regulators in biological processes and diseases, have been implicated in pathogenesis of BA. However, a global view of TFs and how they link to clinical presentations remain explored. Here, we perform a joint transcriptional regulatory network and protein activity inference analysis in order to investigate transcription factor activity in BA. By integration of three independent human BA liver transcriptome datasets, we identify 22 common master regulators, with 14 activated- and 8 repressed TFs. Gene targets of activated TFs are enriched in biological processes of SMAD, NF-kappaB and TGF-beta, while those of repressed TFs are related to lipid metabolism. Mining the clinical association of TFs, we identify inflammation-, fibrosis- and survival associated TFs. In particular, ZNF14 is predictive of poor survival and advanced live fibrosis. Supporting this observation, ZNF14 is positively correlated with T helper cells, cholangiocytes and hepatic stellate cells. In sum, our analysis reveals key clinically associated master regulators for BA.

The role of p.Val444Ala variant in the ABCB11 gene and susceptibility to biliary atresia in Vietnamese patients

Biliary atresia (BA) is the most serious type of obstructive cholangiopathy that occurs in infants. BA can be the cause of death in children under 2 years if untreated early. However, the etiology of the disease is not known. BA is considered to be the result of the destruction of the bile duct system including the accumulation of bile acids. The bile salt export pump, a transporter protein encoded by the ABCB11 gene, plays the main role in the exportation and accumulation of bile acids. The p.Val444Ala variant in this gene is known to be associated with many cholestatic diseases. However, to date no study have been performed to evaluate the association of this variant with susceptibility to the risk of BA. In this study, we aimed to identify the frequency of p.Val444Ala variant and the risk of BA in Vietnamese patients.The polymerase chain reaction (PCR)- restriction fragment length polymorphism method was used to determine the frequency of alleles c.1331T>C (p.Val444Ala, rs2287622) in the ABCB11 gene in 266 Vietnamese patients with BA and 150 healthy people. The gene segment containing the variant was amplified by PCR with specific primers, after that the PCR products were cut by HaeIII restriction enzyme and analyzed on agarose gel to determine the genotypes. The frequency of alleles was assessed statistically to determine the association between these alleles and the risk of disease in patients.In our study, the frequency of alleles c.1331T>C (p.Val444Ala, rs2287622) in the ABCB11 gene was investigated the first time in the patients with BA. The results showed that CC and TC genotypes were significantly different between BA patients and healthy people (P < .01), and the C allele was associated with an increased risk of BA (odds ratio = 2.47; 95% confidence interval: 1.84-3.32; P < .01). The initial results of clinical, biochemical, and genetic analysis in our study suggested that the p.Val444Ala variant in the ABCB11 gene may be a susceptibility factor for the disease in Vietnamese patients with BA. These results provided new insights into the role of this ABCB11 variant in the pathogenesis of BA.

Novel findings from family-based exome sequencing for children with biliary atresia

Biliary atresia (BA) is a progressive inflammation and fibrosis of the biliary tree characterized by the obstruction of bile flow, which results in liver failure, scarring and cirrhosis. This study aimed to explore the elusive aetiology of BA by conducting whole exome sequencing for 41 children with BA and their parents (35 trios, including 1 family with 2 BA-diagnosed children and 5 child-mother cases). We exclusively identified and validated a total of 28 variants (17 X-linked, 6 de novo and 5 homozygous) in 25 candidate genes from our BA cohort. These variants were among the 10% most deleterious and had a low minor allele frequency against the employed databases: Kinh Vietnamese (KHV), GnomAD and 1000 Genome Project. Interestingly, AMER1, INVS and OCRL variants were found in unrelated probands and were first reported in a BA cohort. Liver specimens and blood samples showed identical variants, suggesting that somatic variants were unlikely to occur during morphogenesis. Consistent with earlier attempts, this study implicated genetic heterogeneity and non-Mendelian inheritance of BA.

RNA Splicing: A Versatile Regulatory Mechanism in Pediatric Liver Diseases

With the development of high-throughput sequencing technology, the posttranscriptional mechanism of alternative splicing is becoming better understood. From decades of studies, alternative splicing has been shown to occur in multiple tissues, including the brain, heart, testis, skeletal muscle, and liver. This regulatory mechanism plays an important role in physiological functions in most liver diseases. Currently, due to the absence of symptoms, chronic pediatric liver diseases have a significant impact on public health. Furthermore, the progression of the disease is accelerated in children, leading to severe damage to their liver tissue if no precautions are taken. To this end, this review article summarizes the current knowledge of alternative splicing in pediatric liver diseases, paying special attention to liver damage in the child stage. The discussion of the regulatory role of splicing in liver diseases and its potential as a new therapeutic target is also included.

NCK-associated protein 1 like (nckap1l) minor splice variant regulates intrahepatic biliary network morphogenesis

Impaired formation of the intrahepatic biliary network leads to cholestatic liver diseases, which are frequently associated with autoimmune disorders. Using a chemical mutagenesis strategy in zebrafish combined with computational network analysis, we screened for novel genes involved in intrahepatic biliary network formation. We positionally cloned a mutation in the nckap1l gene, which encodes a cytoplasmic adaptor protein for the WAVE regulatory complex. The mutation is located in the last exon after the stop codon of the primary splice isoform, only disrupting a previously unannotated minor splice isoform, which indicates that the minor splice isoform is responsible for the intrahepatic biliary network phenotype. CRISPR/Cas9-mediated nckap1l deletion, which disrupts both the primary and minor isoforms, showed the same defects. In the liver of nckap1l mutant larvae, WAVE regulatory complex component proteins are degraded specifically in biliary epithelial cells, which line the intrahepatic biliary network, thus disrupting the actin organization of these cells. We further show that nckap1l genetically interacts with the Cdk5 pathway in biliary epithelial cells. These data together indicate that although nckap1l was previously considered to be a hematopoietic cell lineage-specific protein, its minor splice isoform acts in biliary epithelial cells to regulate intrahepatic biliary network formation.

Genetic aspects of biliary atresia etiology

Biliary atresia (BA) is a cholestatic disorder of infancy that is fatal if untreated. Despite years of study the etiology of BA remains unknown. Three etiopathogenic mechanisms may be involved, such as immune dysregulation, environmental factors and genetic susceptibility. Genetic predisposition is being actively studied. Candidate genes associated with BA in certain populations, genes affecting the cholangiocyte cilia function, as well as genes involved in stress responses have been identified. However, the long-term follow-up of twins with BA suggests that genotype is not of paramount importance for the disease development. Both epigenetic patterns and postzygotic somatic mutations may contribute to etiology of the disease. Recently, some evidence is being accumulated on the possible genetic predisposition to certain outcome of Kasai portoenterostomy performed in patients with BA. However, the presence of a number of factors contributing to the development of the disease makes it difficult to identify the genetic markers.

Cholangitis in Patients With Biliary Atresia Receiving Hepatoportoenterostomy: A National Database Study

INTRODUCTION

Biliary atresia (BA) is a progressive form of liver disease in the neonatal period usually requiring hepatoportoenterostomy (HPE). Cholangitis is a common sequelae of HPE but data about which patients are at risk for this complication are limited.

OBJECTIVE

The objective of the study was to determine risk factors associated with cholangitis in a large retrospective cohort after HPE.

METHODS

The Pediatric Health Information System (PHIS) was queried for BA (ICD-9 975.61) and HPE (ICD-9-CM 51.37) admissions from 2004 to 2013. We performed univariate analysis and linear regression with dependent variables of ≥ 2 or ≥ 5 episodes of cholangitis, and independent variables of age at time of HPE, race, ethnicity, gender, insurance, ursodeoxycholic acid (UDCA) use, steroid use, presence of esophageal varices (EV), and portal hypertension (PH).

RESULTS

We identified 1112 subjects with a median age at HPE of 63 days and median number of cholangitis episodes of 2 within 2 years. On multiple regression analysis, black race (odds ratio (OR) 1.51, P = 0.044) and presence of PH (OR 2.24, P < 0.001) were associated with increased risk of ≥ 2 episodes of cholangitis, whereas HPE at >90 days was associated with less risk (OR 0.46, P = 0.001). Among those with ≥5 episodes, Asian race (OR 2.66, P = 0.038), public insurance (OR 1.72, P = 0.043), EV (OR 1.81, P = 0.017), and PH (OR 2.88, P < 0.001) were associated with higher risk.

CONCLUSIONS

Complications, such as cholangitis remain a common problem for patients, after HPE, with median of 2 episodes within 2 years. Higher rates of cholangitis are associated with portal hypertension whereas lower rate is associated with age at HPE of >90 days. Asians, patients with public insurance, and those with portal hypertension are more likely to have recurrent cholangitis.

Impaired Redox and Protein Homeostasis as Risk Factors and Therapeutic Targets in Toxin-Induced Biliary Atresia

BACKGROUND & AIMS

Extrahepatic biliary atresia (BA) is a pediatric liver disease with no approved medical therapy. Recent studies using human samples and experimental modeling suggest that glutathione redox metabolism and heterogeneity play a role in disease pathogenesis. We sought to dissect the mechanistic basis of liver redox variation and explore how other stress responses affect cholangiocyte injury in BA.

METHODS

We performed quantitative in situ hepatic glutathione redox mapping in zebrafish larvae carrying targeted mutations in glutathione metabolism genes and correlated these findings with sensitivity to the plant-derived BA-linked toxin biliatresone. We also determined whether genetic disruption of HSP90 protein quality control pathway genes implicated in human BA altered biliatresone toxicity in zebrafish and human cholangiocytes. An in vivo screening of a known drug library was performed to identify novel modifiers of cholangiocyte injury in the zebrafish experimental BA model, with subsequent validation.

RESULTS

Glutathione metabolism gene mutations caused regionally distinct changes in the redox potential of cholangiocytes that differentially sensitized them to biliatresone. Disruption of human BA-implicated HSP90 pathway genes sensitized zebrafish and human cholangiocytes to biliatresone-induced injury independent of glutathione. Phosphodiesterase-5 inhibitors and other cyclic guanosine monophosphate signaling activators worked synergistically with the glutathione precursor N-acetylcysteine in preventing biliatresone-induced injury in zebrafish and human cholangiocytes. Phosphodiesterase-5 inhibitors enhanced proteasomal degradation and required intact HSP90 chaperone.

CONCLUSION

Regional variation in glutathione metabolism underlies sensitivity to the biliary toxin biliatresone and may account for the reported association between BA transplant-free survival and glutathione metabolism gene expression. Human BA can be causatively linked to genetic modulation of protein quality control. Combined treatment with N-acetylcysteine and cyclic guanosine monophosphate signaling enhancers warrants further investigation as therapy for BA.

Association of common variation in ADD3 and GPC1 with biliary atresia susceptibility

Biliary atresia (BA) is an idiopathic neonatal cholestatic disease. Recent genome-wide association study (GWAS) revealed that common variation of ADD3, GPC1, ARF6, and EFEMP1 gene was associated with BA susceptibility. We aimed to evaluate the association of these genes with BA in Chinese population. Twenty single nucleotide polymorphisms (SNPs) in these four genes were genotyped in 340 BA patients and 1,665 controls. Three SNPs in ADD3 were significantly associated with BA, and rs17095355 was the top SNP (P_Allele = 3.23×10^-6). Meta-analysis of published data and current data indicated that rs17095355 was associated with BA susceptibility in Asians and Caucasians. Three associated SNPs were expression quantitative trait loci (eQTL) for ADD3. Two GPC1 SNPs in high linkage disequilibrium (LD) showed nominal association with BA susceptibility (P_Allele = 0.03 for rs6707262 and P_Allele = 0.04 for rs6750380), and were eQTL of GPC1. Haplotype harboring these two SNPs almost reached the study-wide significance (P = 0.0035). No association for ARF6 and EFEMP1 was found with BA risk in the current population. Our study validated associations of ADD3 and GPC1 SNPs with BA risk in Chinese population and provided evidence of epistatic contributions of genetic factors to BA susceptibility.

Exome Sequencing in Individuals with Isolated Biliary Atresia

Biliary atresia (BA) is a severe pediatric liver disease resulting in necroinflammatory obliteration of the extrahepatic biliary tree. BA presents within the first few months of life as either an isolated finding or with additional syndromic features. The etiology of isolated BA is unknown, with evidence for infectious, environmental, and genetic risk factors described. However, to date, there are no definitive causal genes identified for isolated BA in humans, and the question of whether single gene defects play a major role remains open. We performed exome-sequencing in 101 North American patients of European descent with isolated BA (including 30 parent-child trios) and considered several experimental designs to identify potentially deleterious protein-altering variants that may be involved in the disease. In a case-only analysis, we did not identify genes with variants shared among more than two probands, and burden tests of rare variants using a case-case control design did not yield significant results. In the trio analysis of 30 simplex families (patient and parent trios), we identified 66 de novo variants in 66 genes including potentially deleterious variants in STIP1 and REV1. STIP1 is a co-chaperone for the heat-shock protein, HSP90, and has been shown to have diverse functions in yeast, flies and mammals, including stress-responses. REV1 is known to be a key player in DNA repair pathway and to interact with HSP90. In conclusion, our results do not support the hypothesis that a simple genetic model is responsible for the majority of cases of isolated BA. Our finding of de novo variants in genes linked to evolutionarily conserved stress responses (STIP1 and REV1) suggests that exploration of how genetic susceptibility and environmental exposure may interact to cause BA is warranted.

TGF-β Signaling Plays a Pivotal Role During Developmental Biliary Atresia in Sea Lamprey (Petromyzon marinus)

Biliary atresia (BA) is a rare neonatal disease with unknown causes. Approximately 10% of BA cases develop in utero with other congenital defects that span a large spectrum of disease variations, including degeneration of the gall bladder and bile duct as well as malformation of the liver, intestines, and kidneys. Similar developmental alterations are manifested in a unique animal model, the sea lamprey (Petromyzon marinus), in which BA occurs naturally during metamorphosis. With the likelihood of conserved developmental mechanisms underlying organogenesis and degeneration, lamprey developmental BA may be a useful model to infer mechanisms underlying human embryonic BA. We reasoned that hepatobiliary transcriptomes regulate the transition between landmark stages of BA. Therefore, we examined sea lamprey hepatobiliary transcriptomes at four stages (M0, metamorphic stage 0 or larval stage, no BA; M2, metamorphic stage 2, onset of BA; M5, metamorphic stage 5, BA, and heightened hepatocyte proliferation and reorganization; and JV, juvenile, completion of BA) using messenger RNA sequencing and Kyoto Encyclopedia of Genes and Genomes pathway analyses. We found gene-expression patterns associated with the transition between these stages. In particular, transforming growth factor β (TGF-β), hedgehog, phosphatidylinositol-4,5-bisphosphate 3-kinase-Akt, Wnt, and mitogen-activated protein kinase pathways were involved during biliary degeneration. Furthermore, disrupting the TGF-β signaling pathway with antagonist or small interfering RNA treatments at the onset of BA delayed gall bladder and bile duct degeneration. Conclusion: Distinctive gene-expression patterns are associated with the degeneration of the biliary system during developmental BA. In addition, disrupting TGF-β signaling pathway at the onset of BA delayed biliary degeneration.

Genetics in biliary atresia

PURPOSE OF REVIEW

Biliary atresia is a poorly understood deadly disease. Genetic predisposition factors are suspected albeit not firmly established. This review summarizes recent evidence of genetic alterations in biliary atresia.

RECENT FINDINGS

Whole-genome association studies in biliary atresia patients identified four distinct predisposition loci with four different genes potentially involved in the disease occurrence. Variations in these genes were searched for, but none were found in patients with biliary atresia suggesting complex mechanisms.

SUMMARY

Despite decades since its description and decades of intensive researches, cause of biliary atresia disease remains enigmatic. The inheritance of biliary atresia is not Mendelian. Genetic predisposition factor is one of the explored fields to explain biliary atresia pathogenicity. Biliary atresia has been associated with several inborn syndromes, chromosome anomalies, and gene polymorphisms in specific populations. Four predisposition loci encompassing genes relevant to the disease have been identified, but no pathogenic variations were found in biliary atresia patients. Few reported cases of isolated biliary atresia manifestation in the context of known genetic diseases suggest coincidental findings. Alternatives to classic genetic alterations are proposed to explain genetic predisposition in biliary atresia including noncoding and epigenetic factors. Biliary atresia is most likely related to complex traits making its genetic exploration challenging.

Modeling Heterogeneity in the Genetic Architecture of Ethnically Diverse Groups Using Random Effect Interaction Models

In humans, most genome-wide association studies have been conducted using data from Caucasians and many of the reported findings have not replicated in other populations. This lack of replication may be due to statistical issues (small sample sizes or confounding) or perhaps more fundamentally to differences in the genetic architecture of traits between ethnically diverse subpopulations. What aspects of the genetic architecture of traits vary between subpopulations and how can this be quantified? We consider studying effect heterogeneity using Bayesian random effect interaction models. The proposed methodology can be applied using shrinkage and variable selection methods, and produces useful information about effect heterogeneity in the form of whole-genome summaries (e.g., the proportions of variance of a complex trait explained by a set of SNPs and the average correlation of effects) as well as SNP-specific attributes. Using simulations, we show that the proposed methodology yields (nearly) unbiased estimates when the sample size is not too small relative to the number of SNPs used. Subsequently, we used the methodology for the analyses of four complex human traits (standing height, high-density lipoprotein, low-density lipoprotein, and serum urate levels) in European-Americans (EAs) and African-Americans (AAs). The estimated correlations of effects between the two subpopulations were well below unity for all the traits, ranging from 0.73 to 0.50. The extent of effect heterogeneity varied between traits and SNP sets. Height showed less differences in SNP effects between AAs and EAs whereas HDL, a trait highly influenced by lifestyle, exhibited a greater extent of effect heterogeneity. For all the traits, we observed substantial variability in effect heterogeneity across SNPs, suggesting that effect heterogeneity varies between regions of the genome.

Biliary Atresia Relevant Human Induced Pluripotent Stem Cells Recapitulate Key Disease Features in a Dish

Biliary atresia (BA) is the most common cause of pediatric end-stage liver disease and the etiology is poorly understood. There is no effective therapy for BA partly due to lack of human BA models. Towards developing in vitro human models of BA, disease-specific induced pluripotent stem cells (iPSCs) from 6 BA patients were generated using non-integrating episomal plasmids. In addition, to determine the functional significance of BA-susceptibility genes identified by genome-wide association studies (GWAS) in biliary development, a genome-editing approach was used to create iPSCs with defined mutations in these GWAS BA loci. Using the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 system, isogenic iPSCs deficient in BA-associated genes (GPC1 and ADD3) were created from healthy iPSCs. Both the BA patient-iPSCs and the knock out (KO) iPSCs were studied for their in vitro biliary differentiation potential. These BA-specific iPSCs demonstrated significantly decreased formation of ductal structures, decreased expression of biliary markers including CK7, EpCAM, SOX9, CK19, AE2, and CFTR and increased fibrosis markers such as alpha smooth muscle actin, Loxl2, and Collagen1 compared to controls. Both the patient- and the KO-iPSCs also showed increased yes-associated protein (YAP, a marker of bile duct proliferation/fibrosis). Collagen and YAP were reduced by treatment with the anti-fibrogenic drug pentoxifylline. In summary, these BA-specific human iPSCs showed deficiency in biliary differentiation along with increased fibrosis, the 2 key disease features of BA. These iPSCs can provide new human BA models for understanding the molecular basis of abnormal biliary development and opportunities to identify drugs that have therapeutic effects on BA.

Genetic investigation into an increased susceptibility to biliary atresia in an extended New Zealand Māori family

Background

Biliary atresia (BA), a fibrosing disorder of the developing biliary tract leading to liver failure in infancy, has an elevated incidence in indigenous New Zealand (NZ) Māori. We investigated a high rate of BA in a group of children (n = 12) belonging to a single Māori iwi (or ‘tribe’, related through a remote ancestor).

Methods

Population and geographical data was used to estimate the rate of BA in Māori sub-groups, and a pedigree linking most of the affected children was constructed from oral and documented history. Array genotyping was used to examine hypotheses about the inheritance of a possible genetic risk factor, and the history of the affected population, and Exome Sequencing to search for candidate genes.

Results

Most of these affected children (n = 7) link to a self-reported pedigree and carry a 50-fold increase in BA risk over unrelated Māori (χ² = 296P < 0.001, 95% CI 23–111). Genetic analysis using FEstim and SNP array genotypes revealed no evidence for elevated consanguinity between parents of affected children (FEstim: F (2,21) = 0.469, P > 0.63). Genome-wide quantitation of intervals of contiguous, homozygous-by-state markers reached a similar conclusion (F (2,399) = 1.99, P = 0.138). Principal component analysis and investigation with STRUCTURE found no evidence of increased allele frequency of either a recessive variant, or additive, low-risk variants due to reproductive isolation. To identify candidate causal factors, Exome Sequencing datasets were scrutinised for shared rare coding variants across 8 affected individuals. No rare, non-synonymous, phylogenetically conserved variants were common to 6 or more affected children.

Conclusion

The substantially elevated risk for development of BA in this subgroup could be mediated by genetic factors, but the iwi exhibits no properties indicative of recent or remote reproductive isolation. Resolution of any risk loci may rely on extensive genomic sequencing studies in this iwi or investigation of other mechnaisms such as copy number variation.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0440-0) contains supplementary material, which is available to authorized users.

Biliary Atresia: Clinical and Research Challenges for the Twenty-First Century

Biliary atresia (BA) is a fibroinflammatory disease of the intrahepatic and extrahepatic biliary tree. Surgical hepatic portoenterostomy (HPE) may restore bile drainage, but progression of the intrahepatic disease results in complications of portal hypertension and advanced cirrhosis in most children. Recognizing that further progress in the field is unlikely without a better understanding of the underlying cause(s) and pathogenesis of the disease, the National Institutes of Diabetes and Digestive and Kidney Diseases (NIDDK) sponsored a research workshop focused on innovative and promising approaches and on identifying future areas of research. Investigators discussed recent advances using gestational ultrasound and results of newborn BA screening with serum direct (conjugated) bilirubin that support a prenatal onset of biliary injury. Experimental and human studies implicate the toxic properties of environmental toxins (e.g., biliatresone) and of viruses (e.g., cytomegalovirus) to the biliary system. Among host factors, sequence variants in genes related to biliary development and ciliopathies, a notable lack of a cholangiocyte glycocalyx and of submucosal collagen bundles in the neonatal extrahepatic bile ducts, and an innate proinflammatory bias of the neonatal immune system contribute to an increased susceptibility to damage and obstruction following epithelial injury. These advances form the foundation for a future research agenda focused on identifying the environmental and host factor(s) that cause BA, the potential use of population screening, studies of the mechanisms of prominent fibrosis in young infants, determinations of clinical surrogates of disease progression, and the design of clinical trials that target subgroups of patients with initial drainage following HPE. (Hepatology 2018; 00:000-000).

A genome-wide association study identifies a susceptibility locus for biliary atresia on 2p16.1 within the gene EFEMP1

Biliary atresia (BA) is a rare pediatric cholangiopathy characterized by fibrosclerosing obliteration of the extrahepatic bile ducts, leading to cholestasis, fibrosis, cirrhosis, and eventual liver failure. The etiology of BA remains unknown, although environmental, inflammatory, infectious, and genetic risk factors have been proposed. We performed a genome-wide association study (GWAS) in a European-American cohort of 343 isolated BA patients and 1716 controls to identify genetic loci associated with BA. A second GWAS was performed in an independent European-American cohort of 156 patients with BA and other extrahepatic anomalies and 212 controls to confirm the identified candidate BA-associated SNPs. Meta-analysis revealed three genome-wide significant BA-associated SNPs on 2p16.1 (rs10865291, rs6761893, and rs727878; P < 5 ×10-8), located within the fifth intron of the EFEMP1 gene, which encodes a secreted extracellular protein implicated in extracellular matrix remodeling, cell proliferation, and organogenesis. RNA expression analysis showed an increase in EFEMP1 transcripts from human liver specimens isolated from patients with either BA or other cholestatic diseases when compared to normal control liver samples. Immunohistochemistry demonstrated that EFEMP1 is expressed in cholangiocytes and vascular smooth muscle cells in liver specimens from patients with BA and other cholestatic diseases, but it is absent from cholangiocytes in normal control liver samples. Efemp1 transcripts had higher expression in cholangiocytes and portal fibroblasts as compared with other cell types in normal rat liver. The identification of a novel BA-associated locus, and implication of EFEMP1 as a new BA candidate susceptibility gene, could provide new insights to understanding the mechanisms underlying this severe pediatric disorder.

A Review on Adducin from Functional to Pathological Mechanisms: Future Direction in Cancer

Adducin (ADD) is a family of membrane skeleton proteins including ADD1, ADD2, and ADD3 that are encoded by distinct genes on different chromosomes. Adducin is primarily responsible for the assembly of spectrin-actin network that provides physical support to the plasma membrane and mediates signal transduction in various cellular physiological processes upon regulation by protein kinase C-dependent and calcium/calmodulin-dependent pathways. Abnormal phosphorylation, genetic variations, and alternative splicing of adducin may contribute to alterations in cellular functions involved in pathogenic processes. These alterations are associated with a wide range of diseases including cancer. This paper begins with a discussion on how adducin partakes in the structural formation of membrane skeleton, its regulation, and related functional characteristics, followed by a review on the pathogenesis of hypertension, biliary atresia, and cancer with respect to increased disease susceptibility mediated by adducin polymorphism and/or dysregulation. Given the functional diversity of adducin in different cellular compartments, we aim to provide a knowledge base whereby its pathophysiological roles can be better understood. More importantly, we aim to provide novel insights that may be of significance in turning the adducin model to clinical application.

Single nucleotide polymorphisms within Adducin 3 and Adducin 3 antisense RNA1 genes are associated with biliary atresia in Thai infants

BACKGROUND

A genome-wide association study in East Asians suggested a genetic association between biliary atresia (BA) and a cluster of variants within the Adducin 3 (ADD3) and ADD3 antisense RNA1 (ADD3-AS1) genes. Another study in Thai neonates reported an association between BA and rs17095355. To validate those findings, this study aimed to analyze the BA association with single nucleotide polymorphisms (SNPs) and the additive influence of ADD3 and ADD3-AS1 in Thai neonates.

METHODS

DNAs from 56 BA cases and 166 controls were genotyped for rs2501577, rs11194981, rs12268910 (ADD3) and rs17095355 (ADD3-AS1), using TaqMan PCR. Genotype distributions were compared between the groups, and SNP-SNP interactions were analyzed by combination of allelotypes.

RESULTS

The risk allele frequencies of rs2501577, rs11194981, and rs17095355 in the BA group were significantly higher than in the controls. Univariate analysis showed that recessive variants in the three SNPs were associated with BA risk at ORs of 1.81 (95% CI 1.32-2.50), 1.58 (95% CI 1.14-2.20) and 1.92 (95% CI 1.39-2.66), respectively. SNP-SNP interaction analysis showed that the SNP combination of the two genes rs17095355 and rs2501577 provided an additive increase in BA risk.

CONCLUSION

ADD3 and ADD3-AS1 variants increased susceptibility to BA, suggesting that these genes may play an additive role in the pathogenesis of the disease. In addition, these interactions may give a clue to the overexpression of the ADD3 protein in the liver of BA patients.

Gene-disease associations identify a connectome with shared molecular pathways in human cholangiopathies

Cholangiopathies are a diverse group of progressive diseases whose primary cell targets are cholangiocytes. To identify shared pathogenesis and molecular connectivity among the three main human cholangiopathies (biliary atresia [BA], primary biliary cholangitis [PBC], and primary sclerosing cholangitis [PSC]), we built a comprehensive platform of published data on gene variants, gene expression, and functional studies and applied network-based analytics in the search for shared molecular circuits. Mining the data platform with largest connected component and interactome analyses, we validated previously reported associations and identified essential and hub genes. In addition to disease-specific modules, we found a substantial overlap of disease neighborhoods and uncovered a group of 34 core genes that are enriched for immune processes and abnormal intestine/hepatobiliary mouse phenotypes. Within this core, we identified a gene subcore containing signal transduction and activator of transcription 3, interleukin-6, tumor necrosis factor, and forkhead box P3 prominently placed in a regulatory connectome of genes related to cellular immunity and fibrosis. We also found substantial gene enrichment in the advanced glycation endproduct/receptor for advanced glycation endproducts (RAGE) pathway and showed that RAGE activation induced cholangiocyte proliferation. Conclusion: Human cholangiopathies share pathways enriched by immunity genes and a molecular connectome that links different pathogenic features of BA, PBC, and PSC. (Hepatology 2018;67:676-689).

Update on investigations pertaining to the pathogenesis of biliary atresia

Biliary atresia is a devastating biliary disease of neonates that results in liver transplantation for the vast majority. The etiology of biliary atresia is unknown and is likely multifactorial, with components of genetic predisposition, environmental trigger and autoimmunity contributing to disease pathogenesis. This review highlights recent work related to investigations of disease pathogenesis in biliary atresia.

Three-dimensional structural analysis reveals a Cdk5-mediated kinase cascade regulating hepatic biliary network branching in zebrafish

The intrahepatic biliary network is a highly branched three-dimensional network lined by biliary epithelial cells, but how its branching patterns are precisely established is not clear. We designed a new computer-based algorithm that quantitatively computes the structural differences of the three-dimensional networks. Utilizing the algorithm, we showed that inhibition of Cyclin-dependent kinase 5 (Cdk5) led to reduced branching in the intrahepatic biliary network in zebrafish. Further, we identified a previously unappreciated downstream kinase cascade regulated by Cdk5. Pharmacological manipulations of this downstream kinase cascade produced a crowded branching defect in the intrahepatic biliary network and influenced actin dynamics in biliary epithelial cells. We generated larvae carrying a mutation in cdk5 regulatory subunit 1a (cdk5r1a), an essential activator of Cdk5. cdk5r1a mutant larvae show similar branching defects as those observed in Cdk5 inhibitor-treated larvae. A small-molecule compound that interferes with the downstream kinase cascade rescued the mutant phenotype. These results provide new insights into branching morphogenesis of the intrahepatic biliary network.

Patient complexity and genotype-phenotype correlations in biliary atresia: a cross-sectional analysis

Background

Biliary Atresia (BA) is rare and genetically complex, and the pathogenesis is elusive. The disease course is variable and can represent heterogeneity, which hinders effective disease management. Deciphering the BA phenotypic variance is a priority in clinics and can be achieved by the integrative analysis of genotype and phenotype. We aim to explore the BA phenotypic features and to delineate the source of its variance.

Methods

The study is a cross-sectional observational study collating with case/control association analysis. One-hundred-and-eighty-one type III non-syndromic BA patients and 431 controls were included for case–control association tests, including 89 patients (47.19% males, born June 15th, 1981 to September 17th, 2007) have detailed clinical records with follow-up of the disease course (median ~17.2 years). BA-association genes from the genome-wide gene-based association test on common genetic variants (CV) and rare copy-number-variants (CNVs) from the genome-wide survey, the later comprise only CNVs > 100 kb and found in the BA patients but not in the local population (N = 1,381) or the database (N = 11,943). Hereby comorbidity is defined as a chronic disease that affects the BA patients but has no known relationship with BA or with the BA treatment. We examined genotype-phenotype correlations of CNVs, connectivity of these novel variants with BA-associated CVs, and their role in the BA candidate gene network.

Results

Of the 89 patients, 41.57% have comorbidities, including autoimmune-allergic disorders (22.47%). They carried 29 BA-private CNVs, including 3 CNVs underpinning the carriers’ immunity comorbidity and one JAG1 micro-deletion. The BA-CNV-intersected genes (N = 102) and the CV-tagged genes (N = 103) were both enriched with immune-inflammatory pathway genes (FDR q < 0.20), and the two gene sets were interconnected (permutation p = 0.039). The molecular network representing CVs and rare-CNV association genes fit into a core/periphery structure, the immune genes and their related modules are found at the coherence core of all connections, suggesting its dominant role in the BA pathogenesis pathway.

Conclusions

The study highlights a patient-complexity phenomenon as a novel BA phenotypic feature, which is underpinned by rare-CNVs that biologically converge with CVs into the immune-inflammatory pathway and drives the BA occurrence and the likely BA association with immune diseases in clinics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-017-0259-0) contains supplementary material, which is available to authorized users.

Downregulation of microRNA-145 may contribute to liver fibrosis in biliary atresia by targeting ADD3

BACKGROUND AND OBJECTIVES

Biliary atresia (BA) is a pediatric liver disease characterized by fibro-obliteration and obstruction of the extrahepatic biliary system, that invariably leads to cirrhosis and even death, if left untreated for extended time. However, its pathology and etiology still remained unknown. In this study, we tested the expression of adducin 3 (ADD3), the gene identified as a susceptibility gene in BA by GWAS, and uncovered its upstream regulatory microRNA in the pathogenesis of BA.

METHODS

In this study, 14 infants with BA and 14 infants with choledochal cyst (CC) were enrolled as experimental group and control group, respectively. ADD3 and microRNA-145 (miR-145) expression profiles in liver tissues of BA and CC were determined using qPCR. Luciferase reporter assay was performed to verify the direct interaction between miR-145-5p and ADD3 3' Untranslated Regions (3'UTR). The Lentiviral vectors containing miR-145, miR-145-3p inhibitor, miR-145-5p inhibitor, empty vector were transfected into human hepatic stellate cell line (LX-2) to determine the functional effect of miR-145 on ADD3 expression at both mRNA and protein level.

RESULTS

MiR-145 was shown to be down-regulated in liver tissues of infants with BA compared to CC (p = 0.0267). ADD3, verified as a target of miR-145-5p, was shown to be overexpressed in infants with BA at the mRNA level (p = 0.0118). Transfection of lentiviruses containing miR-145 into LX-2 cells decreased the expression of ADD3 at both mRNA and protein level compared to negative control group, and suppressed the expression of p-Akt at protein level.

CONCLUSIONS

Our study has shown that overexpressed ADD3 and downregulated miR-145 were detected in BA liver tissues. MiR-145-5p was confirmed to target ADD3 by luciferase reporter assay. The downregulation of miR-145 may contribute to liver fibrosis in BA by upregulating the expression of ADD3.

Loss of a Candidate Biliary Atresia Susceptibility Gene, add3a, Causes Biliary Developmental Defects in Zebrafish

OBJECTIVES

Biliary atresia (BA) is a progressive fibroinflammatory cholangiopathy affecting the bile ducts of neonates. Although BA is the leading indication for pediatric liver transplantation, the etiology remains elusive. Adducin 3 (ADD3) and X-prolyl aminopeptidase 1 (XPNPEP1) are 2 genes previously identified in genome-wide association studies as potential BA susceptibility genes. Using zebrafish, we investigated the importance of ADD3 and XPNPEP1 in functional studies.

METHODS

To determine whether loss of either gene leads to biliary defects, we performed morpholino antisense oligonucleotide (MO) knockdown studies targeting add3a and xpnpep1 in zebrafish. Individuals were assessed for decreases in biliary function and the presence of biliary defects. Quantitative polymerase chain reaction was performed on pooled 5 days postfertilization larvae to assess variations in transcriptional expression of genes of interest.

RESULTS

Although both xpnpep1 and add3a are expressed in the developing zebrafish liver, only knockdown of add3a produced intrahepatic defects and decreased biliary function. Similar results were observed in homozygous add3a mutants. MO-mediated knockdown of add3a also showed higher mRNA expression of hedgehog (Hh) targets. Inhibition of Hh signaling rescued biliary defects caused by add3a knockdown. Combined knockdown of add3a and glypican-1 (gpc1), another mediator of Hh activity that is also a BA susceptibility gene, resulted in more severe biliary defects than knockdown of either alone.

CONCLUSIONS

Our results support previous studies identifying ADD3 as a putative genetic risk factor for BA susceptibility. Our results also provide evidence that add3a may be affecting the Hh pathway, an important factor in BA pathogenesis.

Association between rs17095355 polymorphism on 10q24 and susceptibility to biliary atresia: a meta-analysis

OBJECTIVE

Recent studies have identified 10q24-rs17095355 as a susceptibility locus for biliary atresia (BA). To more precisely estimate the association between the rs17095355 polymorphism and BA risk, a meta-analysis was performed.

METHODS

A comprehensive search was conducted to examine all the eligible studies by electronic databases including Elsevier Science Direct, Pubmed, Google Scholar, China National Knowledge Infrastructure (CNKI) and Chinese Biomedical Literature (CBM) up to December 2015. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to assess the strength of the association.

RESULTS

A total of 6 comparisons from 5 relevant studies involving 1000 patients and 3257 controls were included to analyze the association between rs17095355 and BA risk. The pooled OR for T allele of rs17095355 was 1.72 (95%CI 1.53-1.92, p < 0.01) in BA. Stratification by ethnicity indicated the degree of risk of rs17095355 with BA susceptibility was similar in populations of Asian origin. The pooled OR was 1.81 (95%CI 1.60-2.06, p < 0.01).

CONCLUSIONS

This meta-analysis confirms the association of rs17095355 polymorphism and BA development, especially in Asians. More original studies with large sample are needed to replicate this genetic association in different ethnic groups.

The toxin biliatresone causes mouse extrahepatic cholangiocyte damage and fibrosis through decreased glutathione and SOX17

Biliary atresia, the most common indication for pediatric liver transplantation, is a fibrotic disease of unknown etiology affecting the extrahepatic bile ducts of newborns. The recently described toxin biliatresone causes lumen obstruction in mouse cholangiocyte spheroids and represents a new model of biliary atresia. The goal of this study was to determine the cellular changes caused by biliatresone in mammalian cells that ultimately lead to biliary atresia and extrahepatic fibrosis. We treated mouse cholangiocytes in three-dimensional (3D) spheroid culture and neonatal extrahepatic duct explants with biliatresone and compounds that regulate glutathione (GSH). We examined the effects of biliatresone on SOX17 levels and determined the effects of Sox17 knockdown on cholangiocytes in 3D culture. We found that biliatresone caused disruption of cholangiocyte apical polarity and loss of monolayer integrity. Spheroids treated with biliatresone had increased permeability as shown by rhodamine efflux within 5 hours compared with untreated spheroids, which retained rhodamine for longer than 12 hours. Neonatal bile duct explants treated with the toxin showed lumen obstruction with increased subepithelial staining for α-smooth muscle actin and collagen, consistent with fibrosis. Biliatresone caused a rapid and transient decrease in GSH, which was both necessary and sufficient to mediate its effects in cholangiocyte spheroid and bile duct explant systems. It also caused a significant decrease in cholangiocyte levels of SOX17, and Sox17 knockdown in cholangiocyte spheroids mimicked the effects of biliatresone.

CONCLUSION

Biliatresone decreases GSH and SOX17 in mouse cholangiocytes. In 3D cell systems, this leads to cholangiocyte monolayer damage and increased permeability; in extrahepatic bile duct explants, it leads to disruption of the extrahepatic biliary tree and subepithelial fibrosis. This mechanism may be important in understanding human biliary atresia. (Hepatology 2016;64:880-893).

Glutathione antioxidant pathway activity and reserve determine toxicity and specificity of the biliary toxin biliatresone in zebrafish

Biliatresone is an electrophilic isoflavone isolated from Dysphania species plants that has been causatively linked to naturally occurring outbreaks of a biliary atresia (BA)-like disease in livestock. Biliatresone has selective toxicity for extrahepatic cholangiocytes (EHCs) in zebrafish larvae. To better understand its mechanism of toxicity, we performed transcriptional profiling of liver cells isolated from zebrafish larvae at the earliest stage of biliatresone-mediated biliary injury, with subsequent comparison of biliary and hepatocyte gene expression profiles. Transcripts encoded by genes involved in redox stress response, particularly those involved in glutathione (GSH) metabolism, were among the most prominently up-regulated in both cholangiocytes and hepatocytes of biliatresone-treated larvae. Consistent with these findings, hepatic GSH was depleted at the onset of biliary injury, and in situ mapping of the hepatic GSH redox potential using a redox-sensitive green fluorescent protein biosensor showed that it was significantly more oxidized in EHCs both before and after treatment with biliatresone. Pharmacological and genetic manipulation of GSH redox homeostasis confirmed the importance of GSH in modulating biliatresone-induced injury given that GSH depletion sensitized both EHCs and the otherwise resistant intrahepatic cholangiocytes to the toxin, whereas replenishing GSH level by N-acetylcysteine administration or activation of nuclear factor erythroid 2-like 2 (Nrf2), a transcriptional regulator of GSH synthesis, inhibited EHC injury.

CONCLUSION

These findings strongly support redox stress as a critical contributing factor in biliatresone-induced cholangiocyte injury, and suggest that variations in intrinsic stress responses underlie the susceptibility profile. Insufficient antioxidant capacity of EHCs may be critical to early pathogenesis of human BA. (Hepatology 2016;64:894-907).

Identification of a plant isoflavonoid that causes biliary atresia

Biliary atresia (BA) is a rapidly progressive and destructive fibrotic disorder of unknown etiology affecting the extrahepatic biliary tree of neonates. Epidemiological studies suggest that an environmental factor, such as a virus or toxin, is the cause of the disease, although none have been definitively established. Several naturally occurring outbreaks of BA in Australian livestock have been associated with the ingestion of unusual plants by pregnant animals during drought conditions. We used a biliary secretion assay in zebrafish to isolate a previously undescribed isoflavonoid, biliatresone, from Dysphania species implicated in a recent BA outbreak. This compound caused selective destruction of the extrahepatic, but not intrahepatic, biliary system of larval zebrafish. A mutation that enhanced biliatresone toxicity mapped to a region of the zebrafish genome that has conserved synteny with an established human BA susceptibility locus. The toxin also caused loss of cilia in neonatal mouse extrahepatic cholangiocytes in culture and disrupted cell polarity and monolayer integrity in cholangiocyte spheroids. Together, these findings provide direct evidence that BA could be initiated by perinatal exposure to an environmental toxin.

The Role of ARF6 in Biliary Atresia

Background & Aims

Altered extrahepatic bile ducts, gut, and cardiovascular anomalies constitute the variable phenotype of biliary atresia (BA).

Methods

To identify potential susceptibility loci, Caucasian children, normal (controls) and with BA (cases) at two US centers were compared at >550000 SNP loci. Systems biology analysis was carried out on the data. In order to validate a key gene identified in the analysis, biliary morphogenesis was evaluated in 2-5-day post-fertilization zebrafish embryos after morpholino-antisense oligonucleotide knockdown of the candidate gene ADP ribosylation factor-6 (ARF6, Mo-arf6).

Results

Among 39 and 24 cases at centers 1 and 2, respectively, and 1907 controls, which clustered together on principal component analysis, the SNPs rs3126184 and rs10140366 in a 3’ flanking enhancer region for ARF6 demonstrated higher minor allele frequencies (MAF) in each cohort, and 63 combined cases, compared with controls (0.286 vs. 0.131, P = 5.94x10^-7, OR 2.66; 0.286 vs. 0.13, P = 5.57x10^-7, OR 2.66). Significance was enhanced in 77 total cases, which included 14 additional BA genotyped at rs3126184 only (p = 1.58x10^-2, OR = 2.66). Pathway analysis of the 1000 top-ranked SNPs in CHP cases revealed enrichment of genes for EGF regulators (p<1 x10^-7), ERK/MAPK and CREB canonical pathways (p<1 x10^-34), and functional networks for cellular development and proliferation (p<1 x10^-45), further supporting the role of EGFR-ARF6 signaling in BA. In zebrafish embryos, Mo-arf6 injection resulted in a sparse intrahepatic biliary network, several biliary epithelial cell defects, and poor bile excretion to the gall bladder compared with uninjected embryos. Biliary defects were reproduced with the EGFR-blocker AG1478 alone or with Mo-arf6 at lower doses of each agent and rescued with arf6 mRNA.

Conclusions

The BA-associated SNPs identify a chromosome 14q21.3 susceptibility locus encompassing the ARF6 gene. arf6 knockdown in zebrafish implicates early biliary dysgenesis as a basis for BA, and also suggests a role for EGFR signaling in BA pathogenesis.

Biliary Atresia: Clinical Lessons Learned

Biliary atresia is a rare disease of unclear etiology, in which obstruction of the biliary tree causes severe cholestasis leading to cirrhosis and ultimately death if left untreated. Biliary atresia is the leading cause of neonatal cholestasis and the most frequent indication for pediatric liver transplantation. Any infant with persistent jaundice beyond 2 weeks of life needs to be evaluated for biliary atresia with fractionation of the bilirubin into conjugated and unconjugated portions. Early performance of a hepatoportoenterostomy in the first 45 days of life to restore bile flow and lessen further damage to the liver is thought to optimize outcome. Despite surgery, progressive liver scarring occurs, and 80% of patients with biliary atresia will require liver transplantation during childhood.

The Cholangiopathies

Cholangiocytes (ie, the epithelial cells that line the bile ducts) are an important subset of liver cells. They are actively involved in the modification of bile volume and composition, are activated by interactions with endogenous and exogenous stimuli (eg, microorganisms, drugs), and participate in liver injury and repair. The term cholangiopathies refers to a category of chronic liver diseases that share a central target: the cholangiocyte. The cholangiopathies account for substantial morbidity and mortality given their progressive nature, the challenges associated with clinical management, and the lack of effective medical therapies. Thus, cholangiopathies usually result in end-stage liver disease requiring liver transplant to extend survival. Approximately 16% of all liver transplants performed in the United States between 1988 and 2014 were for cholangiopathies. For all these reasons, cholangiopathies are an economic burden on patients, their families, and society. This review offers a concise summary of the biology of cholangiocytes and describes a conceptual framework for development of the cholangiopathies. We also present the recent progress made in understanding the pathogenesis of and how this knowledge has influenced therapies for the 6 common cholangiopathies-primary biliary cirrhosis, primary sclerosing cholangitis, cystic fibrosis involving the liver, biliary atresia, polycystic liver disease, and cholangiocarcinoma-because the latest scientific progress in the field concerns these conditions. We performed a search of the literature in PubMed for published papers using the following terms: cholangiocytes, biliary epithelia, cholestasis, cholangiopathy, and biliary disease. Studies had to be published in the past 5 years (from June 1, 2009, through May 31, 2014), and non-English studies were excluded.