Selecting high-dimensional mixed graphical models using minimal AIC or BIC forests

Background  

Chow and Liu showed that the maximum likelihood tree for multivariate discrete distributions may be found using a maximum weight spanning tree algorithm, for example Kruskal's algorithm. The efficiency of the algorithm makes it tractable for high-dimensional problems.     

Variable generation times and Darwinian fitness measures

Summary Reproductive value (RV) and net reproductive output (R _o) are frequently used fitness measures. We argue that they are only appropriate when intervals between reproductive events are fixed, as they are dimensionless generation-to-generation scalings with units offspring per parent. A fitness measure should account for two different effects of a decrease in generation time: (1) increased survival due to shorter exposure to mortality agents and (2) increased frequency of reproduction.R _o andRV deal with the first of these two effects, while a measure with a physical dimensionper time [T^–1] is needed to account for the second. The Malthusian growth parameter,r, meets this requirement and in situations where time to reproduction is variable, we propose, the instantaneous rate of spread of descendants (from an individual) be used instead ofR _o. As an alternative toRV, we suggest using the instantaneous difference = –r, wherer is the population rate of increase. WhileRV andR _o are dimensionless ratios, , and areper time rates which are appropriate in accounting for alterations in generation time.     

Attachment and Biofilm Forming Capabilities of Staphylococcus epidermidis Strains Isolated from Preterm Infants

Staphylococcus epidermidis, a human commensal, is an important opportunistic, biofilm-forming pathogen and the main cause of late onset sepsis in preterm infants, worldwide. In this study we describe the characteristics of S. epidermidis strains causing late onset (>72 h) bloodstream infection in preterm infants and skin isolates from healthy newborns. Attachment and biofilm formation capability were analyzed in microtiter plates and with transmission electron microscopy (TEM). Clonal relationship among strains was studied with pulsed-field gel electrophoresis. Antimicrobial susceptibility testing was performed, as well as the detection of biofilm-associated genes and of the invasiveness marker IS256 with polymerase chain reaction. Blood and skin isolates had similar attachment and biofilm-forming capabilities and biofilm formation was not related to the presence of specific genes. Filament-like membrane structures were seen by TEM early in the attachment close to the device surface, both in blood and skin strains. Nine of the ten blood isolates contained the IS256 and were also resistant to methicillin and gentamicin in contrast to skin strains. S. epidermidis strains causing bloodstream infection in preterm infants exhibit higher antibiotic resistance and are provided with an invasive genetic equipment compared to skin commensal strains. Adhesion capability to a device surface seems to involve bacterial membrane filaments.     

Enlarging a training set for genomic selection by imputation of un-genotyped animals in populations of varying genetic architecture

Background

The most common application of imputation is to infer genotypes of a high-density panel of markers on animals that are genotyped for a low-density panel. However, the increase in accuracy of genomic predictions resulting from an increase in the number of markers tends to reach a plateau beyond a certain density. Another application of imputation is to increase the size of the training set with un-genotyped animals. This strategy can be particularly successful when a set of closely related individuals are genotyped.

Methods

Imputation on completely un-genotyped dams was performed using known genotypes from the sire of each dam, one offspring and the offspring’s sire. Two methods were applied based on either allele or haplotype frequencies to infer genotypes at ambiguous loci. Results of these methods and of two available software packages were compared. Quality of imputation under different population structures was assessed. The impact of using imputed dams to enlarge training sets on the accuracy of genomic predictions was evaluated for different populations, heritabilities and sizes of training sets.

Results

Imputation accuracy ranged from 0.52 to 0.93 depending on the population structure and the method used. The method that used allele frequencies performed better than the method based on haplotype frequencies. Accuracy of imputation was higher for populations with higher levels of linkage disequilibrium and with larger proportions of markers with more extreme allele frequencies. Inclusion of imputed dams in the training set increased the accuracy of genomic predictions. Gains in accuracy ranged from close to zero to 37.14%, depending on the simulated scenario. Generally, the larger the accuracy already obtained with the genotyped training set, the lower the increase in accuracy achieved by adding imputed dams.

Conclusions

