De novo transcriptome assembly with ABySS

Whole genome and transcriptome integrated analyses guide clinical care of pediatric poor prognosis cancers

The role for routine whole genome and transcriptome analysis (WGTA) for poor prognosis pediatric cancers remains undetermined. Here, we characterize somatic mutations, structural rearrangements, copy number variants, gene expression, immuno-profiles and germline cancer predisposition variants in children and adolescents with relapsed, refractory or poor prognosis malignancies who underwent somatic WGTA and matched germline sequencing. Seventy-nine participants with a median age at enrollment of 8.8 y (range 6 months to 21.2 y) are included. Germline pathogenic/likely pathogenic variants are identified in 12% of participants, of which 60% were not known prior. Therapeutically actionable variants are identified by targeted gene report and whole genome in 32% and 62% of participants, respectively, and increase to 96% after integrating transcriptome analyses. Thirty-two molecularly informed therapies are pursued in 28 participants with 54% achieving a clinical benefit rate; objective response or stable disease ≥6 months. Integrated WGTA identifies therapeutically actionable variants in almost all tumors and are directly translatable to clinical care of children with poor prognosis cancers.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Roast: a tool for reference-free optimization of supertranscriptome assemblies

BACKGROUND

Transcriptomic studies involving organisms for which reference genomes are not available typically start by generating de novo transcriptome or supertranscriptome assembly from the raw RNA-seq reads. Assembling a supertranscriptome is, however, a challenging task due to significantly varying abundance of mRNA transcripts, alternative splicing, and sequencing errors. As a result, popular de novo supertranscriptome assembly tools generate assemblies containing contigs that are partially-assembled, fragmented, false chimeras or have local mis-assemblies leading to decreased assembly accuracy. Commonly available tools for assembly improvement rely primarily on running BLAST using closely related species making their accuracy and reliability conditioned on the availability of the data for closely related organisms.

RESULTS

We present ROAST, a tool for optimization of supertranscriptome assemblies that uses paired-end RNA-seq data from Illumina sequencing platform to iteratively identify and fix assembly errors solely using the error signatures generated by RNA-seq alignment tools including soft-clips, unexpected expression coverage, and reads with mates unmapped or mapped on a different contig to identify and fix various supertranscriptome assembly errors without performing BLAST searches against other organisms. Evaluation results using simulated as well as real datasets show that ROAST significantly improves assembly quality by identifying and fixing various assembly errors.

CONCLUSION

ROAST provides a reference-free approach to optimizing supertranscriptome assemblies highlighting its utility in refining de novo supertranscriptome assemblies of non-model organisms.

Evolution and expression of the duck TRIM gene repertoire

Tripartite motif (TRIM) proteins are involved in development, innate immunity, and viral restriction. TRIM gene repertoires vary between species, likely due to diversification caused by selective pressures from pathogens; however, this has not been explored in birds. We mined a de novo assembled transcriptome for the TRIM gene repertoire of the domestic mallard duck (Anas platyrhynchos), a reservoir host of influenza A viruses. We found 57 TRIM genes in the duck, which represent all 12 subfamilies based on their C-terminal domains. Members of the C-IV subfamily with C-terminal PRY-SPRY domains are known to augment immune responses in mammals. We compared C-IV TRIM proteins between reptiles, birds, and mammals and show that many C-IV subfamily members have arisen independently in these lineages. A comparison of the MHC-linked C-IV TRIM genes reveals expansions in birds and reptiles. The TRIM25 locus with related innate receptor modifiers is adjacent to the MHC in reptile and marsupial genomes, suggesting the ancestral organization. Within the avian lineage, both the MHC and TRIM25 loci have undergone significant TRIM gene reorganizations and divergence, both hallmarks of pathogen-driven selection. To assess the expression of TRIM genes, we aligned RNA-seq reads from duck tissues. C-IV TRIMs had high relative expression in immune relevant sites such as the lung, spleen, kidney, and intestine, and low expression in immune privileged sites such as in the brain or gonads. Gene loss and gain in the evolution of the TRIM repertoire in birds suggests candidate immune genes and potential targets of viral subversion.

Midgut transcriptomic responses to dengue and chikungunya viruses in the vectors Aedes albopictus and Aedes malayensis

Dengue (DENV) and chikungunya (CHIKV) viruses are among the most preponderant arboviruses. Although primarily transmitted through the bite of Aedes aegypti mosquitoes, Aedes albopictus and Aedes malayensis are competent vectors and have an impact on arbovirus epidemiology. Here, to fill the gap in our understanding of the molecular interactions between secondary vectors and arboviruses, we used transcriptomics to profile the whole-genome responses of A. albopictus to CHIKV and of A. malayensis to CHIKV and DENV at 1 and 4 days post-infection (dpi) in midguts. In A. albopictus, 1793 and 339 genes were significantly regulated by CHIKV at 1 and 4 dpi, respectively. In A. malayensis, 943 and 222 genes upon CHIKV infection, and 74 and 69 genes upon DENV infection were significantly regulated at 1 and 4 dpi, respectively. We reported 81 genes that were consistently differentially regulated in all the CHIKV-infected conditions, identifying a CHIKV-induced signature. We identified expressed immune genes in both mosquito species, using a de novo assembled midgut transcriptome for A. malayensis, and described the immune architectures. We found the JNK pathway activated in all conditions, generalizing its antiviral function to Aedines. Our comprehensive study provides insight into arbovirus transmission by multiple Aedes vectors.

Blood-feeding adaptations and virome assessment of the poultry red mite Dermanyssus gallinae guided by RNA-seq

Dermanyssus gallinae is a blood-feeding mite that parasitises wild birds and farmed poultry. Its remarkably swift processing of blood, together with the capacity to blood-feed during most developmental stages, makes this mite a highly debilitating pest. To identify specific adaptations to digestion of a haemoglobin-rich diet, we constructed and compared transcriptomes from starved and blood-fed stages of the parasite and identified midgut-enriched transcripts. We noted that midgut transcripts encoding cysteine proteases were upregulated with a blood meal. Mapping the full proteolytic apparatus, we noted a reduction in the suite of cysteine proteases, missing homologues for Cathepsin B and C. We have further identified and phylogenetically analysed three distinct transcripts encoding vitellogenins that facilitate the reproductive capacity of the mites. We also fully mapped transcripts for haem biosynthesis and the ferritin-based system of iron storage and inter-tissue trafficking. Additionally, we identified transcripts encoding proteins implicated in immune signalling (Toll and IMD pathways) and activity (defensins and thioester-containing proteins), RNAi, and ion channelling (with targets for commercial acaricides such as Fluralaner, Fipronil, and Ivermectin). Viral sequences were filtered from the Illumina reads and we described, in part, the RNA-virome of D. gallinae with identification of a novel virus, Red mite quaranjavirus 1.

Homologous recombination deficiency signatures in gastrointestinal and thoracic cancers correlate with platinum therapy duration

There is emerging evidence about the predictive role of homologous recombination deficiency (HRD), but this is less defined in gastrointestinal (GI) and thoracic malignancies. We reviewed whole genome (WGS) and transcriptomic (RNA-Seq) data from advanced GI and thoracic cancers in the Personalized OncoGenomics trial (NCT02155621) to evaluate HRD scores and single base substitution (SBS)3, which is associated with BRCA1/2 mutations and potentially predictive of defective HRD. HRD scores were calculated by sum of loss of heterozygosity, telomeric allelic imbalance, and large-scale state transitions scores. Regression analyses examined the association between HRD and time to progression on platinum (TTPp). We included 223 patients with GI (n = 154) or thoracic (n = 69) malignancies. TTPp was associated with SBS3 (p < 0.01) but not HRD score in patients with GI malignancies, whereas neither was associated with TTPp in thoracic malignancies. Tumors with gBRCA1/2 mutations and a somatic second alteration exhibited high SBS3 and HRD scores, but these signatures were also present in several tumors with germline but no somatic second alterations, suggesting silencing of the wild-type allele or BRCA1/2 haploinsufficiency. Biallelic inactivation of an HR gene, including loss of XRCC2 and BARD1, was identified in BRCA1/2 wild-type HRD tumors and these patients had prolonged response to platinum. Thoracic cases with high HRD score were associated with high RECQL5 expression (p ≤ 0.025), indicating another potential mechanism of HRD. SBS3 was more strongly associated with TTPp in patients with GI malignancies and may be complementary to using HRD and BRCA status in identifying patients who benefit from platinum therapy.

