Pouya Khankhanian

Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis

Monozygotic or ‘identical’ twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4+ lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among ∼3.6 million single nucleotide polymorphisms (SNPs) or ∼0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of ∼19,000 genes in CD4+ T cells. Only 2 to 176 differences in the methylation of ∼2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to ∼800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

Pathway and network-based analysis of genome-wide association studies in multiple sclerosis

Genome-wide association studies (GWAS) testing several hundred thousand SNPs have been performed in multiple sclerosis (MS) and other complex diseases. Typically, the number of markers in which the evidence for association exceeds the genome-wide significance threshold is very small, and markers that do not exceed this threshold are generally neglected. Classical statistical analysis of these datasets in MS revealed genes with known immunological functions. However, many of the markers showing modest association may represent false negatives. We hypothesize that certain combinations of genes flagged by these markers can be identified if they belong to a common biological pathway. Here we conduct a pathway-oriented analysis of two GWAS in MS that takes into account all SNPs with nominal evidence of association (P < 0.05). Gene-wise P-values were superimposed on a human protein interaction network and searches were conducted to identify sub-networks containing a higher proportion of genes associated with MS than expected by chance. These sub-networks, and others generated at random as a control, were categorized for membership of biological pathways. GWAS from eight other diseases were analyzed to assess the specificity of the pathways identified. In the MS datasets, we identified sub-networks of genes from several immunological pathways including cell adhesion, communication and signaling. Remarkably, neural pathways, namely axon-guidance and synaptic potentiation, were also over-represented in MS. In addition to the immunological pathways previously identified, we report here for the first time the potential involvement of neural pathways in MS susceptibility.

Network-Based Multiple Sclerosis Pathway Analysis with GWAS Data from 15,000 Cases and 30,000 Controls

Multiple sclerosis (MS) is an inflammatory CNS disease with a substantial genetic component, originally mapped to only the human leukocyte antigen (HLA) region. In the last 5 years, a total of seven genome-wide association studies and one meta-analysis successfully identified 57 non-HLA susceptibility loci. Here, we merged nominal statistical evidence of association and physical evidence of interaction to conduct a protein-interaction-network-based pathway analysis (PINBPA) on two large genetic MS studies comprising a total of 15,317 cases and 29,529 controls. The distribution of nominally significant loci at the gene level matched the patterns of extended linkage disequilibrium in regions of interest. We found that products of genome-wide significantly associated genes are more likely to interact physically and belong to the same or related pathways. We next searched for subnetworks (modules) of genes (and their encoded proteins) enriched with nominally associated loci within each study and identified those modules in common between the two studies. We demonstrate that these modules are more likely to contain genes with bona fide susceptibility variants and, in addition, identify several high-confidence candidates (including BCL10, CD48, REL, TRAF3, and TEC). PINBPA is a powerful approach to gaining further insights into the biology of associated genes and to prioritizing candidates for subsequent genetic studies of complex traits.

HLA Diversity in the 1000 Genomes Dataset

The 1000 Genomes Project aims to provide a deep characterization of human genome sequence variation by sequencing at a level that should allow the genome-wide detection of most variants with frequencies as low as 1%. However, in the major histocompatibility complex (MHC), only the top 10 most frequent haplotypes are in the 1% frequency range whereas thousands of haplotypes are present at lower frequencies. Given the limitation of both the coverage and the read length of the sequences generated by the 1000 Genomes Project, the highly variable positions that define HLA alleles may be difficult to identify. We used classical Sanger sequencing techniques to type the HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1 genes in the available 1000 Genomes samples and combined the results with the 103,310 variants in the MHC region genotyped by the 1000 Genomes Project. Using pairwise identity-by-descent distances between individuals and principal component analysis, we established the relationship between ancestry and genetic diversity in the MHC region. As expected, both the MHC variants and the HLA phenotype can identify the major ancestry lineage, informed mainly by the most frequent HLA haplotypes. To some extent, regions of the genome with similar genetic or similar recombination rate have similar properties. An MHC-centric analysis underlines departures between the ancestral background of the MHC and the genome-wide picture. Our analysis of linkage disequilibrium (LD) decay in these samples suggests that overestimation of pairwise LD occurs due to a limited sampling of the MHC diversity. This collection of HLA-specific MHC variants, available on the dbMHC portal, is a valuable resource for future analyses of the role of MHC in population and disease studies.

