Generation Transcriptome Sequencing Research Articles

Transcriptome analysis and exploration of genes involved in the biosynthesis of secoiridoids in Gentiana rhodantha.

Gentiana rhodantha is a medicinally important perennial herb used as traditional Chinese and ethnic medicines. Secoiridoids are one of the major bioactive compounds in G. rhodantha. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from four organs (root, leaf, stem and flower), followed by the de novo sequence assembly. We verified 8-HGO (8-hydroxygeraniol oxidoreductase), which may encode key enzymes of the secoiridoid biosynthesis by qRT-PCR. The mangiferin, swertiamarin and loganic acid contents in root, stem, leaf, and flower were determined by HPLC. The results showed that there were 47,871 unigenes with an average length of 1,107.38 bp. Among them, 1,422 unigenes were involved in 25 standard secondary metabolism-related pathways in the KEGG database. Furthermore, we found that 1,005 unigenes can be divided into 66 transcription factor (TF) families, with no family members exhibiting significant organ-specificity. There were 54 unigenes in G. rhodantha that encoded 17 key enzymes of the secoiridoid biosynthetic pathway. The qRT-PCR of the 8-HGO and HPLC results showed that the relative expression and the mangiferin, swertiamarin, and loganic acid contents of the aerial parts were higher than in the root. Six types of SSR were identified by SSR analysis of unigenes: mono-nucleoside repeat SSR, di-nucleoside repeat SSR, tri-nucleoside repeat SSR, tetra-nucleoside repeat SSR, penta-nucleoside repeat SSR, and hexa-nucleoside repeat SSR. This report not only enriches the Gentiana transcriptome database but helps further study the function and regulation of active component biosynthesis of G. rhodantha.

Chapter 17 - Neonatal and infantile hypotonia

The underlying etiology of neonatal and infantile hypotonia can be divided into primary peripheral and central nervous system and acquired or genetic disorders. The approach to identifying the likeliest cause of hypotonia begins with a bedside assessment followed by a careful review of the birth history and early development and family pedigree and obtaining available genetic studies and age- and disease-appropriate laboratory investigations. Until about a decade ago, the main goal was to identify the clinical signs and a battery of basic investigations including electrophysiology to confirm or exclude a given neuromuscular disorder, however the availability of whole-exome sequencing and next generation sequencing and transcriptome sequencing has simplified the identification of specific underlying genetic defect and improved the accuracy of diagnosis in many related Mendelian disorders.

Transcriptome analysis reveals the spatial and temporal differentiation of gene expression in the sporophyte of Undaria pinnatifida

Multicellular organisms are differentiated into specialized organs having specific functions. Kelps, morphologically diverse brown algae of the order Laminariales, exhibit specialized morphological features with a holdfast, a stipe, and a blade. However, little is known about the expression profiles and their biological functions from the different tissues. Among the kelp species, the edible Pacific kelp, Undaria pinnatifida is an ecologically and economically important species native to the Northeastern Asia. It presents the typical organization of kelps but is characterized by the formation of its sporophyll, the sexually reproductive tissue, on the stipe. To study their respective gene expression profiles, we generated transcriptome sequences from four tissues; the sporophyll, meristem, healthy blade, and decaying apex blade were sampled at the three time points in the development of the sporophyte. We demonstrate an up-regulation of genes related to photosynthetic genes in the healthy blade, of genes related to cell proliferation in the meristem and of genes related to flagella development in the sporophyll. Overall, the study revealed clear spatial and temporal variation in the expression patterns within and between the different tissues that were well correlated with the biological functions in the different tissues.

Tissue-specific transcriptome profiles identify functional differences key to understanding whole plant response to life in variable salinity.