Whenever a reference population resembling the family configuration considered here is available, imputation can be used to achieve an extra increase in accuracy of genomic predictions by enlarging the training set with completely un-genotyped dams. This strategy was shown to be particularly useful for populations with lower levels of linkage disequilibrium, for genomic selection on traits with low heritability, and for species or breeds for which the size of the reference population is limited.     

Bacterial colonization of Hydra hatchlings follows a robust temporal pattern

Animals are colonized by complex bacterial communities. The processes controlling community membership and influencing the establishment of the microbial ecosystem during development are poorly understood. Here we aimed to explore the assembly of bacterial communities in Hydra with the broader goal of elucidating the general rules that determine the temporal progression of bacterial colonization of animal epithelia. We profiled the microbial communities in polyps at various time points after hatching in four replicates. The composition and temporal patterns of the bacterial communities were strikingly similar in all replicates. Distinct features included high diversity of community profiles in the first week, a remarkable but transient adult-like profile 2 weeks after hatching, followed by progressive emergence of a stable adult-like pattern characterized by low species diversity and the preponderance of the Betaproteobacterium Curvibacter. Intriguingly, this process displayed important parallels to the assembly of human fecal communities after birth. In addition, a mathematical modeling approach was used to uncover the organizational principles of this colonization process, suggesting that both, local environmental or host-derived factor(s) modulating the colonization rate, as well as frequency-dependent interactions of individual bacterial community members are important aspects in the emergence of a stable bacterial community at the end of development.     

Dysbiosis of upper respiratory tract microbiota in elderly pneumonia patients

Bacterial pneumonia is a major cause of morbidity and mortality in elderly. We hypothesize that dysbiosis between regular residents of the upper respiratory tract (URT) microbiome, that is balance between commensals and potential pathogens, is involved in pathogen overgrowth and consequently disease. We compared oropharyngeal microbiota of elderly pneumonia patients (n=100) with healthy elderly (n=91) by 16S-rRNA-based sequencing and verified our findings in young adult pneumonia patients (n=27) and young healthy adults (n=187). Microbiota profiles differed significantly between elderly pneumonia patients and healthy elderly (PERMANOVA, P<0.0005). Highly similar differences were observed between microbiota profiles of young adult pneumonia patients and their healthy controls. Clustering resulted in 11 (sub)clusters including 95% (386/405) of samples. We observed three microbiota profiles strongly associated with pneumonia (P<0.05) and either dominated by lactobacilli (n=11), Rothia (n=51) or Streptococcus (pseudo)pneumoniae (n=42). In contrast, three other microbiota clusters (in total n=183) were correlated with health (P<0.05) and were all characterized by more diverse profiles containing higher abundances of especially Prevotella melaninogenica, Veillonella and Leptotrichia. For the remaining clusters (n=99), the association with health or disease was less clear. A decision tree model based on the relative abundance of five bacterial community members in URT microbiota showed high specificity of 95% and sensitivity of 84% (89% and 73%, respectively, after cross-validation) for differentiating pneumonia patients from healthy individuals. These results suggest that pneumonia in elderly and young adults is associated with dysbiosis of the URT microbiome with bacterial overgrowth of single species and absence of distinct anaerobic bacteria. Whether the observed microbiome changes are a cause or a consequence of the development of pneumonia or merely coincide with disease status remains a question for future research.     

Does small-sided-games’ court area influence metabolic,perceptual, and physical performance parameters of young elite basketball players?

The purpose of this study was to investigate the effect of court size on physiological responses and physical performance of young elite basketball players. Twelve male basketball players (18.6 ± 0.5 years; 88.8 ± 14.5 kg; 192.6 ± 6.5 cm) from an under-19 team performed two small-sided games (matches) with different court areas (28x15 m and 28x9 m; 28x15 and 28x9 protocols). The number of players (3x3) was kept the same in each protocol. The players performed a repeated-sprint ability (RSA) test before and after each match. Blood lactate concentration was collected before (pre) and after (post) the matches, and the session rating of perceived exertion (session-RPE) was determined 30 minutes after the match. Best and mean time in the RSA test were not different between the 28x15 and the 28x9 match protocols (p > 0.05). A significant difference was observed for lactate concentration from pre- to post-match (p < 0.05) in both protocols (28x15 and 28x9); however, there was no significant interaction between protocols. A similar session-RPE mean score (28x15: 7.2 ± 1.4 and 28x9: 6.6 ± 1.4) was detected for both protocols (p > 0.05, ES=0.41). In summary, the results of the current study suggest that the different court areas induced similar responses. Although there was no significant difference in effort perception, players tended to perceive a greater effort in the larger court size.     

