FunPat: function-based pattern analysis on RNA-seq time series data

Mechanism of Chinese botanical drug Dizhi pill for myopia: An integrated study based on bioinformatics and network analysis

To identify the active constituents, core targets, immunomodulatory functions and potential mechanisms of Dizhi pill (DZP) in the treatment of myopia. The active constituents and drug targets of DZP were searched in the TCMSP, Herb databases and correlational studies. The targets of myopia were searched in the TTD, Genecards, OMIM and Drugbank databases. Gene expression profile data of GSE136701 were downloaded from the GEO database and subjected to WGCNA and DEG analysis to screen for significant modules and targets of myopia. Intersectional targets of myopia and DZP and core targets of myopia were analyzed through the String database. The GO and KEGG enrichment analyses of the interested targets were conducted. Cibersort algorithm was used for immune infiltration analysis to investigate the immunomodulatory functions of DZP on myopia. Autodock was used to dock the important targets and active constituents. Eight targets (STAT3, PIK3CA, PIK3R1, MAPK1, MAPK3, HSP90AA1, MIP, and LGSN) and 5 active constituents (Quercetin, Beta-sitosterol, Diincarvilone A, Ferulic acid methyl ester, and Naringenin) were identified from DZP. In pathways identified by the GO and KEGG enrichment analyses, "ATP metabolic process" and "AGE-RAGE diabetes complication signaling" pathways were closely related to the mechanisms of DZP in the treatment of myopia. Molecular docking showed that both the intersectional targets and core targets of myopia could bind stably and spontaneously with the active constituents of DZP. This study suggested that the mechanisms of DZP in the treatment of myopia were related to active constituents: Quercetin, Beta-sitosterol, Diincarvilone A, Ferulic acid methyl ester and Naringenin, intersectional targets: STAT3, PIK3CA, PIK3R1, MAPK1, MAPK3, and HSP90AA1, core targets of myopia: MIP and LGSN, AGE-RAGE signaling pathway, positive regulation of ATP metabolic process pathway and immunomodulatory functions.

Learning perturbation-inducible cell states from observability analysis of transcriptome dynamics

A major challenge in biotechnology and biomanufacturing is the identification of a set of biomarkers for perturbations and metabolites of interest. Here, we develop a data-driven, transcriptome-wide approach to rank perturbation-inducible genes from time-series RNA sequencing data for the discovery of analyte-responsive promoters. This provides a set of biomarkers that act as a proxy for the transcriptional state referred to as cell state. We construct low-dimensional models of gene expression dynamics and rank genes by their ability to capture the perturbation-specific cell state using a novel observability analysis. Using this ranking, we extract 15 analyte-responsive promoters for the organophosphate malathion in the underutilized host organism Pseudomonas fluorescens SBW25. We develop synthetic genetic reporters from each analyte-responsive promoter and characterize their response to malathion. Furthermore, we enhance malathion reporting through the aggregation of the response of individual reporters with a synthetic consortium approach, and we exemplify the library's ability to be useful outside the lab by detecting malathion in the environment. The engineered host cell, a living malathion sensor, can be optimized for use in environmental diagnostics while the developed machine learning tool can be applied to discover perturbation-inducible gene expression systems in the compendium of host organisms.

Transcript Assembly and Annotations: Bias and Adjustment

Motivation

Transcript annotations play a critical role in gene expression analysis as they serve as a reference for quantifying isoform-level expression. The two main sources of annotations are RefSeq and Ensembl/GENCODE, but discrepancies between their methodologies and information resources can lead to significant differences. It has been demonstrated that the choice of annotation can have a significant impact on gene expression analysis. Furthermore, transcript assembly is closely linked to annotations, as assembling large-scale available RNA-seq data is an effective data-driven way to construct annotations, and annotations are often served as benchmarks to evaluate the accuracy of assembly methods. However, the influence of different annotations on transcript assembly is not yet fully understood.

Results