ABSTRACTPlants endure environmental stressors via adaptation and phenotypic plasticity. Studying these mechanisms in seagrasses is extremely relevant as they are important primary producers and functionally significant carbon sinks. These mechanisms are not well understood at the tissue level in seagrasses. Using RNA-seq, we generated transcriptome sequences from tissue of leaf, basal leaf meristem and root organs of Posidonia australis, establishing baseline in situ transcriptomic profiles for tissues across a salinity gradient. Samples were collected from four P. australis meadows growing in Shark Bay, Western Australia. Analysis of gene expression showed significant differences between tissue types, with more variation among leaves than meristem or roots. Gene ontology enrichment analysis showed the differences were largely due to the role of photosynthesis, plant growth and nutrient absorption in leaf, meristem and root organs, respectively. Differential gene expression of leaf and meristem showed upregulation of salinity regulation processes in higher salinity meadows. Our study highlights the importance of considering leaf meristem tissue when evaluating whole-plant responses to environmental change.This article has an associated First Person interview with the first author of the paper.

Differential expression of somatostatin receptor (SSTR) subtypes across a spectrum of neuroendocrine neoplasms (NENs).

3071 Background: Targeted therapy of NENs based on the presence of SSTRs fills a unique niche in tumor biology and clinical treatment of patients with solid tumors. SSTRs have multiple isoforms and are collectively expressed in the majority of NENs. However, subtypes are still not routinely tested and thus not assessed for clinical decision-making, especially for patients meriting consideration of targeted radionucleotide therapy. Clarifying the landscape of SSTR subtypes using molecular techniques more sensitive than immunohistochemistry (IHC)-the standard of testing, and identifying associated genomic biomarkers that differ between them, will pave the way for more sophisticated decision-making in the future. Additionally, leveraging transcriptomics to better assess mitotic markers such as Ki-67 to assess tumor grade, would increase diagnostic accuracy. Here we provide initial validation across a spectrum of NENs. Methods:1595 NENs were analyzed using Next Generation Sequencing (592 gene panel, NextSeq), Whole Exome and Transcriptome Sequencing (NovaSeq), and IHC at Caris Life Sciences (Phoenix, AZ). Significance was determined using chi-square, Fisher-Exact or Mann-Whitney U and p-adjusted for multiple comparisons (q<0.05) where applicable. Results: In a subset of 492 NENs with accompanying tumor grading information, a median MKI67 (gene encoding Ki-67) TPM value of 2.27 for low-grade (LG-), and 38.7 for high-grade NENs (HG-NENs) was observed (q<0.05). Using ROC curve analysis, a threshold of MKI67 expression (13.4375 TPM) differentiated LG- from HG-NENs, with a true positive rate of 86.84%, a false positive rate of 11.9% and an AUC of 95% and was subsequently applied to the entire cohort to infer HG/LG. Compared to HG-NENs (n = 862), LG-NENs(n = 733) expressed higher levels of SSTR 1(3.5-fold),2 (2.9-fold) and 5 (1.67-fold) and lower levels of SSTR4 (0.28-fold)(q<0.05). Further, the expression of SSTRs 3 and 4 in HG-NEN (rs= 0.63) and SSTRs 1 and 2 in LG-NENs (rs= 0.64) were positively correlated. Overall, the prevalence of altered TP53, RB1, PIK3CA, APC, KRAS was higher and MEN1 was lower in HG-vs LG-NENs (q<0.05). For each SSTR subtype, we established high and low cohorts based on median expressions. In LG-NENs, increased alterations in TP53 and RB1 were associated with increased expression of SSTRs 1 and 2 and reduced expression of SSTRs 3 and 4. In HG-NENs, increased alterations in APC were associated with increased expression of SSTR 1 and 4 and reduced expression of SSTRs 3 and 5. Additional subtype- and grade-specific alterations were also observed. Conclusions: This study provides evidence that WTS and NGS can be leveraged to predict grade of NENs and define characteristic differences in the genomic landscape across SSTR subtypes in HG and LG NENs. Incorporating the molecular profiling of NENs can thus aid in advancing the development of more tailored therapeutic strategies.

Identification of candidate genes from androgenic gland in Macrobrachium nipponense regulated by eyestalk ablation

