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.     

High mobility group box 1 levels in large vessel vasculitis are not associated with disease activity but are influenced by age and statins

IntroductionTakayasu arteritis (TA) and giant cell arteritis (GCA) are large vessel vasculitides (LVV) that usually present as granulomatous inflammation in arterial walls. High mobility group box 1 (HMGB1) is a nuclear protein that acts as an alarmin when released by dying or activated cells. This study aims to evaluate whether serum HMGB1 can be used as a biomarker in LVV.MethodsTwenty-nine consecutive TA patients with 29 healthy controls (HC) were evaluated in a cross-sectional study. Eighteen consecutive GCA patients with 16 HC were evaluated at the onset of disease and some of them during follow-up. Serum HMGB1 levels were measured by enzyme-linked immunosorbent assay.ResultsIn GCA patients at disease onset mean serum HMGB1 levels did not differ from HC (5.74 ± 4.19 ng/ml vs. 4.17 ± 3.14 ng/ml; p = 0.230). No differences in HMGB1 levels were found between GCA patients with and without polymyalgia rheumatica (p = 0.167), ischemic manifestations (p = 0.873), systemic manifestations (p = 0.474) or relapsing disease (p = 0.608). During follow-up, no significant fluctuations on serum HMGB1 levels were observed from baseline to 3 months (n = 13) (p = 0.075), 12 months (n = 6) (p = 0.093) and at the first relapse (n = 4) (p = 0.202). Serum HMGB1 levels did not differ between TA patients and HC [1.19 (0.45–2.10) ng/ml vs. 1.46 (0.89–3.34) ng/ml; p = 0.181] and no difference was found between TA patients with active disease and in remission [1.31 (0.63–2.16) ng/ml vs. 0.75 (0.39–2.05) ng/ml; p = 0.281]. HMGB1 levels were significantly lower in 16 TA patients on statins compared with 13 patients without statins [0.59 (0.29–1.46) ng/ml vs. 1.93 (0.88–3.34) ng/ml; p = 0.019]. Age was independently associated with higher HMGB1 levels regardless of LVV or control status.ConclusionsPatients with TA and GCA present similar serum HMGB1 levels compared with HC. Serum HMGB1 is not useful to discriminate between active disease and remission. In TA, use of statins was associated with lower HMGB1 levels. HMGB1 is not a biomarker for LVV.     

Reproducibility of MUAP properties in array surface EMG recordings of the upper trapezius and sternocleidomastoid muscle

The use of array surface EMG recordings for detailed assessment of motor control and muscle properties is increasing. Motor unit action potentials (MUAPs) and their properties can be extracted from these recordings. The objective of this study was to determine the reproducibility of variables obtained from array surface EMG recordings of the shoulder and neck muscles during different functional tasks.Eight-channel linear arrays were placed on the upper trapezius (UT) and sternocleidomastoid (SCM) muscles of 12 healthy subjects. Subjects performed 3 tasks: shoulder abduction (90°), ironing (repetitively touching two ends of a horizontal bar in front of the subject), and 90° head turning. The protocol was performed twice while electrodes remained on and repeated a third time a week later.Three global and six MUAP-related variables were calculated. Intra-class correlation coefficients (ICC) were calculated to assess reliability and smallest detectable changes (SDC) were calculated to assess agreement.In general, the EMG variables showed high levels of reliability which suggests they may be effective for differentiating between-subjects. SDC was found to be considerably lower for the frequency-related (5–23%) than for the amplitude-related variables (15–78%), indicating that the frequency-related variables may be more suitable for investigating interventions which aim to modify motor control. There was no difference in reproducibility between global and MUAP-related variables, which justifies their complementary use.     

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.     

Identification of Regions Important for Resistance and Signalling within the Antimicrobial Peptide Transporter BceAB of Bacillus subtilis

In the low-G+C-content Gram-positive bacteria, resistance to antimicrobial peptides is often mediated by so-called resistance modules. These consist of a two-component system and an ATP-binding cassette transporter and are characterized by an unusual mode of signal transduction where the transporter acts as a sensor of antimicrobial peptides, because the histidine kinase alone cannot detect the substrates directly. Thus, the transporters fulfill a dual function as sensors and detoxification systems to confer resistance, but the mechanistic details of these processes are unknown. The paradigm and best-understood example for this is the BceRS-BceAB module of Bacillus subtilis, which mediates resistance to bacitracin, mersacidin, and actagardine. Using a random mutagenesis approach, we here show that mutations that affect specific functions of the transporter BceAB are primarily found in the C-terminal region of the permease, BceB, particularly in the eighth transmembrane helix. Further, we show that while signaling and resistance are functionally interconnected, several mutations could be identified that strongly affected one activity of the transporter but had only minor effects on the other. Thus, a partial genetic separation of the two properties could be achieved by single amino acid replacements, providing first insights into the signaling mechanism of these unusual modules.     

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.     

Stringent and reproducible tetracycline-regulated transgene expression by site-specific insertion at chromosomal loci with pre-characterised induction characteristics

Background  

The ability to regulate transgene expression has many applications, mostly concerning the analysis of gene function. Desirable induction characteristics, such as low un-induced expression, high induced expression and limited cellular heterogeneity, can be seriously impaired by chromosomal position effects at the site of transgene integration. Many clones may therefore need to be screened before one with optimal induction characteristics is identified. Furthermore, such screens must be repeated for each new transgene investigated, and comparisons between clones with different transgenes is complicated by their different integration sites.     

Influence of gravity compensation on muscle activity during reach and retrieval in healthy elderly

IntroductionArm support like gravity compensation may improve arm movements during stroke rehabilitation. It is unknown how gravity compensation affects muscle activation patterns during reach and retrieval movements. Since muscle activity during reach is represented by a component varying with movement velocity and a component supposedly counteracting gravity, we hypothesized that gravity compensation decreases the amplitude of muscle activity, but does not affect the pattern. To examine this, we compared muscle activity during well defined movements with and without gravity compensation in healthy elderly.MethodsTen subjects performed reach and retrieval movements with and without gravity compensation. Muscle activity of biceps, triceps, anterior, middle and posterior parts of deltoid and upper trapezius was compared between the two conditions.ResultsThe level of muscle activity was lower with gravity compensation in all muscles, reaching significance in biceps, anterior deltoid and trapezius (p ? 0.026). The muscle activation pattern did not differ between movements with and without gravity compensation (p ? 0.662).DiscussionGravity compensation only influenced the level of muscle activity but not the muscle activation pattern in terms of timing. Future studies should examine if the influence of gravity compensation is comparable for stroke patients. This may stimulate early and intensive training during rehabilitation.