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.