Full-length transcriptome analysis of shade-induced promotion of tuber production in Pinellia ternata

Abstract Background Pinellia ternata is native to China and has been used as a traditional herb due to its antiemetic, antitussive, analgesic, and anxiolytic effects. When exposed to strong light intensity and high temperature during the reproductive growth process, P. ternata withers in a phenomenon known as “sprout tumble”, which largely limits tuber production. Shade was previously found to delay sprout tumble formation (STF); however, no information exists regarding this process at the molecular level. Hence, we determined the genes involved in tuber development and STF in P. ternata. Results Compared to that with natural sun-light (control), shade significantly induced chlorophyll accumulation, increased chlorophyll fluorescence parameters including initial fluorescence, maximal fluorescence, and qP, and dramatically repressed chlorophyll a:b and NPQ. Catalase (CAT) activity was largely induced by shade, and tuber products were largely increased in this environment. Transcriptome profiles of P. ternata grown in natural sun-light and shaded environments were analyzed by a combination of next generation sequencing (NGS) and third generation single-molecule real-time (SMRT) sequencing. Corrections of SMRT long reads based on NGS short reads yielded 136,163 non-redundant transcripts, with an average N50 length of 2578 bp. In total, 6738 deferentially-expressed genes (DEGs) were obtained from the comparisons, specifically D5S vs D5CK, D20S vs D20CK, D20S vs D5S, and D20CK vs D5CK, of which, 6384 DEGs (94.8%) were generated from the D20S vs D20CK comparison. Gene annotation and functional analyses revealed that these genes were related to auxin signal transduction, polysaccharide and sugar metabolism, phenylpropanoid biosynthesis, and photosynthesis. Moreover, the expression of genes enriched in photosynthesis appeared to be significantly altered by shade. The expression patterns of 16 candidate genes were consistent with changes in their transcript abundance as identified by RNA-Seq, and these might contribute to STF and tuber production. Conclusion The full-length transcripts identified in this study have provided a more accurate depiction of P. ternata gene transcription. Further, we identified potential genes involved in STF and tuber growth. Such data could serve as a genetic resource and a foundation for further research on this important traditional herb.

Download Full-text

PacBio Single-Molecule Long-Read Sequencing Reveals Genes Tolerating Manganese Stress in Schima superba Saplings

Frontiers in Genetics ◽

10.3389/fgene.2021.635043 ◽

2021 ◽

Vol 12 ◽

Author(s):

Fiza Liaquat ◽

Muhammad Farooq Hussain Munis ◽

Samiah Arif ◽

Urooj Haroon ◽

Jianxin Shi ◽

...

Keyword(s):

Single Molecule ◽

Gene Annotation ◽

Treated Group ◽

Full Length ◽

Open Reading Frames ◽

Interacting Protein ◽

Schima Superba ◽

Long Read ◽

First Time ◽

Potential Tool

Schima superba (Theaceae) is a subtropical evergreen tree and is used widely for forest firebreaks and gardening. It is a plant that tolerates salt and typically accumulates elevated amounts of manganese in the leaves. With large ecological amplitude, this tree species grows quickly. Due to its substantial biomass, it has a great potential for soil remediation. To evaluate the thorough framework of the mRNA, we employed PacBio sequencing technology for the first time to generate S. Superba transcriptome. In this analysis, overall, 511,759 full length non-chimeric reads were acquired, and 163,834 high-quality full-length reads were obtained. Overall, 93,362 open reading frames were obtained, of which 78,255 were complete. In gene annotation analyses, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Genes (COG), Gene Ontology (GO), and Non-Redundant (Nr) databases were allocated 91,082, 71,839, 38,914, and 38,376 transcripts, respectively. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT), and Coding Non-Coding Index (CNCI) databases and observed 8,551, 9,174, 20,720, and 18,669 lncRNAs, respectively. Moreover, nine genes were randomly selected for the expression analysis, which showed the highest expression of Gene 6 (Na_Ca_ex gene), and CAX (CAX-interacting protein 4) was higher in manganese (Mn)-treated group. This work provided significant number of full-length transcripts and refined the annotation of the reference genome, which will ease advanced genetic analyses of S. superba.

Download Full-text

Transcriptome analysis and molecular mechanism of linseed (Linum usitatissimum L.) drought tolerance under repeated drought using single-molecule long-read sequencing

BMC Genomics ◽

10.1186/s12864-021-07416-5 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Wei Wang ◽

Lei Wang ◽

Ling Wang ◽

Meilian Tan ◽

