A Genome-Wide Association Study Reveals Loci Influencing Height and Other Conformation Traits in Horses

The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ∼200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ∼40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ∼18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits.

Genome-wide association study of survival from sepsis due to pneumonia: an observational cohort study.

BACKGROUND Sepsis continues to be a major cause of death, disability, and health-care expenditure worldwide. Despite evidence suggesting that host genetics can influence sepsis outcomes, no specific loci have yet been convincingly replicated. The aim of this study was to identify genetic variants that influence sepsis survival. METHODS We did a genome-wide association study in three independent cohorts of white adult patients admitted to intensive care units with sepsis, severe sepsis, or septic shock (as defined by the International Consensus Criteria) due to pneumonia or intra-abdominal infection (cohorts 1-3, n=2534 patients). The primary outcome was 28 day survival. Results for the cohort of patients with sepsis due to pneumonia were combined in a meta-analysis of 1553 patients from all three cohorts, of whom 359 died within 28 days of admission to the intensive-care unit. The most significantly associated single nucleotide polymorphisms (SNPs) were genotyped in a further 538 white patients with sepsis due to pneumonia (cohort 4), of whom 106 died. FINDINGS In the genome-wide meta-analysis of three independent pneumonia cohorts (cohorts 1-3), common variants in the FER gene were strongly associated with survival (p=9·7 × 10(-8)). Further genotyping of the top associated SNP (rs4957796) in the additional cohort (cohort 4) resulted in a combined p value of 5·6 × 10(-8) (odds ratio 0·56, 95% CI 0·45-0·69). In a time-to-event analysis, each allele reduced the mortality over 28 days by 44% (hazard ratio for death 0·56, 95% CI 0·45-0·69; likelihood ratio test p=3·4 × 10(-9), after adjustment for age and stratification by cohort). Mortality was 9·5% in patients carrying the CC genotype, 15·2% in those carrying the TC genotype, and 25·3% in those carrying the TT genotype. No significant genetic associations were identified when patients with sepsis due to pneumonia and intra-abdominal infection were combined. INTERPRETATION We have identified common variants in the FER gene that associate with a reduced risk of death from sepsis due to pneumonia. The FER gene and associated molecular pathways are potential novel targets for therapy or prevention and candidates for the development of biomarkers for risk stratification. FUNDING European Commission and the Wellcome Trust.

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency.

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.

A Multi-Breed Genome-Wide Association Analysis for Canine Hypothyroidism Identifies a Shared Major Risk Locus on CFA12.

Hypothyroidism is a complex clinical condition found in both humans and dogs, thought to be caused by a combination of genetic and environmental factors. In this study we present a multi-breed analysis of predisposing genetic risk factors for hypothyroidism in dogs using three high-risk breeds-the Gordon Setter, Hovawart and the Rhodesian Ridgeback. Using a genome-wide association approach and meta-analysis, we identified a major hypothyroidism risk locus shared by these breeds on chromosome 12 (p = 2.1x10-11). Further characterisation of the candidate region revealed a shared ~167 kb risk haplotype (4,915,018-5,081,823 bp), tagged by two SNPs in almost complete linkage disequilibrium. This breed-shared risk haplotype includes three genes (LHFPL5, SRPK1 and SLC26A8) and does not extend to the dog leukocyte antigen (DLA) class II gene cluster located in the vicinity. These three genes have not been identified as candidate genes for hypothyroid disease previously, but have functions that could potentially contribute to the development of the disease. Our results implicate the potential involvement of novel genes and pathways for the development of canine hypothyroidism, raising new possibilities for screening, breeding programmes and treatments in dogs. This study may also contribute to our understanding of the genetic etiology of human hypothyroid disease, which is one of the most common endocrine disorders in humans.

Salivary Pellicle vs Whole Saliva: The Response of Oral Fibroblasts Based on a Genome-Wide Microarray.

