Ancestry-based stratified analysis of Immunochip data identifies novel associations with celiac disease

m6A Methylated Long Noncoding RNA LOC339803 Regulates Intestinal Inflammatory Response

Cytokine mediated sustained inflammation increases the risk to develop different complex chronic inflammatory diseases, but the implicated mechanisms remain unclear. Increasing evidence shows that long noncoding RNAs (lncRNAs) play key roles in the pathogenesis of inflammatory disorders, while inflammation associated variants are described to affect their function or essential RNA modifications as N6-methyladenosine (m6A) methylation, increasing predisposition to inflammatory diseases. Here, the functional implication of the intestinal inflammation associated lncRNA LOC339803 in the production of cytokines by intestinal epithelial cells is described. Allele-specific m6A methylation is found to affect YTHDC1 mediated protein binding affinity. LOC339803-YTHDC1 interaction dictates chromatin localization of LOC339803 ultimately inducing the expression of NFκB mediated proinflammatory cytokines and contributing to the development of intestinal inflammation. These findings are confirmed using human intestinal biopsy samples from different intestinal inflammatory conditions and controls. Additionally, it is demonstrated that LOC339803 targeting can be a useful strategy for the amelioration of intestinal inflammation in vitro and ex vivo. Overall, the results support the importance of the methylated LOC339803 lncRNA as a mediator of intestinal inflammation, explaining genetic susceptibility and presenting this lncRNA as a potential novel therapeutic target for the treatment of inflammatory intestinal disorders.

Sex bias in celiac disease: XWAS and monocyte eQTLs in women identify TMEM187 as a functional candidate gene

BACKGROUND

Celiac disease (CeD) is an immune-mediated disorder that develops in genetically predisposed individuals upon gluten consumption. HLA risk alleles explain 40% of the genetic component of CeD, so there have been continuing efforts to uncover non-HLA loci that can explain the remaining heritability. As in most autoimmune disorders, the prevalence of CeD is significantly higher in women. Here, we investigated the possible involvement of the X chromosome on the sex bias of CeD.

METHODS

We performed a X chromosome-wide association study (XWAS) and a gene-based association study in women from the CeD Immunochip (7062 cases, 5446 controls). We also constructed a database of X chromosome cis-expression quantitative trait loci (eQTLs) in monocytes from unstimulated (n = 226) and lipopolysaccharide (LPS)-stimulated (n = 130) female donors and performed a Summary-data-based MR (SMR) analysis to integrate XWAS and eQTL information. We interrogated the expression of the potentially causal gene (TMEM187) in peripheral blood mononuclear cells (PBMCs) from celiac patients at onset, on a gluten-free diet, potential celiac patients and non-celiac controls.

RESULTS

The XWAS and gene-based analyses identified 13 SNPs and 25 genes, respectively, 22 of which had not been previously associated with CeD. The X chromosome cis-eQTL analysis found 18 genes with at least one cis-eQTL in naïve female monocytes and 8 genes in LPS-stimulated female monocytes, 2 of which were common to both situations and 6 were unique to LPS stimulation. SMR identified a potentially causal association of TMEM187 expression in naïve monocytes with CeD in women, regulated by CeD-associated, eQTL-SNPs rs7350355 and rs5945386. The CeD-risk alleles were correlated with lower TMEM187 expression. These results were replicated using eQTLs from LPS-stimulated monocytes. We observed higher levels of TMEM187 expression in PBMCs from female CeD patients at onset compared to female non-celiac controls, but not in male CeD individuals.

CONCLUSION

Using X chromosome genotypes and gene expression data from female monocytes, SMR has identified TMEM187 as a potentially causal candidate in CeD. Further studies are needed to understand the implication of the X chromosome in the higher prevalence of CeD in women.

Effect of SARS-CoV-2 infection on asthma patients

Background

SARS-CoV-2 causes coronavirus disease 2019 (COVID-19), a new coronavirus pneumonia, and containing such an international pandemic catastrophe remains exceedingly difficult. Asthma is a severe chronic inflammatory airway disease that is becoming more common around the world. However, the link between asthma and COVID-19 remains unknown. Through bioinformatics analysis, this study attempted to understand the molecular pathways and discover potential medicines for treating COVID-19 and asthma.

Methods

To investigate the relationship between SARS-CoV-2 and asthma patients, a transcriptome analysis was used to discover shared pathways and molecular signatures in asthma and COVID-19. Here, two RNA-seq data (GSE147507 and GSE74986) from the Gene Expression Omnibus were used to detect differentially expressed genes (DEGs) in asthma and COVID-19 patients to find the shared pathways and the potential drug candidates.

