Search

Here we use a top-down and bottom-up approach in landscape ecology to analyze the active microbes processing methane fluxes (FCH4) in seasonally flooded-forest (FOR) and -traditional farming systems (TFS) in Amazonian floodplains flooded with black, white, and clear water. Our results revealed higher CH4 emissions from water-atmosphere interface in clear water floodplain, followed by black and white water floodplain, respectively. Active methanogenic and methanotrophic taxa were ubiquitous at 0–15 and 15–30 cm soil layer in FOR and TFS, with differences among the water types with respect to the richness, evenness and diversity of the methanogenic communities. These ecological results were not generalizable regarding to FOR and TFS sites, soil layers, and non-flooded and flooded periods. Despite the predominant oxidation of CH4 in the non-flooded period, higher richness and diversity of methanotrophs were revealed for FOR and TFS in the flooded period. In turn, the structure of the methanogenic and methanotrophic communities and their variation were influenced mainly by soil physicochemical factors, water type, soil depth and the presence of nitrifiers, as Nitrososphaera and Nitrospira. Our study reveals a signature across methanotrophic communities in soils from Amazon floodplain with different water types, with a putative disproportionate role of NC10 phylum in CH4 mitigation in natural and agricultural Amazonian floodplains. These findings open the possibilities to explore the role of NC10 phylum in the carbon cycling in Amazon. [ABSTRACT FROM AUTHOR]

Disrupted circadian rhythms is a prominent and early feature of neurodegenerative diseases including Huntington’s disease (HD). In HD patients and animal models, striatal and hypothalamic neurons expressing molecular circadian clocks are targets of mutant Huntingtin (mHtt) pathogenicity. Yet how mHtt disrupts circadian rhythms remains unclear. In a genetic screen for modifiers of mHtt effects on circadian behavior in Drosophila, we discovered a role for the neurodegenerative disease gene Ataxin2 (Atx2). Genetic manipulations of Atx2 modify the impact of mHtt on circadian behavior as well as mHtt aggregation and demonstrate a role for Atx2 in promoting mHtt aggregation as well as mHtt-mediated neuronal dysfunction. RNAi knockdown of the Fragile X mental retardation gene, dfmr1, an Atx2 partner, also partially suppresses mHtt effects and Atx2 effects depend on dfmr1. Atx2 knockdown reduces the cAMP response binding protein A (CrebA) transcript at dawn. CrebA transcript level shows a prominent diurnal regulation in clock neurons. Loss of CrebA also partially suppresses mHtt effects on behavior and cell loss and restoration of CrebA can suppress Atx2 effects. Our results indicate a prominent role of Atx2 in mediating mHtt pathology, specifically via its regulation of CrebA, defining a novel molecular pathway in HD pathogenesis. [ABSTRACT FROM AUTHOR]

Development of a robust automated platform for enrichment of extracellular vesicles from low sample volume that matches the needs for next-generation sequencing could remove major hurdles for genomic biomarker discovery. Here, we document a protocol for urinary EVs enrichment by utilizing an automated microfluidic system, termed acoustic trap, followed by next-generation sequencing of microRNAs (miRNAs) for biomarker discovery. Specifically, we compared the sequencing output from two small RNA library preparations, NEXTFlex and CATS, using only 130 pg of input total RNA. The samples prepared using NEXTflex was found to contain larger number of unique miRNAs that was the predominant RNA species whereas rRNA was the dominant RNA species in CATS prepared samples. A strong correlation was found between the miRNA expressions of the acoustic trap technical replicate in the NEXTFlex prepared samples, as well as between the acoustic trap and ultracentrifugation enrichment methods. Together, these results demonstrate a robust and automated strategy for biomarker discovery from small volumes of urine. [ABSTRACT FROM AUTHOR]

Virus discovery based on high-throughput sequencing relies on enrichment for virus sequences prior to library preparation to achieve a sufficient number of viral reads. In general, preparations of double-stranded RNA or total RNA preparations treated to remove rRNA are used for sequence enrichment. We used virus-specific antibodies to immunocapture virions from plant sap to conduct cDNA synthesis, followed by library preparation and HTS. For the four potato viruses PLRV, PVY, PVA and PYV, template preparation by virion immunocapture provided a simpler and less expensive method than the enrichment of total RNA by ribosomal depletion. Specific enrichment of viral sequences without an intermediate amplification step was achieved, and this high coverage of sequences across the viral genomes was important to identify rare sequence variations. Using this approach, the first complete genome sequence of a potato yellowing virus isolate (PYV, DSMZ PV-0706) was determined in this study. PYV can be confidently assigned as a distinct species in the genus Ilarvirus. [ABSTRACT FROM AUTHOR]

Formalin-fixed paraffin-embedded (FFPE) tissues are among the most widely available clinical specimens. Their potential utility as a source of RNA for transcriptome studies would greatly enhance population-based cancer studies. Although preliminary studies suggest FFPE tissue may be used for RNA sequencing, the effect of storage time on these specimens needs to be determined. We conducted this study to determine whether RNA in archived FFPE high-grade ovarian serous adenocarcinomas from Surveillance, Epidemiology and End Results (SEER) registries was present in sufficient quantity and quality for RNA-Seq analysis. FFPE tissues, stored from 7 to 32 years, were obtained from three SEER sites. RNA was extracted, quantified, quality assessed, and subjected to RNA-Seq (a whole transcriptome sequencing technology). FFPE specimens stored for longer periods of time had poorer RNA sample quality as indicated by negative correlations between specimen storage time and fragment distribution values (DV). In addition, sample contamination was a common issue among the RNA, with 41 of 67 samples having 5% to 48% bacterial contamination. However, regardless of specimen storage time and bacterial contamination, 60% of the samples yielded data that enabled gene expression quantification, identifying more than 10,000 genes, with the correlations among most biological replicates above 0.7. This study demonstrates that FFPE high-grade ovarian serous adenocarcinomas specimens stored in repositories for up to 32 years and under varying storage conditions are a promising source of RNA for RNA-Seq. We also describe certain caveats to be considered when designing RNA-Seq studies using archived FFPE tissues. [ABSTRACT FROM AUTHOR]

A fundamental assumption, common to the vast majority of high-throughput transcriptome analyses, is that the expression of most genes is unchanged among samples and that total cellular RNA remains constant. As the number of analyzed experimental systems increases however, different independent studies demonstrate that this assumption is often violated. We present a calibration method using RNA spike-ins that allows for the measurement of absolute cellular abundance of RNA molecules. We apply the method to pooled RNA from cell populations of known sizes. For each transcript, we compute a nominal abundance that can be converted to absolute by dividing by a scale factor determined in separate experiments: the yield coefficient of the transcript relative to that of a reference spike-in measured with the same protocol. The method is derived by maximum likelihood theory in the context is a complete statistical model for sequencing counts contributed by cellular RNA and spike-ins. The counts are based on a sample from a fixed number of cells to which a fixed population of spike-in molecules has been added. We illustrate and evaluate the method with applications to two global expression data sets, one from the model eukaryote Saccharomyces cerevisiae, proliferating at different growth rates, and differentiating cardiopharyngeal cell lineages in the chordate Ciona robusta. We tested the method in a technical replicate dilution study, and in a k-fold validation study. [ABSTRACT FROM AUTHOR]

RNA interference (RNAi) requires RNA-dependent RNA polymerases (RdRPs) in many eukaryotes, and RNAi amplification constitutes the only known function for eukaryotic RdRPs. Yet in animals, classical model organisms can elicit RNAi without possessing RdRPs, and only nematode RNAi was shown to require RdRPs. Here we show that RdRP genes are much more common in animals than previously thought, even in insects, where they had been assumed not to exist. RdRP genes were present in the ancestors of numerous clades, and they were subsequently lost at a high frequency. In order to probe the function of RdRPs in a deuterostome (the cephalochordate Branchiostoma lanceolatum), we performed high-throughput analyses of small RNAs from various Branchiostoma developmental stages. Our results show that Branchiostoma RdRPs do not appear to participate in RNAi: we did not detect any candidate small RNA population exhibiting classical siRNA length or sequence features. Our results show that RdRPs have been independently lost in dozens of animal clades, and even in a clade where they have been conserved (cephalochordates) their function in RNAi amplification is not preserved. Such a dramatic functional variability reveals an unexpected plasticity in RNA silencing pathways. [ABSTRACT FROM AUTHOR]

