Deep Sequencing Approach Research Articles

Genome-wide identification and prediction of long non-coding RNAs in half-smooth tongue sole Cynoglossus semilaevis

At present, the function and profile of long noncoding RNAs (IncRNAs) in half-smooth tongue soles Cynoglossus semilaevis remain poorly understood. Lherefore, we identified IncRNAs in the species using large-scale deep sequencing approaches. A total of 96 726 222 clean reads and 2497 IncRNAs were obtained from the compound samples. In total, 1 694 known IncRNAs and 803 novel IncRNAs were identified. Most of the novel IncRNAs distributed mainly in a range of 200–3 000 nt and contained 2–6 exons. Six novel IncRNAs were selected for expression analysis by quantitative reverse transcription-polymerase chain reaction, of which lnc_770 and lnc_150 were expressed mainly in the ovary. Lhis study provides a valuable resource for IncRNA studies of half-smooth tongue sole and improves our understanding of the function of non-coding RNA in fish.

World Workshop on Oral Medicine VII: Targeting the oral microbiome Part 2: Current knowledge on malignant and potentially malignant oral disorders

The World Workshop on Oral Medicine VII chose the oral microbiome as a focus area. Part 1 presents the methodological state of the science for oral microbiome studies. Part 2 was guided by the question: What is currently known about the microbiome associated with oral squamous cell carcinoma and potentially malignant disorders of the oral mucosa? A scoping review methodology was followed to identify and analyse relevant studies on the composition and potential functions of the oral microbiota using high-throughput sequencing techniques. The authors performed searches in PubMed and EMBASE. After removal of duplicates, a total of 239 potentially studies were identified. Twenty-three studies on oral squamous cell carcinoma, two on oral leukoplakia and four on oral lichen planus were included with substantial differences in diagnostic criteria, sample type, region sequenced and sequencing method utilised. The majority of studies focused on bacterial identification and recorded statistically significant differences in the oral microbiota associated with health and disease. However, even when comparing studies of similar methodology, the microbial differences between health and disease varied considerably. No consensus on the composition of the microbiomes associated with these conditions on genus and species level could be obtained. Six studies on oral squamous cell carcinoma had included in silico predicted microbial functions (genes and/or pathways) and found some similarities between the studies. Attempts to reveal the microbiome associated with oral mucosal diseases are still in its infancy, and the studies demonstrate significant clinical and methodological heterogeneity across disease categories. The immense richness and diversity of the microbiota clearly illustrate that there is a need for additional methodologically comparable studies utilising deep sequencing approaches in significant cohorts of subjects together with functional analyses. Our hope is that following the recipe as outlined in our preceding companion paper, that is Part 1, will enhance achieving this in the future and elucidate the role of the oral microbiome in oral squamous cell carcinoma and potentially malignant disorders of the oral mucosa.

Comparative profile analysis reveals differentially expressed microRNAs regulate anther and pollen development in kenaf cytoplasmic male sterility line.

Cytoplasmic male sterility (CMS) is advantageous in extensive crop breeding and represents a perfect model for understanding anther and pollen development research. MicroRNAs (miRNAs) play key roles in regulating various biological processes. However, the miRNA-mediated regulatory network in kenaf CMS occurrence remains largely unknown. In the present study, a comparative deep sequencing approach was used to investigate the miRNAs and their roles in regulating anther and pollen development during CMS occurrence. We identified 283 known and 46 new candidate miRNAs in kenaf anther. A total of 67 differentially expressed miRNAs (DEMs) were discovered between CMS and its maintainer line. Among them, 40 and 27 miRNAs were up- and downregulated, respectively. These 67 DEMs were predicted to target 189 genes. Validation of DEMs and putative target genes were confirmed by using real-time quantitative PCR. In addition, a potential miRNA-mediated regulatory network, which mainly involves the auxin signaling pathway, signal transduction, glycolysis and energy metabolism, gene expression, transmembrane transport, protein modification and metabolism, and floral development, that mediates anther development during CMS occurrence was proposed. Taken together, our findings provide a better understanding of the molecular mechanism of miRNA regulation in pollen development and CMS occurrence in kenaf.