Results

There were 66 DEGs in all that were classified as common DEGs. Using a protein-protein interaction (PPI) network created using various bioinformatics techniques, five hub genes were found. We found that asthma has some shared links with the progression of COVID-19. Additionally, protein-drug interactions with common DEGs were also identified in the datasets.

Conclusion

We investigated possible links between COVID-19 and asthma using bioinformatics databases, which might be useful in treating COVID-19 patients. More studies on populations affected by these diseases are needed to elucidate the molecular mechanism behind their association.

Genetic variants at the chromosomal region 2q21.3 underlying inhibitor development in patients with severe haemophilia A

Introduction

Inhibitor development affects about 30% of patients with severe haemophilia A (HA) and results from different environmental and genetic risk factors. Previously, we identified the missense variant rs3754689 in the LCT gene linked with this predisposition. Since rs3754689 variant is benign and is located in a conserved haplotype region, we hypothesized that the association signal captured by this variant is located in coinherited, neighbouring genes.

Aim

To identify novel genetic risk factors associated with inhibitor development in coding regions of R3HDM1, UBXN4, CXCR4, MCM6, DARS and miR128‐1 genes.

Methods

Targeted sequencing was performed in 246 severe HA patients (72 with and 174 without inhibitor): 181 previously and 65 newly enrolled.

Results

Forty‐one common and 152 rare variants passed the quality control. Logistic regression analysis of common variants identified rs3754689 and four additional variants (.011 < P < .047; FDR ranging .2‐.38). Logistic regression analysis performed only in the 220 Italian patients showed similar results (.004 < P < .05; FDR ranging .12‐.22). Three of these variants (rs3213892 and rs3816155 in the LCT intron 13 and rs961360 in the R3HDM1 intron10‐exon11 junction) may affect the expression of UBXN4 and R3HDM1, respectively. Rare variants did not show association with inhibitor development. Identified variants were not replicated in the multi‐ethnic SIPPET cohort of 230 severe HA patients.

Conclusion

Due to the limited sample size that may be responsible of the high FDR values, we could not confirm with certainty the analysed association. Further evaluation of the expression levels of analysed genes will confirm or not their role in inhibitor development.

Hearing Function: Identification of New Candidate Genes Further Explaining the Complexity of This Sensory Ability

To date, the knowledge of the genetic determinants behind the modulation of hearing ability is relatively limited. To investigate this trait, we performed Genome-Wide Association Study (GWAS) meta-analysis using genotype and audiometric data (hearing thresholds at 0.25, 0.5, 1, 2, 4, and 8 kHz, and pure-tone averages of thresholds at low, medium, and high frequencies) collected in nine cohorts from Europe, South-Eastern USA, Caucasus, and Central Asia, for an overall number of ~9000 subjects. Three hundred seventy-five genes across all nine analyses were tagged by single nucleotide polymorphisms (SNPs) reaching a suggestive p-value (p < 10⁻⁵). Amongst these, 15 were successfully replicated using a gene-based approach in the independent Italian Salus in the Apulia cohort (n = 1774) at the nominal significance threshold (p < 0.05). In addition, the expression level of the replicated genes was assessed in published human and mouse inner ear datasets. Considering expression patterns in humans and mice, eleven genes were considered particularly promising candidates for the hearing function: BNIP3L, ELP5, MAP3K20, MATN2, MTMR7, MYO1E, PCNT, R3HDM1, SLC9A9, TGFB2, and YTHDC2. These findings represent a further contribution to our understanding of the genetic basis of hearing function and its related diseases.

R3HDM1 haploinsufficiency is associated with mild intellectual disability

R3HDM1 (R3H domain containing 1) is an uncharacterized RNA-binding protein that is highly expressed in the human cerebral cortex. We report the first case of a 12-year-old Japanese male with haploinsufficiency of R3HDM1. He presented with mild intellectual disability (ID) and developmental delay. He had a pericentric inversion of 46,XY,inv(2)(p16.1q21.3)dn with breakpoints in intron 19 of R3HDM1 (2q21.3) and the intergenic region (2p16.1). The R3HDM1 levels in his lymphoblastoid cells were reduced to approximately half that of the healthy controls. However, the expression of MIR128-1, in intron 18 of R3HDM1, was not affected via the pericentric inversion. Knockdown of R3HDM1 in mouse embryonic hippocampal neurons suppressed dendritic growth and branching. Notably, the Database of Genomic Variants reported the case of a healthy control with a 488-kb deletion that included both R3HDM1 and MIR128-1. miR-128 has been reported to inhibit dendritic growth and branching in mouse brain neurons, which directly opposes the novel functions of R3HDM1. These findings suggest that deleting both R3HDM1 and MIR128-1 alleviates the symptoms of the disease caused by loss-of-function mutations in R3HDM1 only. Thus, haploinsufficiency of R3HDM1 in the patient may be the cause of the mild ID due to the genetic imbalance between R3HDM1 and MIR128-1.