A novel splicing variant of DJ-1 in Parkinson's disease induces mitochondrial dysfunction

Several studies have identified mutations in neuroprotective genes in a few cases of Parkinson's disease (PD); however, the role of alternative splicing changes in PD remains unelucidated. Based on the transcriptome analysis of substantia nigra (SN) tissues obtained from PD cases and age-matched healthy controls, we identified a novel alternative splicing variant of DJ-1, lacking exon 6 (DJ-1 ^ΔE6), frequently detected in the SN of patients with PD. We found that the exon 6 skipping of DJ-1 induces mitochondrial dysfunction and impaired antioxidant capability. According to an in silico modeling study, the exon 6 skipping of DJ-1 disrupts the structural state suitable for the oxidation of the cysteine 106 residue that is a prerequisite for activating its neuroprotective roles. Our results suggest that change in DJ-1 alternative splicing may contribute to PD progression and provide an insight for studying PD etiology and its potential therapeutic targets.

Improving Coconut Using Modern Breeding Technologies: Challenges and Opportunities

Coconut (Cocos nucifera L.) is a perennial palm with a wide range of distribution across tropical islands and coastlines. Multitude use of coconut by nature is important in the socio-economic fabric framework among rural smallholders in producing countries. It is a major source of income for 30 million farmers, while 60 million households rely on the coconut industry directly as farm workers and indirectly through the distribution, marketing, and processing of coconut and coconut-based products. Stagnant production, inadequate planting materials, the effects of climate change, as well as pests and diseases are among the key issues that need to be urgently addressed in the global coconut industry. Biotechnology has revolutionized conventional breeding approaches in creating genetic variation for trait improvement in a shorter period of time. In this review, we highlighted the challenges of current breeding strategies and the potential of biotechnological approaches, such as genomic-assisted breeding, next-generation sequencing (NGS)-based genotyping and genome editing tools in improving the coconut. Also, combining these technologies with high-throughput phenotyping approaches and speed breeding could speed up the rate of genetic gain in coconut breeding to solve problems that have been plaguing the industry for decades.

Exceptional response to combination ipilimumab and nivolumab in metastatic uveal melanoma: Insights from genomic analysis

Uveal melanoma is the most common intraocular malignancy and has a poor prognosis compared to other melanoma subtypes with a median overall survival of 6-10 months. With immune checkpoint inhibitor therapy, either PD-1 inhibitor alone or combination ipilimumab/nivolumab (anti-CTLA-4/anti-PD-1), responses are rare and often not durable. We present a case report of a now 66-year-old woman with diffuse metastatic uveal melanoma previously treated with a combination of ipilimumab/nivolumab, followed by maintenance nivolumab. Almost complete resolution of all sites of metastatic disease was observed except for one liver metastasis which regressed partially on immunotherapy. Notably, the patient had a significantly elevated BMI and developed widespread vitiligo on treatment. Whole-genome and transcriptome analysis was performed on the residual liver biopsy and molecular markers that may have contributed to the exceptional response were investigated. Several alterations were observed in genes involved in T-cell responses. Estimates of tumour infiltrating immune cells indicated a high level of plasma cells compared to other uveal melanoma cases, a finding previously associated with indolent disease. The patient also carried several germline SNPs that may have contributed to her treatment response as well as widespread vitiligo. Whole-genome and transcriptome sequencing have provided insight into potential molecular underpinnings of an exceptional treatment response in a tumour type typically associated with poor prognosis. Immunological findings suggest a role for plasma cells in the tumour microenvironment. Elevated BMI and the development of vitiligo may be clinically relevant factors for predicting response to immune checkpoint inhibitor therapy, warranting further studies in patients with uveal melanoma.

The impact of whole genome and transcriptome analysis (WGTA) on predictive biomarker discovery and diagnostic accuracy of advanced malignancies

In this study, we evaluate the impact of whole genome and transcriptome analysis (WGTA) on predictive molecular profiling and histologic diagnosis in a cohort of advanced malignancies. WGTA was used to generate reports including molecular alterations and site/tissue of origin prediction. Two reviewers analyzed genomic reports, clinical history, and tumor pathology. We used National Comprehensive Cancer Network (NCCN) consensus guidelines, Food and Drug Administration (FDA) approvals, and provincially reimbursed treatments to define genomic biomarkers associated with approved targeted therapeutic options (TTOs). Tumor tissue/site of origin was reassessed for most cases using genomic analysis, including a machine learning algorithm (Supervised Cancer Origin Prediction Using Expression [SCOPE]) trained on The Cancer Genome Atlas data. WGTA was performed on 652 cases, including a range of primary tumor types/tumor sites and 15 malignant tumors of uncertain histogenesis (MTUH). At the time WGTA was performed, alterations associated with an approved TTO were identified in 39 (6%) cases; 3 of these were not identified through routine pathology workup. In seven (1%) cases, the pathology workup either failed, was not performed, or gave a different result from the WGTA. Approved TTOs identified by WGTA increased to 103 (16%) when applying 2021 guidelines. The histopathologic diagnosis was reviewed in 389 cases and agreed with the diagnostic consensus after WGTA in 94% of non‐MTUH cases (n = 374). The remainder included situations where the morphologic diagnosis was changed based on WGTA and clinical data (0.5%), or where the WGTA was non‐contributory (5%). The 15 MTUH were all diagnosed as specific tumor types by WGTA. Tumor board reviews including WGTA agreed with almost all initial predictive molecular profile and histopathologic diagnoses. WGTA was a powerful tool to assign site/tissue of origin in MTUH. Current efforts focus on improving therapeutic predictive power and decreasing cost to enhance use of WGTA data as a routine clinical test.

Genes and evolutionary fates of the amanitin biosynthesis pathway in poisonous mushrooms

Why do unrelated poisonous mushrooms (Amanita, Galerina, and Lepiota) make the same deadly toxin, α-amanitin? One of the most effective and fast strategies for organisms to acquire new abilities is through horizontal gene transfer (HGT). With the help of genome sequencing and the finding of two genes for the amanitin biosynthetic pathway, we demonstrate that the pathway distribution resulted from HGT probably through an unknown ancestral fungal donor. In Amanita mushrooms, the pathway evolved, through a series of gene manipulations, to produce very high levels of toxins, generating “the deadliest mushroom known to mankind.”

The deadly toxin α-amanitin is a bicyclic octapeptide biosynthesized on ribosomes. A phylogenetically disjunct group of mushrooms in Agaricales (Amanita, Lepiota, and Galerina) synthesizes α-amanitin. This distribution of the toxin biosynthetic pathway is possibly related to the horizontal transfer of metabolic gene clusters among taxonomically unrelated mushrooms with overlapping habitats. Here, our work confirms that two biosynthetic genes, P450-29 and FMO1, are oxygenases important for amanitin biosynthesis. Phylogenetic and genetic analyses of these genes strongly support their origin through horizontal transfer, as is the case for the previously characterized biosynthetic genes MSDIN and POPB. Our analysis of multiple genomes showed that the evolution of the α-amanitin biosynthetic pathways in the poisonous agarics in the Amanita, Lepiota, and Galerina clades entailed distinct evolutionary pathways including gene family expansion, biosynthetic genes, and genomic rearrangements. Unrelated poisonous fungi produce the same deadly amanitin toxins using variations of the same pathway. Furthermore, the evolution of the amanitin biosynthetic pathway(s) in Amanita species generates a much wider range of toxic cyclic peptides. The results reported here expand our understanding of the genetics, diversity, and evolution of the toxin biosynthetic pathway in fungi.

A platform for oncogenomic reporting and interpretation

Manual interpretation of variants remains rate limiting in precision oncology. The increasing scale and complexity of molecular data generated from comprehensive sequencing of cancer samples requires advanced interpretative platforms as precision oncology expands beyond individual patients to entire populations. To address this unmet need, we introduce a Platform for Oncogenomic Reporting and Interpretation (PORI), comprising an analytic framework that facilitates the interpretation and reporting of somatic variants in cancer. PORI integrates reporting and graph knowledge base tools combined with support for manual curation at the reporting stage. PORI represents an open-source platform alternative to commercial reporting solutions suitable for comprehensive genomic data sets in precision oncology. We demonstrate the utility of PORI by matching 9,961 pan-cancer genome atlas tumours to the graph knowledge base, calculating therapeutically informative alterations, and making available reports describing select individual samples.

