A comprehensive analysis of genetic risk for metabolic syndrome in the Egyptian population via allele frequency investigation and Missense3D predictions.

Diabetes mellitus (DM) represents a major health problem in Egypt and worldwide, with increasing numbers of patients with prediabetes every year. Numerous factors, such as obesity, hyperlipidemia, and hypertension, which have recently become serious concerns, affect the complex pathophysiology of diabetes. These metabolic syndrome diseases are highly linked to genetic variability that drives certain populations, such as Egypt, to be more susceptible to developing DM. Here we conduct a comprehensive analysis to pinpoint the similarities and uniqueness among the Egyptian genome reference and the 1000-genome subpopulations (Europeans, Ad-Mixed Americans, South Asians, East Asians, and Africans), aiming at defining the potential genetic risk of metabolic syndromes. Selected approaches incorporated the analysis of the allele frequency of the different populations' variations, supported by genotypes' principal component analysis. Results show that the Egyptian's reference metabolic genes were clustered together with the Europeans', Ad-Mixed Americans', and South-Asians'. Additionally, 8563 variants were uniquely identified in the Egyptian cohort, from those, two were predicted to cause structural damage, namely, CDKAL1: 6_21065070 (A > T) and PPARG: 3_12351660 (C > T) utilizing the Missense3D database. The former is a protein coding gene associated with Type 2 DM while the latter is a key regulator of adipocyte differentiation and glucose homeostasis. Both variants were detected heterozygous in two different Egyptian individuals from overall 110 sample. This analysis sheds light on the unique genetic traits of the Egyptian population that play a role in the DM high prevalence in Egypt. The proposed analysis pipeline -available through GitHub- could be used to conduct similar analysis for other diseases across populations.

A comprehensive analysis of gut and skin microbiota in canine atopic dermatitis in Shiba Inu dogs.

BACKGROUND: Like its human counterpart, canine atopic dermatitis (cAD) is a chronic relapsing condition; thus, most cAD-affected dogs will require lifelong treatment to maintain an acceptable quality of life. A potential intervention is modulation of the composition of gut microbiota, and in fact, probiotic treatment has been proposed and tried in human atopic dermatitis (AD) patients. Since dogs are currently receiving intensive medical care, this will be the same option for dogs, while evidence of gut dysbiosis in cAD is still missing, although skin microbial profiling in cAD has been conducted in several studies. Therefore, we conducted a comprehensive analysis of both gut and skin microbiota in cAD in one specific cAD-predisposed breed, Shiba Inu. Additionally, we evaluated the impact of commonly used medical management on cAD (Janus kinase; JAK inhibitor, oclacitinib) on the gut and skin microbiota. Furthermore, we genotyped the Shiba Inu dogs according to the mitochondrial DNA haplogroup and assessed its association with the composition of the gut microbiota. RESULTS: Staphylococcus was the most predominant bacterial genus observed in the skin; Escherichia/Shigella and Clostridium sensu stricto were highly abundant in the gut of cAD-affected dogs. In the gut microbiota, Fusobacteria and Megamonas were highly abundant in healthy dogs but significantly reduced in cAD-affected dogs. The abundance of these bacterial taxa was positively correlated with the effect of the treatment and state of the disease. Oclacitinib treatment on cAD-affected dogs shifted the composition of microbiota towards that in healthy dogs, and the latter brought it much closer to healthy microbiota, particularly in the gut. Additionally, even within the same dog breed, the mtDNA haplogroup varied, and there was an association between the mtDNA haplogroup and microbial composition in the gut and skin. CONCLUSIONS: Dysbiosis of both the skin and the gut was observed in cAD in Shiba Inu dogs. Our findings provide a basis for the potential treatment of cAD by manipulating the gut microbiota as well as the skin microbiota. Video Abstract.

North and East African mitochondrial genetic variation needs further characterization towards precision medicine.

