Experimental evolution links post-transcriptional regulation to Leishmania fitness gain.

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

Constrained Fourier estimation of short-term time-series gene expression data reduces noise and improves clustering and gene regulatory network predictions.

BACKGROUND: Biological data suffers from noise that is inherent in the measurements. This is particularly true for time-series gene expression measurements. Nevertheless, in order to to explore cellular dynamics, scientists employ such noisy measurements in predictive and clustering tools. However, noisy data can not only obscure the genes temporal patterns, but applying predictive and clustering tools on noisy data may yield inconsistent, and potentially incorrect, results. RESULTS: To reduce the noise of short-term (< 48 h) time-series expression data, we relied on the three basic temporal patterns of gene expression: waves, impulses and sustained responses. We constrained the estimation of the true signals to these patterns by estimating the parameters of first and second-order Fourier functions and using the nonlinear least-squares trust-region optimization technique. Our approach lowered the noise in at least 85% of synthetic time-series expression data, significantly more than the spline method ([Formula: see text]). When the data contained a higher signal-to-noise ratio, our method allowed downstream network component analyses to calculate consistent and accurate predictions, particularly when the noise variance was high. Conversely, these tools led to erroneous results from untreated noisy data. Our results suggest that at least 5-7 time points are required to efficiently de-noise logarithmic scaled time-series expression data. Investing in sampling additional time points provides little benefit to clustering and prediction accuracy. CONCLUSIONS: Our constrained Fourier de-noising method helps to cluster noisy gene expression and interpret dynamic gene networks more accurately. The benefit of noise reduction is large and can constitute the difference between a successful application and a failing one.

Recombinant production of the lantibiotic nisin using Corynebacterium glutamicum in a two-step process.

BACKGROUND: The bacteriocin nisin is naturally produced by Lactococcus lactis as an inactive prepeptide that is modified posttranslationally resulting in five (methyl-)lanthionine rings characteristic for class Ia bacteriocins. Export and proteolytic cleavage of the leader peptide results in release of active nisin. By targeting the universal peptidoglycan precursor lipid II, nisin has a broad target spectrum including important human pathogens such as Listeria monocytogenes and methicillin-resistant Staphylococcus aureus strains. Industrial nisin production is currently performed using natural producer strains resulting in rather low product purity and limiting its application to preservation of dairy food products. RESULTS: We established heterologous nisin production using the biotechnological workhorse organism Corynebacterium glutamicum in a two-step process. We demonstrate successful biosynthesis and export of fully modified prenisin and its activation to mature nisin by a purified, soluble variant of the nisin protease NisP (sNisP) produced in Escherichia coli. Active nisin was detected by a L. lactis sensor strain with strictly nisin-dependent expression of the fluorescent protein mCherry. Following activation by sNisP, supernatants of the recombinant C. glutamicum producer strain cultivated in standard batch fermentations contained at least 1.25 mg/l active nisin. CONCLUSIONS: We demonstrate successful implementation of a two-step process for recombinant production of active nisin with C. glutamicum. This extends the spectrum of bioactive compounds that may be produced using C. glutamicum to a bacteriocin harboring complex posttranslational modifications. Our results provide a basis for further studies to optimize product yields, transfer production to sustainable substrates and purification of pharmaceutical grade nisin.

Analysis of temporal gene regulation of Listeria monocytogenes revealed distinct regulatory response modes after exposure to high pressure processing.

BACKGROUND: The pathogen Listeria (L.) monocytogenes is known to survive heat, cold, high pressure, and other extreme conditions. Although the response of this pathogen to pH, osmotic, temperature, and oxidative stress has been studied extensively, its reaction to the stress produced by high pressure processing HPP (which is a preservation method in the food industry), and the activated gene regulatory network (GRN) in response to this stress is still largely unknown. RESULTS: We used RNA sequencing transcriptome data of L. monocytogenes (ScottA) treated at 400 MPa and 8∘C, for 8 min and combined it with current information in the literature to create a transcriptional regulation database, depicting the relationship between transcription factors (TFs) and their target genes (TGs) in L. monocytogenes. We then applied network component analysis (NCA), a matrix decomposition method, to reconstruct the activities of the TFs over time. According to our findings, L. monocytogenes responded to the stress applied during HPP by three statistically different gene regulation modes: survival mode during the first 10 min post-treatment, repair mode during 1 h post-treatment, and re-growth mode beyond 6 h after HPP. We identified the TFs and their TGs that were responsible for each of the modes. We developed a plausible model that could explain the regulatory mechanism that L. monocytogenes activated through the well-studied CIRCE operon via the regulator HrcA during the survival mode. CONCLUSIONS: Our findings suggest that the timely activation of TFs associated with an immediate stress response, followed by the expression of genes for repair purposes, and then re-growth and metabolism, could be a strategy of L. monocytogenes to survive and recover extreme HPP conditions. We believe that our results give a better understanding of L. monocytogenes behavior after exposure to high pressure that may lead to the design of a specific knock-out process to target the genes or mechanisms. The results can help the food industry select appropriate HPP conditions to prevent L. monocytogenes recovery during food storage.