Mosquito Small RNA Responses to West Nile and Insect-Specific Virus Infections in Aedes and Culex Mosquito Cells.

Small RNA mediated responses are essential for antiviral defence in mosquitoes, however, they appear to differ per virus-vector combination. To further investigate the diversity of small RNA responses against viruses in mosquitoes, we applied a small RNA deep sequencing approach on five mosquito cell lines: Culex tarsalis CT cells, Aedes albopictus U4.4 and C6/36 cells, Ae. aegypti Aag2 cells (cleared from cell fusing agent virus and Culex Y virus (CYV) by repetitive dsRNA transfections) and Ae. pseudoscutellaris AP-61 cells. De novo assembly of small RNAs revealed the presence of Phasi Charoen-like virus (PCLV), Calbertado virus, Flock House virus and a novel narnavirus in CT cells, CYV in U4.4 cells, and PCLV in Aag2 cells, whereas no insect-specific viruses (ISVs) were detected in C6/36 and AP-61 cells. Next, we investigated the small RNA responses to the identified ISVs and to acute infection with the arthropod-borne West Nile virus (WNV). We demonstrate that AP-61 and C6/36 cells do not produce siRNAs to WNV infection, suggesting that AP-61, like C6/36, are Dicer-2 deficient. CT cells produced a strong siRNA response to the persistent ISVs and acute WNV infection. Interestingly, CT cells also produced viral PIWI-interacting (pi)RNAs to PCLV, but not to WNV or any of the other ISVs. In contrast, in U4.4 and Aag2 cells, WNV siRNAs, and pi-like RNAs without typical ping-pong piRNA signature were observed, while this signature was present in PCLV piRNAs in Aag2 cells. Together, our results demonstrate that mosquito small RNA responses are strongly dependent on both the mosquito cell type and/or the mosquito species and family of the infecting virus.

HIV-1 detection in the olfactory mucosa of HIV-1-infected participants.

HIV infection chronically affects the central nervous system (CNS). Olfactory mucosa is a unique site in the respiratory tract that is directly connected to the CNS; thus we wanted to evaluate olfactory mucosa as a surrogate of CNS sampling. We conducted a preliminary study examining HIV populations and susceptible cells in the olfactory mucosa. Olfactory mucosa was sampled by minimally invasive brushing. Cerebrospinal fluid (CSF) analyses were performed as per routine clinical procedures. Olfactory marker protein, CD4+, CD8+, and trans-activator of transcription (TAT) expressions were assessed by immunohistochemistry. Plasma, CSF, and olfactory mucosa HIV-RNA were quantified using the Cobas AmpliPrep/Cobas TaqMan assay, whereas HIV proviral DNA was evaluated on peripheral blood mononuclear cell and olfactory mucosa. HIV-1 env deep sequencing was performed for phylogenetic analysis. Among ART-naive participants, 88.2% (15/17), and among ART-treated participants, 21.4% (6/28) had detectable HIV-RNA in samples from their olfactory mucosa; CSF escape was more common in patients with olfactory mucosa escape (50 vs. 7.9%; P = 0.010). Olfactory mucosa samples contained few cells positive for CD4, CD8, or HIV-DNA, and no HIV TAT-positive cells, indicating that this approach efficiently samples virions in the olfactory mucosa, but not HIV-infected cells. Yet, using a deep sequencing approach to phylogenetically compare partial HIV env genes in five untreated participants, we identified distinct viral lineages in the OM. The results of this study suggest that nasal brushing is a well tolerated and useful technique for sampling the olfactory mucosa. HIV-RNA was detected in most naïve and in some treated patients, warranting larger longitudinal studies.

Targeted sequencing of a specific gene panel detects a high frequency of ARID1A and PIK3CA mutations in ovarian clear cell carcinoma

