Centralization Within Sub-Experiments Enhances the Biological Relevance of Gene Co-expression Networks: A Plant Mitochondrial Case Study

Gene co-expression networks (GCNs) can be prepared using a variety of mathematical approaches based on data sampled across diverse developmental processes, tissue types, pathologies, mutant backgrounds, and stress conditions. These networks are used to identify genes with similar expression dynamics but are prone to introducing false-positive and false-negative relationships, especially in the instance of large and heterogenous datasets. With the aim of optimizing the relevance of edges in GCNs and enhancing global biological insight, we propose a novel approach that involves a data-centering step performed simultaneously per gene and per sub-experiment, called centralization within sub-experiments (CSE). Using a gene set encoding the plant mitochondrial proteome as a case study, our results show that all CSE-based GCNs assessed had significantly more edges within the majority of the considered functional sub-networks, such as the mitochondrial electron transport chain and its complexes, than GCNs not using CSE; thus demonstrating that CSE-based GCNs are efficient at predicting canonical functions and associated pathways, here referred to as the core gene network. Furthermore, we show that correlation analyses using CSE-processed data can be used to fine-tune prediction of the function of uncharacterized genes; while its use in combination with analyses based on non-CSE data can augment conventional stress analyses with the innate connections underpinning the dynamic system being examined. Therefore, CSE is an effective alternative method to conventional batch correction approaches, particularly when dealing with large and heterogenous datasets. The method is easy to implement into a pre-existing GCN analysis pipeline and can provide enhanced biological relevance to conventional GCNs by allowing users to delineate a core gene network.

AUTHOR SUMMARY

Gene co-expression networks (GCNs) are the product of a variety of mathematical approaches that identify causal relationships in gene expression dynamics but are prone to the misdiagnoses of false-positives and false-negatives, especially in the instance of large and heterogenous datasets. In light of the burgeoning output of next-generation sequencing projects performed on a variety of species, and developmental or clinical conditions; the statistical power and complexity of these networks will undoubtedly increase, while their biological relevance will be fiercely challenged. Here, we propose a novel approach to generate a "core" GCN with enhanced biological relevance. Our method involves a data-centering step that effectively removes all primary treatment/tissue effects, which is simple to employ and can be easily implemented into pre-existing GCN analysis pipelines. The gain in biological relevance resulting from the adoption of this approach was assessed using a plant mitochondrial case study.

A panel of correlates predicts vaccine-induced protection of rats against respiratory challenge with virulent Francisella tularensis

There are no defined correlates of protection for any intracellular pathogen, including the bacterium Francisella tularensis, which causes tularemia. Evaluating vaccine efficacy against sporadic diseases like tularemia using field trials is problematic, and therefore alternative strategies to test vaccine candidates like the Francisella Live Vaccine Strain (LVS), such as testing in animals and applying correlate measurements, are needed. Recently, we described a promising correlate strategy that predicted the degree of vaccine-induced protection in mice given parenteral challenges, primarily when using an attenuated Francisella strain. Here, we demonstrate that using peripheral blood lymphocytes (PBLs) in this approach predicts LVS-mediated protection against respiratory challenge of Fischer 344 rats with fully virulent F. tularensis, with exceptional sensitivity and specificity. Rats were vaccinated with a panel of LVS-derived vaccines and subsequently given lethal respiratory challenges with Type A F. tularensis. In parallel, PBLs from vaccinated rats were evaluated for their functional ability to control intramacrophage Francisella growth in in vitro co-culture assays. PBLs recovered from co-cultures were also evaluated for relative gene expression using a large panel of genes identified in murine studies. In vitro control of LVS intramacrophage replication reflected the hierarchy of protection. Further, despite variability between individuals, 22 genes were significantly more up-regulated in PBLs from rats vaccinated with LVS compared to those from rats vaccinated with the variant LVS-R or heat-killed LVS, which were poorly protective. These genes included IFN-γ, IL-21, NOS2, LTA, T-bet, IL-12rβ2, and CCL5. Most importantly, combining quantifications of intramacrophage growth control with 5-7 gene expression levels using multivariate analyses discriminated protected from non-protected individuals with greater than 95% sensitivity and specificity. The results therefore support translation of this approach to non-human primates and people to evaluate new vaccines against Francisella and other intracellular pathogens.

A statistically compiled test battery for feasible evaluation of knee function after rupture of the Anterior Cruciate Ligament - derived from long-term follow-up data

Purpose

