Early Life and Childhood Environmental Exposures, More Than Genetic Predisposition, Influence Age of Diagnosis in a Diverse Cohort of 2952 Patients With IBD.

BACKGROUND AND AIMS: Inflammatory bowel disease (IBD) develops from a combination of genetic and environmental factors. The aim of this study was to determine the contribution of established environmental risk factors and genetic risk on age of IBD diagnosis in a diverse cohort. METHODS: IBD patients in clinic completed detailed questionnaires. Blood was drawn for genetic analysis. Environmental risk factors and age of diagnosis were analyzed by ethnicity (Hispanic/Latinx or non-Hispanic White [NHW] individuals) and IBD subtype (ulcerative colitis or Crohn's disease [CD]). Weighted genetic risk scores and environmental risk scores were developed. We examined the relationship between environmental risk scores, genetic risk scores, and age of diagnosis. RESULTS: A total of 2952 patients were included: 58.9% had CD. A total of 46.83% were of Hispanic background. Early life exposures like cesarean delivery and being born in a developed country were associated with a younger age of IBD diagnosis. Childhood exposures such as frequent plastic water bottle use and having more than 1 bathroom at home were associated with a younger age of IBD. Hispanic and NHW individuals shared similar susceptibilities to environmental exposures. Environmental factors explained 21% of the variance in age of CD diagnosis and 39% in ulcerative colitis. In models incorporating genetic risk score and environmental risk score, the environment was the only significant factor associated with younger age of IBD diagnosis in all groups. CONCLUSIONS: Early life and childhood exposures impact IBD diagnosis and influence Hispanic and NHW individuals similarly. A cumulative environmental risk score contributes more to age of IBD diagnosis than genetic risk.

Not all roads lead to the immune system: the genetic basis of multiple sclerosis severity.

Multiple sclerosis is a leading cause of neurological disability in adults. Heterogeneity in multiple sclerosis clinical presentation has posed a major challenge for identifying genetic variants associated with disease outcomes. To overcome this challenge, we used prospectively ascertained clinical outcomes data from the largest international multiple sclerosis registry, MSBase. We assembled a cohort of deeply phenotyped individuals of European ancestry with relapse-onset multiple sclerosis. We used unbiased genome-wide association study and machine learning approaches to assess the genetic contribution to longitudinally defined multiple sclerosis severity phenotypes in 1813 individuals. Our primary analyses did not identify any genetic variants of moderate to large effect sizes that met genome-wide significance thresholds. The strongest signal was associated with rs7289446 (ß = -0.4882, P = 2.73 × 10-7), intronic to SEZ6L on chromosome 22. However, we demonstrate that clinical outcomes in relapse-onset multiple sclerosis are associated with multiple genetic loci of small effect sizes. Using a machine learning approach incorporating over 62 000 variants together with clinical and demographic variables available at multiple sclerosis disease onset, we could predict severity with an area under the receiver operator curve of 0.84 (95% CI 0.79-0.88). Our machine learning algorithm achieved positive predictive value for outcome assignation of 80% and negative predictive value of 88%. This outperformed our machine learning algorithm that contained clinical and demographic variables alone (area under the receiver operator curve 0.54, 95% CI 0.48-0.60). Secondary, sex-stratified analyses identified two genetic loci that met genome-wide significance thresholds. One in females (rs10967273; ßfemale = 0.8289, P = 3.52 × 10-8), the other in males (rs698805; ßmale = -1.5395, P = 4.35 × 10-8), providing some evidence for sex dimorphism in multiple sclerosis severity. Tissue enrichment and pathway analyses identified an overrepresentation of genes expressed in CNS compartments generally, and specifically in the cerebellum (P = 0.023). These involved mitochondrial function, synaptic plasticity, oligodendroglial biology, cellular senescence, calcium and G-protein receptor signalling pathways. We further identified six variants with strong evidence for regulating clinical outcomes, the strongest signal again intronic to SEZ6L (adjusted hazard ratio 0.72, P = 4.85 × 10-4). Here we report a milestone in our progress towards understanding the clinical heterogeneity of multiple sclerosis outcomes, implicating functionally distinct mechanisms to multiple sclerosis risk. Importantly, we demonstrate that machine learning using common single nucleotide variant clusters, together with clinical variables readily available at diagnosis can improve prognostic capabilities at diagnosis, and with further validation has the potential to translate to meaningful clinical practice change.

Homozygosity Haplotype and Whole-Exome Sequencing Analysis to Identify Potentially Functional Rare Variants Involved in Multiple Sclerosis among Sardinian Families.

Multiple Sclerosis (MS) is a complex multifactorial autoimmune disease, whose sex- and age-adjusted prevalence in Sardinia (Italy) is among the highest worldwide. To date, 233 loci were associated with MS and almost 20% of risk heritability is attributable to common genetic variants, but many low-frequency and rare variants remain to be discovered. Here, we aimed to contribute to the understanding of the genetic basis of MS by investigating potentially functional rare variants. To this end, we analyzed thirteen multiplex Sardinian families with Immunochip genotyping data. For five families, Whole Exome Sequencing (WES) data were also available. Firstly, we performed a non-parametric Homozygosity Haplotype analysis for identifying the Region from Common Ancestor (RCA). Then, on these potential disease-linked RCA, we searched for the presence of rare variants shared by the affected individuals by analyzing WES data. We found: (i) a variant (43181034 T > G) in the splicing region on exon 27 of CUL9; (ii) a variant (50245517 A > C) in the splicing region on exon 16 of ATP9A; (iii) a non-synonymous variant (43223539 A > C), on exon 9 of TTBK1; (iv) a non-synonymous variant (42976917 A > C) on exon 9 of PPP2R5D; and v) a variant (109859349-109859354) in 3'UTR of MYO16.