We investigate the impact of annotations on transcript assembly. We observe that conflicting conclusions can arise when evaluating assemblers with different annotations. To understand this striking phenomenon, we compare the structural similarity of annotations at various levels and find that the primary structural difference across annotations occurs at the intron-chain level. Next, we examine the biotypes of annotated and assembled transcripts and uncover a significant bias towards annotating and assembling transcripts with intron retentions, which explains above the contradictory conclusions. We develop a standalone tool, available at https://github.com/Shao-Group/irtool, that can be combined with an assembler to generate an assembly without intron retentions. We evaluate the performance of such a pipeline and offer guidance to select appropriate assembling tools for different application scenarios.

Role of the Extracytoplasmic Function Sigma Factor SigE in the Stringent Response of Mycobacterium tuberculosis

Bacteria respond to nutrient starvation implementing the stringent response, a stress signaling system resulting in metabolic remodeling leading to decreased growth rate and energy requirements. A well-characterized model of stringent response in Mycobacterium tuberculosis is the one induced by growth in low phosphate. The extracytoplasmic function (ECF) sigma factor SigE was previously suggested as having a key role in the activation of stringent response. In this study, we challenge this hypothesis by analyzing the temporal dynamics of the transcriptional response of a sigE mutant and its wild-type parental strain to low phosphate using RNA sequencing. We found that both strains responded to low phosphate with a typical stringent response trait, including the downregulation of genes encoding ribosomal proteins and RNA polymerase. We also observed transcriptional changes that support the occurring of an energetics imbalance, compensated by a reduced activity of the electron transport chain, decreased export of protons, and a remodeling of central metabolism. The most striking difference between the two strains was the induction in the sigE mutant of several stress-related genes, in particular, the genes encoding the ECF sigma factor SigH and the transcriptional regulator WhiB6. Since both proteins respond to redox unbalances, their induction suggests that the sigE mutant is not able to maintain redox homeostasis in response to the energetics imbalance induced by low phosphate. In conclusion, our data suggest that SigE is not directly involved in initiating stringent response but in protecting the cell from stress consequent to the low phosphate exposure and activation of stringent response. IMPORTANCE Mycobacterium tuberculosis can enter a dormant state enabling it to establish latent infections and to become tolerant to antibacterial drugs. Dormant bacteria's physiology and the mechanism(s) used by bacteria to enter dormancy during infection are still unknown due to the lack of reliable animal models. However, several in vitro models, mimicking conditions encountered during infection, can reproduce different aspects of dormancy (growth arrest, metabolic slowdown, drug tolerance). The stringent response, a stress response program enabling bacteria to cope with nutrient starvation, is one of them. In this study, we provide evidence suggesting that the sigma factor SigE is not directly involved in the activation of stringent response as previously hypothesized, but it is important to help the bacteria to handle the metabolic stress related to the adaptation to low phosphate and activation of stringent response, thus giving an important contribution to our understanding of the mechanism behind stringent response development.

Advanced variant classification framework reduces the false positive rate of predicted loss of function (pLoF) variants in population sequencing data

Predicted loss of function (pLoF) variants are highly deleterious and play an important role in disease biology, but many of these variants may not actually result in loss-of-function. Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines's PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 autosomal recessive disease-genes from the Genome Aggregation Database (gnomAD, v2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in low per-base expression (pext) score regions, or the presence of cryptic splice rescues. Variants predicted to be potential artifacts or to evade LoF were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of LoF evading variants assessed, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines, and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

The Pancreatic ß-cell Response to Secretory Demands and Adaption to Stress

Pancreatic β cells dedicate much of their protein translation capacity to producing insulin to maintain glucose homeostasis. In response to increased secretory demand, β cells can compensate by increasing insulin production capability even in the face of protracted peripheral insulin resistance. The ability to amplify insulin secretion in response to hyperglycemia is a critical facet of β-cell function, and the exact mechanisms by which this occurs have been studied for decades. To adapt to the constant and fast-changing demands for insulin production, β cells use the unfolded protein response of the endoplasmic reticulum. Failure of these compensatory mechanisms contributes to both type 1 and 2 diabetes. Additionally, studies in which β cells are "rested" by reducing endogenous insulin demand have shown promise as a therapeutic strategy that could be applied more broadly. Here, we review recent findings in β cells pertaining to the metabolic amplifying pathway, the unfolded protein response, and potential advances in therapeutics based on β-cell rest.

