Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted re-sequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrl änder dogs. To our knowledge, this is the first time this family trio WGS-approach, has successfully been used to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrl änder dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G>C (p.R52P) was found to be concordant in eight additional cases and 16 healthy Kromfohrl änder dogs.

H19 lncRNA alters DNA methylation genome wide by regulating S-adenosylhomocysteine hydrolase

DNA methylation is essential for mammalian development and physiology. Here we report that the developmentally regulated H19 lncRNA binds to and inhibits S-adenosylhomocysteine hydrolase (SAHH), the only mammalian enzyme capable of hydrolysing S-adenosylhomocysteine (SAH). SAH is a potent feedback inhibitor of S-adenosylmethionine (SAM)-dependent methyltransferases that methylate diverse cellular components, including DNA, RNA, proteins, lipids and neurotransmitters. We show that H19 knockdown activates SAHH, leading to increased DNMT3B-mediated methylation of an lncRNA-encoding gene Nctc1 within the Igf2-H19-Nctc1 locus. Genome-wide methylation profiling reveals methylation changes at numerous gene loci consistent with SAHH modulation by H19. Our results uncover an unanticipated regulatory circuit involving broad epigenetic alterations by a single abundantly expressed lncRNA that may underlie gene methylation dynamics of development and diseases and suggest that this mode of regulation may extend to other cellular components.

Impact of single nucleotide polymorphisms and of clinical risk factors on new-onset diabetes mellitus in HIV-infected individuals

Background. Metabolic complications, including cardiovascular events and diabetes mellitus (DM), are a major long-term concern in human immunodeficienc virus (HIV)-infected individuals. Recent genome-wide association studies have reliably associated multiple single nucleotide polymorphisms (SNPs) to DM in the general population. Methods. We evaluated the contribution of 22 SNPs identifie in genome-wide association studies and of longitudinally measured clinical factors to DM. We genotyped all 94 white participants in the Swiss HIV Cohort Study who developed DM from 1 January 1999 through 31 August 2009 and 550 participants without DM. Analyses were based on 6054 person-years of follow-up and 13,922 measurements of plasma glucose. Results. The contribution to DM risk explained by SNPs (14% of DM variability) was larger than the contribution to DM risk explained by current or cumulative exposure to different antiretroviral therapy combinations (3% of DM variability). Participants with the most unfavorable genetic score (representing 12% and 19% of the study population, respectively, when applying 2 different genetic scores) had incidence rate ratios for DM of 3.80 (95% confidenc interval [CI], 2.05–7.06) and 2.74 (95% CI, 1.53–4.88), respectively, compared with participants with a favorable genetic score. However, addition of genetic data to clinical risk factors that included body mass index only slightly improved DM prediction. Conclusions. In white HIV-infected persons treated with antiretroviral therapy, the DM effect of genetic variants was larger than the potential toxic effects of antiretroviral therapy. SNPs contributed significantl to DM risk, but their addition to a clinical model improved DM prediction only slightly, similar to studies in the general population.

Child and adolescent psychiatric genetics

The current status of child and adolescent psychiatric genetics appears promising in light of the initiation of genome-wide association studies (GWAS) for diverse polygenic disorders and the molecular elucidation of monogenic Rett syndrome, for which recent functional studies provide hope for pharmacological treatment strategies. Within the last 50 years, tremendous progress has been made in linking genetic variation to behavioral phenotypes and psychiatric disorders. We summarize the major findings of the Human Genome Project and dwell on largely unsuccessful candidate gene and linkage studies. GWAS for the first time offer the possibility to detect single nucleotide polymorphisms and copy number variants without a priori hypotheses as to their molecular etiology. At the same time it is becoming increasingly clear that very large sample sizes are required in order to enable genome wide significant findings, thus necessitating further large-scaled ascertainment schemes for the successful elucidation of the molecular genetics of childhood and adolescent psychiatric disorders. We conclude by reflecting on different scenarios for future research into the molecular basis of early onset psychiatric disorders. This review represents the introductory article of this special issue of the European Child and Adolescent Psychiatry.

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies and candidate gene studies in ulcerative colitis have identified 18 susceptibility loci. We conducted a meta-analysis of six ulcerative colitis genome-wide association study datasets, comprising 6,687 cases and 19,718 controls, and followed up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P < 5 × 10(-8)), increasing the number of ulcerative colitis-associated loci to 47. After annotating associated regions using GRAIL, expression quantitative trait loci data and correlations with non-synonymous SNPs, we identified many candidate genes that provide potentially important insights into disease pathogenesis, including IL1R2, IL8RA-IL8RB, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease risk loci is now 99, including a minimum of 28 shared association signals between Crohn's disease and ulcerative colitis.

Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology

