Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

Genome-wide analysis reveals selection for important traits in domestic horse breeds.

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F(ST)-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse.

Major histocompatibility complex and other allergy-related candidate genes associated with insect bite hypersensitivity in Icelandic horses

Insect bite hypersensitivity (IBH) is an allergic dermatitis of horses caused by bites of insects. IBH is a multifactorial disease with contribution of genetic and environmental factors. Candidate gene association analysis of IBH was performed in a group of 89 Icelandic horses all born in Iceland and imported to Europe. Horses were classified in IBH-affected and non-affected based on clinical signs and history of recurrent dermatitis, and on the results of an in vitro sulfidoleukotriene (sLT)-release assay with Culicoides nubeculosus and Simulium vittatum extract. Different genetic markers were tested for association with IBH by the Fisher's exact test. The effect of the major histocompatibility complex (MHC) gene region was studied by genotyping five microsatellites spanning the MHC region (COR112, COR113, COR114, UM011 and UMN-JH34-2), and exon 2 polymorphisms of the class II Eqca-DRA gene. Associations with Eqca-DRA and COR113 were identified (p < 0.05). In addition, a panel of 20 single nucleotide polymorphisms (SNPs) in 17 candidate allergy-related genes was tested. During the initial screen, no marker from the panel was significantly (p < 0.05) associated with IBH. Five SNPs associated with IBH at p < 0.10 were therefore used for analysis of combined genotypes. Out of them, SNPs located in the genes coding for the CD14 receptor (CD14), interleukin 23 receptor (IL23R), thymic stromal lymphopoietin (TSLP) and transforming growth factor beta 3 (TGFB3) molecules were associated with IBH as parts of complex genotypes. These results are supported by similar associations and by expression data from different horse populations and from human studies.

Surveillance of bovine tuberculosis and risk estimation of a future reservoir formation in wildlife in Switzerland and Liechtenstein.

Bovine tuberculosis (bTB) caused by Mycobacterium bovis or M. caprae has recently (re-) emerged in livestock and wildlife in all countries bordering Switzerland (CH) and the Principality of Liechtenstein (FL). Comprehensive data for Swiss and Liechtenstein wildlife are not available so far, although two native species, wild boar (Sus scrofa) and red deer (Cervus elaphus elaphus), act as bTB reservoirs elsewhere in continental Europe. Our aims were (1) to assess the occurrence of bTB in these wild ungulates in CH/FL and to reinforce scanning surveillance in all wild mammals; (2) to evaluate the risk of a future bTB reservoir formation in wild boar and red deer in CH/FL. Tissue samples collected from 2009 to 2011 from 434 hunted red deer and wild boar and from eight diseased ungulates with tuberculosis-like lesions were tested by direct real-time PCR and culture to detect mycobacteria of the Mycobacterium tuberculosis complex (MTBC). Identification of suspicious colonies was attempted by real-time PCR, genotyping and spoligotyping. Information on risk factors for bTB maintenance within wildlife populations was retrieved from the literature and the situation regarding identified factors was assessed for our study areas. Mycobacteria of the MTBC were detected in six out of 165 wild boar (3.6%; 95% CI: 1.4-7.8) but none of the 269 red deer (0%; 0-1.4). M. microti was identified in two MTBC-positive wild boar, while species identification remained unsuccessful in four cases. Main risk factors for bTB maintenance worldwide, including different causes of aggregation often resulting from intensive wildlife management, are largely absent in CH and FL. In conclusion, M. bovis and M. caprae were not detected but we report for the first time MTBC mycobacteria in Swiss wild boar. Present conditions seem unfavorable for a reservoir emergence, nevertheless increasing population numbers of wild ungulates and offal consumption may represent a risk.

Adaptive evolution during an ongoing range expansion: the invasive bank vole (Myodes glareolus) in Ireland