Kinesin-4 KIF21B limits microtubule growth to allow rapid centrosome polarization in T cells

When a T cell and an antigen-presenting cell form an immunological synapse, rapid dynein-driven translocation of the centrosome toward the contact site leads to reorganization of microtubules and associated organelles. Currently, little is known about how the regulation of microtubule dynamics contributes to this process. Here, we show that the knockout of KIF21B, a kinesin-4 linked to autoimmune disorders, causes microtubule overgrowth and perturbs centrosome translocation. KIF21B restricts microtubule length by inducing microtubule pausing typically followed by catastrophe. Catastrophe induction with vinblastine prevented microtubule overgrowth and was sufficient to rescue centrosome polarization in KIF21B-knockout cells. Biophysical simulations showed that a relatively small number of KIF21B molecules can restrict mirotubule length and promote an imbalance of dynein-mediated pulling forces that allows the centrosome to translocate past the nucleus. We conclude that proper control of microtubule length is important for allowing rapid remodeling of the cytoskeleton and efficient T cell polarization.

The immune system is composed of many types of cells that can recognize foreign molecules and pathogens so they can eliminate them. When cells in the body become infected with a pathogen, they can process the pathogen’s proteins and present them on their own surface. Specialized immune cells can then recognize infected cells and interact with them, forming an ‘immunological synapse’. These synapses play an important role in immune response: they activate the immune system and allow it to kill harmful cells.

To form an immunological synapse, an immune cell must reorganize its internal contents, including an aster-shaped scaffold made of tiny protein tubes called microtubules. The center of this scaffold moves towards the immunological synapse as it forms. This re-orientation of the microtubules towards the immunological synapse is known as 'polarization' and it happens very rapidly, but it is not yet clear how it works.

One molecule involved in the polarization process is called KIF21B, a protein that can walk along microtubules, building up at the ends and affecting their growth. Whether KIF21B makes microtubules grow more quickly, or more slowly, is a matter of debate, and the impact microtubule length has on immunological synapse formation is unknown.

Here, Hooikaas, Damstra et al. deleted the gene for KIF21B from human immune cells called T cells to find out how it affected their ability to form an immunological synapse. Without KIF21B, the T cells grew microtubules that were longer than normal, and had trouble forming immunological synapses. When the T cells were treated with a drug that stops microtubule growth, their ability to form immunological synapses was restored, suggesting a role for KIF21B. To explore this further, Hooikaas, Damstra et al. replaced the missing KIF21B gene with a gene that coded for a version of the protein that could be seen using microscopy. This revealed that, when KIF21B reaches the ends of microtubules, it stops their growth and triggers their disassembly. Computational modelling showed that cells find it hard to reorient their microtubule scaffolding when the individual tubes are too long. It only takes a small number of KIF21B molecules to shorten the microtubules enough to allow the center of the scaffold to move.

Research has linked the KIF21B gene to autoimmune conditions like multiple sclerosis. Microtubules also play an important role in cell division, a critical process driving all types of cancer. Drugs that affect microtubule growth are already available, and a deeper understanding of KIF21B and microtubule regulation in immune cells could help to improve treatments in the future.

Celiac disease susceptibility: The genome and beyond

Celiac Disease (CeD) is an immune-mediated complex disease that is triggered by the ingestion of gluten and develops in genetically susceptible individuals. It has been known for a long time that the Human Leucocyte Antigen (HLA) molecules DQ2 and DQ8 are necessary, although not sufficient, for the disease development, and therefore other susceptibility genes and (epi)genetic events must participate in CeD pathogenesis. The advances in Genomics during the last 15 years have made CeD one of the immune-related disorders with the best-characterized genetic component. In the present work, we will first review the main Genome-Wide Association Studies (GWAS) carried out in the disorder, and emphasize post-GWAS discoveries, including diverse integrative strategies, SNP prioritization approaches, and insights into the Microbiome through the host Genomics. Second, we will explore CeD-related Epigenetics and Epigenomics, mostly focusing on the emerging knowledge of the celiac methylome, and the vast but yet under-explored non-coding RNA (ncRNA) landscape. We conclude that much has been done in the field although there are still completely unvisited areas in the post-Genomics of CeD. Chromatin conformation and accessibility, and Epitranscriptomics are promising domains that need to be unveiled to complete the big picture of the celiac Genome.

