Imputation of non-genotyped individuals based on genotyped relatives: assessing the imputation accuracy of a real case scenario in dairy cattle

Background

Imputation of genotypes for ungenotyped individuals could enable the use of valuable phenotypes created before the genomic era in analyses that require genotypes. The objective of this study was to investigate the accuracy of imputation of non-genotyped individuals using genotype information from relatives.

Methods

Genotypes were simulated for all individuals in the pedigree of a real (historical) dataset of phenotyped dairy cows and with part of the pedigree genotyped. The software AlphaImpute was used for imputation in its standard settings but also without phasing, i.e. using basic inheritance rules and segregation analysis only. Different scenarios were evaluated i.e.: (1) the real data scenario, (2) addition of genotypes of sires and maternal grandsires of the ungenotyped individuals, and (3) addition of one, two, or four genotyped offspring of the ungenotyped individuals to the reference population.

Results

The imputation accuracy using AlphaImpute in its standard settings was lower than without phasing. Including genotypes of sires and maternal grandsires in the reference population improved imputation accuracy, i.e. the correlation of the true genotypes with the imputed genotype dosages, corrected for mean gene content, across all animals increased from 0.47 (real situation) to 0.60. Including one, two and four genotyped offspring increased the accuracy of imputation across all animals from 0.57 (no offspring) to 0.73, 0.82, and 0.92, respectively.

Conclusions

At present, the use of basic inheritance rules and segregation analysis appears to be the best imputation method for ungenotyped individuals. Comparison of our empirical animal-specific imputation accuracies to predictions based on selection index theory suggested that not correcting for mean gene content considerably overestimates the true accuracy. Imputation of ungenotyped individuals can help to include valuable phenotypes for genome-wide association studies or for genomic prediction, especially when the ungenotyped individuals have genotyped offspring.     

Exploring causal networks of bovine milk fatty acids in a multivariate mixed model context

Background

Knowledge regarding causal relationships among traits is important to understand complex biological systems. Structural equation models (SEM) can be used to quantify the causal relations between traits, which allow prediction of outcomes to interventions applied to such a network. Such models are fitted conditionally on a causal structure among traits, represented by a directed acyclic graph and an Inductive Causation (IC) algorithm can be used to search for causal structures. The aim of this study was to explore the space of causal structures involving bovine milk fatty acids and to select a network supported by data as the structure of a SEM.

Results

The IC algorithm adapted to mixed models settings was applied to study 14 correlated bovine milk fatty acids, resulting in an undirected network. The undirected pathway from C4:0 to C12:0 resembled the de novo synthesis pathway of short and medium chain saturated fatty acids. By using prior knowledge, directions were assigned to that part of the network and the resulting structure was used to fit a SEM that led to structural coefficients ranging from 0.85 to 1.05. The deviance information criterion indicated that the SEM was more plausible than the multi-trait model.

Conclusions

The IC algorithm output pointed towards causal relations between the studied traits. This changed the focus from marginal associations between traits to direct relationships, thus towards relationships that may result in changes when external interventions are applied. The causal structure can give more insight into underlying mechanisms and the SEM can predict conditional changes due to such interventions.     

Depression,Anxiety and Glucose Metabolism in the General Dutch Population: The New Hoorn Study

Background

There is a well recognized association between depression and diabetes. However, there is little empirical data about the prevalence of depressive symptoms and anxiety among different groups of glucose metabolism in population based samples. The aim of this study was to determine whether the prevalence of increased levels of depression and anxiety is different between patients with type 2 diabetes and subjects with impaired glucose metabolism (IGM) and normal glucose metabolism (NGM).

Methodology/Principal Findings

Cross-sectional data from a population-based cohort study of 2667 residents, 1261 men and 1406 women aged 40–65 years from the Hoorn region, the Netherlands. Depressive symptoms and anxiety were measured using the Centre for Epidemiologic Studies Depression Scale (CES-D, score ≥16) and the Hospital Anxiety and Depression Scale – Anxiety Subscale (HADS-A, score ≥8), respectively. Glucose metabolism status was determined by oral glucose tolerance test. In the total study population the prevalence of depressive symptoms and anxiety for the NGM, IGM and type 2 diabetes were 12.5, 12.2 and 21.0% (P = 0.004) and 15.0, 15.3 and 19.9% (p = 0.216), respectively. In men, the prevalence of depressive symptoms was 7.7, 9.5 and 19.6% (p<0.001), and in women 16.4, 15.8 and 22.6 (p = 0.318), for participants with NGM, IGM and type 2 diabetes, respectively. Anxiety was not associated with glucose metabolism when stratified for sex. Intergroup differences (NGM vs. IGM and IGM vs. type 2 diabetes) revealed that higher prevalences of depressive symptoms are mainly manifested in participants with type 2 diabetes, and not in participants with IGM.

Conclusions