The eyestalk of crustaceans, such as Macrobrachium nipponense, contains many neurosecretory hormones affecting the process of reproduction, molting, metabolism of glucose, and other functions. In this study, important metabolic pathways and candidate genes involved in male sexual development were selected from M. nipponense. The methodology involved performing long-read and next generation transcriptome sequencing of genes from the androgenic gland after eyestalk ablation. qPCR analysis revealed that the mRNA expression of Mn-IAG was significantly increased after ablation of both the single-side (SS) and double-side (DS) eyestalk, compared with the control group (CG). The long-read transcriptome generated 49,840 non-redundant transcripts. A total of 1319, 2092 and 4351 differentially expressed genes (DEGs) were identified between CG versus SS, SS versus DS and CG versus DS, respectively. These data indicated that ablation of the double-sided eyestalk played stronger regulatory roles than the single-side ablation on male sexual development in M. nipponense. This was consistent with the qPCR analysis. Cell Cycle, Cellular Senescence, Oxidative Phosphorylation, Glycolysis/Gluconeogenesis and Steroid Hormone Biosynthesis were the primary enriched metabolic pathways in all three comparisons, and the important genes from these metabolic pathways were also selected. qPCR permitted secondary confirmation of ten DEGs identified through RNA-seq. RNAi-mediated silencing analyses of Hydroxysteroid dehydrogenase like 1 (HSDL1) revealed that HSDL1 has a positive regulatory effect on testes development. This study provides valuable insight into male sexual development in M. nipponense, including metabolic pathways and genes, paving the way for advanced studies on male sexual development in this species and in other crustaceans.

Abstract 2137: Differential gene expression identifies a transcriptional regulatory network involving ESR1and PITX1 in invasive epithelial ovarian cancer

Abstract The heterogeneous subtypes and stages of epithelial ovarian cancer (EOC) differ in invasive processes, but the molecular factors that drive aggressive behavior remain nebulous. In this study, we analyzed two types of EOC at opposite ends of the invasiveness spectrum, high-grade serous ovarian cancers (HGSOCs) and serous borderline tumors (SBTs). We identified differentially expressed genes to guide follow on studies of regulatory networks. Initially, in a discovery sample set, we generated transcriptome sequences to identify 11 differentially expressed genes in SBTs and HGSOCs, and a nonoverlapping set of 17 differentially expressed genes stage II and stage III HGSOCs. We validated the classification performance of these gene expression signatures on a validation set of SBT and HGSOC with classification rates of 94 and 98% using the 11-gene and 17-gene sets, respectively. Transcription factor binding sites enriched in the differentially expressed gene promoter regions were examined by Chip-seq analysis in ovarian cancer and other cell lines. We identified a suite of transcription factors implicated in the differential regulation associated with invasiveness, including ERα, RARA, PITX1, FOXA1/FOXF1, and BHLHE41/E40. These data implicate estrogen responsive gene regulatory networks in ovarian cancers and suggest that dynamic regulatory interactions determine malignant EOC outcomes. Citation Format: Yichao Li, Sushil Kumar Jaiswal, Rupleen Kaur, Dana Alsaadi, Xiaoyu Liang, Frank Drews, Julie A. DeLoia, Thomas Krivak, Hanna M. Petrykowska, Valer Gotea, Lonnie Welch, Elnitski Laura. Differential gene expression identifies a transcriptional regulatory network involving ESR1and PITX1 in invasive epithelial ovarian cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2137.

Candidate genes involved in the biosynthesis of lignan in Schisandra chinensis fruit based on transcriptome and metabolomes analysis.

Schisandra chinensis Turcz. (Baill.) is a plant species with fruits that have been well known in Far Eastern medicine for a long time. It has traditionally been used as a stimulating and fortifying agent in cases of physical exhaustion and to inhibit fatigue. The major bioactive compounds found in S. chinensis are lignans with a dibenzocyclooctadiene skeleton, but little is known about their biosynthesis in plants. S. chinensis is the ideal medicinal plant for studying the biosynthesis of lignans, especially the dibenzocyclooctadiene skeleton. Genomic information for this important herbal plant is unavailable. To better understand the lignan biosynthesis pathway, we generated transcriptome sequences from the fruit during ripening and performed de novo sequence assembly, yielding 136 843 unique transcripts with N50 of 1778 bp. Putative functions could be assigned to 41 824 transcripts (51.57%) based on BLAST searches against annotation databases including GO (Gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes). Furthermore, 22 candidate cytochrome P450 genes and 15 candidate dirigent proteins genes that were most likely involved in the lignan biosynthesis pathway were discovered based on transcriptome sequencing of S. chinensis. The genomic data obtained from S. chinensis, especially the identification of putative genes involved in the lignan biosynthesis pathway, will facilitate our understanding of lignan biosynthesis at the molecular level. The lignan metabolite profiles were analyzed by metabolomes, the accumulation patterns of 30 metabolites involved in the lignan pathway were studied. Co-expression network of lignan contents and transcriptional changes showed 355 strong correlations (correlation coefficient, R2 > 0.9) between 21 compounds and 153 transcripts. Furthermore, the comprehensive analysis and characterization of the genes involved in lignan pathways and the metabolite profiles of lignans are expected to provide better insight regarding the diversity of the chemical composition, synthetic characteristics, and regulatory mechanisms of this medical herb.