Extreme Phenotype Approach Suggests Taste Transduction Pathway for Carotid Plaque in a Multi-Ethnic Cohort.

BACKGROUND AND PURPOSE: Carotid plaque is a heritable trait and a strong predictor of vascular events. Several loci have been identified for carotid plaque, however, studies in minority populations are lacking. Within a multi-ethnic cohort, we have identified individuals with extreme total carotid plaque area (TCPA), that is, higher or lower TCPA than expected based on traditional vascular risk factors (age, sex, smoking, diabetes mellitus, hypertension, etc). We hypothesized that these individuals are enriched with genetic variants accounting for the plaque burden that cannot be explained by traditional vascular risk factors. Herein, we sought to identify the genetic basis for TCPA using the multi-ethnic cohort. METHODS: Three hundred forty participants (170 from each extreme group) from 3 race/ethnic groups (53% Hispanic, 29% non-Hispanic Black, and 18% non-Hispanic White) were genotyped using a genome-wide single-nucleotide polymorphism (SNP) array and imputed using 1000Genome data. SNP-based analyses using logistic regression and gene-based analyses using VEGAS2 were performed within each race/ethnic group and then meta-analyzed. Genes with P<0.001 were included in an overrepresentation enrichment pathway analysis using WebGestalt. Promising findings were tested for association with ischemic stroke using the MEGASTROKE Consortium data set. RESULTS: No SNP or gene reached genome-wide significance. In the pathway analysis, GO:0050913 (sensory perception of bitter taste) gene set was significantly enriched (P=4.5×10-6, false discovery rate=0.04), which was confirmed in MEGASTROKE (P=0.01). Within the GO:0050913 gene set, 3 genes were associated with extreme TCPA in our study (P<0.001): TAS2R20, TAS2R50, and ITPR3. In TAS2R50, rs1376251 is the top SNP and has been associated with myocardial infarction by others. In ITPR3, a SNP with high regulatory potential (rs3818527, RegulomeScore=1f), and ITPR3 itself were among the top SNP-based and gene-based results and showed consistent evidence for association in all ethnic groups (P<0.05). CONCLUSIONS: Extreme TCPA analysis identified new candidate genes for carotid plaque in understudied populations.

