Direct chromosome-length haplotyping by single-cell sequencing.

Haplotypes are fundamental to fully characterize the diploid genome of an individual, yet methods to directly chart the unique genetic makeup of each parental chromosome are lacking. Here we introduce single-cell DNA template strand sequencing (Strand-seq) as a novel approach to phasing diploid genomes along the entire length of all chromosomes. We demonstrate this by building a complete haplotype for a HapMap individual (NA12878) at high accuracy (concordance 99.3%), without using generational information or statistical inference. By use of this approach, we mapped all meiotic recombination events in a family trio with high resolution (median range ∼14 kb) and phased larger structural variants like deletions, indels, and balanced rearrangements like inversions. Lastly, the single-cell resolution of Strand-seq allowed us to observe loss of heterozygosity regions in a small number of cells, a significant advantage for studies of heterogeneous cell populations, such as cancer cells. We conclude that Strand-seq is a unique and powerful approach to completely phase individual genomes and map inheritance patterns in families, while preserving haplotype differences between single cells.

Mutations in DARS cause hypomyelination with brain stem and spinal cord involvement and leg spasticity.

Inherited white-matter disorders are a broad class of diseases for which treatment and classification are both challenging. Indeed, nearly half of the children presenting with a leukoencephalopathy remain without a specific diagnosis. Here, we report on the application of high-throughput genome and exome sequencing to a cohort of ten individuals with a leukoencephalopathy of unknown etiology and clinically characterized by hypomyelination with brain stem and spinal cord involvement and leg spasticity (HBSL), as well as the identification of compound-heterozygous and homozygous mutations in cytoplasmic aspartyl-tRNA synthetase (DARS). These mutations cause nonsynonymous changes to seven highly conserved amino acids, five of which are unchanged between yeast and man, in the DARS C-terminal lobe adjacent to, or within, the active-site pocket. Intriguingly, HBSL bears a striking resemblance to leukoencephalopathy with brain stem and spinal cord involvement and elevated lactate (LBSL), which is caused by mutations in the mitochondria-specific DARS2, suggesting that these two diseases might share a common underlying molecular pathology. These findings add to the growing body of evidence that mutations in tRNA synthetases can cause a broad range of neurologic disorders.

Deficiency in SLC25A1, encoding the mitochondrial citrate carrier, causes combined D-2- and L-2-hydroxyglutaric aciduria.

The Krebs cycle is of fundamental importance for the generation of the energetic and molecular needs of both prokaryotic and eukaryotic cells. Both enantiomers of metabolite 2-hydroxyglutarate are directly linked to this pivotal biochemical pathway and are found elevated not only in several cancers, but also in different variants of the neurometabolic disease 2-hydroxyglutaric aciduria. Recently we showed that cancer-associated IDH2 germline mutations cause one variant of 2-hydroxyglutaric aciduria. Complementary to these findings, we now report recessive mutations in SLC25A1, the mitochondrial citrate carrier, in 12 out of 12 individuals with combined D-2- and L-2-hydroxyglutaric aciduria. Impaired mitochondrial citrate efflux, demonstrated by stable isotope labeling experiments and the absence of SLC25A1 in fibroblasts harboring certain mutations, suggest that SLC25A1 deficiency is pathogenic. Our results identify defects in SLC25A1 as a cause of combined D-2- and L-2-hydroxyglutaric aciduria.

Myosin IXB variant increases the risk of celiac disease and points toward a primary intestinal barrier defect.

Celiac disease is probably the best-understood immune-related disorder. The disease presents in the small intestine and results from the interplay between multiple genes and gluten, the triggering environmental factor. Although HLA class II genes explain 40% of the heritable risk, non-HLA genes accounting for most of the familial clustering have not yet been identified. Here we report significant and replicable association (P = 2.1 x 10(-6)) to a common variant located in intron 28 of the gene myosin IXB (MYO9B), which encodes an unconventional myosin molecule that has a role in actin remodeling of epithelial enterocytes. Individuals homozygous with respect to the at-risk allele have a 2.3-times higher risk of celiac disease (P = 1.55 x 10(-5)). This result is suggestive of a primary impairment of the intestinal barrier in the etiology of celiac disease, which may explain why immunogenic gluten peptides are able to pass through the epithelial barrier.

Genetic analysis of innate immunity in Crohn's disease and ulcerative colitis identifies two susceptibility loci harboring CARD9 and IL18RAP.

