Predicting expected progeny difference for marbling score in Angus cattle using artificial neural networks and Bayesian regression models

Background

Artificial neural networks (ANN) mimic the function of the human brain and are capable of performing massively parallel computations for data processing and knowledge representation. ANN can capture nonlinear relationships between predictors and responses and can adaptively learn complex functional forms, in particular, for situations where conventional regression models are ineffective. In a previous study, ANN with Bayesian regularization outperformed a benchmark linear model when predicting milk yield in dairy cattle or grain yield of wheat. Although breeding values rely on the assumption of additive inheritance, the predictive capabilities of ANN are of interest from the perspective of their potential to increase the accuracy of prediction of molecular breeding values used for genomic selection. This motivated the present study, in which the aim was to investigate the accuracy of ANN when predicting the expected progeny difference (EPD) of marbling score in Angus cattle. Various ANN architectures were explored, which involved two training algorithms, two types of activation functions, and from 1 to 4 neurons in hidden layers. For comparison, BayesCπ models were used to select a subset of optimal markers (referred to as feature selection), under the assumption of additive inheritance, and then the marker effects were estimated using BayesCπ with π set equal to zero. This procedure is referred to as BayesCpC and was implemented on a high-throughput computing cluster.

Results

The ANN with Bayesian regularization method performed equally well for prediction of EPD as BayesCpC, based on prediction accuracy and sum of squared errors. With the 3K-SNP panel, for example, prediction accuracy was 0.776 using BayesCpC, and ranged from 0.776 to 0.807 using BRANN. With the selected 700-SNP panel, prediction accuracy was 0.863 for BayesCpC and ranged from 0.842 to 0.858 for BRANN. However, prediction accuracy for the ANN with scaled conjugate gradient back-propagation was lower, ranging from 0.653 to 0.689 with the 3K-SNP panel, and from 0.743 to 0.793 with the selected 700-SNP panel.

Conclusions

ANN with Bayesian regularization performed as well as linear Bayesian regression models in predicting additive genetic values, supporting the idea that ANN are useful as universal approximators of functions of interest in breeding contexts.     

What does biologically meaningful mean? A perspective on gene regulatory network validation

Gene regulatory networks (GRNs) are rapidly being delineated, but their quality and biological meaning are often questioned. Here, I argue that biological meaning is challenging to define and discuss reasons why GRN validation should be interpreted cautiously.     

Meta-analysis identified the TNFA -308G > A promoter polymorphism as a risk factor for disease severity in patients with rheumatoid arthritis

Introduction

The goal of this study is to investigate whether the -308G > A promoter polymorphism in the tumor necrosis factor alpha (TNFA) gene is associated with disease severity and radiologic joint damage in a large cohort of patients with rheumatoid arthritis (RA).

Methods

A long-term observational early RA inception cohort (n = 208) with detailed information about disease activity and radiologic damage after 3, 6 and 9 years of disease was genotyped for the TNFA -308G > A promoter polymorphism (rs1800629). A longitudinal regression analysis was performed to assess the effect of genotype on RA disease severity and joint damage. Subsequently, a meta-analysis, including all publically available data, was performed to further test the association between joint erosions and the TNFA polymorphism. To learn more about the mechanism behind the effect of the polymorphism, RNA isolated from peripheral blood from RA patients (n = 66) was used for TNFA gene expression analysis by quantitative PCR.

Results

Longitudinal regression analysis with correction for gender and disease activity showed a significant difference in total joint damage between GG and GA+AA genotype groups (P = 0.002), which was stable over time. The meta-analysis, which included 2,053 patients, confirmed an association of the genetic variant with the development of erosions (odds ratio 0.78, 95% CI 0.62, 0.98). No significant differences in TNFA gene expression were observed for the different genotypes, confirming earlier findings in healthy individuals.

Conclusions

Our data confirm that the TNFA -308G > A promoter polymorphism is associated with joint damage in patients with RA. This is not mediated by differences in TNFA gene expression between genotypes.     

A systematic screen for genes expressed in definitive endoderm by Serial Analysis of Gene Expression (SAGE)

Background  

The embryonic definitive endoderm (DE) gives rise to organs of the gastrointestinal and respiratory tract including the liver, pancreas and epithelia of the lung and colon. Understanding how DE progenitor cells generate these tissues is critical to understanding the cause of visceral organ disorders and cancers, and will ultimately lead to novel therapies including tissue and organ regeneration. However, investigation into the molecular mechanisms of DE differentiation has been hindered by the lack of early DE-specific markers.     

Multidisciplinary transmural rehabilitation for older persons with a stroke: the design of a randomised controlled trial

Background

Stroke is one of the major causes of loss of independence, decreased quality of life and mortality among elderly people. About half of the elderly stroke patients discharged after rehabilitation in a nursing home still experience serious impairments in daily functioning one year post stroke, which can lead to difficulties in picking up and managing their social life. The aim of this study is to evaluate the effectiveness and feasibility of a new multidisciplinary transmural rehabilitation programme for older stroke patients.

Methods