Properties of global- and local-ancestry adjustments in genetic association tests in admixed populations.

Population substructure can lead to confounding in tests for genetic association, and failure to adjust properly can result in spurious findings. Here we address this issue of confounding by considering the impact of global ancestry (average ancestry across the genome) and local ancestry (ancestry at a specific chromosomal location) on regression parameters and relative power in ancestry-adjusted and -unadjusted models. We examine theoretical expectations under different scenarios for population substructure; applying different regression models, verifying and generalizing using simulations, and exploring the findings in real-world admixed populations. We show that admixture does not lead to confounding when the trait locus is tested directly in a single admixed population. However, if there is more complex population structure or a marker locus in linkage disequilibrium (LD) with the trait locus is tested, both global and local ancestry can be confounders. Additionally, we show the genotype parameters of adjusted and unadjusted models all provide tests for LD between the marker and trait locus, but in different contexts. The local ancestry adjusted model tests for LD in the ancestral populations, while tests using the unadjusted and the global ancestry adjusted models depend on LD in the admixed population(s), which may be enriched due to different ancestral allele frequencies. Practically, this implies that global-ancestry adjustment should be used for screening, but local-ancestry adjustment may better inform fine mapping and provide better effect estimates at trait loci.

Investigating multiple sclerosis genetic susceptibility on the founder population of east-central Sardinia via association and linkage analysis of immune-related loci.

BACKGROUND: A wealth of single-nucleotide polymorphisms (SNPs) responsible for multiple sclerosis (MS) susceptibility have been identified; however, they explain only a fraction of MS heritability. OBJECTIVES: We contributed to discovery of new MS susceptibility SNPs by studying a founder population with high MS prevalence. METHODS: We analyzed ImmunoChip data from 15 multiplex families and 94 unrelated controls from the Nuoro Province, Sardinia, Italy. We tested each SNP for both association and linkage with MS, the linkage being explored in terms of identity-by-descent (IBD) sharing excess and using gene dropping to compute a corresponding empirical p-value. By targeting regions that are both associated and in linkage with MS, we increase chances of identifying interesting genomic regions. RESULTS: We identified 486 MS-associated (p < 1 × 10-4) and 18,426 MS-linked (p < 0.05) SNPs. A total of 111 loci were both linked and associated with MS, 18 of them pointing to 14 non-major histocompatibility complex (MHC) genes, and 93 of them located in the MHC region. CONCLUSION: We discovered new suggestive signals and confirmed some previously identified ones. We believe this to represent a significant step toward an understanding of the genetic basis of MS.

European multiple sclerosis risk variants in the south Asian population.

BACKGROUND: In less than a decade, genomewide association studies have identified over 100 single-nucleotide variants that are associated with increased risk of developing multiple sclerosis. However, since these studies have focused almost exclusively on European populations, it is unclear what role these variants might play in determining risk in other ethnic groups. OBJECTIVE: To assess the effects of European multiple sclerosis-associated risk variants in the south Asian population. METHODS: Using a combination of chip-based genotyping and next-generation sequencing, we have assessed 109 European-associated variants in a total of 270 cases and 555 controls from the south Asian population. RESULTS: We found that two-thirds of the tested variants (72/109) showed over representation of the European risk allele in south Asian cases (p < 0.0003). In the rest of the Immunochip array, the most associated variant was rs7318477 which maps close to TNFSF13B, the gene for the B-cell-related protein BAFF. CONCLUSION: Our data indicate substantial overlap in genetic risk architecture between Europeans and south Asians and suggest that the aetiology of the disease may be largely independent of ethnicity.

An ImmunoChip study of multiple sclerosis risk in African Americans.