Genome-wide association studies (GWAS) have defined over 150 genomic regions unequivocally containing variation predisposing to immune-mediated disease. Inferring disease biology from these observations, however, hinges on our ability to discover the molecular processes being perturbed by these risk variants. It has previously been observed that different genes harboring causal mutations for the same Mendelian disease often physically interact. We sought to evaluate the degree to which this is true of genes within strongly associated loci in complex disease. Using sets of loci defined in rheumatoid arthritis (RA) and Crohn's disease (CD) GWAS, we build protein-protein interaction (PPI) networks for genes within associated loci and find abundant physical interactions between protein products of associated genes. We apply multiple permutation approaches to show that these networks are more densely connected than chance expectation. To confirm biological relevance, we show that the components of the networks tend to be expressed in similar tissues relevant to the phenotypes in question, suggesting the network indicates common underlying processes perturbed by risk loci. Furthermore, we show that the RA and CD networks have predictive power by demonstrating that proteins in these networks, not encoded in the confirmed list of disease associated loci, are significantly enriched for association to the phenotypes in question in extended GWAS analysis. Finally, we test our method in 3 non-immune traits to assess its applicability to complex traits in general. We find that genes in loci associated to height and lipid levels assemble into significantly connected networks but did not detect excess connectivity among Type 2 Diabetes (T2D) loci beyond chance. Taken together, our results constitute evidence that, for many of the complex diseases studied here, common genetic associations implicate regions encoding proteins that physically interact in a preferential manner, in line with observations in Mendelian disease.

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.

An unusual splice defect in the mitofusin 2 gene (MFN2) is associated with degenerative axonopathy in Tyrolean Grey cattle

Tyrolean Grey cattle represent a local breed with a population size of approximately 5000 registered cows. In 2003, a previously unknown neurological disorder was recognized in Tyrolean Grey cattle. The clinical signs of the disorder are similar to those of bovine progressive degenerative myeloencephalopathy (weaver syndrome) in Brown Swiss cattle but occur much earlier in life. The neuropathological investigation of an affected calf showed axonal degeneration in the central nervous system (CNS) and femoral nerve. The pedigrees of the affected calves suggested a monogenic autosomal recessive inheritance. We localized the responsible mutation to a 1.9 Mb interval on chromosome 16 by genome-wide association and haplotype mapping. The MFN2 gene located in this interval encodes mitofusin 2, a mitochondrial membrane protein. A heritable human axonal neuropathy, Charcot-Marie-Tooth disease-2A2 (CMT2A2), is caused by MFN2 mutations. Therefore, we considered MFN2 a positional and functional candidate gene and performed mutation analysis in affected and control Tyrolean Grey cattle. We did not find any non-synonymous variants. However, we identified a perfectly associated silent SNP in the coding region of exon 20 of the MFN2 gene. This SNP is located within a putative exonic splice enhancer (ESE) and the variant allele leads to partial retention of the entire intron 19 and a premature stop codon in the aberrant MFN2 transcript. Thus we have identified a highly unusual splicing defect, where an exonic single base exchange leads to the retention of the preceding intron. This splicing defect represents a potential explanation for the observed degenerative axonopathy. Marker assisted selection can now be used to eliminate degenerative axonopathy from Tyrolean Grey cattle.