BackgroundThe objective of this study was to assess the mutational profile in epithelial ovarian cancer using formalin-fixed, paraffin-embedded (FFPE) tumor specimens from a Taiwanese population by performing targeted sequencing of 9 cancer-associated genes. MethodsTargeted sequencing was performed on 32 formalin-fixed, paraffin embedded (FFPE) tumor specimens, consisting of matched samples from 16 epithelial ovarian cancer patients. Genetic alterations in the 9 cancer-associated genes were detected using a deep sequencing (>1000×) approach. ResultsARID1A and PIK3CA were the most frequently mutated genes. Specifically, ARID1A mutations and PIK3CA mutations were detected in 77.8% and 66.7% of ovarian clear cell carcinoma patients, respectively. Mutations in other genes, including MLH1 (6.3%) and CREBBP (6.3%), were detected in the Taiwanese population. We also identified coexisting ARID1A-PIK3CA mutations (43.8%) and ARID1A-KRAS mutations (12.5%) in tumors. It should also be noted that we identified the presence of three coexisting mutations, the ARID1A-KRAS-PIK3CA mutations and the ARID1A-CREBBP-PIK3CA mutations. ConclusionsIn summary, we identified novel genetic alterations in patients with epithelial ovarian carcinoma (EOC) in a Taiwanese populations. Further studies are needed to elucidate the mechanism of chromatin remodeling to examine the role of the PI3K/AKT pathway, to determine the critical roles of these mechanisms in tumor development and the progression of ovarian malignancy and to investigate new targeted therapies. Overall, our findings were reliable and are worthy of further study.

Molecular Characterization of a Novel Endornavirus Conferring Hypovirulence in Rice Sheath Blight Fungus Rhizoctonia solani AG-1 IA Strain GD-2.

The complete sequence and genome organization of a novel Endornavirus from the hypovirulent strain GD-2 of Rhizoctonia solani AG-1 IA, the causal agent of rice sheath blight, were identified using a deep sequencing approach and it was tentatively named as Rhizoctonia solani endornavirus 1 (RsEV1). It was composed of only one segment that was 19,936 bp in length and was found to be the longest endornavirus genome that has been reported so far. The RsEV1 genome contained two open reading frames (ORFs): ORF1 and ORF2. ORF1 contained a glycosyltransferase 1 domain and a conserved RNA-dependent RNA polymerase domain, whereas ORF2 encoded a conserved hypothetical protein. Phylogenetic analysis revealed that RsEV1 was phylogenetically a new endogenous RNA virus. A horizontal transmission experiment indicated that RsEV1 could be transmitted from the host fungal strain GD-2 to a virulent strain GD-118P and resulted in hypovirulence in the derivative isogenic strain GD-118P-V1. Metabolomic analysis showed that 32 metabolites were differentially expressed between GD-118P and its isogenic hypovirulent strain GD-118P-V1. The differential metabolites were mainly classified as organic acids, amino acids, carbohydrates, and the intermediate products of energy metabolism. Pathway annotation revealed that these 32 metabolites were mainly involved in pentose and glucuronate interconversions and glyoxylate, dicarboxylate, starch, and sucrose metabolism, and so on. Taken together, our results showed that RsEV1 is a novel Endornavirus, and the infection of virulent strain GD-118P by RsEV1 caused metabolic disorders and resulted in hypovirulence. The results of this study lay a foundation for the biocontrol of rice sheath blight caused by R. solani AG1-IA.

Inferring the regulatory network of the miRNA-mediated response to biotic and abiotic stress in melon

BackgroundMiRNAs have emerged as key regulators of stress response in plants, suggesting their potential as candidates for knock-in/out to improve stress tolerance in agricultural crops. Although diverse assays have been performed, systematic and detailed studies of miRNA expression and function during exposure to multiple environments in crops are limited.ResultsHere, we present such pioneering analysis in melon plants in response to seven biotic and abiotic stress conditions. Deep-sequencing and computational approaches have identified twenty-four known miRNAs whose expression was significantly altered under at least one stress condition, observing that down-regulation was preponderant. Additionally, miRNA function was characterized by high scale degradome assays and quantitative RNA measurements over the intended target mRNAs, providing mechanistic insight. Clustering analysis provided evidence that eight miRNAs showed a broad response range under the stress conditions analyzed, whereas another eight miRNAs displayed a narrow response range. Transcription factors were predominantly targeted by stress-responsive miRNAs in melon. Furthermore, our results show that the miRNAs that are down-regulated upon stress predominantly have as targets genes that are known to participate in the stress response by the plant, whereas the miRNAs that are up-regulated control genes linked to development.ConclusionAltogether, this high-resolution analysis of miRNA-target interactions, combining experimental and computational work, Illustrates the close interplay between miRNAs and the response to diverse environmental conditions, in melon.