The methylome of the celiac intestinal epithelium harbours genotype-independent alterations in the HLA region

The Human Leucocyte Antigen (HLA) locus and other DNA sequence variants identified in Genome-Wide Association (GWA) studies explain around 50% of the heritability of celiac disease (CD). However, the pathogenesis of CD could be driven by other layers of genomic information independent from sequence variation, such as DNA methylation, and it is possible that allele-specific methylation explains part of the SNP associations. Since the DNA methylation landscape is expected to be different among cell types, we analyzed the methylome of the epithelial and immune cell populations of duodenal biopsies in CD patients and controls separately. We found a cell type-specific methylation signature that includes genes mapping to the HLA region, namely TAP1 and HLA-B. We also performed Immunochip SNP genotyping of the same samples and interrogated the expression of some of the affected genes. Our analysis revealed that the epithelial methylome is characterized by the loss of CpG island (CGI) boundaries, often associated to altered gene expression, and by the increased variability of the methylation across the samples. The overlap between differentially methylated positions (DMPs) and CD-associated SNPs or variants contributing to methylation quantitative trait loci (mQTLs) is minimal. In contrast, there is a notable enrichment of mQTLs among the most significant CD-associated SNPs. Our results support the notion that DNA methylation alterations constitute a genotype-independent event and confirm its role in the HLA region (apart from the well-known, DQ allele-specific effect). Finally, we find that a fraction of the CD-associated variants could exert its phenotypic effect through DNA methylation.

Bioinformatics methodologies for coeliac disease and its comorbidities

Coeliac disease (CD) is a complex, multifactorial pathology caused by different factors, such as nutrition, immunological response and genetic factors. Many autoimmune diseases are comorbidities for CD, and a comprehensive and integrated analysis with bioinformatics approaches can help in evaluating the interconnections among all the selected pathologies. We first performed a detailed survey of gene expression data available in public repositories on CD and less commonly considered comorbidities. Then we developed an innovative pipeline that integrates gene expression, cell-type data and online resources (e.g. a list of comorbidities from the literature), using bioinformatics methods such as gene set enrichment analysis and semantic similarity. Our pipeline is written in R language, available at the following link: http://bioinformatica.isa.cnr.it/COELIAC_DISEASE/SCRIPTS/. We found a list of common differential expressed genes, gene ontology terms and pathways among CD and comorbidities and the closeness among the selected pathologies by means of disease ontology terms. Physicians and other researchers, such as molecular biologists, systems biologists and pharmacologists can use it to analyze pathology in detail, from differential expressed genes to ontologies, performing a comparison with the pathology comorbidities or with other diseases.

Single Nucleotide Polymorphism in the COL11A2 Gene Associated with Heat Pain Sensitivity in Knee Osteoarthritis

Pain is one of the most prominent symptoms of osteoarthritis. However, there is often discordance between the pain experienced by individuals with osteoarthritis and the degree of articular pathology. This suggests that individual differences, including genetic variability in the central processing of nociceptive stimuli, may impact the presentation of osteoarthritis. Here, we show that the single nucleotide polymorphism rs16868943 in the collagen gene COL11A2 is significantly associated with lowered heat pain tolerance on the arm in participants with knee osteoarthritis (P = 1.21 × 10⁻⁶, P = 0.0053 after Bonferroni correction, beta = −3.42). A total of 161 knee osteoarthritis participants were included and evaluated for heat, punctate and pressure pain sensitivity of the affected knee and the ipsilateral arm. Each participant was genotyped for 4392 single nucleotide polymorphisms in genes implicated in pain perception, inflammation and mood and tested for association with pain sensitivity. The minor A allele of single nucleotide polymorphism rs16868943 was significantly associated with lower arm heat pain tolerance after correction for age, gender, race, and study site. This single nucleotide polymorphism was also nominally associated with other measures of heat pain sensitivity, including lowered knee heat pain tolerance (P = 1.14 × 10⁻⁵, P = 0.05 after Bonferroni correction), lowered arm heat pain threshold (P = 0.0039, uncorrected) and lowered knee heat pain threshold (P = 0.003, uncorrected). Addition of genotypes from 91 participants without knee pain produced a significant interaction between knee osteoarthritis status and the rs16868943 single nucleotide polymorphism in heat pain tolerance (P = 1.71 × 10⁻⁵), such that rs16868943 was not associated with heat pain tolerance in participants without knee pain (P = 0.12, beta = 1.3). This is the first study to show genetic association with heat pain tolerance in individuals with osteoarthritis. The association is specific to participants who have already developed knee osteoarthritis, suggesting that the COL11A2 gene, which has previously been associated with familial osteoarthritis, may play a role in pain sensitization after the development of osteoarthritis.

