Transcriptome sequencing and annotation of the polychaete Hermodice carunculata (Annelida, Amphinomidae).

BACKGROUND: The amphinomid polychaete Hermodice carunculata is a cosmopolitan and ecologically important omnivore in coral reef ecosystems, preying on a diverse suite of reef organisms and potentially acting as a vector for coral disease. While amphinomids are a key group for determining the root of the Annelida, their phylogenetic position has been difficult to resolve, and their publically available genomic data was scarce. RESULTS: We performed deep transcriptome sequencing (Illumina HiSeq) and profiling on Hermodice carunculata collected in the Western Atlantic Ocean. We focused this study on 58,454 predicted Open Reading Frames (ORFs) of genes longer than 200 amino acids for our homology search, and Gene Ontology (GO) terms and InterPro IDs were assigned to 32,500 of these ORFs. We used this de novo assembled transcriptome to recover major signaling pathways and housekeeping genes. We also identify a suite of H. carunculata genes related to reproduction and immune response. CONCLUSIONS: We provide a comprehensive catalogue of annotated genes for Hermodice carunculata and expand the knowledge of reproduction and immune response genes in annelids, in general. Overall, this study vastly expands the available genomic data for H. carunculata, of which previously consisted of only 279 nucleotide sequences in NCBI. This underscores the utility of Illumina sequencing for de novo transcriptome assembly in non-model organisms as a cost-effective and efficient tool for gene discovery and downstream applications, such as phylogenetic analysis and gene expression profiling.

Transcriptome deep-sequencing and clustering of expressed isoforms from Favia corals.

BACKGROUND: Genomic and transcriptomic sequence data are essential tools for tackling ecological problems. Using an approach that combines next-generation sequencing, de novo transcriptome assembly, gene annotation and synthetic gene construction, we identify and cluster the protein families from Favia corals from the northern Red Sea. RESULTS: We obtained 80 million 75 bp paired-end cDNA reads from two Favia adult samples collected at 65 m (Fav1, Fav2) on the Illumina GA platform, and generated two de novo assemblies using ABySS and CAP3. After removing redundancy and filtering out low quality reads, our transcriptome datasets contained 58,268 (Fav1) and 62,469 (Fav2) contigs longer than 100 bp, with N50 values of 1,665 bp and 1,439 bp, respectively. Using the proteome of the sea anemone Nematostella vectensis as a reference, we were able to annotate almost 20% of each dataset using reciprocal homology searches. Homologous clustering of these annotated transcripts allowed us to divide them into 7,186 (Fav1) and 6,862 (Fav2) homologous transcript clusters (E-value ≤ 2e(-30)). Functional annotation categories were assigned to homologous clusters using the functional annotation of Nematostella vectensis. General annotation of the assembled transcripts was improved 1-3% using the Acropora digitifera proteome. In addition, we screened these transcript isoform clusters for fluorescent proteins (FPs) homologs and identified seven potential FP homologs in Fav1, and four in Fav2. These transcripts were validated as bona fide FP transcripts via robust fluorescence heterologous expression. Annotation of the assembled contigs revealed that 1.34% and 1.61% (in Fav1 and Fav2, respectively) of the total assembled contigs likely originated from the corals' algal symbiont, Symbiodinium spp. CONCLUSIONS: Here we present a study to identify the homologous transcript isoform clusters from the transcriptome of Favia corals using a far-related reference proteome. Furthermore, the symbiont-derived transcripts were isolated from the datasets and their contribution quantified. This is the first annotated transcriptome of the genus Favia, a major increase in genomics resources available in this important family of corals.

Clinical Outcomes and Evolution of Clonal Hematopoiesis in Patients with Newly Diagnosed Multiple Myeloma.