High-pressure processing-induced transcriptome response during recovery of Listeria monocytogenes.

BACKGROUND: High-pressure processing (HPP) is a commonly used technique in the food industry to inactivate pathogens, including L. monocytogenes. It has been shown that L. monocytogenes is able to recover from HPP injuries and can start to grow again during long-term cold storage. To date, the gene expression profiling of L. monocytogenes during HPP damage recovery at cooling temperature has not been studied. In order identify key genes that play a role in recovery of the damage caused by HPP treatment, we performed RNA-sequencing (RNA-seq) for two L. monocytogenes strains (barotolerant RO15 and barosensitive ScottA) at nine selected time points (up to 48 h) after treatment with two pressure levels (200 and 400 MPa). RESULTS: The results showed that a general stress response was activated by SigB after HPP treatment. In addition, the phosphotransferase system (PTS; mostly fructose-, mannose-, galactitol-, cellobiose-, and ascorbate-specific PTS systems), protein folding, and cobalamin biosynthesis were the most upregulated genes during HPP damage recovery. We observed that cell-division-related genes (divIC, dicIVA, ftsE, and ftsX) were downregulated. By contrast, peptidoglycan-synthesis genes (murG, murC, and pbp2A) were upregulated. This indicates that cell-wall repair occurs as a part of HPP damage recovery. We also observed that prophage genes, including anti-CRISPR genes, were induced by HPP. Interestingly, a large amount of RNA-seq data (up to 85%) was mapped to Rli47, which is a non-coding RNA that is upregulated after HPP. Thus, we predicted that Rli47 plays a role in HPP damage recovery in L. monocytogenes. Moreover, gene-deletion experiments showed that amongst peptidoglycan biosynthesis genes, pbp2A mutants are more sensitive to HPP. CONCLUSIONS: We identified several genes and mechanisms that may play a role in recovery from HPP damage of L. monocytogenes. Our study contributes to new information on pathogen inactivation by HPP.

Establishing recombinant production of pediocin PA-1 in Corynebacterium glutamicum.

Bacteriocins are antimicrobial peptides produced by bacteria to inhibit competitors in their natural environments. Some of these peptides have emerged as commercial food preservatives and, due to the rapid increase in antibiotic resistant bacteria, are also discussed as interesting alternatives to antibiotics for therapeutic purposes. Currently, commercial bacteriocins are produced exclusively with natural producer organisms on complex substrates and are sold as semi-purified preparations or crude fermentates. To allow clinical application, efficacy of production and purity of the product need to be improved. This can be achieved by shifting production to recombinant microorganisms. Here, we identify Corynebacterium glutamicum as a suitable production host for the bacteriocin pediocin PA-1. C. glutamicum CR099 shows resistance to high concentrations of pediocin PA-1 and the bacteriocin was not inactivated when spiked into growing cultures of this bacterium. Recombinant C. glutamicum expressing a synthetic pedACDCgl operon releases a compound that has potent antimicrobial activity against Listeria monocytogenes and Listeria innocua and matches size and mass:charge ratio of commercial pediocin PA-1. Fermentations in shake flasks and bioreactors suggest that low levels of dissolved oxygen are favorable for production of pediocin. Under these conditions, however, reduced activity of the TCA cycle resulted in decreased availability of the important pediocin precursor l-asparagine suggesting options for further improvement. Overall, we demonstrate that C. glutamicum is a suitable host for recombinant production of bacteriocins of the pediocin family.

Factors associated with parental grief reaction following pediatric acquired brain injury.