A locus on chromosome 5 is associated with dilated cardiomyopathy in Doberman Pinschers

Dilated cardiomyopathy (DCM) is a heterogeneous group of heart diseases with a strong genetic background. Currently, many human DCM cases exist where no causative mutation can be identified. DCM also occurs with high prevalence in several large dog breeds. In the Doberman Pinscher a specific DCM form characterized by arrhythmias and/or echocardiographic changes has been intensively studied by veterinary cardiologists. We performed a genome-wide association study in Doberman Pinschers. Using 71 cases and 70 controls collected in Germany we identified a genome-wide significant association to DCM on chromosome 5. We validated the association in an independent cohort collected in the United Kingdom. There is no known DCM candidate gene under the association signal. Therefore, DCM in Doberman Pinschers offers the chance of identifying a novel DCM gene that might also be relevant for human health.

LGI2 truncation causes a remitting focal epilepsy in dogs

One quadrillion synapses are laid in the first two years of postnatal construction of the human brain, which are then pruned until age 10 to 500 trillion synapses composing the final network. Genetic epilepsies are the most common neurological diseases with onset during pruning, affecting 0.5% of 2-10-year-old children, and these epilepsies are often characterized by spontaneous remission. We previously described a remitting epilepsy in the Lagotto romagnolo canine breed. Here, we identify the gene defect and affected neurochemical pathway. We reconstructed a large Lagotto pedigree of around 34 affected animals. Using genome-wide association in 11 discordant sib-pairs from this pedigree, we mapped the disease locus to a 1.7 Mb region of homozygosity in chromosome 3 where we identified a protein-truncating mutation in the Lgi2 gene, a homologue of the human epilepsy gene LGI1. We show that LGI2, like LGI1, is neuronally secreted and acts on metalloproteinase-lacking members of the ADAM family of neuronal receptors, which function in synapse remodeling, and that LGI2 truncation, like LGI1 truncations, prevents secretion and ADAM interaction. The resulting epilepsy onsets at around seven weeks (equivalent to human two years), and remits by four months (human eight years), versus onset after age eight in the majority of human patients with LGI1 mutations. Finally, we show that Lgi2 is expressed highly in the immediate post-natal period until halfway through pruning, unlike Lgi1, which is expressed in the latter part of pruning and beyond. LGI2 acts at least in part through the same ADAM receptors as LGI1, but earlier, ensuring electrical stability (absence of epilepsy) during pruning years, preceding this same function performed by LGI1 in later years. LGI2 should be considered a candidate gene for common remitting childhood epilepsies, and LGI2-to-LGI1 transition for mechanisms of childhood epilepsy remission.

Common obesity risk alleles in childhood attention-deficit/hyperactivity disorder

Children with attention-deficit/hyperactivity disorder (ADHD) have a higher rate of obesity than children without ADHD. Obesity risk alleles may overlap with those relevant for ADHD. We examined whether risk alleles for an increased body mass index (BMI) are associated with ADHD and related quantitative traits (inattention and hyperactivity/impulsivity). We screened 32 obesity risk alleles of single nucleotide polymorphisms (SNPs) in a genome-wide association study (GWAS) for ADHD based on 495 patients and 1,300 population-based controls and performed in silico analyses of the SNPs in an ADHD meta-analysis comprising 2,064 trios, 896 independent cases, and 2,455 controls. In the German sample rs206936 in the NUDT3 gene (nudix; nucleoside diphosphate linked moiety X-type motif 3) was associated with ADHD risk (OR: 1.39; P = 3.4 × 10(-4) ; Pcorr  = 0.01). In the meta-analysis data we found rs6497416 in the intronic region of the GPRC5B gene (G protein-coupled receptor, family C, group 5, member B; P = 7.2 × 10(-4) ; Pcorr  = 0.02) as a risk allele for ADHD. GPRC5B belongs to the metabotropic glutamate receptor family, which has been implicated in the etiology of ADHD. In the German sample rs206936 (NUDT3) and rs10938397 in the glucosamine-6-phosphate deaminase 2 gene (GNPDA2) were associated with inattention, whereas markers in the mitogen-activated protein kinase 5 gene (MAP2K5) and in the cell adhesion molecule 2 gene (CADM2) were associated with hyperactivity. In the meta-analysis data, MAP2K5 was associated with inattention, GPRC5B with hyperactivity/impulsivity and inattention and CADM2 with hyperactivity/impulsivity. Our results justify further research on the elucidation of the common genetic background of ADHD and obesity.