Clonal hematopoiesis (CH) at time of autologous stem cell transplant (ASCT) has been shown to be associated with decreased overall survival (OS) and progression-free survival (PFS) in patients with multiple myeloma not receiving immunomodulatory drugs (IMiD). However, the significance of CH in newly diagnosed patients, including transplant ineligible patients, and its effect on clonal evolution during multiple myeloma therapy in the era of novel agents, has not been well studied. Using our new algorithm to differentiate tumor and germline mutations from CH, we detected CH in approximately 10% of 986 patients with multiple myeloma from the Clinical Outcomes in MM to Personal Assessment of Genetic Profile (CoMMpass) cohort (40/529 transplanted and 59/457 non-transplanted patients). CH was associated with increased age, risk of recurrent bacterial infections and cardiovascular disease. CH at time of multiple myeloma diagnosis was not associated with inferior OS or PFS regardless of undergoing ASCT, and all patients benefited from IMiD-based therapies, irrespective of the presence of CH. Serial sampling of 52 patients revealed the emergence of CH over a median of 3 years of treatment, increasing its prevalence to 25%, mostly with DNMT3A mutations. SIGNIFICANCE: Using our algorithm to differentiate tumor and germline mutations from CH mutations, we detected CH in approximately 10% of patients with newly diagnosed myeloma, including both transplant eligible and ineligible patients. Receiving IMiDs improved outcomes irrespective of CH status, but the prevalence of CH significantly rose throughout myeloma-directed therapy.

Cross center single-cell RNA sequencing study of the immune microenvironment in rapid progressing multiple myeloma.

Despite advancements in understanding the pathophysiology of Multiple Myeloma (MM), the cause of rapid progressing disease in a subset of patients is still unclear. MM's progression is facilitated by complex interactions with the surrounding bone marrow (BM) cells, forming a microenvironment that supports tumor growth and drug resistance. Understanding the immune microenvironment is key to identifying factors that promote rapid progression of MM. To accomplish this, we performed a multi-center single-cell RNA sequencing (scRNA-seq) study on 102,207 cells from 48 CD138- BM samples collected at the time of disease diagnosis from 18 patients with either rapid progressing (progression-free survival (PFS) < 18 months) or non-progressing (PFS > 4 years) disease. Comparative analysis of data from three centers demonstrated similar transcriptome profiles and cell type distributions, indicating subtle technical variation in scRNA-seq, opening avenues for an expanded multicenter trial. Rapid progressors depicted significantly higher enrichment of GZMK+ and TIGIT+ exhausted CD8+ T-cells (P = 0.022) along with decreased expression of cytolytic markers (PRF1, GZMB, GNLY). We also observed a significantly higher enrichment of M2 tolerogenic macrophages in rapid progressors and activation of pro-proliferative signaling pathways, such as BAFF, CCL, and IL16. On the other hand, non-progressive patients depicted higher enrichment for immature B Cells (i.e., Pre/Pro B cells), with elevated expression for markers of B cell development (IGLL1, SOX4, DNTT). This multi-center study identifies the enrichment of various pro-tumorigenic cell populations and pathways in those with rapid progressing disease and further validates the robustness of scRNA-seq data generated at different study centers.

Comprehensive Characterization of the Multiple Myeloma Immune Microenvironment Using Integrated scRNA-seq, CyTOF, and CITE-seq Analysis.

As part of the Multiple Myeloma Research Foundation (MMRF) immune atlas pilot project, we compared immune cells of multiple myeloma bone marrow samples from 18 patients assessed by single-cell RNA sequencing (scRNA-seq), mass cytometry (CyTOF), and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to understand the concordance of measurements among single-cell techniques. Cell type abundances are relatively consistent across the three approaches, while variations are observed in T cells, macrophages, and monocytes. Concordance and correlation analysis of cell type marker gene expression across different modalities highlighted the importance of choosing cell type marker genes best suited to particular modalities. By integrating data from these three assays, we found International Staging System stage 3 patients exhibited decreased CD4+ T/CD8+ T cells ratio. Moreover, we observed upregulation of RAC2 and PSMB9, in natural killer cells of fast progressors compared with those of nonprogressors, as revealed by both scRNA-seq and CITE-seq RNA measurement. This detailed examination of the immune microenvironment in multiple myeloma using multiple single-cell technologies revealed markers associated with multiple myeloma rapid progression which will be further characterized by the full-scale immune atlas project. Significance: scRNA-seq, CyTOF, and CITE-seq are increasingly used for evaluating cellular heterogeneity. Understanding their concordances is of great interest. To date, this study is the most comprehensive examination of the measurement of the immune microenvironment in multiple myeloma using the three techniques. Moreover, we identified markers predicted to be significantly associated with multiple myeloma rapid progression.

Clonal hematopoiesis is associated with adverse outcomes in multiple myeloma patients undergoing transplant.