The interpretation of somatic variants in cancer is challenging due to the scale and complexity of sequencing data. Here, the authors present PORI, an open-source framework for interpreting somatic variants in cancer using graph knowledge base tools, automated reporting, and manual curation.

MultiTrans: An Algorithm for Path Extraction Through Mixed Integer Linear Programming for Transcriptome Assembly

Recent advances in RNA-seq technology have made identification of expressed genes affordable, and thus boosting repaid development of transcriptomic studies. Transcriptome assembly, reconstructing all expressed transcripts from RNA-seq reads, is an essential step to understand genes, proteins, and cell functions. Transcriptome assembly remains a challenging problem due to complications in splicing variants, expression levels, uneven coverage and sequencing errors. Here, we formulate the transcriptome assembly problem as path extraction on splicing graphs (or assembly graphs), and propose a novel algorithm MultiTrans for path extraction using mixed integer linear programming. MultiTrans is able to take into consideration coverage constraints on vertices and edges, the number of paths and the paired-end information simultaneously. We benchmarked MultiTrans against two state-of-the-art transcriptome assemblers, TransLiG and rnaSPAdes. Experimental results show that MultiTrans generates more accurate transcripts compared to TransLiG (using the same splicing graphs) and rnaSPAdes (using the same assembly graphs). MultiTrans is freely available at https://github.com/jzbio/MultiTrans.

Clinical response to nivolumab in an INI1-deficient pediatric chordoma correlates with immunogenic recognition of brachyury

Poorly differentiated chordoma (PDC) is a recently recognized subtype of chordoma characterized by expression of the embryonic transcription factor, brachyury, and loss of INI1. PDC primarily affects children and is associated with a poor prognosis and limited treatment options. Here we describe the molecular and immune tumour microenvironment profiles of two paediatric PDCs produced using whole-genome, transcriptome and whole-genome bisulfite sequencing (WGBS) and multiplex immunohistochemistry. Our analyses revealed the presence of tumour-associated immune cells, including CD8+ T cells, and expression of the immune checkpoint protein, PD-L1, in both patient samples. Molecular profiling provided the rationale for immune checkpoint inhibitor (ICI) therapy, which resulted in a clinical and radiographic response. A dominant T cell receptor (TCR) clone specific for a brachyury peptide-MHC complex was identified from bulk RNA sequencing, suggesting that targeting of the brachyury tumour antigen by tumour-associated T cells may underlie this clinical response to ICI. Correlative analysis with rhabdoid tumours, another INI1-deficient paediatric malignancy, suggests that a subset of tumours may share common immune phenotypes, indicating the potential for a therapeutically targetable subgroup of challenging paediatric cancers.

Anti-oxidant mechanisms of Chlorella pyrenoidosa under acute GenX exposure

GenX, a substitute for perfluorooctanoic acid, has been widely detected in surface water. Due to its bioaccumulation, toxicity and persistence, GenX can cause adverse effects such as oxidative damage on aquatic organisms. To investigate the toxicity of GenX and the anti-oxidant mechanism of algae under acute exposure, the growth, photosynthetic activity and gene expression of Chlorella pyrenoidosa (C. pyrenoidosa) were tested. Results showed that the growth of C. pyrenoidosa was inhibited under acute GenX exposure. The toxicity of GenX increased with time and concentration but was lower than that of the traditional perfluoroalkyl and polyfluoroalkyl substances (PFASs). Furthermore, with the increase of GenX concentration, the production of reactive oxygen species increased, while the level of the anti-oxidant enzyme first increased and then decreased. Changes in photosynthetic parameters also indicated that the photosynthetic system of C. pyrenoidosa was negatively affected by GenX exposure. Transcription analysis revealed that the up-regulation of genes related to the glutathione-ascorbate cycle and photosynthesis is a positive strategy to cope with the oxidative stress caused by acute GenX exposure. Our findings provide new insights into the interactions between emerging PFASs and aquatic organisms at the molecular level.

Independent evolutionary transitions to pueriparity across multiple timescales in the viviparous genus Salamandra

The ability to bear live offspring, viviparity, has evolved multiple times across the tree of life and is a remarkable adaptation with profound life-history and ecological implications. Within amphibians the ancestral reproductive mode is oviparity followed by a larval life stage, but viviparity has evolved independently in all three amphibian orders. Two types of viviparous reproduction can be distinguished in amphibians; larviparity and pueriparity. Larviparous amphibians deliver larvae into nearby ponds and streams, while pueriparous amphibians deliver fully developed juveniles and thus do not require waterbodies for reproduction. Among amphibians, the salamander genus Salamandra is remarkable as it exhibits both inter- and intraspecific variation in the occurrence of larviparity and pueriparity. While the evolutionary relationships among Salamandra lineages have been the focus of several recent studies, our understanding of how often and when transitions between modes occurred is still incomplete. Furthermore, in species with intraspecific variation, the reproductive mode of a given population can only be confirmed by direct observation of births and thus the prevalence of pueriparous populations is also incompletely documented. We used sequence capture to obtain 1,326 loci from 94 individuals from across the geographic range of the genus, focusing on potential reproductive mode transition zones. We also report additional direct observations of pueriparous births for 20 new locations and multiple lineages. We identify at least five independent transitions from the ancestral mode of larviparity to pueriparity among and within species, occurring at different evolutionary timescales ranging from the Pliocene to the Holocene. Four of these transitions occurred within species. Based on a distinct set of markers and analyses, we also confirm previous findings of introgression between species and the need for taxonomic revisions in the genus. We discuss the implications of our findings with respect to the evolution of this complex trait, and the potential of using five independent convergent transitions for further studies on the ecological context in which pueriparity evolves and the genetic architecture of this specialized reproductive mode.

Identification of potential new mosquito-associated viruses of adult Aedes aegypti mosquitoes from Tocantins state, Brazil

Medically important arboviruses such as dengue virus (DENV), Zika virus (ZIKV), and chikungunya virus (CHIKV) are primarily transmitted by the globally distributed mosquito Aedes aegypti. Increasing evidence suggests that the transmission of some viruses can be influenced by mosquito-specific and mosquito-borne viruses. Advancements in high-throughput sequencing (HTS) and bioinformatics have expanded our knowledge on the richness of viruses harbored by mosquitoes. HTS was used to characterize the presence of virus sequences in wild-caught adult Ae. aegypti from Tocantins (TO) state, Brazil. Samples of mosquitoes were collected in four cities of Tocantins state and submitted to RNA isolation, followed by sequencing at an Illumina HiSeq platform. Our results showed initially by Krona the presence of 3% of the sequenced reads belonging to the viral database. After further analysis, the virus sequences were found to have homology to two viral families found in insects Phenuiviridae and Metaviridae. Three possible viral strains including putative new viruses were detected and named Phasi Charoen-like phasivirus isolate To-1 (PCLV To-1), Aedes aegypti To virus 1 (AAToV1), and Aedes aegypti To virus 2 (AAToV2). The results presented in this work contribute to the growing knowledge about the diversity of viruses in mosquitoes and might be useful for future studies on the interaction between insect-specific viruses and arboviruses.

The Localization of Phytohormones within the Gall-inducing Insect Eurosta solidaginis (Diptera: Tephritidae)

The phytohormone production hypothesis suggests that organisms, including insects, induce galls by producing and secreting plant growth hormones. Auxins and cytokinins are classes of phytohormones that induce cell growth and cell division, which could contribute to the plant tissue proliferation which constitutes the covering gall. Bacteria, symbiotic with insects, may also play a part in gall induction by insects through the synthesis of phytohormones or other effectors. Past studies have shown that concentrations of cytokinins and auxins in gall-inducing insects are higher than in their host plants. However, these analyses have involved whole-body extractions. Using immunolocalization of cytokinin and auxin, in the gall inducing stage of Eurosta solidaginis, we found both phytohormones to localize almost exclusively to the salivary glands. Co-localization of phytohormone label with a nucleic acid stain in the salivary glands revealed the absence of Wolbachia sp., the bacterial symbiont of E. solidaginis, which suggests that phytohormone production is symbiont independent. Our findings are consistent with the hypothesis that phytohormones are synthesized in and secreted from the salivary glands of E. solidaginis into host-plant tissues for the purpose of manipulating the host plant.