The present study investigated factors associated with parental grief reaction (PGR) following pediatric acquired brain injury (ABI), and compared PGR to the one exhibited following child death. Fifty-seven parents of 51 children (aged 3-18) whose ABI occurred 1-14 years before participation, completed the multi-scale Two-Track Bereavement Questionnaire; a socio-demographic questionnaire; and a scale assessing perceived behavioural changes in the child. Results from regression analysis indicated that time since injury had no impact on parents' grief other than having an adverse impact on their overall coping and functioning; A higher amount of weekly caring hours predicted only a greater traumatic perception of the loss; Older children's ages but mostly greater parental-perceived behavioural changes, predicted greater PGR on most scales. PGR was compared with the pre-existing data of bereaved parents who completed the same grief questionnaire. Although grief response patterns and intensity were similar in both groups, significant differences were found on scales assessing the continuing bond with the child: relational active grief, close and positive relationship, and conflictual relationship. Our findings indicate that parental grief is multi-dimensional following pediatric ABI and illuminate the interplay between elements characterizing parents' nonfinite vs. finite loss experience.

Genomic characterization of the most barotolerant Listeria monocytogenes RO15 strain compared to reference strains used to evaluate food high pressure processing.

BACKGROUND: High pressure processing (HPP; i.e. 100-600 MPa pressure depending on product) is a non-thermal preservation technique adopted by the food industry to decrease significantly foodborne pathogens, including Listeria monocytogenes, from food. However, susceptibility towards pressure differs among diverse strains of L. monocytogenes and it is unclear if this is due to their intrinsic characteristics related to genomic content. Here, we tested the barotolerance of 10 different L. monocytogenes strains, from food and food processing environments and widely used reference strains including clinical isolate, to pressure treatments with 400 and 600 MPa. Genome sequencing and genome comparison of the tested L. monocytogenes strains were performed to investigate the relation between genomic profile and pressure tolerance. RESULTS: None of the tested strains were tolerant to 600 MPa. A reduction of more than 5 log10 was observed for all strains after 1 min 600 MPa pressure treatment. L. monocytogenes strain RO15 showed no significant reduction in viable cell counts after 400 MPa for 1 min and was therefore defined as barotolerant. Genome analysis of so far unsequenced L. monocytogenes strain RO15, 2HF33, MB5, AB199, AB120, C7, and RO4 allowed us to compare the gene content of all strains tested. This revealed that the three most pressure tolerant strains had more than one CRISPR system with self-targeting spacers. Furthermore, several anti-CRISPR genes were detected in these strains. Pan-genome analysis showed that 10 prophage genes were significantly associated with the three most barotolerant strains. CONCLUSIONS: L. monocytogenes strain RO15 was the most pressure tolerant among the selected strains. Genome comparison suggests that there might be a relationship between prophages and pressure tolerance in L. monocytogenes.

ILDR2 Is a Novel B7-like Protein That Negatively Regulates T Cell Responses.

The B7-like protein family members play critical immunomodulatory roles and constitute attractive targets for the development of novel therapies for human diseases. We identified Ig-like domain-containing receptor (ILDR)2 as a novel B7-like protein with robust T cell inhibitory activity, expressed in immune cells and in immune-privileged and inflamed tissues. A fusion protein, consisting of ILDR2 extracellular domain with an Fc fragment, that binds to a putative counterpart on activated T cells showed a beneficial effect in the collagen-induced arthritis model and abrogated the production of proinflammatory cytokines and chemokines in autologous synovial-like cocultures of macrophages and cytokine-stimulated T cells. Collectively, these findings point to ILDR2 as a novel negative regulator for T cells, with potential roles in the development of immune-related diseases, including autoimmunity and cancer.

Locked-in your heart-shaped box: Familial-role and attachment orientation as predictors of grief in prolonged disorders of consciousness vs. death.

Death or prolonged disorders of consciousness (DOC) of a loved one are both considered relational-losses that severely disrupt attachment-bonds. Grief in both conditions was compared by exploring the impact of familial-role and attachment-orientation. In DOC, caregivers' grief was found significantly intensified relative to Death. Familial-role impacted grief in both conditions alike, with partners' heightened grief in DOC reflecting the complexity of their stagnant bonds. In Death, avoidance-attachment mitigated grief, while in DOC anxiety-attachment accentuated grief, we suggest that while physical-separation in death facilitates the modification of continuing attachment-schema, in DOC, modification may be required while the patient is still alive.

A sensory-motor control model of animal flight explains why bats fly differently in light versus dark.