Collins O. Ogutu ◽

...

Keyword(s):

Dna Repair ◽

Drought Tolerance ◽

Single Molecule ◽

Linum Usitatissimum ◽

Expression Patterns ◽

Transcriptional Response ◽

Enrichment Analysis ◽

Resistant Variety ◽

Proline Biosynthesis ◽

Interacting Protein

Abstract Background Oil flax (linseed, Linum usitatissimum L.) is one of the most important oil crops., However, the increases in drought resulting from climate change have dramatically reduces linseed yield and quality, but very little is known about how linseed coordinates the expression of drought resistance gene in response to different level of drought stress (DS) on the genome-wide level. Results To explore the linseed transcriptional response of DS and repeated drought (RD) stress, we determined the drought tolerance of different linseed varieties. Then we performed full-length transcriptome sequencing of drought-resistant variety (Z141) and drought-sensitive variety (NY-17) under DS and RD stress at the seedling stage using single-molecule real-time sequencing and RNA-sequencing. Gene Ontology (GO) and reduce and visualize GO (REVIGO) enrichment analysis showed that upregulated genes of Z141 were enriched in more functional pathways related to plant drought tolerance than those of NY-17 were under DS. In addition, 4436 linseed transcription factors were identified, and 1190 were responsive to stress treatments. Moreover, protein-protein interaction (PPI) network analysis showed that the proline biosynthesis pathway interacts with stress response genes through RAD50 (DNA repair protein 50) interacting protein 1 (RIN-1). Finally, proline biosynthesis and DNA repair structural gene expression patterns were verified by RT- PCR. Conclusions The drought tolerance of Z141 may be related to its upregulation of drought tolerance genes under DS. Proline may play an important role in linseed drought tolerance by maintaining cell osmotic and protecting DNA from ROS damage. In summary, this study provides a new perspective to understand the drought adaptability of linseed.

Download Full-text

Full-length-transcriptomic analysis in mice and human heart using Single-Molecule Real-time Sequencing (SMRT) identified 15 novel isoforms and a novel promoter region of PGC1-alpha

European Heart Journal ◽

10.1093/ehjci/ehaa946.3582 ◽

2020 ◽

Vol 41 (Supplement_2) ◽

Author(s):

D Oehler ◽

A Goedecke ◽

A Spychala ◽

K Lu ◽

N Gerdes ◽

...

Keyword(s):

Alternative Splicing ◽

Single Molecule ◽

High Fat Diet ◽

Human Heart ◽

Full Length ◽

Funding Source ◽

Smrt Sequencing ◽

High Fat ◽

Exon 1 ◽

Novel Exon

Abstract Background Alternative splicing is a process by which exons within a pre-mRNA are joined or skipped, resulting in isoforms being encoded by a single gene. Alternative Splicing affecting transcription factors may have substantial impact on cellular dynamics. The PPARG Coactivator 1 Alpha (PGC1-α), is a major modulator in energy metabolism. Data from murine skeletal muscle revealed distinctive isoform patterns giving rise to different phenotypes, i.e. mitogenesis and hypertrophy. Here, we aimed to establish a complete dataset of isoforms in murine and human heart applying single-molecule real-time (SMRT)-sequencing as novel approach to identify transcripts without need for assembly, resulting in true full-length sequences. Moreover, we aimed to unravel functional relevance of the various isoforms during experimental ischemia reperfusion (I/R). Methods RNA-Isolation was performed in murine (C57Bl/6J) or human heart tissue (obtained during LVAD-surgery), followed by library preparation and SMRT-Sequencing. Bioinformatic analysis was done using a modified IsoSeq3-Pipeline and OS-tools. Identification of PGC1-α isoforms was fulfilled by similarity search against exonic sequences within the full-length, non-concatemere (FLNC) reads. Isoforms with Open-Reading-Frame (ORF) were manually curated and validated by PCR and Sanger-Sequencing. I/R was induced by ligature of the LAD for 45 min in mice on standard chow as well as on high-fat-high-sucrose diet. Area At Risk (AAR) and remote tissue were collected three and 16 days after I/R or sham-surgery (n=4 per time point). Promotor patterns were analyzed by qPCR. Results Deciphering the full-length transcriptome of murine and human heart resulted in ∼60000 Isoforms with 99% accuracy on mRNA-sequence. Focusing on murine PGC1-α-isoforms we discovered and verified 15 novel transcripts generated by hitherto unknown splicing events. Additionally, we identified a novel Exon 1 originating between the known promoters followed by a valid ORF, suggesting the discovery of a novel promoter. Remarkably, we found a homologous novel Exon1 in human heart, suggesting conservation of the postulated promoter. In I/R the AAR exhibited a significant lower expression of established and novel promoters compared to remote under standard chow 3d post I/R. 16d post I/R, the difference between AAR & Remote equalized in standard chow while remaining under High-Fat-Diet. Conclusion Applying SMRT-technique, we generated the first time a complete full-length-transcriptome of the murine and human heart, identifying 15 novel potentially coding transcripts of PGC1-α and a novel exon 1. These transcripts are differentially regulated in experimental I/R in AAR and remote myocardium, suggesting transcriptional regulation and alternative splicing modulating PGC1-α function in heart. Differences between standard chow and high fat diet suggest impact of impaired glucose metabolism on regulatory processes after myocardial infarction. Funding Acknowledgement Type of funding source: Public grant(s) – National budget only. Main funding source(s): Collaborative Research Centre 1116 (German Research Foundation)