Current progress in the field of next-generation transcriptome sequencing have contributed significantly to the study of various malignancies including glioblastoma multiforme (GBM). Differential sequencing of transcriptomes of patients and non-tumor controls has a potential to reveal novel transcripts with significant role in GBM. One such candidate group of molecules are long non-coding RNAs (lncRNAs) which have been proved to be involved in processes such as carcinogenesis, epigenetic modifications and resistance to various therapeutic approaches. To maximize the value of transcriptome sequencing, a proper protocol for library preparation from tissue-derived RNA needs to be found which would produce high quality transcriptome sequencing data and increase the number of detected lncRNAs. It is important to mention that success of library preparation is determined by the quality of input RNA, which is in case of real-life tissue specimens very often altered in comparison to high quality RNA commonly used by manufacturers for development of library preparation chemistry. In the present study, we used GBM and non-tumor brain tissue specimens and compared three different commercial library preparation kits, namely NEXTflex Rapid Directional qRNA-Seq Kit (Bioo Scientific), SENSE Total RNA-Seq Library Prep Kit (Lexogen) and NEBNext Ultra II Directional RNA Library Prep Kit for Illumina (NEB). Libraries generated using SENSE kit were characterized by the most normal distribution of normalized average GC content, the least amount of over-represented sequences and the percentage of ribosomal RNA reads (0.3–1.5%) and highest numbers of uniquely mapped reads and reads aligning to coding regions. However, NEBNext kit performed better having relatively low duplication rates, even transcript coverage and the highest number of hits in Ensembl database for every biotype of our interest including lncRNAs. Our results indicate that out of three approaches the NEBNext library preparation kit was most suitable for the study of lncRNAs via transcriptome sequencing. This was further confirmed by highly consistent data reached in an independent validation on an expanded cohort. [ABSTRACT FROM AUTHOR]

Preparation of high-quality sequencing libraries is a costly and time-consuming component of metagenomic next generation sequencing (mNGS). While the overall cost of sequencing has dropped significantly over recent years, the reagents needed to prepare sequencing samples are likely to become the dominant expense in the process. Furthermore, libraries prepared by hand are subject to human variability and needless waste due to limitations of manual pipetting volumes. Reduction of reaction volumes, combined with sub-microliter automated dispensing of reagents without consumable pipette tips, has the potential to provide significant advantages. Here, we describe the integration of several instruments, including the Labcyte Echo 525 acoustic liquid handler and the iSeq and NovaSeq Illumina sequencing platforms, to miniaturize and automate mNGS library preparation, significantly reducing the cost and the time required to prepare samples. Through the use of External RNA Controls Consortium (ERCC) spike-in RNAs, we demonstrated the fidelity of the miniaturized preparation to be equivalent to full volume reactions. Furthermore, detection of viral and microbial species from cell culture and patient samples was also maintained in the miniaturized libraries. For 384-well mNGS library preparations, we achieved cost savings of over 80% in materials and reagents alone, and reduced preparation time by 90% compared to manual approaches, without compromising quality or representation within the library. [ABSTRACT FROM AUTHOR]

The binding of transcription factors to short recognition sequences plays a pivotal role in controlling the expression of genes. The sequence and shape characteristics of binding sites influence DNA binding specificity and have also been implicated in modulating the activity of transcription factors downstream of binding. To quantitatively assess the transcriptional activity of tens of thousands of designed synthetic sites in parallel, we developed a synthetic version of STARR-seq (synSTARR-seq). We used the approach to systematically analyze how variations in the recognition sequence of the glucocorticoid receptor (GR) affect transcriptional regulation. Our approach resulted in the identification of a novel highly active functional GR binding sequence and revealed that sequence variation both within and flanking GR’s core binding site can modulate GR activity without apparent changes in DNA binding affinity. Notably, we found that the sequence composition of variants with similar activity profiles was highly diverse. In contrast, groups of variants with similar activity profiles showed specific DNA shape characteristics indicating that DNA shape may be a better predictor of activity than DNA sequence. Finally, using single cell experiments with individual enhancer variants, we obtained clues indicating that the architecture of the response element can independently tune expression mean and cell-to cell variability in gene expression (noise). Together, our studies establish synSTARR as a powerful method to systematically study how DNA sequence and shape modulate transcriptional output and noise. [ABSTRACT FROM AUTHOR]

N6-Methyladenosine (m6A) accounts for approximately 0.2% to 0.6% of all adenosine in mammalian mRNA, representing the most abundant internal mRNA modifications. m6A RNA immunoprecipitation followed by high-throughput sequencing (MeRIP-seq) is a powerful technique to map the m6A location transcriptome-wide. However, this method typically requires 300 μg of total RNA, which limits its application to patient tumors. In this study, we present a refined m6A MeRIP-seq protocol and analysis pipeline that can be applied to profile low-input RNA samples from patient tumors. We optimized the key parameters of m6A MeRIP-seq, including the starting amount of RNA, RNA fragmentation, antibody selection, MeRIP washing/elution conditions, methods for RNA library construction, and the bioinformatics analysis pipeline. With the optimized immunoprecipitation (IP) conditions and a postamplification rRNA depletion strategy, we were able to profile the m6A epitranscriptome using 500 ng of total RNA. We identified approximately 12,000 m6A peaks with a high signal-to-noise (S/N) ratio from 2 lung adenocarcinoma (ADC) patient tumors. Through integrative analysis of the transcriptome, m6A epitranscriptome, and proteome data in the same patient tumors, we identified dynamics at the m6A level that account for the discordance between mRNA and protein levels in these tumors. The refined m6A MeRIP-seq method is suitable for m6A epitranscriptome profiling in a limited amount of patient tumors, setting the ground for unraveling the dynamics of the m6A epitranscriptome and the underlying mechanisms in clinical settings. [ABSTRACT FROM AUTHOR]

Ribosome-binding proteins function broadly in protein synthesis, gene regulation, and cellular homeostasis, but the complete complement of functional ribosome-bound proteins remains unknown. Using quantitative mass spectrometry, we identified late-annotated short open reading frame 2 (Lso2) as a ribosome-associated protein that is broadly conserved in eukaryotes. Genome-wide crosslinking and immunoprecipitation of Lso2 and its human ortholog coiled-coil domain containing 124 (CCDC124) recovered 25S ribosomal RNA in a region near the A site that overlaps the GTPase activation center. Consistent with this location, Lso2 also crosslinked to most tRNAs. Ribosome profiling of yeast lacking LSO2 (lso2Δ) revealed global translation defects during recovery from stationary phase with translation of most genes reduced more than 4-fold. Ribosomes accumulated at start codons, were depleted from stop codons, and showed codon-specific changes in occupancy in lso2Δ. These defects, and the conservation of the specific ribosome-binding activity of Lso2/CCDC124, indicate broadly important functions in translation and physiology. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) are small RNAs that regulate diverse biological processes including multiple aspects of the host-pathogen interface. Consequently, miRNAs are commonly encoded by viruses that undergo long-term persistent infection. Papillomaviruses (PVs) are capable of undergoing persistent infection, but as yet, no widely-accepted PV-encoded miRNAs have been described. The incomplete understanding of PV-encoded miRNAs is due in part to lack of tractable laboratory models for most PV types. To overcome this, we have developed RNA iscovery by forced enome xpression (miDGE), a new wet bench approach to miRNA identification that screens numerous pathogen genomes in parallel. Using miDGE, we screened over 73 different PV genomes for the ability to code for miRNAs. Our results show that most PVs are unlikely to code for miRNAs and we conclusively demonstrate a lack of PV miRNA expression in cancers associated with infections of several high risk HPVs. However, we identified five different high-confidence or highly probable miRNAs encoded by four different PVs (Human PVs 17, 37, 41 and a Fringilla coelebs PV (FcPV1)). Extensive in vitro assays confirm the validity of these miRNAs in cell culture and two FcPV1 miRNAs are further confirmed to be expressed in vivo in a natural host. We show that miRNAs from two PVs (HPV41 & FcPV1) are able to regulate viral transcripts corresponding to the early region of the PV genome. Combined, these findings identify the first canonical PV miRNAs and support that miRNAs of either host or viral origin are important regulators of the PV life cycle. [ABSTRACT FROM AUTHOR]