The two main phenotypes of inflammatory bowel disease (IBD)--Crohn's disease (CD) and ulcerative colitis (UC)--are chronic intestinal inflammatory disorders with a complex genetic background. Using a three-stage design, we performed a functional candidate-gene analysis of innate immune pathway in IBD. In phase I, we typed 354 SNPs from 85 innate immunity genes in 520 Dutch IBD patients (284 CD, 236 UC) and 808 controls. In phase II, ten autosomal SNPs showing association at p < 0.006 in phase I were replicated in a second cohort of 545 IBD patients (326 CD, 219 UC) and 360 controls. In phase III, four SNPs with p < 0.01 in the combined phase I and phase II analysis were genotyped in an additional 786 IBD samples (452 CD, 334 UC) and 768 independent controls. Joint analysis of 1851 IBD patients (1062 CD, 789 UC) and 1936 controls demonstrated strong association to the IL18RAP rs917997 SNP for both CD and UC (p(IBD) 1.9 x 10(-8); OR 1.35). Association in CD is independently supported by the Crohn's disease dataset of the Wellcome Trust Case Control Consortium (imputed SNP rs917997, p = 9.19 x 10(-4)). In addition, an association of the CARD9 rs10870077 SNP to CD and UC was observed (p(IBD) = 3.25 x 10(-5); OR 1.21). Both genes are located in extended haplotype blocks on 2q11-2q12 and 9q34.3, respectively. Our results indicate two IBD loci and further support the importance of the innate immune system in the predisposition to both CD and UC.

Novel association in chromosome 4q27 region with rheumatoid arthritis and confirmation of type 1 diabetes point to a general risk locus for autoimmune diseases.

Recently, association of celiac disease with common single-nucleotide polymorphism (SNP) variants in an extensive linkage-disequilibrium block of 480 kb containing the KIAA1109, Tenr, IL2, and IL21 genes has been demonstrated in three independent populations (rs6822844P combined=1.3 x 10(-14)). The KIAA1109/Tenr/IL2/IL21 block corresponds to the Idd3 locus in the nonobese diabetic mouse model of type 1 diabetes (T1D). This block was recently found to be associated with T1D in a genomewide association study, although this finding lacks unequivocal confirmation. We therefore aimed to investigate whether the KIAA1109/Tenr/IL2/IL21 region is involved in susceptibility to multiple autoimmune diseases. We tested SNP rs6822844 for association with disease in 350 T1D-affected and 1,047 rheumatoid arthritis (RA)-affected Dutch patients and in 929 controls. We replicated the association with T1D (P=.0006; OR 0.64 [95% CI 0.50-0.83]), and revealed a similar novel association with RA (P=.0002; OR 0.72 [95% CI 0.61-0.86]). Our results replicate and extend the association found in the KIAA1109/Tenr/IL2/IL21 gene region with autoimmune diseases, implying that this locus is a general risk factor for multiple autoimmune diseases.

The TRAF1-C5 region on chromosome 9q33 is associated with multiple autoimmune diseases.

OBJECTIVES: The TRAF1-C5 locus has recently been identified as a genetic risk factor for rheumatoid arthritis (RA). Since genetic risk factors tend to overlap with several autoimmune diseases, a study was undertaken to investigate whether this region is associated with type 1 diabetes (TID), celiac disease (CD), systemic sclerosis (SSc) and systemic lupus erythematosus (SLE). METHODS: The most consistently associated SNP, rs10818488, was genotyped in a total of 735 patients with T1D, 1049 with CD, 367 with SSc, 746 with SLE and 3494 ethnically- and geographically-matched healthy individuals. The replication sample set consisted of 99 patients with T1D, 272 with SLE and 482 healthy individuals from Crete. RESULTS: A significant association was detected between the rs10818488 A allele and T1D (OR 1.14, p=0.027) and SLE (OR 1.16, p=0.016), which was replicated in 99 patients with T1D, 272 with SLE and 482 controls from Crete (OR 1.64, p=0.002; OR 1.43, p=0.002, respectively). Joint analysis of all patients with T1D (N=961) and all patients with SLE (N=1018) compared with 3976 healthy individuals yielded an allelic common OR of 1.19 (p=0.002) and 1.22 (p=2.6 x 10(-4)), respectively. However, combining our dataset with the T1D sample set from the WTCCC resulted in a non-significant association (OR 1.06, p=0.087). In contrast, previously unpublished results from the SLEGEN study showed a significant association of the same allele (OR 1.19, p=0.0038) with an overall effect of 1.22 (p=1.02 x 10(-6)) in a total of 1577 patients with SLE and 4215 healthy individuals. CONCLUSION: A significant association was found for the TRAF1-C5 locus in SLE, implying that this region lies in a pathway relevant to multiple autoimmune diseases.

Using mouse models to dissect the genetics of obesity.

Mice have proved to be powerful models for understanding obesity in humans and farm animals. Single-gene mutants and genetically modified mice have been used successfully to discover genes and pathways that can regulate body weight. For polygenic obesity, the most common pattern of inheritance, many quantitative trait loci (QTLs) have been mapped in crosses between selected and inbred mouse lines. Most QTL effects are additive, and diet, age and gender modify the genetic effects. Congenic, recombinant inbred, advanced intercross, and chromosome substitution strains are needed to map QTLs finely, to identify the genes underlying the traits, and to examine interactions between them.

Chromosome-wise dissection of the genome of the extremely big mouse line DU6i.

