A genotype calling algorithm for the Illumina BeadArray platform

Genetic effects on the skin methylome in healthy older twins

Whole-skin DNA methylation variation has been implicated in several diseases, including melanoma, but its genetic basis has not yet been fully characterized. Using bulk skin tissue samples from 414 healthy female UK twins, we performed twin-based heritability and methylation quantitative trait loci (meQTL) analyses for >400,000 DNA methylation sites. We find that the human skin DNA methylome is on average less heritable than previously estimated in blood and other tissues (mean heritability: 10.02%). meQTL analysis identified local genetic effects influencing DNA methylation at 18.8% (76,442) of tested CpG sites, as well as 1,775 CpG sites associated with at least one distal genetic variant. As a functional follow-up, we performed skin expression QTL (eQTL) analyses in a partially overlapping sample of 604 female twins. Colocalization analysis identified over 3,500 shared genetic effects affecting thousands of CpG sites (10,067) and genes (4,475). Mediation analysis of putative colocalized gene-CpG pairs identified 114 genes with evidence for eQTL effects being mediated by DNA methylation in skin, including in genes implicating skin disease such as ALOX12 and CSPG4. We further explored the relevance of skin meQTLs to skin disease and found that skin meQTLs and CpGs under genetic influence were enriched for multiple skin-related genome-wide and epigenome-wide association signals, including for melanoma and psoriasis. Our findings give insights into the regulatory landscape of epigenomic variation in skin.

GenoTriplo: A SNP genotype calling method for triploids

Triploidy is very useful in both aquaculture and some cultivated plants as the induced sterility helps to enhance growth and product quality, as well as acting as a barrier against the contamination of wild populations by escapees. To use genetic information from triploids for academic or breeding purposes, an efficient and robust method to genotype triploids is needed. We developed such a method for genotype calling from SNP arrays, and we implemented it in the R package named GenoTriplo. Our method requires no prior information on cluster positions and remains unaffected by shifted luminescence signals. The method relies on starting the clustering algorithm with an initial higher number of groups than expected from the ploidy level of the samples, followed by merging groups that are too close to each other to be considered as distinct genotypes. Accurate classification of SNPs is achieved through multiple thresholds of quality controls. We compared the performance of GenoTriplo with that of fitPoly, the only published method for triploid SNP genotyping with a free software access. This was assessed by comparing the genotypes generated by both methods for a dataset of 1232 triploid rainbow trout genotyped for 38,033 SNPs. The two methods were consistent for 89% of the genotypes, but for 26% of the SNPs, they exhibited a discrepancy in the number of different genotypes identified. For these SNPs, GenoTriplo had >95% concordance with fitPoly when fitPoly genotyped better. On the contrary, when GenoTriplo genotyped better, fitPoly had less than 50% concordance with GenoTriplo. GenoTriplo was more robust with less genotyping errors. It is also efficient at identifying low-frequency genotypes in the sample set. Finally, we assessed parentage assignment based on GenoTriplo genotyping and observed significant differences in mismatch rates between the best and second-best couples, indicating high confidence in the results. GenoTriplo could also be used to genotype diploids as well as individuals with higher ploidy level by adjusting a few input parameters.

Does social intolerance vary according to cognitive styles, genetic cognitive capacity, or education?

BACKGROUND

Low education, low cognitive abilities, and certain cognitive styles are suggested to predispose to social intolerance and prejudices. Evidence is, however, restricted by comparatively small samples, neglect of confounding variables and genetic factors, and a narrow focus on a single sort of prejudice. We investigated the relationships of education, polygenic cognitive potential, cognitive performance, and cognitive styles with social intolerance in adulthood over a 15-year follow-up.

METHODS

We used data from the prospective population-based Young Finns Study (n = 960-1679). Social intolerance was evaluated with the Social Intolerance Scale in 1997, 2001, and 2011; cognitive performance with the Cambridge Neuropsychological Test Automated Battery in 2011; cognitive styles in 1997; and socioeconomic factors in 1980 (childhood) and 2011 (adulthood); and polygenic cognitive potential was calculated based on genome-wide association studies.

RESULTS

We found that nonrational thinking, polygenic cognitive potential, cognitive performance, or socioeconomic factors were not related to social intolerance. Regarding cognitive styles, low flexibility (B = -0.759, p < .001), high perseverance (B = 1.245, p < .001), and low persistence (B = -0.329, p < .001) predicted higher social intolerance consistently in the analyses.

DISCUSSION

When developing prejudice-reduction interventions, it should be considered that educational level or cognitive performance may not be crucial for development of social intolerance. Adopting certain cognitive styles may play more important roles in development of social intolerance.

Life-Course Associations between Blood Pressure-Related Polygenic Risk Scores and Hypertension in the Bogalusa Heart Study

Genetic information may help to identify individuals at increased risk for hypertension in early life, prior to the manifestation of elevated blood pressure (BP) values. We examined 369 Black and 832 White Bogalusa Heart Study (BHS) participants recruited in childhood and followed for approximately 37 years. The multi-ancestry genome-wide polygenic risk scores (PRSs) for systolic BP (SBP), diastolic BP (DBP), and hypertension were tested for an association with incident hypertension and stage 2 hypertension using Cox proportional hazards models. Race-stratified analyses were adjusted for baseline age, age2, sex, body mass index, genetic principal components, and BP. In Black participants, each standard deviation increase in SBP and DBP PRS conferred a 38% (p = 0.009) and 22% (p = 0.02) increased risk of hypertension and a 74% (p < 0.001) and 50% (p < 0.001) increased risk of stage 2 hypertension, respectively, while no association was observed with the hypertension PRSs. In Whites, each standard deviation increase in SBP, DBP, and hypertension PRS conferred a 24% (p < 0.05), 29% (p = 0.01), and 25% (p < 0.001) increased risk of hypertension, and a 27% (p = 0.08), 29% (0.01), and 42% (p < 0.001) increased risk of stage 2 hypertension, respectively. The addition of BP PRSs to the covariable-only models generally improved the C-statistics (p < 0.05). Multi-ancestry BP PRSs demonstrate the utility of genomic information in the early life prediction of hypertension.

Does better education mitigate risky health behavior? A mendelian randomization study

Education and risky health behaviors are strongly negatively correlated. Education may affect health behaviors by enabling healthier choices through higher disposable income, increasing information about the harmful effects of risky health behaviors, or altering time preferences. Alternatively, the observed negative correlation may stem from reverse causality or unobserved confounders. Based on the data from the Cardiovascular Risk in Young Finns Study linked to register-based information on educational attainment and family background, this paper identifies the causal effect of education on risky health behaviors. To examine causal effects, we used a genetic score as an instrument for years of education. We found that individuals with higher education allocated more attention to healthy habits. In terms of health behaviors, highly educated people were less likely to smoke. Some model specifications also indicated that the highly educated consumed more fruit and vegetables, but the results were imprecise in this regard. No causal effect was found between education and abusive drinking. In brief, inference based on genetic instruments showed that higher education leads to better choices in some but not all dimensions of health behaviors.

The relationship between temperament, polygenic score for intelligence and cognition: A population-based study of middle-aged adults

We investigated whether temperament modifies an association between polygenic intelligence potential and cognitive test performance in midlife. The participants (n = 1647, born between 1962 and 1977) were derived from the Young Finns Study. Temperament was assessed with Temperament and Character Inventory over a 15-year follow-up (1997, 2001, 2007, 2012). Polygenic intelligence potential was assessed with a polygenic score for intelligence. Cognitive performance (visual memory, reaction time, sustained attention, spatial working memory) was assessed with CANTAB in midlife. The PGSI was significantly associated with the overall cognitive performance and performance in visual memory, sustained attention and working memory tests but not reaction time test. Temperament did not correlate with polygenic score for intelligence and did not modify an association between the polygenic score and cognitive performance, either. High persistence was associated with higher visual memory (B = 0.092; FDR-adj. p = 0.007) and low harm avoidance with higher overall cognitive performance, specifically better reaction time (B = -0.102; FDR-adj; p = 0.007). The subscales of harm avoidance had different associations with cognitive performance: higher "anticipatory worry," higher "fatigability," and lower "shyness with strangers" were associated with lower cognitive performance, while the role of "fear of uncertainty" was subtest-related. In conclusion, temperament does not help or hinder one from realizing their genetic potential for intelligence. The overall modest relationships between temperament and cognitive performance advise caution if utilizing temperament-related information e.g. in working-life recruitments. Cognitive abilities may be influenced by temperament variables, such as the drive for achievement and anxiety about test performance, but they involve distinct systems of learning and memory.

Birth weight, adult weight, and cardiovascular biomarkers: Evidence from the Cardiovascular Young Finns Study