IsoDetect: Detection of Splice Isoforms from Third Generation Long Reads Based on Short Feature Sequences

Background:

Transcriptome annotation is the basis for understanding gene structures and analysing gene expression. The transcriptome annotation of many organisms such as humans is far from incomplete, due partly to the challenge in the identification of isoforms that are produced from the same gene through alternative splicing. Third generation sequencing (TGS) reads provide unprecedented opportunity for detecting isoforms due to their long length that exceeds the length of most isoforms. One limitation of current TGS reads-based isoform detection methods is that they are exclusively based on sequence reads, without incorporating the sequence information of annotated isoforms.

Objective:

We aim to develop a method to detect isoforms by incorporating annotated isoforms.

Methods:

Based on annotated isoforms, we propose a splice isoform detection method called IsoDetect. First, the sequence at exon-exon junctions is extracted from annotated isoforms as “short feature sequences”, which is used to distinguish splice isoforms. Second, we align these feature sequences to long reads and partition long reads into groups that contain the same set of feature sequences, thereby avoiding the pair-wise comparison among the large number of long reads. Third, clustering and consensus generation are carried out based on sequence similarity. For the long reads that do not contain any short feature sequence, clustering analysis based on sequence similarity is performed to identify isoforms. Therefore, our method can detect not only known but also novel isoforms.

Result:

Tested on two datasets from Calypte anna and Zebra Finch, IsoDetect shows higher speed and good accuracies compared with four existing methods.

Conclusion:

IsoDetect may become a promising method for isoform detection.

The Central Role of Salivary Metalloproteases in Host Acquired Resistance to Tick Feeding

During feeding on vertebrate hosts, ticks secrete saliva composed of a rich cocktail of bioactive molecules modulating host immune responses. Although most of the proteinaceous fraction of tick saliva is of little immunogenicity, repeated feeding of ticks on mammalian hosts may lead to impairment of tick feeding, preventing full engorgement. Here, we challenged rabbits with repeated feeding of both Ixodes ricinus nymphs and adults and observed the formation of specific antibodies against several tick salivary proteins. Repeated feeding of both I. ricinus stages led to a gradual decrease in engorged weights. To identify the salivary antigens, isolated immunoglobulins from repeatedly infested rabbits were utilized for a protein pull-down from the saliva of pilocarpine-treated ticks. Eluted antigens were first identified by peptide mass fingerprinting with the aid of available I. ricinus salivary gland transcriptomes originating from early phases of tick feeding. To increase the authenticity of immunogens identified, we also performed, for the first time, de novo assembly of the sialome from I. ricinus females fed for six days, a timepoint used for pilocarpine-salivation. The most dominant I. ricinus salivary immunogens identified in our study were zinc-dependent metalloproteases of three different families. To corroborate the role of metalloproteases at the tick/host interface, we fed ticks micro-injected with a zinc metalloprotease inhibitor, phosphoramidon, on a rabbit. These ticks clearly failed to initiate feeding and to engorge. However, neither feeding to ticks immune blood of repeatedly infested rabbits, nor phosphoramidon injection into ticks, prevented their engorgement when fed in vitro on an artificial membrane system. These data show that Zn metalloproteases play a decisive role in the success of tick feeding, mediated by complex molecular interactions between the host immune, inflammatory, and hemostatic processes, which are absent in in vitro feeding. This basic concept warrants further investigation and reconsideration of the current strategies towards the development of an effective “anti-tick” vaccine.

Uncovering Clinically Relevant Gene Fusions with Integrated Genomic and Transcriptomic Profiling of Metastatic Cancers

PURPOSE

Gene fusions are important oncogenic drivers and many are actionable. Whole-genome and transcriptome (WGS and RNA-seq, respectively) sequencing can discover novel clinically relevant fusions.

EXPERIMENTAL DESIGN

Using WGS and RNA-seq, we reviewed the prevalence of fusions in a cohort of 570 patients with cancer, and compared prevalence to that predicted with commercially available panels. Fusions were annotated using a consensus variant calling pipeline (MAVIS) and required that a contig of the breakpoint could be constructed and supported from ≥2 structural variant detection approaches.

RESULTS

In 570 patients with advanced cancer, MAVIS identified 81 recurrent fusions by WGS and 111 by RNA-seq, of which 18 fusions by WGS and 19 by RNA-seq were noted in at least 3 separate patients. The most common fusions were EML4-ALK in thoracic malignancies (9/69, 13%), and CMTM8-CMTM7 in colorectal cancer (4/73, 5.5%). Combined genomic and transcriptomic analysis identified novel fusion partners for clinically relevant genes, such as NTRK2 (novel partners: SHC3, DAPK1), and NTRK3 (novel partners: POLG, PIBF1).

CONCLUSIONS

Utilizing WGS/RNA-seq facilitates identification of novel fusions in clinically relevant genes, and detected a greater proportion than commercially available panels are expected to find. A significant benefit of WGS and RNA-seq is the innate ability to retrospectively identify variants that becomes clinically relevant over time, without the need for additional testing, which is not possible with panel-based approaches.

A decade of de novo transcriptome assembly: Are we there yet?

A decade ago, de novo transcriptome assembly evolved as a versatile and powerful approach to make evolutionary assumptions, analyse gene expression, and annotate novel transcripts, in particular, for non-model organisms lacking an appropriate reference genome. Various tools have been developed to generate a transcriptome assembly, and even more computational methods depend on the results of these tools for further downstream analyses. In this issue of Molecular Ecology Resources, Freedman et al. (Mol Ecol Resourc 2020) present a comprehensive analysis of errors in de novo transcriptome assemblies across public data sets and different assembly methods. They focus on two implicit assumptions that are often violated: First, the assembly presents an unbiased view of the transcriptome. Second, the expression estimates derived from the assembly are reasonable, albeit noisy, approximations of the relative frequency of expressed transcripts. They show that appropriate filtering can reduce this bias but can also lead to the loss of a reasonable number of highly expressed transcripts. Thus, to partly alleviate the noise in expression estimates, they propose a new normalization method called length-rescaled CPM. Remarkably, the authors found considerable distortions at the nucleotide level, which leads to an underestimation of diversity in transcriptome assemblies. The study by Freedman et al. (Mol Ecol Resourc 2020) clearly shows that we have not yet reached "high-quality" in the field of transcriptome assembly. Above all, it helps researchers be aware of these problems and filter and interpret their transcriptome assembly data appropriately and with caution.

RNA-sequencing of the Nyssomyia neivai sialome: a sand fly-vector from a Brazilian endemic area for tegumentary leishmaniasis and pemphigus foliaceus

Leishmaniasis encompasses a spectrum of diseases caused by a protozoan belonging to the genus Leishmania. The parasite is transmitted by the bite of sand flies, which inoculate the promastigote forms into the host’s skin while acquiring a blood meal. Nyssomyia neivai is one of the main vectors of tegumentary leishmaniasis (TL) in Brazil. Southeastern Brazil is an endemic region for TL but also overlaps with an endemic focus for pemphigus foliaceus (PF), also known as Fogo Selvagem. Salivary proteins of sand flies, specifically maxadilan and LJM11, have been related to pemphigus etiopathogenesis in the New World, being proposed as an environmental trigger for autoimmunity. We present a comprehensive description of the salivary transcriptome of the N. neivai, using deep sequencing achieved by the Illumina protocol. In addition, we highlight the abundances of several N. neivai salivary proteins and use phylogenetic analysis to compare with Old- and New-World sand fly salivary proteins. The collection of protein sequences associated with the salivary glands of N. neivai can be useful for monitoring vector control strategies as biomarkers of N. neivai, as well as driving vector-vaccine design for leishmaniasis. Additionally, this catalog will serve as reference to screen for possible antigenic peptide candidates triggering anti-Desmoglein-1 autoantibodies.

iPAC: a genome-guided assembler of isoforms via phasing and combing paths

MOTIVATION

Full-length transcript reconstruction is very important and quite challenging for the widely used RNA-seq data analysis. Currently, available RNA-seq assemblers generally suffered from serious limitations in practical applications, such as low assembly accuracy and incompatibility with latest alignment tools.