DCLK1 variants are associated across schizophrenia and attention deficit/hyperactivity disorder

Doublecortin and calmodulin like kinase 1 (DCLK1) is implicated in synaptic plasticity and neurodevelopment. Genetic variants in DCLK1 are associated with cognitive traits, specifically verbal memory and general cognition. We investigated the role of DCLK1 variants in three psychiatric disorders that have neuro-cognitive dysfunctions: schizophrenia (SCZ), bipolar affective disorder (BP) and attention deficit/hyperactivity disorder (ADHD). We mined six genome wide association studies (GWASs) that were available publically or through collaboration; three for BP, two for SCZ and one for ADHD. We also genotyped the DCLK1 region in additional samples of cases with SCZ, BP or ADHD and controls that had not been whole-genome typed. In total, 9895 subjects were analysed, including 5308 normal controls and 4,587 patients (1,125 with SCZ, 2,496 with BP and 966 with ADHD). Several DCLK1 variants were associated with disease phenotypes in the different samples. The main effect was observed for rs7989807 in intron 3, which was strongly associated with SCZ alone and even more so when cases with SCZ and ADHD were combined (P-value = 4 × 10(-5) and 4 × 10(-6), respectively). Associations were also observed with additional markers in intron 3 (combination of SCZ, ADHD and BP), intron 19 (SCZ+BP) and the 3'UTR (SCZ+BP). Our results suggest that genetic variants in DCLK1 are associated with SCZ and, to a lesser extent, with ADHD and BP. Interestingly the association is strongest when SCZ and ADHD are considered together, suggesting common genetic susceptibility. Given that DCLK1 variants were previously found to be associated with cognitive traits, these results are consistent with the role of DCLK1 in neurodevelopment and synaptic plasticity.

Mapping genetic variants associated with beta-adrenergic responses in inbred mice

β-blockers and β-agonists are primarily used to treat cardiovascular diseases. Inter-individual variability in response to both drug classes is well recognized, yet the identity and relative contribution of the genetic players involved are poorly understood. This work is the first genome-wide association study (GWAS) addressing the values and susceptibility of cardiovascular-related traits to a selective β(1)-blocker, Atenolol (ate), and a β-agonist, Isoproterenol (iso). The phenotypic dataset consisted of 27 highly heritable traits, each measured across 22 inbred mouse strains and four pharmacological conditions. The genotypic panel comprised 79922 informative SNPs of the mouse HapMap resource. Associations were mapped by Efficient Mixed Model Association (EMMA), a method that corrects for the population structure and genetic relatedness of the various strains. A total of 205 separate genome-wide scans were analyzed. The most significant hits include three candidate loci related to cardiac and body weight, three loci for electrocardiographic (ECG) values, two loci for the susceptibility of atrial weight index to iso, four loci for the susceptibility of systolic blood pressure (SBP) to perturbations of the β-adrenergic system, and one locus for the responsiveness of QTc (p<10(-8)). An additional 60 loci were suggestive for one or the other of the 27 traits, while 46 others were suggestive for one or the other drug effects (p<10(-6)). Most hits tagged unexpected regions, yet at least two loci for the susceptibility of SBP to β-adrenergic drugs pointed at members of the hypothalamic-pituitary-thyroid axis. Loci for cardiac-related traits were preferentially enriched in genes expressed in the heart, while 23% of the testable loci were replicated with datasets of the Mouse Phenome Database (MPD). Altogether these data and validation tests indicate that the mapped loci are relevant to the traits and responses studied.

