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.

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.

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.

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.

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.

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.

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.

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.

Two genetic variants explain the association of European ancestry with multiple sclerosis risk in African-Americans.

Epidemiological studies have suggested differences in the rate of multiple sclerosis (MS) in individuals of European ancestry compared to African ancestry, motivating genetic scans to identify variants that could contribute to such patterns. In a whole-genome scan in 899 African-American cases and 1155 African-American controls, we confirm that African-Americans who inherit segments of the genome of European ancestry at a chromosome 1 locus are at increased risk for MS [logarithm of odds (LOD) = 9.8], although the signal weakens when adding an additional 406 cases, reflecting heterogeneity in the two sets of cases [logarithm of odds (LOD) = 2.7]. The association in the 899 individuals can be fully explained by two variants previously associated with MS in European ancestry individuals. These variants tag a MS susceptibility haplotype associated with decreased CD58 gene expression (odds ratio of 1.37; frequency of 84% in Europeans and 22% in West Africans for the tagging variant) as well as another haplotype near the FCRL3 gene (odds ratio of 1.07; frequency of 49% in Europeans and 8% in West Africans). Controlling for all other genetic and environmental factors, the two variants predict a 1.44-fold higher rate of MS in European-Americans compared to African-Americans.

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.

Investigating the Causal Effect of Brain Expression of CCL2, NFKB1, MAPK14, TNFRSF1A, CXCL10 Genes on Multiple Sclerosis: A Two-Sample Mendelian Randomization Approach.

Multiple Sclerosis (MS) exhibits considerable heterogeneity in phenotypic expression, course, prognosis and response to therapy. This suggests this disease involves multiple, as yet poorly understood, causal mechanisms. In this work we assessed the possible causal link between gene expression level of five selected genes related to the pro-inflammatory NF-κB signaling pathway (i.e., CCL2, NFKB1, MAPK14, TNFRSF1A, CXCL10) in ten different brain tissues (i.e., cerebellum, frontal cortex, hippocampus, medulla, occipital cortex, putamen, substantia nigra, thalamus, temporal cortex and intralobular white matter) and MS. We adopted a two-stage Mendelian Randomization (MR) approach for the estimation of the causal effects of interest, based on summary-level data from 20 multiplex Sardinian families and data provided by the United Kingdom Brain Expression Consortium (UKBEC). Through Radial-MR and Cochrane's Q statistics we identified and removed genetic variants which are most likely to be invalid instruments. To estimate the total causal effect, univariable MR was carried out separately for each gene and brain region. We used Inverse-Variance Weighted estimator (IVW) as main analysis and MR-Egger Regression estimator (MR-ER) and Weighted Median Estimator (WME) as sensitivity analysis. As these genes belong to the same pathway and thus they can be closely related, we also estimated their direct causal effects by applying IVW and MR-ER within a multivariable MR (MVMR) approach using set of genetic instruments specific and common (composite) to each multiple exposures represented by the expression of the candidate genes. Univariate MR analysis showed a significant positive total causal effect for CCL2 and NFKB1 respectively in medulla and cerebellum. MVMR showed a direct positive causal effect for NFKB1 and TNFRSF1A, and a direct negative causal effect for CCL2 in cerebellum; while in medulla we observed a direct positive causal effect for CCL2. Since in general we observed a different magnitude for the gene specific causal effect we hypothesize that in cerebellum and medulla the effect of each gene expression is direct but also mediated by the others. These results confirm the importance of the involvement of NF-κB signaling pathway in brain tissue for the development of the disease and improve our understanding in the pathogenesis of MS.

Fine-Mapping Array Design for Multi-Ethnic Studies of Multiple Sclerosis.

While approximately 200 autosomal genetic associations outside of the major histocompatibility complex (MHC) have been identified for multiple sclerosis (MS) risk in European populations, causal variants identified at the majority of these associated loci have been much more elusive. We propose that knowledge gained from replication efforts in Hispanic and African American populations can be utilized to more efficiently fine-map these risk loci. To this end, we have customized a genotyping array by adding ~20,000 bead types (~17,000 variants) to the base content of the Ilumina Infinium expanded multi-ethnic genotyping array and the Infinium ImmunoArray-24 v2 BeadChip. These custom bead types were chosen to allow for the detection of causal variation (1) in the presence of allelic and locus heterogeneity, by incorporating regulatory and coding variation within 1-Mb of previously identified risk variants and (2) in the absence of allelic and locus heterogeneity by incorporation of variants using linkage disequilibrium criteria, which are based on knowledge of replication status in Hispanic and African American study samples. This array has been designed to maximize fine-mapping potential for currently identified MS susceptibility loci, particularly in multi-ethnic populations. The strategies described here could be additionally informative for fine-mapping of other disease phenotypes.

Targeted sequencing of linkage region in Dominican families implicates PRIMA1 and the SPATA7-PTPN21-ZC3H14-EML5-TTC8 locus in carotid-intima media thickness and atherosclerotic events.