Multiple myeloma (MM) is a plasma-cell neoplasm that is treated with high-dose chemotherapy, autologous stem cell transplant (ASCT) and long-term immunomodulatory drug (IMiD) maintenance. The presence of somatic mutations in the peripheral blood is termed clonal hematopoiesis of indeterminate potential (CHIP) and is associated with adverse outcomes. Targeted sequencing of the stem cell product from 629 MM patients treated by ASCT at the Dana-Farber Cancer Institute (2003-2011) detects CHIP in 136/629 patients (21.6%). The most commonly mutated genes are DNMT3A, TET2, TP53, ASXL1 and PPM1D. Twenty-one from fifty-six patients (3.3%) receiving first-line IMiD maintenance develop a therapy-related myeloid neoplasm (TMN). However, regardless of CHIP status, the use of IMiD maintenance associates with improved PFS and OS. In those not receiving IMiD maintenance, CHIP is associated with decreased overall survival (OS) (HR:1.34, p = 0.02) and progression free survival (PFS) (HR:1.45, p < 0.001) due to an increase in MM progression.

Good Citizenship Made Easy: A Step-by-Step Guide to Submitting RNA-Seq Data to NCBI.

The analysis of transcriptome data from non-model organisms contributes to our understanding of diverse aspects of evolutionary biology, including developmental processes, speciation, adaptation, and extinction. Underlying this diversity is one shared feature, the generation of enormous amounts of sequence data. Data availability requirements in most journals oblige researchers to make their raw transcriptome data publicly available, and the databases housed at the National Center for Biotechnology Information (NCBI) are a popular choice for data deposition. Unfortunately, the successful submission of raw sequences to the Sequence Read Archive (SRA) and transcriptome assemblies to the Transcriptome Shotgun Assembly (TSA) can be challenging for novice users, significantly delaying data availability and publication. Here we present two comprehensive protocols for submitting RNA-Seq data to NCBI databases, accompanied by an easy-to-use website that facilitates the timely submission of data by researchers of any experience level. © 2018 by John Wiley & Sons, Inc.

Comparative genomics explains the evolutionary success of reef-forming corals.

Transcriptome and genome data from twenty stony coral species and a selection of reference bilaterians were studied to elucidate coral evolutionary history. We identified genes that encode the proteins responsible for the precipitation and aggregation of the aragonite skeleton on which the organisms live, and revealed a network of environmental sensors that coordinate responses of the host animals to temperature, light, and pH. Furthermore, we describe a variety of stress-related pathways, including apoptotic pathways that allow the host animals to detoxify reactive oxygen and nitrogen species that are generated by their intracellular photosynthetic symbionts, and determine the fate of corals under environmental stress. Some of these genes arose through horizontal gene transfer and comprise at least 0.2% of the animal gene inventory. Our analysis elucidates the evolutionary strategies that have allowed symbiotic corals to adapt and thrive for hundreds of millions of years.

Adaptive Evolution of Eel Fluorescent Proteins from Fatty Acid Binding Proteins Produces Bright Fluorescence in the Marine Environment.

We report the identification and characterization of two new members of a family of bilirubin-inducible fluorescent proteins (FPs) from marine chlopsid eels and demonstrate a key region of the sequence that serves as an evolutionary switch from non-fluorescent to fluorescent fatty acid-binding proteins (FABPs). Using transcriptomic analysis of two species of brightly fluorescent Kaupichthys eels (Kaupichthys hyoproroides and Kaupichthys n. sp.), two new FPs were identified, cloned and characterized (Chlopsid FP I and Chlopsid FP II). We then performed phylogenetic analysis on 210 FABPs, spanning 16 vertebrate orders, and including 163 vertebrate taxa. We show that the fluorescent FPs diverged as a protein family and are the sister group to brain FABPs. Our results indicate that the evolution of this family involved at least three gene duplication events. We show that fluorescent FABPs possess a unique, conserved tripeptide Gly-Pro-Pro sequence motif, which is not found in non-fluorescent fatty acid binding proteins. This motif arose from a duplication event of the FABP brain isoforms and was under strong purifying selection, leading to the classification of this new FP family. Residues adjacent to the motif are under strong positive selection, suggesting a further refinement of the eel protein's fluorescent properties. We present a phylogenetic reconstruction of this emerging FP family and describe additional fluorescent FABP members from groups of distantly related eels. The elucidation of this class of fish FPs with diverse properties provides new templates for the development of protein-based fluorescent tools. The evolutionary adaptation from fatty acid-binding proteins to fluorescent fatty acid-binding proteins raises intrigue as to the functional role of bright green fluorescence in this cryptic genus of reclusive eels that inhabit a blue, nearly monochromatic, marine environment.