Search

Background Count data generated by next-generation sequencing assays do not measure absolute transcript abundances. Instead, the data are constrained to an arbitrary “library size” by the sequencing depth of the assay, and typically must be normalized prior to statistical analysis. The constrained nature of these data means one could alternatively use a log-ratio transformation in lieu of normalization, as often done when testing for differential abundance (DA) of operational taxonomic units (OTUs) in 16S rRNA data. Therefore, we benchmark how well the ALDEx2 package, a transformation-based DA tool, detects differential expression in high-throughput RNA-sequencing data (RNA-Seq), compared to conventional RNA-Seq methods such as edgeR and DESeq2. Results To evaluate the performance of log-ratio transformation-based tools, we apply the ALDEx2 package to two simulated, and two real, RNA-Seq data sets. One of the latter was previously used to benchmark dozens of conventional RNA-Seq differential expression methods, enabling us to directly compare transformation-based approaches. We show that ALDEx2, widely used in meta-genomics research, identifies differentially expressed genes (and transcripts) from RNA-Seq data with high precision and, given sufficient sample sizes, high recall too (regardless of the alignment and quantification procedure used). Although we show that the choice in log-ratio transformation can affect performance, ALDEx2 has high precision (i.e., few false positives) across all transformations. Finally, we present a novel, iterative log-ratio transformation (now implemented in ALDEx2) that further improves performance in simulations. Conclusions Our results suggest that log-ratio transformation-based methods can work to measure differential expression from RNA-Seq data, provided that certain assumptions are met. Moreover, these methods have very high precision (i.e., few false positives) in simulations a

BACKGROUND: Count data generated by next-generation sequencing assays do not measure absolute transcript abundances. Instead, the data are constrained to an arbitrary "library size" by the sequencing depth of the assay, and typically must be normalized prior to statistical analysis. The constrained nature of these data means one could alternatively use a log-ratio transformation in lieu of normalization, as often done when testing for differential abundance (DA) of operational taxonomic units (OTUs) in 16S rRNA data. Therefore, we benchmark how well the ALDEx2 package, a transformation-based DA tool, detects differential expression in high-throughput RNA-sequencing data (RNA-Seq), compared to conventional RNA-Seq methods such as edgeR and DESeq2. RESULTS: To evaluate the performance of log-ratio transformation-based tools, we apply the ALDEx2 package to two simulated, and two real, RNA-Seq data sets. One of the latter was previously used to benchmark dozens of conventional RNA-Seq differential expression methods, enabling us to directly compare transformation-based approaches. We show that ALDEx2, widely used in meta-genomics research, identifies differentially expressed genes (and transcripts) from RNA-Seq data with high precision and, given sufficient sample sizes, high recall too (regardless of the alignment and quantification procedure used). Although we show that the choice in log-ratio transformation can affect performance, ALDEx2 has high precision (i.e., few false positives) across all transformations. Finally, we present a novel, iterative log-ratio transformation (now implemented in ALDEx2) that further improves performance in simulations. CONCLUSIONS: Our results suggest that log-ratio transformation-based methods can work to measure differential expression from RNA-Seq data, provided that certain assumptions are met. Moreover, these methods have very high precision (i.e., few false positives) in simulations and perform w

[Background]: The identification of genes differentially expressed in the skeletal muscle of pigs displaying distinct growth and fatness profiles might contribute to identify the genetic factors that influence the phenotypic variation of such traits. So far, the majority of porcine transcriptomic studies have investigated differences in gene expression at a global scale rather than at the mRNA isoform level. In the current work, we have investigated the differential expression of mRNA isoforms in the gluteus medius (GM) muscle of 52 Duroc HIGH (increased backfat thickness, intramuscular fat and saturated and monounsaturated fatty acids contents) and LOW pigs (opposite phenotype, with an increased polyunsaturated fatty acids content)., [Results]: Our analysis revealed that 10.9% of genes expressed in the GM muscle generate alternative mRNA isoforms, with an average of 2.9 transcripts per gene. By using two different pipelines, one based on the CLC Genomics Workbench and another one on the STAR, RSEM and DESeq2 softwares, we have identified 10 mRNA isoforms that both pipelines categorize as differentially expressed in HIGH vs LOW pigs (P-value < 0.01 and ±0.6 log2fold-change). Only five mRNA isoforms, produced by the ITGA5, SEMA4D, LITAF, TIMP1 and ANXA2 genes, remain significant after correction for multiple testing (q-value < 0.05 and ±0.6 log2fold-change), being upregulated in HIGH pigs., [Conclusions]: The increased levels of specific ITGA5, LITAF, TIMP1 and ANXA2 mRNA isoforms in HIGH pigs is consistent with reports indicating that the overexpression of these four genes is associated with obesity and metabolic disorders in humans. A broader knowledge about the functional attributes of these mRNA variants would be fundamental to elucidate the consequences of transcript diversity on the determinism of porcine phenotypes of economic interest.

Background: Asparagine is a very important nitrogen transport and storage compound in plants due to its high nitrogen/carbon ratio and stability. Asparagine intracellu lar concentration depends on a balance between asparagine biosynthesis and degradation. The main enzymes involved in asparagine metabolism are as paragine synthetase (ASN), asparaginase (NSE) and serine-glyoxylate aminotransfera se (SGAT). The study of the genes encoding for these enzymes in the model legume Lotus japonicus is of particular interest since it has been proposed that asparagine is the principal molecule used to transport reduced nitrogen within the plant in most temperate legumes. Results: A differential expression of genes encoding for seve ral enzymes involved in asparagine metabolism was detected in L. japonicus . ASN is encoded by three genes, LjASN1 was the most highly expressed in mature leaves while LjASN2 expression was negligible and LjASN3 showed a low expression in this organ, suggesting that LjASN1 is the main gene responsible for asparagine synthesis in mature leaves. In young leaves, LjASN3 was the only ASN gene expressed although at low levels, while all the three genes encoding for NSE were highly expressed, especially LjNSE1 .Innodules, LjASN2 and LjNSE2 were the most highly expressed genes, suggesting an important role for these genes in this organ. Several lines of evidence support the connection between asparagine metabolic genes and photorespiration in L. japonicus : a) a mutant plant deficient in LjNSE1 showed a dramatic decrease in the expression of the two genes encoding for SGAT; b) expression of the genes involved in asparagine metabolism is altered in a photorespiratory mutant lacking plastidic glutamine synthetase; c) a clustering analysis indicated a similar pattern of expression among several genes involved in photorespiratory and asparagine metabolism, indicating a clear link between LjASN1 and LjSGAT genes and photorespiration. Conclusions: The results obt

Background: Asparagine is a very important nitrogen transport and storage compound in plants due to its high nitrogen/carbon ratio and stability. Asparagine intracellu lar concentration depends on a balance between asparagine biosynthesis and degradation. The main enzymes involved in asparagine metabolism are as paragine synthetase (ASN), asparaginase (NSE) and serine-glyoxylate aminotransfera se (SGAT). The study of the genes encoding for these enzymes in the model legume Lotus japonicus is of particular interest since it has been proposed that asparagine is the principal molecule used to transport reduced nitrogen within the plant in most temperate legumes. Results: A differential expression of genes encoding for seve ral enzymes involved in asparagine metabolism was detected in L. japonicus . ASN is encoded by three genes, LjASN1 was the most highly expressed in mature leaves while LjASN2 expression was negligible and LjASN3 showed a low expression in this organ, suggesting that LjASN1 is the main gene responsible for asparagine synthesis in mature leaves. In young leaves, LjASN3 was the only ASN gene expressed although at low levels, while all the three genes encoding for NSE were highly expressed, especially LjNSE1 .Innodules, LjASN2 and LjNSE2 were the most highly expressed genes, suggesting an important role for these genes in this organ. Several lines of evidence support the connection between asparagine metabolic genes and photorespiration in L. japonicus : a) a mutant plant deficient in LjNSE1 showed a dramatic decrease in the expression of the two genes encoding for SGAT; b) expression of the genes involved in asparagine metabolism is altered in a photorespiratory mutant lacking plastidic glutamine synthetase; c) a clustering analysis indicated a similar pattern of expression among several genes involved in photorespiratory and asparagine metabolism, indicating a clear link between LjASN1 and LjSGAT genes and photorespiration. Conclusions: The results obt