Range expansions are extremely common, but have only recently begun to attract attention in terms of their genetic consequences. As populations expand, demes at the wave front experience strong genetic drift, which is expected to reduce genetic diversity and potentially cause ‘allele surfing’, where alleles may become fixed over a wide geographical area even if their effects are deleterious. Previous simulation models show that range expansions can generate very strong selective gradients on dispersal, reproduction, competition and immunity. To investigate the effects of range expansion on genetic diversity and adaptation, we studied the population genomics of the bank vole (Myodes glareolus) in Ireland. The bank vole was likely introduced in the late 1920s and is expanding its range at a rate of ~2.5 km/year. Using genotyping-by-sequencing, we genotyped 281 bank voles at 5979 SNP loci. Fourteen sample sites were arranged in three transects running from the introduction site to the wave front of the expansion. We found significant declines in genetic diversity along all three transects. However, there was no evidence that sites at the wave front had accumulated more deleterious mutations. We looked for outlier loci with strong correlations between allele frequency and distance from the introduction site, where the direction of correlation was the same in all three transects. Amongst these outliers, we found significant enrichment for genic SNPs, suggesting the action of selection. Candidates for selection included several genes with immunological functions and several genes that could influence behaviour.

Life-threatening adverse events following therapeutic opioid administration in adults: Is pharmacogenetic analysis useful?

BACKGROUND Systemic approaches are needed to understand how variations in the genes associated with opioid pharmacokinetics and response can be used to predict patient outcome. The application of pharmacogenetic analysis to two cases of life-threatening opioid-induced respiratory depression is presented. The usefulness of genotyping in the context of these cases is discussed. METHODS A panel of 20 functional candidate polymorphisms in genes involved in the opioid biotransformation pathway (CYP2D6, UGT2B7, ABCB1, OPRM1, COMT) were genotyped in these two patients using commercially available genotyping assays. RESULTS In case 1, the patient experienced adverse outcomes when administered codeine and morphine, but not hydromorphone. Genetic test results suggested that this differential response may be due to an inherent propensity to generate active metabolites from both codeine and morphine. These active metabolites are not generated with hydromorphone. In case 2, the patient experienced severe respiratory depression during postoperative recovery following standard doses of morphine. The patient was found to carry genetic variations that result in decreased morphine efflux transporter activity at the blood-brain barrier and increased sensitivity to opioids. CONCLUSIONS Knowledge of the relative contribution of pharmacogenetic biomarkers and their influence on opioid response are continually evolving. Pharmacogenetic analysis, together with clinical history, has the potential to provide mechanistic insight into severe respiratory depressive events in patients who receive opioids at therapeutic doses.

DQB1 locus alone explains most of the risk and protection in narcolepsy with cataplexy in Europe

STUDY OBJECTIVE Prior research has identified five common genetic variants associated with narcolepsy with cataplexy in Caucasian patients. To replicate and/or extend these findings, we have tested HLA-DQB1, the previously identified 5 variants, and 10 other potential variants in a large European sample of narcolepsy with cataplexy subjects. DESIGN Retrospective case-control study. SETTING A recent study showed that over 76% of significant genome-wide association variants lie within DNase I hypersensitive sites (DHSs). From our previous GWAS, we identified 30 single nucleotide polymorphisms (SNPs) with P < 10(-4) mapping to DHSs. Ten SNPs tagging these sites, HLADQB1, and all previously reported SNPs significantly associated with narcolepsy were tested for replication. PATIENTS AND PARTICIPANTS For GWAS, 1,261 narcolepsy patients and 1,422 HLA-DQB1*06:02-matched controls were included. For HLA study, 1,218 patients and 3,541 controls were included. MEASUREMENTS AND RESULTS None of the top variants within DHSs were replicated. Out of the five previously reported SNPs, only rs2858884 within the HLA region (P < 2x10(-9)) and rs1154155 within the TRA locus (P < 2x10(-8)) replicated. DQB1 typing confirmed that DQB1*06:02 confers an extraordinary risk (odds ratio 251). Four protective alleles (DQB1*06:03, odds ratio 0.17, DQB1*05:01, odds ratio 0.56, DQB1*06:09 odds ratio 0.21, DQB1*02 odds ratio 0.76) were also identified. CONCLUSION An overwhelming portion of genetic risk for narcolepsy with cataplexy is found at DQB1 locus. Since DQB1*06:02 positive subjects are at 251-fold increase in risk for narcolepsy, and all recent cases of narcolepsy after H1N1 vaccination are positive for this allele, DQB1 genotyping may be relevant to public health policy.

