Multitrait genome-wide analyses identify new susceptibility loci and candidate drugs to primary sclerosing cholangitis.

Primary sclerosing cholangitis (PSC) is a rare autoimmune bile duct disease that is strongly associated with immune-mediated disorders. In this study, we implemented multitrait joint analyses to genome-wide association summary statistics of PSC and numerous clinical and epidemiological traits to estimate the genetic contribution of each trait and genetic correlations between traits and to identify new lead PSC risk-associated loci. We identified seven new loci that have not been previously reported and one new independent lead variant in the previously reported locus. Functional annotation and fine-mapping nominated several potential susceptibility genes such as MANBA and IRF5. Network-based in silico drug efficacy screening provided candidate agents for further study of pharmacological effect in PSC.

Self-labelled encoder-decoder (SLED) for multi-echo gradient echo-based myelin water imaging.

PURPOSE: Reconstruction of high quality myelin water imaging (MWI) maps is challenging, particularly for data acquired using multi-echo gradient echo (mGRE) sequences. A non-linear least squares fitting (NLLS) approach has often been applied for MWI. However, this approach may produce maps with limited detail and, in some cases, sub-optimal signal to noise ratio (SNR), due to the nature of the voxel-wise fitting. In this study, we developed a novel, unsupervised learning method called self-labelled encoder-decoder (SLED) to improve gradient echo-based MWI data fitting. METHODS: Ultra-high resolution, MWI data was collected from five mouse brains with variable levels of myelination, using a mGRE sequence. Imaging data was acquired using a 7T preclinical MRI system. A self-labelled, encoder-decoder network was implemented in TensorFlow for calculation of myelin water fraction (MWF) based on the mGRE signal decay. A simulated MWI phantom was also created to evaluate the performance of MWF estimation. RESULTS: Compared to NLLS, SLED demonstrated improved MWF estimation, in terms of both stability and accuracy in phantom tests. In addition, SLED produced less noisy MWF maps from high resolution MR microscopy images of mouse brain tissue. It specifically resulted in lower noise amplification for all mouse genotypes that were imaged and yielded mean MWF values in white matter ROIs that were highly correlated with those derived from standard NLLS fitting. Lastly, SLED also exhibited higher tolerance to low SNR data. CONCLUSION: Due to its unsupervised and self-labeling nature, SLED offers a unique alternative to analyze gradient echo-based MWI data, providing accurate and stable MWF estimations.

Novel insights into mouse models of ectopic proplatelet release.

Mature bone marrow (BM) megakaryocytes (MKs) produce platelets by extending proplatelets into sinusoidal blood vessels. Defects in this process can lead to thrombocytopenia and increased risk of bleeding. Mice lacking the actin-regulatory proteins Profilin 1 (PFN1), Wiskott-Aldrich Syndrome protein (WASp), Actin Related Protein 2/3 complex (Arp2/3), or adhesion and degranulation-promoting adapter protein (ADAP) display thrombocytopenia and ectopic release of (pro)platelet-like particles into the BM compartment, pointing to an important axis of actin-mediated directional proplatelet formation. The mechanism underlying ectopic release in these mice is still not completely understood. However, we hypothesized that similar functional defects account for this observation. We analyzed WASp-, ADAP-, PFN1-, and ARPC2-knockout mice to determine the role of actin reorganization and integrin activation in directional proplatelet formation. ADAP-, ARPC2-, and PFN1-deficient MKs displayed reduced adhesion to collagen, defective F-actin organization, and diminished ß1-integrin activation. WASp-deficient MKs showed the strongest reduction in the adhesion assay of collagen and altered F-actin organization with reduced podosome formation. Our results indicate that ADAP, PFN1, WASp, and ARP2/3 are part of the same pathway that regulates polarization processes in MKs and directional proplatelet formation into BM sinusoids.

Genome-wide analysis provides genetic evidence that ACE2 influences COVID-19 risk and yields risk scores associated with severe disease.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters human host cells via angiotensin-converting enzyme 2 (ACE2) and causes coronavirus disease 2019 (COVID-19). Here, through a genome-wide association study, we identify a variant (rs190509934, minor allele frequency 0.2-2%) that downregulates ACE2 expression by 37% (P = 2.7 × 10-8) and reduces the risk of SARS-CoV-2 infection by 40% (odds ratio = 0.60, P = 4.5 × 10-13), providing human genetic evidence that ACE2 expression levels influence COVID-19 risk. We also replicate the associations of six previously reported risk variants, of which four were further associated with worse outcomes in individuals infected with the virus (in/near LZTFL1, MHC, DPP9 and IFNAR2). Lastly, we show that common variants define a risk score that is strongly associated with severe disease among cases and modestly improves the prediction of disease severity relative to demographic and clinical factors alone.