Genetic risk variants in African Americans with multiple sclerosis

Objectives: To assess the association of established multiple sclerosis (MS) risk variants in 3,254 African Americans (1,162 cases and 2,092 controls). Methods: Human leukocyte antigen (HLA)-DRB1, HLA-DQB1, and HLA-A alleles were typed by molecular techniques. Single nucleotide polymorphism (SNP) genotyping was conducted for 76 MS-associated SNPs and 52 ancestry informative marker SNPs selected throughout the genome. Self-declared ancestry was refined by principal component analysis of the ancestry informative marker SNPs. An ancestry-adjusted multivariate model was applied to assess genetic associations. Results: The following major histocompatibility complex risk alleles were replicated: HLA-DRB1*15:01 (odds ratio [OR] = 2.02 [95% confidence interval: 1.54–2.63], p = 2.50e-07), HLA-DRB1*03:01 (OR = 1.58 [1.29–1.94], p = 1.11e-05), as well as HLA-DRB1*04:05 (OR = 2.35 [1.26–4.37], p = 0.007) and the African-specific risk allele of HLA-DRB1*15:03 (OR = 1.26 [1.05–1.51], p = 0.012). The protective association of HLA-A*02:01 was confirmed (OR = 0.72 [0.55–0.93], p = 0.013). None of the HLA-DQB1 alleles were associated with MS. Using a significance threshold of p < 0.01, outside the major histocompatibility complex region, 8 MS SNPs were also found to be associated with MS in African Americans. Conclusion: MS genetic risk in African Americans only partially overlaps with that of Europeans and could explain the difference of MS prevalence between populations.

A genome-wide association study of brain lesion distribution in multiple sclerosis

Brain magnetic resonance imaging is widely used as a diagnostic and monitoring tool in multiple sclerosis and provides a non-invasive, sensitive and reproducible way to track the disease. Topological characteristics relating to the distribution and shape of lesions are recognized as important neuroradiological markers in the diagnosis of multiple sclerosis, although these have been much less well characterized quantitatively than have traditional measures such as T2 hyperintense or T1 hypointense lesion volumes. Here, we used voxel-level 3 T magnetic resonance imaging T1-weighted scans to reconstruct the 3D topology of lesions in 284 subjects with multiple sclerosis and tested whether this is a heritable phenotype. To this end, we extracted the genotypes from a published genome-wide association study on these same individuals and searched for genetic associations with lesion load, shape and topological distribution. Lesion probability maps were created to identify frequently affected areas and to assess the overall distribution of T1 lesions in the subject population as a whole. We then developed an original algorithm to cluster adjacent lesional voxels (cluxels) in each subject and tested whether cluxel topology was significantly associated with any single-nucleotide polymorphism in our data set. To focus on patterns of lesion distribution, we computed the first 10 principal components. Although principal component 1 correlated with lesion load, none of the remaining orthogonal components correlated with any other known variable. We then conducted genome-wide association studies on each of these and found 31 significant associations (false discovery rate <0.01) with principal component 8, which represents a mode of variation of lesion topology in the population. The majority of the loci can be linked to genes related to immune cell function and to myelin and neural growth; some (SYK, MYT1L, TRAPPC9, SLITKR6 and RIC3) have been previously associated with the distribution of white matter lesions in multiple sclerosis. Finally, we used a bioinformatics approach to identify a network of 48 interacting proteins showing genetic associations (P < 0.01) with cluxel topology in multiple sclerosis. This network also contains proteins expressed in immune cells and is enriched in molecules expressed in the central nervous system that contribute to neural development and regeneration. Our results show how quantitative traits derived from brain magnetic resonance images of patients with multiple sclerosis can be used as dependent variables in a genome-wide association study. With the widespread availability of powerful computing and the availability of genotyped populations, integration of imaging and genetic data sets is likely to become a mainstream tool for understanding the complex biological processes of multiple sclerosis and other brain disorders.