Gorillas have been infected with the HERV-K (HML-2) endogenous retrovirus much more recently than humans and chimpanzees

Human endogenous retrovirus-K (HERV-K) human mouse mammary tumor virus-like 2 (HML-2) is the most recently active endogenous retrovirus group in humans, and the only group with human-specific proviruses. HML-2 expression is associated with cancer and other diseases, but extensive searches have failed to reveal any replication-competent proviruses in humans. However, HML-2 proviruses are found throughout the catarrhine primates, and it is possible that they continue to infect some species today. To investigate this possibility, we searched for gorilla-specific HML-2 elements using both in silico data mining and targeted deep-sequencing approaches. We identified 150 gorilla-specific integrations, including 31 2-LTR proviruses. Many of these proviruses have identical LTRs, and are insertionally polymorphic, consistent with very recent integration. One identified provirus has full-length ORFs for all genes, and thus could potentially be replication-competent. We suggest that gorillas may still harbor infectious HML-2 virus and could serve as a model for understanding retrovirus evolution and pathogenesis in humans.

Transcriptome-Wide Mapping 5-Methylcytosine by m5C RNA Immunoprecipitation Followed by Deep Sequencing in Plant.

Transcriptome-wide mapping RNA modification is crucial to understand the distribution and function of RNA modifications. Here, we describe a protocol to transcriptome-wide mapping 5-methylcytosine (m5C) in plant, by a RNA immunoprecipitation followed by deep sequencing (m5C-RIP-seq) approach. The procedure includes RNA extraction, fragmentation, RNA immunoprecipitation, and library construction.

Identification of Novel Genes in Osteoarthritic Fibroblast-Like Synoviocytes Using Next-Generation Sequencing and Bioinformatics Approaches.

Synovitis in osteoarthritis (OA) the consequence of low grade inflammatory process caused by cartilage breakdown products that stimulated the production of pro-inflammatory mediators by fibroblast-like synoviocytes (FLS). FLS participate in joint homeostasis and low grade inflammation in the joint microenvironment triggers FLS transformation. In the current study, we aimed to identify differentially expressed genes and potential miRNA regulations in human OA FLS through deep sequencing and bioinformatics approaches. The 245 differentially expressed genes in OA FLS were identified, and pathway analysis using various bioinformatics databases indicated their enrichment in functions related to altered extracellular matrix organization, cell adhesion and cellular movement. Moreover, among the 14 dysregulated genes with potential miRNA regulations identified, src kinase associated phosphoprotein 2 (SKAP2), adaptor related protein complex 1 sigma 2 subunit (AP1S2), PHD finger protein 21A (PHF21A), lipoma preferred partner (LPP), and transcription factor AP-2 alpha (TFAP2A) showed similar expression patterns in OA FLS and OA synovial tissue datasets in Gene Expression Omnibus database. Ingenuity Pathway Analysis identified the dysregulated LPP participated in cell migration and cell spreading of OA FLS, which was potentially regulated by miR-141-3p. The current findings suggested new perspectives into understanding the novel molecular signatures of FLS involved in the pathogenesis of OA, which may be potential therapeutic targets.

Genetic Heterogeneity and Evolution in Lymphoid Malignancies

