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.

Contribution of Genetic Background and Clinical Risk Factors to Low-Trauma Fractures in Human Immunodeficiency Virus (HIV)-Positive Persons: The Swiss HIV Cohort Study.

Background.  The impact of human genetic background on low-trauma fracture (LTF) risk has not been evaluated in the context of human immunodeficiency virus (HIV) and clinical LTF risk factors. Methods.  In the general population, 6 common single-nucleotide polymorphisms (SNPs) associate with LTF through genome-wide association study. Using genome-wide SNP arrays and imputation, we genotyped these SNPs in HIV-positive, white Swiss HIV Cohort Study participants. We included 103 individuals with a first, physician-validated LTF and 206 controls matched on gender, whose duration of observation and whose antiretroviral therapy start dates were similar using incidence density sampling. Analyses of nongenetic LTF risk factors were based on 158 cases and 788 controls. Results.  A genetic risk score built from the 6 LTF-associated SNPs did not associate with LTF risk, in both models including and not including parental hip fracture history. The contribution of clinical LTF risk factors was limited in our dataset. Conclusions.  Genetic LTF markers with a modest effect size in the general population do not improve fracture prediction in persons with HIV, in whom clinical LTF risk factors are prevalent in both cases and controls.

Bone morphogenetic protein receptor 2 in patients with idiopathic portal hypertension

In idiopathic portal hypertension (IPH) typical vascular lesions are present in the branches of the portal vein or in the perisinusoidal area of the liver. Similar histological alterations have been reported in the pulmonary vasculature of patients with idiopathic pulmonary artery hypertension (IPAH). As IPAH is associated with mutations of the bone morphogenetic protein receptor 2 (BMPR2) gene, the aim of this study was to investigate whether this association might also be found in patients with IPH. Twenty-three samples belonging to 21 unrelated caucasian patients with IPH followed in the hepatic haemodynamic laboratory of the Hospital Clinic in Barcelona were included in the study. All patients were studied for the entire open reading frame and splice site of the BMPR2 gene by direct sequencing and multiple ligation probe amplification (MLPA) in order to detect large deletions/duplications. None of the 23 patients had pulmonary artery hypertension. Four patients presented one single nucleotide polymorphism (SNP) in intron 5, four patients had a SNP in exon 12 and a SNP in exon 1 was found in two cases. Two patients had both intron 5 and exon 12 polymorphisms. All SNPs were previously described. Except for these three SNPs, neither mutations nor rearrangements have been identified in the BMPR2 gene in this population. We did not detect mutations or rearrangements in the coding region of the BMPR2 gene in our patients with IPH. These findings suggest that, in contrast to IPAH, mutations in BMPR2 are not involved in the pathogenesis of IPH.

Genetic variations in bile acid homeostasis are not overrepresented in alcoholic cirrhosis compared to patients with heavy alcohol abuse and absent liver disease

Increased serum bile salt levels have been associated to a single-nucleotide polymorphism in the bile salt export pump (BSEP; ABCB11) in several acquired cholestatic liver diseases but there is little evidence in alcoholic liver disease (ALD). Furthermore, a crosstalk between vitamin D and bile acid synthesis has recently been discovered. Whether this crosstalk has an influence on the course of ALD is unclear to date. Our aim was to analyse the role of genetic polymorphisms in BSEP and the vitamin D receptor gene (NR1I1) on the emergence of cirrhosis in patients with ALD. Therefore, 511 alcoholic patients (131 with cirrhosis and 380 without cirrhosis) underwent ABCB11 genotyping (rs2287622). Of these, 321 (131 with cirrhosis and 190 without cirrhosis) were also tested for NR1I1 polymorphisms (bat-haplotype: BsmI rs1544410, ApaI rs7975232 and TaqI rs731236). Frequencies of ABCB11 and NR1I1 genotypes and haplotypes were compared between alcoholic patients with and without cirrhosis and correlated to serum bile salt, bilirubin and aspartate aminotransferase levels in those with cirrhosis. Frequencies of ABCB11 and NR1I1 genotypes and haplotypes did not differ between the two subgroups and no significant association between genotypes/haplotypes and liver function tests could be determined for neither polymorphism. We conclude that ABCB11 and NR1I1 polymorphisms are obviously not associated with development of cirrhosis in patients with ALD.

Age-modulated association between prefrontal NAA and the BDNF gene