Pooled CRISPR libraries are widely used in high-throughput screening to study various biological processes. Various pooled CRISPR libraries have been shared for CRISPR screens and useful tools have been developed to construct researcher’s own libraries, however, many researchers are struggling to create their own pooled CRISPR libraries: it is a time-consuming, labor-intensive, and expensive process. In this study, we develop a simple method to customize conventional pooled CRISPR libraries using the CRISPR/Cas9 system. We show that conventional pooled CRISPR libraries can be modified by eliminating gRNAs that target positive genes, enabling the identification of unknown target genes in CRISPR screening. CRISPR/Cas9 system can be applied as a precise tool for customizing conventional pooled CRISPR libraries and will broaden the scope of high-throughput screening technology. [ABSTRACT FROM AUTHOR]

Nucleic acid material of adequate quality is crucial for successful high-throughput sequencing (HTS) analysis. DNA and RNA isolated from archival FFPE material are frequently degraded and not readily amplifiable due to chemical damage introduced during fixation. To identify optimal nucleic acid extraction kits, DNA and RNA quantity, quality and performance in HTS applications were evaluated. DNA and RNA were isolated from five sarcoma archival FFPE blocks, using eight extraction protocols from seven kits from three different commercial vendors. For DNA extraction, the truXTRAC FFPE DNA kit from Covaris gave higher yields and better amplifiable DNA, but all protocols gave comparable HTS library yields using Agilent SureSelect XT and performed well in downstream variant calling. For RNA extraction, all protocols gave comparable yields and amplifiable RNA. However, for fusion gene detection using the Archer FusionPlex Sarcoma Assay, the truXTRAC FFPE RNA kit from Covaris and Agencourt FormaPure kit from Beckman Coulter showed the highest percentage of unique read-pairs, providing higher complexity of HTS data and more frequent detection of recurrent fusion genes. truXTRAC simultaneous DNA and RNA extraction gave similar outputs as individual protocols. These findings show that although successful HTS libraries could be generated in most cases, the different protocols gave variable quantity and quality for FFPE nucleic acid extraction. Selecting the optimal procedure is highly valuable and may generate results in borderline quality specimens. [ABSTRACT FROM AUTHOR]

Exosomes are endosome-derived membrane vesicles that contain proteins, lipids, and nucleic acids. The exosomal transcriptome mediates intercellular communication, and represents an understudied reservoir of novel biomarkers for human diseases. Next-generation sequencing enables complex quantitative characterization of exosomal RNAs from diverse sources. However, detailed protocols describing exosome purification for preparation of exosomal RNA-sequence (RNA-Seq) libraries are lacking. Here we compared methods for isolation of exosomes and extraction of exosomal RNA from human cell-free serum, as well as strategies for attaining equal representation of samples within pooled RNA-Seq libraries. We compared commercial precipitation with ultracentrifugation for exosome purification and confirmed the presence of exosomes via both transmission electron microscopy and immunoblotting. Exosomal RNA extraction was compared using four different RNA purification methods. We determined the minimal starting volume of serum required for exosome preparation and showed that high quality exosomal RNA can be isolated from sera stored for over a decade. Finally, RNA-Seq libraries were successfully prepared with exosomal RNAs extracted from human cell-free serum, cataloguing both coding and non-coding exosomal transcripts. This method provides researchers with strategic options to prepare RNA-Seq libraries and compare RNA-Seq data quantitatively from minimal volumes of fresh and archival human cell-free serum for disease biomarker discovery. [ABSTRACT FROM AUTHOR]

The potato aphid, Macrosiphum euphorbiae, is an important agricultural pest that causes economic losses to potato and tomato production. To establish the transcriptome for this aphid, RNA-Seq libraries constructed from aphids maintained on tomato plants were used in Illumina sequencing generating 52.6 million 75–105 bp paired-end reads. The reads were assembled using Velvet/Oases software with SEED preprocessing resulting in 22,137 contigs with an N50 value of 2,003bp. After removal of contigs from tomato host origin, 20,254 contigs were annotated using BLASTx searches against the non-redundant protein database from the National Center for Biotechnology Information (NCBI) as well as IntereProScan. This identified matches for 74% of the potato aphid contigs. The highest ranking hits for over 12,700 contigs were against the related pea aphid, Acyrthosiphon pisum. Gene Ontology (GO) was used to classify the identified M. euphorbiae contigs into biological process, cellular component and molecular function. Among the contigs, sequences of microbial origin were identified. Sixty five contigs were from the aphid bacterial obligate endosymbiont Buchnera aphidicola origin and two contigs had amino acid similarities to viruses. The latter two were named Macrosiphum euphorbiae virus 2 (MeV-2) and Macrosiphum euphorbiae virus 3 (MeV-3). The highest sequence identity to MeV-2 had the Dysaphis plantaginea densovirus, while to MeV-3 is the Hubei sobemo-like virus 49. Characterization of MeV-2 and MeV-3 indicated that both are transmitted vertically from adult aphids to nymphs. MeV-2 peptides were detected in the aphid saliva and only MeV-2 and not MeV-3 nucleic acids were detected inside tomato leaves exposed to virus-infected aphids. However, MeV-2 nucleic acids did not persist in tomato leaf tissues, after clearing the plants from aphids, indicating that MeV-2 is likely an aphid virus. [ABSTRACT FROM AUTHOR]

Decades of research have demonstrated the crucial importance of viruses in freshwater ecosystems. However, few studies have focused on the seasonal dynamics and potential hosts of RNA viruses. We surveyed microbial-sized (i.e. 5–0.2 μm) mixed community plankton transcriptomes for RNA viral genomes and investigated their distribution between microbial and macrobial plankton over a seasonal cycle across three temperate lakes by quantitative reverse transcriptase PCR (qRT-PCR). A total of 30 contigs bearing similarity to RNA viral genomes were recovered from a global assembly of 30 plankton RNA libraries. Of these, only 13 were found in >2 libraries and recruited >100 reads (of 9.13 x 107 total reads), representing several picornaviruses, two tobamoviruses and a reovirus. We quantified the abundance of four picornaviruses and the reovirus monthly from August 2014 to May 2015. Patterns of viral abundance in the >5 μm size fraction and representation in microbial-sized community RNA libraries over time suggest that one picornavirus genotype (TS24835) and the reovirus (TS148892) may infect small (<5 μm) eukaryotic microorganisms, while two other picornaviruses (TS24641 and TS4340) may infect larger (>5 μm) eukaryotic microorganisms or metazoa. Our data also suggest that picornavirus TS152062 may originate from an allochthonous host. All five viral genotypes were present in at least one size fraction across all 3 lakes during the year, suggesting that RNA viruses may easily disperse between adjacent aquatic habitats. Our data therefore demonstrate that RNA viruses are widespread in temperate lacustrine ecosystems, and may provide evidence of viral infection in larger eukaryotes (including metazoa) inhabiting the lakes. [ABSTRACT FROM AUTHOR]

RNA interference (RNAi) is a complex and highly conserved regulatory mechanism mediated via small RNAs (sRNAs). Recent technical advances in high throughput sequencing have enabled an increasingly detailed analysis of sRNA abundances and profiles in specific body parts and tissues. This enables investigations of the localized roles of microRNAs (miRNAs) and small interfering RNAs (siRNAs). However, variation in the proportions of non-coding RNAs in the samples being compared can hinder these analyses. Specific tissues may vary significantly in the proportions of fragments of longer non-coding RNAs (such as ribosomal RNA or transfer RNA) present, potentially reflecting tissue-specific differences in biological functions. For example, in Drosophila, some tissues contain a highly abundant 30nt rRNA fragment (the 2S rRNA) as well as abundant 5’ and 3’ terminal rRNA fragments. These can pose difficulties for the construction of sRNA libraries as they can swamp the sequencing space and obscure sRNA abundances. Here we addressed this problem and present a modified “rRNA blocking” protocol for the construction of high-definition (HD) adapter sRNA libraries, in D. melanogaster reproductive tissues. The results showed that 2S rRNAs targeted by blocking oligos were reduced from >80% to < 0.01% total reads. In addition, the use of multiple rRNA blocking oligos to bind the most abundant rRNA fragments allowed us to reveal the underlying sRNA populations at increased resolution. Side-by-side comparisons of sequencing libraries of blocked and non-blocked samples revealed that rRNA blocking did not change the miRNA populations present, but instead enhanced their abundances. We suggest that this rRNA blocking procedure offers the potential to improve the in-depth analysis of differentially expressed sRNAs within and across different tissues. [ABSTRACT FROM AUTHOR]

