Sequencing Of cDNA Libraries Research Articles

Identification of the LEA family members from Caragana korshinskii (Fabaceae) and functional characterization of CkLEA2-3 in response to abiotic stress in Arabidopsis

Plants are exposed to different abiotic stresses that affect growth, development, and productivity. Previous studies have shown that late embryogenesis abundant (LEA) proteins play important roles in plant adaptation to abiotic stresses. However, reports that demonstrate their biological functions, especially in Caragana korshinskii Kom., are still very limited. In this study, 26 LEA genes were identified from dehydration-treated suppressive subtractive hybridizations cDNA library and transcriptome sequencing data of C. korshinskii, and were classified into seven groups according to their structural features. Quantitative real-time PCR analysis revealed that the CkLEAs were induced by abscisic acid (ABA) and diverse abiotic stresses and widely expressed in various tissues. Moreover, overexpression of CkLEA2-3 in Arabidopsis thaliana L. Heynh resulted in enhanced tolerance to ABA treatment, and osmotic and salt stresses during seed germination. CkLEA2-3 overexpression lines also exhibited resistance to drought stress during seedling development. Taken together, our results indicate that CkLEA2-3 plays positive roles in conferring abiotic stress tolerance in Arabidopsis.

Fluorescence-Activated Cell Sorting-Based Isolation and Characterization of Neural Stem Cells from the Adult Zebrafish Telencephalon.

Adult mammalian brain, including humans, has rather limited addition of new neurons and poor regenerative capacity. In contrast, neural stem cells (NSC) with glial identity and neurogenesis are highly abundant throughout the adult zebrafish brain. Importantly, the activation of NSC and production of new neurons in response to injuries lead to the brain regeneration in zebrafish brain. Therefore, understanding of the molecular pathways regulating NSC behavior in response to injury is crucial in order to set the basis for experimental modification of these pathways in glial cells after injury in the mammalian brain and to elicit neuronal regeneration. Here, we describe the procedure that we successfully used to prospectively isolate NSCs from adult zebrafish telencephalon, extract RNA, and prepare cDNA libraries for next generation sequencing (NGS) and full transcriptome analysis as the first step toward understanding regulatory mechanisms leading to restorative neurogenesis in zebrafish. Moreover, we describe an alternative approach to analyze antigenic properties of NSC in the adult zebrafish brain using intracellular fluorescence activated cell sorting (FACS). We employ this method to analyze the number of proliferating NSCs positive for proliferating cell nuclear antigen (PCNA) in the prospectively isolated population of stem cells.

Strategies for Global RNA Sequencing of the Human Pathogen Neisseria gonorrhoeae.

Over the last several years transcriptomic analysis of bacterial pathogens has become easier and less expensive. This technique is used to determine expression levels for all genes of a particular species or collection of species under a given condition, including genes that are not yet known to exist. While transcriptomics can be a powerful tool to better understand bacterial regulatory responses to specific host environments, the experimental approach and data analysis must be performed correctly to ensure robust, accurate, and translational results. Here, we describe experimental protocols related to transcriptomic analysis of the sexually transmitted disease pathogen Neisseria gonorrhoeae. Methods are described for the extraction of high-quality RNA, examination of RNA to ensure quality, the generation of cDNA libraries for sequencing and the downstream analysis of raw sequencing data to determine gene expression levels. Much of this work can be carried out with equipment and reagents that are readily available, and the methods can be performed by a majority of laboratory groups. RNA-seq and transcriptomic analyses are set to become even more common in the coming years. The protocols described here will provide a standardized set of methods for applying this powerful technique to the study of N. gonorrhoeae under a variety of conditions.

포르말린 폭로에 따른 스쿠티카충의 유전체 발현 특성

Scuticociliate are one of the main parasite affecting the flounder aquaculture industry, causing significant economic losses in Korea. Aquatic formalin has been approved to treat against external parasitic infestation in flounder aquaculture pond. To investigate metabolic effect after treating the parasitic agent aquatic formalin in scuticociliate Miamiensis avidus, the transcriptome analysis was conducted. Through exposing with low concentration(10~50 ppm) of formalin for 96 hr, parasite trancriptomes were sequenced cDNA libraries in the GX Plex sequencer. A total of million read for control and formalin treated group were assembled in 187,210 and 177,959 contigs, respectively. After normalization and assembling with house keeping gene(β-actin), transcripts were assigned to biological processes and functions after annotation in Gene Ontology(GO). Specifically, 10 exclusive transcripts were up- or down-regulated for control and formalin treated group, respectively. Our results provide valuable genetic information for further analysis of the biological responses mechanism of aquatic formalin exposure in M. avidus.