Temporal Transcriptome Analysis Reveals Dynamic Gene Expression Patterns Driving β-Cell Maturation

Newly differentiated pancreatic β cells lack proper insulin secretion profiles of mature functional β cells. The global gene expression differences between paired immature and mature β cells have been studied, but the dynamics of transcriptional events, correlating with temporal development of glucose-stimulated insulin secretion (GSIS), remain to be fully defined. This aspect is important to identify which genes and pathways are necessary for β-cell development or for maturation, as defective insulin secretion is linked with diseases such as diabetes. In this study, we assayed through RNA sequencing the global gene expression across six β-cell developmental stages in mice, spanning from β-cell progenitor to mature β cells. A computational pipeline then selected genes differentially expressed with respect to progenitors and clustered them into groups with distinct temporal patterns associated with biological functions and pathways. These patterns were finally correlated with experimental GSIS, calcium influx, and insulin granule formation data. Gene expression temporal profiling revealed the timing of important biological processes across β-cell maturation, such as the deregulation of β-cell developmental pathways and the activation of molecular machineries for vesicle biosynthesis and transport, signal transduction of transmembrane receptors, and glucose-induced Ca²⁺ influx, which were established over a week before β-cell maturation completes. In particular, β cells developed robust insulin secretion at high glucose several days after birth, coincident with the establishment of glucose-induced calcium influx. Yet the neonatal β cells displayed high basal insulin secretion, which decreased to the low levels found in mature β cells only a week later. Different genes associated with calcium-mediated processes, whose alterations are linked with insulin resistance and deregulation of glucose homeostasis, showed increased expression across β-cell stages, in accordance with the temporal acquisition of proper GSIS. Our temporal gene expression pattern analysis provided a comprehensive database of the underlying molecular components and biological mechanisms driving β-cell maturation at different temporal stages, which are fundamental for better control of the in vitro production of functional β cells from human embryonic stem/induced pluripotent cell for transplantation-based type 1 diabetes therapy.

Temporal Dynamic Methods for Bulk RNA-Seq Time Series Data

Dynamic studies in time course experimental designs and clinical approaches have been widely used by the biomedical community. These applications are particularly relevant in stimuli-response models under environmental conditions, characterization of gradient biological processes in developmental biology, identification of therapeutic effects in clinical trials, disease progressive models, cell-cycle, and circadian periodicity. Despite their feasibility and popularity, sophisticated dynamic methods that are well validated in large-scale comparative studies, in terms of statistical and computational rigor, are less benchmarked, comparing to their static counterparts. To date, a number of novel methods in bulk RNA-Seq data have been developed for the various time-dependent stimuli, circadian rhythms, cell-lineage in differentiation, and disease progression. Here, we comprehensively review a key set of representative dynamic strategies and discuss current issues associated with the detection of dynamically changing genes. We also provide recommendations for future directions for studying non-periodical, periodical time course data, and meta-dynamic datasets.

Pooled analysis of radiation hybrids identifies loci for growth and drug action in mammalian cells

Genetic screens in mammalian cells commonly focus on loss-of-function approaches. To evaluate the phenotypic consequences of extra gene copies, we used bulk segregant analysis (BSA) of radiation hybrid (RH) cells. We constructed six pools of RH cells, each consisting of ∼2500 independent clones, and placed the pools under selection in media with or without paclitaxel. Low pass sequencing identified 859 growth loci, 38 paclitaxel loci, 62 interaction loci, and three loci for mitochondrial abundance at genome-wide significance. Resolution was measured as ∼30 kb, close to single-gene. Divergent properties were displayed by the RH-BSA growth genes compared to those from loss-of-function screens, refuting the balance hypothesis. In addition, enhanced retention of human centromeres in the RH pools suggests a new approach to functional dissection of these chromosomal elements. Pooled analysis of RH cells showed high power and resolution and should be a useful addition to the mammalian genetic toolkit.

High-dimensionality Data Analysis of Pharmacological Systems Associated with Complex Diseases