Animal flight requires fine motor control. However, it is unknown how flying animals rapidly transform noisy sensory information into adequate motor commands. Here we developed a sensorimotor control model that explains vertebrate flight guidance with high fidelity. This simple model accurately reconstructed complex trajectories of bats flying in the dark. The model implies that in order to apply appropriate motor commands, bats have to estimate not only the angle-to-target, as was previously assumed, but also the angular velocity ("proportional-derivative" controller). Next, we conducted experiments in which bats flew in light conditions. When using vision, bats altered their movements, reducing the flight curvature. This change was explained by the model via reduction in sensory noise under vision versus pure echolocation. These results imply a surprising link between sensory noise and movement dynamics. We propose that this sensory-motor link is fundamental to motion control in rapidly moving animals under different sensory conditions, on land, sea, or air.

Landscape of transcription in human cells.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

DTIE, a novel core promoter element that directs start site selection in TATA-less genes.

The transcription start site (TSS) determines the length and composition of the 5' UTR and therefore can have a profound effect on translation. Yet, little is known about the mechanism underlying start site selection, particularly from promoters lacking conventional core elements such as TATA-box and Initiator. Here we report a novel mechanism of start site selection in the TATA- and Initiator-less promoter of miR-22, through a strictly localized downstream element termed DTIE and an upstream distal element. Changing the distance between them reduced promoter strength, altered TSS selection and diminished Pol II recruitment. Biochemical assays suggest that DTIE does not serve as a docking site for TFIID, the major core promoter-binding factor. TFIID is recruited to the promoter through DTIE but is dispensable for TSS selection. We determined DTIE consensus and found it to be remarkably prevalent, present at the same TSS downstream location in ≈20.8% of human promoters, the vast majority of which are TATA-less. Analysis of DTIE in the tumor suppressor p53 confirmed a similar function. Our findings reveal a novel mechanism of transcription initiation from TATA-less promoters.

Iterative sub-network component analysis enables reconstruction of large scale genetic networks.

BACKGROUND: Network component analysis (NCA) became a popular tool to understand complex regulatory networks. The method uses high-throughput gene expression data and a priori topology to reconstruct transcription factor activity profiles. Current NCA algorithms are constrained by several conditions posed on the network topology, to guarantee unique reconstruction (termed compliancy). However, the restrictions these conditions pose are not necessarily true from biological perspective and they force network size reduction, pruning potentially important components. RESULTS: To address this, we developed a novel, Iterative Sub-Network Component Analysis (ISNCA) for reconstructing networks at any size. By dividing the initial network into smaller, compliant subnetworks, the algorithm first predicts the reconstruction of each subnetwork using standard NCA algorithms. It then subtracts from the reconstruction the contribution of the shared components from the other subnetwork. We tested the ISNCA on real, large datasets using various NCA algorithms. The size of the networks we tested and the accuracy of the reconstruction increased significantly. Importantly, FOXA1, ATF2, ATF3 and many other known key regulators in breast cancer could not be incorporated by any NCA algorithm because of the necessary conditions. However, their temporal activities could be reconstructed by our algorithm, and therefore their involvement in breast cancer could be analyzed. CONCLUSIONS: Our framework enables reconstruction of large gene expression data networks, without reducing their size or pruning potentially important components, and at the same time rendering the results more biological plausible. Our ISNCA method is not only suitable for prediction of key regulators in cancer studies, but it can be applied to any high-throughput gene expression data.

Sequence-Dependent Promoter Escape Efficiency Is Strongly Influenced by Bias for the Pretranslocated State during Initial Transcription.

Abortive transcription initiation can be rate-limiting for promoter escape and therefore represents a barrier to productive gene expression. The mechanism for abortive initiation is unknown, but the amount of abortive transcript is known to vary with the composition of the initial transcribed sequence (ITS). Here, we used a thermodynamic model of translocation combined with experimental validation to investigate the relationship between ITS and promoter escape on a set of phage T5 N25 promoters. We found a strong, negative correlation between RNAP's propensity to occupy the pretranslocated state during initial transcription and the efficiency of promoter escape (r = -0.67; p < 10(-6)). This correlation was almost entirely caused by free energy changes due to variation in the RNA 3' dinucleotide sequence at each step, implying that this sequence element controls the disposition of initial transcribing complexes. We tested our model experimentally by constructing a set of novel N25-ITS promoter variants; quantitative transcription analysis again showed a strong correlation (r = -0.81; p < 10(-6)). Our results support a model in which sequence-directed bias for the pretranslocated state during scrunching results in increased backtracking, which limits the efficiency of promoter escape. This provides an answer to the long-standing issue of how sequence composition of the ITS affects promoter escape efficiency.

Reconstruction of temporal activity of microRNAs from gene expression data in breast cancer cell line.