Background: Asparagine is a very important nitrogen transport and storage compound in plants due to its high nitrogen/carbon ratio and stability. Asparagine intracellu lar concentration depends on a balance between asparagine biosynthesis and degradation. The main enzymes involved in asparagine metabolism are as paragine synthetase (ASN), asparaginase (NSE) and serine-glyoxylate aminotransfera se (SGAT). The study of the genes encoding for these enzymes in the model legume Lotus japonicus is of particular interest since it has been proposed that asparagine is the principal molecule used to transport reduced nitrogen within the plant in most temperate legumes. Results: A differential expression of genes encoding for seve ral enzymes involved in asparagine metabolism was detected in L. japonicus . ASN is encoded by three genes, LjASN1 was the most highly expressed in mature leaves while LjASN2 expression was negligible and LjASN3 showed a low expression in this organ, suggesting that LjASN1 is the main gene responsible for asparagine synthesis in mature leaves. In young leaves, LjASN3 was the only ASN gene expressed although at low levels, while all the three genes encoding for NSE were highly expressed, especially LjNSE1 .Innodules, LjASN2 and LjNSE2 were the most highly expressed genes, suggesting an important role for these genes in this organ. Several lines of evidence support the connection between asparagine metabolic genes and photorespiration in L. japonicus : a) a mutant plant deficient in LjNSE1 showed a dramatic decrease in the expression of the two genes encoding for SGAT; b) expression of the genes involved in asparagine metabolism is altered in a photorespiratory mutant lacking plastidic glutamine synthetase; c) a clustering analysis indicated a similar pattern of expression among several genes involved in photorespiratory and asparagine metabolism, indicating a clear link between LjASN1 and LjSGAT genes and photorespiration. Conclusions: The results obt

Background: Edible bird's nest (EBN), produced from solidified saliva secretions of specific swiftlet species during the breeding season, is one of the most valuable animal by-products in the world. The composition and medicinal benefits of EBN have been extensively studied, however, genomic and transcriptomic studies of the salivary glands of these birds have not been conducted. Results: The study described the transcriptomes of salivary glands from three swiftlet species (28 samples) generated by RNASeq. A total of 14,835 annotated genes and 428 unmapped genes were cataloged. The current study investigated the genes and pathways that are associated with the development of salivary gland and EBN composition. Differential expression and pathway enrichment analysis indicated that the expression of CREB3L2 and several signaling pathways involved in salivary gland development, namely, the EGFR, BMP, and MAPK signaling pathways, were up-regulated in swiftlets producing white EBN (Aerodramus fuciphagus) and black EBN (Aerodramus maximus) compared with non-EBN-producing swiftlets (Apus affinis). Furthermore, MGAT, an essential gene for the biosynthesis of N-acetylneuraminic acid (sialic acid), was highly expressed in both white- and black-nest swiftlets compared to non-EBN-producing swiftlets. Interspecies comparison between Aerodramus fuciphagus and Aerodramus maximus indicated that the genes involved in N-acetylneuraminic and fatty acid synthesis were up-regulated in Aerodramus fuciphagus, while alanine and aspartate synthesis pathways were up-regulated in Aerodramus maximus. Furthermore, gender-based analysis revealed that N-glycan trimming pathway was significantly up-regulated in male Aerodramus fuciphagus from its natural habitat (cave) compared to their female counterpart. Conclusions:Transcriptomic analysis of salivary glands of different swiftlet species reveal differential expressions of candidate genes that are involved in salivary gland development and in the

Background: Edible bird's nest (EBN), produced from solidified saliva secretions of specific swiftlet species during the breeding season, is one of the most valuable animal by-products in the world. The composition and medicinal benefits of EBN have been extensively studied, however, genomic and transcriptomic studies of the salivary glands of these birds have not been conducted. Results: The study described the transcriptomes of salivary glands from three swiftlet species (28 samples) generated by RNASeq. A total of 14,835 annotated genes and 428 unmapped genes were cataloged. The current study investigated the genes and pathways that are associated with the development of salivary gland and EBN composition. Differential expression and pathway enrichment analysis indicated that the expression of CREB3L2 and several signaling pathways involved in salivary gland development, namely, the EGFR, BMP, and MAPK signaling pathways, were up-regulated in swiftlets producing white EBN (Aerodramus fuciphagus) and black EBN (Aerodramus maximus) compared with non-EBN-producing swiftlets (Apus affinis). Furthermore, MGAT, an essential gene for the biosynthesis of N-acetylneuraminic acid (sialic acid), was highly expressed in both white- and black-nest swiftlets compared to non-EBN-producing swiftlets. Interspecies comparison between Aerodramus fuciphagus and Aerodramus maximus indicated that the genes involved in N-acetylneuraminic and fatty acid synthesis were up-regulated in Aerodramus fuciphagus, while alanine and aspartate synthesis pathways were up-regulated in Aerodramus maximus. Furthermore, gender-based analysis revealed that N-glycan trimming pathway was significantly up-regulated in male Aerodramus fuciphagus from its natural habitat (cave) compared to their female counterpart. Conclusions:Transcriptomic analysis of salivary glands of different swiftlet species reveal differential expressions of candidate genes that are involved in salivary gland development and in the

Background: Asparagine is a very important nitrogen transport and storage compound in plants due to its high nitrogen/carbon ratio and stability. Asparagine intracellu lar concentration depends on a balance between asparagine biosynthesis and degradation. The main enzymes involved in asparagine metabolism are as paragine synthetase (ASN), asparaginase (NSE) and serine-glyoxylate aminotransfera se (SGAT). The study of the genes encoding for these enzymes in the model legume Lotus japonicus is of particular interest since it has been proposed that asparagine is the principal molecule used to transport reduced nitrogen within the plant in most temperate legumes. Results: A differential expression of genes encoding for seve ral enzymes involved in asparagine metabolism was detected in L. japonicus . ASN is encoded by three genes, LjASN1 was the most highly expressed in mature leaves while LjASN2 expression was negligible and LjASN3 showed a low expression in this organ, suggesting that LjASN1 is the main gene responsible for asparagine synthesis in mature leaves. In young leaves, LjASN3 was the only ASN gene expressed although at low levels, while all the three genes encoding for NSE were highly expressed, especially LjNSE1 .Innodules, LjASN2 and LjNSE2 were the most highly expressed genes, suggesting an important role for these genes in this organ. Several lines of evidence support the connection between asparagine metabolic genes and photorespiration in L. japonicus : a) a mutant plant deficient in LjNSE1 showed a dramatic decrease in the expression of the two genes encoding for SGAT; b) expression of the genes involved in asparagine metabolism is altered in a photorespiratory mutant lacking plastidic glutamine synthetase; c) a clustering analysis indicated a similar pattern of expression among several genes involved in photorespiratory and asparagine metabolism, indicating a clear link between LjASN1 and LjSGAT genes and photorespiration. Conclusions: The results obt

Background: Asparagine is a very important nitrogen transport and storage compound in plants due to its high nitrogen/carbon ratio and stability. Asparagine intracellu lar concentration depends on a balance between asparagine biosynthesis and degradation. The main enzymes involved in asparagine metabolism are as paragine synthetase (ASN), asparaginase (NSE) and serine-glyoxylate aminotransfera se (SGAT). The study of the genes encoding for these enzymes in the model legume Lotus japonicus is of particular interest since it has been proposed that asparagine is the principal molecule used to transport reduced nitrogen within the plant in most temperate legumes. Results: A differential expression of genes encoding for seve ral enzymes involved in asparagine metabolism was detected in L. japonicus . ASN is encoded by three genes, LjASN1 was the most highly expressed in mature leaves while LjASN2 expression was negligible and LjASN3 showed a low expression in this organ, suggesting that LjASN1 is the main gene responsible for asparagine synthesis in mature leaves. In young leaves, LjASN3 was the only ASN gene expressed although at low levels, while all the three genes encoding for NSE were highly expressed, especially LjNSE1 .Innodules, LjASN2 and LjNSE2 were the most highly expressed genes, suggesting an important role for these genes in this organ. Several lines of evidence support the connection between asparagine metabolic genes and photorespiration in L. japonicus : a) a mutant plant deficient in LjNSE1 showed a dramatic decrease in the expression of the two genes encoding for SGAT; b) expression of the genes involved in asparagine metabolism is altered in a photorespiratory mutant lacking plastidic glutamine synthetase; c) a clustering analysis indicated a similar pattern of expression among several genes involved in photorespiratory and asparagine metabolism, indicating a clear link between LjASN1 and LjSGAT genes and photorespiration. Conclusions: The results obt

Background: Asparagine is a very important nitrogen transport and storage compound in plants due to its high nitrogen/carbon ratio and stability. Asparagine intracellu lar concentration depends on a balance between asparagine biosynthesis and degradation. The main enzymes involved in asparagine metabolism are as paragine synthetase (ASN), asparaginase (NSE) and serine-glyoxylate aminotransfera se (SGAT). The study of the genes encoding for these enzymes in the model legume Lotus japonicus is of particular interest since it has been proposed that asparagine is the principal molecule used to transport reduced nitrogen within the plant in most temperate legumes. Results: A differential expression of genes encoding for seve ral enzymes involved in asparagine metabolism was detected in L. japonicus . ASN is encoded by three genes, LjASN1 was the most highly expressed in mature leaves while LjASN2 expression was negligible and LjASN3 showed a low expression in this organ, suggesting that LjASN1 is the main gene responsible for asparagine synthesis in mature leaves. In young leaves, LjASN3 was the only ASN gene expressed although at low levels, while all the three genes encoding for NSE were highly expressed, especially LjNSE1 .Innodules, LjASN2 and LjNSE2 were the most highly expressed genes, suggesting an important role for these genes in this organ. Several lines of evidence support the connection between asparagine metabolic genes and photorespiration in L. japonicus : a) a mutant plant deficient in LjNSE1 showed a dramatic decrease in the expression of the two genes encoding for SGAT; b) expression of the genes involved in asparagine metabolism is altered in a photorespiratory mutant lacking plastidic glutamine synthetase; c) a clustering analysis indicated a similar pattern of expression among several genes involved in photorespiratory and asparagine metabolism, indicating a clear link between LjASN1 and LjSGAT genes and photorespiration. Conclusions: The results obt