Patterns of mutation and clonal evolution in relapsed acute lymphoblastic leukemiaRelapsed acute lymphoblastic leukemia remains a major cause of childhood cancer death, and this remains true despite the advent of new targeted and immunotherapeutic approaches. Recent years have witnessed the use of broad and deep serial genomic profiling approaches to dissect the relationship of genetic variegation to clonal evolution and relapse. Studies of over 90 children treated on St Jude Total Therapy protocols, incorporating genome, exome and transcriptome sequencing, coupled with limiting dilution xenografting to formally elucidate clonal structure have provided multiple key insights. In the majority of cases, the relapse-fated clone is a minor clone at diagnosis, that harbors resistance-enriched (and thus relapse-promoting) mutations at diagnosis, and/or acquires additional mutations the confer resistance after initial therapy. Approximately one third of cases relapse from a major clone, or show polyclonal evolution, and such cases typically have a shorter time to disease recurrence and relapse. A subset of cases exhibit complete discordance for somatic non-silent mutations, DNA copy number alterations and antigen receptor rearrangements between diagnosis and relapse, suggesting relapse represents a second leukemia; however such cases typically preserve the founding chromosomal rearrangement and a subset of non-coding mutations, indicating that relapse arises from an ancestral clone that has undergone divergent evolution early in leukemogenesis. Conversely, a subset of cases relapse with myeloid or lineage ambiguous leukemia but preserve genomic alterations indicating a common clonal origin but lineage plasticity: thus, careful genomic analysis is required to interpret the nature of disease recurrence/relapse. Approximately 15% of cases exhibit hypermutation, particularly in aneuploid leukemia and second or later relapse, associated with distinct mutational signatures and kinetics of hypermutation, thus identifying this process as a driver of treatment failure in a subset of ALL cases. Integrated analysis has identified over 80 recurrent targets of alteration at relapse that show variable patterns of enrichment in rising and falling clones. Importantly, several targets (e.g. NT5C2) are never identified at diagnosis despite deep sequencing approaches, suggesting adverse effects on leukemic fitness, and/or an absolute requirement of prior drug exposure to initiate mutagenesis. Integration of limiting dilution xenografting, coupled with genomic analysis of xenografts and drug exposure has not only formally confirmed and extended inferential clonal structures, but shown that in a subset of cases resistance is present at initial diagnosis, rather than being acquired after drug exposure. Finally, several groups have shown that the relationship of relapse-enriched mutations and relapse by be drug agnostic (e.g. IKZF1) or drug specific (e.g. NT5C2 and thiopurine resistance, and CREBBP and glucocorticoid resistance). As such mutations may now be detected at levels suitable for tracking of minimal residual disease, these insights offer the opportunity to identify the relapse-fated clone early in disease evolution, and modulate therapy accordingly to circumvent relapse. DisclosuresMullighan:Pfizer: Honoraria, Research Funding, Speakers Bureau; Cancer Prevention and Research Institute of Texas: Consultancy; Loxo Oncology: Research Funding; Abbvie: Research Funding; Amgen: Honoraria, Speakers Bureau.

Placebo Controlled Functional Analysis of Eltrombopag in a Preclinical Xenograft Model of Myelodysplastic Syndromes (MDS)