This study quantifies the causal effect of birth weight on cardiovascular biomarkers in adulthood using the Cardiovascular Risk in Young Finns Study (YFS). We apply a multivariable Mendelian randomization (MVMR) method that provides a novel approach to improve inference in causal analysis based on a mediation framework. The results show that birth weight is linked to triglyceride levels (β = -0.294; 95% CI [-0.591, 0.003]) but not to low-density lipoprotein (LDL) cholesterol levels (β = 0.007; 95% CI [-0.168, 0.183]). The total effect of birth weight on triglyceride levels is partly offset by a mediation pathway linking birth weight to adult BMI (β = 0.111; 95% CI [-0.013, 0.234]). The negative total effect is consistent with the fetal programming hypothesis. The positive indirect effect via adult BMI highlights the persistence of body weight throughout a person's life and the adverse effects of high BMI on health. The results are consistent with previous findings that both low birth weight and weight gain increase health risks in adulthood.

Genome-Wide association between EYA1 and Aspirin-induced peptic ulceration

Background

Low-dose aspirin can cause gastric and duodenal ulceration, hereafter called peptic ulcer disease (PUD). Predisposition is thought to be related to clinical and genetic factors; our aim was to identify genetic risk factors associated with aspirin-induced PUD.

Methods

Patients (n=1478) were recruited from 15 UK hospitals. Cases (n=505) were defined as patients with endoscopically confirmed PUD within 2 weeks of using aspirin and non-aspirin Non-Steroidal Anti-Inflammatory Drugs (NSAIDs). They were compared to two control groups: patients with endoscopically confirmed PUD without any history of NSAID use within 3 months of diagnosis (n=495), and patients with no PUD on endoscopy (n=478). A genome-wide association study (GWAS) of aspirin-induced cases (n=247) was compared to 476 controls. The results were validated by replication in another 84 cases and 162 controls.

Findings

The GWAS identified one variant, rs12678747 (p=1·65×10⁻⁷) located in the last intron of EYA1 on chromosome 8. The association was replicated in another sample of 84 PUD patients receiving aspirin (p=0·002). Meta-analysis of discovery and replication cohort data for rs12678747, yielded a genome-wide significant association (p=3·12×10⁻¹¹; OR=2·03; 95% CI 1·65-2·50). Expression of EYA1 was lower at the gastric ulcer edge when compared with the antrum.

Interpretation

Genetic variation in an intron of the EYA1 gene increases the risk of endoscopically confirmed aspirin-induced PUD. Reduced EYA1 expression in the upper gastrointestinal epithelium may modulate risk, but the functional basis of this association will need mechanistic evaluation.

Funding

Department of Health Chair in Pharmacogenetics, MRC Centre for Drug Safety Science and the Barts Cardiovascular NIHR Biomedical Research Centre, British Heart Foundation (BHF)

Birth weight and adult income: An examination of mediation through adult height and body mass

This paper examines the causal links between early human endowments and socioeconomic outcomes in adulthood. We use a genotyped longitudinal survey (Cardiovascular Risk in Young Finns Study) that is linked to the administrative registers of Statistics Finland. We focus on the effect of birth weight on income via two anthropometric mediators: body mass index (BMI) and height in adulthood. We find that (i) the genetic instruments for birth weight, adult height, and adult BMI are statistically powerful; (ii) there is a robust total effect of birth weight on income for men but not for women; (iii) the total effect of birth weight on income for men is partly mediated via height but not via BMI; and (iv) the share of the total effect mediated via height is substantial, of approximately 56%.

Donor genetic variants in interleukin-6 and interleukin-6 receptor associate with biopsy-proven rejection following kidney transplantation

Rejection after kidney transplantation remains an important cause of allograft failure that markedly impacts morbidity. Cytokines are a major player in rejection, and we, therefore, explored the impact of interleukin-6 (IL6) and IL-6 receptor (IL6R) gene polymorphisms on the occurrence of rejection after renal transplantation. We performed an observational cohort study analyzing both donor and recipient DNA in 1271 renal transplant‐pairs from the University Medical Center Groningen in The Netherlands and associated single nucleotide polymorphisms (SNPs) with biopsy-proven rejection after kidney transplantation. The C-allele of the IL6R SNP (Asp358Ala; rs2228145 A > C, formerly rs8192284) in donor kidneys conferred a reduced risk of rejection following renal transplantation (HR 0.78 per C-allele; 95%-CI 0.67–0.90; P = 0.001). On the other hand, the C-allele of the IL6 SNP (at position-174 in the promoter; rs1800795 G > C) in donor kidneys was associated with an increased risk of rejection for male organ donors (HR per C-allele 1.31; 95%-CI 1.08–1.58; P = 0.0006), but not female organ donors (P = 0.33). In contrast, neither the IL6 nor IL6R SNP in the recipient showed an association with renal transplant rejection. In conclusion, donor IL6 and IL6R genotypes but not recipient genotypes represent an independent prognostic marker for biopsy-proven renal allograft rejection.

Cardiovascular Risk Factor Trajectories Since Childhood and Cognitive Performance in Midlife: The Cardiovascular Risk in Young Finns Study

BACKGROUND

Cardiovascular risk factors, such as high blood pressure, adverse serum lipids, and elevated body mass index in midlife, may harm cognitive performance. It is important to note that longitudinal accumulation of cardiovascular risk factors since childhood may be associated with cognitive performance already since childhood, but the previous evidence is scarce. We studied the associations of cardiovascular risk factors from childhood to midlife, their accumulation, and midlife cognitive performance.

METHODS

From 1980, a population-based cohort of 3596 children (3-18 years of age) have been repeatedly followed up for 31 years. Blood pressure, serum lipids, and body mass index were assessed in all follow-ups. Cardiovascular risk factor trajectories from childhood to midlife were identified using latent class growth mixture modeling. Cognitive testing was performed in 2026 participants 34 to 49 years of age using a computerized test. The associations of the cardiovascular risk factor trajectories and cognitive performance were studied for individual cardiovascular risk factors and cardiovascular risk factor accumulation.

RESULTS

Consistently high systolic blood pressure (β=-0.262 SD [95% CI, -0.520 to -0.005]) and serum total cholesterol (β=-0.214 SD [95% CI, -0.365 to -0.064]) were associated with worse midlife episodic memory and associative learning compared with consistently low values. Obesity since childhood was associated with worse visual processing and sustained attention (β=-0.407 SD [95% CI, -0.708 to -0.105]) compared with normal weight. An inverse association was observed for the cardiovascular risk factor accumulation with episodic memory and associative learning (P for trend=0.008; 3 cardiovascular risk factors: β=-0.390 SD [95% CI, -0.691 to -0.088]), with visual processing and sustained attention (P for trend<0.0001; 3 cardiovascular risk factors: β=-0.443 SD [95% CI, -0.730 to -0.157]), and with reaction and movement time (P for trend=0.048; 2 cardiovascular risk factors: β=-0.164 SD [95% CI, -0.318 to -0.010]).

CONCLUSIONS

Longitudinal elevated systolic blood pressure, high serum total cholesterol, and obesity from childhood to midlife were inversely associated with midlife cognitive performance. It is important to note that the higher the number of cardiovascular risk factors, the worse was the observed cognitive performance. Therefore, launching preventive strategies against cardiovascular risk factors beginning from childhood might benefit primordial promotion of cognitive health in adulthood.

Vitamin B12 and transcobalamin in bovine milk: Genetic variation and genome-wide association with loci along the genome

In human nutrition, bovine milk is an essential source of bioavailable vitamin B₁₂ and B₁₂-binding proteins, including transcobalamin. In this study, we estimated genetic parameters for milk content of vitamin B₁₂ and transcobalamin using milk samples from 341 and 663 Danish Holstein cows, respectively. Additionally, we conducted whole-genome association analysis to identify SNP and genes associated with vitamin B₁₂ and transcobalamin. Our results indicated moderate to high heritability for vitamin B₁₂ (0.37 ± 0.18) and transcobalamin (0.61 ± 0.13) content in the Danish Holstein. With a significance threshold of -log₁₀ P-value > 5.87, significant associations were detected between SNP in Bos taurus autosome (BTA)17 and the log-transformed transcobalamin content of milk; no significant association was detected for vitamin B₁₂. The significant region in BTA17 was imputed to full sequence for further fine mapping, and the SNP with the most significant associations to transcobalamin were assigned to the transcobalamin 2 (TCN2) gene.

The Role of Genetic Ancestry as a Risk Factor for Primary Open-angle Glaucoma in African Americans

Purpose

POAG is the leading cause of irreversible blindness in African Americans. In this study, we quantitatively assess the association of autosomal ancestry with POAG risk in a large cohort of self-identified African Americans.

Methods

Subjects recruited to the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study were classified as glaucoma cases or controls by fellowship-trained glaucoma specialists. POAAGG subjects were genotyped using the MEGA Ex array (discovery cohort, n = 3830; replication cohort, n = 2135). Population structure was interrogated using principal component analysis in the context of the 1000 Genomes Project superpopulations.

Results

The majority of POAAGG samples lie on an axis between African and European superpopulations, with great variation in admixture. Cases had a significantly lower mean value of the ancestral component q0 than controls for both cohorts (P = 6.14-4; P = 3-6), consistent with higher degree of African ancestry. Among POAG cases, higher African ancestry was also associated with thinner central corneal thickness (P = 2-4). Admixture mapping showed that local genetic ancestry was not a significant risk factor for POAG. A polygenic risk score, comprised of 23 glaucoma-associated single nucleotide polymorphisms from the NHGRI-EBI genome-wide association study catalog, was significant in both cohorts (P < 0.001), suggesting that both known POAG single nucleotide polymorphisms and an omnigenic ancestry effect influence POAG risk.