Carotid intima-media thickness (cIMT) is a subclinical marker for atherosclerosis. Previously, we reported a quantitative trait locus (QTL) for total cIMT on chromosome 14q and identified PRiMA1, FOXN3 and CCDC88C as candidate genes using a common variants (CVs)-based approach. Herein, we further evaluated the genetic contribution of the QTL to cIMT by resequencing. We sequenced all exons within the QTL and genomic regions of PRiMA1, FOXN3 and CCDC88C in Dominican families with evidence for linkage to the QTL. Unrelated Dominicans from the Northern Manhattan Study (NOMAS) were used for validation. Single-variant-based and gene-based analyses were performed for CVs and rare variants (RVs). The strongest evidence for association with CVs was found in PRiMA1 (p = 8.2 × 10-5 in families, p = 0.01 in NOMAS at rs12587586), and in the five-gene cluster SPATA7-PTPN21-ZC3H14-EML5-TTC8 locus (p = 1.3 × 10-4 in families, p = 0.01 in NOMAS at rs2274736). No evidence for association with RVs was found in PRiMA1. The top marker from previous study in PRiMA1 (rs7152362) was associated with fewer atherosclerotic events (OR = 0.67; p = 0.02 in NOMAS) and smaller cIMT (ß = -0.58, p = 2.8 × 10-4 in Family). Within the five-gene cluster, evidence for association was found for exonic RVs (p = 0.02 in families, p = 0.28 in NOMAS), which was enriched among RVs with higher functional potentials (p = 0.05 in NOMAS for RVs in the top functional tertile). In summary, targeted resequencing provided validation and novel insights into the genetic architecture of cIMT, suggesting stronger effects for RVs with higher functional potentials. Furthermore, our data support the clinical relevance of CVs associated with subclinical atherosclerosis.

Native ancestry is associated with optic neuritis and age of onset in hispanics with multiple sclerosis.

BACKGROUND AND OBJECTIVE: Hispanics with multiple sclerosis (MS) present younger and more often with optic neuritis (ON) as compared to Whites in the western United States. Regional differences related to Hispanic genetic admixture could be responsible. We investigated the association between global genetic ancestry and ON and age at onset of MS in Hispanics. METHODS: Data were obtained for 1033 self-identified Hispanics with MS from four MS-based registries from four academic institutions across the United States January 2016-April 2017. Multivariate regression models, utilizing genetic ancestry estimates for Native American (NA), African, and European ancestry, were used to assess the relationship between genetic ancestry and ON presentation and age of MS onset, defined as age at first symptom. RESULTS: Genetic ancestry and ON proportions varied by region where NA ancestry and ON proportions were highest among Hispanics in the southwestern United States (40% vs. 19% overall for NA and 38% vs. 25% overall for ON). A strong inverse correlation was observed between NA and European ancestry (r = -0.83). ON presentation was associated with younger age of onset (OR: 0.98; 95% CI: 0.96-0.99; P = 7.80 × 10-03) and increased NA ancestry (OR: 2.35 for the highest versus the lowest quartile of NA ancestry; 95% CI: 1.35-4.10; P = 2.60 × 10-03). Younger age of onset was found to be associated with a higher proportion NA (Beta: -5.58; P = 3.49 × 10-02) and African ancestry (Beta: -10.07; P = 1.39 × 10-03). INTERPRETATION: Ethnic differences associated with genetic admixture could influence clinical presentation in Hispanics with MS; underscoring the importance of considering genetic substructure in future clinical, genetic, and epigenetic studies in Hispanics.

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.

Variants in chondroitin sulfate metabolism genes in thrombotic storm.

INTRODUCTION: Thrombotic storm (TS) presents as a severe, acute thrombotic phenotype, characterized by multiple clotting events and frequently affecting younger adults. Understanding the extensive hypercoagulation of an extreme phenotype as TS will also provide insight into the pathogenesis of a wider spectrum of thrombotic disorders. MATERIAL AND METHODS: We completed whole exome sequencing on 26 TS patients, including 1 multiplex family, 13 trios and 12 isolated TS patients. We examined both dominant and recessive inheritance models for known thrombotic factors as well as performed a genome-wide screen. Identified genes of interest in the family and trios were screened in the remaining TS patients. Variants were filtered on frequency (<5% in 1000 genomes), conservation and function in gene and were annotated for effect on protein and overall functionality. RESULTS: We observed an accumulation of variants in genes linked to chondroitin sulfate (CS), but not heparan sulfate metabolism. Sixteen conserved, rare missense and nonsense variants in genes involved in CS metabolism (CHPF, CHPF2, CHST3, CHST12, CHST15, SLC26A2, PAPSS2, STAB2) were identified in over one-third of the TS patients. In contrast, we identified only seven variants in known thrombosis genes (including FV Leiden). CONCLUSIONS: As CS has multiple functions in the glycocalyx protecting the endothelial cells, reduced availability of CS could diminish the normal control mechanisms for blood coagulation, making these CS metabolism genes strong potential risk factors for TS. Overall, no single gene was identified with strong evidence for TS causality; however, our data suggest TS is mediated by an accumulation of rare pro-thrombotic risk factors.

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.