iCTNet: a Cytoscape plugin to produce and analyze integrative complex traits networks.

BackgroundThe speed at which biological datasets are being accumulated stands in contrast to our ability to integrate them meaningfully. Large-scale biological databases containing datasets of genes, proteins, cells, organs, and diseases are being created but they are not connected. Integration of these vast but heterogeneous sources of information will allow the systematic and comprehensive analysis of molecular and clinical datasets, spanning hundreds of dimensions and thousands of individuals. This integration is essential to capitalize on the value of current and future molecular- and cellular-level data on humans to gain novel insights about health and disease.ResultsWe describe a new open-source Cytoscape plugin named iCTNet (i ntegrated C omplex T raits Net works). iCTNet integrates several data sources to allow automated and systematic creation of networks with up to five layers of omics information: phenotype-SNP association, protein-protein interaction, disease-tissue, tissue-gene, and drug-gene relationships. It facilitates the generation of general or specific network views with diverse options for more than 200 diseases. Built-in tools are provided to prioritize candidate genes and create modules of specific phenotypes.ConclusionsiCTNet provides a user-friendly interface to search, integrate, visualize, and analyze genome-scale biological networks for human complex traits. We argue this tool is a key instrument that facilitates systematic integration of disparate large-scale data through network visualization, ultimately allowing the identification of disease similarities and the design of novel therapeutic approaches.The online database and Cytoscape plugin are freely available for academic use at: http://www.cs.queensu.ca/ictnet

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.

In depth comparison of an individual’s DNA and its lymphoblastoid cell line using whole genome sequencing

BackgroundA detailed analysis of whole genomes can be now achieved with next generation sequencing. Epstein Barr Virus (EBV) transformation is a widely used strategy in clinical research to obtain an unlimited source of a subject’s DNA. Although the mechanism of transformation and immortalization by EBV is relatively well known at the transcriptional and proteomic level, the genetic consequences of EBV transformation are less well understood. A detailed analysis of the genetic alterations introduced by EBV transformation is highly relevant, as it will inform on the usefulness and limitations of this approach.ResultsWe used whole genome sequencing to assess the genomic signature of a low-passage lymphoblastoid cell line (LCL). Specifically, we sequenced the full genome (40X) of an individual using DNA purified from fresh whole blood as well as DNA from his LCL. A total of 217.33 Gb of sequence were generated from the cell line and 238.95 Gb from the normal genomic DNA. We determined with high confidence that 99.2% of the genomes were identical, with no reproducible changes in structural variation (chromosomal rearrangements and copy number variations) or insertion/deletion polymorphisms (indels).ConclusionsOur results suggest that, at this level of resolution, the LCL is genetically indistinguishable from its genomic counterpart and therefore their use in clinical research is not likely to introduce a significant bias.

Genetic associations with brain cortical thickness in multiple sclerosis

Multiple sclerosis (MS) is characterized by temporal and spatial dissemination of demyelinating lesions in the central nervous system. Associated neurodegenerative changes contributing to disability have been recognized even at early disease stages. Recent studies show the importance of gray matter damage for the accrual of clinical disability rather than white matter where demyelination is easily visualized by magnetic resonance imaging (MRI). The susceptibility to MS is influenced by genetic risk, but genetic factors associated with the disability are not known. We used MRI data to determine cortical thickness in 557 MS cases and 75 controls and in another cohort of 219 cases. We identified nine areas showing different thickness between cases and controls (regions of interest, ROI) (eight of them were negatively correlated with Kurtzkes expanded disability status scale, EDSS) and conducted genome‐wide association studies (GWAS) in 464 and 211 cases available from the two data sets. No marker exceeded genome‐wide significance in the discovery cohort. We next combined nominal statistical evidence of association with physical evidence of interaction from a curated human protein interaction network, and searched for subnetworks enriched with nominally associated genes and for commonalities between the two data sets. This network‐based pathway analysis of GWAS detected gene sets involved in glutamate signaling, neural development and an adjustment of intracellular calcium concentration. We report here for the first time gene sets associated with cortical thinning of MS. These genes are potentially correlated with disability of MS.