Download Full-text

Spatially resolved transcriptomics reveals unique gene signatures associated with human temporal cortical architecture and Alzheimer's pathology

10.1101/2021.07.07.451554 ◽

2021 ◽

Author(s):

Shuo Chen ◽

Yuzhou Chang ◽

Liangping Li ◽

Diana Acosta ◽

Cody Morrison ◽

...

Keyword(s):

Gene Expression ◽

Single Molecule ◽

Brain Region ◽

Expression Patterns ◽

Middle Temporal Gyrus ◽

Marker Genes ◽

Gene Signatures ◽

Unique Gene ◽

Spatially Resolved ◽

And Gender

Alzheimer's disease (AD) is pathologically characterized by amyloid beta (Aβ) plaques, neurofibrillary tangles (tau aggregates), and alterations in microglia, astrocytes and oligodendrocytes. The mesial temporal lobe is a vulnerable brain region in early AD; however, little is known about the transcriptome-scale gene expression in this region and its relation to AD pathology. Here we use the 10x Genomics Visium platform in combination with co-immunofluorescence staining of AD-associated pathological markers to define the spatial topography of gene expression in the middle temporal gyrus (MTG) from both early AD and age- and gender-matched control cases. We identify unique marker genes for six cortical layers and the adjacent white matter as well as gene expression patterns and alterations that showcase unique gene signatures and pathways associated with a range of AD pathology. Also, gene co-expression analyses of differentially expressed genes (DEGs) between AD and controls reveal four unique gene modules, which significantly change their co-expression patterns in the presence of variations of AD pathology. Furthermore, we validate the changes of key representative DEGs that are associated with AD pathology in neurons, microglia, astrocytes and oligodendrocytes using single-molecule fluorescent in situ hybridization. In summary, we provide a rich resource for the spatial transcriptomic profile of the human MTG, which will contribute to our understanding of the complex architecture and AD pathology of this vulnerable brain region.

Download Full-text

Effects of Pathological Mutations on the Prion-Like Polymerisation of MyD88

10.1101/351726 ◽

2018 ◽

Author(s):

Ailís O’Carroll ◽

Brieuc Chauvin ◽

James Brown ◽

Ava Meagher ◽

Joanne Coyle ◽

...

Keyword(s):

Single Molecule ◽

Self Assembly ◽

Point Mutations ◽

Adaptor Protein ◽

Full Length ◽

Binary Response ◽

Death Domain ◽

Tir Domain ◽

The Impact

AbstractA novel concept has emerged whereby the higher-order self-assembly of proteins provides a simple and robust mechanism for signal amplification. This appears to be a universal signalling mechanism within the innate immune system, where the recognition of pathogens or danger-associated molecular patterns need to trigger a strong, binary response within cells. Previously, multiple structural studies have been limited to single domains, expressed and assembled at high protein concentrations. We therefore set out to develop new in vitro strategies to characterise the behaviour of full-length proteins at physiological levels. In this study we focus on the adaptor protein MyD88, which contains two domains with different self-assembly properties: a TIR domain that can polymerise similarly to the TIR domain of Mal, and a Death Domain that has been shown to oligomerise with helical symmetry in the Myddosome complex. To visualize the behaviour of full-length MyD88 without purification steps, we use single-molecule fluorescence coupled to eukaryotic cell-free protein expression. These experiments demonstrate that at low protein concentration, only full-length MyD88 forms prion-like polymers. We also demonstrate that the metastability of MyD88 polymerisation creates the perfect binary response required in innate signalling: the system is silenced at normal concentrations but upstream signalling creates a “seed” that triggers polymerisation and amplification of the response. These findings pushed us to re-interpret the role of polymerisation in MyD88-related diseases and we studied the impact of disease-associated point mutations L93P, R196C and L252P/L265P at the molecular level. We discovered that all mutations completely block the ability of MyD88 to polymerise. We also confirm that L252P, a gain-of-function mutation, allows the MyD88 mutant to form extremely stable oligomers, even when expressed at low nanomolar concentrations. Thus, our results are consistent with and greatly add to the findings on the Myddosomes digital ‘all-or-none’ responses and the behaviour of the oncogenic mutation of MyD88.

