Fine-scale population structure analysis of seven local Swiss sheep breeds using genome-wide SNP data

As part of the global sheep Hapmap project, 24 individuals from each of seven indigenous Swiss sheep breeds (Bundner Oberländer sheep (BOS), Engadine Red sheep (ERS), Swiss Black-Brown Mountain sheep (SBS), Swiss Mirror sheep (SMS), Swiss White Alpine (SWA) sheep, Valais Blacknose sheep (VBS) and Valais Red sheep (VRS)), were genotyped using Illumina’s Ovine SNP50 BeadChip. In total, 167 animals were subjected to a detailed analysis for genetic diversity using 45 193 informative single nucleotide polymorphisms. The results of the phylogenetic analyses supported the known proximity between populations such as VBS and VRS or SMS and SWA. Average genomic relatedness within a breed was found to be 12 percent (BOS), 5 percent (ERS), 9 percent (SBS), 10 percent (SMS), 9 percent (SWA), 12 percent (VBS) and 20 percent (VRS). Furthermore, genomic relationships between breeds were found for single individuals from SWA and SMS, VRS and VBS as well as VRS and BOS. In addition, seven out of 40 indicated parent–offspring pairs could not be confirmed. These results were further supported by results from the genome-wide population cluster analysis. This study provides a better understanding of fine-scale population structures within and between Swiss sheep breeds. This relevant information will help to increase the conservation activities of the local Swiss sheep breeds.

A Deletion in the VLDLR Gene in Eurasier Dogs with Cerebellar Hypoplasia Resembling a Dandy-Walker-Like Malformation (DWLM).

Dandy-Walker-like malformation (DWLM) is the result of aberrant brain development and mainly characterized by cerebellar hypoplasia. DWLM affected dogs display a non-progressive cerebellar ataxia. Several DWLM cases were recently observed in the Eurasier dog breed, which strongly suggested a monogenic autosomal recessive inheritance in this breed. We performed a genome-wide association study (GWAS) with 9 cases and 11 controls and found the best association of DWLM with markers on chromosome 1. Subsequent homozygosity mapping confirmed that all 9 cases were homozygous for a shared haplotype in this region, which delineated a critical interval of 3.35 Mb. We sequenced the genome of an affected Eurasier and compared it with the Boxer reference genome and 47 control genomes of dogs from other breeds. This analysis revealed 4 private non-synonymous variants in the critical interval of the affected Eurasier. We genotyped these variants in additional dogs and found perfect association for only one of these variants, a single base deletion in the VLDLR gene encoding the very low density lipoprotein receptor. This variant, VLDLR:c.1713delC is predicted to cause a frameshift and premature stop codon (p.W572Gfs*10). Variants in the VLDLR gene have been shown to cause congenital cerebellar ataxia and mental retardation in human patients and Vldlr knockout mice also display an ataxia phenotype. Our combined genetic data together with the functional knowledge on the VLDLR gene from other species thus strongly suggest that VLDLR:c.1713delC is indeed causing DWLM in Eurasier dogs.

VKORC1 and CYP2C9 genotypes are predictors of warfarin-related outcomes in children

BACKGROUND Despite substantial evidence supporting a pharmacogenetic approach to warfarin therapy in adults, evidence on the importance of genetics in warfarin therapy in children is limited, particularly for clinical outcomes. We assessed the contribution of CYP2C9/VKORC1/CYP4F2 genotypes and variation in other genes involved in vitamin K and coagulation pathways to warfarin dose and related clinical outcomes in children. PROCEDURE Clinical and genetic data for 93 children (age ≤ 18 years) who received warfarin therapy were obtained. DNA was genotyped for 93 selected single nucleotide polymorphisms using a custom assay. RESULTS With a median age of 4.8 years, our cohort included more young children than most previous studies. Overall, 76.3% of dose variability was explained by weight, indication, VKORC1-1639G/A and CYP2C9 *2/*3, with genotypes accounting for 21.1% of variability. There was a strong correlation (R(2) = 0.68; P < 0.001) between actual and predicted warfarin dose using a pediatric genotype-based dosing model. VKORC1 genotype had a significant impact on time to therapeutic international normalized ratio (INR) (P = 0.047) and time to over-anticoagulation (INR > 4; P = 0.024) during the initiation of therapy. CYP2C9*3 carriers were also at increased risk of major bleeding while receiving warfarin (adjusted OR = 11.28). An additional variant in CYP2C9 (rs7089580) was significantly associated with warfarin dose (P = 0.020) in a multivariate clinical and genetic model. CONCLUSIONS This study confirms the importance of VKORC1/CYP2C9 genotypes for warfarin dosing in a young pediatric cohort and demonstrates an impact of genetic factors on clinical outcomes in children. Furthermore, we identified an additional variant in CYP2C9 of potential relevance for warfarin dosing in children.

Identification of amino acid substitutions with compensational effects in the attachment protein of canine distemper virus.

The hemagglutinin (H) gene of canine distemper virus (CDV) encodes the receptor-binding protein. This protein, together with the fusion (F) protein, is pivotal for infectivity since it contributes to the fusion of the viral envelope with the host cell membrane. Of the two receptors currently known for CDV (nectin-4 and the signaling lymphocyte activation molecule [SLAM]), SLAM is considered the most relevant for host susceptibility. To investigate how evolution might have impacted the host-CDV interaction, we examined the functional properties of a series of missense single nucleotide polymorphisms (SNPs) naturally accumulating within the H-gene sequences during the transition between two distinct but related strains. The two strains, a wild-type strain and a consensus strain, were part of a single continental outbreak in European wildlife and occurred in distinct geographical areas 2 years apart. The deduced amino acid sequence of the two H genes differed at 5 residues. A panel of mutants carrying all the combinations of the SNPs was obtained by site-directed mutagenesis. The selected mutant, wild type, and consensus H proteins were functionally evaluated according to their surface expression, SLAM binding, fusion protein interaction, and cell fusion efficiencies. The results highlight that the most detrimental functional effects are associated with specific sets of SNPs. Strikingly, an efficient compensational system driven by additional SNPs appears to come into play, virtually neutralizing the negative functional effects. This system seems to contribute to the maintenance of the tightly regulated function of the H-gene-encoded attachment protein. Importance: To investigate how evolution might have impacted the host-canine distemper virus (CDV) interaction, we examined the functional properties of naturally occurring single nucleotide polymorphisms (SNPs) in the hemagglutinin gene of two related but distinct strains of CDV. The hemagglutinin gene encodes the attachment protein, which is pivotal for infection. Our results show that few SNPs have a relevant detrimental impact and they generally appear in specific combinations (molecular signatures). These drastic negative changes are neutralized by compensatory mutations, which contribute to maintenance of an overall constant bioactivity of the attachment protein. This compensational mechanism might reflect the reaction of the CDV machinery to the changes occurring in the virus following antigenic variations critical for virulence.

Impact of common risk factors of fibrosis progression in chronic hepatitis C.

OBJECTIVE The natural course of chronic hepatitis C varies widely. To improve the profiling of patients at risk of developing advanced liver disease, we assessed the relative contribution of factors for liver fibrosis progression in hepatitis C. DESIGN We analysed 1461 patients with chronic hepatitis C with an estimated date of infection and at least one liver biopsy. Risk factors for accelerated fibrosis progression rate (FPR), defined as ≥0.13 Metavir fibrosis units per year, were identified by logistic regression. Examined factors included age at infection, sex, route of infection, HCV genotype, body mass index (BMI), significant alcohol drinking (≥20 g/day for ≥5 years), HIV coinfection and diabetes. In a subgroup of 575 patients, we assessed the impact of single nucleotide polymorphisms previously associated with fibrosis progression in genome-wide association studies. Results were expressed as attributable fraction (AF) of risk for accelerated FPR. RESULTS Age at infection (AF 28.7%), sex (AF 8.2%), route of infection (AF 16.5%) and HCV genotype (AF 7.9%) contributed to accelerated FPR in the Swiss Hepatitis C Cohort Study, whereas significant alcohol drinking, anti-HIV, diabetes and BMI did not. In genotyped patients, variants at rs9380516 (TULP1), rs738409 (PNPLA3), rs4374383 (MERTK) (AF 19.2%) and rs910049 (major histocompatibility complex region) significantly added to the risk of accelerated FPR. Results were replicated in three additional independent cohorts, and a meta-analysis confirmed the role of age at infection, sex, route of infection, HCV genotype, rs738409, rs4374383 and rs910049 in accelerating FPR. CONCLUSIONS Most factors accelerating liver fibrosis progression in chronic hepatitis C are unmodifiable.

A conserved role for p48 homologs in protecting dopaminergic neurons from oxidative stress.

Parkinson's disease (PD) is the most common neurodegenerative movement disorder characterized by the progressive loss of dopaminergic (DA) neurons. Both environmental and genetic factors are thought to contribute to the pathogenesis of PD. Although several genes linked to rare familial PD have been identified, endogenous risk factors for sporadic PD, which account for the majority of PD cases, remain largely unknown. Genome-wide association studies have identified many single nucleotide polymorphisms associated with sporadic PD in neurodevelopmental genes including the transcription factor p48/ptf1a. Here we investigate whether p48 plays a role in the survival of DA neurons in Drosophila melanogaster and Caenorhabditis elegans. We show that a Drosophila p48 homolog, 48-related-2 (Fer2), is expressed in and required for the development and survival of DA neurons in the protocerebral anterior medial (PAM) cluster. Loss of Fer2 expression in adulthood causes progressive PAM neuron degeneration in aging flies along with mitochondrial dysfunction and elevated reactive oxygen species (ROS) production, leading to the progressive locomotor deficits. The oxidative stress challenge upregulates Fer2 expression and exacerbates the PAM neuron degeneration in Fer2 loss-of-function mutants. hlh-13, the worm homolog of p48, is also expressed in DA neurons. Unlike the fly counterpart, hlh-13 loss-of-function does not impair development or survival of DA neurons under normal growth conditions. Yet, similar to Fer2, hlh-13 expression is upregulated upon an acute oxidative challenge and is required for the survival of DA neurons under oxidative stress in adult worms. Taken together, our results indicate that p48 homologs share a role in protecting DA neurons from oxidative stress and degeneration, and suggest that loss-of-function of p48 homologs in flies and worms provides novel tools to study gene-environmental interactions affecting DA neuron survival.

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.

Variants in the Mannose-binding Lectin Gene MBL2 do not Associate With Sepsis Susceptibility or Survival in a Large European Cohort.

BACKGROUND  Sepsis is an increasingly common condition, which continues to be associated with unacceptably high mortality. A large number of association studies have investigated susceptibility to, or mortality from, sepsis for variants in the functionally important immune-related gene MBL2. These studies have largely been underpowered and contradictory. METHODS  We genotyped and analyzed 4 important MBL2 single nucleotide polymorphisms (SNPs; rs5030737, rs1800450, rs1800451, and rs7096206) in 1839 European community-acquired pneumonia (CAP) and peritonitis sepsis cases, and 477 controls from the United Kingdom. We analyzed the following predefined subgroups and outcomes: 28-day and 6 month mortality from sepsis due to CAP or peritonitis combined, 28-day mortality from CAP sepsis, peritonitis sepsis, pneumococcal sepsis or sepsis in younger patients, and susceptibility to CAP sepsis or pneumococcal sepsis in the United Kingdom. RESULTS  There were no significant associations (all P-values were greater than .05 after correction for multiple testing) between MBL2 genotypes and any of our predefined analyses. CONCLUSIONS  In this large, well-defined cohort of immune competent adult patients, no associations between MBL2 genotype and sepsis susceptibility or outcome were identified.

Molecular characterization and SNP development for the porcine IL6 and IL10 genes

Different cytokines are secreted in response to specific microbial molecules referred to as pathogen associated molecular patterns (PAMPs). Interleukin 6 (IL6) and interleukin 10 (IL10), both secreted by macrophages and lymphocytes, play a central role in the immunological response. In this work we obtained the genomic structure and complete DNA sequence of the porcine IL6 and IL10 genes and identified polymorphisms in the genomic sequences of these genes on a panel of ten different pig breeds. Comparative intra- and interbreed sequence analysis revealed a total of eight polymorphisms in the porcine IL6 gene and 21 in the porcine IL10 gene, which include single nucleotide polymorphisms (SNPs) and insertion deletion polymorphisms (indels). Additionally, the chromosomal localization of the IL10 gene was determined by FISH and RH mapping.

Molecular characterization and chromosomal assignment of the bovine glycinamide ribonucleotide formyltransferase (GART) gene on cattle chromosome 1q12.1-q12.2

The mammalian glycinamide ribonucleotide formyltransferase (GART) genes encode a trifunctional polypeptide involved in the de novo purine biosynthesis. We isolated a bacterial artificial chromosome (BAC) clone containing the bovine GART gene and determined the complete DNA sequence of the BAC clone. Cloning and characterization of the bovine GART gene revealed that the bovine gene consists of 23 exons spanning approximately 27 kb. RT-PCR amplification of bovine GART in different organs showed the expression of two GART transcripts in cattle similar to human and mouse. The GART transcripts encode two proteins of 1010 and 433 amino acids, respectively. Eleven single nucleotide polymorphisms (SNPs) were detected in a mutation scan of 24 unrelated animals of three different cattle breeds, including one SNP that affects the amino acid sequence of GART. The chromosomal localization of the gene was determined by fluorescence in situ hybridization. Comparative genome analysis between cattle, human and mouse indicates that the chromosomal location of the bovine GART gene is in agreement with a previously published mapping report.

Imputation of sequence level genotypes in the Franches-Montagnes horse breed

BACKGROUND A cost-effective strategy to increase the density of available markers within a population is to sequence a small proportion of the population and impute whole-genome sequence data for the remaining population. Increased densities of typed markers are advantageous for genome-wide association studies (GWAS) and genomic predictions. METHODS We obtained genotypes for 54 602 SNPs (single nucleotide polymorphisms) in 1077 Franches-Montagnes (FM) horses and Illumina paired-end whole-genome sequencing data for 30 FM horses and 14 Warmblood horses. After variant calling, the sequence-derived SNP genotypes (~13 million SNPs) were used for genotype imputation with the software programs Beagle, Impute2 and FImpute. RESULTS The mean imputation accuracy of FM horses using Impute2 was 92.0%. Imputation accuracy using Beagle and FImpute was 74.3% and 77.2%, respectively. In addition, for Impute2 we determined the imputation accuracy of all individual horses in the validation population, which ranged from 85.7% to 99.8%. The subsequent inclusion of Warmblood sequence data further increased the correlation between true and imputed genotypes for most horses, especially for horses with a high level of admixture. The final imputation accuracy of the horses ranged from 91.2% to 99.5%. CONCLUSIONS Using Impute2, the imputation accuracy was higher than 91% for all horses in the validation population, which indicates that direct imputation of 50k SNP-chip data to sequence level genotypes is feasible in the FM population. The individual imputation accuracy depended mainly on the applied software and the level of admixture.

Gene flow and the maintenance of species boundaries

Hybrid zones are regions where individuals from genetically differentiated populations meet and mate, resulting in at least some offspring of mixed ancestry. Patterns of gene flow (introgression) in hybrid zones vary across the genome, allowing assessment of the role of individual genes or genome regions in reproductive isolation. Here, we document patterns of introgression between two recently diverged species of field crickets. We sampled at a very fine spatial scale and genotyped crickets for 110 highly differentiated single nucleotide polymorphisms (SNPs) identified through transcriptome scans. Using both genomic and geographic cline analysis, we document remarkably abrupt transitions (<100 m) in allele frequencies for 50 loci, despite high levels of gene flow at other loci. These are among the steepest clines documented for any hybridizing taxa. Furthermore, the cricket hybrid zone provides one of the clearest examples of the semi-permeability of species boundaries. Comparisons between data from the fine-scale transect and data (for the same set of markers) from sampling a much larger area in a different region of the cricket hybrid zone reveal consistent patterns of introgression for individual loci. The consistency in patterns of introgression between these two distant and distinct regions of the hybrid zone suggests that strong selection is acting to maintain abrupt discontinuities within the hybrid zone and that genomic regions with restricted introgression likely include genes that contribute to nonecological prezygotic barriers.

Association of lectin pathway proteins with intra-abdominal Candida infection in high-risk surgical intensive-care unit patients. A prospective cohort study within the fungal infection network of Switzerland

OBJECTIVES Human studies on the role of mannose-binding lectin (MBL) in patients with invasive candidiasis have yielded conflicting results. We investigated the influence of MBL and other lectin pathway proteins on Candida colonization and intra-abdominal candidiasis (IAC) in a cohort of high-risk patients. METHODS Prospective observational cohort study of 89 high-risk intensive-care unit (ICU) patients. Levels of lectin pathway proteins at study entry and six MBL2 single-nucleotide polymorphisms were analyzed by sandwich-type immunoassays and genotyping, respectively, and correlated with development of heavy Candida colonization (corrected colonization index (CCI) ≥0.4) and occurrence of IAC during a 4-week period. RESULTS Within 4 weeks after inclusion a CCI ≥0.4 and IAC was observed in 47% and 38% of patients respectively. Neither serum levels of MBL, ficolin-1, -2, -3, MASP-2 or collectin liver 1 nor MBL2 genotypes were associated with a CCI ≥0.4. Similarly, none of the analyzed proteins was found to be associated with IAC with the exception of lower MBL levels (HR 0.74, p = 0.02) at study entry. However, there was no association of MBL deficiency (<0.5 μg/ml), MBL2 haplo- or genotypes with IAC. CONCLUSION Lectin pathway protein levels and MBL2 genotype investigated in this study were not associated with heavy Candida colonization or IAC in a cohort of high-risk ICU patients.

Significant lethality following liver resection in A20 heterozygous knockout mice uncovers a key role for A20 in liver regeneration

Hepatic expression of A20, including in hepatocytes, increases in response to injury, inflammation and resection. This increase likely serves a hepatoprotective purpose. The characteristic unfettered liver inflammation and necrosis in A20 knockout mice established physiologic upregulation of A20 as integral to the anti-inflammatory and anti-apoptotic armamentarium of hepatocytes. However, the implication of physiologic upregulation of A20 in modulating hepatocytes' proliferative responses following liver resection remains controversial. To resolve the impact of A20 on hepatocyte proliferation and the liver's regenerative capacity, we examined whether decreased A20 expression, as in A20 heterozygous knockout mice, affects outcome following two-third partial hepatectomy. A20 heterozygous mice do not demonstrate a striking liver phenotype, indicating that their A20 expression levels are still sufficient to contain inflammation and cell death at baseline. However, usually benign partial hepatectomy provoked a staggering lethality (>40%) in these mice, uncovering an unsuspected phenotype. Heightened lethality in A20 heterozygous mice following partial hepatectomy resulted from impaired hepatocyte proliferation due to heightened levels of cyclin-dependent kinase inhibitor, p21, and deficient upregulation of cyclins D1, E and A, in the context of worsened liver steatosis. A20 heterozygous knockout minimally affected baseline liver transcriptome, mostly circadian rhythm genes. Nevertheless, this caused differential expression of >1000 genes post hepatectomy, hindering lipid metabolism, bile acid biosynthesis, insulin signaling and cell cycle, all critical cellular processes for liver regeneration. These results demonstrate that mere reduction of A20 levels causes worse outcome post hepatectomy than full knockout of bona fide liver pro-regenerative players such as IL-6, clearly ascertaining A20's primordial role in enabling liver regeneration. Clinical implications of these data are of utmost importance as they caution safety of extensive hepatectomy for donation or tumor in carriers of A20/TNFAIP3 single nucleotide polymorphisms alleles that decrease A20 expression or function, and prompt the development of A20-based liver pro-regenerative therapies.

Uric Acid and Cardiovascular Events: A Mendelian Randomization Study

Obesity and diets rich in uric acid-raising components appear to account for the increased prevalence of hyperuricemia in Westernized populations. Prevalence rates of hypertension, diabetes mellitus, CKD, and cardiovascular disease are also increasing. We used Mendelian randomization to examine whether uric acid is an independent and causal cardiovascular risk factor. Serum uric acid was measured in 3315 patients of the Ludwigshafen Risk and Cardiovascular Health Study. We calculated a weighted genetic risk score (GRS) for uric acid concentration based on eight uric acid-regulating single nucleotide polymorphisms. Causal odds ratios and causal hazard ratios (HRs) were calculated using a two-stage regression estimate with the GRS as the instrumental variable to examine associations with cardiometabolic phenotypes (cross-sectional) and mortality (prospectively) by logistic regression and Cox regression, respectively. Our GRS was not consistently associated with any biochemical marker except for uric acid, arguing against pleiotropy. Uric acid was associated with a range of prevalent diseases, including coronary artery disease. Uric acid and the GRS were both associated with cardiovascular death and sudden cardiac death. In a multivariate model adjusted for factors including medication, causal HRs corresponding to each 1-mg/dl increase in genetically predicted uric acid concentration were significant for cardiovascular death (HR, 1.77; 95% confidence interval, 1.12 to 2.81) and sudden cardiac death (HR, 2.41; 95% confidence interval, 1.16 to 5.00). These results suggest that high uric acid is causally related to adverse cardiovascular outcomes, especially sudden cardiac death.