Depressive symptoms, but not anxiety are associated with glucose metabolism. This association is mainly determined by a higher prevalence of depressive symptoms in participants with type 2 diabetes and not in participants with IGM.     

Plasma globotriaosylsphingosine: Diagnostic value and relation to clinical manifestations of Fabry disease

Fabry disease is an X-linked lysosomal storage disorder due to deficiency of alpha-Galactosidase A, causing accumulation of globotriaosylceramide and elevated plasma globotriaosylsphingosine (lysoGb3). The diagnostic value and clinical relevance of plasma lysoGb3 concentration was investigated. All male and adult female patients with classical Fabry disease could be discerned by an elevated plasma lysoGb3. In young pre-symptomatic Fabry heterozygotes, lysoGb3 levels can be normal. Individuals carrying the R112H and P60L mutations, without classical Fabry symptoms, showed no elevated plasma lysoGb3. Multiple regression analysis showed that there is no correlation of plasma lysoGb3 concentration with total disease severity score in Fabry males. However, plasma lysoGb3 concentration did correlate with white matter lesions (odds ratio: 6.1 per 100 nM lysoGb3 increase (95% CI: 1.4–25.9, p = 0.015). In females, plasma lysoGb3 concentration correlated with overall disease severity. Furthermore, plasma lysoGb3 level was related to left ventricular mass (19.5 ± 5.5 g increase per 10 nM lysoGb3 increase; p = 0.001). In addition, it was assessed whether lifetime exposure to lysoGb3 correlates with disease manifestations. Male Fabry patients with a high lysoGb3 exposure (> 10,000 U), were moderately or severely affected, only one mildly. Female patients with a low exposure (< 1000 U) were asymptomatic or mildly affected. A large proportion of the females with an exposure > 1000 U showed disease complications. Plasma lysoGb3 is useful for the diagnosis of Fabry disease. LysoGb3 is an independent risk factor for development of cerebrovascular white matter lesions in male patients and left ventricular hypertrophy in females. Disease severity correlates with exposure to plasma lysoGb3.     

Drivers of recombinant soluble influenza A virus hemagglutinin and neuraminidase expression in mammalian cells

Recombinant soluble trimeric influenza A virus hemagglutinins (HA) and tetrameric neuraminidases (NAs) have proven to be excellent tools to decipher biological properties. Receptor binding and sialic acid cleavage by recombinant proteins correlate satisfactorily compared to whole viruses. Expression of HA and NA can be achieved in a plethora of different laboratory hosts. For immunological and receptor interaction studies however, insect and mammalian cell expressed proteins are preferred due to the presence of N‐linked glycosylation and disulfide bond formation. Because mammalian‐cell expression is widely applied, an increased expression yield is an important goal. Here we report that using codon‐optimized genes and sfGFP fusions, the expression yield of HA can be significantly improved. sfGFP also significantly increased expression yields when fused to the N‐terminus of NA. In this study, a suite of different hemagglutinin and neuraminidase constructs are described, which can be valuable tools to study a wide array of different HAs, NAs and their mutants.     

The phenotype and genotype of rheumatoid arthritis in the Democratic Republic of Congo

Introduction

Little is known about rheumatoid arthritis in the black, particularly in Congolese, populations. Our objective was to describe the phenotype and genotype of rheumatoid arthritis (RA) in Congolese.

Methods

All consecutive rheumatoid arthritis (RA) patients attending Kinshasa University Hospital in a three-year time period were included. Demographics, clinical features and tobacco consumption were noted. Disease Activity Score (DAS)-28 based on the erythrocyte sedimentation rate (ESR), Health Assessment Questionnaire (HAQ), anti-citrullinated peptide antibodies (CCP) antibodies and rheumatoid factor (RF) were determined. Radiographs were scored according to Sharp-van der Heijde. On a subset of patients and controls HLA-DRB1 typing was performed.

Results

A total of 114 females and 14 males aged 51.2 ± 14.9 were included. Mean duration of symptoms was four years. Moderate tobacco consumption was reported in a minority of patients. DAS-28 at first visit was >5.1 and HAQ ≥0.5 in all patients. X-rays showed joint erosions and/or joint space narrowing, mostly of a moderate grade in 55.8% of patients. Anti-CCP and/or RF were present in 48.6% of patients with available data (n = 72) and in 3.0% of controls (n = 67). Radiographic changes and nodules were more frequent in RF or anti-CCP positive patients. One copy of the shared epitope was found in 13 patients (35.1%) and 3 controls (12.5%). Two copies were found in one patient (2.7%) and in one control (4.2%).

Conclusion

Congolese patients with RA consult long after disease onset. Despite this delay, the majority presents without major damage and is RF, anti-CCP and SE negative. We put forward the hypothesis that besides different environmental factors there is probably also a particular genetic risk profile in Congolese patients, different from the HLA-DRB1 shared epitope.