Transcription Factor Binding Site Enrichment Analysis in Co-Expression Modules in Celiac Disease

The aim of this study was to construct celiac co-expression patterns at a whole genome level and to identify transcription factors (TFs) that could drive the gliadin-related changes in coordination of gene expression observed in celiac disease (CD). Differential co-expression modules were identified in the acute and chronic responses to gliadin using expression data from a previous microarray study in duodenal biopsies. Transcription factor binding site (TFBS) and Gene Ontology (GO) annotation enrichment analyses were performed in differentially co-expressed genes (DCGs) and selection of candidate regulators was performed. Expression of candidates was measured in clinical samples and the activation of the TFs was further characterized in C2BBe1 cells upon gliadin challenge. Enrichment analyses of the DCGs identified 10 TFs and five were selected for further investigation. Expression changes related to active CD were detected in four TFs, as well as in several of their in silico predicted targets. The activation of TFs was further characterized in C2BBe1 cells upon gliadin challenge, and an increase in nuclear translocation of CAMP Responsive Element Binding Protein 1 (CREB1) and IFN regulatory factor-1 (IRF1) in response to gliadin was observed. Using transcriptome-wide co-expression analyses we are able to propose novel genes involved in CD pathogenesis that respond upon gliadin stimulation, also in non-celiac models.

Nuclear Factor-kappaB in Autoimmunity: Man and Mouse

NF-κB (nuclear factor-kappa B) is a transcription complex crucial for host defense mediated by innate and adaptive immunity, where canonical NF-κB signaling, mediated by nuclear translocation of RelA, c-Rel, and p50, is important for immune cell activation, differentiation, and survival. Non-canonical signaling mediated by nuclear translocation of p52 and RelB contributes to lymphocyte maturation and survival and is also crucial for lymphoid organogenesis. We outline NF-κB signaling and regulation, then summarize important molecular contributions of NF-κB to mechanisms of self-tolerance. We relate these mechanisms to autoimmune phenotypes described in what is now a substantial catalog of immune defects conferred by mutations in NF-κB pathways in mouse models. Finally, we describe Mendelian autoimmune syndromes arising from human NF-κB mutations, and speculate on implications for understanding sporadic autoimmune disease.

Whole exome sequencing of a consanguineous family identifies the possible modifying effect of a globally rare AK5 allelic variant in celiac disease development among Saudi patients

Celiac disease (CD), a multi-factorial auto-inflammatory disease of the small intestine, is known to occur in both sporadic and familial forms. Together HLA and Non-HLA genes can explain up to 50% of CD’s heritability. In order to discover the missing heritability due to rare variants, we have exome sequenced a consanguineous Saudi family presenting CD in an autosomal recessive (AR) pattern. We have identified a rare homozygous insertion c.1683_1684insATT, in the conserved coding region of AK5 gene that showed classical AR model segregation in this family. Sequence validation of 200 chromosomes each of sporadic CD cases and controls, revealed that this extremely rare (EXac MAF 0.000008) mutation is highly penetrant among general Saudi populations (MAF is 0.62). Genotype and allelic distribution analysis have indicated that this AK5 (c.1683_1684insATT) mutation is negatively selected among patient groups and positively selected in the control group, in whom it may modify the risk against CD development [p<0.002]. Our observation gains additional support from computational analysis which predicted that Iso561 insertion shifts the existing H-bonds between 400^th and 556^th amino acid residues lying near the functional domain of adenylate kinase. This shuffling of amino acids and their H-bond interactions is likely to disturb the secondary structure orientation of the polypeptide and induces the gain-of-function in nucleoside phosphate kinase activity of AK5, which may eventually down-regulates the reactivity potential of CD4⁺ T-cells against gluten antigens. Our study underlines the need to have population-specific genome databases to avoid false leads and to identify true candidate causal genes for the familial form of celiac disease.