Conclusions

In sum, the POAAGG study population is very admixed, with a higher degree of African ancestry associated with an increased POAG risk. Further analyses should consider social and environmental factors as possible confounding factors for disease predisposition.

Diabetes and Genetics: A Relationship Between Genetic Risk Alleles, Clinical Phenotypes and Therapeutic Approaches

Unveiling human genome through successful completion of Human Genome Project and International HapMap Projects with the advent of state of art technologies has shed light on diseases associated genetic determinants. Identification of mutational landscapes such as copy number variation, single nucleotide polymorphisms or variants in different genes and loci have revealed not only genetic risk factors responsible for diseases but also region(s) playing protective roles. Diabetes is a global health concern with two major types - type 1 diabetes (T1D) and type 2 diabetes (T2D). Great progress in understanding the underlying genetic predisposition to T1D and T2D have been made by candidate gene studies, genetic linkage studies, genome wide association studies with substantial number of samples. Genetic information has importance in predicting clinical outcomes. In this review, we focus on recent advancement regarding candidate gene(s) associated with these two traits along with their clinical parameters as well as therapeutic approaches perceived. Understanding genetic architecture of these disease traits relating clinical phenotypes would certainly facilitate population stratification in diagnosing and treating T1D/T2D considering the doses and toxicity of specific drugs.

Childhood Exposure to Parental Smoking and Midlife Cognitive Function

We studied whether exposure to parental smoking in childhood/adolescence is associated with midlife cognitive function, leveraging data from the Cardiovascular Risk in Young Finns Study. A population-based cohort of 3,596 children/adolescents aged 3-18 years was followed between 1980 and 2011. In 2011, cognitive testing was performed on 2,026 participants aged 34-49 years using computerized testing. Measures of secondhand smoke exposure in childhood/adolescence consisted of parental self-reports of smoking and participants' serum cotinine levels. Participants were classified into 3 exposure groups: 1) no exposure (nonsmoking parents, cotinine <1.0 ng/mL); 2) hygienic parental smoking (1-2 smoking parents, cotinine <1.0 ng/mL); and 3) nonhygienic parental smoking (1-2 smoking parents, cotinine ≥1.0 ng/mL). Analyses adjusted for sex, age, family socioeconomic status, polygenic risk score for cognitive function, adolescent/adult smoking, blood pressure, and serum total cholesterol level. Compared with the nonexposed, participants exposed to nonhygienic parental smoking were at higher risk of poor (lowest quartile) midlife episodic memory and associative learning (relative risk (RR) = 1.38, 95% confidence interval (CI): 1.08, 1.75), and a weak association was found for short-term and spatial working memory (RR = 1.25, 95% CI: 0.98, 1.58). Associations for those exposed to hygienic parental smoking were nonsignificant (episodic memory and associative learning: RR = 1.19, 95% CI: 0.92, 1.54; short-term and spatial working memory: RR = 1.10, 95% CI: 0.85, 1.34). We conclude that avoiding childhood/adolescence secondhand smoke exposure promotes adulthood cognitive function.

Accurate calling of homeoallelic genotypes of iSelect markers using inbred structured populations

MOTIVATION

Genetic map construction is a foundational step in the analysis of structured experimental populations. For markers that hybridize to several genetically similar locations, or where several alleles are present (such as in multiparental populations), current methods often discard the marker or incorrectly call the genotypes. These errors result in information loss, or incorrect genotypes that can corrupt map construction.

RESULTS

We present a new approach for simultaneously performing genetic map construction and marker calling. Our new approach allows the calling of a larger number of markers, a larger number of unique alleles per marker and the correct use of markers which hybridize to multiple genetically similar locations. We demonstrate our new approach using simulations, a biparental wheat population and an eight-parent population of spring bread wheat. Applying our method to the eight-parent population increased the number of mapped markers by 71%. We show that the new genetic map allows the investigation of synteny in ways that were not previously possible in that dataset.

AVAILABILITY AND IMPLEMENTATION

The method described in this article has been incorporated into R package mpMap2. It is available from CRAN and also from https://github.com/rohan-shah/mpMap2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genetic determinants of the humoral immune response in MS

OBJECTIVE

In this observational study, we investigated the impact of genetic factors at the immunoglobulin heavy chain constant locus on chromosome 14 and the major histocompatibility complex region on intrathecal immunoglobulin G, A, and M levels as well as on B cells and plasmablasts in the CSF and blood of patients with multiple sclerosis (MS).

METHODS

Using regression analyses, we tested genetic variants on chromosome 14 and imputed human leukocyte antigen (HLA) alleles for associations with intrathecal immunoglobulins in 1,279 patients with MS or clinically isolated syndrome and with blood and CSF B cells and plasmablasts in 301 and 348 patients, respectively.

RESULTS

The minor alleles of variants on chromosome 14 were associated with higher intrathecal immunoglobulin G levels (β = 0.58 [0.47 to 0.68], lowest adjusted p = 2.32 × 10^-23), and lower intrathecal immunoglobulin M (β = -0.56 [-0.67 to -0.46], p = 2.06 × 10^-24) and A (β = -0.42 [-0.54 to -0.31], p = 7.48 × 10^-11) levels. Alleles from the HLA-B*07:02-DRB1*15:01-DQA1*01:02-DQB1*06:02 haplotype were associated with higher (lowest p = 2.14 × 10^-7) and HLA-B*44:02 with lower (β = -0.35 [-0.54 to -0.17], p = 1.38 × 10^-2) immunoglobulin G levels. Of interest, different HLA alleles were associated with lower intrathecal immunoglobulin M (HLA-C*02:02, β = -0.45 [-0.61 to -0.28], p = 1.01 × 10^-5) and higher immunoglobulin A levels (HLA-DQA1*01:03-DQB1*06:03-DRB1*13:01 haplotype, β = 0.40 [0.21 to 0.60], p = 4.46 × 10^-3). The impact of HLA alleles on intrathecal immunoglobulin G and M levels could mostly be explained by associations with CSF B cells and plasmablasts.

CONCLUSION

Although some HLA alleles seem to primarily drive the extent of humoral immune responses in the CNS by increasing CSF B cells and plasmablasts, genetic variants at the immunoglobulin heavy chain constant locus might regulate intrathecal immunoglobulins levels via different mechanisms.

SureTypeSC-a Random Forest and Gaussian mixture predictor of high confidence genotypes in single-cell data

MOTIVATION

Accurate genotyping of DNA from a single cell is required for applications such as de novo mutation detection, linkage analysis and lineage tracing. However, achieving high precision genotyping in the single-cell environment is challenging due to the errors caused by whole-genome amplification. Two factors make genotyping from single cells using single nucleotide polymorphism (SNP) arrays challenging. The lack of a comprehensive single-cell dataset with a reference genotype and the absence of genotyping tools specifically designed to detect noise from the whole-genome amplification step. Algorithms designed for bulk DNA genotyping cause significant data loss when used for single-cell applications.

RESULTS

In this study, we have created a resource of 28.7 million SNPs, typed at high confidence from whole-genome amplified DNA from single cells using the Illumina SNP bead array technology. The resource is generated from 104 single cells from two cell lines that are available from the Coriell repository. We used mother-father-proband (trio) information from multiple technical replicates of bulk DNA to establish a high quality reference genotype for the two cell lines on the SNP array. This enabled us to develop SureTypeSC-a two-stage machine learning algorithm that filters a substantial part of the noise, thereby retaining the majority of the high quality SNPs. SureTypeSC also provides a simple statistical output to show the confidence of a particular single-cell genotype using Bayesian statistics.

AVAILABILITY AND IMPLEMENTATION

The implementation of SureTypeSC in Python and sample data are available in the GitHub repository: https://github.com/puko818/SureTypeSC.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Longitudinal association of a body mass index (BMI) genetic risk score with growth and BMI changes across the life course: The Cardiovascular Risk in Young Finns Study

BACKGROUND

The role of genetic risk scores associated with adult body mass index (BMI) on BMI levels across the life course is unclear. We examined if a 97 single nucleotide polymorphism weighted genetic risk score (wGRS97) associated with age-related progression in BMI at different life stages and distinct developmental trajectories of BMI across the early life course.

METHODS

2188 Cardiovascular Risk in Young Finns Study participants born pre-1980 who had genotype data and objective measurements of height and weight collected up to 8 times from age 6 to 49 years. Associations were examined using Individual Growth Curve analysis, Latent Class Growth Mixture Modelling, and Poisson modified regression.

RESULTS