Deep sequencing and transcriptome analyses to identify genes involved in iridoid biosynthesis in the medicinal plant Valeriana jatamansi Jones

Valeriana jatamansi Jones is utilized for medicinal purposes in China, and is also an important substitute for European Valeriana officinalis. The major active principles are generally called valepotriates, which belong to iridoids compounds. To better understand the iridoid biosynthesis pathway in V. jatamansi, we generated transcriptome sequences from the leaf and root tissues, and performed de novo sequence assembly, a total of 183,524,060 transcripts and 61,876 unigenes for V. jatamansi were obtained from 13.28 Gb clean reads. 56,641 unigenes were annotated by public databases, while 5,235 unigenes remained unannotated. Different unigenes in V. jatamansi were identified by MISA analysis, and 5,195 unigenes containing Simple Sequence Repeat (SSR) were identified. When examining the annotation of transcriptome contigs against the KEGG database, we identified 24 unigenes that could be classified into 24 enzyme categories associated with three metabolic pathways leading to iridoid biosynthesis, 6 genes of MVA pathways, 9 genes of MEP pathways and 9 genes of iridoids pathways. We selected 9 genes encoding key enzymes in the iridoid pathway of V. jatamansi to examine their organ specificity of expression using quantitative real-time PCR (qPCR). In conclusion, we generated a comprehensive transcriptome assembly representing the gene space in V. jatamansi, and the genomic dataset and analyses presented here lay the foundation for further research on this important medicinal plant.

Nuclear loci developed from multiple transcriptomes yield high resolution in phylogeny of scaly tree ferns (Cyatheaceae) from China and Vietnam

The scaly tree ferns, Cyatheaceae, form a large natural group of ca. 640 species. They comprise an ideal model for studying the biogeography of plants due to their broad distribution across the tropical to south-temperate areas of the world. However, tracing the biogeographic history of this highly diversified group has been hampered by unresolved evolutionary relationships among the major clades. Here, we generated transcriptome sequences of five species in three genera of Cyatheaceae (Alsophila, Gymnosphaera, and Sphaeropteris) and used them to search for single-copy nuclear loci for phylogenetic reconstruction. We identified a total of 818 candidate single-copy loci across multiple Cyatheaceae species. To test their phylogenetic utility, we further obtained sequence data of 12 of these loci for 76 samples representing all 13 known species of scaly tree ferns in China and Vietnam. Phylogenetic analyses based on multispecies coalescent and, alternatively, concatenation models yielded congruent results with high resolution. Additionally, we used the 12 loci to identify genetic signals of hybridization. Overall, our results demonstrated that multiple, single-copy loci are informative and efficient tools for phylogenetic or evolutionary studies of scaly tree ferns.

Identification of candidate genes involved in the response to different abiotic stresses in potato (Solanum tuberosum L.)

Plants growing in natural habitats are exposed to multiple environmental stresses resulting from abiotic factors such as heat, drought, and cold, which have a significant impact on cultivated potato. We have evaluated in two Solanum tuberosum varieties, Soprano and Kondor, the adaptation to different abiotic stresses as heat, cold, and drought. For this purpose plants of both varieties were stressed, and when they showed symptoms, RNA extraction was carried out and a cDNA library for each sample was constructed. The objective of this study was to detect and analyse the genes involved in the responses to abiotic stresses in potato. The assay generated transcriptome sequences from both varieties, and a total of 5.579.655 reads and 8420 putative candidate genes were generated. 4.027 of the candidate genes were polymorphic and presented a different number of patterns defined by a varying number of SNPs. Many of the generated candidate genes showed differential expression, since the candidate gene was present in the stressed plant, but not in the control plant. The application of this methodology allows us to detect numerous candidate genes or specific alleles/allele combinations, which are differentially expressed in specific samples after the application of different abiotic stresses. This will be useful to identify superior alleles which can be used in Marker Assisted Selection for resistance and tolerance to abiotic stresses.