[Background]: African swine fever (ASF) is a re-expanding devastating viral disease currently threatening the pig industry worldwide. MicroRNAs are a class of 17–25 nucleotide non- coding RNAs that have been shown to have critical functions in a wide variety of biological processes, such as cell differentiation, cell cycle regulation, carcinogenesis, apoptosis, regulation of immunity as well as in viral infections by cleavage or translational repression of mRNAs. Nevertheless, there is no information about miRNA expression in an ASFV infection., [Methods]: In this proof-of-concept study, we have analyzed miRNAs expressed in spleen and submandibular lymph node of experimentally infected pigs with a virulent (E75) or its derived attenuated (E75CV1) ASFV strain, as well as, at different times post-infection with the virulent strain, by high throughput sequencing of small RNA libraries., [Results]: Spleen presented a more differential expression pattern than lymph nodes in an ASFV infection. Of the most abundant miRNAs, 12 were differentially expressed in both tissues at two different times in infected animals with the virulent strain. Of these, miR-451, miR-145-5p, miR-181a and miR-122 presented up-regulation at late times post-infection while miR-92a, miR-23a, miR-92b-3p, miR-126-5p, miR-126-3p, miR-30d, miR-23b and miR-92c showed down-regulation. Of the 8 differentially expressed miRNAs identified at the same time post-infection in infected animals with the virulent strain compared with animals infected with its attenuated strain, miR-126-5p, miR-92c, miR-92a, miR-30e-5p and miR-500a-5p presented up-regulation whereas miR-125b, miR-451 and miR-125a were down-regulated. All these miRNAs have been shown to be associated with cellular genes involved in pathways related to the immune response, virus-host interactions as well as with several viral genes., [Conclusion]: The study of miRNA expression will contribute to a better understanding of African swine fever virus pathogenesis, essential in the development of any disease control strategy.

[Backgroud] Enteromyxum leei is a myxozoan parasite that produces a slow-progressing intestinal disease. This parasite invades the paracellular space of the intestinal epithelium and progresses from the posterior to the anterior intestine. The aim of the present study was to gain insights into fish T cell responses in the gilthead sea bream-E. leei infection model using a PCR-array with 30 signature molecules for different leukocyte responses in head kidney, spleen, anterior and posterior intestine., [Results] The PCR-array results suggest that E. leei induced migration of T cells from head kidney to intestines where TH1, CTL and TH17 profiles were activated and kept in balance by the upregulation of regulatory cytokines. These results were partially validated by the use of cross-reacting antibodies and BrdU immunostaining to monitor proliferation. Zap70 immunostaining supported the increased number of T cells in the anterior intestine detected by gene expression, but double staining with BrdU did not show active proliferation of this cell type at a local level, supporting the migration from lymphohaematopoietic tissues to the site of infection. Global analyses of the expression profiles revealed a clear separation between infected and exposed, but non-infected fish, more evident in the target organ. Exposed, non-infected animals showed an intermediate phenotype closer to the control fish., [Conclusions] These results evidence a clear modulation of the T cell response of gilthead sea bream upon E. leei infection. The effects occurred both at local and systemic levels, but the response was stronger and more specific at the site of infection, the intestine. Altogether, this research poses a promising basis to understand the response against this important parasite and establish effective preventive or palliative measures.

RNA sequencing studies with complex designs and transcript-resolution analyses involve multiple hypotheses per gene; however, conventional approaches fail to control the false discovery rate (FDR) at gene level. We propose stageR, a two-stage testing paradigm that leverages the increased power of aggregated gene-level tests and allows post hoc assessment for significant genes. This method provides gene-level FDR control and boosts power for testing interaction effects. In transcript-level analysis, it provides a framework that performs powerful gene-level tests while maintaining biological interpretation at transcript-level resolution. The procedure is applicable whenever individual hypotheses can be aggregated, providing a unified framework for complex high-throughput experiments.

BACKGROUND RNA-seq is now widely used to quantitatively assess gene expression, expression differences and isoform switching, and promises to deliver results for the entire transcriptome. However, whether the transcriptional state of a gene can be captured accurately depends critically on library preparation, read alignment, expression estimation and the tests for differential expression and isoform switching. There are comparisons available for the individual steps but there is not yet a systematic investigation which specific genes are impacted by biases throughout the entire analysis workflow. It is especially unclear whether for a given gene, with current methods and protocols, expression changes and isoform switches can be detected. RESULTS For the human genes, we report their detectability under various conditions using different approaches. Overall, we find that the input material has the biggest influence and may, depending on the protocol and RNA degradation, exhibit already strong length-dependent over- and underrepresentation of transcripts. The alignment step aligns for 50% of the isoforms up to 99% of the reads correctly; only in the presence of transcript modifications mainly short isoforms will have a low alignment rate. In our dataset, we found that, depending on the aligner and the input material used, the expression estimation of up to 93% of the genes being accurate within a factor of two; with the deviations being due to ambiguous alignments. Detection of differential expression using a negative-binomial count model works reliably for our simulated data but is dependent on the count accuracy. Interestingly, using the fold-change instead of the p-value as a score for differential expression yields the same performance in the situation of three replicates and the true change being two-fold. Isoform switching is harder to detect and for at least 109 genes the isoform differences evade detection independent of the method used. CONCLUSIONS RNA-seq is

Background: RNA sequencing (RNA-Seq) has emerged as a powerful approach for the detection of differential gene expression with both high-throughput and high resolution capabilities possible depending upon the experimental design chosen. Multiplex experimental designs are now readily available, these can be utilised to increase the numbers of samples or replicates profiled at the cost of decreased sequencing depth generated per sample. These strategies impact on the power of the approach to accurately identify differential expression. This study presents a detailed analysis of the power to detect differential expression in a range of scenarios including simulated null and differential expression distributions with varying numbers of biological or technical replicates, sequencing depths and analysis methods. Results: Differential and non-differential expression datasets were simulated using a combination of negative binomial and exponential distributions derived from real RNA-Seq data. These datasets were used to evaluate the performance of three commonly used differential expression analysis algorithms and to quantify the changes in power with respect to true and false positive rates when simulating variations in sequencing depth, biological replication and multiplex experimental design choices. Conclusions: This work quantitatively explores comparisons between contemporary analysis tools and experimental design choices for the detection of differential expression using RNA-Seq. We found that the DESeq algorithm performs more conservatively than edgeR and NBPSeq. With regard to testing of various experimental designs, this work strongly suggests that greater power is gained through the use of biological replicates relative to library (technical) replicates and sequencing depth. Strikingly, sequencing depth could be reduced as low as 15% without substantial impacts on false positive or true positive rates.

[Background]: There has been an emergence and expansion of tick-borne diseases in Europe, Asia and North America in recent years, including Lyme disease, tick-borne encephalitis and human anaplasmosis. The primary vectors implicated are hard ticks of the genus Ixodes. Although much is known about the host response to these bacterial and viral pathogens, there is limited knowledge of the cellular responses to infection within the tick vector. The bacterium Anaplasma phagocytophilum is able to bypass apoptotic processes in ticks, enabling infection to proceed. However, the tick cellular responses to infection with the flaviviruses tick-borne encephalitis virus (TBEV) and louping ill virus (LIV), which cause tick-borne encephalitis and louping ill respectively, are less clear. [Results]: Infection and transcriptional analysis of the Ixodes ricinus tick cell line IRE/CTVM20 with the viruses LIV and TBEV, and the bacterium A. phagocytophilum, identified activation of common and distinct cellular pathways. In particular, commonly-upregulated genes included those that modulate apoptotic pathways, putative anti-pathogen genes, and genes that influence the tick innate immune response, including selective activation of toll genes. [Conclusion]: These data provide an insight into potential key genes involved in the tick cellular response to viral or bacterial infection, which may promote cell survival and host resistance.

[Background]: There has been an emergence and expansion of tick-borne diseases in Europe, Asia and North America in recent years, including Lyme disease, tick-borne encephalitis and human anaplasmosis. The primary vectors implicated are hard ticks of the genus Ixodes. Although much is known about the host response to these bacterial and viral pathogens, there is limited knowledge of the cellular responses to infection within the tick vector. The bacterium Anaplasma phagocytophilum is able to bypass apoptotic processes in ticks, enabling infection to proceed. However, the tick cellular responses to infection with the flaviviruses tick-borne encephalitis virus (TBEV) and louping ill virus (LIV), which cause tick-borne encephalitis and louping ill respectively, are less clear. [Results]: Infection and transcriptional analysis of the Ixodes ricinus tick cell line IRE/CTVM20 with the viruses LIV and TBEV, and the bacterium A. phagocytophilum, identified activation of common and distinct cellular pathways. In particular, commonly-upregulated genes included those that modulate apoptotic pathways, putative anti-pathogen genes, and genes that influence the tick innate immune response, including selective activation of toll genes. [Conclusion]: These data provide an insight into potential key genes involved in the tick cellular response to viral or bacterial infection, which may promote cell survival and host resistance.