A two group multicentre randomised controlled trial is used to evaluate the effects of the rehabilitation programme. The programme consists of three care modules: 1) neurorehabilitation treatment for elderly stroke patients; 2) empowerment training for patient and informal caregiver; and 3) stroke education for patient and informal caregiver. The total programme has a duration of between two and six months, depending on the individual problems of the patient and informal caregiver. The control group receives usual care in the nursing home and after discharge.Patients aged 65 years and over are eligible for study participation when they are admitted to a geriatric rehabilitation unit in a nursing home due to a recent stroke and are expected to be able to return to their original home environment after discharge. Data are gathered by face-to-face interviews, self-administered questionnaires, focus groups and registration forms. Primary outcomes for patients are activity level after stroke, functional dependence, perceived quality of life and social participation. Outcomes for informal caregivers are perceived care burden, objective care burden, quality of life and perceived health. Outcome measures of the process evaluation are implementation fidelity, programme deliverance and the opinion of the stroke professionals, patients and informal caregivers about the programme. Outcome measures of the economic evaluation are the healthcare utilisation and associated costs. Data are collected at baseline, and after six and 12 months. The first results of the study will be expected in 2014.

Trial registration

International Standard Randomised Controlled Trial Register Number ISRCTN62286281, The Dutch Trial Register NTR2412

     

Mitochondrial growth and division during the cell cycle in HeLa cells

The growth and division of mitochondria during the cell cycle was investigated by a morphometric analysis of electron micrographs of synchronized HeLa cells. The ratio of total outer membrane contour length to cytoplasmic area did not vary significantly during the cell cycle, implying a continuous growth of the mitochondrial outer membrane. The mean fraction of cytoplasmic area occupied by mitochondrial profiles was likewise found to remain constant, indicating that the increase in total mitochondrial volume per cell occurs continuously during interphase, in such a way that the mitochondrial complement occupies a constant fraction( approximately 10-11(percent)) of the volume of the cytoplasm. The mean area, outer membrane contour length, and axis ratio of the mitochondrial profiles also did not vary appreciably during the cell cycle; furthermore, the close similarity of the frequency distributions of these parameters for the six experimental time-points suggested a stable mitochondrial shape distribution. The constancy of both the mean mitochondrial profile area and the number of mitochondrial profiles per unit of cytoplasmic area was interpreted to indicate the continuous division of mitochondria at the level of the cell population. Furthermore, no evidence was found for the occurrence of synchronous mitochondrial growth and division within individual cells. Thus, it appears that, in HeLa cells, there is no fixed temporal relationship between the growth and division of mitochondria and the events of the cell cycle. A number of statistical methods were developed for the purpose of making numerical estimates of certain three-dimensional cellular and mitochondrial parameters. Mean cellular and cytoplasmic volumes were calculated for the six time-points; both exhibited a nonlinear, approx. twofold increase. A comparison of the axis ratio distributions of the mitochondrial profiles with theoretical distributions expected from random sectioning of bodies of various three-dimensional shapes allowed the derivation of an "average" mitochondrial shape. This, in turn, permitted calculations to be made which expressed the two-dimensional results in three-dimensional terms. Thus, the estimated values for the number of mitochondria per unit of cytoplasmic volume and for the mean mitochondrial volume were found to remain constant during the cell cycle, while the estimated number of mitochondria per cell increase approx. twofold in an essentially continuous manner.     

Contribution of exogenous genetic elements to the group A Streptococcus metagenome

Variation in gene content among strains of a bacterial species contributes to biomedically relevant differences in phenotypes such as virulence and antimicrobial resistance. Group A Streptococcus (GAS) causes a diverse array of human infections and sequelae, and exhibits a complex pathogenic behavior. To enhance our understanding of genotype-phenotype relationships in this important pathogen, we determined the complete genome sequences of four GAS strains expressing M protein serotypes (M2, M4, and 2 M12) that commonly cause noninvasive and invasive infections. These sequences were compared with eight previously determined GAS genomes and regions of variably present gene content were assessed. Consistent with the previously determined genomes, each of the new genomes is approximately 1.9 Mb in size, with approximately 10% of the gene content of each encoded on variably present exogenous genetic elements. Like the other GAS genomes, these four genomes are polylysogenic and prophage encode the majority of the variably present gene content of each. In contrast to most of the previously determined genomes, multiple exogenous integrated conjugative elements (ICEs) with characteristics of conjugative transposons and plasmids are present in these new genomes. Cumulatively, 242 new GAS metagenome genes were identified that were not present in the previously sequenced genomes. Importantly, ICEs accounted for 41% of the new GAS metagenome gene content identified in these four genomes. Two large ICEs, designated 2096-RD.2 (63 kb) and 10750-RD.2 (49 kb), have multiple genes encoding resistance to antimicrobial agents, including tetracycline and erythromycin, respectively. Also resident on these ICEs are three genes encoding inferred extracellular proteins of unknown function, including a predicted cell surface protein that is only present in the genome of the serotype M12 strain cultured from a patient with acute poststreptococcal glomerulonephritis. The data provide new information about the GAS metagenome and will assist studies of pathogenesis, antimicrobial resistance, and population genomics.