An optimized MALDI MSI protocol for spatial detection of tryptic peptides in fresh frozen prostate tissue

MALDI MS imaging (MSI) is a powerful analytical tool for spatial peptide detection in heterogeneous tissues. Proper sample preparation is crucial to achieve high quality, reproducible measurements. Here we developed an optimized protocol for spatially resolved proteolytic peptide detection with MALDI time-of-flight MSI of fresh frozen prostate tissue sections. The parameters tested included four different tissue washes, four methods of protein denaturation, four methods of trypsin digestion (different trypsin densities, sprayers, and incubation times), and five matrix deposition methods (different sprayers, settings, and matrix concentrations). Evaluation criteria were the number of detected and excluded peaks, percentage of high mass peaks, signal-to-noise ratio, spatial localization, and average intensities of identified peptides, all of which were integrated into a weighted quality evaluation scoring system. Based on these scores, the optimized protocol included an ice-cold EtOH+H₂O wash, a 5 min heating step at 95°C, tryptic digestion incubated for 17h at 37°C and CHCA matrix deposited at a final amount of 1.8 μg/mm². Including a heat-induced protein denaturation step after tissue wash is a new methodological approach that could be useful also for other tissue types. This optimized protocol for spatial peptide detection using MALDI MSI facilitates future biomarker discovery in prostate cancer and may be useful in studies of other tissue types.     

Modelling spatial dynamics of fish

Our ability to model spatial distributions of fish populations is reviewed by describing the available modelling tools. Ultimate models of the individual's motivation for behavioural decisions are derived from evolutionary ecology. Mechanistic models for how fish sense and may respond to their surroundings are presented for vision, olfaction, hearing, the lateral line and other sensory organs. Models for learning and memory are presented, based both upon evolutionary optimization premises and upon neurological information processing and decision making. Functional tools for modelling behaviour and life histories can be categorized as belonging to an optimization or an adaptation approach. Among optimization tools, optimal foraging theory, life history theory, ideal free distribution, game theory and stochastic dynamic programming are presented. Among adaptation tools, genetic algorithms and the combination with artificial neural networks are described. The review advocates the combination of evolutionary and neurological approaches to modelling spatial dynamics of fish.     

The DiversiLab system versus pulsed-field gel electrophoresis: Characterisation of extended spectrum β-lactamase producing Escherichia coli and Klebsiella pneumoniae

Fast and reliable epidemiological typing methods for identifying outbreaks and epidemic strains of extended spectrum β-lactamase (ESBL) producing Enterobacteriaceae are urgently needed. The DiversiLab system (DL) has been proposed for these purposes. We compared DL to pulsed-field gel electrophoresis (PFGE) on a national collection of ESBL-producing Escherichia coli (n = 258; of which 226 isolates were typeable with PFGE) and Klebsiella pneumoniae (n = 48) isolated in 2007. For E. coli the Wallace coefficients showed that the probability of two isolates of the same DL type having the same PFGE type was only 19.8% and the probability of two isolates of the same PFGE type having the same DL type was 90.4%. For K. pneumoniae the Wallace coefficients showed that the probability of two isolates of the same DL type having the same PFGE type was 100% and the probability of two isolates of the same PFGE type having the same DL type was 79%, indicating that for this K. pneumoniae strain collection DL was slightly more discriminatory. Only four of 48 isolates had discordant results with the two methods. In E. coli 42% of the isolates were sequence type 131 and these isolates were related at > 95% similarity with DL and at ≥ 60% similarity with PFGE. In summary, for E. coli DL performed well in identifying isolates related by PFGE, but overestimated the genetic relatedness in the studied collection. This indicates that DL could be a primary screening method for excluding unrelated isolates. Isolates shown to be related must be confirmed with a more discriminatory method. For K. pneumoniae, DL discriminated well but overestimated the diversity of the isolates compared to PFGE, assuming a risk of missing possible genetic relatedness.