It is widely accepted that molecular reductionist views of highly complex human physiologic activity, e.g., the aging process, as well as therapeutic drug efficacy are largely oversimplifications. Currently some of the most effective appreciation of biologic disease and drug response complexity is achieved using high-dimensionality (H-D) data streams from transcriptomic, proteomic, metabolomics, or epigenomic pipelines. Multiple H-D data sets are now common and freely accessible for complex diseases such as metabolic syndrome, cardiovascular disease, and neurodegenerative conditions such as Alzheimer's disease. Over the last decade our ability to interrogate these high-dimensionality data streams has been profoundly enhanced through the development and implementation of highly effective bioinformatic platforms. Employing these computational approaches to understand the complexity of age-related diseases provides a facile mechanism to then synergize this pathologic appreciation with a similar level of understanding of therapeutic-mediated signaling. For informative pathology and drug-based analytics that are able to generate meaningful therapeutic insight across diverse data streams, novel informatics processes such as latent semantic indexing and topological data analyses will likely be important. Elucidation of H-D molecular disease signatures from diverse data streams will likely generate and refine new therapeutic strategies that will be designed with a cognizance of a realistic appreciation of the complexity of human age-related disease and drug effects. We contend that informatic platforms should be synergistic with more advanced chemical/drug and phenotypic cellular/tissue-based analytical predictive models to assist in either de novo drug prioritization or effective repurposing for the intervention of aging-related diseases. SIGNIFICANCE STATEMENT: All diseases, as well as pharmacological mechanisms, are far more complex than previously thought a decade ago. With the advent of commonplace access to technologies that produce large volumes of high-dimensionality data (e.g., transcriptomics, proteomics, metabolomics), it is now imperative that effective tools to appreciate this highly nuanced data are developed. Being able to appreciate the subtleties of high-dimensionality data will allow molecular pharmacologists to develop the most effective multidimensional therapeutics with effectively engineered efficacy profiles.

Improving Mycoplasma ovipneumoniae culture medium by a comparative transcriptome method

Mycoplasma ovipneumoniae (Mo) is difficult to culture, resulting in many difficulties in related research and application. Since nucleotide metabolism is a basic metabolism affects growth, this study conducted a “point-to-point” comparison of the corresponding growth phases between the Mo NM151 strain and the Mycoplasma mycoides subsp. capri (Mmc) PG3 strain. The results showed that the largest difference in nucleotide metabolism was found in the stationary phase. Nucleotide synthesis in PG3 was mostly de novo, while nucleotide synthesis in NM151 was primarily based on salvage synthesis. Compared with PG3, the missing reactions of NM151 referred to the synthesis of deoxythymine monophosphate. We proposed and validated a culture medium with added serine to fill this gap and prolong the stationary phase of NM151. This solved the problem of the fast death of Mo, which is significant for related research and application.

TimeMeter assesses temporal gene expression similarity and identifies differentially progressing genes

Comparative time series transcriptome analysis is a powerful tool to study development, evolution, aging, disease progression and cancer prognosis. We develop TimeMeter, a statistical method and tool to assess temporal gene expression similarity, and identify differentially progressing genes where one pattern is more temporally advanced than the other. We apply TimeMeter to several datasets, and show that TimeMeter is capable of characterizing complicated temporal gene expression associations. Interestingly, we find: (i) the measurement of differential progression provides a novel feature in addition to pattern similarity that can characterize early developmental divergence between two species; (ii) genes exhibiting similar temporal patterns between human and mouse during neural differentiation are under strong negative (purifying) selection during evolution; (iii) analysis of genes with similar temporal patterns in mouse digit regeneration and axolotl blastema differentiation reveals common gene groups for appendage regeneration with potential implications in regenerative medicine.

State-Transition Analysis of Time-Sequential Gene Expression Identifies Critical Points That Predict Development of Acute Myeloid Leukemia