Download Full-text

Single molecule real time (SMRT) full length RNA-sequencing reveals novel and distinct mRNA isoforms in human bone marrow cell subpopulations

10.1101/555326 ◽

2019 ◽

Cited By ~ 1

Author(s):

Anne Deslattes Mays ◽

Marcel O. Schmidt ◽

Garrett T. Graham ◽

Elizabeth Tseng ◽

Primo Baybayan ◽

...

Keyword(s):

Bone Marrow ◽

Real Time ◽

Rna Sequencing ◽

Single Molecule ◽

Marrow Cell ◽

Full Length ◽

Protein Isoforms ◽

Mrna Isoforms ◽

Transcript Isoforms ◽

Cell Subpopulations

AbstractHematopoietic cells are continuously replenished from progenitor cells that reside in the bone marrow. To evaluate molecular changes during this process, we analyzed the transcriptomes of freshly harvested human bone marrow progenitor (lineage-negative) and differentiated (lineage-positive) cells by single molecule, real time (SMRT) full length RNA sequencing. This analysis revealed a ∼5-fold higher number of transcript isoforms than previously detected and showed a distinct composition of individual transcript isoforms characteristic for bone marrow subpopulations. A detailed analysis of mRNA isoforms transcribed from the ANXA1 and EEF1A1 loci confirmed their distinct composition. The expression of proteins predicted from the transcriptome analysis was validated by mass spectrometry and validated previously unknown protein isoforms predicted e.g. for EEF1A1. These protein isoforms distinguished the lineage negative cell population from the lineage positive cell population. Finally, transcript isoforms expressed from paralogous gene loci (e.g. CFD, GATA2, HLA-A, B & C) also distinguished cell subpopulations but were only detectable by full length RNA sequencing. Thus, qualitatively distinct transcript isoforms from individual genomic loci separate bone marrow cell subpopulations indicating complex transcriptional regulation and protein isoform generation during hematopoiesis.

Download Full-text

Two classes of active transcription sites and their roles in developmental regulation

10.1101/2020.06.23.167338 ◽

2020 ◽

Author(s):

Sarah Robinson-Thiewes ◽

John McCloskey ◽

Judith Kimble

Keyword(s):

Single Molecule ◽

Developmental Regulation ◽

Full Length ◽

Genes Encoding ◽

Transcription Sites ◽

Core Components ◽

Active Transcription ◽

Multiple Levels ◽

Nascent Transcripts ◽

Transcriptional Output

AbstractGenes encoding powerful developmental regulators are exquisitely controlled, often at multiple levels. Here, we use single molecule FISH (smFISH) to investigate nuclear active transcription sites (ATS) and cytoplasmic mRNAs of three key regulatory genes along the C. elegans germline developmental axis. The genes encode ERK/MAP kinase and core components of the Notch-dependent transcription complex. Using differentially-labeled probes spanning either a long first intron or downstream exons, we identify two ATS classes that differ in transcriptional progression: iATS harbor partial nascent transcripts while cATS harbor full-length nascent transcripts. Remarkably, the frequencies of iATS and cATS are patterned along the germline axis in a gene-, stage- and sex-specific manner. Moreover, regions with more frequent iATS make fewer full-length nascent transcripts and mRNAs, whereas those with more frequent cATS produce more of them. We propose that the regulated balance of these two ATS classes has a major impact on transcriptional output during development.

Download Full-text

Allosteric activation of preformed EGF receptor dimers by binding of a single ligand

10.21203/rs.3.rs-182557/v1 ◽

2021 ◽

Author(s):

Endang Purba ◽

Ei-ichiro Saita ◽

Reetesh Akhouri ◽

Lars-Göran Öfverstedt ◽

Gunnar Wilken ◽

...

Keyword(s):

Ligand Binding ◽

Single Molecule ◽

