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.     

The challenge-phage assay reveals differences in the binding equilibria of mutant Escherichia coli Trp super-repressors in vivo.

The phenotypes of four mutant Escherichia coli Trp repressor proteins with increased activities have been examined in vivo using the challenge-phage assay, an assay based on a positive genetic selection for DNA binding. These proteins, which differ by single amino acid changes from the wild type (Glu13-->Lys, Glu18-->Lys, Glu49-->Lys and Ala77-->Val), require less L-tryptophan than wild-type repressor for activation in vivo, and are super-aporepressors. However, none of the four mutant repressors binds DNA in a corepressor-independent manner. Three of the four mutant repressors (with Glu-->Lys changes) are more active when complexed with tryptophan, and are superholorepressors. Challenge-phage assays with excess tryptophan rank the mutant holorepressors in the same order as determined by binding studies in vitro. Challenge-phage assays with limiting tryptophan reveal additional phenotypic differences among the mutant proteins. These results show that the challenge-phage assay is a robust assay for measuring the relative affinities of specific protein-DNA interactions in vivo.     

Incomplete activation of Escherichia coli hemolysin (HlyA) due to mutations in the 3' region of hlyC.

Mutational analysis of the carboxy-terminal region of Escherichia coli HlyC was performed by site-directed mutagenesis. Replacement of residue Val-127 or Lys-129 reduced the activity of HlyC to about 30 or 60%, respectively, of that of the wild type, while replacement of Gly-128 reduced the activity to less than 1% of the wild-type level. Complete inactivation of HlyC was caused by a double mutation, replacement of Gly-128 with valine and of Lys-129 with isoleucine. Analysis of culture supernatants from mutants with reduced hemolytic activity by two-dimensional gel electrophoresis revealed the production and simultaneous secretion of nonacylated, monoacylated, and fully acylated HlyA forms, demonstrating impairment of the acylation reaction, possibly due to a decreased affinity of HlyC for the individual HlyA acylation sites.     

Proton NMR bandshape studies of lamellar liquid crystals and gel phases containing lecithins and cholesterol.

Proton NMR spectra for gel and liquid crystalline samples, composed of dimyristoyl and/or dipalmitoyl lecithin, cholesterol and water, can be consistently interpreted in terms of mesophase symmetry and molecular diffusion according to a model proposed by Wennerstrom (Wennerstrom, H. (1973) Chem. Phys. Lett. 18, 41-44). It is shown by computer simulation that the characteristic "super-lorentzian" bandshape of the lamellar mesophase can be described by the superposition of three gaussian curves. The NMR signal of the gel phase can be simulated by the superposition of two gaussian curves with widths at half height of 2.5 kHz and 19 kHz. An upper limit of the lateral diffusion coefficient of the lecithin molecules in the gel phase is calculated to be about 5-10(-15) m-2/s. It is therefore concluded that the static intermolecular dipolar couplings average to zero in the lamellar mesophase. An estimation of the order parameter of the liquid crystalline phase is made from experimental data and a calculated "rigid lattice" linewidth. A two phase system is shown to exist in the temperature range 28-34 degrees C for a mesophase of a mixture of dimyristoyl and dipalmitoyl lecithin. The presence of cholesterol results in enhanced lateral diffusion of the lecithin molecules at temperatures below the Chapman transition point.     

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.     

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.