Metagenomics and next generation sequencing were employed to discover new virus natural enemies of the fire ant, Solenopsis invicta Buren in its native range (i.e., Formosa, Argentina) with the ultimate goal of testing and releasing new viral pathogens into U.S. S. invicta populations to provide natural, sustainable control of this ant. RNA was purified from worker ants from 182 S. invicta colonies, which was pooled into 4 groups according to location. A library was created from each group and sequenced using Illumina Miseq technology. After a series of winnowing methods to remove S. invicta genes, known S. invicta virus genes, and all other non-virus gene sequences, 61,944 unique singletons were identified with virus identity. These were assembled de novo yielding 171 contiguous sequences with significant identity to non-plant virus genes. Fifteen contiguous sequences exhibited very high expression rates and were detected in all four gene libraries. One contig (Contig_29) exhibited the highest expression level overall and across all four gene libraries. Random amplification of cDNA ends analyses expanded this contiguous sequence yielding a complete virus genome, which we have provisionally named Solenopsis invicta virus 5 (SINV-5). SINV-5 is a positive-sense, single-stranded RNA virus with genome characteristics consistent with insect-infecting viruses from the family Dicistroviridae. Moreover, the replicative genome strand of SINV-5 was detected in worker ants indicating that S. invicta serves as host for the virus. Many additional sequences were identified that are likely of viral origin. These sequences await further investigation to determine their origins and relationship with S. invicta. This study expands knowledge of the RNA virome diversity found within S. invicta populations. [ABSTRACT FROM AUTHOR]

In contrast to Trypanosoma brucei gambiense and T. b. rhodesiense (the causative agents of human African trypanosomiasis), T. b. brucei is lysed by apolipoprotein-L1 (apoL1)-containing human serum trypanolytic factors (TLF), rendering it non-infectious to humans. While the mechanisms of TLF1 uptake, apoL1 membrane integration, and T. b. gambiense and T. b. rhodesiense apoL1-resistance have been extensively characterised, our understanding of the range of factors that drive apoL1 action in T. b. brucei is limited. Selecting our bloodstream-form T. b. brucei RNAi library with recombinant apoL1 identified an array of factors that supports the trypanocidal action of apoL1, including six putative ubiquitin modifiers and several proteins putatively involved in membrane trafficking; we also identified the known apoL1 sensitivity determinants, TbKIFC1 and the V-ATPase. Most prominent amongst the novel apoL1 sensitivity determinants was a putative ubiquitin ligase. Intriguingly, while loss of this ubiquitin ligase reduces parasite sensitivity to apoL1, its loss enhances parasite sensitivity to TLF1-dominated normal human serum, indicating that free and TLF1-bound apoL1 have contrasting modes-of-action. Indeed, loss of the known human serum sensitivity determinants, p67 (lysosomal associated membrane protein) and the cathepsin-L regulator, ‘inhibitor of cysteine peptidase’, had no effect on sensitivity to free apoL1. Our findings highlight a complex network of proteins that influences apoL1 action, with implications for our understanding of the anti-trypanosomal action of human serum. [ABSTRACT FROM AUTHOR]

microRNAs are small non-coding RNAs that are important regulators of gene expression in a range of animals, including nematodes. We have analysed a cluster of four miRNAs from the pathogenic nematode species Haemonchus contortus that are closely linked in the genome. We find that the cluster is conserved only in clade V parasitic nematodes and in some ascarids, but not in other clade III species nor in clade V free-living nematodes. Members of the cluster are present in parasite excretory-secretory products and can be detected in the abomasum and draining lymph nodes of infected sheep, indicating their release in vitro and in vivo. As observed for other parasitic nematodes, H. contortus adult worms release extracellular vesicles (EV). Small RNA libraries were prepared from vesicle-enriched and vesicle-depleted supernatants from both adult worms and L4 stage larvae. Comparison of the miRNA species in the different fractions indicated that specific miRNAs are packaged within vesicles, while others are more abundant in vesicle-depleted supernatant. Hierarchical clustering analysis indicated that the gut is the likely source of vesicle-associated miRNAs in the L4 stage, but not in the adult worm. These findings add to the growing body of work demonstrating that miRNAs released from parasitic helminths may play an important role in host-parasite interactions. [ABSTRACT FROM AUTHOR]

A noticeable proportion of low transcribed genes involved in wood formation in conifers may have been missed in previous transcriptomic studies. This could be the case for genes related to less abundant cell types, such as axial parenchyma and resin ducts, and genes related to juvenile wood. In this study, two normalized libraries have been obtained from the cambial zone of young individuals of Pinus canariensis C. Sm. ex DC, a species in which such cells are comparatively abundant. These two libraries cover earlywood (EW) and latewood (LW) differentiation, and reads have been de novo meta-assembled into one transcriptome. A high number of previously undescribed genes have been found. The transcriptional profiles during the growing season have been analyzed and several noticeable differences with respect to previous studies have been found. This work contributes to a more complete picture of wood formation in conifers. The genes and their transcription profiles described here provide a useful molecular tool for further studies focused on relevant developmental issues, such as wound response and the formation of traumatic wood, re-sprouting, etc., presumably related to those cells. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules that play important roles in the innate immunity system of invertebrates, especially in the aspect of antivirus. In the present study, high-throughput small RNA Illumina sequencing systems were used to identify differentially expressed miRNAs (DEMs) from the intestines of Procambarus clarkii that were infected with white spot syndrome virus (WSSV). As a result, 39 known and 12 novel miRNAs were identified in both NG and WG small RNA libraries. Seven DEMs were determined to be involved in the antiviral innate immunity in the intestines of P. clarkii. The results of the target gene predictions of the DEMs showed that the putative target genes of these 7 DEMs are related to tight junctions, vascular smooth muscle contraction regulation of the actin cytoskeleton, focal adhesion, RNA transport, mRNA surveillance, viral carcinogenesis, and Salmonella infection. These results provide theoretical insights for future studies on the antiviral immunity of crustaceans. [ABSTRACT FROM AUTHOR]

Cellular permissiveness to HIV infection is highly heterogeneous across individuals. Heterogeneity is also found across CD4+ T cells from the same individual, where only a fraction of cells gets infected. To explore the basis of permissiveness, we performed single-cell RNA-seq analysis of non-infected CD4+ T cells from high and low permissive individuals. Transcriptional heterogeneity translated in a continuum of cell states, driven by T-cell receptor-mediated cell activation and was strongly linked to permissiveness. Proteins expressed at the cell surface and displaying the highest correlation with T cell activation were tested as biomarkers of cellular permissiveness to HIV. FACS sorting using antibodies against several biomarkers of permissiveness led to an increase of HIV cellular infection rates. Top candidate biomarkers included CD25, a canonical activation marker. The combination of CD25 high expression with other candidate biomarkers led to the identification of CD298, CD63 and CD317 as the best biomarkers for permissiveness. CD25highCD298highCD63highCD317high cell population showed an enrichment of HIV-infection of up to 28 fold as compared to the unsorted cell population. The purified hyper-permissive cell subpopulation was characterized by a downregulation of interferon-induced genes and several known restriction factors. Single-cell RNA-seq analysis coupled with functional characterization of cell biomarkers provides signatures of the “HIV-permissive cell”. [ABSTRACT FROM AUTHOR]