Usage of Single-Cell RNA Sequencing to Unveil Immune Lymphoid Cell Precursors.

Single-cell RNA sequencing (scRNA-seq) is a powerful tool to study immune cells, which enables an unbiased way to discover novel cell populations, biological meaningful cellular heterogeneity, and cell lineage development trajectories. Advances in scRNA-seq technologies and computational data analysis have driven a revolution in our understanding of the immune system in health and disease. Technically, the key step for scRNA-seq analysis is making a high-quality cDNA library for sequencing. Here, we describe a plate-based protocol to prepare single-cell cDNA library of bone marrow innate lymphoid precursors for next generation sequencing-based transcriptome analysis.

First Transcriptome Analysis of Iranian Scorpion, Mesobuthus Eupeus Venom Gland.

Scorpions are generally an important source of bioactive components, including toxins and other small peptides as attractive molecules for new drug development. Mesobuthus eupeus, from medically important and widely distributed Buthidae family, is the most abundant species in Iran. Researchers are interesting on the gland of this scorpion due to the complexity of its venom. Here, we have analyzed the transcriptome based on expressed sequence tag (EST) database from the venom tissue of Iranian M. eupeus by constructing a cDNA library and subsequent Sanger sequencing of obtained inserts. Sixty-three unique transcripts were identified, which were grouped in different categories, including Toxins (44 transcripts), Cell Proteins (9 transcripts), Antimicrobial Peptides (4 transcripts) and Unknown Peptides (3 transcripts). The analysis of the ESTs revealed several new components categorized among various toxin families with effect on ion channels. Sequence analysis of a new precursor provides evidence to validate the first CaTxs from M. eupeus. The results are exploration of the diversity of precursors expressed of Iranian M. eupeus venom gland. We further described comparative analysis of venom components of Iranian M. eupeus with other sibling species.

Cloning and characterization of MdGST1 from red apple leaves

Glutathione S-transferase (GST) is involved in the downstream steps of the anthocyanin biosynthesis pathway in plants. However, the gene(s) encoding this enzyme have not been isolated from apple (Malus × domestica Borkh.) yet. We isolated a gene encoding GST from leaves of the red-fleshed apple variety ‘Royalty’ by full cDNA library sequencing and the 3′ rapid-amplification of cDNA ends method, and designated it MdGST1. In total, seven different MdGST1 transcripts were found. These had three different untranslated but identical protein-coding regions. Phylogenetic analysis showed that MdGST1 is a TT19-type GST, which is involved in anthocyanin transport. qRT-PCR analyses showed that the transcript level of MdGST1 was much higher in red leaves than in bagged or green leaves. When MdGST1 was introduced into a non-pigmented mutant of Arabidopsis, the transformants showed a visible purple phenotype in leaves and stems. Our results suggest that MdGST1 plays a role in anthocyanin biosynthesis. This information will help to improve the understanding of the mechanism of apple coloration.

Characterization of transcriptional responses mediated by benzo[a]pyrene stress in a new marine fish model of goby, Mugilogobius chulae.