Serum ferritin levels are associated with a distinct phenotype of chronic hepatitis C poorly responding to pegylated interferon-alpha and ribavirin therapy

Elevated serum ferritin levels may reflect a systemic inflammatory state as well as increased iron storage, both of which may contribute to an unfavorable outcome of chronic hepatitis C (CHC). We therefore performed a comprehensive analysis of the role of serum ferritin and its genetic determinants in the pathogenesis and treatment of CHC. To this end, serum ferritin levels at baseline of therapy with pegylated interferon-alpha and ribavirin or before biopsy were correlated with clinical and histological features of chronic hepatitis C virus (HCV) infection, including necroinflammatory activity (N = 970), fibrosis (N = 980), steatosis (N = 886), and response to treatment (N = 876). The association between high serum ferritin levels (> median) and the endpoints was assessed by logistic regression. Moreover, a candidate gene as well as a genome-wide association study of serum ferritin were performed. We found that serum ferritin ≥ the sex-specific median was one of the strongest pretreatment predictors of treatment failure (univariate P < 0.0001, odds ratio [OR] = 0.45, 95% confidence interval [CI] = 0.34-0.60). This association remained highly significant in a multivariate analysis (P = 0.0002, OR = 0.35, 95% CI = 0.20-0.61), with an OR comparable to that of interleukin (IL)28B genotype. When patients with the unfavorable IL28B genotypes were stratified according to high versus low ferritin levels, SVR rates differed by > 30% in both HCV genotype 1- and genotype 3-infected patients (P < 0.001). Serum ferritin levels were also independently associated with severe liver fibrosis (P < 0.0001, OR = 2.67, 95% CI = 1.68-4.25) and steatosis (P = 0.002, OR = 2.29, 95% CI = 1.35-3.91), but not with necroinflammatory activity (P = 0.3). Genetic variations had only a limited impact on serum ferritin levels. Conclusion: In patients with CHC, elevated serum ferritin levels are independently associated with advanced liver fibrosis, hepatic steatosis, and poor response to interferon-alpha-based therapy.

Bipolar disorder risk alleles in children with ADHD

Bipolar disorder (BD) and attention deficit/hyperactivity disorder (ADHD) may share common genetic risk factors as indicated by the high co-morbidity of BD and ADHD, their phenotypic overlap especially in pediatric populations, the high heritability of both disorders, and the co-occurrence in families. We therefore examined whether known polygenic BD risk alleles are associated with ADHD. We chose the eight best SNPs of the recent genome-wide association study (GWAS) of BD patients of German ancestry and the nine SNPs from international GWAS meeting a 'genome-wide significance' level of α = 5 × 10(-8). A GWAS was performed in 495 ADHD children and 1,300 population-based controls using HumanHap550v3 and Human660 W-Quadv1 BeadArrays. We found no significant association of childhood ADHD with single BD risk alleles surviving adjustment for multiple testing. Yet, risk alleles for BD and ADHD were directionally consistent at eight of nine loci with the strongest support for three SNPs in or near NCAN, BRE, and LMAN2L. The polygene analysis for the BP risk alleles at all 14 loci indicated a higher probability of being a BD risk allele carrier in the ADHD cases as compared to the controls. At a moderate power to detect association with ADHD, if true effects were close to estimates from GWAS for BD, our results suggest that the possible contribution of BD risk variants to childhood ADHD risk is considerably lower than for BD. Yet, our findings should encourage researchers to search for common genetic risk factors in BD and childhood ADHD in future studies.