Comparative transcriptome analysis of roots, stems and leaves of Isodon amethystoides reveals candidate genes involved in Wangzaozins biosynthesis

BackgroundIsodon amethystoides (Ben-th) Cy Wu et Hsuan is an important traditional medicinal plant endowed with pharmacological properties effective in the treatment of various diseases, including pulmonary tuberculosis. The tetracyclic diterpenoids, Wangzaozins (Wangzaozin A, glaucocalyxin A, glaucocalyxin B), are the major bioactive compounds of I. amethystoides. However, the molecular information about the biosynthesis of these compounds still remains unclear.ResultsAn examination of the accumulated levels of Wangzaozins in I. amethystoides revealed considerable variations in the root, stem, and leaf tissues of this plant, indicating possible differences in metabolite biosynthesis and accumulation among various tissues. To better elucidate the tetracyclic diterpenoid biosynthesis pathway, we generated transcriptome sequences from the root, stem, and leaf tissues, and performed de novo sequence assembly, yielding 230,974 transcripts and 114,488 unigenes, with average N50 lengths of 1914 and 1241 bp, respectively. Putative functions could be assigned to 73,693 transcripts (31.9%) based on BLAST searches against annotation databases, including GO, KEGG, Swiss-Prot, NR, and Pfam. Moreover, the candidate genes involving in the diterpenoid biosynthesis, such as CPS, KSL, were also analyzed. The expression profiles of eight transcripts, involving the tetracyclic diterpenoid biosynthesis, were validated in different I. amethystoides tissues by qRT-PCR, unraveling the gene expression profile of the pathway. The differential expressions of ISPD, ISPF and ISPH (MEP pathway), and IaCPS and IaKSL (diterpenoid pathway) candidate genes in leaves and roots, may contribute to the high accumulation of Wangzaozins in I. amethystoides leaves.ConclusionThe genomic dataset and analyses reported here lay the foundations for further research on this important medicinal plant.

Transcriptome Analysis of Bronchoalveolar Lavage Fluid From Children With Mycoplasma pneumoniae Pneumonia Reveals Natural Killer and T Cell-Proliferation Responses.

Mycoplasma pneumoniae pneumonia (MPP) is one of the most common community-acquired pneumonia; this study is to explore the immune-pathogenesis of children MPP. Next-generation transcriptome sequencing was performed on the bronchoalveolar lavage fluid cells from six children with MPP and three children with foreign body aspiration as control. Some of the results had been validated by quantitative real-time PCR in an expanded group of children. Results revealed 810 differentially expressed genes in MPP group comparing to control group, of which 412 genes including RLTPR, CARD11 and RASAL3 were upregulated. These upregulated genes were mainly enriched in mononuclear cell proliferation and signaling biological processes. Kyoto encyclopedia of genes and genomes pathway analysis revealed that hematopoietic cell linage pathway, natural killer cell-mediated cytotoxicity pathway, and T cell receptor signaling pathway were significantly upregulated in MPP children. In addition, significant alternative splicing events were found in GNLY and SLC11A1 genes, which may cause the differential expressions of these genes. Our results suggest that NK and CD8+ T cells are over activated and proliferated in MPP children; the upregulated IFNγ, PRF1, GZMB, FASL, and GNLY may play important roles in the pathogenesis of children MPP.

The Inhibitive Effect of Salidroside on Apoptosis of EA.hy926 Cells Induced by H2O2 under Simulated Microgravity