Introduction: Thrombocytopenia is a common complication among MDS patients. Thus, many patients are dependent on platelet (PLT) transfusions, which give short-term therapeutic relief but are also associated with considerable clinical risks. In this context, thrombopoietin receptor agonists (TRAs) are under investigation as alternative treatment option, albeit with the concern that these substances may promote adverse events in MDS. However, beside potential positive effects on thrombopoiesis in MDS patients the TRA Eltrombopag (EPAG) has also been shown to exert positive disease modifying effects in vitro (Roth et al., Blood 2012). Using a MDS xenograft model, we here investigate the efficacy of EPAG and its influence on clonal composition on primary patient derived MDS xenografts and present data from an ongoing study.Methods: Currently, samples from n=18 MDS patients (MDS del(5q)=2, MDS-MLD=6, MDS-RS-MLD=1 MDS-EB-1=2, MDS-EB-2=7) have been xenografted into NSG mice by intrafemoral co-injection of CD34+ hematopoietic stem cells and mesenchymal stromal cells using a modified protocol according to Medyouf et al., Cell Stem Cell 2014. Long term engraftment is assessed 12 weeks post-transplant by intrafemoral bone marrow (BM) biopsy and mice with positive human engraftment are subsequently treated with either EPAG (50mg/kg) or vehicle control for 18 weeks. During that time, the mice are bled every two weeks and BM aspiration is performed every six weeks. Human hematopoietic cells are FACS sorted. In peripheral blood, human PLTs are specifically and absolutely counted with a FACS assay based on hCD41+ cells and beads. To track clonal composition of MDS samples upon xenografting and EPAG treatment in comparison to placebo control, the original patient sample and the final MDS xenograft sample are being whole exome sequenced (WES). Interspersed time points are analyzed with a patient individual amplicon based deep sequencing approach (Mossner et al., Blood 2016) to calculate dynamics of variant allele frequencies (VAF) in dependency of treatment.Results: To date, n=12 patient samples have been analyzed for human engraftment after 12 weeks post-transplant. Of these, n=7 (58%) have shown positive human engraftment and are being treated with EPAG versus placebo. To this end, one case has been completely followed up, including final molecular analysis. This MDS high risk case (MDS-EB-2) with a clinical PLT count of 29x109 PLT/L was transplanted into n=3 NSG mice. While two mice treated with EPAG survived the complete duration of the experiment, the placebo mouse died prematurely due to severe weight loss after 6 weeks of treatment. Further, EPAG treatment led to an initial rise of human PLT levels, while the placebo treated mouse presented a continuous decline of human PLTs, showing the efficacy of EPAG on human xenografts in the model. This observation has been confirmed in another case currently still under treatment. Molecular tracking by WES confirmed MDS patient specific molecular lesions in the MDS xenograft such as monosomy 7 and the disease related mutations CBL, DNMT3A and EZH2 with VAFs of 83%/43%/23% respectively. The monosomy 7 was detectable in all mice. CBL and DNMT3A exhibited similar VAFs in mouse EPAG1 (VAF=100%/54%), EPAG2 (VAF=100%/34%) and placebo (VAF=100%/50%). The EZH2 mutation was only detected in mouse EPAG2 (VAF=11%). Interestingly, the placebo mouse acquired a de novo mutation of U2AF1 (VAF=10%), which was not detectable in the initial patient sample or the EPAG treated mice. This spliceosomal mutation is associated with a higher risk of transformation to AML and shorter survival (Graubert et al., Nat Genet 2012; Makishima et al., Blood 2012).Conclusions: Our data show first proof of principle results that new treatment options can be tested successfully in a preclinical murine xenograft model of primary MDS patient samples in a placebo controlled experimental setting. This approach allows the performance of patient individual substance testing that can segregate substance specific effects from natural disease progression in the same patient. Clinical parameters such as human PLT production and molecular clonal composition can be measured with a high confidence in vivo. Our current data show preliminary support for the hypothesis that EPAG may be efficacious in increasing PLT production in MDS patients without adversely influencing the underlying clonal composition. DisclosuresNowak:Novartis: Research Funding.

Genome-wide transcriptional response of microRNAs to the benzo(a)pyrene stress in amphioxus Branchiostoma belcheri

Amphioxus, a cephalochordate found in sand habitats in shallow in-shore seawaters, has been widely used as a model in comparative immunology of chordates. However, the role of microRNAs (miRNAs) in amphioxus under abiotic stress, particularly xenobiotics with strong toxicity, remains largely unknown. Here, a widespread marine contaminant, benzo(a)pyrene (BaP) is used to evaluate its toxic effects on miRNA expression of amphioxus. Six small RNA libraries were sequenced from Branchiostoma belcheri. A total of 144 known and 157 novel miRNAs were identified using deep sequencing and bioinformatics approaches. A total of 58 differentially expressed miRNAs (DEMs) were screened, including 25 up- and 33 down-regulated DEMs under BaP stress. Target genes possibly regulated by DEMs were predicted, and their functional enrichment analyses were performed. Targets of DEMs are primarily involved in xenobiotic and cellular homeostasis, catabolic and transport process. They could be largely linked to nine immune- and toxin detoxification-related pathways, including metabolism of xenobiotics by cytochrome P450, drug metabolism-other enzymes, and drug metabolism-cytochrome P450, etc. Furthermore, quantitative real-time PCR (qRT-PCR) analysis for 12 key BaP-responsive DEMs validates the accuracy of deep sequencing. Experiments were then conducted to investigate their expression responses to BaP stress at different time intervals in detail to further determine their expression dynamic in responses of B. belcheri towards BaP exposure. This study, to the best of our knowledge, investigates the regulatory roles of miRNAs in the toxicological response of amphioxus for the first time, providing valuable information for the protection of lone existing cephalochordate amphioxus.