RESULTS

We introduce iPAC, a new genome-guided assembler for reconstruction of isoforms, which revolutionizes the usage of paired-end and sequencing depth information via phasing and combing paths over a newly designed phasing graph. Tested on both simulated and real datasets, it is to some extent superior to all the salient assemblers of the same kind. Especially, iPAC is significantly powerful in recovery of lowly expressed transcripts while others are not.

AVAILABILITY AND IMPLEMENTATION

iPAC is freely available at http://sourceforge.net/projects/transassembly/files.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The sialotranscriptome of the gopher-tortoise tick, Amblyomma tuberculatum

The gopher tortoise tick, Amblyomma tuberculatum, is known to parasitize keystone ectotherm reptile species. The biological success of ticks requires precise mechanisms to evade host hemostatic and immune responses. Acquisition of a full blood meal requires attachment, establishment of the blood pool, and engorgement of the tick. Tick saliva contains molecules which counter the host responses to allow uninterrupted feeding on the host. RNASeq of the salivary glands of Amblyomma tuberculatum ticks were sequenced resulting in 138,030 pyrosequencing reads which were assembled into 29,991 contigs. A total of 1875 coding sequences were deduced from the transcriptome assembly, including 602 putative secretory and 982 putative housekeeping proteins. The annotated data sets are available as a hyperlinked spreadsheet. The sialotranscriptome assembled for this tick species made available a valuable resource for mining novel pharmacological activities and comparative analysis.

Integrated analysis of sialotranscriptome and sialoproteome of the brown dog tick Rhipicephalus sanguineus (s.l.): Insights into gene expression during blood feeding

Tick salivary glands secrete a complex saliva into their hosts which modulates vertebrate hemostasis, immunity and tissue repair mechanisms. Transcriptomic studies revealed a large number of transcripts coding for structural and secreted protein products in a single tick species. These transcripts are organized in several large families according to their products. Not all transcripts are expressed at the same time, transcription profile switches at intervals, characterizing the phenomenon of "sialome switching". In this work, using transcriptomic and proteomic analysis we explored the sialome of Rhipicephalus sanguineus (s.l.) adult female ticks feeding on a rabbit. The correlations between transcriptional and translational results in the different groups were evaluated, confirming the "sialome switching" and validating the idea that the expression switch may serve as a mechanism of escape from the host immunity. Recombination breakpoints were identified in lipocalin and metalloprotease families, indicating this mechanism could be a possible source of diversity in the tick sialome. Another remarkable observation was the identification of host-derived proteins as a component of tick salivary gland content. These results and disclosed sequences contribute to our understanding of tick feeding biology, to the development of novel anti-tick methods, and to the discovery of novel pharmacologically active products. SIGNIFICANCE: Ticks are a burden by themselves to humans and animals, and vectors of viral, bacterial, protozoal and helminthic diseases. Their saliva has anti-clotting, anti-platelet, vasodilatory and immunomodulatory activities that allows successful feeding and pathogen transmission. Previous transcriptomic studies indicate ticks to have over one thousand transcripts coding for secreted salivary proteins. These transcripts code for proteins of diverse families, but not all are transcribed simultaneously, but rather transiently, in a succession. Here we explored the salivary transcriptome and proteome of the brown dog tick, Rhipicephalus sanguineus. A protein database of over 20 thousand sequences was "de novo" assembled from over 600 million nucleotide reads, from where over two thousand polypeptides were identified by mass spectrometry. The proteomic data was shown to vary in time with the transcription profiles, validating the idea that the expression switch may serve as a mechanism of escape from the host immunity. Analysis of the transcripts coding for lipocalin and metalloproteases indicate their genes to contain signals of breakpoint recombination suggesting a new mechanism responsible for the large diversity in tick salivary proteins. These results and the disclosed sequences contribute to our understanding of the success ticks enjoy as ectoparasites, to the development of novel anti-tick methods, and to the discovery of novel pharmacologically active products.

Transcriptome profiling of venom gland from wasp species: de novo assembly, functional annotation, and discovery of molecular markers

Background

Vespa velutina, one of the most aggressive and fearful wasps in China, can cause grievous allergies and toxic reactions, leading to organ failure and even death. However, there is little evidence on molecular data regarding wasps. Therefore, we aimed to provide an insight into the transcripts expressed in the venom gland of wasps.

Results

In our study, high-throughput RNA sequencing was performed using the venom glands of four wasp species. First, the mitochondrial cytochrome C oxidase submit I (COI) barcoding and the neighbor joining (NJ) tree were used to validate the unique identity and lineage of each individual species. After sequencing, a total of 127,630 contigs were generated and 98,716 coding domain sequences (CDS) were predicted from the four species. The Gene ontology (GO) enrichment analysis of unigenes revealed their functional role in important biological processes (BP), molecular functions (MF) and cellular components (CC). In addition, c-type, p1 type, p2 type and p3 type were the most commonly found simple sequence repeat (SSR) types in the four species of wasp transcriptome. There were differences in the distribution of SSRs and single nucleotide polymorphisms (SNPs) among the four wasp species.

Conclusions

The transcriptome data generated in this study will improve our understanding on bioactive proteins and venom-related genes in wasp venom gland and provide a basis for pests control and other applications. To our knowledge, this is the first study on the identification of large-scale genomic data and the discovery of microsatellite markers from V. tropica ducalis and V. analis fabricius.

Molecular mechanisms underlying milk production and viviparity in the cockroach, Diploptera punctata

Viviparous reproduction is characterized by maternal retention of developing offspring within the reproductive tract during gestation, culminating in live birth. In some cases, a mother will provide nutrition beyond that present in the yolk; this is known as matrotrophic viviparity. While this phenomenon is best associated with mammals, it is observed in insects such as the viviparous cockroach, Diploptera punctata. Female D. punctata carry developing embryos in the brood sac, a reproductive organ that acts as both a uterus and a placenta by protecting and providing a nutritive secretion to the intrauterine developing progeny. While the basic physiology of D. punctata pregnancy has been characterized, little is known about the molecular mechanisms underlying this phenomenon. This study combined RNA-seq analysis, RNA interference, and other assays to characterize molecular and physiological changes associated with D. punctata reproduction. A comparison of four stages of the female reproductive cycle and males revealed unique gene expression profiles corresponding to each stage and between sexes. Differentially regulated transcripts of interest include the previously identified family of milk proteins and transcripts associated with juvenile hormone metabolism. RNA interference and methoprene application experiments established the potential impacts of bothbreakdown and synthesis reduction of juvenile hormone in maintaining pregnancy in D. punctata. These studies provide the comprehensive molecular mechanisms associated with cockroach viviparity, which will be a critical resource for comparative purposes among viviparity in insect systems.

Error, noise and bias in de novo transcriptome assemblies

De novo transcriptome assembly is a powerful tool, and has been widely used over the last decade for making evolutionary inferences. However, it relies on two implicit assumptions: that the assembled transcriptome is an unbiased representation of the underlying expressed transcriptome, and that expression estimates from the assembly are good, if noisy approximations of the relative abundance of expressed transcripts. Using publicly available data for model organisms, we demonstrate that, across assembly algorithms and data sets, these assumptions are consistently violated. Bias exists at the nucleotide level, with genotyping error rates ranging from 30% to 83%. As a result, diversity is underestimated in transcriptome assemblies, with consistent underestimation of heterozygosity in all but the most inbred samples. Even at the gene level, expression estimates show wide deviations from map-to-reference estimates, and positive bias at lower expression levels. Standard filtering of transcriptome assemblies improves the robustness of gene expression estimates but leads to the loss of a meaningful number of protein-coding genes, including many that are highly expressed. We demonstrate a computational method, length-rescaled CPM, to partly alleviate noise and bias in expression estimates. Researchers should consider ways to minimize the impact of bias in transcriptome assemblies.

Tick-Borne Encephalitis Virus Infection Alters the Sialome of Ixodes ricinus Ticks During the Earliest Stages of Feeding