The aims of this study were: (i) to determine to what degree multiple sclerosis-associated loci discovered in European populations also influence susceptibility in African Americans; (ii) to assess the extent to which the unique linkage disequilibrium patterns in African Americans can contribute to localizing the functionally relevant regions or genes; and (iii) to search for novel African American multiple sclerosis-associated loci. Using the ImmunoChip custom array we genotyped 803 African American cases with multiple sclerosis and 1516 African American control subjects at 130 135 autosomal single nucleotide polymorphisms. We conducted association analysis with rigorous adjustments for population stratification and admixture. Of the 110 non-major histocompatibility complex multiple sclerosis-associated variants identified in Europeans, 96 passed stringent quality control in our African American data set and of these, >70% (69) showed over-representation of the same allele amongst cases, including 21 with nominally significant evidence for association (one-tailed test P < 0.05). At a further eight loci we found nominally significant association with an alternate correlated risk-tagging single nucleotide polymorphism from the same region. Outside the regions known to be associated in Europeans, we found seven potentially associated novel candidate multiple sclerosis variants (P < 10(-4)), one of which (rs2702180) also showed nominally significant evidence for association (one-tailed test P = 0.034) in an independent second cohort of 620 African American cases and 1565 control subjects. However, none of these novel associations reached genome-wide significance (combined P = 6.3 × 10(-5)). Our data demonstrate substantial overlap between African American and European multiple sclerosis variants, indicating common genetic contributions to multiple sclerosis risk.

Sequencing of candidate genes in Dominican families implicates both rare exonic and common non-exonic variants for carotid intima-media thickness at bifurcation.

Through linkage and tagSNP-based association studies in 100 Dominican Republic (DR) families, we previously identified ANLN and AOAH (7p14.3) as candidate genes for carotid intima-media thickness at bifurcation (bIMT). Introns, exons, and flanking regions of ANLN and AOAH were re-sequenced in 151 individuals from nine families with evidence for linkage at 7p14.3. For common variants [CV, minor allele frequency (MAF) ≥5 %], single variant-based analysis was performed. For rare variants (RV, MAF < 5 %), gene-based analysis aggregating all RVs within a gene was performed. CV analysis revealed the strongest signal at rs3815483 (P = 0.0003) in ANLN and rs60023210 (P = 0.00005) in AOAH. In ANLN, RV analysis found suggestive evidence for association with exonic RVs (P = 0.08), and in particular non-synonymous RVs (P = 0.04) but not with all RVs (P = 0.15). The variant alleles of all non-synonymous RVs segregated with the major allele of rs3815483 and were associated with lower bIMT while a novel synonymous RV segregated with the minor allele of rs3815483 and was associated with greater bIMT. Additional analysis in 561 DR individuals found suggestive evidence for association with all ANLN non-synonymous RVs (P = 0.08). In AOAH, no evidence for association with RVs was detected. Instead, conditional analysis revealed that multiple independent intronic CVs are associated with bIMT in addition to rs60023210. We demonstrate the utility of using family-based studies to evaluate the contribution of RVs. Our data suggest two modes of genetic architecture underlying the linkage and association at ANLN (multiple exonic RVs) and AOAH (multiple intronic CVs with uncharacterized functions).

Clinical Expression of Multiple Sclerosis in Hispanic Whites of Primarily Caribbean Ancestry.

OBJECTIVE: The clinical characteristics of multiple sclerosis (MS) are not well defined in Hispanic populations. We hypothesized that disease presentation in Hispanic white (HW) patients will be different from non-Hispanic white (NHW) patients given their ancestral background and reported lower disease prevalence. This study was undertaken to compare HW of primarily Caribbean ancestry to NHW on clinical characteristics of MS. METHODS: We assessed 312 HW and 312 NHW patients with definite MS for clinical disease characteristics obtained through consented review of medical records. In order to assess the relationship between age-related phenotypes and ethnicity, linear regression was used. Logistic regression was used to assess the relationship between ethnicity and descriptors of disease presentation and severity as well as presence of neurological symptoms. RESULTS: We observed a significantly younger age at diagnosis (p = 1.38E-02) and age at exam (p = 2.36E-05) in HW. However, age at first symptom did not differ significantly between the two groups. Furthermore, within HW, the mean age at first symptom and age at diagnosis was significantly younger in those born in the United States (p < 1.00E-03 for both). Interestingly, we noted an increase in ambulatory disability in HW patients, primarily among those with relapsing disease (p = 4.18E-03). CONCLUSIONS: We found several differences in age-related phenotypes and disease severity between HW of primarily Caribbean origin and NHW patients. To our knowledge, this is the largest study to date that examined the clinical characteristics of MS in Hispanic patients of largely Caribbean origin.

