Transcriptome assembly
Dotaz
Zobrazit nápovědu
The hop plant (Humulus lupulus L.) produces several valuable secondary metabolites, such as prenylflavonoid, bitter acids, and essential oils. These compounds are biosynthesized in glandular trichomes (lupulin glands) endowed with pharmacological properties and widely implicated in the beer brewing industry. The present study is an attempt to generate exhaustive information of transcriptome dynamics and gene regulatory mechanisms involved in biosynthesis and regulation of these compounds, developmental changes including trichome development at three development stages, namely leaf, bract, and mature lupulin glands. Using high-throughput RNA-Seq technology, a total of 61.13, 50.01, and 20.18 Mb clean reads in the leaf, bract, and lupulin gland libraries, respectively, were obtained and assembled into 43,550 unigenes. The putative functions were assigned to 30,996 transcripts (71.17%) based on basic local alignment search tool similarity searches against public sequence databases, including GO, KEGG, NR, and COG families, which indicated that genes are principally involved in fundamental cellular and molecular functions, and biosynthesis of secondary metabolites. The expression levels of all unigenes were analyzed in leaf, bract, and lupulin glands tissues of hop. The expression profile of transcript encoding enzymes of BCAA metabolism, MEP, and shikimate pathway was most up-regulated in lupulin glands compared with leaves and bracts. Similarly, the expression levels of the transcription factors and structural genes that directly encode enzymes involved in xanthohumol, bitter acids, and terpenoids biosynthesis pathway were found to be significantly enhanced in lupulin glands, suggesting that production of these metabolites increases after the leaf development. In addition, numerous genes involved in primary metabolism, lipid metabolism, photosynthesis, generation of precursor metabolites/energy, protein modification, transporter activity, and cell wall component biogenesis were differentially regulated in three developmental stages, suggesting their involvement in the dynamics of the lupulin gland development. The identification of differentially regulated trichome-related genes provided a new foundation for molecular research on trichome development and differentiation in hop. In conclusion, the reported results provide directions for future functional genomics studies for genetic engineering or molecular breeding for augmentation of secondary metabolite content in hop.
- MeSH
- flavonoidy biosyntéza chemie metabolismus MeSH
- genová ontologie MeSH
- Humulus chemie metabolismus MeSH
- listy rostlin genetika metabolismus MeSH
- propiofenony chemie metabolismus MeSH
- regulace genové exprese u rostlin MeSH
- rostlinné proteiny genetika metabolismus MeSH
- sekvenování transkriptomu MeSH
- terpeny chemie metabolismus MeSH
- transkripční faktory metabolismus MeSH
- transkriptom genetika MeSH
- trichomy genetika metabolismus ultrastruktura MeSH
- Publikační typ
- časopisecké články MeSH
The major histocompatibility complex (MHC) plays a central role in the adaptive immune response and is the most polymorphic gene family in vertebrates. Although high-throughput sequencing has increasingly been used for genotyping families of co-amplifying MHC genes, its potential to facilitate early steps in the characterisation of MHC variation in nonmodel organism has not been fully explored. In this study we evaluated the usefulness of de novo transcriptome assembly in characterisation of MHC sequence diversity. We found that although de novo transcriptome assembly of MHC I genes does not reconstruct sequences of individual alleles, it does allow the identification of conserved regions for PCR primer design. Using the newly designed primers, we characterised MHC I sequences in the bank vole. Phylogenetic analysis of the partial MHC I coding sequence (2-4 exons) of the bank vole revealed a lack of orthology to MHC I of other Cricetidae, consistent with the high gene turnover of this region. The diversity of expressed alleles was characterised using ultra-deep sequencing of the third exon that codes for the peptide-binding region of the MHC molecule. High allelic diversity was demonstrated, with 72 alleles found in 29 individuals. Interindividual variation in the number of expressed loci was found, with the number of alleles per individual ranging from 5 to 14. Strong signatures of positive selection were found for 8 amino acid sites, most of which are inferred to bind antigens in human MHC, indicating conservation of structure despite rapid sequence evolution.
- MeSH
- alely MeSH
- Arvicolinae genetika MeSH
- DNA primery MeSH
- exony MeSH
- fylogeneze MeSH
- genetická variace MeSH
- genotyp MeSH
- geny MHC třídy I * MeSH
- hlavní histokompatibilní komplex genetika MeSH
- multigenová rodina MeSH
- myši MeSH
- transkriptom * MeSH
- vysoce účinné nukleotidové sekvenování MeSH
- zvířata MeSH
- Check Tag
- myši MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
BACKGROUND: The phlebotomine sand fly Phlebotomus perniciosus (Diptera: Psychodidae, Phlebotominae) is a major Old World vector of the protozoan Leishmania infantum, the etiological agent of visceral and cutaneous leishmaniases in humans and dogs, a worldwide re-emerging diseases of great public health concern, affecting 101 countries. Despite the growing interest in the study of this sand fly species in the last years, the development of genomic resources has been limited so far. To increase the available sequence data for P. perniciosus and to start studying the molecular basis of the sexual differentiation in sand flies, we performed whole transcriptome Illumina RNA sequencing (RNA-seq) of adult males and females and de novo transcriptome assembly. RESULTS: We assembled 55,393 high quality transcripts, of which 29,292 were unique, starting from adult whole body male and female pools. 11,736 transcripts had at least one functional annotation, including full-length low abundance salivary transcripts, 981 transcripts were classified as putative long non-coding RNAs and 244 transcripts encoded for putative novel proteins specific of the Phlebotominae sub-family. Differential expression analysis identified 8590 transcripts significantly biased between sexes. Among them, some show relaxation of selective constraints when compared to their orthologs of the New World sand fly species Lutzomyia longipalpis. CONCLUSIONS: In this paper, we present a comprehensive transcriptome resource for the sand fly species P. perniciosus built from short-read RNA-seq and we provide insights into sex-specific gene expression at adult stage. Our analysis represents a first step towards the identification of sex-specific genes and pathways and a foundation for forthcoming investigations into this important vector species, including the study of the evolution of sex-biased genes and of the sexual differentiation in phlebotomine sand flies.
- MeSH
- hmyz - vektory genetika MeSH
- Leishmania infantum genetika patogenita MeSH
- leishmanióza viscerální genetika parazitologie MeSH
- lidé MeSH
- Phlebotomus genetika parazitologie MeSH
- pohlavní dimorfismus MeSH
- psi MeSH
- sekvence aminokyselin MeSH
- transkriptom genetika MeSH
- vysoce účinné nukleotidové sekvenování MeSH
- zvířata MeSH
- Check Tag
- lidé MeSH
- mužské pohlaví MeSH
- psi MeSH
- ženské pohlaví MeSH
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Hazelnut (Corylus), which has high commercial and nutritional benefits, is an important tree for producing nuts and nut oil consumed as ingredient especially in chocolate. While Corylus avellana L. (Euro-pean hazelnut, Betulaceae) and Corylus colurna L. (Turkish hazelnut, Betulaceae) are the two common hazelnut species in Europe, C. avellana L. (Tombul hazelnut) is grown as the most widespread hazelnut species in Turkey, and C. colurna L., which is the most important genetic resource for hazelnut breeding, exists naturally in Anatolia. We generated the transcriptome data of these two Corylus species and used these data for gene discovery and gene expression profiling. Total RNA from young leaves, flowers (male and female), buds, and husk shoots of C. avellana and C. colurna were used for two different libraries and were sequenced using Illumina HiSeq4000 with 100 bp paired-end reads. The transcriptome data 10.48 and 10.30 Gb of C. avellana and C. colurna, respectively, were assembled into 70,265 and 88,343 unigenes, respectively. These unigenes were functionally annotated using the TRAPID platform. We identified 25,312 and 27,051 simple sequen-ce repeats (SSRs) for C. avellana and C. colurna, respectively. TL1, GMPM1, N, 2MMP, At1g29670, CHIB1 unigenes were selected for validation with qPCR. The first de novo transcriptome data of C. co-lurna were used to compare data of C. avellana of commercial importance. These data constitute a valuable extension of the publicly available transcriptomic resource aimed at breeding, medicinal, and industrial research studies.
- MeSH
- líska * genetika metabolismus MeSH
- ořechy MeSH
- stanovení celkové genové exprese MeSH
- Publikační typ
- časopisecké články MeSH
- Geografické názvy
- Turecko MeSH
Ectropis oblique Prout (Lepidoptera: Geometridae) is one of the main pests that damages the tea crop in Southeast Asia. To understand the molecular mechanisms of its feeding biology, transcriptomes of the alimentary tract (AT) and of the body minus the AT of E. oblique were successfully sequenced and analyzed in this study. A total of 36,950 unigenes from de novo sequences were assembled. After analysis using six annotation databases (e.g., Gene Ontology, Kyoto Encyclopedia of Genes and Genome, and NCBI nr), a series of putative genes were found for this insect species that were related to digestion, detoxification, the immune system, and Bacillus thuringiensis (Bt) receptors. From this series of genes, 21 were randomly selected to verify the relative expression levels of transcripts using quantitative real-time polymerase chain reaction. These results will provide an invaluable genomic resource for future studies on the molecular mechanisms of E. oblique, which will be useful in developing biological control strategies for this pest.
BACKGROUND: RNA-sequencing analysis is increasingly utilized to study gene expression in non-model organisms without sequenced genomes. Aethionema arabicum (Brassicaceae) exhibits seed dimorphism as a bet-hedging strategy - producing both a less dormant mucilaginous (M+) seed morph and a more dormant non-mucilaginous (NM) seed morph. Here, we compared de novo and reference-genome based transcriptome assemblies to investigate Ae. arabicum seed dimorphism and to evaluate the reference-free versus -dependent approach for identifying differentially expressed genes (DEGs). RESULTS: A de novo transcriptome assembly was generated using sequences from M+ and NM Ae. arabicum dry seed morphs. The transcripts of the de novo assembly contained 63.1% complete Benchmarking Universal Single-Copy Orthologs (BUSCO) compared to 90.9% for the transcripts of the reference genome. DEG detection used the strict consensus of three methods (DESeq2, edgeR and NOISeq). Only 37% of 1533 differentially expressed de novo assembled transcripts paired with 1876 genome-derived DEGs. Gene Ontology (GO) terms distinguished the seed morphs: the terms translation and nucleosome assembly were overrepresented in DEGs higher in abundance in M+ dry seeds, whereas terms related to mRNA processing and transcription were overrepresented in DEGs higher in abundance in NM dry seeds. DEGs amongst these GO terms included ribosomal proteins and histones (higher in M+), RNA polymerase II subunits and related transcription and elongation factors (higher in NM). Expression of the inferred DEGs and other genes associated with seed maturation (e.g. those encoding late embryogenesis abundant proteins and transcription factors regulating seed development and maturation such as ABI3, FUS3, LEC1 and WRI1 homologs) were put in context with Arabidopsis thaliana seed maturation and indicated that M+ seeds may desiccate and mature faster than NM. The 1901 transcriptomic DEG set GO-terms had almost 90% overlap with the 2191 genome-derived DEG GO-terms. CONCLUSIONS: Whilst there was only modest overlap of DEGs identified in reference-free versus -dependent approaches, the resulting GO analysis was concordant in both approaches. The identified differences in dry seed transcriptomes suggest mechanisms underpinning previously identified contrasts between morphology and germination behaviour of M+ and NM seeds.
- MeSH
- anotace sekvence MeSH
- Brassicaceae genetika růst a vývoj MeSH
- genom rostlinný MeSH
- genová ontologie MeSH
- klíčení MeSH
- regulace genové exprese u rostlin * MeSH
- rostlinné proteiny genetika MeSH
- semena rostlinná genetika růst a vývoj MeSH
- stanovení celkové genové exprese MeSH
- transkriptom * MeSH
- vysoce účinné nukleotidové sekvenování MeSH
- Publikační typ
- časopisecké články MeSH
Mitochondrial complex I (NADH:ubiquinone oxidoreductase) must be assembled precisely from 45 protein subunits for it to function correctly. One of its mitochondrial DNA (mtDNA) encoded subunits, ND1, is incorporated during the early stages of complex I assembly. However, little is known about how mutations in ND1 affect this assembly process. We found that in human 143B cybrid cells carrying a homoplasmic MT-ND1 mutation, ND1 protein could not be translated. As a result, the early stages of complex I assembly were disrupted, with mature complex I undetectable and complex I-linked respiration severely reduced to 2.0% of control levels. Interestingly, complex IV (ferrocytochrome c:oxygen oxidoreductase) steady-state levels were also reduced to 40.3%, possibly due to its diminished stability in the absence of respiratory supercomplex formation. This was in comparison with 143B cybrid controls (that contained wild-type mtDNA on the same nuclear background), which exhibited normal complex I, complex IV, and supercomplex assembly. We conclude that the loss of ND1 stalls complex I assembly during the early stages of its biogenesis, which not only results in the loss of mature complex I but also disrupts the stability of complex IV and the respiratory supercomplex to cause mitochondrial dysfunction.-Lim, S. C., Hroudová, J., Van Bergen, N. J., Lopez Sanchez, M. I. G., Trounce, I. A., McKenzie, M. Loss of mitochondrial DNA-encoded protein ND1 results in disruption of complex I biogenesis during early stages of assembly.
BACKGROUND: Parasites employ proteases to evade host immune systems, feed and replicate and are often the target of anti-parasite strategies to disrupt these interactions. Myxozoans are obligate cnidarian parasites, alternating between invertebrate and fish hosts. Their genes are highly divergent from other metazoans, and available genomic and transcriptomic datasets are limited. Some myxozoans are important aquaculture pathogens such as Sphaerospora molnari replicating in the blood of farmed carp before reaching the gills for sporogenesis and transmission. Proliferative stages cause a massive systemic lymphocyte response and the disruption of the gill epithelia by spore-forming stages leads to respiratory problems and mortalities. In the absence of a S. molnari genome, we utilized a de novo approach to assemble the first transcriptome of proliferative myxozoan stages to identify S. molnari proteases that are upregulated during the first stages of infection when the parasite multiplies massively, rather than in late spore-forming plasmodia. Furthermore, a subset of orthologs was used to characterize 3D structures and putative druggable targets. RESULTS: An assembled and host filtered transcriptome containing 9436 proteins, mapping to 29,560 contigs was mined for protease virulence factors and revealed that cysteine proteases were most common (38%), at a higher percentage than other myxozoans or cnidarians (25-30%). Two cathepsin Ls that were found upregulated in spore-forming stages with a presenilin like aspartic protease and a dipeptidyl peptidase. We also identified downregulated proteases in the spore-forming development when compared with proliferative stages including an astacin metallopeptidase and lipases (qPCR). In total, 235 transcripts were identified as putative proteases using a MEROPS database. In silico analysis of highly transcribed cathepsins revealed potential drug targets within this data set that should be prioritised for development. CONCLUSIONS: In silico surveys for proteins are essential in drug discovery and understanding host-parasite interactions in non-model systems. The present study of S. molnari's protease arsenal reveals previously unknown proteases potentially used for host exploitation and immune evasion. The pioneering dataset serves as a model for myxozoan virulence research, which is of particular importance as myxozoan diseases have recently been shown to emerge and expand geographically, due to climate change.
- MeSH
- antiparazitární látky farmakologie terapeutické užití MeSH
- faktory virulence MeSH
- fylogeneze MeSH
- kapři mikrobiologie MeSH
- Myxozoa genetika růst a vývoj MeSH
- nemoci ryb parazitologie terapie MeSH
- objevování léků MeSH
- parazitární nemoci u zvířat parazitologie terapie MeSH
- proteasy genetika MeSH
- transkriptom MeSH
- zvířata MeSH
- Check Tag
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.
- MeSH
- anotace sekvence metody MeSH
- data mining metody MeSH
- fylogeneze MeSH
- genom mitochondriální * MeSH
- Isoptera genetika MeSH
- reprodukovatelnost výsledků MeSH
- sekvence nukleotidů genetika MeSH
- sekvenční analýza DNA MeSH
- sekvenování transkriptomu MeSH
- transkriptom genetika MeSH
- vysoce účinné nukleotidové sekvenování MeSH
- zvířata MeSH
- Check Tag
- zvířata MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- validační studie MeSH
We report a series of microarray-based leaf and crown transcriptome comparisons involving three barley cultivars (cvs. Luxor, Igri and Atlas 68) which express differing degrees of frost tolerance. The transcripts were obtained following the exposure of seedlings to low (above and below zero) temperatures, aiming to identify those genes and signalling/metabolic pathways which are associated with frost tolerance. Both the leaves and the crowns responded to low temperature by the up-regulation of a suite of abscisic acid (ABA)-responsive genes, most of which have already been recognized as components of the plant low temperature response. The inter-cultivar comparison indicated that genes involved in maintaining the leaf's capacity to synthesize protein and to retain chloroplast activity were important for the expression of frost tolerance. In the crown, the repression of genes associated with nucleosome assembly and transposon regulation were the most relevant transcriptional changes associated with frost tolerance, highlighting the role of gene repression in the cold acclimation response.
- MeSH
- down regulace MeSH
- ječmen (rod) genetika metabolismus MeSH
- listy rostlin genetika metabolismus MeSH
- mapování chromozomů MeSH
- nukleozomy genetika metabolismus MeSH
- reakce na chladový šok MeSH
- regulace genové exprese u rostlin MeSH
- rostlinné geny MeSH
- rostlinné proteiny genetika metabolismus MeSH
- transkripční elongační faktory genetika metabolismus MeSH
- transkriptom MeSH
- upregulace MeSH
- zmrazování MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH