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.

Standard genotyping overestimates transmission of Mycobacterium tuberculosis among immigrants in a low incidence country.

Immigrants from high tuberculosis (TB) incidence regions are a risk group for TB in low-incidence countries such as Switzerland. In a previous analysis of a nationwide collection of 520 Mycobacterium tuberculosis isolates from 2000-2008, we identified 35 clusters comprising 90 patients based on standard genotyping (24-loci MIRU-VNTR and spoligotyping). Here, we used whole genome sequencing (WGS) to revisit these transmission clusters. Genome-based transmission clusters were defined as isolate pairs separated by ≤12 single nucleotide polymorphisms (SNPs). WGS confirmed 17/35 (49%) MIRU-VNTR clusters; the other 18 clusters contained pairs separated by >12 SNPs. Most transmission clusters (3/4) of Swiss-born patients were confirmed by WGS, as opposed to 25% (4/16) of clusters involving only foreign-born patients. The overall clustering proportion using standard genotyping was 17% (90 patients, 95% confidence interval [CI]: 14-21%), but only 8% (43 patients, 95% CI: 6-11%) using WGS. The clustering proportion was 17% (67/401, 95% CI: 13-21%) using standard genotyping and 7% (26/401, 95% CI: 4-9%) using WGS among foreign-born patients, and 19% (23/119, 95% CI: 13-28%) and 14% (17/119, 95% CI: 9-22%), respectively, among Swiss-born patients. Using weighted logistic regression, we found weak evidence for an association between birth origin and transmission (aOR 2.2, 95% CI: 0.9-5.5, comparing Swiss-born patients to others). In conclusion, standard genotyping overestimated recent TB transmission in Switzerland when compared to WGS, particularly among immigrants from high TB incidence regions, where genetically closely related strains often predominate. We recommend the use of WGS to identify transmission clusters in low TB incidence settings.

Evaluation of candidate genes for susceptibility to spina bifida

Neural tube defects (NTDs) are the most common severely disabling birth defects in the United States, with a frequency of approximately 1–2 of every 1,000 births. This text includes the identification and evaluation of candidate susceptibility genes that confer risk for the development of neural tube defects (NTDs). The project focused on isolated meningomyelocele, also termed spina bifida (SB). ^ Spina bifida is a complex disease with multifactorial inheritance, therefore the subject population (consisting of North American Caucasians and Hispanics of Mexicali-American descent) was composed of 459 simplex SB families who were tested for genetic associations utilizing the transmission disequilibrium test (TDT), a nonparametric linkage technique. Three categories of candidate genes were studied, including (1) human equivalents of genes determined in mouse models to cause NTDs, (2) HOX and PAX genes, and (3) the MTHFR gene involved in the metabolic pathway of folate. ^ The C677T variant of the 5,10-methylenetetrahydrofolate reductase (MTHFR) gene was the first mutation in this gene to be implicated as a risk factor for NTDs. Our evaluation of the MTHFR gene provides evidence that maternal C677T homozygosity is a risk factor for upper level spina bifida defects in Hispanics [OR = 2.3, P = 0.02]. This observed risk factor is of great importance due to the high prevalence of this homozygous genotype in the Hispanic population. Additionally, maternal C677T/A1298C compound heterozygosity is a risk factor for upper level spina bifida defects in non-Hispanic whites [OR = 3.6, P = 0.03]. ^ For TDT analysis, our total population of 1128 subjects were genotyped for 54 markers from within and/or flanking the 20 candidate genes/gene regions of interest. Significant TDT findings were obtained for 3 of the 54 analyzed markers: d20s101 flanking the PAX1 gene (P = 0.019), d1s228 within the PAX7 gene (P = 0.011), and d2s110 within the PAX8 gene (P = 0.013). These results were followed-up by testing the genes directly for mutations utilizing single-strand conformational analysis (SSCA) and direct sequencing. Multiple variations were detected in each of these PAX genes; however, these variations were not passed from parent to child in phase with the positively transmitted alleles. Therefore, these variations do not contribute to the susceptibility of spina bifida, but rather are previously unreported single nucleotide polymorphisms. ^