Benzo[a]pyrene (BaP) is one of the most studied targets among polycyclic aromatic hydrocarbons (PAHs). Because of the complexity of the toxicity mechanism in BaP, little is known about the molecular mechanism at the level of transcription of BaP in marine fishes. The primary objective of this study was to investigate the molecular basis of the effects of BaP on marine fish, using Mugilogobius chulae (Smith 1932) as the model. A closed colony of M. chulae was used for the BaP toxicity test. Two fish liver samples per replicate from each group were excised and blended into one sample by pooling an equal amount of liver tissue. Total RNA of all samples was extracted separately. Equal quantities of total RNA from the three replicates of the two groups were pooled for sequencing. The sequencing cDNA libraries were sequenced using Illumina HiSeq 2000 system. Differentially expressed genes were detected with the DEGSeq R package. In total, 52,364,032 and 53,771,748 clean nucleotide reads were obtained in the control and BaP-exposed libraries, respectively, with N50 lengths of 1277 and 1288bp, respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed a significant enrichment of genes related to detoxification, transportation, and lipid metabolism. We also identified, for the first time, an association between endoplasmic reticulum dysfunction and lipid metabolism resulting from BaP exposure. Using quantitative real-time PCR, some effective molecular biomarkers for monitoring of BaP-polluted seawater were identified. The results demonstrate that BaP enhanced the expression of genes involved in detoxification in M. chulae and inhibited that of genes related to lipid metabolism, possibly by suppressing the expression of numerous ER-related genes involved in fat digestion and absorption.

WED 139 Gene expression regulation in CD4+ t-cell dysfunction in ms

BackgroundGenetic risk variants for complex disorders such as multiple sclerosis (MS) rarely encode protein. Growing evidence suggests risk variants regulate gene expression in specific disease-relevant cells by altering regulatory genetic sequences. This can be explored by correlating genotype at risk loci with expression of nearby genes (cis-regulation), with greatest sensitivity demonstrated in the most relevant cell types. Identifying genes under regulation by the MS risk variants could be a crucial step in identifying novel therapeutic targets.Aims1) Interrogate MS risk loci for cis-regulatory effects using allele-specific expression (ASE) analysis. 2) Establish a publicly available resource of genome-wide RNA-seq expression data in CSF-derived CD4 +T lymphocytes. 3) Correlate gene expression data with longitudinal clinical data to seek biomarkers of prognosis.Method and progressWe have extracted RNA from FACS-sorted CSF-derived CD4 +T lymphocytes from over 100 individuals. cDNA library preparation and sequencing for each sample is underway.Data analysisASE analysis makes use of a pre-selected transcribed SNP to differentiate maternal and paternal transcripts of each gene of interest. When an individual is heterozygous at an additional functional regulatory locus the ratio of maternal: paternal expression will deviate from 1:1. Each MS risk locus will be interrogated for this effect.

Identification of key genes and pathways involved in response to pain in goat and sheep by transcriptome sequencing

PurposeThis aim of this study was to investigate the key genes and pathways involved in the response to pain in goat and sheep by transcriptome sequencing.MethodsChronic pain was induced with the injection of the complete Freund’s adjuvant (CFA) in sheep and goats. The animals were divided into four groups: CFA-treated sheep, control sheep, CFA-treated goat, and control goat groups (n = 3 in each group). The dorsal root ganglions of these animals were isolated and used for the construction of a cDNA library and transcriptome sequencing. Differentially expressed genes (DEGs) were identified in CFA-induced sheep and goats and gene ontology (GO) enrichment analysis was performed.ResultsIn total, 1748 and 2441 DEGs were identified in CFA-treated goat and sheep, respectively. The DEGs identified in CFA-treated goats, such as C-C motif chemokine ligand 27 (CCL27), glutamate receptor 2 (GRIA2), and sodium voltage-gated channel alpha subunit 3 (SCN3A), were mainly enriched in GO functions associated with N-methyl-d-aspartate (NMDA) receptor, inflammatory response, and immune response. The DEGs identified in CFA-treated sheep, such as gamma-aminobutyric acid (GABA)-related DEGs (gamma-aminobutyric acid type A receptor gamma 3 subunit [GABRG3], GABRB2, and GABRB1), SCN9A, and transient receptor potential cation channel subfamily V member 1 (TRPV1), were mainly enriched in GO functions related to neuroactive ligand-receptor interaction, NMDA receptor, and defense response.ConclusionsOur data indicate that NMDA receptor, inflammatory response, and immune response as well as key DEGs such as CCL27, GRIA2, and SCN3A may regulate the process of pain response during chronic pain in goats. Neuroactive ligand-receptor interaction and NMDA receptor as well as GABA-related DEGs, SCN9A, and TRPV1 may modulate the process of response to pain in sheep. These DEGs may serve as drug targets for preventing chronic pain.

The Future of Cross-Linking and Immunoprecipitation (CLIP).

