Search

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. ‘Population-based linkage analysis’ (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10−6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority.

Multiple sclerosis (MS) is a debilitating, chronic demyelinating disease of the central nervous system affecting over 2 million people worldwide. The TAM family of receptor tyrosine kinases (TYRO3, AXL and MERTK) have been implicated as important players during demyelination in both animal models of MS and in the human disease. We therefore conducted an association study to identify single nucleotide polymorphisms (SNPs) within genes encoding the TAM receptors and their ligands associated with MS. Analysis of genotype data from a genome-wide association study which consisted of 1618 MS cases and 3413 healthy controls conducted by the Australia and New Zealand Multiple Sclerosis Genetics Consortium (ANZgene) revealed several SNPs within the MERTK gene (Chromosome 2q14.1, Accession Number NG_011607.1) that showed suggestive association with MS. We therefore interrogated 28 SNPs in MERTK in an independent replication cohort of 1140 MS cases and 1140 healthy controls. We found 12 SNPs that replicated, with 7 SNPs showing p-values of less than 10-5 when the discovery and replication cohorts were combined. All 12 replicated SNPs were in strong linkage disequilibrium with each other. In combination, these data suggest the MERTK gene is a novel risk gene for MS susceptibility.

To identify multiple sclerosis (MS) susceptibility loci, we conducted a genome-wide association study (GWAS) in 1,618 cases and used shared data for 3,413 controls. We performed replication in an independent set of 2,256 cases and 2,310 controls, for a total of 3,874 cases and 5,723 controls. We identified risk-associated SNPs on chromosome 12q13-14 (rs703842, P = 5.4 x 10(-11); rs10876994, P = 2.7 x 10(-10); rs12368653, P = 1.0 x 10(-7)) and upstream of CD40 on chromosome 20q13 (rs6074022, P = 1.3 x 10(-7); rs1569723, P = 2.9 x 10(-7)). Both loci are also associated with other autoimmune diseases. We also replicated several known MS associations (HLA-DR15, P = 7.0 x 10(-184); CD58, P = 9.6 x 10(-8); EVI5-RPL5, P = 2.5 x 10(-6); IL2RA, P = 7.4 x 10(-6); CLEC16A, P = 1.1 x 10(-4); IL7R, P = 1.3 x 10(-3); TYK2, P = 3.5 x 10(-3)) and observed a statistical interaction between SNPs in EVI5-RPL5 and HLA-DR15 (P = 0.001).

To identify multiple sclerosis (MS) susceptibility loci, we conducted a genome-wide association study (GWAS) in 1,618 cases and used shared data for 3,413 controls. We performed replication in an independent set of 2,256 cases and 2,310 controls, for a total of 3,874 cases and 5,723 controls. We identified risk-associated SNPs on chromosome 12q13-14 (rs703842, P = 5.4 x 10(-11); rs10876994, P = 2.7 x 10(-10); rs12368653, P = 1.0 x 10(-7)) and upstream of CD40 on chromosome 20q13 (rs6074022, P = 1.3 x 10(-7); rs1569723, P = 2.9 x 10(-7)). Both loci are also associated with other autoimmune diseases. We also replicated several known MS associations (HLA-DR15, P = 7.0 x 10(-184); CD58, P = 9.6 x 10(-8); EVI5-RPL5, P = 2.5 x 10(-6); IL2RA, P = 7.4 x 10(-6); CLEC16A, P = 1.1 x 10(-4); IL7R, P = 1.3 x 10(-3); TYK2, P = 3.5 x 10(-3)) and observed a statistical interaction between SNPs in EVI5-RPL5 and HLA-DR15 (P = 0.001).

In many Candida species, the leucine CUG codon is decoded by a tRNA with two unusual properties: it is a ser-tRNA and, uniquely, has guanosine at position 33 (G33). Using a combination of enzymatic (V-1, RNase, RnI nuclease) and chemical (Pb2+, imidazole) probing of the native Candida albicans ser-tRNA(CAG), we demonstrate that the overall tertiary structure of this tRNA resembles that of a ser-tRNA rather than a leu-tRNA, except within the anticodon arm where there is considerable disruption of the anticodon stem. Using non-modified in vitro transcripts of the C. albicans ser-tRNA(CAG), carrying G, C, U or A at position 33, we demonstrate that it is specifically a G residue at this position that induces the atypical anticodon stem structure. Further quantitative evidence for an unusual structure in the anticodon arm of the G33-tRNA is provided by the observed change in kinetics of methylation of the G at position 37, by purified Escherichia coli m(1)G37 methyltransferase. We conclude that the anticodon arm distortion, induced by a guanosine base at position 33 in the anticodon loop of this novel tRNA, results in reduced decoding ability which has facilitated the evolution of this tRNA without extinction of the species encoding it.

The human pathogenic yeast Candida albicans and a number of other Candida species translate the standard leucine CUG codon as serine, This is the latest addition to an increasing number of alterations to the standard genetic code which invalidate the theory that the code is frozen and universal. The unexpected finding that some organisms evolved alternative genetic codes raises two important questions: how have these alternative codes evolved and what evolutionary advantages could they create to allow for their selection? To address these questions in the context of serine CUG translation in C.albicans, we have searched for unique structural features in seryl-tRNA(CAG), which translates the leucine CUG codon as serine, and attempted to reconstruct the early stages of this genetic code switch in the closely related yeast species Saccharomyces cerevisiae. We show that a purine at position 33 (G(33)) in the C.albicans Ser-tRNA(CAG) anticodon loop, which replaces a conserved pyrimidine found in all other tRNAs, is a key structural element in the reassignment of the CUG codon from leucine to serine in that it decreases the decoding efficiency of the tRNA, thereby allowing cells to survive low level serine CUG translation. Expression of this tRNA in S.cerevisiae induces the stress response which allows cells to acquire thermotolerance. We argue that acquisition of thermotolerance may represent a positive selection for this genetic code change by allowing yeasts to adapt to sudden changes in environmental conditions and therefore colonize new ecological niches.