The wGRS97 associated with BMI from age 6 years with peak effect sizes observed at age 30 years (females: 1.14 kg/m²; males: 1.09 kg/m² higher BMI per standard deviation increase in wGRS97). The association between wGRS97 and BMI became stronger with age in childhood but slowed in adolescence, especially in females, and weakened at age 35-40 years. A higher wGRS97 associated with an increased BMI velocity in childhood and adulthood, but not with BMI change in adulthood. Compared with belonging to a 'normal stable' life-course trajectory group (normal BMI from childhood to adulthood), a one standard deviation higher wGRS97 associated with a 13-127% increased risk of belonging to a less favourable life-course BMI trajectory group.

CONCLUSIONS

Individuals with genetic susceptibility to higher adult BMI have higher levels and accelerated rates of increase in BMI in childhood/adolescence, and are at increased risk of having a less favourable life-course BMI trajectory.

Family-based genome-wide association study of leprosy in Vietnam

Leprosy is a chronic infectious disease of the skin and peripheral nerves with a strong genetic predisposition. Recent genome-wide approaches have identified numerous common variants associated with leprosy, almost all in the Chinese population. We conducted the first family-based genome-wide association study of leprosy in 622 affected offspring from Vietnam, followed by replication in an independent sample of 1181 leprosy cases and 668 controls of the same ethnic origin. The most significant results were observed within the HLA region, in which six SNPs displayed genome-wide significant associations, all of which were replicated in the independent case/control sample. We investigated the signal in the HLA region in more detail, by conducting a multivariate analysis on the case/control sample of 319 GWAS-suggestive HLA hits for which evidence for replication was obtained. We identified three independently associated SNPs, two located in the HLA class I region (rs1265048: OR = 0.69 [0.58-0.80], combined p-value = 5.53x10-11; and rs114598080: OR = 1.47 [1.46-1.48], combined p-value = 8.77x10-13), and one located in the HLA class II region (rs3187964 (OR = 1.67 [1.55-1.80], combined p-value = 8.35x10-16). We also validated two previously identified risk factors for leprosy: the missense variant rs3764147 in the LACC1 gene (OR = 1.52 [1.41-1.63], combined p-value = 5.06x10-14), and the intergenic variant rs6871626 located close to the IL12B gene (OR = 0.73 [0.61-0.84], combined p-value = 6.44x10-8). These results shed new light on the genetic control of leprosy, by dissecting the influence of HLA SNPs, and validating the independent role of two additional variants in a large Vietnamese sample.

Education leads to a more physically active lifestyle: Evidence based on Mendelian randomization

Physical inactivity is a major health risk worldwide. Observational studies suggest that higher education is positively related to physical activity, but it is not clear whether this relationship constitutes a causal effect. Using participants (N = 1651) drawn from the Cardiovascular Risk in Young Finns Study linked to nationwide administrative data from Statistics Finland, this study examined whether educational attainment, measured by years of education, is related to adulthood physical activity in terms of overall physical activity, weekly hours of intensive activity, total steps per day, and aerobic steps per day. We employed ordinary least squares (OLS) models and extended the analysis using an instrumental variables approach (Mendelian randomization, MR) with a genetic risk score as an instrument for years of education. Based on the MR results, it was found that years of education is positively related to physical activity. On average, one additional year of education leads to a 0.62-unit higher overall physical activity (P < .01), 0.26 more hours of weekly intensive activity (P < .05), 560 more steps per day (P < .10), and 390 more aerobic steps per day (P < .09). The findings indicate that education may be a factor leading to higher leisure-time physical activity and thus promoting global health.

Insights into malaria susceptibility using genome-wide data on 17,000 individuals from Africa, Asia and Oceania

The human genetic factors that affect resistance to infectious disease are poorly understood. Here we report a genome-wide association study in 17,000 severe malaria cases and population controls from 11 countries, informed by sequencing of family trios and by direct typing of candidate loci in an additional 15,000 samples. We identify five replicable associations with genome-wide levels of evidence including a newly implicated variant on chromosome 6. Jointly, these variants account for around one-tenth of the heritability of severe malaria, which we estimate as ~23% using genome-wide genotypes. We interrogate available functional data and discover an erythroid-specific transcription start site underlying the known association in ATP2B4, but are unable to identify a likely causal mechanism at the chromosome 6 locus.  Previously reported HLA associations do not replicate in these samples. This large dataset will provide a foundation for further research on the genetic determinants of malaria resistance in diverse populations.

Four genome-wide associated loci are currently known for malaria susceptibility. Here, the authors expand on earlier work by combining data from 11 malaria-endemic countries and additional population sequencing informing an African-enriched imputation reference panel, with findings including a previously unreported association on chromosome 6.

Multivariate Genome-wide Association Analysis of a Cytokine Network Reveals Variants with Widespread Immune, Haematological, and Cardiometabolic Pleiotropy

Cytokines are essential regulatory components of the immune system, and their aberrant levels have been linked to many disease states. Despite increasing evidence that cytokines operate in concert, many of the physiological interactions between cytokines, and the shared genetic architecture that underlies them, remain unknown. Here, we aimed to identify and characterize genetic variants with pleiotropic effects on cytokines. Using three population-based cohorts (n = 9,263), we performed multivariate genome-wide association studies (GWAS) for a correlation network of 11 circulating cytokines, then combined our results in meta-analysis. We identified a total of eight loci significantly associated with the cytokine network, of which two (PDGFRB and ABO) had not been detected previously. In addition, conditional analyses revealed a further four secondary signals at three known cytokine loci. Integration, through the use of Bayesian colocalization analysis, of publicly available GWAS summary statistics with the cytokine network associations revealed shared causal variants between the eight cytokine loci and other traits; in particular, cytokine network variants at the ABO, SERPINE2, and ZFPM2 loci showed pleiotropic effects on the production of immune-related proteins, on metabolic traits such as lipoprotein and lipid levels, on blood-cell-related traits such as platelet count, and on disease traits such as coronary artery disease and type 2 diabetes.

Abdominal adiposity and cardiometabolic risk factors in children and adolescents: a Mendelian randomization analysis

BACKGROUND

Mendelian randomization studies in adults suggest that abdominal adiposity is causally associated with increased risk of type 2 diabetes and coronary artery disease in adults, but its causal effect on cardiometabolic risk in children remains unclear.

OBJECTIVE

We aimed to study the causal relation of abdominal adiposity with cardiometabolic risk factors in children by applying Mendelian randomization.

METHODS

We constructed a genetic risk score (GRS) using variants previously associated with waist-to-hip ratio adjusted for BMI (WHRadjBMI) and examined its associations with cardiometabolic factors by linear regression and Mendelian randomization in a meta-analysis of 6 cohorts, including 9895 European children and adolescents aged 3-17 y.

RESULTS

WHRadjBMI GRS was associated with higher WHRadjBMI (β = 0.021 SD/allele; 95% CI: 0.016, 0.026 SD/allele; P = 3 × 10-15) and with unfavorable concentrations of blood lipids (higher LDL cholesterol: β = 0.006 SD/allele; 95% CI: 0.001, 0.011 SD/allele; P = 0.025; lower HDL cholesterol: β = -0.007 SD/allele; 95% CI: -0.012, -0.002 SD/allele; P = 0.009; higher triglycerides: β = 0.007 SD/allele; 95% CI: 0.002, 0.012 SD/allele; P = 0.006). No differences were detected between prepubertal and pubertal/postpubertal children. The WHRadjBMI GRS had a stronger association with fasting insulin in children and adolescents with overweight/obesity (β = 0.016 SD/allele; 95% CI: 0.001, 0.032 SD/allele; P = 0.037) than in those with normal weight (β = -0.002 SD/allele; 95% CI: -0.010, 0.006 SD/allele; P = 0.605) (P for difference = 0.034). In a 2-stage least-squares regression analysis, each genetically instrumented 1-SD increase in WHRadjBMI increased circulating triglycerides by 0.17 mmol/L (0.35 SD, P = 0.040), suggesting that the relation between abdominal adiposity and circulating triglycerides may be causal.

CONCLUSIONS

Abdominal adiposity may have a causal, unfavorable effect on plasma triglycerides and potentially other cardiometabolic risk factors starting in childhood. The results highlight the importance of early weight management through healthy dietary habits and physically active lifestyle among children with a tendency for abdominal adiposity.

Impact of the rumen microbiome on milk fatty acid composition of Holstein cattle

Background

Fatty acids (FA) in bovine milk derive through body mobilization, de novo synthesis or from the feed via the blood stream. To be able to digest feedstuff, the cow depends on its rumen microbiome. The relative abundance of the microbes has been shown to differ between cows. To date, there is little information on the impact of the microbiome on the formation of specific milk FA. Therefore, in this study, our aim was to investigate the impact of the rumen bacterial microbiome on milk FA composition. Furthermore, we evaluated the predictive value of the rumen microbiome and the host genetics on the composition of individual FA in milk.

Results

Our results show that the proportion of variance explained by the rumen bacteria composition (termed microbiability or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h_{B}^{2}$$\end{document}hB2) was generally smaller than that of the genetic component (heritability), and that rumen bacteria influenced most C15:0, C17:0, C18:2 n-6, C18:3 n-3 and CLA cis-9, trans-11 with estimated \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h_{B}^{2}$$\end{document}hB2 ranging from 0.26 to 0.42. For C6:0, C8:0, C10:0, C12:0, C16:0, C16:1 cis-9 and C18:1 cis-9, the variance explained by the rumen bacteria component was close to 0. In general, both the rumen microbiome and the host genetics had little value for predicting FA phenotype. Compared to genetic information only, adding rumen bacteria information resulted in a significant improvement of the predictive value for C15:0 from 0.22 to 0.38 (P = 9.50e−07) and C18:3 n-3 from 0 to 0.29 (P = 8.81e−18).

Conclusions

The rumen microbiome has a pronounced influence on the content of odd chain FA and polyunsaturated C18 FA, and to a lesser extent, on the content of the short- and medium-chain FA in the milk of Holstein cattle. The accuracy of prediction of FA phenotypes in milk based on information from either the animal’s genotypes or rumen bacteria composition was very low.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0464-8) contains supplementary material, which is available to authorized users.