A cross between IR64 (high-yielding but drought-susceptible) and Aday Sel (drought-tolerant) rice cultivars yielded a stable line with enhanced grain yield under drought screening field trials at International Rice Research Institute. The major effect qDTY4.1 drought tolerance and yield QTL was detected in the IR77298-14-1-2-10 Backcrossed Inbred Line (BIL) and its IR87705-7-15-B Near Isogenic Line (NIL) with 93.9% genetic similarity to IR64. Although rice yield is extremely susceptible to water stress at reproductive stage, currently, there is only one report on the detection of drought-responsive microRNAs in inflorescence tissue of a Japonica rice line. In this study, more drought-responsive microRNAs were identified in the inflorescence tissues of IR64, IR77298-14-1-2-10 and IR87705-7-15-B via next-generation sequencing. Among the 32 families of inflorescence-specific non-conserved microRNAs that were identified, 22 families were up-regulated in IR87705-7-15-B. Overall 9 conserved and 34 non-conserved microRNA families were found as drought-responsive in rice inflorescence with 5 conserved and 30 non-conserved families induced in the IR87705-7-15-B. The observation of more drought-responsive non-conserved microRNAs may imply their prominence over conserved microRNAs in drought response mechanisms of rice inflorescence. Gene Ontology annotation analysis on the target genes of drought-responsive microRNAs identified in IR87705-7-15-B revealed over-representation of biological processes including development, signalling and response to stimulus. Particularly, four inflorescence-specific microRNAs viz. osa-miR5485, osa-miR5487, osa-miR5492 and osa-miR5517, and two non-inflorescence specific microRNAs viz. osa-miR169d and osa-miR169f.2 target genes that are involved in flower or embryonic development. Among them, osa-miR169d, osa-miR5492 and osa-miR5517 are related to flowering time control. It is also worth mentioning that osa-miR2118 and osa-miR2275, which are implicated in the biosynthesis of rice inflorescence-specific small interfering RNAs, were induced in IR87705-7-15-B but repressed in IR77298-14-1-2-10. Further, gene search within qDTY4.1 QTL region had identified multiple copies of NBS-LRR resistance genes (potential target of osa-miR2118), subtilisins and genes implicated in stomatal movement, ABA metabolism and cuticular wax biosynthesis. [ABSTRACT FROM AUTHOR]

Plant microRNAs (miRNAs) are a class of non-coding RNAs that play important roles in plant development, defense, and symptom development. Here, 547 known miRNAs representing 129 miRNA families, and 282 potential novel miRNAs were identified in Beta macrocarpa using small RNA deep sequencing. A phylogenetic analysis was performed, and 8 Beta lineage-specific miRNAs were identified. Through a differential expression analysis, miRNAs associated with Beet necrotic yellow vein virus (BNYVV) infection were identified and confirmed using a microarray analysis and stem-loop RT-qPCR. In total, 103 known miRNAs representing 38 miRNA families, and 45 potential novel miRNAs were differentially regulated, with at least a two-fold change, in BNYVV-infected plants compared with that of the mock-inoculated control. Targets of these differentially expressed miRNAs were also predicted by degradome sequencing. These differentially expressed miRNAs were involved in hormone biosynthesis and signal transduction pathways, and enhanced axillary bud development and plant defenses. This work is the first to describe miRNAs of the plant genus Beta and may offer a reference for miRNA research in other species in the genus. It provides valuable information on the pathogenicity mechanisms of BNYVV. [ABSTRACT FROM AUTHOR]

RNA-seq-based SSU (small subunit) rRNA (ribosomal RNA) analysis has provided a better understanding of potentially active microbial community within environments. However, for RNA-seq library construction, high quantities of purified RNA are typically required. We propose a modified RNA-seq method for SSU rRNA-based microbial community analysis that depends on the direct ligation of a 5’ adaptor to RNA before reverse-transcription. The method requires only a low-input quantity of RNA (10–100 ng) and does not require a DNA removal step. The method was initially tested on three mock communities synthesized with enriched SSU rRNA of archaeal, bacterial and fungal isolates at different ratios, and was subsequently used for environmental samples of high or low biomass. For high-biomass salt-marsh sediments, enriched SSU rRNA and total nucleic acid-derived RNA-seq datasets revealed highly consistent community compositions for all of the SSU rRNA sequences, and as much as 46.4%-59.5% of 16S rRNA sequences were suitable for OTU (operational taxonomic unit)-based community and diversity analyses with complete coverage of V1-V2 regions. OTU-based community structures for the two datasets were also highly consistent with those determined by all of the 16S rRNA reads. For low-biomass samples, total nucleic acid-derived RNA-seq datasets were analyzed, and highly active bacterial taxa were also identified by the OTU-based method, notably including members of the previously underestimated genus Nitrospira and phylum Acidobacteria in tap water, members of the phylum Actinobacteria on a shower curtain, and members of the phylum Cyanobacteria on leaf surfaces. More than half of the bacterial 16S rRNA sequences covered the complete region of primer 8F, and non-coverage rates as high as 38.7% were obtained for phylum-unclassified sequences, providing many opportunities to identify novel bacterial taxa. This modified RNA-seq method will provide a better snapshot of diverse microbial communities, most notably by OTU-based analysis, even communities with low-biomass samples. [ABSTRACT FROM AUTHOR]

Sheep are a key source of meat, milk and fibre for the global livestock sector, and an important biomedical model. Global analysis of gene expression across multiple tissues has aided genome annotation and supported functional annotation of mammalian genes. We present a large-scale RNA-Seq dataset representing all the major organ systems from adult sheep and from several juvenile, neonatal and prenatal developmental time points. The Ovis aries reference genome (Oar v3.1) includes 27,504 genes (20,921 protein coding), of which 25,350 (19,921 protein coding) had detectable expression in at least one tissue in the sheep gene expression atlas dataset. Network-based cluster analysis of this dataset grouped genes according to their expression pattern. The principle of ‘guilt by association’ was used to infer the function of uncharacterised genes from their co-expression with genes of known function. We describe the overall transcriptional signatures present in the sheep gene expression atlas and assign those signatures, where possible, to specific cell populations or pathways. The findings are related to innate immunity by focusing on clusters with an immune signature, and to the advantages of cross-breeding by examining the patterns of genes exhibiting the greatest expression differences between purebred and crossbred animals. This high-resolution gene expression atlas for sheep is, to our knowledge, the largest transcriptomic dataset from any livestock species to date. It provides a resource to improve the annotation of the current reference genome for sheep, presenting a model transcriptome for ruminants and insight into gene, cell and tissue function at multiple developmental stages. [ABSTRACT FROM AUTHOR]

Albinism, a phenomenon characterized by pigmentation deficiency on the ocular side of Japanese flounder (Paralichthys olivaceus), has caused significant damage. Limited mRNA and microRNA (miRNA) information is available on fish pigmentation deficiency. In this study, a high-throughput sequencing strategy was employed to identify the mRNA and miRNAs involved in P. olivaceus albinism. Based on P. olivaceus genome, RNA-seq identified 21,787 know genes and 711 new genes by transcripts assembly. Of those, 235 genes exhibited significantly different expression pattern (fold change ≥2 or ≤0.5 and q-value≤0.05), including 194 down-regulated genes and 41 up-regulated genes in albino versus normally pigmented individuals. These genes were enriched to 81 GO terms and 9 KEGG pathways (p≤0.05). Among those, the pigmentation related pathways-Melanogenesis and tyrosine metabolism were contained. High-throughput miRNA sequencing identified a total of 475 miRNAs, including 64 novel miRNAs. Furthermore, 33 differentially expressed miRNAs containing 13 up-regulated and 20 down-regulated miRNAs were identified in albino versus normally pigmented individuals (fold change ≥1.5 or ≤0.67 and p≤0.05). The next target prediction discovered a variety of putative target genes, of which, 134 genes including Tyrosinase (TYR), Tyrosinase-related protein 1 (TYRP1), Microphthalmia-associated transcription factor (MITF) were overlapped with differentially expressed genes derived from RNA-seq. These target genes were significantly enriched to 254 GO terms and 103 KEGG pathways (p<0.001). Of those, tyrosine metabolism, lysosomes, phototransduction pathways, etc., attracted considerable attention due to their involvement in regulating skin pigmentation. Expression patterns of differentially expressed mRNA and miRNAs were validated in 10 mRNA and 10 miRNAs by qRT-PCR. With high-throughput mRNA and miRNA sequencing and analysis, a series of interested mRNA and miRNAs involved in fish pigmentation are identified. And the miRNA-mRNA regulatory network also provides a solid starting point for further elucidation of fish pigmentation deficiency. [ABSTRACT FROM AUTHOR]