BACKGROUND: MicroRNAs (miRNAs) are small non-coding RNAs that regulate genes at the post-transcriptional level in spatiotemporal manner. Several miRNAs are identified as prognostic and diagnostic markers in many human cancers. Estimation of the temporal activities of the miRNAs is an important step in the way to understand the complex interactions of these important regulatory elements with transcription factors (TFs) and target genes (TGs). However, current research on miRNA activities excludes network dynamics from the studies, disregarding the important element of time in the regulatory network analysis. RESULTS: In the current study, we combined experimentally verified miRNA-TG interactions with breast cancer microarray TG expression data to identify key miRNAs and compute their temporal activity using network component analysis (NCA). The computed activities showed that miRNAs were regulated in a time dependent manner. Our results allowed constructing a synergistic network of miRNAs using the computed miRNA activities and their shared regulation of TGs. We further extended this network by incorporating miRNA-TG, miRNA-TF, TF-miRNA and TF-TG regulations in the context of breast cancer. Our integrated network identified several miRNAs known to be involved in breast cancer regulation and revealed several novel miRNAs. Our further analysis detected substantial involvement of the miRNAs miR-324, miR-93, miR-615 and miR-1 in breast cancer, which was not known previously. Next, combining our integrated networks with functional annotation of differentially expressed genes resulted in new sub-networks. These sub-networks allowed us to identify the key miRNAs and their interactions with TFs and TGs of several biological processes involved in breast cancer. The identified markers are validated for their potential as prognostic markers for breast cancer through survival analysis. CONCLUSIONS: Our dynamical analysis of the miRNA interactions greatly helps to discover new network based markers, and is highly applicable (but not limited) to cancer research.

Robust, fully quantifiable and scalable bioprocess utilizing spent sulfite liquor with Corynebacterium glutamicum.

In this study, a bioprocessing strategy was designed to valorize ultra-filtered spent sulfite liquor (UF-SSL) without prior detoxification steps as well as using it purely as a carbon source supplement to defined or complex media. Hence, a minimal medium for the bioconversion of UF-SSL with Corynebacterium glutamicum was developed and process robustness and reproducibility were validated. Process quantifiability was ensured by development of a biomass measurement technique for matrices with high water-insoluble solids and verified using elemental balancing. Mechanistic modeling based on Monod equations was used to identify batch kinetics. In a final step, scale-up of the developed process was performed to showcase process maturity towards commercialisation.

Genome-assisted Identification, Purification, and Characterization of Bacteriocins.

Bacteriocins are antimicrobial peptides with activity against antibiotic resistant bacterial pathogens. Here, we describe a set of methods aimed at purifying, identifying, and characterizing new bacteriocins. The purification consists of ammonium sulphate precipitation, cation-exchange chromatography, and reversed-phase chromatography. The yield of the bacteriocin is quantified by bacteriocin antimicrobial activity in a microtiter plate assay after each purification step. The mass of the purified bacteriocin is assessed by MALDI TOF MS analysis of the active fractions after reversed-phase chromatography. The mass is compared with the theoretical mass based on genetic information from the whole genome sequencing of the bacteriocin producer strain. Physicochemical characterization is performed by assessing antimicrobial activity following heat and protease treatments. Fluorescent techniques are used to examine the capacity of the bacteriocin to disrupt membrane integrity. Herein a set of protocols for purification and characterization of the bacteriocin nisin Z is used as a typical example in this paper.

Novel computational analysis of large transcriptome datasets identifies sets of genes distinguishing chronic obstructive pulmonary disease from healthy lung samples.

Chronic obstructive pulmonary disease (COPD) kills over three million people worldwide every year. Despite its high global impact, the knowledge about the underlying molecular mechanisms is still limited. In this study, we aimed to extend the available knowledge by identifying a small set of COPD-associated genes. We analysed different publicly available gene expression datasets containing whole lung tissue (WLT) and airway epithelium (AE) samples from over 400 human subjects for differentially expressed genes (DEGs). We reduced the resulting sets of 436 and 663 DEGs using a novel computational approach that utilises a random depth-first search to identify genes which improve the distinction between COPD patients and controls along the first principle component of the data. Our method identified small sets of 10 and 15 genes in the WLT and AE, respectively. These sets of genes significantly (p < 10-20) distinguish COPD patients from controls with high fidelity. The final sets revealed novel genes like cysteine rich protein 1 (CRIP1) or secretoglobin family 3A member 2 (SCGB3A2) that may underlie fundamental molecular mechanisms of COPD in these tissues.