A genome-wide association study reveals specific transferases as candidate loci for bovine milk oligosaccharides synthesis

BACKGROUND

Human milk oligosaccharides (OS) play a key role in brain and gut microbiota development of the neonate, but the underlying biosynthetic steps of OS in the mammary gland are still largely unknown. As bovine milk contains OS with somewhat similar structures and functionalities there is increased interest in further understanding the genetic basis underlying the OS content of milk for eventual extraction and generation of value-added ingredients for infant formulas and nutraceuticals. The present study is the first to report on genetic parameter estimation as well as on a genome wide association study (GWAS) from the largest bovine milk OS dataset analyzed to date.

RESULTS

In total 15 different bovine milk OS were monitored. Heritabilities ranged from 0 to 0.68 in Danish Holstein and from 0 to 0.92 in Danish Jersey. The GWAS identified in total 1770 SNPs (FDR < 0.10) for five different OS in Danish Holstein and 6913 SNPs (FDR < 0.10) for 11 OS in Danish Jersey. In Danish Holstein, a major overlapping QTL was identified on BTA1 for LNH and LNT explaining 24% of the variation in these OS. The most significant SNPs were associated with B3GNT5, a gene encoding a glycosyltransferase involved in glycan synthesis. In Danish Jersey, a very strong QTL was detected for the OS with composition 2 Hex 1 HexNAc (isomer 1) on BTA11. The most significant SNP had -log10(P-value) of 52.88 (BOVINEHD1100030300) and was assigned to ABO, a gene encoding ABO blood group glycosyltransferases. This SNP has been reported to be a missense mutation and explains 56% of the OS variation. Other candidate genes of interest identified for milk OS were ALG3, B3GALNT2, LOC520336, PIGV, MAN1C1, ST6GALNAC6, GLT6D1, GALNT14, GALNT17, COLGALT2, LFNG and SIGLEC.

CONCLUSION

To our knowledge, this is the first study documenting a solid breeding potential for bovine milk OS and a strong indication of specific candidate genes related to OS synthesis underlying this genetic influence. This new information has the potential to guide breeding strategies to achieve production of milk with higher diversity and concentration of OS and ultimately facilitate large-scale extraction of bovine milk OS.

Health endowment and later-life outcomes in the labour market: Evidence using genetic risk scores and reduced-form models

This paper examines the relationship between health endowment and later-life outcomes in the labour market. The analysis is based on reduced-form models in which labour market outcomes are regressed on genetic variants related to the increased risk of cardiovascular diseases. We use linked Finnish data that have many strengths. Genetic risk scores constitute exogenous measures for health endowment, and accurate administrative tax records on earnings, employment and social income transfers provide a comprehensive account of an individual's long-term performance in the labour market. The results show that although the direction of an effect is generally consistent with theoretical reasoning, the effects of health endowment on outcomes are statistically weak, and the hypothesis of no effect can be rejected only in one case: genetic endowment related to obesity influences male earnings and employment in prime age. Due to the sample size (N = 1651), the results should be interpreted with caution and should be confirmed in larger samples and in other institutional settings.

A genome-wide association study of mitochondrial DNA copy number in two population-based cohorts

BACKGROUND

Mitochondrial DNA copy number (mtDNA CN) exhibits interindividual and intercellular variation, but few genome-wide association studies (GWAS) of directly assayed mtDNA CN exist. We undertook a GWAS of qPCR-assayed mtDNA CN in the Avon Longitudinal Study of Parents and Children (ALSPAC) and the UK Blood Service (UKBS) cohort. After validating and harmonising data, 5461 ALSPAC mothers (16-43 years at mtDNA CN assay) and 1338 UKBS females (17-69 years) were included in a meta-analysis. Sensitivity analyses restricted to females with white cell-extracted DNA and adjusted for estimated or assayed cell proportions. Associations were also explored in ALSPAC children and UKBS males.

RESULTS

A neutrophil-associated locus approached genome-wide significance (rs709591 [MED24], β (change in SD units of mtDNA CN per allele) [SE] - 0.084 [0.016], p = 1.54e-07) in the main meta-analysis of adult females. This association was concordant in magnitude and direction in UKBS males and ALSPAC neonates. SNPs in and around ABHD8 were associated with mtDNA CN in ALSPAC neonates (rs10424198, β [SE] 0.262 [0.034], p = 1.40e-14), but not other study groups. In a meta-analysis of unrelated individuals (N = 11,253), we replicated a published association in TFAM (β [SE] 0.046 [0.017], p = 0.006), with an effect size much smaller than that observed in the replication analysis of a previous in silico GWAS.

CONCLUSIONS

In a hypothesis-generating GWAS, we confirm an association between TFAM and mtDNA CN and present putative loci requiring replication in much larger samples. We discuss the limitations of our work, in terms of measurement error and cellular heterogeneity, and highlight the need for larger studies to better understand nuclear genomic control of mtDNA copy number.

Does education protect against depression? Evidence from the Young Finns Study using Mendelian randomization

Using participants (N = 1733) drawn from the nationally representative longitudinal Young Finns Study (YFS) we estimate the effect of education on depressive symptoms. In 2007, when the participants were between 30 and 45 years old, they reported their depressive symptoms using a revised version of Beck's Depression Inventory. Education was measured using register information on the highest completed level of education in 2007, which was converted to years of education. To identify a causal relationship between education and depressive symptoms we use an instrumental variables approach (Mendelian randomization, MR) with a genetic risk score as an instrument for years of education. The genetic risk score was based on 74 genetic variants, which were associated with years of education in a genome-wide association study (GWAS). Because the genetic variants are randomly assigned at conception, they induce exogenous variation in years of education and thus identify a causal effect if the assumptions of the MR approach are met. In Ordinary Least Squares (OLS) estimation years of education in 2007 were negatively associated with depressive symptoms in 2007 (b = -0.027, 95% Confidence Interval (CI) = -0.040, -0.015). However, the results based on Mendelian randomization suggested that the effect is not causal (b = 0.017; 95% CI = -0.144, 0.178). This indicates that omitted variables correlated with education and depression may bias the linear regression coefficients and exogenous variation in education caused by differences in genetic make-up does not seem to protect against depressive symptoms.

Genetic polymorphism of sterol transporters in children with future gallstones

BACKGROUND & AIMS

Gallstone disease is related to hypersecretion of cholesterol in bile, and low serum phytosterol levels. We examined how genetic polymorphisms of sterol transporters affect childhood cholesterol metabolism trait predicting adult gallstone disease.

PATIENTS AND METHODS

In retrospective controlled study, we determined D19H polymorphism of ABCG8 gene, genetic variation at Niemann-Pick C1-like 1 (NPC1L1) gene locus (rs41279633, rs17655652, rs2072183, rs217434 and rs2073548), and serum cholesterol, noncholesterol sterols and lipids in children affected by gallstones decades later (n = 66) and controls (n = 126).

RESULTS

In childhood, phytosterols were lower (9.7%-23.4%) in carriers of risk allele 19H compared to 19D homozygotes. Lowest campesterol/cholesterol tertile consisted of 1.9-times more future gallstone subjects, and 3.7-times more 19H carriers than highest one. Campesterol/cholesterol-ratio was highest in 19D homozygote controls, but ∼11% lower in gallstone 19D homozygotes and ∼25% lower among gallstone and control carriers of 19H. Gallstone subjects with alleles CC of rs41279633 and TT of rs217434 of NPC1L1 had ∼18% lower campesterol/cholesterol-ratio compared to mutation carriers.

CONCLUSIONS

Risk trait of cholesterol metabolism (low phytosterols) in childhood favouring cholesterol gallstone disease later in adulthood is influenced by risk variant 19H of ABCG8 and obviously also other factors. NPC1L1 variants have minor influence on noncholesterol sterols.

Pathway-Wide Genetic Risks in Chlamydial Infections Overlap between Tissue Tropisms: A Genome-Wide Association Scan