Microgravity or simulated microgravity(SM) can induce cardiovascular dysfunction in which vascular endothelial cells play an important role. Studies showed that oxidative stress was common in vascular tissues under simulated microgravity and was also increased in vascular endothelial cells during spaceflight. The purpose of this experiment is to study whether vascular endothelial cells have different biological response to oxidative stress under simulated microgravity with new generation transcriptome sequencing technology and discuss the cell protection effect of salidroside. EA.hy926 cells, a fusion of HUVECs and the lung carcinoma cell line A549, were randomly divided into control group(Con72) and clinostat group(Z72), and after clinostating for 72 hours, the two groups were randomly given the following treatment at the same time: stimulated by H2O2, stimulation by H2O2 and salidroside intervention. Apoptosis was measured by flow cytometry and transcriptome sequencing was performed by new generation sequencing technology, RT‐qPCR was used to verify the differentially expressed genes and detect the effect of salidroside on the expression of these genes. The results showed that there were obvious apoptosis in EA.hy926 cells treated by H2O2 under simulated microgravity, which was significantly inhibited by salidroside. There were six significant differentially expressed genes screened by transcriptome sequencing, including BCL‐2A1, FAM196B, TMEM158, PPP1R16B, CXCL8, PPP1R3B. RT‐qPCR results showed that BCL‐2A1 and FAM196B were consistent with the transciptome sequencing results, meanwhile, BCL2A1 and FAM196B genes were significantly increased after salidroside intervention in H2O2‐treated EA.hy926 cells under simulated microgravity. These results showed that apoptosis was induced significantly in H2O2‐treated EA.hy926 cells under simulated microgravity, which can be inhibited by salidroside, the active ingredient of Rhodiola rosea, through the regulation of BCL‐2A1 and FAM196B gene expression.Support or Funding InformationFunded by State Key Laboratory of Space Medicine Fundamentals and Application, China Astronaut Research and Training Center, grant No.SMFA14B01This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Establishment of molecular genetic approaches to study gene expression and function in an invasive hemipteran, Halyomorpha halys

Hemiptera is a large clade of insects understudied in terms of developmental biology. Halyomorpha halys, the Brown Marmorated Stink Bug (BMSB, referred to throughout as H. halys), is an invasive hemipteran pest of the mid-Atlantic region of the USA that has rapidly spread to other regions in recent years, devastating a wide range of crops using a piercing and sucking mechanism. Its phylogenetic position, polyphagous habits, and rapid spread in the USA suggested that H. halys would be an ideal system to broaden our knowledge of developmental mechanisms in insects. We and others previously generated transcriptome sequences from different life stages of this insect. Here, we describe tools to examine gene expression patterns in whole-mount H. halys embryos and to test the response of H. halys to RNA interference (RNAi). We show that spatial and temporal patterns of gene expression in H. halys can be effectively monitored by both immunostaining and in situ hybridization. We also show that delivery of dsRNA to adult females knocks down gene function in offspring, using the homeotic gene Sex combs reduced (Scr). Knockdown of Hh-Scr resulted in dramatic malformations of the mouthparts, demonstrating for the first time that RNAi is effective in this species. Our results suggest that, despite difficulties with long-term laboratory culture of H. halys, this species shows promise as a developmental system.

Deep sequencing and transcriptome analyses to identify genes involved in secoiridoid biosynthesis in the Tibetan medicinal plant Swertia mussotii

Swertia mussotii Franch. is an important traditional Tibetan medicinal plant with pharmacological properties effective in the treatment of various ailments including hepatitis. Secoiridoids are the major bioactive compounds in S. mussotii. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from the root, leaf, stem, and flower tissues, and performed de novo sequence assembly, yielding 98,613 unique transcripts with an N50 of 1,085 bp. Putative functions could be assigned to 35,029 transcripts (35.52%) based on BLAST searches against annotation databases including GO and KEGG. The expression profiles of 39 candidate transcripts encoding the key enzymes for secoiridoid biosynthesis were examined in different S. mussotii tissues, validated by qRT-PCR, and compared with the homologous genes from S. japonica, a species in the same family, unveiling the gene expression, regulation, and conservation of the pathway. The examination of the accumulated levels of three bioactive compounds, sweroside, swertiamarin, and gentiopicroside, revealed their considerable variations in different tissues, with no significant correlation with the expression profiles of key genes in the pathway, suggesting complex biological behaviours in the coordination of metabolite biosynthesis and accumulation. The genomic dataset and analyses presented here lay the foundation for further research on this important medicinal plant.