An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs.

BACKGROUNDS & AIMS: Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intrahepatic bile ducts eventually leads to cirrhosis. Many patients have inadequate response to licensed medications, motivating the search for novel therapies. Previous genome-wide association studies (GWAS) and meta-analyses (GWMA) of PBC have identified numerous risk loci for this condition, providing insight into its aetiology. We undertook the largest GWMA of PBC to date, aiming to identify additional risk loci and prioritise candidate genes for in silico drug efficacy screening. METHODS: We combined new and existing genotype data for 10,516 cases and 20,772 controls from 5 European and 2 East Asian cohorts. RESULTS: We identified 56 genome-wide significant loci (20 novel) including 46 in European, 13 in Asian, and 41 in combined cohorts; and a 57th genome-wide significant locus (also novel) in conditional analysis of the European cohorts. Candidate genes at newly identified loci include FCRL3, INAVA, PRDM1, IRF7, CCR6, CD226, and IL12RB1, which each play key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, JAK-STAT signalling, and differentiation of T helper (TH)1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some of which are well-established in the treatment of other autoimmune disorders. CONCLUSIONS: This study has identified additional risk loci for PBC, provided a hierarchy of agents that could be trialled in this condition, and emphasised the value of genetic and genomic approaches to drug discovery in complex disorders. LAY SUMMARY: Primary biliary cholangitis (PBC) is a chronic liver disease that eventually leads to cirrhosis. In this study, we analysed genetic information from 10,516 people with PBC and 20,772 healthy individuals recruited in Canada, China, Italy, Japan, the UK, or the USA. We identified several genetic regions associated with PBC. Each of these regions contains several genes. For each region, we used diverse sources of evidence to help us choose the gene most likely to be involved in causing PBC. We used these 'candidate genes' to help us identify medications that are currently used for treatment of other conditions, which might also be useful for treatment of PBC.

WAVE2 suppresses mTOR activation to maintain T cell homeostasis and prevent autoimmunity.

Cytoskeletal regulatory protein dysfunction has been etiologically linked to inherited diseases associated with immunodeficiency and autoimmunity, but the mechanisms involved are incompletely understood. Here, we show that conditional Wave2 ablation in T cells causes severe autoimmunity associated with increased mammalian target of rapamycin (mTOR) activation and metabolic reprogramming that engender spontaneous activation and accelerated differentiation of peripheral T cells. These mice also manifest diminished antigen-specific T cell responses associated with increased inhibitory receptor expression, dysregulated mitochondrial function, and reduced cell survival upon activation. Mechanistically, WAVE2 directly bound mTOR and inhibited its activation by impeding mTOR interactions with RAPTOR (regulatory-associated protein of mTOR) and RICTOR (rapamycin-insensitive companion of mTOR). Both the T cell defects and immunodysregulatory disease were ameliorated by pharmacological mTOR inhibitors. Thus, WAVE2 restraint of mTOR activation is an absolute requirement for maintaining the T cell homeostasis supporting adaptive immune responses and preventing autoimmunity.

X Chromosome Contribution to the Genetic Architecture of Primary Biliary Cholangitis.

BACKGROUND & AIMS: Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease. METHODS: We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals). RESULTS: Single-marker association analyses found approximately 100 loci displaying P < 5 × 10-4, with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 × 10-6; odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 × 10-8), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 × 10-9; OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively). CONCLUSIONS: This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus.

Transcriptional regulators of the Golli/myelin basic protein locus integrate additive and stealth activities.

Myelin is composed of plasma membrane spirally wrapped around axons and compacted into dense sheaths by myelin-associated proteins. Myelin is elaborated by neuroepithelial derived oligodendrocytes in the central nervous system (CNS) and by neural crest derived Schwann cells in the peripheral nervous system (PNS). While some myelin proteins accumulate in only one lineage, myelin basic protein (Mbp) is expressed in both. Overlapping the Mbp gene is Golli, a transcriptional unit that is expressed widely both within and beyond the nervous system. A super-enhancer domain within the Golli/Mbp locus contains multiple enhancers shown previously to drive reporter construct expression specifically in oligodendrocytes or Schwann cells. In order to determine the contribution of each enhancer to the Golli/Mbp expression program, and to reveal if functional interactions occur among them, we derived mouse lines in which they were deleted, either singly or in different combinations, and relative mRNA accumulation was measured at key stages of early development and at maturity. Although super-enhancers have been shown previously to facilitate interaction among their component enhancers, the enhancers investigated here demonstrated largely additive relationships. However, enhancers demonstrating autonomous activity strictly in one lineage, when missing, were found to significantly reduce output in the other, thus revealing cryptic "stealth" activity. Further, in the absence of a key oligodendrocyte enhancer, Golli accumulation was markedly and uniformly attenuated in all cell types investigated. Our observations suggest a model in which enhancer-mediated DNA-looping and potential super-enhancer properties underlie Golli/Mbp regulatory organization.