PURPOSE Whole saliva comprises components of the salivary pellicle that spontaneously forms on surfaces of implants and teeth. However, there are no studies that functionally link the salivary pellicle with a possible change in gene expression. MATERIALS AND METHODS This study examined the genetic response of oral fibroblasts exposed to the salivary pellicle and whole saliva. Oral fibroblasts were seeded onto a salivary pellicle and the respective untreated surface. Oral fibroblasts were also exposed to freshly harvested sterile-filtered whole saliva. A genome-wide microarray of oral fibroblasts was performed, followed by gene ontology screening with DAVID functional annotation clustering, KEGG pathway analysis, and the STRING functional protein association network. RESULTS Exposure of oral fibroblasts to saliva caused 61 genes to be differentially expressed (P < .05). Gene ontology screening assigned the respective genes into 262 biologic processes, 3 cellular components, 13 molecular functions, and 7 pathways. Most remarkable was the enrichment in the inflammatory response. None of the genes regulated by whole saliva was significantly changed when cells were placed onto a salivary pellicle. CONCLUSION The salivary pellicle per se does not provoke a significant inflammatory response of oral fibroblasts in vitro, whereas sterile-filtered whole saliva does produce a strong inflammatory response.

Detection of convergent genome-wide signals of adaptation to tropical forests in humans

Tropical forests are believed to be very harsh environments for human life. It is unclear whether human beings would have ever subsisted in those environments without external resources. It is therefore possible that humans have developed recent biological adaptations in response to specific selective pressures to cope with this challenge. To understand such biological adaptations we analyzed genome-wide SNP data under a Bayesian statistics framework, looking for outlier markers with an overly large extent of differentiation between populations living in a tropical forest, as compared to genetically related populations living outside the forest in Africa and the Americas. The most significant positive selection signals were found in genes related to lipid metabolism, the immune system, body development, and RNA Polymerase III transcription initiation. The results are discussed in the light of putative tropical forest selective pressures, namely food scarcity, high prevalence of pathogens, difficulty to move, and inefficient thermoregulation. Agreement between our results and previous studies on the pygmy phenotype, a putative prototype of forest adaptation, were found, suggesting that a few genetic regions previously described as associated with short stature may be evolving under similar positive selection in Africa and the Americas. In general, convergent evolution was less pervasive than local adaptation in one single continent, suggesting that Africans and Amerindians may have followed different routes to adapt to similar environmental selective pressures.

Genome-wide association studies based on sequence-derived genotypes reveal new QTL associated with conformation and performance traits in the Franches-Montagnes horse breed

To identify novel quantitative trait loci (QTL) within horses, we performed genome-wide association studies (GWAS) based on sequence-level genotypes for conformation and performance traits in the Franches-Montagnes (FM) horse breed. Sequence-level genotypes of FM horses were derived by re-sequencing 30 key founders and imputing 50K data of genotyped horses. In total, we included 1077 FM horses genotyped for ~4 million SNPs and their respective de-regressed breeding values of the traits in the analysis. Based on this dataset, we identified a total of 14 QTL associated with 18 conformation traits and one performance trait. Therefore, our results suggest that the application of sequence-derived genotypes increases the power to identify novel QTL which were not identified previously based on 50K SNP chip data.

Genome-wide Association Studies Identify Genetic Loci Associated with Albuminuria in Diabetes

Elevated concentrations of albumin in the urine, albuminuria, are a hallmark of diabetic kidney disease and associate with increased risk for end-stage renal disease and cardiovascular events. To gain insight into the pathophysiological mechanisms underlying albuminuria, we conducted meta-analyses of genome-wide association studies and independent replication in up to 5,825 individuals of European ancestry with diabetes mellitus and up to 46,061 without diabetes, followed by functional studies. Known associations of variants in CUBN, encoding cubilin, with the urinary albumin-to-creatinine ratio (UACR) were confirmed in the overall sample (p=2.4*10(-10)). Gene-by-diabetes interactions were detected and confirmed for variants in HS6ST1 and near RAB38/CTSC. SNPs at these loci demonstrated a genetic effect on UACR in individuals with but not without diabetes. The change in average UACR per minor allele was 21% for HS6ST1 and 13% for RAB38/CTSC (p=6.3*10(-7) and 5.8*10(-7), respectively). Experiments using streptozotocin-treated diabetic Rab38 knockout and control rats showed higher urinary albumin concentrations and reduced amounts of megalin and cubilin at the proximal tubule cell surface in Rab38 knockout vs. control rats. Relative expression of RAB38 was higher in tubuli of patients with diabetic kidney disease compared to controls. The loci identified here confirm known and highlight novel pathways influencing albuminuria.

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.