[Background] DNA microarray technology allows the analysis of genome structure and dynamics at genome-wide scale. Expression microarrays (EMA) contain probes for annotated open reading frames (ORF) and are widely used for the analysis of differential gene expression. By contrast, tiling microarrays (TMA) have a much higher probe density and provide unbiased genome-wide coverage. The purpose of this study was to develop a protocol to exploit the high resolution of TMAs for quantitative measurement of DNA strand-specific differential expression of annotated and non-annotated transcripts., [Results] We extensively filtered probes present in Affymetrix Genechip Yeast Genome 2.0 expression and GeneChip S. pombe 1.0FR tiling microarrays to generate custom Chip Description Files (CDF) in order to compare their efficiency. We experimentally tested the potential of our approach by measuring the differential expression of 4904 genes in the yeast Schizosaccharomyces pombe growing under conditions of oxidative stress. The results showed a Pearson correlation coefficient of 0.943 between both platforms, indicating that TMAs are as reliable as EMAs for quantitative expression analysis. A significant advantage of TMAs over EMAs is the possibility of detecting non-annotated transcripts generated only under specific physiological conditions. To take full advantage of this property, we have used a target-labelling protocol that preserves the original polarity of the transcripts and, therefore, allows the strand-specific differential expression of non-annotated transcripts to be determined. By using a segmentation algorithm prior to generating the corresponding custom CDFs, we identified and quantitatively measured the expression of 510 transcripts longer than 180 nucleotides and not overlapping previously annotated ORFs that were differentially expressed at least 2-fold under oxidative stress., [Conclusions] We show that the information derived from TMA hybridization can be processed simultaneously for high-resolution qualitative and quantitative analysis of the differential expression of well-characterized genes and of previously non-annotated and antisense transcripts. The consistency of the performance of TMA, their genome-wide coverage and adaptability to updated genome annotations, and the possibility of measuring strand-specific differential expression makes them a tool of choice for the analysis of gene expression in any organism for which TMA platforms are available.

A response to 'Comprehensive evaluation of differential gene expression analysis methods for RNA-seq data' by Rapaport F, Khanin R, Liang Y, Pirun M, Krek A, Zumbo P, Mason CE, Socci ND and Betel D in Genome Biology, 2013, 14:R95.

Background: Exoskeletal hardening in crustaceans can be attributed to mineralization and sclerotization of the organic matrix. Glycoproteins have been implicated in the calcification process of many matrices. Sclerotization, on the other hand, is catalysed by phenoloxidases, which also play a role in melanization and the immunological response in arthropods. Custom cDNA microarrays from Portunus pelagicus were used to identify genes possibly associated with the activation pathways involved in these processes. Results: Two genes potentially involved in the recognition of glycosylation, the C-type lectin receptor and the mannose-binding protein, were found to display molt cycle-related differential expression profiles. C-type lectin receptor up-regulation was found to coincide with periods associated with new uncalcified cuticle formation, while the up-regulation of mannose-binding protein occurred only in the post-molt stage, during which calcification takes place, implicating both in the regulation of calcification. Genes presumed to be involved in the phenoloxidase activation pathway that facilitates sclerotization also displayed molt cycle-related differential expression profiles. Members of the serine protease superfamily, trypsin-like and chymotrypsin-like, were up-regulated in the intermolt stage when compared to post-molt, while trypsin-like was also up-regulated in pre-molt compared to ecdysis. Additionally, up-regulation in pre- and intermolt stages was observed by transcripts encoding other phenoloxidase activators including the putative antibacterial protein carcinin-like, and clotting protein precursor-like. Furthermore, hemocyanin, itself with phenoloxidase activity, displayed an identical expression pattern to that of the phenoloxidase activators, i.e. up-regulation in pre- and intermolt. Conclusion: Cuticle hardening in crustaceans is a complex process that is precisely timed to occur in the post-molt stage of the molt cycle. We have identified differe

Background: Exoskeletal hardening in crustaceans can be attributed to mineralization and sclerotization of the organic matrix. Glycoproteins have been implicated in the calcification process of many matrices. Sclerotization, on the other hand, is catalysed by phenoloxidases, which also play a role in melanization and the immunological response in arthropods. Custom cDNA microarrays from Portunus pelagicus were used to identify genes possibly associated with the activation pathways involved in these processes. Results: Two genes potentially involved in the recognition of glycosylation, the C-type lectin receptor and the mannose-binding protein, were found to display molt cycle-related differential expression profiles. C-type lectin receptor up-regulation was found to coincide with periods associated with new uncalcified cuticle formation, while the up-regulation of mannose-binding protein occurred only in the post-molt stage, during which calcification takes place, implicating both in the regulation of calcification. Genes presumed to be involved in the phenoloxidase activation pathway that facilitates sclerotization also displayed molt cycle-related differential expression profiles. Members of the serine protease superfamily, trypsin-like and chymotrypsin-like, were up-regulated in the intermolt stage when compared to post-molt, while trypsin-like was also up-regulated in pre-molt compared to ecdysis. Additionally, up-regulation in pre- and intermolt stages was observed by transcripts encoding other phenoloxidase activators including the putative antibacterial protein carcinin-like, and clotting protein precursor-like. Furthermore, hemocyanin, itself with phenoloxidase activity, displayed an identical expression pattern to that of the phenoloxidase activators, i.e. up-regulation in pre- and intermolt. Conclusion: Cuticle hardening in crustaceans is a complex process that is precisely timed to occur in the post-molt stage of the molt cycle. We have identified differe

Background: Exoskeletal hardening in crustaceans can be attributed to mineralization and sclerotization of the organic matrix. Glycoproteins have been implicated in the calcification process of many matrices. Sclerotization, on the other hand, is catalysed by phenoloxidases, which also play a role in melanization and the immunological response in arthropods. Custom cDNA microarrays from Portunus pelagicus were used to identify genes possibly associated with the activation pathways involved in these processes. Results: Two genes potentially involved in the recognition of glycosylation, the C-type lectin receptor and the mannose-binding protein, were found to display molt cycle-related differential expression profiles. C-type lectin receptor up-regulation was found to coincide with periods associated with new uncalcified cuticle formation, while the up-regulation of mannose-binding protein occurred only in the post-molt stage, during which calcification takes place, implicating both in the regulation of calcification. Genes presumed to be involved in the phenoloxidase activation pathway that facilitates sclerotization also displayed molt cycle-related differential expression profiles. Members of the serine protease superfamily, trypsin-like and chymotrypsin-like, were up-regulated in the intermolt stage when compared to post-molt, while trypsin-like was also up-regulated in pre-molt compared to ecdysis. Additionally, up-regulation in pre- and intermolt stages was observed by transcripts encoding other phenoloxidase activators including the putative antibacterial protein carcinin-like, and clotting protein precursor-like. Furthermore, hemocyanin, itself with phenoloxidase activity, displayed an identical expression pattern to that of the phenoloxidase activators, i.e. up-regulation in pre- and intermolt. Conclusion: Cuticle hardening in crustaceans is a complex process that is precisely timed to occur in the post-molt stage of the molt cycle. We have identified differe

[Background]: Neuroblastic tumours (NBTs) represent a heterogeneous spectrum of neoplastic diseases associated with multiple genetic alterations. Structural and numerical chromosomal changes are frequent and are predictive parameters of NBTs outcome. We performed a comparative analysis of the biological entities constituted by NBTs with different ploidy status. [Methods]: Gene expression profiling of 49 diagnostic primary NBTs with ploidy data was performed using oligonucleotide microarray. Further analyses using Quantitative Real-Time Polymerase Chain Reaction (Q-PCR); array-Comparative Genomic Hybridization (aCGH); and Fluorescent in situ Hybridization (FISH) were performed to investigate the correlation between aneuploidy, chromosomal changes and gene expression profiles. [Results]: Gene expression profiling of 49 primary near-triploid and near-diploid/tetraploid NBTs revealed distinct expression profiles associated with each NBT subgroup. A statistically significant portion of genes mapped to 1p36 (P = 0.01) and 17p13-q21 (P < 0.0001), described as recurrently altered in NBTs. Over 90% of these genes showed higher expression in near-triploid NBTs and the majority are involved in cell differentiation pathways. Specific chromosomal abnormalities observed in NBTs, 1p loss, 17q and whole chromosome 17 gains, were reflected in the gene expression profiles. Comparison between gene copy number and expression levels suggests that differential expression might be only partly dependent on gene copy number. Intratumoural clonal heterogeneity was observed in all NBTs, with marked interclonal variability in near-diploid/tetraploid tumours. [Conclusion]: NBTs with different cellular DNA content display distinct transcriptional profiles with a significant portion of differentially expressed genes mapping to specific chromosomal regions known to be associated with outcome. Furthermore, our results demonstrate that these specific genetic abnormalities are highly heterogeneous in