Comparative transcriptome analysis reveals candidate genes for the biosynthesis of natural insecticide in Tanacetum cinerariifolium

BackgroundPyrethrins are monoterpenoids and consist of either a chrysanthemic acid or pyrethric acid with a rethrolone moiety. Natural pyrethrins are safe and eco-friendly while possessing strong insecticidal properties. Despite such advantages of commercial value coming with the eco-friendly tag, most enzymes/genes involved in the pyrethrin biosynthesis pathway remain unidentified and uncharacterized. Since the flowers of Tanacetum cinerariifolium are rich in major pyrethrins, next generation transcriptome sequencing was undertaken to compare the flowers and the leaves of the plant de novo to identify differentially expressed transcripts and ascertain which among them might be involved in and responsible for the differential accumulation of pyrethrins in T. cinerariifolium flowers.ResultsIn this first tissue specific transcriptome analysis of the non-model plant T. cinerariifolium, a total of 23,200,000 and 28,500,110 high quality Illumina next generation sequence reads, with a length of 101 bp, were generated for the flower and leaf tissue respectively. After functional enrichment analysis and GO based annotation using public protein databases such as UniRef, PFAM, SMART, KEGG and NR, 4443 and 8901 unigenes were identified in the flower and leaf tissue respectively. These could be assigned to 13344 KEGG pathways and the pyrethrin biosynthesis contextualized. The 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway was involved in the biosynthesis of acid moiety of pyrethrin and this pathway predominated in the flowers as compared to the leaves. However, enzymes related to oxylipin biosynthesis were found predominantly in the leaf tissue, which suggested that major steps of pyrethrin biosynthesis occurred in the flowers.ConclusionsTranscriptome comparison between the flower and leaf tissue of T. cinerariifolium provided an elaborate list of tissue specific transcripts that was useful in elucidating the differences in the expression of the biosynthetic pathways leading to differential presence of pyrethrin in the flowers. The information generated on genes, pathways and markers related to pyrethrin biosynthesis in this study will be helpful in enhancing the production of these useful compounds for value added breeding programs. Related proteome comparison to overlay our transcriptome comparison can generate more relevant information to better understand flower specific accumulation of secondary metabolites in general and pyrethrin accumulation in particular.

Assessing the Gene Content of the Megagenome: Sugar Pine (Pinus lambertiana).

Sugar pine (Pinus lambertiana Douglas) is within the subgenus Strobus with an estimated genome size of 31 Gbp. Transcriptomic resources are of particular interest in conifers due to the challenges presented in their megagenomes for gene identification. In this study, we present the first comprehensive survey of the P. lambertiana transcriptome through deep sequencing of a variety of tissue types to generate more than 2.5 billion short reads. Third generation, long reads generated through PacBio Iso-Seq have been included for the first time in conifers to combat the challenges associated with de novo transcriptome assembly. A technology comparison is provided here to contribute to the otherwise scarce comparisons of second and third generation transcriptome sequencing approaches in plant species. In addition, the transcriptome reference was essential for gene model identification and quality assessment in the parallel project responsible for sequencing and assembly of the entire genome. In this study, the transcriptomic data were also used to address questions surrounding lineage-specific Dicer-like proteins in conifers. These proteins play a role in the control of transposable element proliferation and the related genome expansion in conifers.

The Clonal Architecture of CXCR4mutations in Waldenstrom's Macroglobulinemia Shows Highly Variable Subclonal Distribution, and Multiple Mutations within Individual Patients Indicative of Targeted Genomic Instability