INTRODUCTION: Mitochondria are maternally inherited cell organelles with their own genome, and perform various functions in eukaryotic cells such as energy production and cellular homeostasis. Due to their inheritance and manifold biological roles in health and disease, mitochondrial genetics serves a dual purpose of tracing the history as well as disease susceptibility of human populations across the globe. This work requires a comprehensive catalogue of commonly observed genetic variations in the mitochondrial DNAs for all regions throughout the world. So far, however, certain regions, such as North and East Africa have been understudied. OBJECTIVES: To address this shortcoming, we have created the most comprehensive quality-controlled North and East African mitochondrial data set to date and use it for characterizing mitochondrial genetic variation in this region. METHODS: We compiled 11 published cohorts with novel data for mitochondrial genomes from 159 Sudanese individuals. We combined these 641 mitochondrial sequences with sequences from the 1000 Genomes (n = 2504) and the Human Genome Diversity Project (n = 828) and used the tool haplocheck for extensive quality control and detection of in-sample contamination, as well as Nanopore long read sequencing for haplogroup validation of 18 samples. RESULTS: Using a subset of high-coverage mitochondrial sequences, we predict 15 potentially novel haplogroups in North and East African subjects and observe likely phylogenetic deviations from the established PhyloTree reference for haplogroups L0a1 and L2a1. CONCLUSION: Our findings demonstrate common hitherto unexplored variants in mitochondrial genomes of North and East Africa that lead to novel phylogenetic relationships between haplogroups present in these regions. These observations call for further in-depth population genetic studies in that region to enable the prospective use of mitochondrial genetic variation for precision medicine.

Genetic analysis of CFH and MCP in Egyptian patients with immune-complex proliferative glomerulonephritis.

Glomerulonephritis (GN) is a complex disease with intricate underlying pathogenic mechanisms. The possible role of underlying complement dysregulation is not fully elucidated in some GN subsets, especially in the setting of autoimmunity or infection. In the current study, diagnosed cases of lupus nephritis (LN) and post-infectious GN (PIGN) were recruited for molecular genetic analysis and targeted next-generation DNA sequencing was performed for two main complement regulating genes: in the fluid phase; CFH, and on tissue surfaces; MCP. Three heterozygous pathogenic variants in CFH (Q172*, W701*, and W1096*) and one likely pathogenic heterozygous variant in MCP (C223R) have been identified in four of the studied LN cases. Additionally, among the several detected variants of uncertain significance, one novel variant (CFH:F614S) was identified in 74% of the studied LN cases and in 65% of the studied PIGN cases. This variant was detected for the first time in the Egyptian population. These findings suggest that subtle mutations may be present in complement regulating genes in patients with immune-complex mediated category of GN that may add to the disease pathogenesis. These findings also call for further studies to delineate the impact of these gene variants on the protein function, the disease course, and outcome.

A survey of genome-wide association studies, polygenic scores and UK Biobank highlights resources for autoimmune disease genetics.