IL28B expression depends on a novel TT/-G polymorphism which improves HCV clearance prediction

Approximately 3% of the world population is chronically infected with the hepatitis C virus (HCV), with potential development of cirrhosis and hepatocellular carcinoma. Despite the availability of new antiviral agents, treatment remains suboptimal. Genome-wide association studies (GWAS) identified rs12979860, a polymorphism nearby IL28B, as an important predictor of HCV clearance. We report the identification of a novel TT/-G polymorphism in the CpG region upstream of IL28B, which is a better predictor of HCV clearance than rs12979860. By using peripheral blood mononuclear cells (PBMCs) from individuals carrying different allelic combinations of the TT/-G and rs12979860 polymorphisms, we show that induction of IL28B and IFN-γ–inducible protein 10 (IP-10) mRNA relies on TT/-G, but not rs12979860, making TT/-G the only functional variant identified so far. This novel step in understanding the genetic regulation of IL28B may have important implications for clinical practice, as the use of TT/G genotyping instead of rs12979860 would improve patient management.

Tracking a tuberculosis outbreak over 21 years: strain-specific single nucleotide polymorphism-typing combined with targeted whole genome sequencing.

BACKGROUND  Whole genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single nucleotide polymorphism (SNP)-typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. METHODS  Based on genome sequences of three historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1,642 patient isolates, and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. RESULTS  We identified 68 patients associated with the outbreak strain. Most were diagnosed in 1991-1995, but cases were observed until 2011. Two thirds belonged to the homeless and substance abuser milieu. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into three sub-clusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. CONCLUSIONS  Strain-specific SNP-genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real-time and at high resolution.

The brown coat colour of Coppernecked goats is associated with a non-synonymous variant at the TYRP1 locus on chromosome 8

The recent development of a goat SNP genotyping microarray enables genome-wide association studies in this important livestock species. We investigated the genetic basis of the black and brown coat colour in Valais Blacknecked and Coppernecked goats. A genome-wide association analysis using goat SNP50 BeadChip genotypes of 22 cases and 23 controls allowed us to map the locus for the brown coat colour to goat chromosome 8. The TYRP1 gene is located within the associated chromosomal region, and TYRP1 variants cause similar coat colour phenotypes in different species. We thus considered TYRP1 as a strong positional and functional candidate. We resequenced the caprine TYRP1 gene by Sanger and Illumina sequencing and identified two non-synonymous variants, p.Ile478Thr and p.Gly496Asp, that might have a functional impact on the TYRP1 protein. However, based on the obtained pedigree and genotype data, the brown coat colour in these goats is not due to a single recessive loss-of-function allele. Surprisingly, the genotype distribution and the pedigree data suggest that the (496) Asp allele might possibly act in a dominant manner. The (496) Asp allele was present in 77 of 81 investigated Coppernecked goats and did not occur in black goats. This strongly suggests heterogeneity underlying the brown coat colour in Coppernecked goats. Functional experiments or targeted matings will be required to verify the unexpected preliminary findings.

Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

Bacteriological diagnosis and molecular strain typing of Mycobacterium bovis and Mycobacterium caprae

The primary isolation of a Mycobacterium sp. of the Mycobacterium tuberculosis complex from an infected animal provides a definitive diagnosis of tuberculosis. However, as Mycobacterium bovis and Mycobacterium caprae are difficult to isolate, particularly for animals in the early stages of disease, success is dependent on the optimal performance of all aspects of the bacteriological process, from the initial choice of tissue samples at post-mortem examination or clinical samples, to the type of media and conditions used to cultivate the microorganism. Each step has its own performance characteristics, which can contribute to sensitivity and specificity of the procedure, and may need to be optimized in order to achieve the gold standard diagnosis. Having isolated the slow-growing mycobacteria, species identification and fine resolution strain typing are keys to understanding the epidemiology of the disease and to devise strategies to limit transmission of infection. New technologies have emerged that can now even discriminate different isolates from the same animal. In this review we highlight the key factors that contribute to the accuracy of bacteriological diagnosis of M. bovis and M. caprae, and describe the development of advanced genotyping techniques that are increasingly used in diagnostic laboratories for the purpose of supporting detailed epidemiological investigations.