Genome-wide interaction study identifies RCBTB1 as a modifier for smoking effect on carotid intima-media thickness.

OBJECTIVE: Carotid intima-media thickness (cIMT), a marker for atherosclerosis, is affected by smoking and has substantial interindividual variation. We sought to identify the genetic moderators influencing the effect of smoking on cIMT. APPROACH AND RESULTS: With a multistage design using 722 379 single nucleotide polymorphisms (SNP), a genome-wide interaction study was performed in a discovery sample of 669 Hispanics, followed by replication in 589 subjects (264 Hispanics, 172 non-Hispanic blacks, 153 non-Hispanic whites). Assuming an additive genetic model, regression analysis was performed to test for smoking-SNP interaction on cIMT while controlling for age, sex, and the top 3 principal components of ancestry. The strongest interaction in Hispanics was found with a synonymous splicing SNP (rs3751383) in exon 9 of RCBTB1 (P=2.5e(-6) in discovery sample; P=0.01 in the Hispanic replication sample; P<8.8e(-9) in the combined Hispanic sample). Stratification analysis in the combined Hispanic sample showed that smoking had no effect on cIMT among rs3751383 G homozygote (P=0.15), a moderate effect among rs3751383 heterozygote (P=0.01), and a strong effect among rs3751383 A homozygote (P=2.1e(-7)). A consistent trend was observed in the non-Hispanic white and black data sets, leading to an interaction effect of P<2.9e(-9) in the meta-analysis of all 1258 subjects. CONCLUSIONS: Our study represents the first genome-wide smoking-SNP interaction study of cIMT and identifies RCBTB1 as a modifier of the smoking effect on cIMT. Testing for gene-environment interactions can help uncover genetic factors that contribute to the interindividual variation in response to the same environmental exposure.

Interrogating the complex role of chromosome 16p13.13 in multiple sclerosis susceptibility: independent genetic signals in the CIITA-CLEC16A-SOCS1 gene complex.

Multiple sclerosis (MS) is a neurodegenerative, autoimmune disease of the central nervous system, and numerous studies have shown that MS has a strong genetic component. Independent studies to identify MS-associated genes have often indicated multiple signals in physically close genomic regions, although by their proximity it is not always clear if these data indicate redundant or truly independent genetic signals. Recently, three MS study samples were genotyped in parallel using an Illumina Custom BeadChip. These revealed multiple significantly associated single-nucleotide polymorphisms within a 600 kb stretch on chromosome 16p13. Here we present a detailed analysis of variants in this region that clarifies the independent nature of these signals. The linkage disequilibrium patterns in the region and logistic regression analysis of the associations suggest that this region likely harbors three independent MS disease loci. Further, we examined cis-expression QTLs, histone modifications and CCCTC-binding factor (CTCF) binding data in the region. We also tested for correlated expression of the genes from the region using whole-genome expression array data from lymphoblastoid cell lines. Three of the genes show expression correlations across loci. Furthermore, in the GM12878 lymphoblastoid cell line, these three genes are in a continuous region devoid of H3K27 methylation, suggesting an open chromatin configuration. This region likely only contributes minimal risk to MS; however, investigation of this region will undoubtedly provide insight into the functional mechanisms of these genes. These data highlight the importance of taking a closer look at the expression and function of chromosome 16p13 in the pathogenesis of MS.