Curation and storage of formalin-fixed, paraffin-embedded (FFPE) samples are standard procedures in hospital pathology laboratories around the world. Many thousands of such samples exist and could be used for next generation sequencing analysis. Retrospective analyses of such samples are important for identifying molecular correlates of carcinogenesis, treatment history and disease outcomes. Two major hurdles in using FFPE material for sequencing are the damaged nature of the nucleic acids and the labor-intensive nature of nucleic acid purification. These limitations and a number of other issues that span multiple steps from nucleic acid purification to library construction are addressed here. We optimized and automated a 96-well magnetic bead-based extraction protocol that can be scaled to large cohorts and is compatible with automation. Using sets of 32 and 91 individual FFPE samples respectively, we generated libraries from 100 ng of total RNA and DNA starting amounts with 95–100% success rate. The use of the resulting RNA in micro-RNA sequencing was also demonstrated. In addition to offering the potential of scalability and rapid throughput, the yield obtained with lower input requirements makes these methods applicable to clinical samples where tissue abundance is limiting. [ABSTRACT FROM AUTHOR]

Background: Biochemical methods are available for enriching 5′ ends of RNAs in prokaryotes, which are employed in the differential RNA-seq (dRNA-seq) and the more recent Cappable-seq protocols. Computational methods are needed to locate RNA 5′ ends from these data by statistical analysis of the enrichment. Although statistical-based analysis methods have been developed for dRNA-seq, they may not be suitable for Cappable-seq data. The more efficient enrichment method employed in Cappable-seq compared with dRNA-seq could affect data distribution and thus algorithm performance. Results: We present Transformation of Nucleotide Enrichment Ratios (ToNER), a tool for statistical modeling of enrichment from RNA-seq data obtained from enriched and unenriched libraries. The tool calculates nucleotide enrichment scores and determines the global transformation for fitting to the normal distribution using the Box-Cox procedure. From the transformed distribution, sites of significant enrichment are identified. To increase power of detection, meta-analysis across experimental replicates is offered. We tested the tool on Cappable-seq and dRNA-seq data for identifying Escherichia coli transcript 5′ ends and compared the results with those from the TSSAR tool, which is designed for analyzing dRNA-seq data. When combining results across Cappable-seq replicates, ToNER detects more known transcript 5′ ends than TSSAR. In general, the transcript 5′ ends detected by ToNER but not TSSAR occur in regions which cannot be locally modeled by TSSAR. Conclusion: ToNER uses a simple yet robust statistical modeling approach, which can be used for detecting RNA 5′ends from Cappable-seq data, in particular when combining information from experimental replicates. The ToNER tool could potentially be applied for analyzing other RNA-seq datasets in which enrichment for other structural features of RNA is employed. The program is freely available for download at ToNER webpage () and GitHub repository (). [ABSTRACT FROM AUTHOR]

Spermatogenesis is a complex process that involves many elements. However, until now, little is known at the molecular level about spermatogenesis in poultry. Here we investigated microRNAs and their target genes that may be involved in germ cell development and spermatogonial in chicken. We used next-generation sequencing to analyze miRNA expression profiles in three types of germline cells: primordial germ cells (PGCs), spermatogonial stem cells (SSCs), and spermatogonia (Sp) during early stage of spermatogenesis. After validated the candidate miRNAs and corresponding genes’ expression in three types of cells, we found 15 miRNAs that were enriched 21 target genes that may be involved in spermatogenesis. Among the enriched miRNAs, miR-202-5p/3p were up-regulated in the Sp library and down-regulated in the PGCs library. Through RT-qPCR and Dual-Luciferase reporter assay, we confirmed that miR-202-5p bind to LIMK2 and involved in germ cell development. Collectively, we firstly discover the novel miRNAs, like miR-202-5p, and its related genes and pathways, expressed during the early spermatogonial stage in chicken, which will provide new clues for deciphering the molecular mechanism of the miRNAs regulating germline stem cell differentiation and spermatogenesis in chicken. [ABSTRACT FROM AUTHOR]

Antibody libraries are important resources to derive antibodies to be used for a wide range of applications, from structural and functional studies to intracellular protein interference studies to developing new diagnostics and therapeutics. Whatever the goal, the key parameter for an antibody library is its complexity (also known as diversity), i.e. the number of distinct elements in the collection, which directly reflects the probability of finding in the library an antibody against a given antigen, of sufficiently high affinity. Quantitative evaluation of antibody library complexity and quality has been for a long time inadequately addressed, due to the high similarity and length of the sequences of the library. Complexity was usually inferred by the transformation efficiency and tested either by fingerprinting and/or sequencing of a few hundred random library elements. Inferring complexity from such a small sampling is, however, very rudimental and gives limited information about the real diversity, because complexity does not scale linearly with sample size. Next-generation sequencing (NGS) has opened new ways to tackle the antibody library complexity quality assessment. However, much remains to be done to fully exploit the potential of NGS for the quantitative analysis of antibody repertoires and to overcome current limitations. To obtain a more reliable antibody library complexity estimate here we show a new, PCR-free, NGS approach to sequence antibody libraries on Illumina platform, coupled to a new bioinformatic analysis and software (Diversity Estimator of Antibody Library, DEAL) that allows to reliably estimate the complexity, taking in consideration the sequencing error. [ABSTRACT FROM AUTHOR]

Normalization is an essential step with considerable impact on high-throughput RNA sequencing (RNA-seq) data analysis. Although there are numerous methods for read count normalization, it remains a challenge to choose an optimal method due to multiple factors contributing to read count variability that affects the overall sensitivity and specificity. In order to properly determine the most appropriate normalization methods, it is critical to compare the performance and shortcomings of a representative set of normalization routines based on different dataset characteristics. Therefore, we set out to evaluate the performance of the commonly used methods (DESeq, TMM-edgeR, FPKM-CuffDiff, TC, Med UQ and FQ) and two new methods we propose: Med-pgQ2 and UQ-pgQ2 (per-gene normalization after per-sample median or upper-quartile global scaling). Our per-gene normalization approach allows for comparisons between conditions based on similar count levels. Using the benchmark Microarray Quality Control Project (MAQC) and simulated datasets, we performed differential gene expression analysis to evaluate these methods. When evaluating MAQC2 with two replicates, we observed that Med-pgQ2 and UQ-pgQ2 achieved a slightly higher area under the Receiver Operating Characteristic Curve (AUC), a specificity rate > 85%, the detection power > 92% and an actual false discovery rate (FDR) under 0.06 given the nominal FDR (≤0.05). Although the top commonly used methods (DESeq and TMM-edgeR) yield a higher power (>93%) for MAQC2 data, they trade off with a reduced specificity (<70%) and a slightly higher actual FDR than our proposed methods. In addition, the results from an analysis based on the qualitative characteristics of sample distribution for MAQC2 and human breast cancer datasets show that only our gene-wise normalization methods corrected data skewed towards lower read counts. However, when we evaluated MAQC3 with less variation in five replicates, all methods performed similarly. Thus, our proposed Med-pgQ2 and UQ-pgQ2 methods perform slightly better for differential gene analysis of RNA-seq data skewed towards lowly expressed read counts with high variation by improving specificity while maintaining a good detection power with a control of the nominal FDR level. [ABSTRACT FROM AUTHOR]

RNA-sequencing (RNA-seq) has become the standard method for unbiased analysis of gene expression but also provides access to more complex transcriptome features, including alternative RNA splicing, RNA editing, and even detection of fusion transcripts formed through chromosomal translocations. However, differences in library methods can adversely affect the ability to recover these different types of transcriptome data. For example, some methods have bias for one end of transcripts or rely on low-efficiency steps that limit the complexity of the resulting library, making detection of rare transcripts less likely. We tested several commonly used methods of RNA-seq library preparation and found vast differences in the detection of advanced transcriptome features, such as alternatively spliced isoforms and RNA editing sites. By comparing several different protocols available for the Ion Proton sequencer and by utilizing detailed bioinformatics analysis tools, we were able to develop an optimized random primer based RNA-seq technique that is reliable at uncovering rare transcript isoforms and RNA editing features, as well as fusion reads from oncogenic chromosome rearrangements. The combination of optimized libraries and rapid Ion Proton sequencing provides a powerful platform for the transcriptome analysis of research and clinical samples. [ABSTRACT FROM AUTHOR]