[Background]: One key step in gene expression is the biogenesis of mRNA ribonucleoparticle complexes (mRNPs). Formation of the mRNP requires the participation of a number of conserved factors such as the THO complex. THO interacts physically and functionally with the Sub2/UAP56 RNA-dependent ATPase, and the Yra1/REF1/ALY RNA-binding protein linking transcription, mRNA export and genome integrity. Given the link between genome instability and cancer, we have performed a comparative analysis of the expression patterns of THOC1, a THO complex subunit, and ALY in tumor samples. [Methods]: The mRNA levels were measured by quantitative real-time PCR and hybridization of a tumor tissue cDNA array; and the protein levels and distribution by immunostaining of a custom tissue array containing a set of paraffin-embedded samples of different tumor and normal tissues followed by statistical analysis. [Results]: We show that the expression of two mRNP factors, THOC1 and ALY are altered in several tumor tissues. THOC1 mRNA and protein levels are up-regulated in ovarian and lung tumors and down-regulated in those of testis and skin, whereas ALY is altered in a wide variety of tumors. In contrast to THOC1, ALY protein is highly detected in normal proliferative cells, but poorly in high-grade cancers. [Conclusions]: These results suggest a differential connection between tumorogenesis and the expression levels of human THO and ALY. This study opens the possibility of defining mRNP biogenesis factors as putative players in cell proliferation that could contribute to tumor development.

[Background]: Bovine anaplasmosis, caused by the rickettsial tick-borne pathogen Anaplasma marginale (Rickettsiales: Anaplasmataceae), is vectored by Rhipicephalus (Boophilus) microplus in many tropical and subtropical regions of the world. A. marginale undergoes a complex developmental cycle in ticks which results in infection of salivary glands from where the pathogen is transmitted to cattle. In previous studies, we reported modification of gene expression in Dermacentor variabilis and cultured Ixodes scapularis tick cells in response to infection with A. marginale. In these studies, we extended these findings by use of a functional genomics approach to identify genes differentially expressed in R. microplus male salivary glands in response to A. marginale infection. Additionally, a R. microplus-derived cell line, BME26, was used for the first time to also study tick cell gene expression in response to A. marginale infection., [Results]: Suppression subtractive hybridization libraries were constructed from infected and uninfected ticks and used to identify genes differentially expressed in male R. microplus salivary glands infected with A. marginale. A total of 279 ESTs were identified as candidate differentially expressed genes. Of these, five genes encoding for putative histamine-binding protein (22Hbp), von Willebrand factor (94Will), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz) and proline-rich protein BstNI subfamily 3 precursor (7BstNI3) were confirmed by real-time RT-PCR to be down-regulated in tick salivary glands infected with A. marginale. The impact of selected tick genes on A. marginale infections in tick salivary glands and BME26 cells was characterized by RNA interference. Silencing of the gene encoding for putative flagelliform silk protein (100Silk) resulted in reduced A. marginale infection in both tick salivary glands and cultured BME26 cells, while silencing of the gene encoding for subolesin (4D8) significantly reduced infection only in cultured BME26 cells. The knockdown of the gene encoding for putative metallothionein (93Meth), significantly up-regulated in infected cultured BME26 cells, resulted in higher A. marginale infection levels in tick cells., [Conclusions]: Characterization of differential gene expression in salivary glands of R. microplus in response to A. marginale infection expands our understanding of the molecular mechanisms at the tick-pathogen interface. Functional studies suggested that differentially expressed genes encoding for subolesin, putative von Willebrand factor and flagelliform silk protein could play a role in A. marginale infection and multiplication in ticks. These tick genes found to be functionally relevant for tick-pathogen interactions will likely be candidates for development of vaccines designed for control of both ticks and tick-borne pathogens.

[Background] Parasitic angiosperm Orobanche crenata infection represents a major constraint for the cultivation of legumes worldwide. The level of protection achieved to date is either incomplete or ephemeral. Hence, an efficient control of the parasite requires a better understanding of its interaction and associated resistance mechanisms at molecular levels., [Results] In order to study the plant response to this parasitic plant and the molecular basis of the resistance we have used a proteomic approach. The root proteome of two accessions of the model legume Medicago truncatula displaying differences in their resistance phenotype, in control as well as in inoculated plants, over two time points (21 and 25 days post infection), has been compared. We report quantitative as well as qualitative differences in the 2-DE maps between early- (SA 27774) and late-resistant (SA 4087) genotypes after Coomassie and silver-staining: 69 differential spots were observed between non-inoculated genotypes, and 42 and 25 spots for SA 4087 and SA 27774 non-inoculated and inoculated plants, respectively. In all, 49 differential spots were identified by peptide mass fingerprinting (PMF) following MALDI-TOF/TOF mass spectrometry. Many of the proteins showing significant differences between genotypes and after parasitic infection belong to the functional category of defense and stress-related proteins. A number of spots correspond to proteins with the same function, and might represent members of a multigenic family or post-transcriptional forms of the same protein., [Conclusion] The results obtained suggest the existence of a generic defense mechanism operating during the early stages of infection and differing in both genotypes. The faster response to the infection observed in the SA 27774 genotype might be due to the action of proteins targeted against key elements needed for the parasite's successful infection, such as protease inhibitors. Our data are discussed and compared with those previously obtained with pea [1] and transcriptomic analysis of other plant-pathogen and plant-parasitic plant systems.

To access full text version of this article. Please click on the hyperlink "Full Text" at the bottom of this page, Chemokines and their receptors are essential in the recruitment and positioning of lymphocytes. To address the question of B cell migration into the inflamed synovial tissue of patients with rheumatoid arthritis (RA), peripheral blood naive B cells, memory B cells and plasma cells were analyzed for cell surface expression of the chemokine receptors CXCR3, CXCR4, CXCR5, CCR5, CCR6, CCR7 and CCR9. For comparison, B cells in the peripheral blood of patients with the autoimmune disease systemic lupus erythematosus (SLE) or with the degenerative disease osteoarthritis (OA) were analyzed. Expression levels of chemokine receptors were measured by flow cytometry and were compared between the different patient groups and healthy individuals. The analysis of chemokine receptor expression showed that the majority of peripheral blood B cells is positive for CXCR3, CXCR4, CXCR5, CCR6 and CCR7. Whereas a small fraction of B cells were positive for CCR5, practically no expression of CCR9 was found. In comparison with healthy individuals, in patients with RA a significant fraction of B cells showed a decreased expression of CXCR5 and CCR6 and increased levels of CXCR3. The downregulation of CXCR5 correlated with an upregulation of CXCR3. In patients with SLE, significant changes in CXCR5 expression were seen. The functionality of the chemokine receptors CXCR3 and CXCR4 was demonstrated by transmigration assays with the chemokines CXCL10 and CXCL12, respectively. Our results suggest that chronic inflammation leads to modulation of chemokine receptor expression on peripheral blood B cells. However, differences between patients with RA and patients with SLE point toward a disease-specific regulation of receptor expression. These differences may influence the migrational behavior of B cells.

Background: Sponges (phylum Porifera) constantly interact with microbes. They graze on microbes from the water column by filter-feeding and they harbor symbiotic partners within their bodies. In experimental setups, sponges take up symbionts at lower rates compared with seawater microbes. This suggests that sponges have the capacity to differentiate between microbes and preferentially graze in non-symbiotic microbes, although the underlying mechanisms of discrimination are still poorly understood. Genomic studies showed that, compared to other animal groups, sponges present an extended repertoire of immune receptors, in particular NLRs, SRCRs, and GPCRs, and a handful of experiments showed that sponges regulate the expression of these receptors upon encounter with microbial elicitors. We hypothesize that sponges may rely on differential expression of their diverse repertoire of poriferan immune receptors to sense different microbial consortia while filter-feeding. To test this, we characterized the transcriptomic response of two sponge species, Aplysina aerophoba and Dysidea avara, upon incubation with microbial consortia extracted from A. aerophoba in comparison with incubation with seawater microbes. The sponges were sampled after 1 h, 3 h, and 5 h for RNA-Seq differential gene expression analysis. Results: D. avara incubated with A. aerophoba-symbionts regulated the expression of genes related to immunity, ubiquitination, and signaling. Within the set of differentially-expressed immune genes we identified different families of Nucleotide Oligomerization Domain (NOD)-Like Receptors (NLRs). These results represent the first experimental evidence that different types of NLRs are involved in microbial discrimination in a sponge. In contrast, the transcriptomic response of A. aerophoba to its own symbionts involved comparatively fewer genes and lacked genes encoding for immune receptors. Conclusion: Our work suggests that: (i) the transcriptomic response of sponges up