Genome-wide association study of executive function in a multi-ethnic cohort implicates LINC01362: Results from the northern Manhattan study.

Executive function is a cognitive domain with sizable heritability representing higher-order cognitive abilities. Genome-wide association studies (GWAS) of executive function are sparse, particularly in populations underrepresented in medical research. We performed a GWAS on a composite measure of executive function that included measures of mental flexibility and reasoning using data from the Northern Manhattan Study, a racially and ethnically diverse cohort (N = 1077, 69% Hispanic, 17% non-Hispanic Black and 14% non-Hispanic White). Four SNPs located in the long intergenic non-protein coding RNA 1362 gene, LINC01362, on chromosome 1p31.1, were significantly associated with the composite measure of executive function in this cohort (top SNP rs2788328, ß = 0.22, p = 3.1 × 10-10). The associated SNPs have been shown to influence expression of the tubulin tyrosine ligase like 7 gene, TTLL7 and the protein kinase CAMP-activated catalytic subunit beta gene, PRKACB, in several regions of the brain involved in executive function. Together, these findings present new insight into the genetic underpinnings of executive function in an understudied population.

Heritability Estimation of Multiple Sclerosis Related Plasma Protein Levels in Sardinian Families with Immunochip Genotyping Data.

This work aimed at estimating narrow-sense heritability, defined as the proportion of the phenotypic variance explained by the sum of additive genetic effects, via Haseman-Elston regression for a subset of 56 plasma protein levels related to Multiple Sclerosis (MS). These were measured in 212 related individuals (with 69 MS cases and 143 healthy controls) obtained from 20 Sardinian families with MS history. Using pedigree information, we found seven statistically significant heritable plasma protein levels (after multiple testing correction), i.e., Gc (h2 = 0.77; 95%CI: 0.36, 1.00), Plat (h2 = 0.70; 95%CI: 0.27, 0.95), Anxa1 (h2 = 0.68; 95%CI: 0.27, 1.00), Sod1 (h2 = 0.58; 95%CI: 0.18, 0.96), Irf8 (h2 = 0.56; 95%CI: 0.19, 0.99), Ptger4 (h2 = 0.45; 95%CI: 0.10, 0.96), and Fadd (h2 = 0.41; 95%CI: 0.06, 0.84). A subsequent analysis was performed on these statistically significant heritable plasma protein levels employing Immunochip genotyping data obtained in 155 healthy controls (92 related and 63 unrelated); we found a meaningful proportion of heritable plasma protein levels' variability explained by a small set of SNPs. Overall, the results obtained, for these seven MS-related proteins, emphasized a high additive genetic variance component explaining plasma levels' variability.

Plasma Protein Levels Analysis in Multiple Sclerosis Sardinian Families Identified C9 and CYP24A1 as Candidate Biomarkers.

Here we investigate protein levels in 69 multiple sclerosis (MS) cases and 143 healthy controls (HC) from twenty Sardinian families to search for promising biomarkers in plasma. Using antibody suspension bead array technology, the plasma levels of 56 MS-related proteins were obtained. Differences between MS cases and HC were estimated using Linear Mixed Models or Linear Quantile Mixed Models. The proportion of proteins level variability, explained by a set of 119 MS-risk SNPs as to the literature, was also quantified. Higher plasma C9 and CYP24A1 levels were found in MS cases compared to HC (p < 0.05 after Holm multiple testing correction), with protein level differences estimated as, respectively, 0.53 (95% CI: 0.25, 0.81) and 0.42 (95% CI: 0.19, 0.65) times plasma level standard deviation measured in HC. Furthermore, C9 resulted in both statistically significantly higher relapsing-remitting MS (RRMS) and secondary-progressive MS (SPMS) compared to HC, with SPMS showing the highest differences. Instead, CYP24A1 was statistically significantly higher only in RRMS as compared to HC. Respectively, 26% (95% CI: 10%, 44%) and 16% (95% CI: 9%, 39%) of CYP24A1 and C9 plasma level variability was explained by known MS-risk SNPs. Our results highlight C9 and CYP24A1 as potential biomarkers in plasma for MS and allow us to gain insight into molecular disease mechanisms.