To understand the assembly and functional outcomes of protein-RNA regulation, it is crucial to precisely identify the positions of such interactions. Cross-linking and immunoprecipitation (CLIP) serves this purpose by exploiting covalent protein-RNA cross-linking and RNA fragmentation, along with a series of stringent purification and quality control steps to prepare complementary DNA (cDNA) libraries for sequencing. Here we describe the core steps of CLIP, its primary variations, and the approaches to data analysis. We present the application of CLIP to studies of specific cell types in genetically engineered mice and discuss the mechanistic and physiologic insights that have already been gained from studies using CLIP. We conclude by discussing the future opportunities for CLIP, including studies of human postmortem tissues from disease patients and controls, RNA epigenetic modifications, and RNA structure. These and other applications of CLIP will continue to unravel fundamental gene regulatory mechanisms while providing important biologic and clinically relevant insights.

Human plasma and serum extracellular small RNA reference profiles and their clinical utility

Circulating extracellular RNAs (exRNAs) have the potential to serve as biomarkers for a wide range of medical conditions. However, limitations in existing exRNA isolation methods and a lack of knowledge on parameters affecting exRNA variability in human samples may hinder their successful discovery and clinical implementation. Using combinations of denaturants, reducing agents, proteolysis, and revised organic extraction, we developed an automated, high-throughput approach for recovery of exRNAs and exDNA from the same biofluid sample. We applied this method to characterize exRNAs from 312 plasma and serum samples collected from 13 healthy volunteers at 12 time points over a 2-month period. Small RNA cDNA library sequencing identified nearly twofold increased epithelial-, muscle-, and neuroendocrine-cell-specific miRNAs in females, while fasting and hormonal cycle showed little effect. External standardization helped to detect quantitative differences in erythrocyte and platelet-specific miRNA contributions and in miRNA concentrations between biofluids. It also helped to identify a study participant with a unique exRNA phenotype featuring a miRNA signature of up to 20-fold elevated endocrine-cell-specific miRNAs and twofold elevated total miRNA concentrations stable for over 1 year. Collectively, these results demonstrate an efficient and quantitative method to discern exRNA phenotypes and suggest that plasma and serum RNA profiles are stable over months and can be routinely monitored in long-term clinical studies.

Using Fluidigm C1 to Generate Single-Cell Full-Length cDNA Libraries for mRNA Sequencing.

Single-cell RNA sequencing has evolved into a benchmark application to study cellular heterogeneity, advancing our understanding of cellular differentiation, disease progression, and gene regulation in a multitude of research areas. The generation of high-quality cDNA, an important step in the experimental workflow when generating sequence-ready libraries, is critical to maximizing data quality. Here we describe a strategy that uses a microfluidic device (i.e., the C1™ IFC) to synthesize full-length cDNA from single cells in a fully automated, nanoliter-scale format. The device also facilitates confirmation of the presence of a single, viable cell and recording of phenotypic information, quality control measures that are crucial for streamlining downstream data processing and enhancing overall data validity.

MicroRNA Expression Analysis: Next-Generation Sequencing.

MicroRNAs are a class of small noncoding RNAs that function as regulators involving in many biological processes. The evaluation of miRNAs and their targets has been aided by miRNA expression profiling studies including multiplex PCR, microarrays, and recent next-generation sequencing tools. Next-generation sequencing has enabled us to profile thousands of genes in a single experiment and overcome the background signal and cross-hybridization issues of microarrays. Next-generation sequencing also allows for the simultaneous confirmation of known miRNAs and discovery of new miRNAs, and significantly reduces costs while providing billions of nucleotide information within a single experiment. Here we describe a detailed procedure of generation of miRNA library for next-generation sequencing to increase the efficiency of adapter ligation and finally construct a more specific cDNA library for sequencing and analyses for miRNA expression profiling.

Individual Nucleotide Resolution UV Cross-Linking and Immunoprecipitation (iCLIP) to Determine Protein-RNA Interactions.