Temporal dynamics of gene expression inform cellular and molecular perturbations associated with disease development and evolution. Given the complexity of high-dimensional temporal genomic data, an analytic framework guided by a robust theory is needed to interpret time-sequential changes and to predict system dynamics. Here we model temporal dynamics of the transcriptome of peripheral blood mononuclear cells in a two-dimensional state-space representing states of health and leukemia using time-sequential bulk RNA-seq data from a murine model of acute myeloid leukemia (AML). The state-transition model identified critical points that accurately predict AML development and identifies stepwise transcriptomic perturbations that drive leukemia progression. The geometry of the transcriptome state-space provided a biological interpretation of gene dynamics, aligned gene signals that are not synchronized in time across mice, and allowed quantification of gene and pathway contributions to leukemia development. Our state-transition model synthesizes information from multiple cell types in the peripheral blood and identifies critical points in the transition from health to leukemia to guide interpretation of changes in the transcriptome as a whole to predict disease progression. SIGNIFICANCE: These findings apply the theory of state transitions to model the initiation and development of acute myeloid leukemia, identifying transcriptomic perturbations that accurately predict time to disease development.See related commentary by Kuijjer, p. 3072 GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/15/3157/F1.large.jpg.

Trajectory-based differential expression analysis for single-cell sequencing data

Trajectory inference has radically enhanced single-cell RNA-seq research by enabling the study of dynamic changes in gene expression. Downstream of trajectory inference, it is vital to discover genes that are (i) associated with the lineages in the trajectory, or (ii) differentially expressed between lineages, to illuminate the underlying biological processes. Current data analysis procedures, however, either fail to exploit the continuous resolution provided by trajectory inference, or fail to pinpoint the exact types of differential expression. We introduce tradeSeq, a powerful generalized additive model framework based on the negative binomial distribution that allows flexible inference of both within-lineage and between-lineage differential expression. By incorporating observation-level weights, the model additionally allows to account for zero inflation. We evaluate the method on simulated datasets and on real datasets from droplet-based and full-length protocols, and show that it yields biological insights through a clear interpretation of the data.

Comparative analysis of differential gene expression tools for RNA sequencing time course data

RNA sequencing (RNA-seq) has become a standard procedure to investigate transcriptional changes between conditions and is routinely used in research and clinics. While standard differential expression (DE) analysis between two conditions has been extensively studied, and improved over the past decades, RNA-seq time course (TC) DE analysis algorithms are still in their early stages. In this study, we compare, for the first time, existing TC RNA-seq tools on an extensive simulation data set and validated the best performing tools on published data. Surprisingly, TC tools were outperformed by the classical pairwise comparison approach on short time series (<8 time points) in terms of overall performance and robustness to noise, mostly because of high number of false positives, with the exception of ImpulseDE2. Overlapping of candidate lists between tools improved this shortcoming, as the majority of false-positive, but not true-positive, candidates were unique for each method. On longer time series, pairwise approach was less efficient on the overall performance compared with splineTC and maSigPro, which did not identify any false-positive candidate.

Temporal dynamics in meta longitudinal RNA-Seq data

Identification of differentially expressed genes has been a high priority task of downstream analyses to further advances in biomedical research. Investigators have been faced with an array of issues in dealing with more complicated experiments and metadata, including batch effects, normalization, temporal dynamics (temporally differential expression), and isoform diversity (isoform-level quantification and differential splicing events). To date, there are currently no standard approaches to precisely and efficiently analyze these moderate or large-scale experimental designs, especially with combined metadata. In this report, we propose comprehensive analytical pipelines to precisely characterize temporal dynamics in differential expression of genes and other genomic features, i.e., the variability of transcripts, isoforms and exons, by controlling batch effects and other nuisance factors that could have significant confounding effects on the main effects of interest in comparative models and may result in misleading interpretations.

Subject level clustering using a negative binomial model for small transcriptomic studies

Background

Unsupervised clustering represents one of the most widely applied methods in analysis of high-throughput ‘omics data. A variety of unsupervised model-based or parametric clustering methods and non-parametric clustering methods have been proposed for RNA-seq count data, most of which perform well for large samples, e.g. N ≥ 500. A common issue when analyzing limited samples of RNA-seq count data is that the data follows an over-dispersed distribution, and thus a Negative Binomial likelihood model is often used. Thus, we have developed a Negative Binomial model-based (NBMB) clustering approach for application to RNA-seq studies.