Chlamydia trachomatis is the most commonly diagnosed bacterial sexually transmitted infection and can lead to tubal factor infertility, a disease characterised by fibrosis of the fallopian tubes. Genetic polymorphisms in molecular pathways involving G protein-coupled receptor signalling, the Akt/PI3K cascade, the mitotic cell cycle, and immune response have been identified in association with the development of trachomatous scarring, an ocular form of chlamydia-related fibrotic pathology. In this case-control study, we performed genome-wide association and pathways-based analysis in a sample of 71 Dutch women who attended an STI clinic who were seropositive for Chlamydia trachomatis antibodies and 169 high-risk Dutch women who sought similar health services but who were seronegative. We identified two regions of within-gene SNP association with Chlamydia trachomatis serological response and found that GPCR signalling and cell cycle pathways were also associated with the trait. These pathway-level associations appear to be common to immunological sequelae of chlamydial infections in both ocular and urogenital tropisms. These pathways may be central mediators of human refractoriness to chlamydial diseases.

Prediction of Adulthood Obesity Using Genetic and Childhood Clinical Risk Factors in the Cardiovascular Risk in Young Finns Study

BACKGROUND

Obesity is a known risk factor for cardiovascular disease. Early prediction of obesity is essential for prevention. The aim of this study is to assess the use of childhood clinical factors and the genetic risk factors in predicting adulthood obesity using machine learning methods.

METHODS AND RESULTS

A total of 2262 participants from the Cardiovascular Risk in YFS (Young Finns Study) were followed up from childhood (age 3-18 years) to adulthood for 31 years. The data were divided into training (n=1625) and validation (n=637) set. The effect of known genetic risk factors (97 single-nucleotide polymorphisms) was investigated as a weighted genetic risk score of all 97 single-nucleotide polymorphisms (WGRS97) or a subset of 19 most significant single-nucleotide polymorphisms (WGRS19) using boosting machine learning technique. WGRS97 and WGRS19 were validated using external data (n=369) from BHS (Bogalusa Heart Study). WGRS19 improved the accuracy of predicting adulthood obesity in training (area under the curve [AUC=0.787 versus AUC=0.744, P<0.0001) and validation data (AUC=0.769 versus AUC=0.747, P=0.026). WGRS97 improved the accuracy in training (AUC=0.782 versus AUC=0.744, P<0.0001) but not in validation data (AUC=0.749 versus AUC=0.747, P=0.785). Higher WGRS19 associated with higher body mass index at 9 years and WGRS97 at 6 years. Replication in BHS confirmed our findings that WGRS19 and WGRS97 are associated with body mass index.

CONCLUSIONS

WGRS19 improves prediction of adulthood obesity. Predictive accuracy is highest among young children (3-6 years), whereas among older children (9-18 years) the risk can be identified using childhood clinical factors. The model is helpful in screening children with high risk of developing obesity.

2SNP heritability and effects of genetic variants for neutrophil-to-lymphocyte and platelet-to-lymphocyte ratio

Neutrophil-to-lymphocyte ratio (NLR) and platelet-to-lymphocyte ratio (PLR) are important biomarkers for disease development and progression. To gain insight into the genetic causes of variance in NLR and PLR in the general population, we conducted genome-wide association (GWA) analyses and estimated SNP heritability in a sample of 5901 related healthy Dutch individuals. GWA analyses identified a new genome-wide significant locus on the HBS1L-MYB intergenic region for PLR, which replicated in a sample of 2538 British twins. For platelet count, we replicated three known genome-wide significant loci in our cohort (at CCDC71L-PIK3CG, BAK1 and ARHGEF3). For neutrophil count, we replicated the PSMD3 locus. For the identified top SNPs, we found significant cis and trans expression quantitative trait loci effects for several loci involved in hematological and immunological pathways. Linkage Disequilibrium score (LD) regression analyses for PLR and NLR confirmed that both traits are heritable, with a polygenetic SNP heritability for PLR of 14.1%, and for NLR of 2.4%. Genetic correlations were present between ratios and the constituent counts, with the genetic correlation (r=0.45) of PLR with platelet count reaching statistical significance. In conclusion, we established that two important biomarkers have a significant heritable SNP component, and identified the first genome-wide locus for PLR.

An interaction map of circulating metabolites, immune gene networks, and their genetic regulation

BACKGROUND

Immunometabolism plays a central role in many cardiometabolic diseases. However, a robust map of immune-related gene networks in circulating human cells, their interactions with metabolites, and their genetic control is still lacking. Here, we integrate blood transcriptomic, metabolomic, and genomic profiles from two population-based cohorts (total N = 2168), including a subset of individuals with matched multi-omic data at 7-year follow-up.

RESULTS

We identify topologically replicable gene networks enriched for diverse immune functions including cytotoxicity, viral response, B cell, platelet, neutrophil, and mast cell/basophil activity. These immune gene modules show complex patterns of association with 158 circulating metabolites, including lipoprotein subclasses, lipids, fatty acids, amino acids, small molecules, and CRP. Genome-wide scans for module expression quantitative trait loci (mQTLs) reveal five modules with mQTLs that have both cis and trans effects. The strongest mQTL is in ARHGEF3 (rs1354034) and affects a module enriched for platelet function, independent of platelet counts. Modules of mast cell/basophil and neutrophil function show temporally stable metabolite associations over 7-year follow-up, providing evidence that these modules and their constituent gene products may play central roles in metabolic inflammation. Furthermore, the strongest mQTL in ARHGEF3 also displays clear temporal stability, supporting widespread trans effects at this locus.

CONCLUSIONS

This study provides a detailed map of natural variation at the blood immunometabolic interface and its genetic basis, and may facilitate subsequent studies to explain inter-individual variation in cardiometabolic disease.

Does higher education protect against obesity? Evidence using Mendelian randomization

OBJECTIVES

The aim of this explorative study was to examine the effect of education on obesity using Mendelian randomization.

METHODS

Participants (N=2011) were from the on-going nationally representative Young Finns Study (YFS) that began in 1980 when six cohorts (aged 30, 33, 36, 39, 42 and 45 in 2007) were recruited. The average value of BMI (kg/m²) measurements in 2007 and 2011 and genetic information were linked to comprehensive register-based information on the years of education in 2007. We first used a linear regression (Ordinary Least Squares, OLS) to estimate the relationship between education and BMI. To identify a causal relationship, we exploited Mendelian randomization and used a genetic score as an instrument for education. The genetic score was based on 74 genetic variants that genome-wide association studies (GWASs) have found to be associated with the years of education. Because the genotypes are randomly assigned at conception, the instrument causes exogenous variation in the years of education and thus enables identification of causal effects.

RESULTS

The years of education in 2007 were associated with lower BMI in 2007/2011 (regression coefficient (b)=-0.22; 95% Confidence Intervals [CI]=-0.29, -0.14) according to the linear regression results. The results based on Mendelian randomization suggests that there may be a negative causal effect of education on BMI (b=-0.84; 95% CI=-1.77, 0.09).

CONCLUSION

The findings indicate that education could be a protective factor against obesity in advanced countries.

A fast and accurate method for detection of IBD shared haplotypes in genome-wide SNP data

Identical by descent (IBD) segments are used to understand a number of fundamental issues in genetics. IBD segments are typically detected using long stretches of identical alleles between haplotypes in phased, whole-genome SNP data. Phase or SNP call errors in genomic data can degrade accuracy of IBD detection and lead to false-positive/negative calls and to under/overextension of true IBD segments. Furthermore, the number of comparisons increases quadratically with sample size, requiring high computational efficiency. We developed a new IBD segment detection program, FISHR (Find IBD Shared Haplotypes Rapidly), in an attempt to accurately detect IBD segments and to better estimate their endpoints using an algorithm that is fast enough to be deployed on very large whole-genome SNP data sets. We compared the performance of FISHR to three leading IBD segment detection programs: GERMLINE, refined IBD, and HaploScore. Using simulated and real genomic sequence data, we show that FISHR is slightly more accurate than all programs at detecting long (>3 cm) IBD segments but slightly less accurate than refined IBD at detecting short (~1 cm) IBD segments. More centrally, FISHR outperforms all programs in determining the true endpoints of IBD segments, which is crucial for several applications of IBD information. FISHR takes two to three times longer than GERMLINE to run, whereas both GERMLINE and FISHR were orders of magnitude faster than refined IBD and HaploScore. Overall, FISHR provides accurate IBD detection in unrelated individuals and is computationally efficient enough to be utilized on large SNP data sets >60 000 individuals.

Estimates of linkage disequilibrium and effective population sizes in Chinese Merino (Xinjiang type) sheep by genome-wide SNPs

Knowledge of linkage disequilibrium (LD) is important for effective genome-wide association studies and accurate genomic prediction. Chinese Merino (Xinjiang type) is well-known fine wool sheep breed. However, the extent of LD across the genome remains unexplored. In this study, we calculated autosomal LD based on genome-wide SNPs of 635 Chinese Merino (Xinjiang type) sheep by Illumina Ovine SNP50 BeadChip. A moderate level of LD (r² ≥ 0.25) across the whole genome was observed at short distances of 0–10 kb. Further, the ancestral effective population size (N_e) was analyzed by extent of LD and found that N_e increased with the increase of generations and declined rapidly within the most recent 50 generations, which is consistent with the history of Chinese Merino sheep breeding, initiated in 1971. We also noted that even when the effective population size was estimated across different single chromosomes, N_e only ranged from 140.36 to 183.33 at five generations in the past, exhibiting a rapid decrease compared with that at ten generations in the past. These results indicated that the genetic diversity in Chinese Merino sheep recently decreased and proper protective measures should be taken to maintain the diversity. Our datasets provided essential genetic information to track molecular variations which potentially contribute to phenotypic variation in Chinese Merino sheep.

Stature and long-term labor market outcomes: Evidence using Mendelian randomization

We use the Young Finns Study (N=∼2000) on the measured height linked to register-based long-term labor market outcomes. The data contain six age cohorts (ages 3, 6, 9, 12, 15 and 18, in 1980) with the average age of 31.7, in 2001, and with the female share of 54.7. We find that taller people earn higher earnings according to the ordinary least squares (OLS) estimation. The OLS models show that 10cm of extra height is associated with 13% higher earnings. We use Mendelian randomization, with the genetic score as an instrumental variable (IV) for height to account for potential confounders that are related to socioeconomic background, early life conditions and parental investments, which are otherwise very difficult to fully account for when using covariates in observational studies. The IV point estimate is much lower and not statistically significant, suggesting that the OLS estimation provides an upward biased estimate for the height premium. Our results show the potential value of using genetic information to gain new insights into the determinants of long-term labor market success.

Genetic variants alter T-bet binding and gene expression in mucosal inflammatory disease

The polarization of CD4+ T cells into distinct T helper cell lineages is essential for protective immunity against infection, but aberrant T cell polarization can cause autoimmunity. The transcription factor T-bet (TBX21) specifies the Th1 lineage and represses alternative T cell fates. Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) that may be causative for autoimmune diseases. The majority of these polymorphisms are located within non-coding distal regulatory elements. It is considered that these genetic variants contribute to disease by altering the binding of regulatory proteins and thus gene expression, but whether these variants alter the binding of lineage-specifying transcription factors has not been determined. Here, we show that SNPs associated with the mucosal inflammatory diseases Crohn's disease, ulcerative colitis (UC) and celiac disease, but not rheumatoid arthritis or psoriasis, are enriched at T-bet binding sites. Furthermore, we identify disease-associated variants that alter T-bet binding in vitro and in vivo. ChIP-seq for T-bet in individuals heterozygous for the celiac disease-associated SNPs rs1465321 and rs2058622 and the IBD-associated SNPs rs1551398 and rs1551399, reveals decreased binding to the minor disease-associated alleles. Furthermore, we show that rs1465321 is an expression quantitative trait locus (eQTL) for the neighboring gene IL18RAP, with decreased T-bet binding associated with decreased expression of this gene. These results suggest that genetic polymorphisms may predispose individuals to mucosal autoimmune disease through alterations in T-bet binding. Other disease-associated variants may similarly act by modulating the binding of lineage-specifying transcription factors in a tissue-selective and disease-specific manner.