Clinical test batteries for evaluation of knee function after injury to the Anterior Cruciate Ligament (ACL) should be valid and feasible, while reliably capturing the outcome of rehabilitation. There is currently a lack of consensus as to which of the many available assessment tools for knee function that should be included. The present aim was to use a statistical approach to investigate the contribution of frequently used tests to avoid redundancy, and filter them down to a proposed comprehensive and yet feasible test battery for long-term evaluation after ACL injury.

Methods

In total 48 outcome variables related to knee function, all potentially relevant for a long-term follow-up, were included from a cross-sectional study where 70 ACL-injured (17–28 years post injury) individuals were compared to 33 controls. Cluster analysis and logistic regression were used to group variables and identify an optimal test battery, from which a summarized estimator of knee function representing various functional aspects was derived.

Results

As expected, several variables were strongly correlated, and the variables also fell into logical clusters with higher within-correlation (max ρ = 0.61) than between clusters (max ρ = 0.19). An extracted test battery with just four variables assessing one-leg balance, isokinetic knee extension strength and hop performance (one-leg hop, side hop) were mathematically combined to an estimator of knee function, which acceptably classified ACL-injured individuals and controls. This estimator, derived from objective measures, correlated significantly with self-reported function, e.g. Lysholm score (ρ = 0.66; p<0.001).

Conclusions

The proposed test battery, based on a solid statistical approach, includes assessments which are all clinically feasible, while also covering complementary aspects of knee function. Similar test batteries could be determined for earlier phases of ACL rehabilitation or to enable longitudinal monitoring. Such developments, established on a well-grounded consensus of measurements, would facilitate comparisons of studies and enable evidence-based rehabilitation.

Models derived from in vitro analyses of spleen, liver, and lung leukocyte functions predict vaccine efficacy against the Francisella tularensis Live Vaccine Strain (LVS)

Currently, there are no licensed vaccines and no correlates of protection against Francisella tularensis, which causes tularemia. We recently demonstrated that measuring in vitro control of intramacrophage bacterial growth by murine F. tularensis-immune splenocytes, as well as transcriptional analyses, discriminated Francisella vaccines of different efficacies. Further, we identified potential correlates of protection against systemic challenge. Here, we extended this approach by studying leukocytes derived from lungs and livers of mice immunized by parenteral and respiratory routes with F. tularensis vaccines. Liver and lung leukocytes derived from intradermally and intranasally vaccinated mice controlled in vitro Francisella Live Vaccine Strain (LVS) intramacrophage replication in patterns similar to those of splenocytes. Gene expression analyses of potential correlates also revealed similar patterns in liver cells and splenocytes. In some cases (e.g., tumor necrosis factor alpha [TNF-α], interleukin 22 [IL-22], and granulocyte-macrophage colony-stimulating factor [GM-CSF]), liver cells exhibited even higher relative gene expression, whereas fewer genes exhibited differential expression in lung cells. In contrast with their strong ability to control LVS replication, splenocytes from intranasally vaccinated mice expressed few genes with a hierarchy of expression similar to that of splenocytes from intradermally vaccinated mice. Thus, the relative levels of gene expression vary between cell types from different organs and by vaccination route. Most importantly, because studies comparing cell sources and routes of vaccination supported the predictive validity of this coculture and gene quantification approach, we combined in vitro LVS replication with gene expression data to develop analytical models that discriminated between vaccine groups and successfully predicted the degree of vaccine efficacy. Thus, this strategy remains a promising means of identifying and quantifying correlative T cell responses. IMPORTANCE Identifying and quantifying correlates of protection is especially challenging for intracellular bacteria, including Francisella tularensis. F. tularensis is classified as a category A bioterrorism agent, and no vaccines have been licensed in the United States, but tularemia is a rare disease. Therefore, clinical trials to test promising vaccines are impractical. In this report, we further evaluated a novel approach to developing correlates by assessing T cell immune responses in lungs and livers of differentially vaccinated mice; these nonprofessional immune tissues are colonized by Francisella. The relative degree of vaccine efficacy against systemic challenge was reflected by the ability of immune T cells, particularly liver T cells, to control the intramacrophage replication of bacteria in vitro and by relative gene expression of several immunological mediators. We therefore developed analytical models that combined bacterial replication data and gene expression data. Several resulting models provided excellent discrimination between vaccines of different efficacies.

Pr is a member of a restricted class of σ70-dependent promoters that lack a recognizable -10 element