Background: A recent study of the gene expression patterns of Zika virus (ZIKV) infected human neural progenitor cells (hNPCs) revealed transcriptional dysregulation and identified cell cycle-related pathways that are affected by infection. However deeper exploration of the information present in the RNA-Seq data can be used to further elucidate the manner in which Zika infection of hNPCs affects the transcriptome, refining pathway predictions and revealing isoform-specific dynamics. Methodology/Principal findings: We analyzed data published by Tang et al. using state-of-the-art tools for transcriptome analysis. By accounting for the experimental design and estimation of technical and inferential variance we were able to pinpoint Zika infection affected pathways that highlight Zika’s neural tropism. The examination of differential genes reveals cases of isoform divergence. Conclusions: Transcriptome analysis of Zika infected hNPCs has the potential to identify the molecular signatures of Zika infected neural cells. These signatures may be useful for diagnostics and for the resolution of infection pathways that can be used to harvest specific targets for further study. [ABSTRACT FROM AUTHOR]

Chagas disease is one of the main parasitic diseases found in Latin America and it is estimated that between six and seven million people are infected worldwide. Its etiologic agent, the protozoan Trypanosoma cruzi, is transmitted by triatomines, some of which from the genus Rhodnius. Twenty species are currently recognized in this genus, including some closely related species with low levels of morphological differentiation, such as Rhodnius montenegrensis and Rhodnius robustus. In order to investigate genetic differences between these two species, we generated large-scale RNA-sequencing data (consisting of four RNA-seq libraries) from the heads and salivary glands of males of R. montenegrensis and R. robustus. Transcriptome assemblies produced for each species resulted in 64,952 contigs for R. montenegrensis and 70,894 contigs for R. robustus, with N50 of approximately 2,100 for both species. SNP calling based on the more complete R. robustus assembly revealed 3,055 fixed interspecific differences and 216 transcripts with high levels of divergence which contained only fixed differences between the two species. A gene ontology enrichment analysis revealed that these highly differentiated transcripts were enriched for eight GO terms related to AP-2 adaptor complex, as well as other interesting genes that could be involved in their differentiation. The results show that R. montenegrensis and R. robustus have a substantial quantity of fixed interspecific polymorphisms, which suggests a high degree of genetic divergence between the two species and likely corroborates the species status of R. montenegrensis. [ABSTRACT FROM AUTHOR]

Lead (Pb) is one of the nonessential and toxic metals that threaten the environment and human health. Medicago sativa L. is a legume with high salt tolerance and high biomass production. It is not only a globally important forage crop but is also an ideal plant for phytoremediation. However, the biological and molecular mechanisms that respond to heavy metals are still not well defined in M. sativa. In this study, de novo and strand-specific RNA-sequencing was performed to identify genes involved in the Pb stress response in M. sativa roots. A total of 415,350 unigenes were obtained from the assembled cDNA libraries, among which 5,416 were identified as significantly differentially expressed genes (DEGs) (false discovery rate < 0.005) between cDNA libraries from control and Pb-treated plants. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs showed they mainly clustered with terms associated with binding, transport, membranes, and the pathways related to signal and energy metabolism. Moreover, a number of candidate genes included antioxidant enzymes, metal transporters, and transcription factors involved in heavy metal response were upregulated under Pb stress. Quantitative real-time PCR(qRT-PCR) validation of the expression patterns of 10 randomly selected candidate DEGs were consistent with the transcriptome analysis results. Thus, this study offers new information towards the investigation of biological changes and molecular mechanisms related to Pb stress response in plants. [ABSTRACT FROM AUTHOR]

Kelps are founding species of temperate marine ecosystems, living in intertidal coastal areas where they are often challenged by generalist and specialist herbivores. As most sessile organisms, kelps develop defensive strategies to restrain grazing damage and preserve their own fitness during interactions with herbivores. To decipher some inducible defense and signaling mechanisms, we carried out metabolome and transcriptome analyses in two emblematic kelp species, Lessonia spicata from South Pacific coasts and Laminaria digitata from North Atlantic, when challenged with their main specialist herbivores. Mass spectrometry based metabolomics revealed large metabolic changes induced in these two brown algae following challenges with their own specialist herbivores. Targeted metabolic profiling of L. spicata further showed that free fatty acid (FFA) and amino acid (AA) metabolisms were particularly regulated under grazing. An early stress response was illustrated by the accumulation of Sulphur containing amino acids in the first twelve hours of herbivory pressure. At latter time periods (after 24 hours), we observed FFA liberation and eicosanoid oxylipins synthesis likely representing metabolites related to stress. Global transcriptomic analysis identified sets of candidate genes specifically induced by grazing in both kelps. qPCR analysis of the top candidate genes during a 48-hours time course validated the results. Most of these genes were particularly activated by herbivore challenge after 24 hours, suggesting that transcriptional reprogramming could be operated at this time period. We demonstrated the potential utility of these genes as molecular markers for herbivory by measuring their inductions in grazed individuals of field harvested L. digitata and L. spicata. By unravelling the regulation of some metabolites and genes following grazing pressure in two kelps representative of the two hemispheres, this work contributes to provide a set of herbivore-induced chemical and molecular responses in kelp species, showing similar inducible responses upon specialist herbivores in their respective ecosystems. [ABSTRACT FROM AUTHOR]

Avoiding biases in next generation sequencing (NGS) library preparation is crucial for obtaining reliable sequencing data. Recently, a new library preparation method has been introduced which has eliminated the need for the ligation step. This method, termed SMART (switching mechanism at the 5′ end of the RNA transcript), is based on template switching reverse transcription. To date, there has been no systematic analysis of the additional biases introduced by this method. We analysed the genomic distribution of sequenced reads prepared from genomic DNA using the SMART methodology and found a strong bias toward long (≥12bp) poly dA/dT containing genomic loci. This bias is unique to the SMART-based library preparation and does not appear when libraries are prepared with conventional ligation based methods. Although this bias is obvious only when performing paired end sequencing, it affects single end sequenced samples as well. Our analysis demonstrates that sequenced reads originating from SMART-DNA libraries are heavily skewed toward genomic poly dA/dT tracts. This bias needs to be considered when deciding to use SMART based technology for library preparation. [ABSTRACT FROM AUTHOR]

The molecular classification of glioblastoma (GBM) based on gene expression might better explain outcome and response to treatment than clinical factors. Whole transcriptome sequencing using next-generation sequencing platforms is rapidly becoming accepted as a tool for measuring gene expression for both research and clinical use. Fresh frozen (FF) tissue specimens of GBM are difficult to obtain since tumor tissue obtained at surgery is often scarce and necrotic and diagnosis is prioritized over freezing. After diagnosis, leftover tissue is usually stored as formalin-fixed paraffin-embedded (FFPE) tissue. However, RNA from FFPE tissues is usually degraded, which could hamper gene expression analysis. We compared RNA-Seq data obtained from matched pairs of FF and FFPE GBM specimens. Only three FFPE out of eleven FFPE-FF matched samples yielded informative results. Several quality-control measurements showed that RNA from FFPE samples was highly degraded but maintained transcriptomic similarities to RNA from FF samples. Certain issues regarding mutation analysis and subtype prediction were detected. Nevertheless, our results suggest that RNA-Seq of FFPE GBM specimens provides reliable gene expression data that can be used in molecular studies of GBM if the RNA is sufficiently preserved. [ABSTRACT FROM AUTHOR]

The presence and relative stability of extracellular RNAs (exRNAs) in biofluids has led to an emerging recognition of their promise as ‘liquid biopsies’ for diseases. Most prior studies on discovery of exRNAs as disease-specific biomarkers have focused on microRNAs (miRNAs) using technologies such as qRT-PCR and microarrays. The recent application of next-generation sequencing to discovery of exRNA biomarkers has revealed the presence of potential novel miRNAs as well as other RNA species such as tRNAs, snoRNAs, piRNAs and lncRNAs in biofluids. At the same time, the use of RNA sequencing for biofluids poses unique challenges, including low amounts of input RNAs, the presence of exRNAs in different compartments with varying degrees of vulnerability to isolation techniques, and the high abundance of specific RNA species (thereby limiting the sensitivity of detection of less abundant species). Moreover, discovery in human diseases often relies on archival biospecimens of varying age and limiting amounts of samples. In this study, we have tested RNA isolation methods to optimize profiling exRNAs by RNA sequencing in individuals without any known diseases. Our findings are consistent with other recent studies that detect microRNAs and ribosomal RNAs as the major exRNA species in plasma. Similar to other recent studies, we found that the landscape of biofluid microRNA transcriptome is dominated by several abundant microRNAs that appear to comprise conserved extracellular miRNAs. There is reasonable correlation of sets of conserved miRNAs across biological replicates, and even across other data sets obtained at different investigative sites. Conversely, the detection of less abundant miRNAs is far more dependent on the exact methodology of RNA isolation and profiling. This study highlights the challenges in detecting and quantifying less abundant plasma miRNAs in health and disease using RNA sequencing platforms. [ABSTRACT FROM AUTHOR]