A history and overview of phenotypic variability in CYP2D6 activity

CYP2D6 is a human cytochrome P450 that is responsible for the metabolism of a large number of drugs and chemicals. Interest in CYP2D6 has largely centered on the wide interindividual variability in its catalytic activity that stems from a common genetic polymorphism in the CYP2D6 gene. Two major phenotypes exist, extensive metabolizer (EM) and poor metabolizer (PM), together with the two less studied phenotypes of ultrarapid metabolizer (UM) and intermediate metabolizer. These phenotypes are the expression of an underlying allelomorphism in CYP2D6 and are also context dependent. Several drugs that are CYP2D6 substrates display polymorphic metabolism, that is, the existence in the population of multiple phenotypes, in particular EM and PM. The most notable drugs in this regard are debrisoquine and sparteine, although there are also data for a few others, in particular, dextromethorphan and metoprolol. Many nongenetic factors can alter the expression of CYP2D6 phenotypes, the most significant of which is the presence of other drugs. In this context, the EM phenotype may not be immutable, with potential conversion into a PM phenocopy, due to significantly impaired CYP2D6 metabolism in the presence of other CYP2D6 substrates and inhibitors. This phenotype interconversion generated great concern and helped drive the movement away from phenotyping based upon drug administration to genotyping of acquired DNA samples. However, ascertaining the presence of CYP2D6 alleles in a DNA sample does not determine the metabolism and pharmacokinetics of CYP2D6 substrates in that subject: it is a forecast, much like the weather forecast and, as we all know regarding the weather, the forecast can be inaccurate at times.

Carriage of the PNPLA3 rs738409 C >G polymorphism confers an increased risk of non-alcoholic fatty liver disease associated hepatocellular carcinoma.

BACKGROUND & AIMS Subtle inter-patient genetic variation and environmental factors combine to determine disease progression in non-alcoholic fatty liver disease (NAFLD). Carriage of the PNPLA3 rs738409 c.444C >G minor allele (encoding the I148M variant) has been robustly associated with advanced NAFLD. Although most hepatocellular carcinoma (HCC) is related to chronic viral hepatitis or alcoholic liver disease, the incidence of NAFLD-related HCC is increasing. We examined whether rs738409 C >G was associated with HCC-risk in patients with NAFLD. METHODS PNPLA3 rs738409 genotype was determined by allelic discrimination in 100 European Caucasians with NAFLD-related HCC and 275 controls with histologically characterised NAFLD. RESULTS Genotype frequencies were significantly different between NAFLD-HCC cases (CC=28, CG=43, GG=29) and NAFLD-controls (CC=125, CG=117, GG=33) (p=0.0001). In multivariate analysis adjusted for age, gender, diabetes, BMI, and presence of cirrhosis, carriage of each copy of the rs738409 minor (G) allele conferred an additive risk for HCC (adjusted OR 2.26 [95% CI 1.23-4.14], p=0.0082), with GG homozygotes exhibiting a 5-fold [1.47-17.29], p=0.01 increased risk over CC. When compared to the UK general population (1958 British Birth Cohort, n=1476), the risk-effect was more pronounced (GC vs. CC: unadjusted OR 2.52 [1.55-4.10], p=0.0002; GG vs. CC: OR 12.19 [6.89-21.58], p<0.0001). CONCLUSIONS Carriage of the PNPLA3 rs738409 C >G polymorphism is not only associated with greater risk of progressive steatohepatitis and fibrosis but also of HCC. If validated, these findings suggest that PNPLA3 genotyping has the potential to contribute to multi-factorial patient-risk stratification, identifying those to whom HCC surveillance may be targeted.

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.