Autoimmune diseases share a general mechanism of auto-antigens harming tissues. Still. they are phenotypically diverse, with genetic as well as environmental factors contributing to their etiology at varying degrees. Associated genomic loci and variants have been identified in numerous genome-wide association studies (GWAS), whose results are increasingly used for polygenic scores (PGS) that are used to predict disease risk. At the same time, a technological shift from genotyping arrays to next generation sequencing (NGS) is ongoing. NGS allows the identification of virtually all - including rare - genetic variants, which in combination with methodological developments promises to improve the prediction of disease risk and elucidate molecular mechanisms underlying disease. Here we review current, publicly available autoimmune disease GWAS and PGS data based on information from the GWAS and PGS catalog, respectively. We summarize autoimmune diseases investigated, respective studies conducted and their results. Further, we review genetic data and autoimmune disease patients in the UK Biobank (UKB), the largest resource for genetic and phenotypic data available for academic research. We find that only comparably prevalent autoimmune diseases are covered by the UKB and at the same time assessed by both GWAS and PGS catalogs. These are systemic (systemic lupus erythematosus) as well as organ-specific, affecting the gastrointestinal tract (inflammatory bowel disease as well as specifically Crohn's disease and ulcerative colitis), joints (juvenile ideopathic arthritis, psoriatic arthritis, rheumatoid arthritis, ankylosing spondylitis), glands (Sjögren syndrome), the nervous system (multiple sclerosis), and the skin (vitiligo).

Elucidating Hexanucleotide Repeat Number and Methylation within the X-Linked Dystonia-Parkinsonism (XDP)-Related SVA Retrotransposon in TAF1 with Nanopore Sequencing.

BACKGROUND: X-linked dystonia-parkinsonism (XDP) is an adult-onset neurodegenerative disorder characterized by progressive dystonia and parkinsonism. It is caused by a SINE-VNTR-Alu (SVA) retrotransposon insertion in the TAF1 gene with a polymorphic (CCCTCT)n domain that acts as a genetic modifier of disease onset and expressivity. METHODS: Herein, we used Nanopore sequencing to investigate SVA genetic variability and methylation. We used blood-derived DNA from 96 XDP patients for amplicon-based deep Nanopore sequencing and validated it with fragment analysis which was performed using fluorescence-based PCR. To detect methylation from blood- and brain-derived DNA, we used a Cas9-targeted approach. RESULTS: High concordance was observed for hexanucleotide repeat numbers detected with Nanopore sequencing and fragment analysis. Within the SVA locus, there was no difference in genetic variability other than variations of the repeat motif between patients. We detected high CpG methylation frequency (MF) of the SVA and flanking regions (mean MF = 0.94, SD = ±0.12). Our preliminary results suggest only subtle differences between the XDP patient and the control in predicted enhancer sites directly flanking the SVA locus. CONCLUSIONS: Nanopore sequencing can reliably detect SVA hexanucleotide repeat numbers, methylation and, lastly, variation in the repeat motif.

Enhanced Glycolysis Is Required for Antileishmanial Functions of Neutrophils Upon Infection With Leishmania donovani.

Visceral leishmaniasis (VL) is a fatal parasitic disease if untreated. Treatment options of VL diminish due to emerging drug resistance. Although the principal host cells for the multiplication of Leishmania are macrophages, neutrophils are the first cells infected with the parasites rapidly after parasite inoculation. Leishmania can survive in neutrophils despite the potent antimicrobial effector functions of neutrophils that can eliminate the parasites. Recently, the growing field of immunometabolism provided strong evidence for the therapeutic potential in targeting metabolic processes as a means of controlling immune effector functions. Therefore, the understanding of the immunometabolic profile of neutrophils during Leishmania infection could provide new promising targets for host-directed therapies against VL. To our knowledge, this is the first study addressing the bioenergetics profile of L. donovani-infected primary human neutrophils. Transcriptome analysis of L. donovani-infected neutrophils revealed an early significant upregulation of several glycolytic enzymes. Extracellular flux analysis showed that glycolysis and glycolytic capacity were upregulated in L. donovani-infected neutrophils at 6 h post infection. An increased glucose uptake and accumulation of glycolytic end products were further signs for an elevated glycolytic metabolism in L. donovani-infected neutrophils. At the same time point, oxidative phosphorylation provided NADPH for the oxidative burst but did not contribute to ATP production. Inhibition of glycolysis with 2-DG significantly reduced the survival of L. donovani promastigotes in neutrophils and in culture. However, this reduction was due to a direct antileishmanial effect of 2-DG and not a consequence of enhanced antileishmanial activity of neutrophils. To further address the impact of glucose metabolism during the first days of infection in vivo, we treated C57BL/6 mice with 2-DG prior to infection with L. donovani and assessed the parasite load one day and seven days post infection. Our results show, that seven days post-infection the parasite load of 2-DG treated animals was significantly higher than in mock treated animals. This data indicates that glycolysis serves as major energy source for antimicrobial effector functions against L. donovani. Inhibition of glycolysis abrogates important neutrophil effector functions that are necessary the initial control of Leishmania infection.

Whole-genome sequencing reveals new Alzheimer's disease-associated rare variants in loci related to synaptic function and neuronal development.

INTRODUCTION: Genome-wide association studies have led to numerous genetic loci associated with Alzheimer's disease (AD). Whole-genome sequencing (WGS) now permits genome-wide analyses to identify rare variants contributing to AD risk. METHODS: We performed single-variant and spatial clustering-based testing on rare variants (minor allele frequency [MAF] ≤1%) in a family-based WGS-based association study of 2247 subjects from 605 multiplex AD families, followed by replication in 1669 unrelated individuals. RESULTS: We identified 13 new AD candidate loci that yielded consistent rare-variant signals in discovery and replication cohorts (4 from single-variant, 9 from spatial-clustering), implicating these genes: FNBP1L, SEL1L, LINC00298, PRKCH, C15ORF41, C2CD3, KIF2A, APC, LHX9, NALCN, CTNNA2, SYTL3, and CLSTN2. DISCUSSION: Downstream analyses of these novel loci highlight synaptic function, in contrast to common AD-associated variants, which implicate innate immunity and amyloid processing. These loci have not been associated previously with AD, emphasizing the ability of WGS to identify AD-associated rare variants, particularly outside of the exome.

In silico candidate variant and gene identification using inbred mouse strains.

Mice are the most widely used animal model to study genotype to phenotype relationships. Inbred mice are genetically identical, which eliminates genetic heterogeneity and makes them particularly useful for genetic studies. Many different strains have been bred over decades and a vast amount of phenotypic data has been generated. In addition, recently whole genome sequencing-based genome-wide genotype data for many widely used inbred strains has been released. Here, we present an approach for in silico fine-mapping that uses genotypic data of 37 inbred mouse strains together with phenotypic data provided by the user to propose candidate variants and genes for the phenotype under study. Public genome-wide genotype data covering more than 74 million variant sites is queried efficiently in real-time to provide those variants that are compatible with the observed phenotype differences between strains. Variants can be filtered by molecular consequences and by corresponding molecular impact. Candidate gene lists can be generated from variant lists on the fly. Fine-mapping together with annotation or filtering of results is provided in a Bioconductor package called MouseFM. In order to characterize candidate variant lists under various settings, MouseFM was applied to two expression data sets across 20 inbred mouse strains, one from neutrophils and one from CD4+ T cells. Fine-mapping was assessed for about 10,000 genes, respectively, and identified candidate variants and haplotypes for many expression quantitative trait loci (eQTLs) reported previously based on these data. For albinism, MouseFM reports only one variant allele of moderate or high molecular impact that only albino mice share: a missense variant in the Tyr gene, reported previously to be causal for this phenotype. Performing in silico fine-mapping for interfrontal bone formation in mice using four strains with and five strains without interfrontal bone results in 12 genes. Of these, three are related to skull shaping abnormality. Finally performing fine-mapping for dystrophic cardiac calcification by comparing 9 strains showing the phenotype with eight strains lacking it, we identify only one moderate impact variant in the known causal gene Abcc6. In summary, this illustrates the benefit of using MouseFM for candidate variant and gene identification.

Whole-genome sequencing reveals new Alzheimer's disease-associated rare variants in loci related to synaptic function and neuronal development.

INTRODUCTION: Genome-wide association studies have led to numerous genetic loci associated with Alzheimer's disease (AD). Whole-genome sequencing (WGS) now permit genome-wide analyses to identify rare variants contributing to AD risk. METHODS: We performed single-variant and spatial clustering-based testing on rare variants (minor allele frequency ≤1%) in a family-based WGS-based association study of 2,247 subjects from 605 multiplex AD families, followed by replication in 1,669 unrelated individuals. RESULTS: We identified 13 new AD candidate loci that yielded consistent rare-variant signals in discovery and replication cohorts (4 from single-variant, 9 from spatial-clustering), implicating these genes: FNBP1L, SEL1L, LINC00298, PRKCH, C15ORF41, C2CD3, KIF2A, APC, LHX9, NALCN, CTNNA2, SYTL3, CLSTN2. DISCUSSION: Downstream analyses of these novel loci highlight synaptic function, in contrast to common AD-associated variants, which implicate innate immunity. These loci have not been previously associated with AD, emphasizing the ability of WGS to identify AD-associated rare variants, particularly outside of coding regions.

Qtlizer: comprehensive QTL annotation of GWAS results.

Exploration of genetic variant-to-gene relationships by quantitative trait loci such as expression QTLs is a frequently used tool in genome-wide association studies. However, the wide range of public QTL databases and the lack of batch annotation features complicate a comprehensive annotation of GWAS results. In this work, we introduce the tool "Qtlizer" for annotating lists of variants in human with associated changes in gene expression and protein abundance using an integrated database of published QTLs. Features include incorporation of variants in linkage disequilibrium and reverse search by gene names. Analyzing the database for base pair distances between best significant eQTLs and their affected genes suggests that the commonly used cis-distance limit of 1,000,000 base pairs might be too restrictive, implicating a substantial amount of wrongly and yet undetected eQTLs. We also ranked genes with respect to the maximum number of tissue-specific eQTL studies in which a most significant eQTL signal was consistent. For the top 100 genes we observed the strongest enrichment with housekeeping genes (P = 2 × 10-6) and with the 10% highest expressed genes (P = 0.005) after grouping eQTLs by r2 > 0.95, underlining the relevance of LD information in eQTL analyses. Qtlizer can be accessed via https://genehopper.de/qtlizer or by using the respective Bioconductor R-package ( https://doi.org/10.18129/B9.bioc.Qtlizer ).

An integrated personal and population-based Egyptian genome reference.

A small number of de novo assembled human genomes have been reported to date, and few have been complemented with population-based genetic variation, which is particularly important for North Africa, a region underrepresented in current genome-wide references. Here, we combine long- and short-read whole-genome sequencing data with recent assembly approaches into a de novo assembly of an Egyptian genome. The assembly demonstrates well-balanced quality metrics and is complemented with variant phasing via linked reads into haploblocks, which we associate with gene expression changes in blood. To construct an Egyptian genome reference, we identify genome-wide genetic variation within a cohort of 110 Egyptian individuals. We show that differences in allele frequencies and linkage disequilibrium between Egyptians and Europeans may compromise the transferability of European ancestry-based genetic disease risk and polygenic scores, substantiating the need for multi-ethnic genome references. Thus, the Egyptian genome reference will be a valuable resource for precision medicine.

Alzheimer's disease risk SNPs show no strong effect on miRNA expression in human lymphoblastoid cell lines.

The role of microRNAs (miRNAs) in the pathogenesis of Alzheimer's disease (AD) is currently extensively investigated. In this study, we assessed the potential impact of AD genetic risk variants on miRNA expression by performing large-scale bioinformatic data integration. Our analysis was based on genetic variants from 3 AD genome-wide association studies (GWASs). Association with miRNA expression was tested by expression quantitative trait locus analysis using next-generation miRNA sequencing data generated in lymphoblastoid cell lines. Although, overall, we did not identify a strong effect of AD GWAS variants on miRNA expression in this cell type, we highlight 2 notable outliers, that is, miR-29c-5p and miR-6840-5p. MiR-29c-5p was recently reported to be involved in the regulation of BACE1 and SORL1 expression. In conclusion, despite 2 exceptions, our large-scale assessment provides only limited support for the hypothesis that AD GWAS variants act as miRNA expression quantitative trait loci.

Improved risk-stratification for posterior fossa ependymoma of childhood considering clinical, histological and genetic features - a retrospective analysis of the HIT ependymoma trial cohort.

INTRODUCTION: Risk stratification of children with ependymomas of the posterior fossa in current therapeutic protocols is mainly based on clinical criteria. We aimed to identify independent outcome predictors for this disease entity by a systematic integrated analysis of clinical, histological and genetic information in a defined cohort of patients treated according to the German HIT protocols. METHODS: Tumor samples of 134 patients aged 0.2-15.9 years treated between 1999 and 2010 according to HIT protocols were analyzed for histological features including mitotic activity, necrosis and vascular proliferation and genomic alterations by SNP and molecular inversion probe analysis. Survival analysis was performed by Kaplan-Meier method with log rank test and multivariate Cox regression analysis. RESULTS: Residual tumor after surgery, chromosome 1q gain and structural genomic alterations were identified as predictors of significantly shorter event-free (EFS) and overall survival (OS). Furthermore, specific histological features including vascular proliferation, necrosis and high mitotic activity were predictive for shorter OS. Multivariate Cox regression revealed residual tumor, chromosome 1q gain and mitotic activity as independent predictors of both EFS and OS. Using these independent predictors of outcome, we were able to build a 3-tiered risk stratification model that separates patients with standard, intermediate and high risk, and which outperforms current stratification procedures. CONCLUSION: The integration of defined clinical, histological and genetic parameters led to an improved risk-stratification model for posterior fossa ependymoma of childhood. After validation in independent cohorts this model may provide the basis for risk-adapted treatment of children with ependymomas of the posterior fossa.

Differential expression of microRNAs in Alzheimer's disease brain, blood, and cerebrospinal fluid.

INTRODUCTION: Several microRNAs (miRNAs) have been implicated in Alzheimer's disease pathogenesis, but the evidence from individual case-control studies remains inconclusive. METHODS: A systematic literature review was performed, followed by standardized multistage data extraction, quality control, and meta-analyses on eligible data for brain, blood, and cerebrospinal fluid specimens. Results were compared with miRNAs reported in the abstracts of eligible studies or recent qualitative reviews to assess novelty. RESULTS: Data from 147 independent data sets across 107 publications were quantitatively assessed in 461 meta-analyses. Twenty-five, five, and 32 miRNAs showed studywide significant differential expression (α < 1·08 × 10-4) in brain, cerebrospinal fluid, and blood-derived specimens, respectively, with 5 miRNAs showing differential expression in both brain and blood. Of these 57 miRNAs, 13 had not been reported in the abstracts of previous original or review articles. DISCUSSION: Our systematic assessment of differential miRNA expression is the first of its kind in Alzheimer's disease and highlights several miRNAs of potential relevance.

Meta-analyses identify differentially expressed micrornas in Parkinson's disease.

OBJECTIVE: MicroRNA (miRNA)-mediated (dys)regulation of gene expression has been implicated in Parkinson's disease (PD), although results of miRNA expression studies remain inconclusive. We aimed to identify miRNAs that show consistent differential expression across all published expression studies in PD. METHODS: We performed a systematic literature search on miRNA expression studies in PD and extracted data from eligible publications. After stratification for brain, blood, and cerebrospinal fluid (CSF)-derived specimen, we performed meta-analyses across miRNAs assessed in three or more independent data sets. Meta-analyses were performed using effect-size- and p-value-based methods, as applicable. RESULTS: After screening 599 publications, we identified 47 data sets eligible for meta-analysis. On these, we performed 160 meta-analyses on miRNAs quantified in brain (n = 125), blood (n = 31), or CSF (n = 4). Twenty-one meta-analyses were performed using effect sizes. We identified 13 significantly (Bonferroni-adjusted α = 3.13 × 10-4 ) differentially expressed miRNAs in brain (n = 3) and blood (n = 10) with consistent effect directions across studies. The most compelling findings were with hsa-miR-132-3p (p = 6.37 × 10-5 ), hsa-miR-497-5p (p = 1.35 × 10-4 ), and hsa-miR-133b (p = 1.90 × 10-4 ) in brain and with hsa-miR-221-3p (p = 4.49 × 10-35 ), hsa-miR-214-3p (p = 2.00 × 10-34 ), and hsa-miR-29c-3p (p = 3.00 × 10-12 ) in blood. No significant signals were found in CSF. Analyses of genome-wide association study data for target genes of brain miRNAs showed significant association (α = 9.40 × 10-5 ) of genetic variants in nine loci. INTERPRETATION: We identified several miRNAs that showed highly significant differential expression in PD. Future studies may assess the possible role of the identified brain miRNAs in pathogenesis and disease progression as well as the potential of the top blood miRNAs as biomarkers for diagnosis, progression, or prediction of PD. ANN NEUROL 2019;85:835-851.

Evidence for a potential role of miR-1908-5p and miR-3614-5p in autoimmune disease risk using integrative bioinformatics.

Genome-wide association studies (GWAS) have identified a large number of genetic risk loci for autoimmune diseases. However, the functional variants underlying these disease associations remain largely unknown. There is evidence that microRNA-mediated regulation may play an important role in this context. Therefore, we assessed whether autoimmune disease loci unfold their effects via altering microRNA expression in relevant immune cells. To this end, we performed comprehensive data integration of many large and publicly available datasets to combine information on autoimmune disease risk loci with RNA-Seq-based microRNA expression data. Specifically, we carried out microRNA expression quantitative trait loci (eQTL) analyses across 115 GWAS regions associated with 12 autoimmune diseases using next-generation sequencing data of 345 lymphoblastoid cell lines. Statistical analyses included the application and extension of a recently proposed framework (joint likelihood mapping) to microRNA expression data and microRNA target gene enrichment analyses of relevant GWAS data. Overall, only a minority of autoimmune disease risk loci may exert their pathophysiologic effects by altering microRNA expression based on JLIM. However, detailed functional fine-mapping revealed two independent GWAS regions harboring autoimmune disease risk SNPs with significant effects on microRNA expression. These relate to SNPs associated with Crohn's disease (CD; rs102275) and rheumatoid arthritis (RA; rs968567), which affect the expression of miR-1908-5p (prs102275 = 1.44e-20, prs968567 = 2.54e-14). In addition, an independent CD risk SNP, rs3853824, was found to alter the expression of miR-3614-5p (p = 5.70e-7). To support these findings, we demonstrate that GWAS signals for RA and CD were enriched in genes predicted to be targeted by both microRNAs (all with p < 0.05). In summary, our study points towards a potential pathophysiological role of miR-1908-5p and miR-3614-5p in autoimmunity.