Deep sequencing prompts the modification of a real-time RT-PCR for the serotype-specific detection of polioviruses

Polioviruses are members of the Enterovirus C species and asymptomatic fecal shedding allows for their transmission and persistence in a community, as well as the emergence of vaccine-derived polioviruses. Using three serotype-specific real-time RT-PCR (rRT-PCR) assays, the shedding and circulation of oral poliovirus vaccine (OPV) strains was previously investigated in a prospective cohort of Mexican children, their contacts, and nearby sewage. Subsequently, a deep sequencing approach targeting the P1 genomic region was applied to characterize OPV strains previously detected by rRT-PCR. Amplifiable RNA was obtained for sequencing from 40.3% (58/144) of stool samples and 71.4% (15/21) of sewage using nucleic acids extracted directly from primary rRT-PCR-positive specimens. Sequencing detected one or more OPV serotypes in 62.1% (36/58) of stool and 53.3% (8/15) of sewage samples. All stool and sewage samples in which poliovirus was not detected by deep sequencing contained at least one non-polio enterovirus C (NPEV-C) strain. To improve screening specificity, a modified, two-step, OPV serotype-specific multiplex rRT-PCR was evaluated. In stool specimens, the overall agreement between the original assays and the multiplex was 70.3%. By serotype, the overall agreement was 95.7% for OPV serotype-1 (S1), 65.6% for S2, and 96.1% for S3. Furthermore, most original rRT-PCR positive/multiplex rRT-PCR negative results were collected in the summer and fall months, consistent with NPEV-C circulation patterns. In conclusion, this deep sequencing approach allowed for the characterization of OPV sequences directly from clinical samples and facilitated the implementation of a more specific multiplex rRT-PCR for OPV detection and serotyping.

The spatial and temporal dynamics of nuclear RNAi-targeted retrotransposon transcripts in Caenorhabditis elegans.

Nuclear RNA interference provides a unique approach to the study of RNA-mediated transgenerational epigenetic inheritance. A paradox in the field is that expression of target loci is necessary for the initiation and maintenance of their silencing. How expression and repression are coordinated during animal development is poorly understood. To resolve this gap, we took imaging, deep-sequencing and genetic approaches towards delineating the developmental regulation and subcellular localization of RNA transcripts of two representative endogenous targets, the LTR retrotransposons Cer3 and Cer8. By examining wild-type worms and a collection of mutant strains, we found that the expression and silencing cycle of Cer3 and Cer8 is coupled with embryonic and germline development. Strikingly, endogenous targets exhibit a hallmark of nuclear enrichment of their RNA transcripts. In addition, germline and somatic repressions of Cer3 have different genetic requirements for three heterochromatin enzymes, MET-2, SET-25 and SET-32, in conjunction with the nuclear Argonaute protein HRDE-1. These results provide the first comprehensive cellular and developmental characterization of nuclear RNAi activities throughout the animal reproductive cycle.

Identification and characterization of microRNAs in a cestode Hydatigera taeniaeformis using deep sequencing approach

Hydatigera taeniaeformis (formerly known as Taenia taeniaeformis) is a parasitic tapeworm that has a worldwide distribution. H. taeniaeformis is naturally transmitted between mice and cats and threatens to human health, especially those who are in close contact with pets. MicroRNAs (miRNAs) are a class of small regulatory non-coding RNAs involved in the regulation of parasite growth and development, parasite infection and immunology, and host-pathogen interactions. The miRNA profile of H. taeniaeformis remains to be elucidated. Herein, 47 conserved miRNAs (grouped into 34 miRNA families) and 4 novel miRNAs were identified in H. taeniaeformis metacestodes using deep sequencing approach. Among them, hta-miR-71, -let-7, and-miR-87 was absolutely predominant in H. taeniaeformis metacestodes. Moreover, comparative analysis revealed the presence of miR-71/2 and miR-4989/277 clusters in H. taeniaeformis. Nucleotide bias analysis of identified miRNAs showed that the adenine (A) was the dominant nucleotide at the beginning of the miRNAs, particularly at the positions of third and 7th nucleotides. The study provides rich data for further understandings of H. taeniaeformis biology.

