Set-base dynamical parameter estimation and model invalidation for biochemical reaction networks

Background  

Mathematical modeling and analysis have become, for the study of biological and cellular processes, an important complement to experimental research. However, the structural and quantitative knowledge available for such processes is frequently limited, and measurements are often subject to inherent and possibly large uncertainties. This results in competing model hypotheses, whose kinetic parameters may not be experimentally determinable. Discriminating among these alternatives and estimating their kinetic parameters is crucial to improve the understanding of the considered process, and to benefit from the analytical tools at hand.     

Differential Expression Levels of Integrin α6 Enable the Selective Identification and Isolation of Atrial and Ventricular Cardiomyocytes

Rationale

Central questions such as cardiomyocyte subtype emergence during cardiogenesis or the availability of cardiomyocyte subtypes for cell replacement therapy require selective identification and purification of atrial and ventricular cardiomyocytes. However, current methodologies do not allow for a transgene-free selective isolation of atrial or ventricular cardiomyocytes due to the lack of subtype specific cell surface markers.

Methods and Results

In order to develop cell surface marker-based isolation procedures for cardiomyocyte subtypes, we performed an antibody-based screening on embryonic mouse hearts. Our data indicate that atrial and ventricular cardiomyocytes are characterized by differential expression of integrin α6 (ITGA6) throughout development and in the adult heart. We discovered that the expression level of this surface marker correlates with the intracellular subtype-specific expression of MLC-2a and MLC-2v on the single cell level and thereby enables the discrimination of cardiomyocyte subtypes by flow cytometry. Based on the differential expression of ITGA6 in atria and ventricles during cardiogenesis, we developed purification protocols for atrial and ventricular cardiomyocytes from mouse hearts. Atrial and ventricular identities of sorted cells were confirmed by expression profiling and patch clamp analysis.

Conclusion

Here, we introduce a non-genetic, antibody-based approach to specifically isolate highly pure and viable atrial and ventricular cardiomyocytes from mouse hearts of various developmental stages. This will facilitate in-depth characterization of the individual cellular subsets and support translational research applications.     

Reference genes for quantitative reverse transcription-polymerase chain reaction expression studies in wild and cultivated peanut

Background

Molecular genetic studies on rare tumour entities, such as bone tumours, often require the use of decalcified, formalin-fixed, paraffin-embedded tissue (dFFPE) samples. Regardless of which decalcification procedure is used, this introduces a vast breakdown of DNA that precludes the possibility of further molecular genetic testing. We set out to establish a robust protocol that would overcome these intrinsic hurdles for bone tumour research.

Findings

The goal of our study was to establish a protocol, using a modified DNA isolation procedure and quality controls, to select decalcified samples suitable for array-CGH testing. Archival paraffin blocks were obtained from 9 different pathology departments throughout Europe, using different fixation, embedding and decalcification procedures, in order to preclude a bias for certain lab protocols. Isolated DNA samples were subjected to direct chemical labelling and enzymatic labelling systems and were hybridised on a high resolution oligonucleotide chip containing 44,000 reporter elements. Genomic alterations (gains and losses) were readily detected in most of the samples analysed. For example, both homozygous deletions of 0.6 Mb and high level of amplifications of 0.7 Mb were identified.

Conclusions

We established a robust protocol for molecular genetic testing of dFFPE derived DNA, irrespective of fixation, decalcification or sample type used. This approach may greatly facilitate further genetic testing on rare tumour entities where archival decalcified, formalin fixed samples are the only source.     

Genomic selection of purebred animals for crossbred performance in the presence of dominant gene action

Background

Genomic selection is an appealing method to select purebreds for crossbred performance. In the case of crossbred records, single nucleotide polymorphism (SNP) effects can be estimated using an additive model or a breed-specific allele model. In most studies, additive gene action is assumed. However, dominance is the likely genetic basis of heterosis. Advantages of incorporating dominance in genomic selection were investigated in a two-way crossbreeding program for a trait with different magnitudes of dominance. Training was carried out only once in the simulation.

Results

When the dominance variance and heterosis were large and overdominance was present, a dominance model including both additive and dominance SNP effects gave substantially greater cumulative response to selection than the additive model. Extra response was the result of an increase in heterosis but at a cost of reduced purebred performance. When the dominance variance and heterosis were realistic but with overdominance, the advantage of the dominance model decreased but was still significant. When overdominance was absent, the dominance model was slightly favored over the additive model, but the difference in response between the models increased as the number of quantitative trait loci increased. This reveals the importance of exploiting dominance even in the absence of overdominance. When there was no dominance, response to selection for the dominance model was as high as for the additive model, indicating robustness of the dominance model. The breed-specific allele model was inferior to the dominance model in all cases and to the additive model except when the dominance variance and heterosis were large and with overdominance. However, the advantage of the dominance model over the breed-specific allele model may decrease as differences in linkage disequilibrium between the breeds increase. Retraining is expected to reduce the advantage of the dominance model over the alternatives, because in general, the advantage becomes important only after five or six generations post-training.