Comprehensive Profiling of the Rheumatoid Arthritis Antibody Repertoire.

OBJECTIVE: Autoantibodies against citrullinated proteins are found in 64-89% of rheumatoid arthritis (RA) patients, with 88-99% specificity. This study was undertaken to create an unbiased, comprehensive profile of serum antibodies against the human proteome, including the citrullinome and the homocitrullinome, in RA patients, using a high-density peptide array. METHODS: Our high-density peptide array, consisting of >4.6 million peptides, contained the entire annotated human proteome. The 20,246 proteins were represented as overlapping 16-mer peptides. In addition to native peptides, citrullinated and homocitrullinated peptides were included, as substitutions for arginine and lysine, and provided a comprehensive screen against all possible epitopes. Twenty-six serum samples (from 8 controls and 18 RA patients) were profiled on the high-density peptide array. Using RA-specific epitopes, we constructed an 8-epitope diagnostic biomarker on a Gyrolab xPlore instrument with a cohort of 92 serum samples (from 29 controls and 63 RA patients). The diagnostic biomarker was further validated with an independent cohort of 181 serum samples (from 54 controls and 127 RA patients). RESULTS: In the initial cohort the diagnostic performance of the 8-epitope biomarker yielded 96.6% specificity and 92.1% sensitivity. The overall diagnostic performance in the validation cohort was 94.4% specificity and 85% sensitivity. In both cohorts, the performance of the 8-epitope diagnostic biomarker compared favorably against the Abnova cyclic citrullinated peptide 2 (CCP2) assay. Using data from the peptide array, we identified novel RA-specific epitopes and formed the basis of a new RA diagnostic assay. CONCLUSION: Comprehensive antibody profiling using a high-density peptide array not only identified novel RA-specific epitopes but also allowed us to construct a novel diagnostic biomarker that is as specific as and more sensitive than the Abnova CCP2 assay.

A modifier in the 129S2/SvPasCrl genome is responsible for the viability of Notch1[12f/12f] mice.

BACKGROUND: Mouse NOTCH1 carries a highly conserved O-fucose glycan at Thr466 in epidermal growth factor-like repeat 12 (EGF12) of the extracellular domain. O-Fucose at this site has been shown by X-ray crystallography to be recognized by both DLL4 and JAG1 Notch ligands. We previously showed that a Notch1 Thr466Ala mutant exhibits very little ligand-induced NOTCH1 signaling in a reporter assay, whereas a Thr466Ser mutation enables the transfer of O-fucose and reverts the NOTCH1 signaling defect. We subsequently generated a mutant mouse with the Thr466Ala mutation termed Notch1[12f](Notch1tm2Pst). Surprisingly, homozygous Notch1[12f/12f] mutants on a mixed background were viable and fertile. RESULTS: We now report that after backcrossing to C57BL/6 J mice for 11-15 generations, few homozygous Notch1[12f/12f] embryos were born. Timed mating showed that embryonic lethality occurred by embryonic day (E) ~E11.5, somewhat delayed compared to mice lacking Notch1 or Pofut1 (the O-fucosyltransferase that adds O-fucose to Notch receptors), which die at ~E9.5. The phenotype of C57BL/6 J Notch1[12f/12f] embryos was milder than mutants affected by loss of a canonical Notch pathway member, but disorganized vasculogenesis in the yolk sac, delayed somitogenesis and development were characteristic. In situ hybridization of Notch target genes Uncx4.1 and Dll3 or western blot analysis of NOTCH1 cleavage did not reveal significant differences at E9.5. However, qRT-PCR of head cDNA showed increased expression of Dll3, Uncx4.1 and Notch1 in E9.5 Notch1[12f/12f] embryos. Sequencing of cDNA from Notch1[12f/12f] embryo heads and Southern analysis showed that the Notch1[12f] locus was intact following backcrossing. We therefore looked for evidence of modifying gene(s) by crossing C57BL/6 J Notch1 [12f/+] mice to 129S2/SvPasCrl mice. Intercrosses of the F1 progeny gave viable F2 Notch1[12f/12f] mice. CONCLUSION: We conclude that the 129S2/SvPasCrl genome contains a dominant modifying gene that rescues the functions of NOTCH1[12f] in signaling. Identification of the modifying gene has the potential to illuminate novel factor(s) that promote Notch signaling when an O-fucose glycan is absent from EGF12 of NOTCH1.