Identification and expression analysis of miRNA in hybrid snakehead by deep sequencing approach and their targets prediction

MicroRNAs (miRNAs) play important regulatory roles in numerous biological processes, but there is no report on miRNAs of hybrid snakehead. In this study, four independent small RNA libraries were constructed from the spleen, liver kidney and muscle of hybrid snakehead. These libraries were sequenced using deep sequencing technology, as result, a total of 1,067,172, 1,152,002, 1,653,941 and 970,866 clean reads from these four libraries were obtained. 252 known miRNAs and 63 putative novel miRNAs in these small RNA dataset were identified. The stem-loop RT-qPCR analysis indicated that eight known miRNAs and two novel miRNAs show different expression in eight different kinds of tissues. For better understanding the functions of miRNAs, 95,947 predicated target genes were analyzed by GO and their pathways, the results indicated that these targets of the identified miRNAs are involved in a broad range of physiological functions.

Species distinctions among closely related strains of Eustigmatophyceae (Stramenopiles) emphasizing ITS2 sequence-structure data: Eustigmatos and Vischeria

ABSTRACTThere is an increasing interest in the Eustigmatophyceae, a class of stramenopile microalgae, because they offer a variety of high-value health-beneficial compounds, e.g. polyunsaturated fatty acids (PUFAs), while concomitantly producing high biomass. Clarification of the taxonomy of these organisms at the species level is important in order to achieve reproducible results and constant yields of valuable compounds in their exploitation. Here the distinction of the, so far exclusively, morphologically defined species of the genera Eustigmatos and Vischeria was tested. Distinctions inferred from almost full 18S and ITS2 rRNA as well as plastid-encoded rbcL gene sequences were evaluated following a morphological investigation. The ITS2 secondary-structure-based phylogenies separated independent lineages (species) with long internal branches. This recommends ITS2 as a promising marker for a DNA metabarcoding approach (culture-independent biodiversity assessment). In contrast, the 18S V4 region which is commonly used in metabarcoding was almost invariant, whereas the almost full length sequences distinguished eight groups/types of strains. Monophyly of the species was supported by shared ITS2 secondary structure features, making them distinct from other eustigmatophyte lineages in concordance with phylogenetic analyses. No groups of strains were congruently supported by all three markers. Consequently, the previous distinction of two genera on the basis of morphology cannot be retained and the species should be accommodated in a single genus, Vischeria. Taxonomic changes among the species with the definition of epitypes, on the basis of cryopreserved strains, are recommended. Two findings point to a more complex evolutionary history of the species. The rbcL and nuclear markers resulted in disparate groupings of strains. In three species divergent intragenomic ITS2 paralogues were revealed. Therefore, a still broader taxon sampling, in conjunction with a deep sequencing approach, is needed for a more comprehensive understanding of the complex evolution of eustigmatophyte species.

HHV-6 encoded small non-coding RNAs define an intermediate and early stage in viral reactivation

Human herpesvirus 6A and 6B frequently acquires latency. HHV-6 activation has been associated with various human diseases. Germ line inheritance of chromosomally integrated HHV-6 makes viral DNA-based analysis difficult for determination of early stages of viral activation. We characterized early stages of HHV-6 activation using high throughput transcriptomics studies and applied the results to understand virus activation under clinical conditions. Using a latent HHV-6A cell culture model in U2OS cells, we identified an early stage of viral reactivation, which we define as transactivation that is marked by transcription of several viral small non-coding RNAs (sncRNAs) in the absence of detectable increase in viral replication and proteome. Using deep sequencing approaches, we detected previously known as well as a new viral sncRNAs that characterized viral transactivation and differentiated it from latency. Here we show changes in human transcriptome upon viral transactivation that reflect multiple alterations in mitochondria-associated pathways, which was supported by observation of increased mitochondrial fragmentation in virus reactivated cells. Furthermore, we present here a unique clinical case of DIHS/DRESS associated death where HHV-6 sncRNA-U14 was abundantly detected throughout the body of the patient in the presence of low viral DNA. In this study, we have identified a unique and early stage of viral activation that is characterized by abundant transcription of viral sncRNAs, which can serve as an ideal biomarker under clinical conditions.