The Pr promoter is the first verified member of a class of bacterial σ⁷⁰-promoters that only possess a single match to consensus within its −10 element. In its native context, the activity of this promoter determines the ability of Pseudomonas putida CF600 to degrade phenolic compounds, which provides proof-of-principle for the significance of such promoters. Lack of identity within the −10 element leads to non-detection of Pr-like promoters by current search engines, because of their bias for detection of the −10 motif. Here, we report a mutagenesis analysis of Pr that reveals strict sequence requirements for its activity that includes an essential −15 element and preservation of non-consensus bases within its −35 and −10 elements. We found that highly similar promoters control plasmid- and chromosomally- encoded phenol degradative systems in various Pseudomonads. However, using a purpose-designed promoter-search algorithm and activity analysis of potential candidate promoters, no bona fide Pr-like promoter could be found in the entire genome of P. putida KT2440. Hence, Pr-like σ⁷⁰-promoters, which have the potential to be a widely distributed class of previously unrecognized promoters, are in fact highly restricted and remain in a class of their own.

BALB/c mice, but not C57BL/6 mice immunized with a ΔclpB mutant of Francisella tularensis subspecies tularensis are protected against respiratory challenge with wild-type bacteria: association of protection with post-vaccination and post-challenge immune responses

Francisella tularensis subspecies tularensis is highly virulent for humans especially when it is inhaled. Therefore, it has the potential to be used as a biothreat agent. Vaccines against F. tularensis will need to be approved in accordance with the FDA Animal Rule. This will require identification of robust correlates of protection in experimental animals and the demonstration that similar immune responses are generated in vaccinated humans. Towards this goal, we have developed an experimental live vaccine strain by deleting the gene, clpB, encoding a heat shock protein from virulent subsp. tularensis strain, SCHU S4. SCHU S4ΔclpB administered intradermally protects BALB/c, but not C57BL/6 mice from subsequent respiratory challenge with wildtype SCHU S4. A comparison of post-vaccination and post-challenge immune responses in these two mouse strains shows an association between several antibody and cytokine responses and protection. In particular, elevated IFNγ levels in the skin 2 days after vaccination, sero-conversion to hypothetical membrane protein FTT_1778c, and to 30S ribosomal protein S1 (FTT_0183c) of F. tularensis after 30 days of vaccination, and elevated levels of pulmonary IL-17 on day 7 after respiratory challenge with SCHU S4 were all associated with protection.

Molecular immune responses to aerosol challenge with Francisella tularensis in mice inoculated with live vaccine candidates of varying efficacy

Background

Francisella tularensis is a facultative intracellular bacterial pathogen and the etiological agent of tularemia. The subspecies F. tularensis tularensis is especially virulent for humans when inhaled and respiratory tularemia is associated with high mortality if not promptly treated. A live vaccine strain (LVS) derived from the less virulent holarctica subspecies confers incomplete protection against aerosol challenge with subsp. tularensis. Moreover, correlates of protection have not been established for LVS.

Methodology/Principal Findings

In the present study we compare molecular immune responses elicited by LVS and two defined deletion mutants of clinical subsp. tularensis strain, SCHU S4, that confer enhanced protection in a mouse model. BALB/c mice were immunized intradermally then challenged with an aerosol of SCHU S4 six weeks later. Changes in the levels of a selected panel of cytokines and chemokines were examined in the lungs, spleens, and sera of vaccinated and challenged mice. Mostly, increased cytokine and chemokine levels correlated with increased bacterial burden. However, after adjusting for this variable, immunization with either of the two Schu S4 mutants resulted in higher levels of several pulmonary cytokines, versus those resulting after LVS immunization, including IL-17. Moreover, treatment of mice immunized with ΔclpB with anti-IL-17 antibodies post-challenge enhanced lung infection.

Conclusions/Significance

This is the first report characterizing local and systemic cytokine and chemokine responses in mice immunized with vaccines with different efficacies against aerosol challenge with virulent F. tularensis subsp. tularensis. It shows that increases in the levels of most of these immunomodulators, including those known to be critical for protective immunity, do not superficially correlate with protection unless adjusted for the effects of bacterial burden. Additionally, several cytokines were selectively suppressed in the lungs of naïve mice, suggesting that one mechanism of vaccine action is to overcome this pathogen-induced immunosuppression.

Swelling behavior of the tomato cell wall network

The swelling of tomato pectin and isolated tomato pericarp cell wall material was investigated in aqueous media under different ionic conditions, pH, and external osmotic stress. Conditions were chosen to include those that would be encountered in vivo. Swelling in these systems was strongly influenced by the polyelectrolyte nature of the polymer and the extent of cross-linking with divalent counterions.

Hydration of pectic polysaccharides