RNA-binding proteins (RBPs) interact with and determine the fate of many cellular RNA transcripts. In doing so they help direct many essential roles in cellular physiology, while their perturbed activity can contribute to disease etiology. In this chapter we detail a functional genomics approach, termed individual nucleotide resolution UV cross-linking and immunoprecipitation (iCLIP), that can determine the interactions of RBPs with their RNA targets in high throughput and at nucleotide resolution. iCLIP achieves this by exploiting UV-induced covalent cross-links formed between RBPs and their target RNAs to both purify the RBP-RNA complexes under stringent conditions, and to cause reverse transcription stalling that then identifies the direct cross-link sites in the high throughput sequenced cDNA libraries.

A Comparison of mRNA Sequencing with Random Primed and 3\u2032-Directed Libraries

Creating a cDNA library for deep mRNA sequencing (mRNAseq) is generally done by random priming, creating multiple sequencing fragments along each transcript. A 3′-end-focused library approach cannot detect differential splicing, but has potentially higher throughput at a lower cost, along with the ability to improve quantification by using transcript molecule counting with unique molecular identifiers (UMI) that correct PCR bias. Here, we compare an implementation of such a 3′-digital gene expression (3′-DGE) approach with “conventional” random primed mRNAseq. Given our particular datasets on cultured human cardiomyocyte cell lines, we find that, while conventional mRNAseq detects ~15% more genes and needs ~500,000 fewer reads per sample for equivalent statistical power, the resulting differentially expressed genes, biological conclusions, and gene signatures are highly concordant between two techniques. We also find good quantitative agreement at the level of individual genes between two techniques for both read counts and fold changes between given conditions. We conclude that, for high-throughput applications, the potential cost savings associated with 3′-DGE approach are likely a reasonable tradeoff for modest reduction in sensitivity and inability to observe alternative splicing, and should enable many larger scale studies focusing on not only differential expression analysis, but also quantitative transcriptome profiling.

Overexpression of CtCHS1 Increases Accumulation of Quinochalcone in Safflower.

Carthami flos, the dried petal of safflower (Carthamus tinctorius L.) has been widely used in traditional Chinese medicine to treat cardiovascular and cerebrovascular diseases, in which quinochalcone glucosides such as hydrosafflower yellow A (HSYA), carthamin are uniquely present and have been identified as active compounds. In the present study, through sequencing of a safflower floret cDNA library and subsequent microarray analysis, we found 23 unigenes (5 PALs, 1 C4Hs, 5 4CLs, 6 CHSs, 2 CHIs, 2 DFRs, 2 FLSs) involved in flavonoid pathway, of which 4 were up-regulated differentially during quinochalcone glucosides accumulation with the floret developing stage. The up-regulated genes were verified by PCR methods. Considering chalcone synthase are entry enzyme in flavonoid biosynthesis, CHS1 was focused on target gene to verify its function furtherly. Bioinformation analysis showed that CHS1 shared 86.94% conserved residues with CHS from other plants. Subcellular localization showed that CtCHS1 was localized in cytoplasm in onion epidermal cells. The transgenic safflower plant with overexpression CtCHS1 by Agrobacterium-mediated pollen-tube pathway method was firstly generated. The results present that expression of PAL2, PAL3, CHS1, CHS4, CHS6 increased and expression of CHI1 and CHI2 decreased in the transgenic plant floret. Meanwhile, the accumulation of quinochalcone glucosides increased by ∼20–30% and accumulation of quercetin-3-β-D-glucoside and quercetin decreased by 48 and 63% in the transgenic plant floret. These results suggested that CtCHS1 played an important role in quinochalcone glucosides biosynthesis rather than flavonol biosynthesis. These results also demonstrated that the pollen-tube pathway method was an efficient method for gene transformation in safflower. Our study will provide a deep understanding of potential synthetic genes involved in quinochalcone biosynthetic pathway.

Identification of genes that potentially affect social dominance hierarchy in adult male giant freshwater prawns (Macrobrachium rosenbergii)