Rare variants in previously identified linkage regions associated with carotid plaque in Dominican Republic families.

Carotid plaque is a subclinical measure of atherosclerosis. We have previously shown measures of carotid plaque to be heritable in a sample of 100 Dominican families and found evidence for linkage and association of common variants (CVs) on 7q36, 11p15, 14q32 and 15q23 with plaque presence. Our current study aimed to refine these regions further and identify rare variants (RVs) influencing plaque presence. Therefore, we performed targeted sequencing of the one LOD unit down region on 7q36, 11p15, 14q32 and 15q23 in 12 Dominican families with evidence for linkage to plaque presence. Gene-based RV analyses were performed using the Sequence Association Test for familial data (F-SKAT) under two filtering algorithms; 1. all exonic RVs and 2. non-synonymous RVs. Replication analyses were performed using a sample of 22 Dominican families and 556 unrelated Dominicans with Exome Array data. To identify additional non-synonymous RVs influencing plaque, we looked for co-segregation of RVs with plaque in each of the sequenced families. Our most strongly associated gene with evidence for replication was AMPD3 which showed suggestive association with plaque presence in the sequenced families (exonic RV p = 0.003, nonsynonymous RV p = 0.005) and replication families (exonic RV p = 0.04, nonsynonymous RV p = 0.02). Examination of the sequenced family pedigrees revealed two missense variants on chromosome 11 which co-segregated with plaque presence in one of our families; rs61751342 (located in DENND2B), and rs61760882 (located in RNF141). The rs61751342 missense variant is an eQTL for SCUBE2 in the atrial appendage. Notably, SCUBE2 encodes a protein which interacts with vascular endothelial growth factor (VEGF) receptor 2 to regulate VEGF-induced angiogenesis, thus providing biologic plausibility for this gene in atherosclerosis. In conclusion, using targeted sequencing of previously-identified linkage regions, we have identified suggestive evidence for the role of RVs in carotid plaque pathogenesis.

Lamina Propria Phagocyte Profiling Reveals Targetable Signaling Pathways in Refractory Inflammatory Bowel Disease.