Ticks are hematophagous arthropods that transmit a number of pathogens while feeding. Among these is tick-borne encephalitis virus (TBEV), a flavivirus transmitted by Ixodes ricinus ticks in the temperate zone of Europe. The infection results in febrile illness progressing to encephalitis and meningitis with a possibility of fatality or long-term neurological sequelae. The composition of tick saliva plays an essential role in the initial virus transmission during tick feeding. Ticks secrete a diverse range of salivary proteins to modulate the host response, such as lipocalins to control the itch and inflammatory response, and both proteases and protease inhibitors to prevent blood coagulation. Here, the effect of viral infection of adult females of Ixodes ricinus was studied with the goal of determining how the virus alters the tick sialome to modulate host tissue response at the site of infection. Uninfected ticks or those infected with TBEV were fed on mice and removed and dissected one- and 3-h post-attachment. RNA from the salivary glands of these ticks, as well as from unfed ticks, was extracted and subjected to next-generation sequencing to determine the expression of key secreted proteins at each timepoint. Genes showing statistically significant up- or down-regulation between infected and control ticks were selected and compared to published literature to ascertain their function. From this, the effect of tick viral infection on the modulation of the tick-host interface was determined. Infected ticks were found to differentially express a number of uncategorized genes, proteases, Kunitz-type serine protease inhibitors, cytotoxins, and lipocalins at different timepoints. These virus-induced changes to the tick sialome may play a significant role in facilitating virus transmission during the early stages of tick feeding.

Compacta: a fast contig clustering tool for de novo assembled transcriptomes

Background

RNA-Seq is the preferred method to explore transcriptomes and to estimate differential gene expression. When an organism has a well-characterized and annotated genome, reads obtained from RNA-Seq experiments can be directly mapped to that genome to estimate the number of transcripts present and relative expression levels of these transcripts. However, for unknown genomes, de novo assembly of RNA-Seq reads must be performed to generate a set of contigs that represents the transcriptome. These contig sets contain multiple transcripts, including immature mRNAs, spliced transcripts and allele variants, as well as products of close paralogs or gene families that can be difficult to distinguish. Thus, tools are needed to select a set of less redundant contigs to represent the transcriptome for downstream analyses. Here we describe the development of Compacta to produce contig sets from de novo assemblies.

Results

Compacta is a fast and flexible computational tool that allows selection of a representative set of contigs from de novo assemblies. Using a graph-based algorithm, Compacta groups contigs into clusters based on the proportion of shared reads. The user can determine the minimum coverage of the contigs to be clustered, as well as a threshold for the proportion of shared reads in the clustered contigs, thus providing a dynamic range of transcriptome compression that can be adapted according to experimental aims. We compared the performance of Compacta against state of the art clustering algorithms on assemblies from Arabidopsis, mouse and mango, and found that Compacta yielded more rapid results and had competitive precision and recall ratios. We describe and demonstrate a pipeline to tailor Compacta parameters to specific experimental aims.

Conclusions

Compacta is a fast and flexible algorithm for the determination of optimum contig sets that represent the transcriptome for downstream analyses.

Transcriptional profiling and physiological roles of Aedes aegypti spermathecal-related genes

BACKGROUND

Successful mating of female mosquitoes typically occurs once, with the male sperm being stored in the female spermatheca for every subsequent oviposition event. The female spermatheca is responsible for the maintenance, nourishment, and protection of the male sperm against damage during storage. Aedes aegypti is a major vector of arboviruses, including Yellow Fever, Dengue, Chikungunya, and Zika. Vector control is difficult due to this mosquito high reproductive capacity.

RESULTS

Following comparative RNA-seq analyses of spermathecae obtained from virgin and inseminated females, eight transcripts were selected based on their putative roles in sperm maintenance and survival, including energy metabolism, chitin components, transcriptional regulation, hormonal signaling, enzymatic activity, antimicrobial activity, and ionic homeostasis. In situ RNA hybridization confirmed tissue-specific expression of the eight transcripts. Following RNA interference (RNAi), observed outcomes varied between targeted transcripts, affecting mosquito survival, egg morphology, fecundity, and sperm motility within the spermathecae.

CONCLUSIONS

This study identified spermatheca-specific transcripts associated with sperm storage in Ae. aegypti. Using RNAi we characterized the role of eight spermathecal transcripts on various aspects of female fecundity and offspring survival. RNAi-induced knockdown of transcript AeSigP-66,427, coding for a Na+/Ca2+ protein exchanger, specifically interfered with egg production and reduced sperm motility. Our results bring new insights into the molecular basis of sperm storage and identify potential targets for Ae. aegypti control.

Exposure to the Florida red tide dinoflagellate, Karenia brevis, and its associated brevetoxins induces ecophysiological and proteomic alterations in Porites astreoides

As reef-building corals are increasingly being exposed to persistent threats that operate on both regional and global scales, there is a pressing need to better understand the complex processes that diminish coral populations. This study investigated the impacts of the Florida red tide dinoflagellate Karenia brevis and associated brevetoxins on selected facets of coral biology using Porites astreoides as a model system. When provided with choice assays, P. astreoides larvae were shown to actively avoid seawater containing red tide (5×10⁵ cells L^-1–7.6×10⁶ cells L^-1) or purified brevetoxins (0.018 μg mL^-1 brevetoxin-2 and 0.0018 μg mL^-1 brevetoxin-3). However, forced exposure to similar treatments induced time-dependent physiological and behavioral changes that were captured by PAM fluorometry and settlement and survival assays, respectively. Adult fragments of P. astreoides exposed to red tide or associated brevetoxins displayed signs of proteomic alterations that were characterized by the use of an iTRAQ-based quantitative proteomic analysis. The novel use of this technique with P. astreoides demonstrated that protein regulation was highly contingent upon biological versus chemical treatment (i.e. live K. brevis vs. solely brevetoxin exposure) and that several broad pathways associated with cell stress were affected including redox homeostasis, protein folding, energy metabolism and reactive oxygen species production. The results herein provide new insight into the ecology, behavior and sublethal stress of reef-building corals in response to K. brevis exposure and underscore the importance of recognizing the potential of red tide to act as a regional stressor to these important foundation species.

Fluorouracil sensitivity in a head and neck squamous cell carcinoma with a somatic DPYD structural variant

Head and neck squamous cell carcinoma (HNSCC) is one of the most common cancers worldwide and represents a heterogeneous group of tumors, the majority of which are treated with a combination of surgery, radiation, and chemotherapy. Fluoropyrimidine (5-FU) and its oral prodrug, capecitabine, are commonly prescribed treatments for several solid tumor types including HNSCC. 5-FU-associated toxicity is observed in ∼30% of treated patients and is largely caused by germline polymorphisms in DPYD, which encodes dihydropyrimidine dehydrogenase, a key enzyme of 5-FU catabolism and deactivation. Although the association of germline DPYD alterations with toxicity is well-described, the potential contribution of somatic DPYD alterations to 5-FU sensitivity has not been explored. In a patient with metastatic HNSCC, in-depth genomic and transcriptomic integrative analysis on a biopsy from a metastatic neck lesion revealed alterations in genes that are associated with 5-FU uptake and metabolism. These included a novel somatic structural variant resulting in a partial deletion affecting DPYD, a variant of unknown significance affecting SLC29A1, and homozygous deletion of MTAP. There was no evidence of deleterious germline polymorphisms that have been associated with 5-FU toxicity, indicating a potential vulnerability of the tumor to 5-FU therapy. The discovery of the novel DPYD variant led to the initiation of 5-FU treatment that resulted in a rapid response lasting 17 wk, with subsequent relapse due to unknown resistance mechanisms. This suggests that somatic alterations present in this tumor may serve as markers for tumor sensitivity to 5-FU, aiding in the selection of personalized treatment strategies.

De novo assembly of the Indian blue peacock (Pavo cristatus) genome using Oxford Nanopore technology and Illumina sequencing

Background

The Indian peafowl (Pavo cristanus) is native to South Asia and is the national bird of India. Here we present a draft genome sequence of the male blue peacock using Illumina and Oxford Nanopore technology (ONT).

Results

ONT sequencing gave ∼2.3-fold sequencing coverage, whereas Illumina generated 150–base pair paired-end sequence data at 284.6-fold coverage from 5 libraries. Subsequently, we generated a 0.915-gigabase pair de novo assembly of the peacock genome with a scaffold N50 of 0.23 megabase pairs (Mb). We predict that the peacock genome contains 23,153 protein-coding genes and 75.3 Mb (7.33%) of repetitive sequences.