Molecular Mechanisms ◽

Egf Receptor ◽

Electron Tomography ◽

Growth Factor Receptor ◽

Full Length ◽

Cancer Therapies ◽

Receptor Dimer

Abstract Aberrant activation of the epidermal growth factor receptor (EGFR) by mutations has been implicated in a variety of human cancers. Elucidation of the structure of the full-length receptor is essential to understand the molecular mechanisms underlying its activation. Unlike previously anticipated, here, we report that purified full-length EGFR adopts a homodimeric form in vitro before and after activation. Cryo-electron tomography analysis of the purified receptor also showed that the extracellular domains of the receptor dimer, which are conformationally flexible before activation, are stabilised by ligand binding. Consistently, optical single-molecule observation also demonstrated that binding of only one ligand activates the receptor dimer on the cell surface. Based on these results, we propose an allosteric model for the activation of EGFR dimers by ligand binding. Our results demonstrate how oncogenic mutations spontaneously activate the receptor and shed light on the development of novel cancer therapies.

Download Full-text

Metagenomics reveals niche partitioning within the phototrophic zone of a microbial mat

10.1101/151704 ◽

2017 ◽

Author(s):

Jackson Z Lee ◽

R Craig Everroad ◽

Ulas Karaoz ◽

Angela M Detweiler ◽

Jennifer Pett-Ridge ◽

...

Keyword(s):

Microbial Mats ◽

Environmental Gradients ◽

Gene Annotation ◽

Niche Partitioning ◽

Metabolic Diversity ◽

Sequence Coverage ◽

Carbohydrate Utilization ◽

Strong Light ◽

Elkhorn Slough ◽

Conserved Gene

AbstractHypersaline photosynthetic microbial mats are stratified microbial communities known for their taxonomic and metabolic diversity and strong light-driven day-night environmental gradients. In this study of the upper photosynthetic zone of hypersaline microbial mats of Elkhorn Slough, California (USA), we show how reference-based and reference-free methods can be used to meaningfully assess microbial ecology and genetic partitioning in these complex microbial systems. Mapping of metagenome reads to the dominantCyanobacteriaobserved in the system,Coleofasciculus (Microcoleus) chthonoplastes, was used to examine strain variants within these metagenomes. Highly conserved gene subsystems indicate a core genome for the species, and a number of variant genes and subsystems suggest strain level differentiation, especially for carbohydrate utilization. Metagenome sequence coverage binning was used to assess ecosystem partitioning of remaining microbes. Functional gene annotation of these bins (primarily ofProteobacteria, Bacteroidetes,andCyanobacteria) recapitulated the known biogeochemical functions in microbial mats using a genetic basis, and also revealed evidence of novel functional diversity within theGemmatimonadetesandGammaproteobacteria. Combined, these two approaches show how genetic partitioning can inform biogeochemical partitioning of the metabolic diversity within microbial ecosystems.

Download Full-text

A method for high precision sequencing of near full-length 16S rRNA genes on an Illumina MiSeq

PeerJ ◽

10.7717/peerj.2492 ◽

2016 ◽

Vol 4 ◽

pp. e2492 ◽

Cited By ~ 29

Author(s):

Catherine M. Burke ◽

Aaron E. Darling

Keyword(s):

16S Rrna ◽

16S Rrna Gene ◽

High Throughput ◽

Single Molecule ◽

Illumina Miseq ◽

Full Length ◽

16S Rrna Genes ◽

Rrna Genes ◽

Rrna Gene ◽

Bacterial Taxonomy

BackgroundThe bacterial 16S rRNA gene has historically been used in defining bacterial taxonomy and phylogeny. However, there are currently no high-throughput methods to sequence full-length 16S rRNA genes present in a sample with precision.ResultsWe describe a method for sequencing near full-length 16S rRNA gene amplicons using the high throughput Illumina MiSeq platform and test it using DNA from human skin swab samples. Proof of principle of the approach is demonstrated, with the generation of 1,604 sequences greater than 1,300 nt from a single Nano MiSeq run, with accuracy estimated to be 100-fold higher than standard Illumina reads. The reads were chimera filtered using information from a single molecule dual tagging scheme that boosts the signal available for chimera detection.ConclusionsThis method could be scaled up to generate many thousands of sequences per MiSeq run and could be applied to other sequencing platforms. This has great potential for populating databases with high quality, near full-length 16S rRNA gene sequences from under-represented taxa and environments and facilitates analyses of microbial communities at higher resolution.

Download Full-text