The extreme high-body-weight-selected mouse line DU6i is a polygenic model for growth research, harboring many small-effect QTL. We dissected the genome of this line into 19 autosomes and the Y chromosome by the construction of a new panel of chromosome substitution strains (CSS). The DU6i chromosomes were transferred to a DBA/2 mice genetic background by marker-assisted recurrent backcrossing. Mitochondria and the X chromosome were of DBA/2 origin in the backcross. During the construction of these novel strains, >4000 animals were generated, phenotyped, and genotyped. Using these data, we studied the genetic control of variation in body weight and weight gain at 21, 42, and 63 days. The unique data set facilitated the analysis of chromosomal interaction with sex and parent-of-origin effects. All analyzed chromosomes affected body weight and weight gain either directly or in interaction with sex or parent of origin. The effects were age specific, with some chromosomes showing opposite effects at different stages of development.

Resequencing three candidate genes for major depressive disorder in a Dutch cohort.

Major depressive disorder (MDD) is a psychiatric disorder, characterized by periods of low mood of more than two weeks, loss of interest in normally enjoyable activities and behavioral changes. MDD is a complex disorder and does not have a single genetic cause. In 2009 a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. Many of the top signals of this GWAS mapped to a region spanning the gene PCLO, and the non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became genome wide significant after post-hoc analysis. We performed resequencing of PCLO, GRM7, and SLC6A4 in 50 control samples from the GAIN-MDD cohort, to detect new genomic variants. Subsequently, we genotyped these variants in the entire GAIN-MDD cohort and performed association analysis to investigate if rs2522833 is the causal variant or simply in linkage disequilibrium with a more associated variant. GRM7 and SLC6A4 are both candidate genes for MDD from literature. We aimed to gather more evidence that rs2522833 is indeed the causal variant in the GAIN-MDD cohort or to find a previously undetected common variant in either PCLO, GRM7, or SLC6A4 with a higher association in this cohort. After next generation sequencing and association analysis we excluded the possibility of an undetected common variant to be more associated. For neither PCLO nor GRM7 we found a more associated variant. For SLC6A4, we found a new SNP that showed a lower P-value (P = 0.07) than in the GAIN-MDD GWAS (P = 0.09). However, no evidence for genome-wide significance was found. Although we did not take into account rare variants, we conclude that our results provide further support for the hypothesis that the non-synonymous coding SNP rs2522833 in the PCLO gene is indeed likely to be the causal variant in the GAIN-MDD cohort.

A fine-mapping study of 7 top scoring genes from a GWAS for major depressive disorder.

Major depressive disorder (MDD) is a psychiatric disorder that is characterized--amongst others--by persistent depressed mood, loss of interest and pleasure and psychomotor retardation. Environmental circumstances have proven to influence the aetiology of the disease, but MDD also has an estimated 40% heritability, probably with a polygenic background. In 2009, a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. A non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became only nominally significant after post-hoc analysis with an Australian cohort which used similar ascertainment. The absence of genome-wide significance may be caused by low SNP coverage of genes. To increase SNP coverage to 100% for common variants (m.a.f.>0.1, r(2)>0.8), we selected seven genes from the GAIN-MDD GWAS: PCLO, GZMK, ANPEP, AFAP1L1, ST3GAL6, FGF14 and PTK2B. We genotyped 349 SNPs and obtained the lowest P-value for rs2715147 in PCLO at P = 6.8E-7. We imputed, filling in missing genotypes, after which rs2715147 and rs2715148 showed the lowest P-value at P = 1.2E-6. When we created a haplotype of these SNPs together with the non-synonymous coding SNP rs2522833, the P-value decreased to P = 9.9E-7 but was not genome wide significant. Although our study did not identify a more strongly associated variant, the results for PCLO suggest that the causal variant is in high LD with rs2715147, rs2715148 and rs2522833.

Multiple independent variants in 6q21-22 associated with susceptibility to celiac disease in the Dutch, Finnish and Hungarian populations.

Celiac disease is an inflammatory enteropathy caused by intolerance to gluten. Previous linkage studies in the Dutch, Finnish and Hungarian populations have revealed a locus on chromosome 6q21-22 conferring susceptibility to celiac disease. This locus has previously been implicated in susceptibility to other autoimmune diseases such as Crohn's disease and type 1 diabetes. We performed fine mapping on 446 independent individuals with celiac disease and 641 controls of Dutch origin, testing 872 tagging SNPs in a 22 Mb region of chromosome 6. The 12 most promising SNPs were followed up in 2071 individuals from 284 Finnish and 357 Hungarian celiac disease families to identify risk variants in this region. Multiple markers in the region were significantly associated with celiac disease in the Dutch material. Two SNPs, rs9391227 and rs4946111, were significantly associated with celiac disease in the Finnish population. The association to rs9391227 represents the strongest association signal found in the Finnish (P = 0.003, OR 0.66) as well as the combined Dutch, Finnish and Hungarian populations (P = 3.6 × 10(-5), OR 0.76). The rs9391227 is situated downstream of the HECT domain and ankyrin repeat containing, E3 ubiquitin protein ligase 1 (HACE1) gene and is contained within a region of strong linkage disequilibrium enclosing HACE1. Two additional, independent, susceptibility variants in the 6q21-22 region were also found in a meta-analysis of the three populations. The 6q21-22 region was confirmed as a celiac disease susceptibility locus and harbors multiple independent associations, some of which may implicate ubiquitin-pathways in celiac disease susceptibility.