Background: Whole genome sequencing (WGS) identified activating CXCR4WHIM somatic mutations in nearly 30% of patients with Waldenstrom's Macroglobulinemia (WM) (Blood 123(11):1637-46). Both nonsense and frameshift CXCR4WHIM mutations occur in WM, with over 30 different types of mutations described within the regulatory carboxyl-terminal domain of CXCR4. CXCR4WHIM mutations almost always occur with activating MYD88 mutations, and impact both disease presentation and treatment outcome (Blood 123(18):2791-6; NEJM 372(15):1430-40.). The clonal architecture of CXCR4WHIM mutations relative to MYD88 mutations and their role in disease evolution remains to be clarified.Methods: We used Sanger sequencing and highly sensitive AS-PCR assays that we developed for the most common CXCR4WHIM mutations (S338X C>A and C>G) to evaluate for CXCR4WHIM mutations. In conjunction with an AS-PCR MYD88L265P assay that we previously developed (Leukemia 28(8):1698-707), we also profiled tumor samples for MYD88L265P and CXCR4S338X mutations in 164 WM, 12 IgM MGUS, 20 MZL, 32 CLL, 14 MM, 7 non-IGM MGUS patients, and 32 healthy donors. Next generation transcriptome sequencing data was also performed for validation in select cases.Results: AS-PCR detected CXCR4S338X mutations in WM and IgM MGUS patients not revealed by Sanger sequencing. By combined AS-PCR and Sanger sequencing, CXCR4WHIM mutations were identified in 44/102 (43%), 21/62 (34%), 2/12 (17%), and 1/20 (5%) untreated WM, previously treated WM, IgM MGUS and MZL patients, respectively, but not in CLL, MM, non-IGM MGUS patients or healthy donors. Cancer cell fraction analysis in WM and IgM MGUS patients showed CXCR4S338X mutations were primarily subclonal, with highly variable clonal distribution (median 35.1%, range 1.2%-97.5%; Figure 1). Sanger sequencing identified 3 patients with multiple CXCR4 mutations, which were shown to be compound heterozygous by TA cloning and sequencing of at least 40 clones. The addition of AS-PCR to the Sanger sequencing results also revealed multiple CXCR4WHIM mutations in many individual patients that included homozygous and compound heterozygous mutations that were validated by next generation sequencing that offered a median of 5,819 (range 5,217-12,235) reads that overlapped the mutated loci.Conclusions: Taken together, these findings show that CXCR4WHIM mutations are more common in WM patients than previously revealed by WGS or Sanger sequencing. Moreover, CXCR4 mutations are primarily subclonal supporting their acquisition after MYD88L265P in WM oncogenesis. The exclusive finding of frameshift and nonsense but not missense variants within the carboxyl-terminal domain of CXCR4 suggests that significant selection pressures exist for activating mutations within the WM clone. Lastly, multiple CXCR4WHIM mutations are common in WM patients indicative of targeted genomic instability within the carboxyl-terminal domain of CXCR4. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Transcriptome analysis of the adult rumen fluke Paramphistomum cervi following next generation sequencing

Rumen flukes are parasitic trematodes (Platyhelminthes: Digenea) of major socioeconomic importance in many countries. Key representatives, such as Paramphistomum cervi, can cause “Rumen fluke disease” or paramphistomosis and undermine economic animal productivity and welfare. P. cervi is primarily a problem in sheep, goat and buffalo production as a consequence of reduced weight gain and milk production, clinical disease or death. Recent technological advances in genomics and bioinformatics now provide unique opportunities for the identification and pre-validation of drug targets and vaccines through improved understanding of the biology of pathogens such as P. cervi and their relationship with their hosts at the molecular level. Here, we report next generation transcriptome sequencing analysis for P. cervi. RNAseq libraries were generated from RNA extracted from 15 adult P. cervi parasites sampled from each of three different host species (sheep, goat and buffalo) and a reference transcriptome was generated by assembly of all Ion Torrent PGM sequencing data. Raw reads (7,433,721 in total) were initially filtered for host nucleotide contamination and ribosomal RNAs and the remaining reads were assembled into 43,753 high confidence transcript contigs. In excess of 50% of the assembled transcripts were annotated with domain- or protein sequence similarity derived functional information. The reference adult P. cervi transcriptome will serve as a basis for future work on the biology of this important parasite. Using the widely investigated trematode virulence factor and vaccine candidate Cathepsin L as an example, the epitope GPISIAINA was found to be conserved in P. cervi isolated from three different host species supporting its candidacy for vaccine development and illustrating the utility of the adult P. cervi transcriptome.