Background: Sponges (phylum Porifera) constantly interact with microbes. They graze on microbes from the water column by filter-feeding and they harbor symbiotic partners within their bodies. In experimental setups, sponges take up symbionts at lower rates compared with seawater microbes. This suggests that sponges have the capacity to differentiate between microbes and preferentially graze in non-symbiotic microbes, although the underlying mechanisms of discrimination are still poorly understood. Genomic studies showed that, compared to other animal groups, sponges present an extended repertoire of immune receptors, in particular NLRs, SRCRs, and GPCRs, and a handful of experiments showed that sponges regulate the expression of these receptors upon encounter with microbial elicitors. We hypothesize that sponges may rely on differential expression of their diverse repertoire of poriferan immune receptors to sense different microbial consortia while filter-feeding. To test this, we characterized the transcriptomic response of two sponge species, Aplysina aerophoba and Dysidea avara, upon incubation with microbial consortia extracted from A. aerophoba in comparison with incubation with seawater microbes. The sponges were sampled after 1 h, 3 h, and 5 h for RNA-Seq differential gene expression analysis. Results: D. avara incubated with A. aerophoba-symbionts regulated the expression of genes related to immunity, ubiquitination, and signaling. Within the set of differentially-expressed immune genes we identified different families of Nucleotide Oligomerization Domain (NOD)-Like Receptors (NLRs). These results represent the first experimental evidence that different types of NLRs are involved in microbial discrimination in a sponge. In contrast, the transcriptomic response of A. aerophoba to its own symbionts involved comparatively fewer genes and lacked genes encoding for immune receptors. Conclusion: Our work suggests that: (i) the transcriptomic response of sponges up

Background: Sponges (phylum Porifera) constantly interact with microbes. They graze on microbes from the water column by filter-feeding and they harbor symbiotic partners within their bodies. In experimental setups, sponges take up symbionts at lower rates compared with seawater microbes. This suggests that sponges have the capacity to differentiate between microbes and preferentially graze in non-symbiotic microbes, although the underlying mechanisms of discrimination are still poorly understood. Genomic studies showed that, compared to other animal groups, sponges present an extended repertoire of immune receptors, in particular NLRs, SRCRs, and GPCRs, and a handful of experiments showed that sponges regulate the expression of these receptors upon encounter with microbial elicitors. We hypothesize that sponges may rely on differential expression of their diverse repertoire of poriferan immune receptors to sense different microbial consortia while filter-feeding. To test this, we characterized the transcriptomic response of two sponge species, Aplysina aerophoba and Dysidea avara, upon incubation with microbial consortia extracted from A. aerophoba in comparison with incubation with seawater microbes. The sponges were sampled after 1 h, 3 h, and 5 h for RNA-Seq differential gene expression analysis. Results: D. avara incubated with A. aerophoba-symbionts regulated the expression of genes related to immunity, ubiquitination, and signaling. Within the set of differentially-expressed immune genes we identified different families of Nucleotide Oligomerization Domain (NOD)-Like Receptors (NLRs). These results represent the first experimental evidence that different types of NLRs are involved in microbial discrimination in a sponge. In contrast, the transcriptomic response of A. aerophoba to its own symbionts involved comparatively fewer genes and lacked genes encoding for immune receptors. Conclusion: Our work suggests that: (i) the transcriptomic response of sponges up

[Background]: The role of non-coding RNAs in the porcine muscle metabolism is poorly understood, with few studies investigating their expression patterns in response to nutrient supply. Therefore, we aimed to investigate the changes in microRNAs (miRNAs), long intergenic non-coding RNAs (lincRNAs) and mRNAs muscle expression before and after food intake., [Results]: We measured the miRNA, lincRNA and mRNA expression levels in the gluteus medius muscle of 12 gilts in a fasting condition (AL-T0) and 24 gilts fed ad libitum during either 5 h. (AL-T1, N = 12) or 7 h. (AL-T2, N = 12) prior to slaughter. The small RNA fraction was extracted from muscle samples retrieved from the 36 gilts and sequenced, whereas lincRNA and mRNA expression data were already available. In terms of mean and variance, the expression profiles of miRNAs and lincRNAs in the porcine muscle were quite different than those of mRNAs. Food intake induced the differential expression of 149 (AL-T0/AL-T1) and 435 (AL-T0/AL-T2) mRNAs, 6 (AL-T0/AL-T1) and 28 (AL-T0/AL-T2) miRNAs and none lincRNAs, while the number of differentially dispersed genes was much lower. Among the set of differentially expressed miRNAs, we identified ssc-miR-148a-3p, ssc-miR-22-3p and ssc-miR-1, which play key roles in the regulation of glucose and lipid metabolism. Besides, co-expression network analyses revealed several miRNAs that putatively interact with mRNAs playing key metabolic roles and that also showed differential expression before and after feeding. One case example was represented by seven miRNAs (ssc-miR-148a-3p, ssc-miR-151-3p, ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-421-5p, ssc-miR-493-5p and ssc-miR-503) which putatively interact with the PDK4 mRNA, one of the master regulators of glucose utilization and fatty acid oxidation., [Conclusions]: As a whole, our results evidence that microRNAs are likely to play an important role in the porcine skeletal muscle metabolic adaptation to nutrient availability.

Background One strategy for identifying targets of a regulatory factor is to perturb the factor and use high-throughput RNA sequencing to examine the consequences. However, distinguishing direct targets from secondary effects and experimental noise can be challenging when confounding signal is present in the background at varying levels. Results Here, we present a statistical modeling strategy to identify microRNAs that are primary substrates of target-directed miRNA degradation (TDMD) mediated by ZSWIM8. This method uses a bi-beta-uniform mixture (BBUM) model to separate primary from background signal components, leveraging the expectation that primary signal is restricted to upregulation and not downregulation upon loss of ZSWIM8. The BBUM model strategy retained the apparent sensitivity and specificity of the previous ad hoc approach but was more robust against outliers, achieved a more consistent stringency, and could be performed using a single cutoff of false discovery rate (FDR). Conclusions We developed the BBUM model, a robust statistical modeling strategy to account for background secondary signal in differential expression data. It performed well for identifying primary substrates of TDMD and should be useful for other applications in which the primary regulatory targets are only upregulated or only downregulated. The BBUM model, FDR-correction algorithm, and significance-testing methods are available as an R package at https://github.com/wyppeter/bbum .

Background Mesenchymal stromal cells (MSCs) have broad potential as a cell therapy including for the treatment of drug-resistant inflammatory conditions with abnormal T cell proliferation such as graft-versus-host disease (GVHD). Clinical success, however, has been complicated by the heterogeneity of culture-expanded MSCs as well as donor variability. Here, we devise culture conditions that promote expansion of MSCs with enhanced immunomodulatory functions both in vitro and in animal models of GVHD. Methods Human bone marrow-derived MSCs were expanded at high-confluency (MSCHC) and low-confluency state (MSCLC). Their immunomodulatory properties were evaluated with in vitro co-culture assays based on suppression of activated T cell proliferation and secretion of pro-inflammatory cytokines from activated T cells. Metabolic state of these cells was determined, while RNA sequencing was performed to explore transcriptome of these MSCs. Ex vivo expanded MSCHC or MSCLC was injected into human peripheral blood mononuclear cells (PBMC)-induced GVHD mouse model to determine their in vivo therapeutic efficacy based on clinical grade scoring, human CD45+ blood count and histopathological examination. Results As compared to MSCLC, MSCHC significantly reduced both the proliferation of anti-CD3/CD28-activated T cells and secretion of pro-inflammatory cytokines upon MSCHC co-culture across several donors even in the absence of cytokine priming. Mechanistically, metabolic analysis of MSCHC prior to co-culture with activated T cells showed increased glycolytic metabolism and lactate secretion compared to MSCLC, consistent with their ability to inhibit T cell proliferation. Transcriptome analysis further revealed differential expression of immunomodulatory genes including TRIM29, BPIFB4, MMP3 and SPP1 in MSCHC as well as enriched pathways including cytokine–cytokine receptor interactions, cell adhesion and PI3K-AKT signalling. Lastly, we demonstrate in a human PBMC-induced GV

BackgroundMicroRNAs (miRNAs) are short endogenous, single-stranded, noncoding small RNA molecules of approximately 22 nucleotides in length. They regulate gene expression posttranscriptionally by silencing mRNA expression, thus orchestrating many physiological processes. The Small Ruminant Lentiviruses (SRLV) group includes the Visna Maedi Virus (VMV) and Caprine Arthritis Encephalitis (CAEV) viruses, which cause a disease in sheep and goats characterized by pneumonia, mastitis, arthritis and encephalitis. Their main target cells are from the monocyte/macrophage lineage. To date, there are no studies on the role of miRNAs in this viral disease.ResultsUsing RNA-seq technology and bioinformatics analysis, the expression levels of miRNAs during different clinical stages of infection were studied. A total of 212 miRNAs were identified, of which 46 were conserved sequences in other species but found for the first time in sheep, and 12 were completely novel. Differential expression analysis comparing the uninfected and seropositive groups showed changes in several miRNAs; however, no significant differences were detected between seropositive asymptomatic and diseased sheep. The robust increase in the expression level of oar-miR-21 is consistent with its increased expression in other viral diseases. Furthermore, the target prediction of the dysregulated miRNAs revealed that they control genes involved in proliferation-related signalling pathways, such as the PI3K-Akt, AMPK and ErbB pathways.ConclusionsTo the best of our knowledge, this is the first study reporting miRNA profiling in sheep in response to SRLV infection. The known functions of oar-miR-21 as a regulator of inflammation and proliferation appear to be a possible cause of the lesions caused in the sheep's lungs. This miRNA could be an indicator for the severity of the lung lesions, or a putative target for therapeutic intervention.