For most sample types, the automation of RNA and DNA sample preparation workflows enables high throughput next-generation sequencing (NGS) library preparation. Greater adoption of small RNA (sRNA) sequencing has been hindered by high sample input requirements and inherent ligation side products formed during library preparation. These side products, known as adapter dimer, are very similar in size to the tagged library. Most sRNA library preparation strategies thus employ a gel purification step to isolate tagged library from adapter dimer contaminants. At very low sample inputs, adapter dimer side products dominate the reaction and limit the sensitivity of this technique. Here we address the need for improved specificity of sRNA library preparation workflows with a novel library preparation approach that uses modified adapters to suppress adapter dimer formation. This workflow allows for lower sample inputs and elimination of the gel purification step, which in turn allows for an automatable sRNA library preparation protocol. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) play an important role in responding to biotic and abiotic stresses in plants. Jujube witches’-broom a phytoplasma disease of Ziziphus jujuba is prevalent in China and is a serious problem to the industry. However, the molecular mechanism of the disease is poorly understood. In this study, genome-wide identification and analysis of microRNAs in response to witches’-broom was performed. A total of 85 conserved miRNA unique sequences belonging to 32 miRNA families and 24 novel miRNA unique sequences, including their complementary miRNA* strands were identified from small RNA libraries derived from a uninfected and witches’-broom infected Z. jujuba plant. Differentially expressed miRNAs associated with Jujube witches’-broom disease were investigated between the two libraries, and 12 up-regulated miRNAs and 10 down- regulated miRNAs identified with more than 2 fold changes. Additionally, 40 target genes of 85 conserved miRNAs and 49 target genes of 24 novel miRNAs were predicted and their putative functions assigned. Using the modified 5’-RACE method, we confirmed that SPL and MYB were cleaved by miR156 and miR159, respectively. Our results provide insight into the molecular mechanisms of witches’-broom disease in Z. jujuba. [ABSTRACT FROM AUTHOR]

Bombyx mori nucleopolyhedrosis virus (BmNPV) is a major pathogen that threatens the growth and sustainability of the sericulture industry. Since microRNAs (miRNAs) have been shown to play important roles in host-pathogen interactions, in this study we investigated the effects of BmNPV infection on silkworm microRNAs expression profile. To achieve this, we constructed and deep-sequenced two small RNA libraries generated from BmNPV infected and un-infected larvae. The results revealed that 38 silkworm miRNAs were differentially expressed after BmNPV infection. Based on the GO analysis, their predicted target genes were found to be involved in diverse functions such as binding, catalytic, virion and immune response to stimulus suggesting their potential roles in host-virus interactions. Using the dual-luciferase reporter assay, we confirmed that Bmo-miR-277-5p, up-regulated in BmNPV-infected larvae, targeted the B. mori DNA cytosine-5 methyltransferase (Dnmt2) gene which may play potential role in silkworm-BmNPV interaction. These results provide new insights into exploring the interaction mechanism between silkworm and BmNPV. [ABSTRACT FROM AUTHOR]

Diabetic nephropathy (DN) is a major complication of diabetes and is caused by an imbalance in the expression of certain genes that activate or inhibit vital cellular functions of kidney. Despite several recent advances, the pathogenesis of DN remains far from clear, suggesting the need to carry out studies identifying molecular aspects, such as gene expression, that could play a key role in the development of DN. There are several techniques to analyze transcriptome in living organisms. In this study, the suppression subtractive hybridization (SSH) method was used to generate up- and down-regulated subtracted cDNA libraries in the kidney of streptozotocin (STZ)-induced diabetic rats. Northern-blot analysis was used to confirm differential expression ratios from the obtained SSH clones to identify genes related to DN. 400 unique SSH clones were randomly chosen from the two subtraction libraries (200 of each) and verified as differentially expressed. According to blast screening and functional annotation, 20.2% and 20.9% of genes were related to metabolism proteins, 9% and 3.6% to transporters and channels, 16% and 6.3% to transcription factors, 19% and 17.2% to hypothetical proteins, and finally 24.1 and 17.2% to unknown genes, from the down- and up-regulated libraries, respectively. The down- and up-regulated cDNA libraries differentially expressed in the kidney of STZ diabetic rats have been successfully constructed and some identified genes could be highly important in DN.

The objective of this work was to determine the transcript profile of tomato plants (Lycopersicon esculentum Mill.), during Fusarium oxysporum f. sp. lycopersici infection and after foliar application of salicylic acid. The suppression subtractive hybridization (SSH) technique was used to generate a cDNA library enriched for transcripts differentially expressed. A total of 307 clones was identified in two subtractive libraries, which allowed the isolation of several defense-related genes that play roles in different mechanisms of plant resistance to phytopathogens. Genes with unknown roles were also isolated from the two libraries, which indicates the possibility of identifying new genes not yet reported in studies of stress/defense response. The SSH technique is effective for identification of resistance genes activated by salicylic acid and F. oxysporum f. sp. lycopersici infection. Not only the application of this technique enables a cost effective isolation of differentially expressed sequences, but also it allows the identification of novel sequences in tomato from a relative small number of sequences.O objetivo deste trabalho foi determinar o perfil de transcritos em plantas de tomate (Lycopersicon esculentum Mill.), durante a infecção com Fusarium oxysporum f. sp. lycopersici e após a aplicação foliar de ácido salicílico. A técnica de hibridização subtrativa por supressão (SSH) foi utilizada para gerar uma biblioteca de cDNA enriquecida por transcritos diferencialmente expressos. Foram identificados 307 clones, em duas bibliotecas subtrativas, que permitiram o isolamento de diversos genes de defesa com função em diferentes processos relacionados à resistência vegetal contra patógenos. Também foram isolados, nas duas bibliotecas, genes com função desconhecida, o que indica a possibilidade de identificação de novos genes que ainda não tenham sido relatados em estudos anteriores de resposta a estresses e defesa, em plantas. A técnica SSH é eficiente em identificar genes de resistência, ativados pelo ácido salicílico e pela infecção com Fusarium oxysporum f. sp. lycopersici. A aplicação dessa técnica não apenas possibilita isolar seqüências diferencialmente expressas, a baixo custo, como também permite a identificação de novas seqüências, em tomate, a partir de um número relativamente pequeno de seqüências.

With the rapid accumulation of publicly available small RNA sequencing datasets, third-party meta-analysis across many datasets is becoming increasingly powerful. Although removing the 3´ adapter is an essential step for small RNA sequencing analysis, the adapter sequence information is not always available in the metadata. The information can be also erroneous even when it is available. In this study, we developed DNApi, a lightweight Python software package that predicts the 3´ adapter sequence de novo and provides the user with cleansed small RNA sequences ready for down stream analysis. Tested on 539 publicly available small RNA libraries accompanied with 3´ adapter sequences in their metadata, DNApi shows near-perfect accuracy (98.5%) with fast runtime (~2.85 seconds per library) and efficient memory usage (~43 MB on average). In addition to 3´ adapter prediction, it is also important to classify whether the input small RNA libraries were already processed, i.e. the 3´ adapters were removed. DNApi perfectly judged that given another batch of datasets, 192 publicly available processed libraries were “ready-to-map” small RNA sequence. DNApi is compatible with Python 2 and 3, and is available at . The 731 small RNA libraries used for DNApi evaluation were from human tissues and were carefully and manually collected. This study also provides readers with the curated datasets that can be integrated into their studies. [ABSTRACT FROM AUTHOR]