BACKGROUND AND AIMS: Lamina propria phagocytes are key mediators of inflammatory bowel disease (IBD). We aimed to understand the transcriptomic and functional differences in these cells based on location, disease type, inflammation state, and medication use in patients with IBD. METHODS: Phagocytic immune cells in the lamina propria, as defined by the marker CD11b, were isolated from 54 unique patients (n = 111 gut mucosal biopsies). We performed flow cytometry for cell phenotyping (n = 30) and RNA sequencing with differential gene expression analysis (n = 58). We further cultured these cells in vitro and exposed them to janus kinase inhibitors to measure cytokine output (n = 27). Finally, we matched patient genomic data to our RNA sequencing data to perform candidate gene expression quantitative trait locus analysis (n = 34). RESULTS: We found distinct differences in gene expression between CD11b+ cells from the colon vs ileum, as well as in different inflammatory states and, to a lesser degree, IBD types (Crohn's disease or ulcerative colitis). These genes mapped to targetable immune pathways and metabolic and cancer pathways. We further explored the janus kinase-signal transducer and activator of transcription pathway, which was upregulated across many comparisons including in biopsies from anti-tumor necrosis factor refractory patients. We found that isolated CD11b+ cells treated with janus kinase inhibitors had decreased secretion of cytokines tumor necrosis factora and interleukin-8 (P ≤ .05). We also found 3 genetic variants acting as expression quantitative trait loci (P ≤ .1) within our CD11b+ data set. CONCLUSIONS: Lamina propria phagocytes from IBD mucosa provide pathogenetic clues on the nature of treatment refractoriness and inform new targets for therapy.

Ancestral risk modification for multiple sclerosis susceptibility detected across the Major Histocompatibility Complex in a multi-ethnic population.

The Major Histocompatibility Complex (MHC) makes the largest genetic contribution to multiple sclerosis (MS) susceptibility, with 32 independent effects across the region explaining 20% of the heritability in European populations. Variation is high across populations with allele frequency differences and population-specific risk alleles identified. We sought to identify MHC-specific MS susceptibility variants and assess the effect of ancestral risk modification within 2652 Latinx and Hispanic individuals as well as 2435 Black and African American individuals. We have identified several novel susceptibility alleles which are rare in European populations including HLA-B*53:01, and we have utilized the differing linkage disequilibrium patterns inherent to these populations to identify an independent role for HLA-DRB1*15:01 and HLA-DQB1*06:02 on MS risk. We found a decrease in Native American ancestry in MS cases vs controls across the MHC, peaking near the previously identified MICB locus with a decrease of ~5.5% in Hispanics and ~0.4% in African Americans. We have identified several susceptibility variants, including within the MICB gene region, which show global ancestry risk modification and indicate ancestral differences which may be due in part to correlated environmental factors. We have also identified several susceptibility variants for which MS risk is modified by local ancestry and indicate true ancestral genetic differences; including HLA-DQB1*06:02 for which MS risk for European allele carriers is almost two times the risk for African allele carriers. These results validate the importance of investigating MS susceptibility at an ancestral level and offer insight into the epidemiology of MS phenotypic diversity.

Large-scale sequencing identifies multiple genes and rare variants associated with Crohn's disease susceptibility.

Genome-wide association studies (GWASs) have identified hundreds of loci associated with Crohn's disease (CD). However, as with all complex diseases, robust identification of the genes dysregulated by noncoding variants typically driving GWAS discoveries has been challenging. Here, to complement GWASs and better define actionable biological targets, we analyzed sequence data from more than 30,000 patients with CD and 80,000 population controls. We directly implicate ten genes in general onset CD for the first time to our knowledge via association to coding variation, four of which lie within established CD GWAS loci. In nine instances, a single coding variant is significantly associated, and in the tenth, ATG4C, we see additionally a significantly increased burden of very rare coding variants in CD cases. In addition to reiterating the central role of innate and adaptive immune cells as well as autophagy in CD pathogenesis, these newly associated genes highlight the emerging role of mesenchymal cells in the development and maintenance of intestinal inflammation.

Genetic loci for blood lipid levels identified by linkage and association analyses in Caribbean Hispanics.

To identify genetic loci influencing blood lipid levels in Caribbean Hispanics, we first conducted a genome-wide linkage scan in 1,211 subjects from 100 Dominican families on five lipid quantitative traits: total cholesterol (TC), low density lipoprotein cholesterol (LDL-C), high density lipoprotein cholesterol (HDL-C), triglycerides (TG), and LDL-C/HDL-C ratio. We then investigated the association between blood lipid levels and 21,361 single nucleotide polymorphisms (SNP) under the 1-logarithm of odds (LOD) unit down regions of linkage peaks in an independent community-based subcohort (N = 814, 42% Dominican) from the Northern Manhattan Study (NOMAS). We found significant linkage evidence for LDL-C/HDL-C on 7p12 (multipoint LOD = 3.91) and for TC on 16q23 (LOD = 3.35). In addition, we identified suggestive linkage evidence of LOD > 2.0 on 15q23 for TG, 16q23 for LDL-C, 19q12 for TC and LDL-C, and 20p12 for LDL-C. In the association analysis of the linkage peaks, we found that seven SNPs near FLJ45974 were associated with LDL-C/HDL-C with a nominal P < 3.5 × 10(-5), in addition to associations (P < 0.0001) for other lipid traits with SNPs in or near CDH13, SUMF2, TLE3, FAH, ARNT2, TSHZ3, ZNF343, RPL7AL2, and TMC3. Further studies are warranted to perform in-depth investigations of functional genetic variants in these regions.