Conclusions

We report a high-quality assembly of the peacock genome using a hybrid approach of sequences generated by both Illumina and ONT. The long-read chemistry generated by ONT was useful for addressing challenges related to de novo assembly, particularly at regions containing repetitive sequences spanning longer than the read length, and which could not be resolved with only short-read–based assembly. Contig assembly of Illumina short reads gave an N50 of 1,639 bases, whereas with ONT, the N50 increased by >9-fold to 14,749 bases. The initial contig assembly based on Illumina sequencing reads alone gave 685,241 contigs. Further scaffolding on assembled contigs using both Illumina and ONT sequencing reads resulted in a final assembly of 15,025 super-scaffolds, with an N50 of ∼0.23 Mb. Ninety-five percent of proteins predicted by homology matched with those in a public repository, verifying the completeness of our assembly. Like other phylogenetic studies of avian conserved genes, we found P. cristatus to be most closely related to Gallus gallus, followed by Meleagris gallopavo and Anas platyrhynchos. Compared with the recently published peacock genome assembly, the current, superior, hybrid assembly has greater sequencing depth, fewer non-ATGC sequences, and fewer scaffolds.

Comparative Analysis of Strategies for De Novo Transcriptome Assembly in Prokaryotes: Streptomyces clavuligerus as a Case Study

The performance of software tools for de novo transcriptome assembly greatly depends on the selection of software parameters. Up to now, the development of de novo transcriptome assembly for prokaryotes has not been as remarkable as that for eukaryotes. In this contribution, Rockhopper2 was used to perform a comparative transcriptome analysis of Streptomyces clavuligerus exposed to diverse environmental conditions. The study focused on assessing the incidence of software parameters on software performance for the identification of differentially expressed genes as a final goal. For this, a statistical optimization was performed using the Transrate Assembly Score (TAS). TAS was also used for evaluating the software performance and for comparing it with related tools, e.g., Trinity. Transcriptome redundancy and completeness were also considered for this analysis. Rockhopper2 and Trinity reached a TAS value of 0.55092 and 0.58337, respectively. Trinity assembles transcriptomes with high redundancy, with 55.6% of transcripts having some duplicates. Additionally, we observed that the total number of differentially expressed genes (DEG) and their annotation greatly depends on the method used for removing redundancy and the tools used for transcript quantification. To our knowledge, this is the first work aimed at assessing de novo assembly software for prokaryotic organisms.

Comprehensive genomic profiling of glioblastoma tumors, BTICs, and xenografts reveals stability and adaptation to growth environments

Glioblastoma multiforme (GBM) is the most deadly brain tumor, and currently lacks effective treatment options. Brain tumor-initiating cells (BTICs) and orthotopic xenografts are widely used in investigating GBM biology and new therapies for this aggressive disease. However, the genomic characteristics and molecular resemblance of these models to GBM tumors remain undetermined. We used massively parallel sequencing technology to decode the genomes and transcriptomes of BTICs and xenografts and their matched tumors in order to delineate the potential impacts of the distinct growth environments. Using data generated from whole-genome sequencing of 201 samples and RNA sequencing of 118 samples, we show that BTICs and xenografts resemble their parental tumor at the genomic level but differ at the mRNA expression and epigenomic levels, likely due to the different growth environment for each sample type. These findings suggest that a comprehensive genomic understanding of in vitro and in vivo GBM model systems is crucial for interpreting data from drug screens, and can help control for biases introduced by cell-culture conditions and the microenvironment in mouse models. We also found that lack of MGMT expression in pretreated GBM is linked to hypermutation, which in turn contributes to increased genomic heterogeneity and requires new strategies for GBM treatment.

An insight into the sialome, mialome and virome of the horn fly, Haematobia irritans

BACKGROUND

The horn fly (Haematobia irritans) is an obligate blood feeder that causes considerable economic losses in livestock industries worldwide. The control of this cattle pest is mainly based on insecticides; unfortunately, in many regions, horn flies have developed resistance. Vaccines or biological control have been proposed as alternative control methods, but the available information about the biology or physiology of this parasite is rather scarce.

RESULTS

We present a comprehensive description of the salivary and midgut transcriptomes of the horn fly (Haematobia irritans), using deep sequencing achieved by the Illumina protocol, as well as exploring the virome of this fly. Comparison of the two transcriptomes allow for identification of uniquely salivary or uniquely midgut transcripts, as identified by statistically differential transcript expression at a level of 16 x or more. In addition, we provide genomic highlights and phylogenetic insights of Haematobia irritans Nora virus and present evidence of a novel densovirus, both associated to midgut libraries of H. irritans.

CONCLUSIONS

We provide a catalog of protein sequences associated with the salivary glands and midgut of the horn fly that will be useful for vaccine design. Additionally, we discover two midgut-associated viruses that infect these flies in nature. Future studies should address the prevalence, biological effects and life cycles of these viruses, which could eventually lead to translational work oriented to the control of this economically important cattle pest.

Sketching and Sublinear Data Structures in Genomics

Large-scale genomics demands computational methods that scale sublinearly with the growth of data. We review several data structures and sketching techniques that have been used in genomic analysis methods. Specifically, we focus on four key ideas that take different approaches to achieve sublinear space usage and processing time: compressed full-text indices, approximate membership query data structures, locality-sensitive hashing, and minimizers schemes. We describe these techniques at a high level and give several representative applications of each.

Nanopore sequencing data analysis: state of the art, applications and challenges

The nanopore sequencing process is based on the transit of a DNA molecule through a nanoscopic pore, and since the 90s is considered as one of the most promising approaches to detect polymeric molecules. In 2014, Oxford Nanopore Technologies (ONT) launched a beta-testing program that supplied the scientific community with the first prototype of a nanopore sequencer: the MinION. Thanks to this program, several research groups had the opportunity to evaluate the performance of this novel instrument and develop novel computational approaches for analyzing this new generation of data. Despite the short period of time from the release of the MinION, a large number of algorithms and tools have been developed for base calling, data handling, read mapping, de novo assembly and variant discovery. Here, we face the main computational challenges related to the analysis of nanopore data, and we carry out a comprehensive and up-to-date survey of the algorithmic solutions adopted by the bioinformatic community comparing performance and reporting limits and advantages of using this new generation of sequences for genomic analyses. Our analyses demonstrate that the use of nanopore data dramatically improves the de novo assembly of genomes and allows for the exploration of structural variants with an unprecedented accuracy and resolution. However, despite the impressive improvements reached by ONT in the past 2 years, the use of these data for small-variant calling is still challenging, and at present, it needs to be coupled with complementary short sequences for mitigating the intrinsic biases of nanopore sequencing technology.

GAAP: A Genome Assembly + Annotation Pipeline

Genomic analysis begins with de novo assembly of short-read fragments in order to reconstruct full-length base sequences without exploiting a reference genome sequence. Then, in the annotation step, gene locations are identified within the base sequences, and the structures and functions of these genes are determined. Recently, a wide range of powerful tools have been developed and published for whole-genome analysis, enabling even individual researchers in small laboratories to perform whole-genome analyses on their objects of interest. However, these analytical tools are generally complex and use diverse algorithms, parameter setting methods, and input formats; thus, it remains difficult for individual researchers to select, utilize, and combine these tools to obtain their final results. To resolve these issues, we have developed a genome analysis pipeline (GAAP) for semiautomated, iterative, and high-throughput analysis of whole-genome data. This pipeline is designed to perform read correction, de novo genome (transcriptome) assembly, gene prediction, and functional annotation using a range of proven tools and databases. We aim to assist non-IT researchers by describing each stage of analysis in detail and discussing current approaches. We also provide practical advice on how to access and use the bioinformatics tools and databases and how to implement the provided suggestions. Whole-genome analysis of Toxocara canis is used as case study to show intermediate results at each stage, demonstrating the practicality of the proposed method.

Identifying similar transcripts in a related organism from de Bruijn graphs of RNA-Seq data, with applications to the study of salt and waterlogging tolerance in Melilotus

Background

A popular strategy to study alternative splicing in non-model organisms starts from sequencing the entire transcriptome, then assembling the reads by using de novo transcriptome assembly algorithms to obtain predicted transcripts. A similarity search algorithm is then applied to a related organism to infer possible function of these predicted transcripts. While some of these predictions may be inaccurate and transcripts with low coverage are often missed, we observe that it is possible to obtain a more complete set of transcripts to facilitate possible functional assignments by starting the search from the intermediate de Bruijn graph that contains all branching possibilities.