Variation at DENND1B and Asthma on the Island of Tristan da Cunha.

A high prevalence of asthma has been documented among the inhabitants of Tristan da Cunha, an isolated island in the South Atlantic. The population derives from just 28 founders. We performed lung function testing, including methacholine inhalation challenge, allergen skin prick testing, and collected DNA from essentially all of the current island population (269 individuals), and genotyped a panel of 43 single-nucleotide polymorphisms (SNPs) reported as associated with asthma and atopy. We carried out a mixed-model association analysis using the known pedigree. There were 96 individuals diagnosed as asthmatic (36%), and heritability estimates were similar to those from nonisolated population samples (multifactorial threshold model, h2 = 48%). The first component from a genetic principal components analysis using the entire SNP panel was nonlinearly associated with asthma, with the maximum risk to those intermediate to reference (Human Genome Diversity Project) European and African samples means. The single most strongly associated SNP was rs2786098 (p = 5.5 × 10-5), known to regulate the gene DENND1B. This explained approximately one-third of the trait heritability, with an allelic odds ratio for the A allele of 2.6. Among A/A carriers, 10 out of 12 individuals were asthmatic. The rs2786098*A variant was initially reported to decrease the risk of childhood (atopic) asthma in European but slightly increase the risk in African-descended populations, and does significantly alter Th2 cell function. Despite an absence of overall association with this variant in recent asthma genome wide association studies meta-analyses, an effect may exist on the particular genetic background of the Tristan da Cunha population.

Genetic interaction analysis among oncogenesis-related genes revealed novel genes and networks in lung cancer development.

The development of cancer is driven by the accumulation of many oncogenesis-related genetic alterations and tumorigenesis is triggered by complex networks of involved genes rather than independent actions. To explore the epistasis existing among oncogenesis-related genes in lung cancer development, we conducted pairwise genetic interaction analyses among 35,031 SNPs from 2027 oncogenesis-related genes. The genotypes from three independent genome-wide association studies including a total of 24,037 lung cancer patients and 20,401 healthy controls with Caucasian ancestry were analyzed in the study. Using a two-stage study design including discovery and replication studies, and stringent Bonferroni correction for multiple statistical analysis, we identified significant genetic interactions between SNPs in RGL1:RAD51B (OR=0.44, p value=3.27x10-11 in overall lung cancer and OR=0.41, p value=9.71x10-11 in non-small cell lung cancer), SYNE1:RNF43 (OR=0.73, p value=1.01x10-12 in adenocarcinoma) and FHIT:TSPAN8 (OR=1.82, p value=7.62x10-11 in squamous cell carcinoma) in our analysis. None of these genes have been identified from previous main effect association studies in lung cancer. Further eQTL gene expression analysis in lung tissues provided information supporting the functional role of the identified epistasis in lung tumorigenesis. Gene set enrichment analysis revealed potential pathways and gene networks underlying molecular mechanisms in overall lung cancer as well as histology subtypes development. Our results provide evidence that genetic interactions between oncogenesis-related genes play an important role in lung tumorigenesis and epistasis analysis, combined with functional annotation, provides a valuable tool for uncovering functional novel susceptibility genes that contribute to lung cancer development by interacting with other modifier genes.

Genetic Variants in Antineutrophil Cytoplasmic Antibody-Associated Vasculitis: A Bayesian Approach and Systematic Review.

A number of genome-wide association studies (GWASs) and meta-analyses of genetic variants have been performed in antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis. We reinterpreted previous studies using false-positive report probability (FPRP) and Bayesian false discovery probability (BFDP). This study searched publications in PubMed and Excerpta Medica Database (EMBASE) up to February 2018. Identification of noteworthy associations were analyzed using FPRP and BFDP, and data (i.e., odds ratio (OR), 95% confidence interval (CI), p-value) related to significant associations were separately extracted. Using filtered gene variants, gene ontology (GO) enrichment analysis and protein⁻protein interaction (PPI) networks were performed. Overall, 241 articles were identified, and 7 were selected for analysis. Single nucleotide polymorphisms (SNPs) discovered by GWASs were shown to be noteworthy, whereas only 27% of significant results from meta-analyses of observational studies were noteworthy. Eighty-five percent of SNPs with borderline p-values (5.0 × 10-8 < p < 0.05) in GWASs were found to be noteworthy. No overlapping SNPs were found between PR3-ANCA and MPO-ANCA vasculitis. GO analysis revealed immune-related GO terms, including "antigen processing and presentation of peptide or polysaccharide antigen via major histocompatibility complex (MHC) class II", "interferon-gamma-mediated (IFN-γ) signaling pathway". By using FPRP and BFDP, network analysis of noteworthy genetic variants discovered genetic risk factors associated with the IFN-γ pathway as novel mechanisms potentially implicated in the complex pathogenesis of ANCA-associated vasculitis.