Conclusion

Under dominance and without retraining, genomic selection based on the dominance model is superior to the additive model and the breed-specific allele model to maximize crossbred performance through purebred selection.     

Accuracies of genomic breeding values in American Angus beef cattle using K-means clustering for cross-validation

Background

Genomic selection is a recently developed technology that is beginning to revolutionize animal breeding. The objective of this study was to estimate marker effects to derive prediction equations for direct genomic values for 16 routinely recorded traits of American Angus beef cattle and quantify corresponding accuracies of prediction.

Methods

Deregressed estimated breeding values were used as observations in a weighted analysis to derive direct genomic values for 3570 sires genotyped using the Illumina BovineSNP50 BeadChip. These bulls were clustered into five groups using K-means clustering on pedigree estimates of additive genetic relationships between animals, with the aim of increasing within-group and decreasing between-group relationships. All five combinations of four groups were used for model training, with cross-validation performed in the group not used in training. Bivariate animal models were used for each trait to estimate the genetic correlation between deregressed estimated breeding values and direct genomic values.

Results

Accuracies of direct genomic values ranged from 0.22 to 0.69 for the studied traits, with an average of 0.44. Predictions were more accurate when animals within the validation group were more closely related to animals in the training set. When training and validation sets were formed by random allocation, the accuracies of direct genomic values ranged from 0.38 to 0.85, with an average of 0.65, reflecting the greater relationship between animals in training and validation. The accuracies of direct genomic values obtained from training on older animals and validating in younger animals were intermediate to the accuracies obtained from K-means clustering and random clustering for most traits. The genetic correlation between deregressed estimated breeding values and direct genomic values ranged from 0.15 to 0.80 for the traits studied.

Conclusions

These results suggest that genomic estimates of genetic merit can be produced in beef cattle at a young age but the recurrent inclusion of genotyped sires in retraining analyses will be necessary to routinely produce for the industry the direct genomic values with the highest accuracy.     

Breeding value prediction for production traits in layer chickens using pedigree or genomic relationships in a reduced animal model

Background

Genomic selection involves breeding value estimation of selection candidates based on high-density SNP genotypes. To quantify the potential benefit of genomic selection, accuracies of estimated breeding values (EBV) obtained with different methods using pedigree or high-density SNP genotypes were evaluated and compared in a commercial layer chicken breeding line.

Methods

The following traits were analyzed: egg production, egg weight, egg color, shell strength, age at sexual maturity, body weight, albumen height, and yolk weight. Predictions appropriate for early or late selection were compared. A total of 2,708 birds were genotyped for 23,356 segregating SNP, including 1,563 females with records. Phenotypes on relatives without genotypes were incorporated in the analysis (in total 13,049 production records).The data were analyzed with a Reduced Animal Model using a relationship matrix based on pedigree data or on marker genotypes and with a Bayesian method using model averaging. Using a validation set that consisted of individuals from the generation following training, these methods were compared by correlating EBV with phenotypes corrected for fixed effects, selecting the top 30 individuals based on EBV and evaluating their mean phenotype, and by regressing phenotypes on EBV.

Results

Using high-density SNP genotypes increased accuracies of EBV up to two-fold for selection at an early age and by up to 88% for selection at a later age. Accuracy increases at an early age can be mostly attributed to improved estimates of parental EBV for shell quality and egg production, while for other egg quality traits it is mostly due to improved estimates of Mendelian sampling effects. A relatively small number of markers was sufficient to explain most of the genetic variation for egg weight and body weight.     

Sensitive gene expression profiling of human T cell subsets reveals parallel post-thymic differentiation for CD4+ and CD8+ lineages