The hydration and swelling of pectic polysaccharides was examined at different pHs and ionic strengths as a function of osmotic stress. For weakly charged pectic polysaccharides at low concentrations of a monovalent salt (20 mM), the main driving force for swelling originates from a polyelectrolyte effect due to the translational entropy of ions within the film. Swelling is reduced at higher salt concentrations and lower pHs. Polyelectrolyte collapse and minimal swelling is observed for more highly charged pectic polysaccharides. Replacement of the Na(+) counterion with Ca(2+) results in minimal swelling and the formation of network structures even for the weakly charged pectic polysaccharides.

Hydration of pectic polysaccharides

The hydration and swelling of pectic polysaccharides was examined at different pHs and ionic strengths as a function of osmotic stress. For weakly charged pectic polysaccharides at low concentrations of a monovalent salt (20 mM), the main driving force for swelling originates from a polyelectrolyte effect due to the translational entropy of ions within the film. Swelling is reduced at higher salt concentrations and lower pHs. Polyelectrolyte collapse and minimal swelling is observed for more highly charged pectic polysaccharides. Replacement of the Na(+) counterion with Ca(2+) results in minimal swelling and the formation of network structures even for the weakly charged pectic polysaccharides.

13C MAS NMR studies of the effects of hydration on the cell walls of potatoes and Chinese water chestnuts

13C NMR with magic angle spinning (MAS) has been employed to investigate the cell walls of potatoes and Chinese water chestnuts over a range of hydration levels. Both single-pulse excitation (SPEMAS) and cross-polarization (CPMAS) experiments were carried out. Hydration led to a substantial increase in signal intensities of galactan and galacturonan in the SPEMAS spectra and a decrease in line width, implying mobilization in the backbone and side chains of pectin. In CPMAS spectra of both samples, noncellulose components showed signal loss as hydration increased. However, the signals of some galacturonan in the 3(1) helix configuration remained in the spectra even when the water content was as high as 110%. Cellulose was unaffected. It is concluded that the pectic polysaccharides experience a distribution of molecular conformations and mobility, whereas cellulose remained as typical rigid solid.

Ileal effluent as a fermentation substrate: implications for butyrate production in the colon

Fermentation of fiber can lead to an enhanced production of short chain fatty acids (SCFA) and, hence, contribute to the proposed anticarcinogenic properties of butyrate in the colon. The fermentation of fiber isolates and the corresponding ileal effluents has been compared under in vitro conditions. Yield of SCFA per gram of substrate fermented was similar for isolates and fiber-enriched effluents (approximately 4.9 mmol/g) and it could be inferred that nonfiber components of effluent also generated SCFA. Butyrate production was highest for glucan-based polymers (approximately 30% total SCFA) and, from the measured acidogenic profile, production of SCFA will occur mainly in the proximal colon. The buffering capacity of ileal effluents during fermentation restrict the potential for a reduction in pH during acidogenesis compared to fiber isolates. This buffering capacity could limit the bioavailability of butyrate in the colon and, hence, the ability to satisfy the proposed antineoplastic properties of butyrate in the colon.

Characterisation of the binding capacities and affinities of wheat bran, fruit and vegetable fibres for MeIQx, before and after fermentation

Binding capacities and affinities of dietary fibres and fermented fibre residues for MeIQx are much higher for brans than for fruit and vegetable fibres. Bran residues concentrated distally in the colon would provide a substantial binding matrix. Small intestinal conditions are not so conducive to hydrophobic adsorption.

The effects of pH and bile salts on the binding of MeIQx to wheat bran fibre

The effects of pH and of the bile salts, sodium cholate, chenodeoxycholate, taurocholate, deoxycholate and glycodeoxycholate on the adsorption of the heterocyclic aromatic amine mutagen and carcinogen 2-amino-3,8-dimethylimidazo[4,5-f]-quinoxaline (MeIQx) onto wheat bran cell wall material have been examined. The highest binding affinity of MeIQx for bran was at pH 5.5. Binding affinity declined more rapidly with pH > 5.5 than with pH < 5.5. Bile salts in solution did not appreciably affect the binding of MeIQx to bran, but where the bile salts formed a suspension then the adsorption was reduced. Where precipitation of bile salts occurred then the free concentration of MeIQx was also reduced, indicating that MeIQx binding/immobilization could be enhanced through interaction with bile salts.

The consequences of fruit and vegetable fibre fermentation on their binding capacity for MeIQx and the effects of soluble fibre sources on the binding affinity of wheat bran preparations

