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.

Mechanism of KMT5B haploinsufficiency in neurodevelopment in humans and mice.

Pathogenic variants in KMT5B, a lysine methyltransferase, are associated with global developmental delay, macrocephaly, autism, and congenital anomalies (OMIM# 617788). Given the relatively recent discovery of this disorder, it has not been fully characterized. Deep phenotyping of the largest (n = 43) patient cohort to date identified that hypotonia and congenital heart defects are prominent features that were previously not associated with this syndrome. Both missense variants and putative loss-of-function variants resulted in slow growth in patient-derived cell lines. KMT5B homozygous knockout mice were smaller in size than their wild-type littermates but did not have significantly smaller brains, suggesting relative macrocephaly, also noted as a prominent clinical feature. RNA sequencing of patient lymphoblasts and Kmt5b haploinsufficient mouse brains identified differentially expressed pathways associated with nervous system development and function including axon guidance signaling. Overall, we identified additional pathogenic variants and clinical features in KMT5B-related neurodevelopmental disorder and provide insights into the molecular mechanisms of the disorder using multiple model systems.

A mutational hotspot in AMOTL1 defines a new syndrome of orofacial clefting, cardiac anomalies, and tall stature.

AMOTL1 encodes angiomotin-like protein 1, an actin-binding protein that regulates cell polarity, adhesion, and migration. The role of AMOTL1 in human disease is equivocal. We report a large cohort of individuals harboring heterozygous AMOTL1 variants and define a core phenotype of orofacial clefting, congenital heart disease, tall stature, auricular anomalies, and gastrointestinal manifestations in individuals with variants in AMOTL1 affecting amino acids 157-161, a functionally undefined but highly conserved region. Three individuals with AMOTL1 variants outside this region are also described who had variable presentations with orofacial clefting and multi-organ disease. Our case cohort suggests that heterozygous missense variants in AMOTL1, most commonly affecting amino acid residues 157-161, define a new orofacial clefting syndrome, and indicates an important functional role for this undefined region.

GATA1 controls numbers of hematopoietic progenitors and their response to autoimmune neuroinflammation.

GATA-binding factor 1 (GATA1) is a transcription factor that governs the development and function of multiple hematopoietic cell lineages. GATA1 is expressed in hematopoietic stem and progenitor cells (HSPCs) and is essential for erythroid lineage commitment; however, whether it plays a role in hematopoietic stem cell (HSC) biology and the development of myeloid cells, and what that role might be, remains unclear. We initially set out to test the role of eosinophils in experimental autoimmune encephalomyelitis (EAE), a model of central nervous system autoimmunity, using mice lacking a double GATA-site (ΔdblGATA), which lacks eosinophils due to the deletion of the dblGATA enhancer to Gata1, which alters its expression. ΔdblGATA mice were resistant to EAE, but not because of a lack of eosinophils, suggesting that these mice have an additional defect. ΔdblGATA mice with EAE had fewer inflammatory myeloid cells than the control mice, suggesting that resistance to EAE is caused by a defect in myeloid cells. Naïve ΔdblGATA mice also showed reduced frequency of CD11b+ myeloid cells in the blood, indicating a defect in myeloid cell production. Examination of HSPCs revealed fewer HSCs and myeloid cell progenitors in the ΔdblGATA bone marrow (BM), and competitive BM chimera experiments showed a reduced capacity of the ΔdblGATA BM to reconstitute immune cells, suggesting that reduced numbers of ΔdblGATA HSPCs cause a functional deficit during inflammation. Taken together, our data show that GATA1 regulates the number of HSPCs and that reduced GATA1 expression due to dblGATA deletion results in a diminished immune response following the inflammatory challenge.

TBX6 as a cause of a combined skeletal-kidney dysplasia syndrome.

TBX6 encodes transcription-factor box 6, a transcription factor critical to paraxial mesoderm segmentation and somitogenesis during embryonic development. TBX6 haploinsufficiency is believed to drive the skeletal and kidney phenotypes associated with the 16p11.2 deletion syndrome. Heterozygous and biallelic variants in TBX6 are associated with vertebral and rib malformations (TBX6-associated congenital scoliosis) and spondylocostal dysostosis, and heterozygous TBX6 variants are associated with increased risk of genitourinary tract malformations. Combined skeletal and kidney phenotypes in individuals harboring heterozygous or biallelic TBX6 variants are rare. Here, we present seven individuals with vertebral and rib malformations and structural kidney differences associated with heterozygous TBX6 gene deletion in trans with a hypomorphic TBX6 allele or biallelic TBX6 variants. Our case series highlights the association between TBX6 and both skeletal and kidney disease.

Identification of Novel Loci Shared by Juvenile Idiopathic Arthritis Subtypes Through Integrative Genetic Analysis.

OBJECTIVE: Juvenile idiopathic arthritis (JIA) is the most common chronic immune-mediated joint disease among children and encompasses a heterogeneous group of immune-mediated joint disorders classified into 7 subtypes according to clinical presentation. However, phenotype overlap and biologic evidence suggest a shared mechanistic basis between subtypes. This study was undertaken to systematically investigate shared genetic underpinnings of JIA subtypes. METHODS: We performed a heterogeneity-sensitive genome-wide association study encompassing a total of 1,245 JIA cases (classified into 7 subtypes) and 9,250 controls, followed by fine-mapping of candidate causal variants at each genome-wide significant locus, functional annotation, and pathway and network analysis. We further identified candidate drug targets and drug repurposing opportunities by in silico analyses. RESULTS: In addition to the major histocompatibility complex locus, we identified 15 genome-wide significant loci shared between at least 2 JIA subtypes, including 10 novel loci. Functional annotation indicated that candidate genes at these loci were expressed in diverse immune cell types. CONCLUSION: This study identified novel genetic loci shared by JIA subtypes. Our findings identified candidate mechanisms underlying JIA subtypes and candidate targets with drug repurposing opportunities for JIA treatment.

Lymphocytes Utilize Somatic Mutations, Epigenetic Silencing, and the Proteasome to Escape Truncated WASP Expression.

Wiskott-Aldrich Syndrome Protein (WASP) deficiency causes Wiskott-Aldrich Syndrome (WAS), a sex-linked disorder characterized by combined immunodeficiency, microthrombocytopenia, and eczema. Like WASP-deficient humans, WASP-deficient mice produce normal numbers of functionally defective T cells. Here, we report a WAS patient with a novel germline frameshifting WAS mutation encoding a truncated form of WASP lacking the C-terminal cofilin homology (C) and the acidic region (A) domains (WASPΔCA). Although stably overexpressed in embryonic kidney cell lines, WASPΔCA was undetectable in circulating patient leukocytes. Deep sequencing, transcript profiling, and protein degradation analyses demonstrated patient lymphocytes employ an array of genetic, epigenetic, and proteasomal strategies to avoid expressing WASPΔCA.

Genetic analysis for type 1 diabetes genes in juvenile dermatomyositis unveils genetic disease overlap.

OBJECTIVES: JDM is a serious autoimmune and complex genetic disease. Another autoimmune genetic disease, type 1 diabetes (T1D), has been observed for significantly increased prevalence in families with JDM, while increased JDM risk has also been observed in T1D cases. This study aimed to study whether these two autoimmune diseases, JDM and T1D, share common genetic susceptibility. METHODS: From 169 JDM families, 121 unrelated cases with European ancestry (EA) were identified by genome-wide genotyping, principal component analysis and identical-by-descent (IBD) analysis. T1D genetic risk score (GRS) were calculated in these cases and were compared with 361 EA T1D cases and 1943 non-diabetes EA controls. A total of 113 cases of the 121 unrelated European cases were sequenced by whole exome sequencing. RESULTS: We observed increased T1D GRS in JDM cases (P = 9.42E-05). Using whole exome sequencing, we uncovered the T1D genes, phospholipase B1, cystic fibrosis transmembrane conductance regulator, tyrosine hydroxylase, CD6 molecule, perforin 1 and dynein axonemal heavy chain 2, potentially associated with JDM by the burden test of rare functional coding variants. CONCLUSION: Novel mechanisms of JDM related to these T1D genes are suggested by this study, which may imply novel therapeutic targets for JDM and warrant further study.

A novel MBTPS2 variant associated with BRESHECK syndrome impairs sterol-regulated transcription and the endoplasmic reticulum stress response.

Ichthyosis follicularis, atrichia, and photophobia syndrome (IFAP syndrome) is a rare, X-linked disorder caused by pathogenic variants in membrane-bound transcription factor protease, site 2 (MBTPS2). Pathogenic MBTPS2 variants also cause BRESHECK syndrome, characterized by the IFAP triad plus intellectual disability and multiple congenital anomalies. Here we present a patient with ichthyosis, sparse hair, pulmonic stenosis, kidney dysplasia, hypospadias, growth failure, thrombocytopenia, anemia, bone marrow fibrosis, and chronic diarrhea found by research-based exome sequencing to harbor a novel, maternally inherited MBTPS2 missense variant (c.766 G>A; (p.Val256Leu)). In vitro modeling supports variant pathogenicity, with impaired cell growth in cholesterol-depleted media, attenuated activation of the sterol regulatory element-binding protein pathway, and failure to activate the endoplasmic reticulum stress response pathway. Our case expands both the genetic and phenotypic spectrum of BRESHECK syndrome to include a novel MBTPS2 variant and cytopenias, bone marrow fibrosis, and chronic diarrhea.

Expanding the genetic landscape of oral-facial-digital syndrome with two novel genes.

Oral-facial-digital syndromes (OFDS) are a heterogeneous and rare group of Mendelian disorders characterized by developmental abnormalities of the oral cavity, face, and digits caused by dysfunction of the primary cilium, a mechanosensory organelle that exists atop most cell types that facilitates organ patterning and growth. OFDS is inherited both in an X-linked dominant, X-linked recessive, and autosomal recessive manner. Importantly, though many of the causal genes for OFDS have been identified, up to 40% of OFD syndromes are of unknown genetic basis. Here we describe three children with classical presentations of OFDS including lingual hamartomas, polydactyly, and characteristic facial features found by exome sequencing to harbor variants in causal genes not previously associated with OFDS. We describe a female with hypothalamic hamartoma, urogenital sinus, polysyndactyly, and multiple lingual hamartomas consistent with OFDVI with biallelic pathogenic variants in CEP164, a gene associated with ciliopathy-spectrum disease, but never before with OFDS. We additionally describe two unrelated probands with postaxial polydactyly, multiple lingual hamartomas, and dysmorphic features both found to be homozygous for an identical TOPORS missense variant, c.29 C>A; (p.Pro10Gln). Heterozygous TOPORS pathogenic gene variants are associated with autosomal dominant retinitis pigmentosa, but never before with syndromic ciliopathy. Of note, both probands are of Dominican ancestry, suggesting a possible founder allele.

A new syndrome of moyamoya disease, kidney dysplasia, aminotransferase elevation, and skin disease associated with de novo variants in RNF213.

Ring-finger protein 213 (RNF213) encodes a protein of unknown function believed to play a role in cellular metabolism and angiogenesis. Gene variants are associated with susceptibility to moyamoya disease. Here, we describe two children with moyamoya disease who also demonstrated kidney disease, elevated aminotransferases, and recurrent skin lesions found by exome sequencing to have de novo missense variants in RNF213. These cases highlight the ability of RNF213 to cause Mendelian moyamoya disease in addition to acting as a genetic susceptibility locus. The cases also suggest a new, multi-organ RNF213-spectrum disease characterized by liver, skin, and kidney pathology in addition to severe moyamoya disease caused by heterozygous, de novo C-terminal RNF213 missense variants.

Evaluating sequence data quality from the Swift Accel-Amplicon CFTR Panel.

Cystic fibrosis (CF) is one of the most common genetic diseases worldwide with high carrier frequencies across different ethnicities. Next generation sequencing of the cystic fibrosis transmembrane conductance regulator (CFTR) gene has proven to be an effective screening tool to determine carrier status with high detection rates. Here, we evaluate the performance of the Swift Biosciences Accel-Amplicon CFTR Capture Panel using CFTR-positive DNA samples. This assay is a one-day protocol that allows for one-tube reaction of 87 amplicons that span all coding regions, 5' and 3'UTR, as well as four intronic regions. In this study, we provide the FASTQ, BAM, and VCF files on seven unique CFTR-positive samples and one normal control sample (14 samples processed including repeated samples). This method generated sequencing data with high coverage and near 100% on-target reads. We found that coverage depth was correlated with the GC content of each exon. This dataset is instrumental for clinical laboratories that are evaluating this technology as part of their carrier screening program.