The giant freshwater prawn (GFP) or Macrobrachium rosenbergii is one of the most important freshwater aquaculture crustaceans worldwide. An enduring problem with the improvement of this species in culture is the social dominance hierarchy that is manifested in male morphotypes where a few fast growing males (blue claw) individuals inhibit the growth of subdominant males. Here, we applied a transcriptomics approach to identify genes that contribute to expression of dominant phenotypes. Six different GFP tissues (abdominal muscle, hepatopancreas, stomach, claw, eyestalk and testis/ovary) were dissected for RNA extraction to prepare cDNA libraries for sequencing in Illumina NextSeq 500 platform. Illumina sequence reads were assembled and annotated for identifying candidate genes. This dataset currently provides the most comprehensive transcriptome data set available for GFP (147,507,623 reads assembled into 30,522 annotated contigs). 40 candidate genes were identified that could potentially play a role in determining adult male morphotype. Candidate genes identified include; Red pigment concentrating hormone (RPCH), Adipokinetic hormone receptor, Pigment dispersing hormone precursor (PDH), Pheromone and odorant, Visual pigment like receptor, Opsin, Glycoprotein 93 and Lola protein isoform a. The genes identified here provide a resource for understanding adult male morphotype development and the hierarchical dominance system in adult male GFP. In the future this information can be applied potentially to manipulate male morphotype frequencies in culture to bias adult male phenotypic ratios to improve overall productivity in practical genetic improvement programs.

The Accumulation of miRNAs Differentially Modulated by Drought Stress Is Affected by Grafting in Grapevine.

Grapevine (Vitis vinifera) is routinely grafted, and rootstocks inducing drought tolerance represent a source for adapting vineyards to climate change in temperate areas. Our goal was to investigate drought stress effects on microRNA (miRNA) abundance in a drought-resistant grapevine rootstock, M4 (Vitis vinifera × Vitis berlandieri), compared with a commercial cultivar, Cabernet Sauvignon, using their autografts and reciprocal grafts. RNA extracted from roots and leaves of droughted and irrigated plants of different graft combinations was used to prepare cDNA libraries for small RNA sequencing and to analyze miRNAs by quantitative real-time polymerase chain reaction (RT-qPCR). Measurements of leaf water potential, leaf gas exchange, and root hydraulic conductance attested that, under irrigation, M4 reduced water loss in comparison with cultivar Cabernet Sauvignon mostly through nonhydraulic, root-specific mechanisms. Under drought, stomatal conductance decreased at similar levels in the two genotypes. Small RNA sequencing allowed the identification of 70 conserved miRNAs and the prediction of 28 novel miRNAs. Different accumulation trends of miRNAs, observed upon drought and in different genotypes and organs, were confirmed by RT-qPCR Corresponding target transcripts, predicted in silico and validated by RT-qPCR, often showed opposite expression profiles than the related miRNAs. Drought effects on miRNA abundance differed between the two genotypes. Furthermore, the concentration of drought-responsive miRNAs in each genotype was affected by reciprocal grafting, suggesting either the movement of signals inducing miRNA expression in the graft partner or, possibly, miRNA transport between scion and rootstock. These results open new perspectives in the selection of rootstocks for improving grapevine adaptation to drought.

De novo transcriptome of the muga silkworm, Antheraea assamensis (Helfer)

Antheraea assamensis (Lepidoptera: Saturniidae), is a semi-domesticated silkworm known to be endemic to Assam and the adjoining hilly areas of Northeast India. It is the only producer of a unique, commercially important variety of golden silk called “muga silk”. Herein, we report the de novo transcriptome of A. assamensis reared on Machilus bombycina leaves for the first time. Short reads generated by high throughput sequencing of cDNA libraries from multiple tissues, viz. alimentary canal, silk gland and residual body of the 5th instar of muga silkworm were assembled into transcripts via a de novo assembly pipeline followed by functional annotation and classification. A total of 1,21,433 transcripts were generated from ~231 million raw reads of which ~74% (89,583) were either allocated a functional annotation or categorized under Pfam/COG/KEGG categories. Identification of differentially expressed transcripts and their comparative sequence analysis revealed candidate genes related to silk synthesis, viz. silk gland factor-1 and 3, sericin-like transcript, etc. with conserved forkhead, homeo- and POU domains. Several candidate anti-microbial peptides which may have potential anti-bacterial, anti-fungal or anti-parasitic activity in A. assamensis were also identified. T/A and AT/TA were predicted to be the most abundant mono- and di-nucleotide simple sequence repeat markers in the transcriptome. Transcriptome validation was carried out by quantitative real-time PCR (qPCR) amplification of eight transcripts. The resources generated by this study will expand the periphery of existing genomic data on A. assamensis facilitating future in-depth studies on its unknown aspects.