Fruits and vegetables provide dietary fibre some of which is partly soluble in the upper gut; in the colon it is highly fermentable. Using alcohol-insoluble residues prepared from a range of fruits and vegetables the effects of fermentation on the changes in composition and binding capacity have been assessed for the hydrophobic mutagen 2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline (MeIQx). Fermentation was extensive and resulted in destruction of most of the pectic polysaccharides. Of the unfermented vegetable fibre only cabbage had a measurable hydrophobic binding capacity. The binding capacities of unfermented apple, carrot and sugar beet were negligible. After fermentation, binding capacities, (per mg of fermented residue), increased. Although fermented cabbage was found to have the highest capacity of the fruit and vegetable fibres this remained less than the least effective of unfermented wheat bran samples which had a relatively high affinity for MeIQx. Mucin inhibited the binding of MeIQx to wheat bran fibre but apple fibre did not. The results show that the contribution of fruit and vegetable fibre to a hydrophobic binding matrix in the colon is insignificant and the suggested harmful effect of fruit and vegetable fibre, maintaining hydrophobic mutagens in solution, can be prevented by the presence of wheat bran fibre.

The effect of fibre source and fermentation on the apparent hydrophobic binding properties of wheat bran preparations for the mutagen 2-amino-3,8- dimethylimidazo[4,5-f]quinoxaline (MeIQx)

Preparations of dietary fibre have an ability to bind potentially toxic compounds in vitro, but it is unclear how the measured binding properties are a reflection of either the relative affinity of such compounds for fibre or the saturation binding capacity of the fibre. Cooking and processing of foods and fermentation activity in the colon can result in significant modification of the structure of the fibre matrix. Hence, binding properties measured in vitro may be significantly altered from those at the proposed site of fibre activity. Using cell wall material (CWM) prepared from three distinct wheat fibre sources, the effects of fermentation have been monitored and the consequences on the binding properties assessed. The CWM preparations were feremented for either 0, 6, 18 or 24 h in vitro, using a human faecal inoculum. The coarse bran CWM was also subjected to a simulated gastric treatment. Fermentation resulted in a loss of material from each CWM preparation, the loss being consistent with the bran source, and the maximum extent of fermentation was reached within 18 h. The in vitro binding of the mutagen 2-amino-3,8- dimethylimidazo[4,5-f]quinoxaline (MeIQx) to the unfermented and fermented fibre fractions showed fine bran CWM had a much higher affinity for MeIQx than coarse bran. A 6 h fermentation had little effect on the binding affinity to fine bran, but the affinity for coarse bran became similar to that of the fine bran. Further fermentation for 18 h or 24 h resulted in a further slight increase in the binding affinity for coarse and fine bran, but, more significantly, the binding capacity of each was increased. The binding properties of the beeswing bran were not significantly affected by fermentation. Wheat bran has the potential to bind hydrophobic mutagens in the diet and this potential can be enhanced after fermentation under colonic conditions.

Cell-wall polysaccharides and glycoproteins of parenchymatous tissues of runner bean (Phaseolus coccineus)

1. Polymers were solubilized from the cell walls of parenchyma from mature runner-bean pods with minimum degradation by successive extractions with cyclohexane-trans-1,2-diamine-NNN'N'-tetra-acetate (CDTA), Na2CO3 and KOH to leave the alpha-cellulose residue, which contained cross-linked pectic polysaccharides and Hyp-rich glycoproteins. These were solubilized with chlorite/acetic acid and cellulase. The polymers were fractionated by anion-exchange chromatography, and fractions were subjected to methylation analysis. 2. The pectic polysaccharides differed in their ease of extraction, and a small proportion were highly cross-linked. The bulk of the pectic polysaccharides solubilized by CDTA and Na2CO3 were less branched than those solubilized by KOH. There was good evidence that most of the pectic polysaccharides were not degraded during extraction. 3. The protein-containing fractions included Hyp-rich and Hyp-poor glycoproteins associated with easily extractable pectic polysaccharides, Hyp-rich glycoproteins solubilized with 4M-KOH+borate, the bulk of which were not associated with pectic polysaccharides, and highly cross-linked Hyp-rich glycoproteins. 4. Isodityrosine was not detected, suggesting that it does not have a (major) cross-linking role in these walls. Instead, it is suggested that phenolics, presumably linked to C-5 of 3,5-linked Araf residues of Hyp-rich glycoproteins, serve to cross-link some of the polymers. 5. There were two main types of xyloglucan, with different degrees of branching. The bulk of the less branched xyloglucans were solubilized by more-concentrated alkali. The anomeric configurations of the sugars in one of the highly branched xyloglucans were determined by 13C-n.m.r. spectroscopy. 6. The structural features of the cell-wall polymers and complexes are discussed in relation to the structure of the cell walls of parenchyma tissues.