Results

We develop an algorithm to extract similar transcripts in a related organism by starting the search from the de Bruijn graph that represents the transcriptome instead of from predicted transcripts. We show that our algorithm is able to recover more similar transcripts than existing algorithms, with large improvements in obtaining longer transcripts and a finer resolution of isoforms. We apply our algorithm to study salt and waterlogging tolerance in two Melilotus species by constructing new RNA-Seq libraries.

Conclusions

We have developed an algorithm to identify paths in the de Bruijn graph that correspond to similar transcripts in a related organism directly. Our strategy bypasses the transcript prediction step in RNA-Seq data and makes use of support from evolutionary information.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5702-5) contains supplementary material, which is available to authorized users.

Genomic characterization of a well-differentiated grade 3 pancreatic neuroendocrine tumor

Pancreatic neuroendocrine neoplasms (PanNENs) represent a minority of pancreatic neoplasms that exhibit variability in prognosis. Ongoing mutational analyses of PanNENs have found recurrent abnormalities in chromatin remodeling genes (e.g., DAXX and ATRX), and mTOR pathway genes (e.g., TSC2, PTEN PIK3CA, and MEN1), some of which have relevance to patients with related familial syndromes. Most recently, grade 3 PanNENs have been divided into two groups based on differentiation, creating a new group of well-differentiated grade 3 neuroendocrine tumors (PanNETs) that have had a limited whole-genome level characterization to date. In a patient with a metastatic well-differentiated grade 3 PanNET, our study utilized whole-genome sequencing of liver metastases for the comparative analysis and detection of single-nucleotide variants, insertions and deletions, structural variants, and copy-number variants, with their biologic relevance confirmed by RNA sequencing. We found that this tumor most notably exhibited a TSC1-disrupting fusion, showed a novel CHD7–BEND2 fusion, and lacked any somatic variants in ATRX, DAXX, and MEN1.

RNA-seq analysis of the salivary glands and midgut of the Argasid tick Ornithodoros rostratus

Ornithodoros rostratus is a South American argasid tick which importance relies on its itchy bite and potential as disease vector. They feed on a wide variety of hosts and secrete different molecules in their saliva and intestinal content that counteract host defences and help to accommodate and metabolize the relatively large quantity of blood upon feeding. The present work describes the transcriptome profile of salivary gland (SG) and midgut (MG) of O. rostratus using Illumina sequencing. A total of 8,031 contigs were assembled and assigned to different functional classes. Secreted proteins were the most abundant in the SG and accounted for ~67% of all expressed transcripts with contigs with identity to lipocalins and acid tail proteins being the most representative. On the other hand, immunity genes were upregulated in MG with a predominance of defensins and lysozymes. Only 10 transcripts in SG and 8 in MG represented ~30% of all RNA expressed in each tissue and one single contig (the acid tail protein ORN-9707) represented ~7% of all expressed contigs in SG. Results highlight the functional difference of each organ and identified the most expressed classes and contigs of O. rostratus SG and MG.

An insight into the sialotranscriptome and virome of Amazonian anophelines

Background

Saliva of mosquitoes contains anti-platelet, anti-clotting, vasodilatory, anti-complement and anti-inflammatory substances that help the blood feeding process. The salivary polypeptides are at a fast pace of evolution possibly due to their relative lack of structural constraint and possibly also by positive selection on their genes leading to evasion of host immune pressure.

Results

In this study, we used deep mRNA sequence to uncover for the first time the sialomes of four Amazonian anophelines species (Anopheles braziliensis, A. marajorara, A. nuneztovari and A. triannulatus) and extend the knowledge of the A. darlingi sialome. Two libraries were generated from A. darlingi mosquitoes, sampled from two localities separated ~ 1100 km apart. A total of 60,016 sequences were submitted to GenBank, which will help discovery of novel pharmacologically active polypeptides and the design of specific immunological markers of mosquito exposure. Additionally, in these analyses we identified and characterized novel phasmaviruses and anpheviruses associated to the sialomes of A. triannulatus, A. marajorara and A. darlingi species.

Conclusions

Besides their pharmacological properties, which may be exploited for the development of new drugs (e.g. anti-thrombotics), salivary proteins of blood feeding arthropods may be turned into tools to prevent and/or better control vector borne diseases; for example, through the development of vaccines or biomarkers to evaluate human exposure to vector bites. The sialotranscriptome study reported here provided novel data on four New World anopheline species and allowed to extend our knowledge on the salivary repertoire of A. darlingi. Additionally, we discovered novel viruses following analysis of the transcriptomes, a procedure that should become standard within future RNAseq studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5545-0) contains supplementary material, which is available to authorized users.

Genomic resources for the Neotropical tree genus Cedrela (Meliaceae) and its relatives

BACKGROUND

Tree species in the genus Cedrela P. Browne are threatened by timber overexploitation across the Neotropics. Genetic identification of processed timber can be used to supplement wood anatomy to assist in the taxonomic and source validation of protected species and populations of Cedrela. However, few genetic resources exist that enable both species and source identification of Cedrela timber products. We developed several 'omic resources including a leaf transcriptome, organelle genome (cpDNA), and diagnostic single nucleotide polymorphisms (SNPs) that may assist the classification of Cedrela specimens to species and geographic origin and enable future research on this widespread Neotropical tree genus.

RESULTS

We designed hybridization capture probes to enrich for thousands of genes from both freshly preserved leaf tissue and from herbarium specimens across eight Meliaceae species. We first assembled a draft de novo transcriptome for C. odorata, and then identified putatively low-copy genes. Hybridization probes for 10,001 transcript models successfully enriched 9795 (98%) of these targets, and analysis of target capture efficiency showed that probes worked effectively for five Cedrela species, with each species showing similar mean on-target sequence yield and depth. The probes showed greater enrichment efficiency for Cedrela species relative to the other three distantly related Meliaceae species. We provide a set of candidate SNPs for species identification of four of the Cedrela species included in this analysis, and present draft chloroplast genomes for multiple individuals of eight species from four genera in the Meliaceae.

CONCLUSIONS

Deforestation and illegal logging threaten forest biodiversity globally, and wood screening tools offer enforcement agencies new approaches to identify illegally harvested timber. The genomic resources described here provide the foundation required to develop genetic screening methods for Cedrela species identification and source validation. Due to their transferability across the genus and family as well as demonstrated applicability for both fresh leaves and herbarium specimens, the genomic resources described here provide additional tools for studies examining the ecology and evolutionary history of Cedrela and related species in the Meliaceae.

The relapsing fever spirochete Borrelia turicatae persists in the highly oxidative environment of its soft-bodied tick vector

The relapsing fever spirochete Borrelia turicatae possesses a complex life cycle in its soft-bodied tick vector, Ornithodoros turicata. Spirochetes enter the tick midgut during a blood meal, and, during the following weeks, spirochetes disseminate throughout O. turicata. A population persists in the salivary glands allowing for rapid transmission to the mammalian hosts during tick feeding. Little is known about the physiological environment within the salivary glands acini in which B. turicatae persists. In this study, we examined the salivary gland transcriptome of O. turicata ticks and detected the expression of 57 genes involved in oxidant metabolism or antioxidant defences. We confirmed the expression of five of the most highly expressed genes, including glutathione peroxidase (gpx), thioredoxin peroxidase (tpx), manganese superoxide dismutase (sod-1), copper-zinc superoxide dismutase (sod-2), and catalase (cat) by reverse-transcriptase droplet digital polymerase chain reaction (RT-ddPCR). We also found distinct differences in the expression of these genes when comparing the salivary glands and midguts of unfed O. turicata ticks. Our results indicate that the salivary glands of unfed O. turicata nymphs are highly oxidative environments where reactive oxygen species (ROS) predominate, whereas midgut tissues comprise a primarily nitrosative environment where nitric oxide synthase is highly expressed. Additionally, B. turicatae was found to be hyperresistant to ROS compared with the Lyme disease spirochete Borrelia burgdorferi, suggesting it is uniquely adapted to the highly oxidative environment of O. turicata salivary gland acini.