Results

We have developed a Negative Binomial Model-Based (NBMB) method to cluster samples using a stochastic version of the expectation-maximization algorithm. A simulation study involving various scenarios was completed to compare the performance of NBMB to Gaussian model-based or Gaussian mixture modeling (GMM). NBMB was also applied for the clustering of two RNA-seq studies; type 2 diabetes study (N = 96) and TCGA study of ovarian cancer (N = 295). Simulation results showed that NBMB outperforms GMM applied with different transformations in majority of scenarios with limited sample size. Additionally, we found that NBMB outperformed GMM for small clusters distance regardless of sample size. Increasing total number of genes with fixed proportion of differentially expressed genes does not change the outperformance of NBMB, but improves the overall performance of GMM. Analysis of type 2 diabetes and ovarian cancer tumor data with NBMB found good agreement with the reported disease subtypes and the gene expression patterns. This method is available in an R package on CRAN named NB.MClust.

Conclusion

Use of Negative Binomial model based clustering is advisable when clustering over dispersed RNA-seq count data.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2556-9) contains supplementary material, which is available to authorized users.

The soil organic matter decomposition mechanisms in ectomycorrhizal fungi are tuned for liberating soil organic nitrogen

Many trees form ectomycorrhizal symbiosis with fungi. During symbiosis, the tree roots supply sugar to the fungi in exchange for nitrogen, and this process is critical for the nitrogen and carbon cycles in forest ecosystems. However, the extents to which ectomycorrhizal fungi can liberate nitrogen and modify the soil organic matter and the mechanisms by which they do so remain unclear since they have lost many enzymes for litter decomposition that were present in their free-living, saprotrophic ancestors. Using time-series spectroscopy and transcriptomics, we examined the ability of two ectomycorrhizal fungi from two independently evolved ectomycorrhizal lineages to mobilize soil organic nitrogen. Both species oxidized the organic matter and accessed the organic nitrogen. The expression of those events was controlled by the availability of glucose and inorganic nitrogen. Despite those similarities, the decomposition mechanisms, including the type of genes involved as well as the patterns of their expression, differed markedly between the two species. Our results suggest that in agreement with their diverse evolutionary origins, ectomycorrhizal fungi use different decomposition mechanisms to access organic nitrogen entrapped in soil organic matter. The timing and magnitude of the expression of the decomposition activity can be controlled by the below-ground nitrogen quality and the above-ground carbon supply.

Trendy: segmented regression analysis of expression dynamics in high-throughput ordered profiling experiments

BACKGROUND

High-throughput expression profiling experiments with ordered conditions (e.g. time-course or spatial-course) are becoming more common for studying detailed differentiation processes or spatial patterns. Identifying dynamic changes at both the individual gene and whole transcriptome level can provide important insights about genes, pathways, and critical time points.

RESULTS

We present an R package, Trendy, which utilizes segmented regression models to simultaneously characterize each gene's expression pattern and summarize overall dynamic activity in ordered condition experiments. For each gene, Trendy finds the optimal segmented regression model and provides the location and direction of dynamic changes in expression. We demonstrate the utility of Trendy to provide biologically relevant results on both microarray and RNA-sequencing (RNA-seq) datasets.

CONCLUSIONS

Trendy is a flexible R package which characterizes gene-specific expression patterns and summarizes changes of global dynamics over ordered conditions. Trendy is freely available on Bioconductor with a full vignette at https://bioconductor.org/packages/release/bioc/html/Trendy.html .

HiMAP: Robust phylogenomics from highly multiplexed amplicon sequencing