A genome-wide association study confirms PNPLA3 and identifies TM6SF2 and MBOAT7 as risk loci for alcohol-related cirrhosis

Alcohol misuse is the leading cause of cirrhosis and the second most common indication for liver transplantation in the Western world. We performed a genome-wide association study for alcohol-related cirrhosis in individuals of European descent (712 cases and 1,426 controls) with subsequent validation in two independent European cohorts (1,148 cases and 922 controls). We identified variants in the MBOAT7 (P = 1.03 × 10(-9)) and TM6SF2 (P = 7.89 × 10(-10)) genes as new risk loci and confirmed rs738409 in PNPLA3 as an important risk locus for alcohol-related cirrhosis (P = 1.54 × 10(-48)) at a genome-wide level of significance. These three loci have a role in lipid processing, suggesting that lipid turnover is important in the pathogenesis of alcohol-related cirrhosis.

Mapping of a new locus for congenital anomalies of the kidney and urinary tract on chromosome 8q24

Congenital anomalies of the kidney and urinary tract (CAKUT) account for the majority of end-stage renal disease in children (50%). Previous studies have mapped autosomal dominant loci for CAKUT. We here report a genome-wide search for linkage in a large pedigree of Somalian descent containing eight affected individuals with a non-syndromic form of CAKUT.

Characterization of the CopR regulon of Lactococcus lactis IL1403

To identify components of the copper homeostatic mechanism of Lactococcus lactis, we employed two-dimensional gel electrophoresis to detect changes in the proteome in response to copper. Three proteins upregulated by copper were identified: glyoxylase I (YaiA), a nitroreductase (YtjD), and lactate oxidase (LctO). The promoter regions of these genes feature cop boxes of consensus TACAnnTGTA, which are the binding site of CopY-type copper-responsive repressors. A genome-wide search for cop boxes revealed 28 such sequence motifs. They were tested by electrophoretic mobility shift assays for the interaction with purified CopR, the CopY-type repressor of L. lactis. Seven of the cop boxes interacted with CopR in a copper-sensitive manner. They were present in the promoter region of five genes, lctO, ytjD, copB, ydiD, and yahC; and two polycistronic operons, yahCD-yaiAB and copRZA. Induction of these genes by copper was confirmed by real-time quantitative PCR. The copRZA operon encodes the CopR repressor of the regulon; a copper chaperone, CopZ; and a putative copper ATPase, CopA. When expressed in Escherichia coli, the copRZA operon conferred copper resistance, suggesting that it functions in copper export from the cytoplasm. Other member genes of the CopR regulon may similarly be involved in copper metabolism.

LUPA: a European initiative taking advantage of the canine genome architecture for unravelling complex disorders in both human and dogs

The domestic dog offers a unique opportunity to explore the genetic basis of disease, morphology and behaviour. Humans share many diseases with our canine companions, making dogs an ideal model organism for comparative disease genetics. Using newly developed resources, genome-wide association studies in dog breeds are proving to be exceptionally powerful. Towards this aim, veterinarians and geneticists from 12 European countries are collaborating to collect and analyse the DNA from large cohorts of dogs suffering from a range of carefully defined diseases of relevance to human health. This project, named LUPA, has already delivered considerable results. The consortium has collaborated to develop a new high density single nucleotide polymorphism (SNP) array. Mutations for four monogenic diseases have been identified and the information has been utilised to find mutations in human patients. Several complex diseases have been mapped and fine mapping is underway. These findings should ultimately lead to a better understanding of the molecular mechanisms underlying complex diseases in both humans and their best friend.

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.