BackgroundMicroRNAs (miRNAs) are short endogenous, single-stranded, noncoding small RNA molecules of approximately 22 nucleotides in length. They regulate gene expression posttranscriptionally by silencing mRNA expression, thus orchestrating many physiological processes. The Small Ruminant Lentiviruses (SRLV) group includes the Visna Maedi Virus (VMV) and Caprine Arthritis Encephalitis (CAEV) viruses, which cause a disease in sheep and goats characterized by pneumonia, mastitis, arthritis and encephalitis. Their main target cells are from the monocyte/macrophage lineage. To date, there are no studies on the role of miRNAs in this viral disease.ResultsUsing RNA-seq technology and bioinformatics analysis, the expression levels of miRNAs during different clinical stages of infection were studied. A total of 212 miRNAs were identified, of which 46 were conserved sequences in other species but found for the first time in sheep, and 12 were completely novel. Differential expression analysis comparing the uninfected and seropositive groups showed changes in several miRNAs; however, no significant differences were detected between seropositive asymptomatic and diseased sheep. The robust increase in the expression level of oar-miR-21 is consistent with its increased expression in other viral diseases. Furthermore, the target prediction of the dysregulated miRNAs revealed that they control genes involved in proliferation-related signalling pathways, such as the PI3K-Akt, AMPK and ErbB pathways.ConclusionsTo the best of our knowledge, this is the first study reporting miRNA profiling in sheep in response to SRLV infection. The known functions of oar-miR-21 as a regulator of inflammation and proliferation appear to be a possible cause of the lesions caused in the sheep's lungs. This miRNA could be an indicator for the severity of the lung lesions, or a putative target for therapeutic intervention.

BackgroundMicroRNAs (miRNAs) are short endogenous, single-stranded, noncoding small RNA molecules of approximately 22 nucleotides in length. They regulate gene expression posttranscriptionally by silencing mRNA expression, thus orchestrating many physiological processes. The Small Ruminant Lentiviruses (SRLV) group includes the Visna Maedi Virus (VMV) and Caprine Arthritis Encephalitis (CAEV) viruses, which cause a disease in sheep and goats characterized by pneumonia, mastitis, arthritis and encephalitis. Their main target cells are from the monocyte/macrophage lineage. To date, there are no studies on the role of miRNAs in this viral disease.ResultsUsing RNA-seq technology and bioinformatics analysis, the expression levels of miRNAs during different clinical stages of infection were studied. A total of 212 miRNAs were identified, of which 46 were conserved sequences in other species but found for the first time in sheep, and 12 were completely novel. Differential expression analysis comparing the uninfected and seropositive groups showed changes in several miRNAs; however, no significant differences were detected between seropositive asymptomatic and diseased sheep. The robust increase in the expression level of oar-miR-21 is consistent with its increased expression in other viral diseases. Furthermore, the target prediction of the dysregulated miRNAs revealed that they control genes involved in proliferation-related signalling pathways, such as the PI3K-Akt, AMPK and ErbB pathways.ConclusionsTo the best of our knowledge, this is the first study reporting miRNA profiling in sheep in response to SRLV infection. The known functions of oar-miR-21 as a regulator of inflammation and proliferation appear to be a possible cause of the lesions caused in the sheep's lungs. This miRNA could be an indicator for the severity of the lung lesions, or a putative target for therapeutic intervention.

BackgroundMicroRNAs (miRNAs) are short endogenous, single-stranded, noncoding small RNA molecules of approximately 22 nucleotides in length. They regulate gene expression posttranscriptionally by silencing mRNA expression, thus orchestrating many physiological processes. The Small Ruminant Lentiviruses (SRLV) group includes the Visna Maedi Virus (VMV) and Caprine Arthritis Encephalitis (CAEV) viruses, which cause a disease in sheep and goats characterized by pneumonia, mastitis, arthritis and encephalitis. Their main target cells are from the monocyte/macrophage lineage. To date, there are no studies on the role of miRNAs in this viral disease.ResultsUsing RNA-seq technology and bioinformatics analysis, the expression levels of miRNAs during different clinical stages of infection were studied. A total of 212 miRNAs were identified, of which 46 were conserved sequences in other species but found for the first time in sheep, and 12 were completely novel. Differential expression analysis comparing the uninfected and seropositive groups showed changes in several miRNAs; however, no significant differences were detected between seropositive asymptomatic and diseased sheep. The robust increase in the expression level of oar-miR-21 is consistent with its increased expression in other viral diseases. Furthermore, the target prediction of the dysregulated miRNAs revealed that they control genes involved in proliferation-related signalling pathways, such as the PI3K-Akt, AMPK and ErbB pathways.ConclusionsTo the best of our knowledge, this is the first study reporting miRNA profiling in sheep in response to SRLV infection. The known functions of oar-miR-21 as a regulator of inflammation and proliferation appear to be a possible cause of the lesions caused in the sheep's lungs. This miRNA could be an indicator for the severity of the lung lesions, or a putative target for therapeutic intervention.

[Background]: RNA sequencing is a widely used technology for differential expression analysis. However, the RNA-Seq do not provide accurate absolute measurements and the results can be different for each pipeline used. The major problem in statistical analysis of RNA-Seq and in the omics data in general, is the small sample size with respect to the large number of variables. In addition, experimental design must be taken into account and few tools consider it. [Results]: We propose OMICfpp, a method for the statistical analysis of RNA-Seq paired design data. First, we obtain a p-value for each case-control pair using a binomial test. These p-values are aggregated using an ordered weighted average (OWA) with a given orness previously chosen. The aggregated p-value from the original data is compared with the aggregated p-value obtained using the same method applied to random pairs. These new pairs are generated using between-pairs and complete randomization distributions. This randomization p-value is used as a raw p-value to test the differential expression of each gene. The OMICfpp method is evaluated using public data sets of 68 sample pairs from patients with colorectal cancer. We validate our results through bibliographic search of the reported genes and using simulated data set. Furthermore, we compared our results with those obtained by the methods edgeR and DESeq2 for paired samples. Finally, we propose new target genes to validate these as gene expression signatures in colorectal cancer. OMICfpp is available at http://www.uv.es/ayala/software/OMICfpp_0.2.tar.gz. [Conclusions]: Our study shows that OMICfpp is an accurate method for differential expression analysis in RNA-Seq data with paired design. In addition, we propose the use of randomized p-values pattern graphic as a powerful and robust method to select the target genes for experimental validation.

[Background] MicroRNAs (miRNAs) are short endogenous, single-stranded, noncoding small RNA molecules of approximately 22 nucleotides in length. They regulate gene expression posttranscriptionally by silencing mRNA expression, thus orchestrating many physiological processes. The Small Ruminant Lentiviruses (SRLV) group includes the Visna Maedi Virus (VMV) and Caprine Arthritis Encephalitis (CAEV) viruses, which cause a disease in sheep and goats characterized by pneumonia, mastitis, arthritis and encephalitis. Their main target cells are from the monocyte/macrophage lineage. To date, there are no studies on the role of miRNAs in this viral disease., [Results] Using RNA-seq technology and bioinformatics analysis, the expression levels of miRNAs during different clinical stages of infection were studied. A total of 212 miRNAs were identified, of which 46 were conserved sequences in other species but found for the first time in sheep, and 12 were completely novel. Differential expression analysis comparing the uninfected and seropositive groups showed changes in several miRNAs; however, no significant differences were detected between seropositive asymptomatic and diseased sheep. The robust increase in the expression level of oar-miR-21 is consistent with its increased expression in other viral diseases. Furthermore, the target prediction of the dysregulated miRNAs revealed that they control genes involved in proliferation-related signalling pathways, such as the PI3K-Akt, AMPK and ErbB pathways., [Conclusions] To the best of our knowledge, this is the first study reporting miRNA profiling in sheep in response to SRLV infection. The known functions of oar-miR-21 as a regulator of inflammation and proliferation appear to be a possible cause of the lesions caused in the sheep’s lungs. This miRNA could be an indicator for the severity of the lung lesions, or a putative target for therapeutic intervention.