The differentiation of CD4(+) or CD8(+) T cells following priming of naive cells is central in the establishment of the immune response against pathogens or tumors. However, our understanding of this complex process and the significance of the multiple subsets of differentiation remains controversial. Gene expression profiling has opened new directions of investigation in immunobiology. Nonetheless, the need for substantial amount of biological material often limits its application range. In this study, we have developed procedures to perform microarray analysis on amplified cDNA from low numbers of cells, including primary T lymphocytes, and applied this technology to the study of CD4 and CD8 lineage differentiation. Gene expression profiling was performed on samples of 1000 cells from 10 different subpopulations, defining the major stages of post-thymic CD4(+) or CD8(+) T cell differentiation. Surprisingly, our data revealed that while CD4(+) and CD8(+) T cell gene expression programs diverge at early stages of differentiation, they become increasingly similar as cells reach a late differentiation stage. This suggests that functional heterogeneity between Ag experienced CD4(+) and CD8(+) T cells is more likely to be located early during post-thymic differentiation, and that late stages of differentiation may represent a common end in the development of T-lymphocytes.     

Influence of Clonal Variation on the Efficiency of Arundo donax Propagation Methods

The recent growth of interest in Arundo donax L. (Giant reed) as an energy crop is due to its great vegetative vigor and high biomass productivity. This perennial rhizomatous plant is able to grow in a wide range of pedo-climatic conditions and it has been employed by man for a great number of purposes. This has promoted its worldwide diffusion, despite its sexual sterility. Field establishment represents the most expensive phase of the whole cultivation cycle, because of the high cost of the propagules. In this work three agamic propagation methods, rhizome transplantation, hydroponic cultivation, and in vitro propagation, were compared with the aim of improving their efficiency, and reducing costs and time required. We focused our attention on in vitro and hydroponic cultivation, finding clones characterized by a high response to these propagation methods. In vitro propagation is the most widespread technique adopted in the large-scale production of propagules; according to our results hydroponic cultivation also appears to be a valid and cheap propagation method, making it possible to obtain in about 2 months, vigorous plants ready to be transplanted into the field. We found a strong effect of the clone used on the success rate of the specific propagation method, suggesting that clonal selection could be particularly useful in improving the efficiency of propagation techniques, contributing to cost reduction. Evaluation of field performance of plants propagated by tissue culture is an important matter, considering the possible residual effects of hormones on plants’ growth, and the importance of first year growth for field establishment. These experiments highlighted the lower biomass production from plants obtained by tissue culture. A novel cytokinin-like hormone, meta-topolin, was tested in vitro, in a preliminary study to evaluate its possible use for A. donax L. propagation.

     

Evaluation of coagulation activation after Rhinovirus infection in patients with asthma and healthy control subjects: an observational study

Background

Asthma exacerbations are frequently triggered by rhinovirus infections. Both asthma and respiratory tract infection can activate haemostasis. Therefore we hypothesized that experimental rhinovirus-16 infection and asthmatic airway inflammation act in synergy on the haemostatic balance.

Methods

28 patients (14 patients with mild allergic asthma and 14 healthy non-allergic controls) were infected with low-dose rhinovirus type 16. Venous plasma and bronchoalveolar lavage fluid (BAL fluid) were obtained before and 6 days after infection to evaluate markers of coagulation activation, thrombin-antithrombin complexes, von Willebrand factor, plasmin-antiplasmin complexes, plasminogen activator inhibitor type-1, endogenous thrombin potential and tissue factor-exposing microparticles by fibrin generation test, in plasma and/or BAL fluid. Data were analysed by nonparametric tests (Wilcoxon, Mann Whitney and Spearman correlation).

Results

13 patients with mild asthma (6 females, 19-29 y) and 11 healthy controls (10 females, 19-31 y) had a documented Rhinovirus-16 infection. Rhinovirus-16 challenge resulted in a shortening of the fibrin generation test in BAL fluid of asthma patients (t = -1: 706 s vs. t = 6: 498 s; p = 0.02), but not of controls (t = -1: 693 s vs. t = 6: 636 s; p = 0.65). The fold change in tissue factor-exposing microparticles in BAL fluid inversely correlated with the fold changes in eosinophil cationic protein and myeloperoxidase in BAL fluid after virus infection (r = -0.517 and -0.528 resp., both p = 0.01).Rhinovirus-16 challenge led to increased plasminogen activator inhibitor type-1 levels in plasma in patients with asthma (26.0 ng/mL vs. 11.5 ng/mL in healthy controls, p = 0.04). Rhinovirus-16 load in BAL showed a linear correlation with the fold change in endogenous thrombin potential, plasmin-antiplasmin complexes and plasminogen activator inhibitor type-1.

Conclusions

Experimental rhinovirus infection induces procoagulant changes in the airways of patients with asthma through increased activity of tissue factor-exposing microparticles. These microparticle-associated procoagulant changes are associated with both neutrophilic and eosinophilic inflammation. Systemic activation of haemostasis increases with Rhinoviral load.

Trial registration

This trial was registered at the Dutch trial registry (http://www.trialregister.nl): NTR1677.