TIE: A Method to Electroporate Long DNA Templates into Preimplantation Embryos for CRISPR-Cas9 Gene Editing.

Precise genome editing using CRISPR typically requires delivery of guide RNAs, Cas9 endonuclease, and DNA repair templates. Both microinjection and electroporation effectively deliver these components into mouse zygotes provided the DNA template is an oligonucleotide of only a few hundred base pairs. However, electroporation completely fails with longer double-stranded DNAs leaving microinjection as the only delivery option. Here, we overcome this limitation by first injecting all CRISPR components, including long plasmid-sized DNA templates, into the sub-zona pellucida space. There they are retained, supporting subsequent electroporation. We show that this simple and well-tolerated method achieves intracellular reagent concentrations sufficient to effect precise gene edits.

Prevalence and Clinical Implication of Double Mutations in Hypertrophic Cardiomyopathy: Revisiting the Gene-Dose Effect.

BACKGROUND: Available data suggests that double mutations in patients with hypertrophic cardiomyopathy are not rare and are associated with a more severe phenotype. Most of this data, however, is based on noncontemporary variant classification. METHODS AND RESULTS: Clinical data of all hypertrophic cardiomyopathy patients with 2 rare genetic variants were retrospectively reviewed and compared with a group of patients with a single disease-causing variant. Furthermore, a literature search was performed for all studies with information on prevalence and outcome of patients with double mutations. Classification of genetic variants was reanalyzed according to current guidelines. In our cohort (n=1411), 9% of gene-positive patients had 2 rare variants in sarcomeric genes but only in 1 case (0.4%) were both variants classified as pathogenic. Patients with 2 rare variants had a trend toward younger age at presentation when compared with patients with a single mutation. All other clinical variables were similar. In data pooled from cohort studies in the literature, 8% of gene-positive patients were published to have double mutations. However, after reanalysis of reported variants, this prevalence diminished to 0.4%. All patients with 2 radical mutations in MYBPC3 in the literature had severe disease with death or heart transplant during the first year of life. Data on other specific genotype-phenotype correlations were scarce. CONCLUSIONS: Double mutations in patients with hypertrophic cardiomyopathy are much less common than previously estimated. With the exception of double radical MYBPC3 mutations, there is little data to guide clinical decision making in cases with double mutations.

Identification of Functional and Expression Polymorphisms Associated With Risk for Antineutrophil Cytoplasmic Autoantibody-Associated Vasculitis.

OBJECTIVE: To identify risk alleles relevant to the causal and biologic mechanisms of antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV). METHODS: A genome-wide association study and subsequent replication study were conducted in a total cohort of 1,986 cases of AAV (patients with granulomatosis with polyangiitis [Wegener's] [GPA] or microscopic polyangiitis [MPA]) and 4,723 healthy controls. Meta-analysis of these data sets and functional annotation of identified risk loci were performed, and candidate disease variants with unknown functional effects were investigated for their impact on gene expression and/or protein function. RESULTS: Among the genome-wide significant associations identified, the largest effect on risk of AAV came from the single-nucleotide polymorphism variants rs141530233 and rs1042169 at the HLA-DPB1 locus (odds ratio [OR] 2.99 and OR 2.82, respectively) which, together with a third variant, rs386699872, constitute a triallelic risk haplotype associated with reduced expression of the HLA-DPB1 gene and HLA-DP protein in B cells and monocytes and with increased frequency of complementary proteinase 3 (PR3)-reactive T cells relative to that in carriers of the protective haplotype. Significant associations were also observed at the SERPINA1 and PTPN22 loci, the peak signals arising from functionally relevant missense variants, and at PRTN3, in which the top-scoring variant correlated with increased PRTN3 expression in neutrophils. Effects of individual loci on AAV risk differed between patients with GPA and those with MPA or between patients with PR3-ANCAs and those with myeloperoxidase-ANCAs, but the collective population attributable fraction for these variants was substantive, at 77%. CONCLUSION: This study reveals the association of susceptibility to GPA and MPA with functional gene variants that explain much of the genetic etiology of AAV, could influence and possibly be predictors of the clinical presentation, and appear to alter immune cell proteins and responses likely to be key factors in the pathogenesis of AAV.