© The Author(s), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in BMC Genomics 14 (2013): 412, doi:10.1186/1471-2164-14-412., Sexual reproduction is a widely studied biological process because it is critically important to the genetics, evolution, and ecology of eukaryotes. Despite decades of study on this topic, no comprehensive explanation has been accepted that explains the evolutionary forces underlying its prevalence and persistence in nature. Monogonont rotifers offer a useful system for experimental studies relating to the evolution of sexual reproduction due to their rapid reproductive rate and close relationship to the putatively ancient asexual bdelloid rotifers. However, little is known about the molecular underpinnings of sex in any rotifer species. We generated mRNA-seq libraries for obligate parthenogenetic (OP) and cyclical parthenogenetic (CP) strains of the monogonont rotifer, Brachionus calyciflorus, to identify genes specific to both modes of reproduction. Our differential expression analysis identified receptors with putative roles in signaling pathways responsible for the transition from asexual to sexual reproduction. Differential expression of a specific copy of the duplicated cell cycle regulatory gene CDC20 and specific copies of histone H2A suggest that such duplications may underlie the phenotypic plasticity required for reproductive mode switch in monogononts. We further identified differential expression of genes involved in the formation of resting eggs, a process linked exclusively to sex in this species. Finally, we identified transcripts from the bdelloid rotifer Adineta ricciae that have significant sequence similarity to genes with higher expression in CP strains of B. calyciflorus. Our analysis of global gene expression differences between facultatively sexual and exclusively asexual populations of B. calyciflorus provides insights into the molecular nature of sexual reproduction in rotifers. Furthermore, our results offer insight into the evolution of obligate asexuality in bdelloid rotifers and provide indicators important for the use of monogononts as, This work was funded by National Institutes of Health Institute of General Medical Sciences (grant number 5R01GM079484, to JML and DMW).

© The Author(s), 2013. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in BMC Genomics 14 (2013): 412, doi:10.1186/1471-2164-14-412., Sexual reproduction is a widely studied biological process because it is critically important to the genetics, evolution, and ecology of eukaryotes. Despite decades of study on this topic, no comprehensive explanation has been accepted that explains the evolutionary forces underlying its prevalence and persistence in nature. Monogonont rotifers offer a useful system for experimental studies relating to the evolution of sexual reproduction due to their rapid reproductive rate and close relationship to the putatively ancient asexual bdelloid rotifers. However, little is known about the molecular underpinnings of sex in any rotifer species. We generated mRNA-seq libraries for obligate parthenogenetic (OP) and cyclical parthenogenetic (CP) strains of the monogonont rotifer, Brachionus calyciflorus, to identify genes specific to both modes of reproduction. Our differential expression analysis identified receptors with putative roles in signaling pathways responsible for the transition from asexual to sexual reproduction. Differential expression of a specific copy of the duplicated cell cycle regulatory gene CDC20 and specific copies of histone H2A suggest that such duplications may underlie the phenotypic plasticity required for reproductive mode switch in monogononts. We further identified differential expression of genes involved in the formation of resting eggs, a process linked exclusively to sex in this species. Finally, we identified transcripts from the bdelloid rotifer Adineta ricciae that have significant sequence similarity to genes with higher expression in CP strains of B. calyciflorus. Our analysis of global gene expression differences between facultatively sexual and exclusively asexual populations of B. calyciflorus provides insights into the molecular nature of sexual reproduction in rotifers. Furthermore, our results offer insight into the evolution of obligate asexuality in bdelloid rotifers and provide indicators important for the use of monogononts as, This work was funded by National Institutes of Health Institute of General Medical Sciences (grant number 5R01GM079484, to JML and DMW).

Molecular gene signatures are useful tools to characterize the physiological state of cell populations, but most have developed under a narrow range of conditions and cell types and are often restricted to a set of gene identities. Focusing on the transcriptional response to hypoxia, we aimed to generate widely applicable classifiers sourced from the results of a meta-analysis of 69 differential expression datasets which included 425 individual RNA-seq experiments from 33 different human cell types exposed to different degrees of hypoxia (0.1–5%O2) for 2–48 h. The resulting decision trees include both gene identities and quantitative boundaries, allowing for easy classification of individual samples without control or normoxic reference. Each tree is composed of 3–5 genes mostly drawn from a small set of just 8 genes (EGLN1, MIR210HG, NDRG1, ANKRD37, TCAF2, PFKFB3, BHLHE40, and MAFF). In spite of their simplicity, these classifiers achieve over 95% accuracy in cross validation and over 80% accuracy when applied to additional challenging datasets. Our results indicate that the classifiers are able to identify hypoxic tumor samples from bulk RNAseq and hypoxic regions within tumor from spatially resolved transcriptomics datasets. Moreover, application of the classifiers to histological sections from normal tissues suggest the presence of a hypoxic gene expression pattern in the kidney cortex not observed in other normoxic organs. Finally, tree classifiers described herein outperform traditional hypoxic gene signatures when compared against a wide range of datasets. This work describes a set of hypoxic gene signatures, structured as simple decision tress, that identify hypoxic samples and regions with high accuracy and can be applied to a broad variety of gene expression datasets and formats.

Background: Type VI secretion systems (T6SS) are widespread among Gram-negative bacteria and have a potential role as essential virulence factors or to maintain symbiotic interactions. Three T6SS gene clusters were identified in the genome of E. amylovora CFBP 1430, of which T6SS-1 and T6SS-3 represent complete T6SS machineries, while T6SS-2 is reduced in its gene content. Results: To assess the contribution of T6SSs to virulence and potential transcriptomic changes of E. amylovora CFBP 1430, single and double mutants in two structural genes were generated for T6SS-1 and T6SS-3. Plant assays showed that mutants in T6SS-3 were slightly more virulent in apple shoots while inducing less disease symptoms on apple flowers, indicating that T6SSs have only a minor effect on virulence of E. amylovora CFBP 1430. The mutations led under in vitro conditions to the differential expression of type III secretion systems, iron acquisition, chemotaxis, flagellar, and fimbrial genes. Comparison of the in planta and in vitro transcriptome data sets revealed a common differential expression of three processes and a set of chemotaxis and motility genes. Additional experiments proved that T6SS mutants are impaired in their motility. Conclusion: These results suggest that the deletion of T6SSs alters metabolic and motility processes. Nevertheless, the difference in lesion development in apple shoots and flower necrosis of T6SS mutants was indicative that T6SSs influences the disease progression and the establishment of the pathogen on host plants.

Background: Type VI secretion systems (T6SS) are widespread among Gram-negative bacteria and have a potential role as essential virulence factors or to maintain symbiotic interactions. Three T6SS gene clusters were identified in the genome of E. amylovora CFBP 1430, of which T6SS-1 and T6SS-3 represent complete T6SS machineries, while T6SS-2 is reduced in its gene content. Results: To assess the contribution of T6SSs to virulence and potential transcriptomic changes of E. amylovora CFBP 1430, single and double mutants in two structural genes were generated for T6SS-1 and T6SS-3. Plant assays showed that mutants in T6SS-3 were slightly more virulent in apple shoots while inducing less disease symptoms on apple flowers, indicating that T6SSs have only a minor effect on virulence of E. amylovora CFBP 1430. The mutations led under in vitro conditions to the differential expression of type III secretion systems, iron acquisition, chemotaxis, flagellar, and fimbrial genes. Comparison of the in planta and in vitro transcriptome data sets revealed a common differential expression of three processes and a set of chemotaxis and motility genes. Additional experiments proved that T6SS mutants are impaired in their motility. Conclusion: These results suggest that the deletion of T6SSs alters metabolic and motility processes. Nevertheless, the difference in lesion development in apple shoots and flower necrosis of T6SS mutants was indicative that T6SSs influences the disease progression and the establishment of the pathogen on host plants.

[Background] The Norway lobster, Nephrops norvegicus, is economically important in European fisheries and is a key organism in local marine ecosystems. Despite multi-faceted scientific interest in this species, our current knowledge of genetic resources in this species remains very limited. Here, we generated a reference de novo transcriptome for N. norvegicus from multiple tissues in both sexes. Bioinformatic analyses were conducted to detect transcripts that were expressed exclusively in either males or females. Patterns were validated via RT-PCR., [Results] Sixteen N. norvegicus libraries were sequenced from immature and mature ovary, testis and vas deferens (including the masculinizing androgenic gland). In addition, eyestalk, brain, thoracic ganglia and hepatopancreas tissues were screened in males and both immature and mature females. RNA-Sequencing resulted in >600 million reads. De novo assembly that combined the current dataset with two previously published libraries from eyestalk tissue, yielded a reference transcriptome of 333,225 transcripts with an average size of 708 base pairs (bp), with an N50 of 1272 bp. Sex-specific transcripts were detected primarily in gonads followed by hepatopancreas, brain, thoracic ganglia, and eyestalk, respectively. Candidate transcripts that were expressed exclusively either in males or females were highlighted and the 10 most abundant ones were validated via RT-PCR. Among the most highly expressed genes were Serine threonine protein kinase in testis and Vitellogenin in female hepatopancreas. These results align closely with gene annotation results. Moreover, a differential expression heatmap showed that the majority of differentially expressed transcripts were identified in gonad and eyestalk tissues. Results indicate that sex-specific gene expression patterns in Norway lobster are controlled by differences in gene regulation pattern between males and females in somatic tissues., [Conclusions] The current study presents the first multi-tissue reference transcriptome for the Norway lobster that can be applied to future biological, wild restocking and fisheries studies. Sex-specific markers were mainly expressed in males implying that males may experience stronger selection than females. It is apparent that differential expression is due to sex-specific gene regulatory pathways that are present in somatic tissues and not from effects of genes located on heterogametic sex chromosomes. The N. norvegicus data provide a foundation for future gene-based reproductive studies.