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.

Mechanisms of Impaired Lung Development and Ciliation in Mannosidase-1-Alpha-2 (Man1a2) Mutants.

BACKGROUND: Ciliary defects cause heterogenous phenotypes related to mutation burden which lead to impaired development. A previously reported homozygous deletion in the Man1a2 gene causes lethal respiratory failure in newborn pups and decreased lung ciliation compared with wild type (WT) pups. The effects of heterozygous mutation, and the potential for rescue are not known. PURPOSE: We hypothesized that survival and lung ciliation, (a) would decrease progressively in Man1a2 +/- heterozygous and Man1a2 -/- null newborn pups compared with WT, and (b) could be enhanced by gestational treatment with N-Acetyl-cysteine (NAC), an antioxidant. METHODS: Man1a2+/- adult mice were fed NAC or placebo from a week before breeding through gestation. Survival of newborn pups was monitored for 24 h. Lungs, liver and tails were harvested for morphology, genotyping, and transcriptional profiling. RESULTS: Survival (p = 0.0001, Kaplan-Meier) and percent lung ciliation (p = 0.0001, ANOVA) measured by frequency of Arl13b+ respiratory epithelial cells decreased progressively, as hypothesized. Compared with placebo, gestational NAC treatment enhanced (a) lung ciliation in pups with each genotype, (b) survival in heterozygous pups (p = 0.017) but not in WT or null pups. Whole transcriptome of lung but not liver demonstrated patterns of up- and down-regulated genes that were identical in living heterozygous and WT pups, and completely opposite to those in dead heterozygous and null pups. Systems biology analysis enabled reconstruction of protein interaction networks that yielded functionally relevant modules and their interactions. In these networks, the mutant Man1a2 enzyme contributes to abnormal synthesis of proteins essential for lung development. The associated unfolded protein, hypoxic and oxidative stress responses can be mitigated with NAC. Comparisons with the developing human fetal lung transcriptome show that NAC likely restores normal vascular and epithelial tube morphogenesis in Man1a2 mutant mice. CONCLUSION: Survival and lung ciliation in the Man1a2 mutant mouse, and its improvement with N-Acetyl cysteine is genotype-dependent. NAC-mediated rescue depends on the central role for oxidative and hypoxic stress in regulating ciliary function and organogenesis during development.

A Newtonian framework for community detection in undirected biological networks.

Community detection is a key problem of interest in network analysis, with applications in a variety of domains such as biological networks, social network modeling, and communication pattern analysis. In this paper, we present a novel framework for community detection that is motivated by a physical system analogy. We model a network as a system of point masses, and drive the process of community detection, by leveraging the Newtonian interactions between the point masses. Our framework is designed to be generic and extensible relative to the model parameters that are most suited for the problem domain. We illustrate the applicability of our approach by applying the Newtonian Community Detection algorithm on protein-protein interaction networks of E. coli , C. elegans, and S. cerevisiae. We obtain results that are comparable in quality to those obtained from the Newman-Girvan algorithm, a widely employed divisive algorithm for community detection. We also present a detailed analysis of the structural properties of the communities produced by our proposed algorithm, together with a biological interpretation using E. coli protein network as a case study. A functional enrichment heat map is constructed with the Gene Ontology functional mapping, in addition to a pathway analysis for each community. The analysis illustrates that the proposed algorithm elicits communities that are not only meaningful from a topological standpoint, but also possess biological relevance. We believe that our algorithm has the potential to serve as a key computational tool for driving therapeutic applications involving targeted drug development for personalized care delivery.

Community Structure Analysis of Gene Interaction Networks in Duchenne Muscular Dystrophy.

Duchenne Muscular Dystrophy (DMD) is an important pathology associated with the human skeletal muscle and has been studied extensively. Gene expression measurements on skeletal muscle of patients afflicted with DMD provides the opportunity to understand the underlying mechanisms that lead to the pathology. Community structure analysis is a useful computational technique for understanding and modeling genetic interaction networks. In this paper, we leverage this technique in combination with gene expression measurements from normal and DMD patient skeletal muscle tissue to study the structure of genetic interactions in the context of DMD. We define a novel framework for transforming a raw dataset of gene expression measurements into an interaction network, and subsequently apply algorithms for community structure analysis for the extraction of topological communities. The emergent communities are analyzed from a biological standpoint in terms of their constituent biological pathways, and an interpretation that draws correlations between functional and structural organization of the genetic interactions is presented. We also compare these communities and associated functions in pathology against those in normal human skeletal muscle. In particular, differential enhancements are observed in the following pathways between pathological and normal cases: Metabolic, Focal adhesion, Regulation of actin cytoskeleton and Cell adhesion, and implication of these mechanisms are supported by prior work. Furthermore, our study also includes a gene-level analysis to identify genes that are involved in the coupling between the pathways of interest. We believe that our results serve to highlight important distinguishing features in the structural/functional organization of constituent biological pathways, as it relates to normal and DMD cases, and provide the mechanistic basis for further biological investigations into specific pathways differently regulated between normal and DMD patients. These findings have the potential to serve as fertile ground for therapeutic applications involving targeted drug development for DMD.

Modularity detection in protein-protein interaction networks.

BACKGROUND: Many recent studies have investigated modularity in biological networks, and its role in functional and structural characterization of constituent biomolecules. A technique that has shown considerable promise in the domain of modularity detection is the Newman and Girvan (NG) algorithm, which relies on the number of shortest-paths across pairs of vertices in the network traversing a given edge, referred to as the betweenness of that edge. The edge with the highest betweenness is iteratively eliminated from the network, with the betweenness of the remaining edges recalculated in every iteration. This generates a complete dendrogram, from which modules are extracted by applying a quality metric called modularity denoted by Q. This exhaustive computation can be prohibitively expensive for large networks such as Protein-Protein Interaction Networks. In this paper, we present a novel optimization to the modularity detection algorithm, in terms of an efficient termination criterion based on a target edge betweenness value, using which the process of iterative edge removal may be terminated. RESULTS: We validate the robustness of our approach by applying our algorithm on real-world protein-protein interaction networks of Yeast, C.Elegans and Drosophila, and demonstrate that our algorithm consistently has significant computational gains in terms of reduced runtime, when compared to the NG algorithm. Furthermore, our algorithm produces modules comparable to those from the NG algorithm, qualitatively and quantitatively. We illustrate this using comparison metrics such as module distribution, module membership cardinality, modularity Q, and Jaccard Similarity Coefficient. CONCLUSIONS: We have presented an optimized approach for efficient modularity detection in networks. The intuition driving our approach is the extraction of holistic measures of centrality from graphs, which are representative of inherent modular structure of the underlying network, and the application of those measures to efficiently guide the modularity detection process. We have empirically evaluated our approach in the specific context of real-world large scale biological networks, and have demonstrated significant savings in computational time while maintaining comparable quality of detected modules.