High-throughput sequencing has fundamentally changed how molecular phylogenetic data sets are assembled, and phylogenomic data sets commonly contain 50- to 100-fold more loci than those generated using traditional Sanger sequencing-based approaches. Here, we demonstrate a new approach for building phylogenomic data sets using single-tube, highly multiplexed amplicon sequencing, which we name HiMAP (highly multiplexed amplicon-based phylogenomics) and present bioinformatic pipelines for locus selection based on genomic and transcriptomic data resources and postsequencing consensus calling and alignment. This method is inexpensive and amenable to sequencing a large number (hundreds) of taxa simultaneously and requires minimal hands-on time at the bench (<1/2 day), and data analysis can be accomplished without the need for read mapping or assembly. We demonstrate this approach by sequencing 878 amplicons in single reactions for 82 species of tephritid fruit flies across seven genera (384 individuals), including some of the most economically important agricultural insect pests. The resulting filtered data set (>150,000-bp concatenated alignment, ~20% missing character sites across all individuals and amplicons) contained >40,000 phylogenetically informative characters, and although some discordance was observed between analyses, it provided unparalleled resolution of many phylogenetic relationships in this group. Most notably, we found high support for the generic status of Zeugodacus and the sister relationship between Dacus and Zeugodacus. We discuss HiMAP, with regard to its molecular and bioinformatic strengths, and the insight the resulting data set provides into relationships of this diverse insect group.

Insight into Genes Regulating Postharvest Aflatoxin Contamination of Tetraploid Peanut from Transcriptional Profiling

Postharvest aflatoxin contamination is a challenging issue that affects peanut quality. Aflatoxin is produced by fungi belonging to the Aspergilli group, and is known as an acutely toxic, carcinogenic, and immune-suppressing class of mycotoxins. Evidence for several host genetic factors that may impact aflatoxin contamination has been reported, e.g., genes for lipoxygenase (PnLOX1 and PnLOX2/PnLOX3 that showed either positive or negative regulation with Aspergillus infection), reactive oxygen species, and WRKY (highly associated with or differentially expressed upon infection of maize with Aspergillus flavus); however, their roles remain unclear. Therefore, we conducted an RNA-sequencing experiment to differentiate gene response to the infection by A. flavus between resistant (ICG 1471) and susceptible (Florida-07) cultivated peanut genotypes. The gene expression profiling analysis was designed to reveal differentially expressed genes in response to the infection (infected vs. mock-treated seeds). In addition, the differential expression of the fungal genes was profiled. The study revealed the complexity of the interaction between the fungus and peanut seeds as the expression of a large number of genes was altered, including some in the process of plant defense to aflatoxin accumulation. Analysis of the experimental data with "keggseq," a novel designed tool for Kyoto Encyclopedia of Genes and Genomes enrichment analysis, showed the importance of α-linolenic acid metabolism, protein processing in the endoplasmic reticulum, spliceosome, and carbon fixation and metabolism pathways in conditioning resistance to aflatoxin accumulation. In addition, coexpression network analysis was carried out to reveal the correlation of gene expression among peanut and fungal genes. The results showed the importance of WRKY, toll/Interleukin1 receptor-nucleotide binding site leucine-rich repeat (TIR-NBS-LRR), ethylene, and heat shock proteins in the resistance mechanism.

Mycobacterium tuberculosis-induced miR-155 subverts autophagy by targeting ATG3 in human dendritic cells

Autophagy is a primordial eukaryotic pathway, which provides the immune system with multiple mechanisms for the elimination of invading pathogens including Mycobacterium tuberculosis (Mtb). As a consequence, Mtb has evolved different strategies to hijack the autophagy process. Given the crucial role of human primary dendritic cells (DC) in host immunity control, we characterized Mtb-DC interplay by studying the contribution of cellular microRNAs (miRNAs) in the post-transcriptional regulation of autophagy related genes. From the expression profile of de-regulated miRNAs obtained in Mtb-infected human DC, we identified 7 miRNAs whose expression was previously found to be altered in specimens of TB patients. Among them, gene ontology analysis showed that miR-155, miR-155* and miR-146a target mRNAs with a significant enrichment in biological processes linked to autophagy. Interestingly, miR-155 was significantly stimulated by live and virulent Mtb and enriched in polysome-associated RNA fraction, where actively translated mRNAs reside. The putative pair interaction among the E2 conjugating enzyme involved in LC3-lipidation and autophagosome formation-ATG3-and miR-155 arose by target prediction analysis, was confirmed by both luciferase reporter assay and Atg3 immunoblotting analysis of miR-155-transfected DC, which showed also a consistent Atg3 protein and LC3 lipidated form reduction. Late in infection, when miR-155 expression peaked, both the level of Atg3 and the number of LC3 puncta per cell (autophagosomes) decreased dramatically. In accordance, miR-155 silencing rescued autophagosome number in Mtb infected DC and enhanced autolysosome fusion, thereby supporting a previously unidentified role of the miR-155 as inhibitor of ATG3 expression. Taken together, our findings suggest how Mtb can manipulate cellular miRNA expression to regulate Atg3 for its own survival, and highlight the importance to develop novel therapeutic strategies against tuberculosis that would boost autophagy.