High resolution HLA haplotyping by imputation for a British population bioresource

This study aimed to establish the occurrence and frequency of HLA alleles and haplotypes for a healthy British Caucasian population bioresource from Oxfordshire. We present the results of imputation from HLA SNP genotyping data using SNP2HLA for 5553 individuals from Oxford Biobank, defining one- and two-field alleles together with amino acid polymorphisms. We show that this achieves a high level of accuracy with validation using sequence-specific primer amplification PCR. We define six- and eight-locus HLA haplotypes for this population by Bayesian methods implemented using PHASE. We determine patterns of linkage disequilibrium and recombination for these individuals involving classical HLA loci and show how analysis within a haplotype block structure may be more tractable for imputed data. Our findings contribute to knowledge of HLA diversity in healthy populations and further validate future large-scale use of HLA imputation as an informative approach in population bioresources.

Analysis of Genome-Wide Association Data

The last decade has seen substantial advances in the understanding of the genetics of complex traits and disease. This has been largely driven by genome-wide association studies (GWAS), which have identified thousands of genetic loci associated with these traits and disease. This chapter provides a guide on how to perform GWAS on both binary (case-control) and quantitative traits. As poor data quality, through both genotyping failures and unobserved population structure, is a major cause of false-positive genetic associations, there is a particular focus on the crucial steps required to prepare the SNP data prior to analysis. This is followed by the methods used to perform the actual GWAS and visualization of the results.

Estimation of genetic parameters and detection of chromosomal regions affecting the major milk proteins and their post translational modifications in Danish Holstein and Danish Jersey cattle

Background

In the Western world bovine milk products are an important protein source in human diet. The major proteins in bovine milk are the four caseins (CN), α_S1-, α_S2-, β-, and k-CN and the two whey proteins, β-LG and α-LA. It has been shown that both the amount of specific CN and their isoforms including post-translational modifications (PTM) influence technological properties of milk. Therefore, the aim of this study was to 1) estimate genetic parameters for individual proteins in Danish Holstein (DH) (n = 371) and Danish Jersey (DJ) (n = 321) milk, and 2) detect genomic regions associated with specific milk protein and their different PTM forms using a genome-wide association study (GWAS) approach.

Results

For DH, high heritability estimates were found for protein percentage (0.47), casein percentage (0.43), k-CN (0.77), β-LG (0.58), and α-LA (0.40). For DJ, high heritability estimates were found for protein percentage (0.70), casein percentage (0.52), and α-LA (0.44). The heritability for G-k-CN, U-k-CN and GD was higher in the DH compared to the DJ, whereas the heritability for the PD of α_S1-CN was lower in DH compared to DJ, whereas the PD for α_S2-CN was higher in DH compared to DJ. The GWAS results for the main milk proteins were in line what has been earlier published. However, we showed that there were SNPs specifically regulating G-k-CN in DH. Some of these SNPs were assigned to casein protein kinase genes (CSNK1G3, PRKCQ).

Conclusion

The genetic analysis of the major milk proteins and their PTM forms revealed that these were heritable in both DH and DJ. In DH, genomic regions specific for glycosylation of k-CN were detected. Furthermore, genomic regions for the major milk proteins confirmed the regions on BTA6 (casein cluster), BTA11 (PEAP), and BTA14 (DGAT1) as important regions influencing protein composition in milk. The results from this study provide confidence that it is possible to breed for specific milk protein including the different PTM forms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0421-2) contains supplementary material, which is available to authorized users.

Role of Conventional Childhood Risk Factors Versus Genetic Risk in the Development of Type 2 Diabetes and Impaired Fasting Glucose in Adulthood: The Cardiovascular Risk in Young Finns Study

OBJECTIVE

We examined whether the addition of novel genetic risk variant data to conventional childhood risk factors improves risk assessment of impaired fasting glucose (IFG) and type 2 diabetes in adulthood.

RESEARCH DESIGN AND METHODS

An association of a weighted genetic risk score (wGRS) based on 73 risk variants with IFG and type 2 diabetes was analyzed in 2,298 participants of the Cardiovascular Risk in Young Finns Study who were followed for 24-31 years from childhood to adulthood. In addition, the value of the wGRS in pediatric prediction of type 2 diabetes was examined.

RESULTS

Of the 2,298 participants, 484 (21.8%) and 79 (3.4%) had IFG or type 2 diabetes in adulthood, respectively. Adjusting for age, sex, baseline BMI, parental diabetes, mother's BMI, fasting insulin concentration, systolic blood pressure, and smoking status, wGRS was associated with an increased risk of IFG (odds ratio 1.64 [95% CI 1.33-2.01] per unit increase in the wGRS) and type 2 diabetes (2.22 [1.43-3.44]). Incorporating wGRS into pediatric risk models improved model discrimination and reclassification properties. Area under the receiver operating curve improved for IFG (from 0.678 to 0.691, P = 0.015), combined IFG and type 2 diabetes outcome (from 0.678 to 0.692, P = 0.007), and type 2 diabetes (from 0.728 to 0.749, P = 0.158). The net reclassification improvement and integrated discrimination improvement were significant for all outcomes.

CONCLUSIONS

A multifactorial approach combining genetic and clinical risk factors may be useful in identifying children at high risk for adult IFG and type 2 diabetes.

A new model calling procedure for Illumina BeadArray data

Background

Accurate genotype calling for high throughput Illumina data is an important step to extract more genetic information for a large scale genome wide association studies. Many popular calling algorithms use mixture models to infer genotypes of a large number of single nucleotide polymorphisms in a fast and efficient way. In practice, mixture models are mostly restricted to infer genotypes for common SNPs where their minor allele frequencies are quite large. However, it is still challenging to accurately genotype rare variants, especially for some rare variants where the boundaries of their genotypes are not clearly defined.

Results

To further improve the call accuracy and the quality of genotypes on rare variants, a new model calling procedure, named M-D, is proposed to infer genotypes for the Illumina BeadArray data. In this calling procedure, a Gaussian Mixture Model and a Dirichlet Process Gaussian Mixture Model are integrated to infer genotypes.