Brain-derived neurotrophic factor (BDNF) has been implicated in the pathophysiology of psychiatric and neurological disorders and in the mechanisms of antidepressant pharmacotherapy. Psychiatric and neurological conditions have also been associated with reduced brain levels of N-acetyl-aspartate (NAA), which has been used as a putative marker of neural integrity. However, few studies have explored the relationship between BDNF polymorphisms and NAA levels directly. Here, we present data from a single-voxel proton magnetic resonance spectroscopy study of 64 individuals and explore the relationship between BDNF polymorphisms and prefrontal NAA level. Our results indicate an association between a single nucleotide polymorphism (SNP) within BDNF, known as rs1519480, and reduced NAA level (p = 0.023). NAA levels were further predicted by age and Asian ancestry. There was a significant rs1519480 × age interaction on NAA level (p = 0.031). Specifically, the effect of rs1519480 on NAA level became significant at age ⩾34.17 yr. NAA level decreased with advancing age for genotype TT (p = 0.001) but not for genotype CT (p = 0.82) or CC (p = 0.34). Additional in silico analysis of 142 post-mortem brain samples revealed an association between the same SNP and reduced BDNF mRNA expression in the prefrontal cortex. The rs1519480 SNP influences BDNF mRNA expression and has an impact on prefrontal NAA level over time. This genetic mechanism may contribute to inter-individual variation in cognitive performance seen during normal ageing, as well as contributing to the risk for developing psychiatric and neurological conditions.

A noncoding melanophilin gene (MLPH) SNP at the splice donor of exon 1 represents a candidate causal mutation for coat color dilution in dogs

Coat color dilution in several breeds of dog is characterized by a specific pigmentation phenotype and sometimes accompanied by hair loss and recurrent skin inflammation, the so-called color dilution alopecia or black hair follicular dysplasia. Coat color dilution (d) is inherited as a Mendelian autosomal recessive trait. In a previous study, MLPH polymorphisms showed perfect cosegregation with the dilute phenotype within breeds. However, different dilute haplotypes were found in different breeds, and no single polymorphism was identified in the coding sequence that was likely to be causative for the dilute phenotype. We resequenced the 5'-region of the canine MLPH gene and identified a strong candidate single nucleotide polymorphism within the nontranslated exon 1, which showed perfect association to the dilute phenotype in 65 dilute dogs from 7 different breeds. The A/G polymorphism is located at the last nucleotide of exon 1 and the mutant A-allele is predicted to reduce splicing efficiency 8-fold. An MLPH mRNA expression study using quantitative reverse transcriptase-polymerase chain reaction confirmed that dd animals had only about approximately 25% of the MLPH transcript compared with DD animals. These results provide preliminary evidence that the reported regulatory MLPH mutation might represent a causal mutation for coat color dilution in dogs.

A MIF haplotype is associated with the outcome of patients with severe sepsis: a case control study

BACKGROUND: Macrophage migration inhibitory factor (MIF) plays an important regulatory role in sepsis. In the promoter region a C/G single nucleotide polymorphism (SNP) at position -173 (rs755622) and a CATT5-8 microsatellite at position -794 are related to modified promoter activity. The purpose of the study was to analyze their association with the incidence and outcome of severe sepsis. METHODS: Genotype distributions and allele frequencies in 169 patients with severe sepsis, 94 healthy blood donors and 183 postoperative patients without signs of infection or inflammation were analyzed by real time PCR and Sequence analysis. All included individuals were Caucasians. RESULTS: Genotype distribution and allele frequencies of severe sepsis patients were comparable to both control groups. However, the genotype and allele frequencies of both polymorphisms were associated significantly with the outcome of severe sepsis. The highest risk of dying from severe sepsis was detectable in patients carrying a haplotype with the alleles -173 C and CATT7 (p = 0.0005, fisher exact test, RR = 1,806, CI: 1.337 to 2.439). CONCLUSION: The haplotype with the combination of the -173 C allele and the -794 CATT7 allele may not serve as a marker for susceptibility to sepsis, but may help identify septic patients at risk of dying.

SNPs for parentage testing and traceability in globally diverse breeds of sheep

DNA-based parentage determination accelerates genetic improvement in sheep by increasing pedigree accuracy. Single nucleotide polymorphism (SNP) markers can be used for determining parentage and to provide unique molecular identifiers for tracing sheep products to their source. However, the utility of a particular "parentage SNP" varies by breed depending on its minor allele frequency (MAF) and its sequence context. Our aims were to identify parentage SNPs with exceptional qualities for use in globally diverse breeds and to develop a subset for use in North American sheep. Starting with genotypes from 2,915 sheep and 74 breed groups provided by the International Sheep Genomics Consortium (ISGC), we analyzed 47,693 autosomal SNPs by multiple criteria and selected 163 with desirable properties for parentage testing. On average, each of the 163 SNPs was highly informative (MAF≥0.3) in 48±5 breed groups. Nearby polymorphisms that could otherwise confound genetic testing were identified by whole genome and Sanger sequencing of 166 sheep from 54 breed groups. A genetic test with 109 of the 163 parentage SNPs was developed for matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry. The scoring rates and accuracies for these 109 SNPs were greater than 99% in a panel of North American sheep. In a blinded set of 96 families (sire, dam, and non-identical twin lambs), each parent of every lamb was identified without using the other parent's genotype. In 74 ISGC breed groups, the median estimates for probability of a coincidental match between two animals (PI), and the fraction of potential adults excluded from parentage (PE) were 1.1×10(-39) and 0.999987, respectively, for the 109 SNPs combined. The availability of a well-characterized set of 163 parentage SNPs facilitates the development of high-throughput genetic technologies for implementing accurate and economical parentage testing and traceability in many of the world's sheep breeds.