Mycobacterium tuberculosis (Mtb) is one of the most successful pathogens in human history and remains the second leading cause of death from an infectious agent worldwide. The major reason of Mtb success relies on its ability to evade host immunity. Autophagy, a cellular mechanism involved in intracellular pathogen elimination, is one of the pathways hijacked by Mtb to elude the control of dendritic cells (DC), major cellular effectors of immune response. Recently, it has become clear that Mtb infection not only alters cellular gene expression, but also controls the level of small RNA molecules, namely microRNAs (miRNAs), which function as negative regulators of mRNA translation into protein. In the present study, we observed that the infection of human DC with Mtb leads to a strong induction of host miR-155, a critical regulator of host immune response. By mean of miR-155 induction, Mtb reduces Atg3 protein content, a crucial enzyme needed for the initial phase of the autophagic process. Interestingly, miR-155 silencing during Mtb infection restores Atg3 level and rescues autophagy. These findings contribute to better elucidate Mtb-triggered escape mechanisms and highlight the importance to develop host-directed therapies to combat tuberculosis based on autophagy boosting.

De novo assembly and comparative transcriptome characterization of Poecilobdella javanica provide insight into blood feeding of medicinal leeches

Leeches (family Hirudinidae) are classic model invertebrates used in diverse clinical treatments, such as reconstructive microsurgery, hypertension, and gangrene treatment.

Comparative transcriptional analysis of hop responses to infection with Verticillium nonalfalfae

KEY MESSAGE

Dynamic transcriptome profiling revealed excessive, yet ineffective, immune response to V. nonalfalfae infection in susceptible hop, global gene downregulation in shoots of resistant hop and only a few infection-associated genes in roots. Hop (Humulus lupulus L.) production is hampered by Verticillium wilt, a disease predominantly caused by the soil-borne fungus Verticillium nonalfalfae. Only a few hop cultivars exhibit resistance towards it and mechanisms of this resistance have not been discovered. In this study, we compared global transcriptional responses in roots and shoots of resistant and susceptible hop plants infected by a lethal strain of V. nonalfalfae. Time-series differential gene expression profiles between infected and mock inoculated plants were determined and subjected to network-based analysis of functional enrichment. In the resistant hop cultivar, a remarkably low number of genes were differentially expressed in roots in response to V. nonalfalfae infection, while the majority of differentially expressed genes were down-regulated in shoots. The most significantly affected genes were related to cutin biosynthesis, cell wall biogenesis, lateral root development and terpenoid biosynthesis. On the other hand, susceptible hop exhibited a strong defence response in shoots and roots, including increased expression of genes associated with plant responses, such as innate immunity, wounding, jasmonic acid pathway and chitinase activity. Strong induction of defence-associated genes in susceptible hop and a low number of infection-responsive genes in the roots of resistant hop are consistent with previous findings, confirming the pattern of excessive response of the susceptible cultivar, which ultimately fails to protect the plant from V. nonalfalfae. This research offers a multifaceted overview of transcriptional responses of susceptible and resistant hop cultivars to V. nonalfalfae infection and represents a valuable resource in the study of this plant-pathogen interaction.

Preface: BITS2014, the annual meeting of the Italian Society of Bioinformatics

This Preface introduces the content of the BioMed Central journal Supplements related to BITS2014 meeting, held in Rome, Italy, from the 26th to the 28th of February, 2014.