Conclusions

Applications to Illumina data illustrate that this new approach can improve calling performance compared to other popular genotyping algorithms.

Five endometrial cancer risk loci identified through genome-wide association analysis

We conducted a meta-analysis of three endometrial cancer genome-wide association studies (GWAS) and two follow-up phases totaling 7,737 endometrial cancer cases and 37,144 controls of European ancestry. Genome-wide imputation and meta-analysis identified five new risk loci of genome-wide significance at likely regulatory regions on chromosomes 13q22.1 (rs11841589, near KLF5), 6q22.31 (rs13328298, in LOC643623 and near HEY2 and NCOA7), 8q24.21 (rs4733613, telomeric to MYC), 15q15.1 (rs937213, in EIF2AK4, near BMF) and 14q32.33 (rs2498796, in AKT1, near SIVA1). We also found a second independent 8q24.21 signal (rs17232730). Functional studies of the 13q22.1 locus showed that rs9600103 (pairwise r(2) = 0.98 with rs11841589) is located in a region of active chromatin that interacts with the KLF5 promoter region. The rs9600103[T] allele that is protective in endometrial cancer suppressed gene expression in vitro, suggesting that regulation of the expression of KLF5, a gene linked to uterine development, is implicated in tumorigenesis. These findings provide enhanced insight into the genetic and biological basis of endometrial cancer.

Genetic diversity, linkage disequilibrium, population structure and construction of a core collection of Prunus avium L. landraces and bred cultivars

BACKGROUND

Depiction of the genetic diversity, linkage disequilibrium (LD) and population structure is essential for the efficient organization and exploitation of genetic resources. The objectives of this study were to (i) to evaluate the genetic diversity and to detect the patterns of LD, (ii) to estimate the levels of population structure and (iii) to identify a 'core collection' suitable for association genetic studies in sweet cherry.

RESULTS

A total of 210 genotypes including modern cultivars and landraces from 16 countries were genotyped using the RosBREED cherry 6 K SNP array v1. Two groups, mainly bred cultivars and landraces, respectively, were first detected using STRUCTURE software and confirmed by Principal Coordinate Analysis (PCoA). Further analyses identified nine subgroups using STRUCTURE and Discriminant Analysis of Principal Components (DAPC). Several sub-groups correspond to different eco-geographic regions of landraces distribution. Linkage disequilibrium was evaluated showing lower values than in peach, the reference Prunus species. A 'core collection' containing 156 accessions was selected using the maximum length sub tree method.

CONCLUSION

The present study constitutes the first population genetics analysis in cultivated sweet cherry using a medium-density SNP (single nucleotide polymorphism) marker array. We provided estimations of linkage disequilibrium, genetic structure and the definition of a first INRA's Sweet Cherry core collection useful for breeding programs, germplasm management and association genetics studies.

A proposal for improvement of genotyping performance for ethnically homogeneous population using DNA microarray

DNA microarray is used to determine the genotypes of several hundred thousand to several million SNPs (Single Nucleotide Polymorphisms) on multiple samples at a time. In the conventional method of genotyping using DNA microarray, it is assumed that each SNP has three types of genotypes: two homozygous and one heterozygous genotypes. However, in an ethnically homogeneous population, there are cases when all the samples of a SNP belong to one homozygous genotype, and there are some other cases, especially in the SNPs of low MAFs (Minor Allele Frequencies), each sample belongs to either of the two genotypes: one homozygous and one heterozygous genotypes. In those cases, the conventional method of genotyping may fail to properly determine the genotypes of the samples. In this paper, we propose a new genotyping method, which can be used as a post-processing technique of the conventional genotyping method, for re-judgment of the SNPs having one or two types of genotypes. The proposed method takes fluctuations of the fluorescence intensities of the signals of DNA microarray into account, assigns genotypes to samples from those genotype patterns that may occur under natural mating conditions and applies different genotype judgment methods depending on the number of genotype clusters of a SNP. We evaluate our proposed method using the data of 1000 Genome Project and have found that our proposed method is able to improve the genotyping performance of the conventional method.

A multi-factorial analysis of response to warfarin in a UK prospective cohort

Background

Warfarin is the most widely used oral anticoagulant worldwide, but it has a narrow therapeutic index which necessitates constant monitoring of anticoagulation response. Previous genome-wide studies have focused on identifying factors explaining variance in stable dose, but have not explored the initial patient response to warfarin, and a wider range of clinical and biochemical factors affecting both initial and stable dosing with warfarin.

Methods

A prospective cohort of 711 patients starting warfarin was followed up for 6 months with analyses focusing on both non-genetic and genetic factors. The outcome measures used were mean weekly warfarin dose (MWD), stable mean weekly dose (SMWD) and international normalised ratio (INR) > 4 during the first week. Samples were genotyped on the Illumina Human610-Quad chip. Statistical analyses were performed using Plink and R.

Results

VKORC1 and CYP2C9 were the major genetic determinants of warfarin MWD and SMWD, with CYP4F2 having a smaller effect. Age, height, weight, cigarette smoking and interacting medications accounted for less than 20 % of the variance. Our multifactorial analysis explained 57.89 % and 56.97 % of the variation for MWD and SMWD, respectively. Genotypes for VKORC1 and CYP2C9*3, age, height and weight, as well as other clinical factors such as alcohol consumption, loading dose and concomitant drugs were important for the initial INR response to warfarin. In a small subset of patients for whom data were available, levels of the coagulation factors VII and IX (highly correlated) also played a role.

Conclusion

Our multifactorial analysis in a prospectively recruited cohort has shown that multiple factors, genetic and clinical, are important in determining the response to warfarin. VKORC1 and CYP2C9 genetic polymorphisms are the most important determinants of warfarin dosing, and it is highly unlikely that other common variants of clinical importance influencing warfarin dosage will be found. Both VKORC1 and CYP2C9*3 are important determinants of the initial INR response to warfarin. Other novel variants, which did not reach genome-wide significance, were identified for the different outcome measures, but need replication.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-015-0255-y) contains supplementary material, which is available to authorized users.

Genetic and genomic relationship between methane production measured in breath and fatty acid content in milk samples from Danish Holsteins

In this study the objective was to estimate the genetic and genomic relationship between methane-related traits and milk fatty acid profiles. This was done using two different estimation procedures: a single nucleotide polymorphism-based genomic relationship matrix and a classical pedigree-based relationship matrix. Data was generated on three Danish Holstein herds and a total of 339 cows were available for the study. Methane phenotypes were generated in milking robots during milking over a weekly period and the milk phenotypes were quantified from milk from one milking. Genetic and genomic parameters were estimated using a mixed linear model. Results showed that heritability estimates were comparable between models, but the standard error was lower for genomic heritabilities compared with genetic heritabilities. Genetic as well as genomic correlations were highly variable and had high standard errors, reflecting a similar pattern as for the heritability estimates with lower standard errors for the genomic correlations compared with the pedigree-based genetic correlations. Many of the correlations though had a magnitude that makes further studies on larger datasets worthwhile. The results indicate that genotypes are highly valuable in studies where limited number of phenotypes can be recorded. Also it shows that there is some significant genetic association between methane in the breath of the cow and milk fatty acids profiles.

Novel 6p21.3 Risk Haplotype Predisposes to Acute Coronary Syndrome

BACKGROUND

The HLA-DRB1*01 allele of the human leukocyte antigen has been associated with acute coronary syndrome. Genome-wide association studies have revealed associations with human leukocyte antigen and non-human leukocyte antigen genes of 3 major histocompatibility complex gene classes but not at allelic level.

METHODS AND RESULTS

We conducted a large-scale genetic analysis on a case-control cohort comprising 5376 acute coronary syndrome cases and 4852 unrelated controls from 4 populations of 2 European countries. We analyzed the risk candidate allele of HLA-DRB1*01 by genomic real-time polymerase chain reaction together with high-density single nucleotide polymorphisms of the major histocompatibility complex to precisely identify risk loci for acute coronary syndrome with effective clinical implications. We found a risk haplotype for the disease containing single nucleotide polymorphisms from BTNL2 and HLA-DRA genes and the HLA-DRB1*01 allele. The association of the haplotype appeared in 3 of the 4 populations, and the direction of the effect was consistent in the fourth. Coronary samples from subjects homozygous for the disease-associated haplotype showed higher BTNL2 mRNA levels (r=0.760; P<0.00001).We localized, with immunofluorescence staining, BTNL2 in CD68-positive macrophages of the coronary artery plaques. In homozygous cases, BTNL2 blocking, in T-cell stimulation assays, enhanced CD4(+)FOXP3(+) regulatory T cell proliferation significantly (blocking versus nonblocking; P<0.05).

CONCLUSIONS

In cases with the risk haplotype for acute coronary syndrome, these results suggest involvement of enhanced immune reactions. BTNL2 may have an inhibitory effect on FOXP3(+) T cell proliferation, especially in patients homozygous for the risk alleles.

CLINICAL TRIAL REGISTRATION

https://www.clinicaltrials.gov; Unique Identifier: NCT00417534.