Reviving the Transcriptome Studies: An Insight Into the Emergence of Single-Molecule Transcriptome Sequencing

The cacao gene atlas: a transcriptome developmental atlas reveals highly tissue-specific and dynamically-regulated gene networks in Theobroma cacao L

BACKGROUND

Theobroma cacao, the cocoa tree, is a tropical crop grown for its highly valuable cocoa solids and fat which are the basis of a 200-billion-dollar annual chocolate industry. However, the long generation time and difficulties associated with breeding a tropical tree crop have limited the progress of breeders to develop high-yielding disease-resistant varieties. Development of marker-assisted breeding methods for cacao requires discovery of genomic regions and specific alleles of genes encoding important traits of interest. To accelerate gene discovery, we developed a gene atlas composed of a large dataset of replicated transcriptomes with the long-term goal of progressing breeding towards developing high-yielding elite varieties of cacao.

RESULTS

We describe the creation of the Cacao Transcriptome Atlas, its global characterization and define sets of genes co-regulated in highly organ- and temporally-specific manners. RNAs were extracted and transcriptomes sequenced from 123 different tissues and stages of development representing major organs and developmental stages of the cacao lifecycle. In addition, several experimental treatments and time courses were performed to measure gene expression in tissues responding to biotic and abiotic stressors. Samples were collected in replicates (3-5) to enable statistical analysis of gene expression levels for a total of 390 transcriptomes. To promote wide use of these data, all raw sequencing data, expression read mapping matrices, scripts, and other information used to create the resource are freely available online. We verified our atlas by analyzing the expression of genes with known functions and expression patterns in Arabidopsis (ACT7, LEA19, AGL16, TIP13, LHY, MYB2) and found their expression profiles to be generally similar between both species. We also successfully identified tissue-specific genes at two thresholds in many tissue types represented and a set of genes highly conserved across all tissues.

CONCLUSION

The Cacao Gene Atlas consists of a gene expression browser with graphical user interface and open access to raw sequencing data files as well as the unnormalized and CPM normalized read count data mapped to several cacao genomes. The gene atlas is a publicly available resource to allow rapid mining of cacao gene expression profiles. We hope this resource will be used to help accelerate the discovery of important genes for key cacao traits such as disease resistance and contribute to the breeding of elite varieties to help farmers increase yields.

A multi-tissue de novo transcriptome assembly and relative gene expression of the vulnerable freshwater salmonid Thymallus ligericus

Freshwater ecosystems are among the most endangered ecosystems worldwide. While numerous taxa are on the verge of extinction as a result of global changes and direct or indirect anthropogenic activity, genomic and transcriptomic resources represent a key tool for comprehending species' adaptability and serve as the foundation for conservation initiatives. The Loire grayling, Thymallus ligericus, is a freshwater European salmonid endemic to the upper Loire River basin. The species is comprised of fragmented populations that are dispersed over a small area and it has been identified as a vulnerable species. Here, we provide a multi-tissue de novo transcriptome assembly of T. ligericus. The completeness and integrity of the transcriptome were assessed before and after redundancy removal with lineage-specific libraries from Eukaryota, Metazoa, Vertebrata, and Actinopterygii. Relative gene expression was assessed for each of the analyzed tissues, using the de novo assembled transcriptome and a genome-based analysis using the available T. thymallus genome as a reference. The final assembly, with a contig N50 of 1221 and Benchmarking Universal Single-Copy Orthologs (BUSCO) scores above 94%, is made accessible along with structural and functional annotations and relative gene expression of the five tissues (NCBI SRA and FigShare databases). This is the first transcriptomic resource for this species, which provides a foundation for future research on this and other salmonid species that are increasingly exposed to environmental stressors.

Status of single-cell RNA sequencing for reproductive toxicology in zebrafish and the transcriptomic trade-off

The utilization of transcriptomic studies identifying profiles of gene expression, especially in toxicogenomics, has catapulted next-generation sequencing to the forefront of reproductive toxicology. An innovative yet underutilized RNA sequencing technique emerging into this field is single-cell RNA sequencing (scRNA-seq), which provides sequencing at the individual cellular level of gonad tissue. ScRNA-seq provides a novel and unique perspective for identifying distinct cellular profiles, including identification of rare cell subtypes. The specificity of scRNA-seq is a powerful tool for reproductive toxicity research, especially for translational animal models including zebrafish. Studies to date not only have focused on 'tissue atlassing' or characterizing what cell types make up different tissues but have also begun to include toxicant exposure as a factor that this review aims to explore. Future scRNA-seq studies will contribute to understanding exposure-induced outcomes; however, the trade-offs with traditional methods need to be considered.

Identification of testis development-related genes by combining Iso-Seq and RNA-Seq in Zeugodacus tau

Introduction: Zeugodacus tau (Walker) is an invasive pest. An effective method to control this pest is the sterile insect technique (SIT). To better apply this technique, it is necessary to understand testis development progression. Methods: Differentially expressed genes (DEGs) during testis development were analyzed by PacBio Iso-Seq and RNA-seq. Results: RNA-Seq library of Z. tau testes on day 1, 6, and 11 post eclosion were constructed. We identified 755 and 865 differentially expressed genes in the comparisons of T6 (testes on day 6) vs. T1 and T11 vs. T1, respectively. The KEGG pathway analysis showed that the DEGs were significantly enriched in retinol metabolism, vitamin B6 metabolism, and ascorbate and aldarate metabolism pathways. Knockdown of retinol dehydrogenase 12-like (rdh12-like), pyridoxal kinase (pdxk) and regucalcin (rgn), the representative gene in each of the above 3 pathways, reduced the hatching rate of Z. tau offspring. In addition, we identified 107 Drosophila spermatogenesis-related orthologous genes in Z. tau, of which innexin 2 (inx2) exhibited significantly up-regulated expression throughout testis development, and the knockdown of this gene reduced offspring hatching rate. Discussion: Our data indicated that rdh12-like, pdxk, rgn, and inx2 genes were related to testis development, and they were conserved in tephritid species. These results suggested that this gene might have the same function in tephritid. The findings provide an insight into testis development and spermatogenesis in tephritid species.

Full-Length Transcriptome Sequencing and Comparative Transcriptomic Analyses Provide Comprehensive Insight into Molecular Mechanisms of Flavonoid Metabolites Biosynthesis in Styphnolobium japonicum

Styphnolobium japonicum L. is a commonly consumed plant in China, known for its medicinal and nutritional benefits. This study focuses on the medicinal properties influenced by flavonoid metabolites, which vary during flower development. Utilizing full-length transcriptome sequencing on S. japonicum flowers, we observed changes in gene expression levels as the flowers progressed through growth stages. During stages S1 and S2, key genes related to flavonoid synthesis (PAL, 4CL, CHS, F3H, etc.) exhibited heightened expression. A weighted gene co-expression network analysis (WGCNA) identified regulatory genes (MYB, bHLH, WRKY) potentially involved in the regulatory network with flavonoid biosynthesis-related genes. Our findings propose a regulatory mechanism for flavonoid synthesis in S. japonicum flowers, elucidating the genetic underpinnings of this process. The identified candidate genes present opportunities for genetic enhancements in S. japonicum, offering insights into potential applications for improving its medicinal attributes.

iFLAS: positive-unlabeled learning facilitates full-length transcriptome-based identification and functional exploration of alternatively spliced isoforms in maize

The advent of full-length transcriptome sequencing technologies has accelerated the discovery of novel splicing isoforms. However, existing alternative splicing (AS) tools are either tailored for short-read RNA-Seq data or designed for human and animal studies. The disparities in AS patterns between plants and animals still pose a challenge to the reliable identification and functional exploration of novel isoforms in plants. Here, we developed integrated full-length alternative splicing analysis (iFLAS), a plant-optimized AS toolkit that introduced a semi-supervised machine learning method known as positive-unlabeled (PU) learning to accurately identify novel isoforms. iFLAS also enables the investigation of AS functions from various perspectives, such as differential AS, poly(A) tail length, and allele-specific AS (ASAS) analyses. By applying iFLAS to three full-length transcriptome sequencing datasets, we systematically identified and functionally characterized maize (Zea mays) AS patterns. We found intron retention not only introduces premature termination codons, resulting in lower expression levels of isoforms, but may also regulate the length of 3'UTR and poly(A) tail, thereby affecting the functional differentiation of isoforms. Moreover, we observed distinct ASAS patterns in two genes within heterosis offspring, highlighting their potential value in breeding. These results underscore the broad applicability of iFLAS in plant full-length transcriptome-based AS research.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing

BACKGROUND

Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates.

METHODS

In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species.

RESULTS

A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

Characterization and analysis of multi-organ full-length transcriptomes in Sphaeropteris brunoniana and Alsophila latebrosa highlight secondary metabolism and chloroplast RNA editing pattern of tree ferns

BACKGROUND

Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events.

RESULTS

In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species.

CONCLUSIONS

Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Full-length transcriptome characterization of Platycladus orientalis based on the PacBio platform

As a unique and native conifer in China, Platycladus orientalis is widely used in soil erosion control, garden landscapes, timber, and traditional Chinese medicine. However, due to the lack of reference genome and transcriptome, it is limited to the further molecular mechanism research and gene function mining. To develop a full-length reference transcriptome, tissues from five different parts of P. orientalis and four cone developmental stages were sequenced and analyzed by single-molecule real-time (SMRT) sequencing through the PacBio platform in this study. Overall, 37,111 isoforms were detected by PacBio with an N50 length of 2,317 nt, an average length of 1,999 bp, and the GC content of 41.81%. Meanwhile, 36,120 coding sequences, 5,645 simple sequence repeats (SSRs), 1,201 non-coding RNAs (lncRNAs), and 182 alternative splicing (AS) events with five types were identified using the results obtained from the PacBio transcript isoforms. Furthermore, 1,659 transcription factors (TFs) were detected and belonged to 51 TF families. A total of 35,689 transcripts (96.17%) were annotated through the NCBI nr, KOG, Swiss-Prot and KEGG databases, and 385 transcript isoforms related to 8 types of hormones were identified incorporated into plant hormone signal transduction pathways. The assembly and revelation of the full-length transcriptome of P. orientalis offer a pioneering insight for future investigations into gene function and genetic breeding within Platycladus species.

PerFSeeB: designing long high-weight single spaced seeds for full sensitivity alignment with a given number of mismatches

BACKGROUND

Technical progress in computational hardware allows researchers to use new approaches for sequence alignment problems. For a given sequence, we usually use smaller subsequences (anchors) to find possible candidate positions within a reference sequence. We may create pairs ("position", "subsequence") for the reference sequence and keep all such records without compression, even on a budget computer. As sequences for new and reference genomes differ, the goal is to find anchors, so we tolerate differences and keep the number of candidate positions with the same anchors to a minimum. Spaced seeds (masks ignoring symbols at specific locations) are a way to approach the task. An ideal (full sensitivity) spaced seed should enable us to find all such positions subject to a given maximum number of mismatches permitted.

RESULTS

Several algorithms to assist seed generation are presented. The first one finds all permitted spaced seeds iteratively. We observe specific patterns for the seeds of the highest weight. There are often periodic seeds with a simple relation between block size, length of the seed and read. The second algorithm produces blocks for periodic seeds for blocks of up to 50 symbols and up to nine mismatches. The third algorithm uses those lists to find spaced seeds for reads of an arbitrary length. Finally, we apply seeds to a real dataset and compare results for other popular seeds.

CONCLUSIONS

PerFSeeB approach helps to significantly reduce the number of reads' possible alignment positions for a known number of mismatches. Lists of long, high-weight spaced seeds are available in Additional file 1. The seeds are best in weight compared to seeds from other papers and can usually be applied to shorter reads. Codes for all algorithms and periodic blocks can be found at https://github.com/vtman/PerFSeeB .

Characterizing the tumor immune microenvironment of ependymomas using targeted gene expression profiles and RNA sequencing

BACKGROUND

Defining the tumor immune microenvironment (TIME) of patients using transcriptome analysis is gaining more popularity. Here, we examined and discussed the pros and cons of using RNA sequencing for fresh frozen samples and targeted gene expression immune profiles (NanoString) for formalin-fixed, paraffin-embedded (FFPE) samples to characterize the TIME of ependymoma samples.

RESULTS

Our results showed a stable expression of the 40 housekeeping genes throughout all samples. The Pearson correlation of the endogenous genes was high. To define the TIME, we first checked the expression of the PTPRC gene, known as CD45, and found it was above the detection limit in all samples by both techniques. T cells were identified consistently using the two types of data. In addition, both techniques showed that the immune landscape was heterogeneous in the 6 ependymoma samples used for this study.

CONCLUSIONS

The low-abundant genes were detected in higher quantities using the NanoString technique, even when FFPE samples were used. RNA sequencing is better suited for biomarker discovery, fusion gene detection, and getting a broader overview of the TIME. The technique that was used to measure the samples had a considerable effect on the type of immune cells that were identified. The limited number of tumor-infiltrating immune cells compared to the high density of tumor cells in ependymoma can limit the sensitivity of RNA expression techniques regarding the identification of the infiltrating immune cells.

Full-Length Transcriptome Sequencing Reveals the Molecular Mechanism of Metasequoia glyptostroboides Seed Responding to Aging

Metasequoia glyptostroboides, Hu and W. C. Cheng, as the only surviving relict species of the Taxodiaceae Metasequoia genus, is a critically endangered and protected species in China. There is a risk of extinction due to the low vigor of M. glyptostroboides seeds, and the physiological mechanism of seed aging in M. glyptostroboides is not yet clear. In order to investigate the physiological and molecular mechanisms underlying the aging process of M. glyptostroboides seeds, we analyzed the antioxidant system and transcriptome at 0, 2, 4, 6, and 8 days after artificial accelerated aging treatment at 40 °C and 100% relative humidity. It was found that the germination percentage of fresh dried M. glyptostroboides seeds was 54 ± 5.29%, and significantly declined to 9.33 ± 1.88% after 6 days of aging, and then gradually decreased until the seed died on day 8. Superoxide dismutase (SOD) activity, ascorbic acid (AsA), glutathione (GSH) content and superoxide anion (O₂^·-) content and production rate significantly decreased, while malondialdehyde (MDA) and hydrogen peroxide (H₂O₂) content and glutathione peroxidase (GPX) and catalase (CAT) activity gradually increased during the aging process. A total of 42,189 unigenes were identified in the whole transcriptome, and 40,446 (95.86%) unigenes were annotated in at least one protein database. A total of 15,376 differentially expressed genes (DEGs) were obtained; KEGG enrichment analysis results revealed that seed aging may be mainly involved in the protein-processing pathways in endoplasmic reticulum, oxidative phosphorylation, and ascorbate and aldarate metabolism. Weighted gene co-expression network analysis (WGCNA) revealed that the dark magenta, orange, and medium purple modules were highly correlated with physiological indicators such as SOD, CAT, and GSH and further identified 40 hub genes such as Rboh, ACO, HSF, and CML as playing important roles in the antioxidant network of M. glyptostroboides seeds. These findings provide a broader perspective for studying the regulatory mechanism of seed aging and a large number of potential target genes for the breeding of other endangered gymnosperms.

Mitochondrial dysfunction and oxidative stress in Alzheimer's disease, and Parkinson's disease, Huntington's disease and Amyotrophic Lateral Sclerosis -An updated review

Misfolded proteins in the central nervous system can induce oxidative damage, which can contribute to neurodegenerative diseases in the mitochondria. Neurodegenerative patients face early mitochondrial dysfunction, impacting energy utilization. Amyloid-ß and tau problems both have an effect on mitochondria, which leads to mitochondrial malfunction and, ultimately, the onset of Alzheimer's disease. Cellular oxygen interaction yields reactive oxygen species within mitochondria, initiating oxidative damage to mitochondrial constituents. Parkinson's disease, linked to oxidative stress, α-synuclein aggregation, and inflammation, results from reduced brain mitochondria activity. Mitochondrial dynamics profoundly influence cellular apoptosis via distinct causative mechanisms. The condition known as Huntington's disease is characterized by an expansion of polyglutamine, primarily impactingthe cerebral cortex and striatum. Research has identified mitochondrial failure as an early pathogenic mechanism contributing to HD's selective neurodegeneration. The mitochondria are organelles that exhibit dynamism by undergoing fragmentation and fusion processes to attain optimal bioenergetic efficiency. They can also be transported along microtubules and regulateintracellular calcium homeostasis through their interaction with the endoplasmic reticulum. Additionally, the mitochondria produce free radicals. The functions of eukaryotic cells, particularly in neurons, have significantly deviated from the traditionally assigned role of cellular energy production. Most of them areimpaired in HD, which may lead to neuronal dysfunction before symptoms manifest. This article summarises the most important changes in mitochondrial dynamics that come from neurodegenerative diseases including Alzheimer's, Parkinson's, Huntington's and Amyotrophic Lateral Sclerosis. Finally, we discussed about novel techniques that can potentially treat mitochondrial malfunction and oxidative stress in four most dominating neuro disorders.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Neuropharmacology of Alcohol Addiction with Special Emphasis on Proteomic Approaches for Identification of Novel Therapeutic Targets

Alcohol is a generic pharmacological agent with only a few recognized primary targets. Nmethyl- D-aspartate, gamma-aminobutyric acid, glycine, 5-hydroxytryptamine 3 (serotonin), nicotinic acetylcholine receptors, and L-type Ca2+ channels and G-protein-activated inwardly rectifying K channels are all involved. Following the first hit of alcohol on specific brain targets, the second wave of indirect effects on various neurotransmitter/neuropeptide systems begins, leading to the typical acute behavioral effects of alcohol, which range from disinhibition to sedation and even hypnosis as alcohol concentrations rise. Recent research has revealed that gene regulation is significantly more complex than previously thought and does not fully explain changes in protein levels. As a result, studying the proteome directly, which differs from the genome/transcriptome in terms of complexity and dynamicity, has provided unique insights into extraordinary advances in proteomic techniques that have changed the way we can analyze the composition, regulation, and function of protein complexes and pathways underlying altered neurobiological conditions. Neuroproteomics has the potential to revolutionize alcohol research by allowing researchers to gain a better knowledge of how alcohol impacts protein structure, function, connections, and networks on a global scale. The amount of information collected from these breakthroughs can aid in identifying valuable biomarkers for early detection and improved prognosis of an alcohol use disorder and future pharmaceutical targets for the treatment of alcoholism.

Engineering Resistance against Sclerotinia sclerotiorum Using a Truncated NLR (TNx) and a Defense-Priming Gene

The association of both cell-surface PRRs (Pattern Recognition Receptors) and intracellular receptor NLRs (Nucleotide-Binding Leucine-Rich Repeat) in engineered plants have the potential to activate strong defenses against a broad range of pathogens. Here, we describe the identification, characterization, and in planta functional analysis of a novel truncated NLR (TNx) gene from the wild species Arachis stenosperma (AsTIR19), with a protein structure lacking the C-terminal LRR (Leucine Rich Repeat) domain involved in pathogen perception. Overexpression of AsTIR19 in tobacco plants led to a significant reduction in infection caused by Sclerotinia sclerotiorum, with a further reduction in pyramid lines containing an expansin-like B gene (AdEXLB8) potentially involved in defense priming. Transcription analysis of tobacco transgenic lines revealed induction of hormone defense pathways (SA; JA-ET) and PRs (Pathogenesis-Related proteins) production. The strong upregulation of the respiratory burst oxidase homolog D (RbohD) gene in the pyramid lines suggests its central role in mediating immune responses in plants co-expressing the two transgenes, with reactive oxygen species (ROS) production enhanced by AdEXLB8 cues leading to stronger defense response. Here, we demonstrate that the association of potential priming elicitors and truncated NLRs can produce a synergistic effect on fungal resistance, constituting a promising strategy for improved, non-specific resistance to plant pathogens.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

AnnotaPipeline: An integrated tool to annotate eukaryotic proteins using multi-omics data

Assignment of gene function has been a crucial, laborious, and time-consuming step in genomics. Due to a variety of sequencing platforms that generates increasing amounts of data, manual annotation is no longer feasible. Thus, the need for an integrated, automated pipeline allowing the use of experimental data towards validation of in silico prediction of gene function is of utmost relevance. Here, we present a computational workflow named AnnotaPipeline that integrates distinct software and data types on a proteogenomic approach to annotate and validate predicted features in genomic sequences. Based on FASTA (i) nucleotide or (ii) protein sequences or (iii) structural annotation files (GFF3), users can input FASTQ RNA-seq data, MS/MS data from mzXML or similar formats, as the pipeline uses both transcriptomic and proteomic information to corroborate annotations and validate gene prediction, providing transcription and expression evidence for functional annotation. Reannotation of the available Arabidopsis thaliana, Caenorhabditis elegans, Candida albicans, Trypanosoma cruzi, and Trypanosoma rangeli genomes was performed using the AnnotaPipeline, resulting in a higher proportion of annotated proteins and a reduced proportion of hypothetical proteins when compared to the annotations publicly available for these organisms. AnnotaPipeline is a Unix-based pipeline developed using Python and is available at: https://github.com/bioinformatics-ufsc/AnnotaPipeline.

Molecular Breeding and Drought Tolerance in Chickpea

Cicer arietinum L. is the third greatest widely planted imperative pulse crop worldwide, and it belongs to the Leguminosae family. Drought is the utmost common abiotic factor on plants, distressing their water status and limiting their growth and development. Chickpea genotypes have the natural ability to fight drought stress using certain strategies viz., escape, avoidance and tolerance. Assorted breeding methods, including hybridization, mutation, and marker-aided breeding, genome sequencing along with omics approaches, could be used to improve the chickpea germplasm lines(s) against drought stress. Root features, for instance depth and root biomass, have been recognized as the greatest beneficial morphological factors for managing terminal drought tolerance in the chickpea. Marker-aided selection, for example, is a genomics-assisted breeding (GAB) strategy that can considerably increase crop breeding accuracy and competence. These breeding technologies, notably marker-assisted breeding, omics, and plant physiology knowledge, underlined the importance of chickpea breeding and can be used in future crop improvement programmes to generate drought-tolerant cultivars(s).

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

Normalized Workflow to Optimize Hybrid De Novo Transcriptome Assembly for Non-Model Species: A Case Study in Lilium ledebourii (Baker) Boiss

A high-quality transcriptome is required to advance numerous bioinformatics workflows. Nevertheless, the effectuality of tools for de novo assembly and real precision assembled transcriptomes looks somewhat unexplored, particularly for non-model organisms with complicated (very long, heterozygous, polyploid) genomes. To disclose the performance of various transcriptome assembly programs, this study built 11 single assemblies and analyzed their performance on some significant reference-free and reference-based criteria. As well as to reconfirm the outputs of benchmarks, 55 BLAST were performed and compared using 11 constructed transcriptomes. Concisely, normalized benchmarking demonstrated that Velvet–Oases suffer from the worst results, while the EvidentialGene strategy can provide the most comprehensive and accurate transcriptome of Lilium ledebourii (Baker) Boiss. The BLAST results also confirmed the superiority of EvidentialGene, so it could capture even up to 59% more (than Velvet–Oases) unique gene hits. To promote assembly optimization, with the help of normalized benchmarking, PCA and AHC, it is emphasized that each metric can only provide part of the transcriptome status, and one should never settle for just a few evaluation criteria. This study supplies a framework for benchmarking and optimizing the efficiency of assembly approaches to analyze RNA-Seq data and reveals that selecting an inefficient assembly strategy might result in less identification of unique gene hits.

A companion to the preclinical common data elements for genomics, transcriptomics, and epigenomics data in rodent epilepsy models. A report of the TASK3-WG4 omics working group of the ILAE/AES joint translational TASK force

The International League Against Epilepsy/American Epilepsy Society (ILAE/AES) Joint Translational Task Force established the TASK3 working groups to create common data elements (CDEs) for various preclinical epilepsy research disciplines. The aim of the CDEs is to improve the standardization of experimental designs across a range of epilepsy research-related methods. Here, we have generated CDE tables with key parameters and case report forms (CRFs) containing the essential contents of the study protocols for genomics, transcriptomics, and epigenomics in rodent models of epilepsy, with a specific focus on adult rats and mice. We discuss the important elements that need to be considered for genomics, transcriptomics, and epigenomics methodologies, providing a rationale for the parameters that should be collected. This is the first in a two-part series of omics papers with the second installment to cover proteomics, lipidomics, and metabolomics in adult rodents.

Time Domains of Hypoxia Responses and -Omics Insights

The ability to respond rapidly to changes in oxygen tension is critical for many forms of life. Challenges to oxygen homeostasis, specifically in the contexts of evolutionary biology and biomedicine, provide important insights into mechanisms of hypoxia adaptation and tolerance. Here we synthesize findings across varying time domains of hypoxia in terms of oxygen delivery, ranging from early animal to modern human evolution and examine the potential impacts of environmental and clinical challenges through emerging multi-omics approaches. We discuss how diverse animal species have adapted to hypoxic environments, how humans vary in their responses to hypoxia (i.e., in the context of high-altitude exposure, cardiopulmonary disease, and sleep apnea), and how findings from each of these fields inform the other and lead to promising new directions in basic and clinical hypoxia research.

Full-length transcriptome sequencing analysis and characterization, development and validation of microsatellite markers in Kengyilia melanthera

As a typical psammophyte of the Triticeae, Kengyilia melanthera possesses high feeding potential and great utilization values in desertification control in the Qinghai-Tibet Plateau. However, few gene function and genetic studies have been performed in K. melanthera. In this study, single-molecule real-time sequencing technology was used to obtain the full-length transcriptome sequence of K. melanthera, following the functional annotation of transcripts and prediction of coding sequences (CDSs), transcription factors (TFs), and long noncoding RNA (lncRNA) sequences. Meanwhile, a total of 42,433 SSR loci were detected, with 5'-UTRs having the most SSR loci and trinucleotide being the most abundant type. In total, 108,399 SSR markers were designed, and 300 SSR markers were randomly selected for diversity verification of K. melanthera. A total of 49 polymorphic SSR markers were used to construct the genetic relationships of 56 K. melanthera accessions, among which 21 SSR markers showed good cross-species transferability among the related species. In conclusion, the full-length transcriptome sequence of the K. melanthera will assist gene prediction and promote molecular biology and genomics research, and the polymorphic SSR markers will promote molecular-assisted breeding and related research of K. melanthera and its relatives.

Identifying Terpenoid Biosynthesis Genes in Euphorbia maculata via Full-Length cDNA Sequencing

The annual herb Euphorbia maculata L. produces anti-inflammatory and biologically active substances such as triterpenoids, tannins, and polyphenols, and it is used in traditional Chinese medicine. Of these bioactive compounds, terpenoids, also called isoprenoids, are major secondary metabolites in E. maculata. Full-length cDNA sequencing was carried out to characterize the transcripts of terpenoid biosynthesis reference genes and determine the copy numbers of their isoforms using PacBio SMRT sequencing technology. The Illumina short-read sequencing platform was also employed to identify differentially expressed genes (DEGs) in the secondary metabolite pathways from leaves, roots, and stems. PacBio generated 62 million polymerase reads, resulting in 81,433 high-quality reads. From these high-quality reads, we reconstructed a genome of 20,722 genes, in which 20,246 genes (97.8%) did not have paralogs. About 33% of the identified genes had two or more isoforms. DEG analysis revealed that the expression level differed among gene paralogs in the leaf, stem, and root. Whole sets of paralogs and isoforms were identified in the mevalonic acid (MVA), methylerythritol phosphate (MEP), and terpenoid biosynthesis pathways in the E. maculata L. The nucleotide information will be useful for identifying orthologous genes in other terpenoid-producing medicinal plants.

Full-Length Transcriptome Reconstruction Reveals the Genetic Mechanisms of Eyestalk Displacement and Its Potential Implications on the Interspecific Hybrid Crab (Scylla serrata ♀ × S. paramamosain ♂)

Simple Summary

The eyestalk is a key organ in crustaceans that produces neurohormones and regulates a range of physiological functions. Eyestalk displacement was discovered in some first-generation (F1) offspring of the novel interspecific hybrid crab (Scylla serrata ♀ × S. paramamosain ♂). To uncover the genetic mechanism underlying eyestalk displacement and its potential implications, high-quality transcriptome was reconstructed using single-molecule real-time (SMRT) sequencing. A total of 37 significantly differential alternative splicing (DAS) events (17 up-regulated and 20 down-regulated) and 1475 significantly differential expressed transcripts (DETs) (492 up-regulated and 983 down-regulated) were detected in hybrid crabs with displaced eyestalks (DH). The most significant DAS events and DETs were annotated as being endoplasmic reticulum chaperone BiP and leucine-rich repeat protein lrrA-like isoform X2. In addition, the top ten significant gene ontology (GO) terms were related to the cuticle or chitin. Overall, this study highlights the underlying genetic mechanisms of eyestalk displacement and provide useful knowledge for mud crab (Scylla spp.) crossbreeding.

Abstract

The lack of high-quality juvenile crabs is the greatest impediment to the growth of the mud crab (Scylla paramamosain) industry. To obtain high-quality hybrid offspring, a novel hybrid mud crab (S. serrata ♀ × S. paramamosain ♂) was successfully produced in our previous study. Meanwhile, an interesting phenomenon was discovered, that some first-generation (F1) hybrid offspring’s eyestalks were displaced during the crablet stage I. To uncover the genetic mechanism underlying eyestalk displacement and its potential implications, both single-molecule real-time (SMRT) and Illumina RNA sequencing were implemented. Using a two-step collapsing strategy, three high-quality reconstructed transcriptomes were obtained from purebred mud crabs (S. paramamosain) with normal eyestalks (SPA), hybrid crabs with normal eyestalks (NH), and hybrid crabs with displaced eyestalks (DH). In total, 37 significantly differential alternative splicing (DAS) events (17 up-regulated and 20 down-regulated) and 1475 significantly differential expressed transcripts (DETs) (492 up-regulated and 983 down-regulated) were detected in DH. The most significant DAS events and DETs were annotated as being endoplasmic reticulum chaperone BiP and leucine-rich repeat protein lrrA-like isoform X2. In addition, the top ten significant GO terms were related to the cuticle or chitin. Overall, high-quality reconstructed transcriptomes were obtained for the novel interspecific hybrid crab and provided valuable insights into the genetic mechanisms of eyestalk displacement in mud crab (Scylla spp.) crossbreeding.

The Long Read Transcriptome of Rice (Oryza sativa ssp. japonica var. Nipponbare) Reveals Novel Transcripts

BACKGROUND

High-throughput next-generation sequencing technologies offer a powerful approach to characterizing the transcriptomes of plants. Long read sequencing has been shown to support the discovery of novel isoforms of transcripts. This approach enables the generation of full-length sequences revealing splice variants that may be important in regulating gene action. Investigation of the diversity of transcripts in the rice transcriptome including splice variants was conducted using PacBio long-read sequence data to improve the annotation of the rice genome.

RESULTS

A cDNA library was prepared from RNA extracted from leaves, roots, seeds, inflorescences, and panicles of O. sativa ssp. japonica var Nipponbare and sequenced on a PacBio Sequel platform. This produced 346,190 non-redundant full-length non-chimeric reads (FLNC) resulting in 33,504 high-quality transcripts. Half of the transcripts were multi-exonic and entirely matched with the reference transcripts. However, 14,874 novel isoforms were also identified resulting predominantly from intron retention and at least one novel splice site. Intron retention was the prevalent alternative splicing event and exon skipping was the least observed. Of 73,659 splice junctions, 12,755 (17%) represented novel splice junctions with canonical and non-canonical intron boundaries. The complexity of the transcriptome was examined in detail for 19 starch synthesis-related genes, defining 276 spliced isoforms of which 94 splice variants were novel.

CONCLUSION

The data reveal the great complexity of the rice transcriptome. The novel transcripts provide new insights that may be a key input in future research to improve the annotation of the rice genome.

Identification of genes involved in drought tolerance in seedlings of the desert grass, Psammochloa villosa (Poaceae), based on full-length isoform sequencing and de novo assembly from short reads

Psammochloa villosa is a perennial herbaceous plant that is dominant within arid regions of the Inner Mongolian Plateau and the Qinghai-Tibet Plateau in China, where it is an endemic species and exhibits strong drought tolerance and wind resistance. To study drought tolerance in P. villosa and determine its molecular basis, we simulated high and moderate drought stress in a controlled environment and then analyzed transcriptome sequences by combining long-read sequences from a representative, wild-grown individual with short reads from the treatment groups. We obtained 184,076 high-quality isoforms as a reference and 168,650 genes (91.6%), which we were able to annotate according to public databases. Ultimately, we obtained 119,005 unigenes representing the transcriptome of P. villosa under drought stress and, among these, we identified 3089 differentially expressed genes and 1484 transcription factors. Physiologically, P. villosa that was exposed to high and moderate drought stress had reduced germination rates and shorter buds but generated more chlorophyll, which is atypical under drought stress and possibly reflects an adaptation of these plants to their arid environment. We inferred that significantly upregulated genes were annotated as 'Chlorophyll a-b binding protein' and 'Light-harvesting chlorophyll-protein' among drought and control groups. Broadly, our analyses revealed that drought stress triggered many genome-level responses, especially related to mitigation of radical oxygen species (ROS), which increase in concentration under drought stress. In particular, in the high drought stress group compared with the control, GO enrichment analysis revealed a significant enrichment of upregulated genes (n = 10) involved in mitigation of oxidative stress. Similarly, using KEGG we found significant enrichment of genes in the phenylpropanoid biosynthesis pathway (11 genes), which yields phenols that scavenge ROS. We also inferred that many genes involved in metabolism of arginine and proline, which may serve as both scavengers of ROS and osmoprotectants that interact with stress response genes based on our protein-protein interaction network analysis. We verified the relative expression levels of eight genes associated with mitigation of ROS, DNA repair, and transmembrane transporter activity using qRT-PCR, and the results were consistent with our inferences from transcriptomes. This study provides insights into the genomic and physiological basis of drought tolerance in P. villosa and represents a resource for development of the species as a forage crop via molecular breeding within arid lands.

Full-Length Transcriptome Characterization and Comparative Analysis of Chosenia arbutifolia

As a unique tree species in the Salicaceae family, Chosenia arbutifolia is used primarily for construction materials and landscape planting in China. Compared with other Salicaceae species members, the genomic resources of C. arbutifolia are extremely scarce. Thus, in the present study, the full-length transcriptome of C. arbutifolia was sequenced by single-molecular real-time sequencing (SMRT) technology based on the PacBio platform. Then, it was compared against those of other Salicaceae species. We generated 17,397,064 subreads and 95,940 polished reads with an average length of 1812 bp, which were acquired through calibration, clustering, and polishing. In total, 50,073 genes were reconstructed, of which 48,174 open reading frames, 4281 long non-coding RNAs, and 3121 transcription factors were discovered. Functional annotation revealed that 47,717 genes had a hit in at least one of five reference databases. Moreover, a set of 12,332 putative SSR markers were screened among the reconstructed genes. Single-copy and special orthogroups, and divergent and conserved genes, were identified and analyzed to find divergence among C. arbutifolia and the five Salicaceae species. To reveal genes involved in a specific function and pathway, enrichment analyses for GO and KEGG were also performed. In conclusion, the present study empirically confirmed that SMRT sequencing realistically depicted the C. arbutifolia transcriptome and provided a comprehensive reference for functional genomic research on Salicaceae species.

Sequencing the Origins of Life

One goal of origins of life research is to understand how primitive informational and catalytic biopolymers emerged and evolved. Recently, a number of sequencing techniques have been applied to analysis of replicating and evolving primitive biopolymer systems, providing a sequence-specific and high-resolution view of primitive chemical processes. Here, we review application of sequencing techniques to analysis of synthetic and primitive nucleic acids and polypeptides. This includes next-generation sequencing of primitive polymerization and evolution processes, followed by discussion of other novel biochemical techniques that could contribute to sequence analysis of primitive biopolymer driven chemical systems. Further application of sequencing to origins of life research, perhaps as a life detection technology, could provide insight into the origin and evolution of informational and catalytic biopolymers on early Earth or elsewhere.

RNA sequencing and its applications in cancer and rare diseases

With the invention of RNA sequencing over a decade ago, diagnosis and identification of the gene-related diseases entered a new phase that enabled more accurate analysis of the diseases that are difficult to approach and analyze. RNA sequencing has availed in-depth study of transcriptomes in different species and provided better understanding of rare diseases and taxonomical classifications of various eukaryotic organisms. Development of single-cell, short-read, long-read and direct RNA sequencing using both blood and biopsy specimens of the organism together with recent advancement in computational analysis programs has made the medical professional’s ability in identifying the origin and cause of genetic disorders indispensable. Altogether, such advantages have evolved the treatment design since RNA sequencing can detect the resistant genes against the existing therapies and help medical professions to take a further step in improving methods of treatments towards higher effectiveness and less side effects. Therefore, it is of essence to all researchers and scientists to have deeper insight in all available methods of RNA sequencing while taking a step-in therapy design.

Analysis of Transcriptomic Data Generated from Drug-Treated Cancer Cell Line

Resistance to chemotherapy is one major obstacle in current cancer treatment. Therefore, understanding the molecular basis of the acquisition of resistance is vital for the design and development of appropriate cancer therapy. Importantly, acquisition of resistance is not a single-step process, and the molecular signature of cells dynamically changes during this process. With the advent of next-generation omic technologies, today one can precisely map the molecular alterations not only in a population of tumor cells but also at the single-cell level as they attain chemo-resistance. In this chapter, we describe a detailed transcriptomic pipeline following next-generation sequencing for mapping alteration in expression during the process of attainment of resistance. We provide comprehensive information on the process to (1) track the differential expression of transcripts, (2) understand the gene ontology functions, (3) filter out candidate key genes, (4) identify the pathways regulated by them, and (5) generate a map of their probable interactions. We assume that our analytical method will be useful for research in this direction.

Improved High-Quality Genome Assembly and Annotation of Pineapple (Ananas comosus) Cultivar MD2 Revealed Extensive Haplotype Diversity and Diversified FRS/FRF Gene Family

Pineapple (Ananas comosus (L.) Merr.) is the second most important tropical fruit crop globally, and ‘MD2’ is the most important cultivated variety. A high-quality genome is important for molecular-based breeding, but available pineapple genomes still have some quality limitations. Here, PacBio and Hi-C data were used to develop a new high-quality MD2 assembly and gene prediction. Compared to the previous MD2 assembly, major improvements included a 26.6-fold increase in contig N50 length, phased chromosomes, and >6000 new genes. The new MD2 assembly also included 161.6 Mb additional sequences and >3000 extra genes compared to the F153 genome. Over 48% of the predicted genes harbored potential deleterious mutations, indicating that the high level of heterozygosity in this species contributes to maintaining functional alleles. The genome was used to characterize the FAR1-RELATED SEQUENCE (FRS) genes that were expanded in pineapple and rice. Transposed and dispersed duplications contributed to expanding the numbers of these genes in the pineapple lineage. Several AcFRS genes were differentially expressed among tissue-types and stages of flower development, suggesting that their expansion contributed to evolving specialized functions in reproductive tissues. The new MD2 assembly will serve as a new reference for genetic and genomic studies in pineapple.

Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing

Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community.

•

There is widespread transcript diversity in the cortex and many novel transcripts

•

Some genes display big differences in isoform number between human and mouse cortex

•

There is evidence of differential transcript usage between human fetal and adult cortex

•

There are many novel isoforms of genes associated with human brain disease

Bioinformatic Tools for the Analysis and Prediction of ncRNA Interactions

Noncoding RNAs (ncRNAs) play prominent roles in the regulation of gene expression via their interactions with other biological molecules such as proteins and nucleic acids. Although much of our knowledge about how these ncRNAs operate in different biological processes has been obtained from experimental findings, computational biology can also clearly substantially boost this knowledge by suggesting possible novel interactions of these ncRNAs with other molecules. Computational predictions are thus used as an alternative source of new insights through a process of mutual enrichment because the information obtained through experiments continuously feeds through into computational methods. The results of these predictions in turn shed light on possible interactions that are subsequently validated experimentally. This review describes the latest advances in databases, bioinformatic tools, and new in silico strategies that allow the establishment or prediction of biological interactions of ncRNAs, particularly miRNAs and lncRNAs. The ncRNA species described in this work have a special emphasis on those found in humans, but information on ncRNA of other species is also included.

Hybrid Sequencing in Different Types of Goat Skeletal Muscles Reveals Genes Regulating Muscle Development and Meat Quality

Domestic goats are commonly reared for meat and milk production in several regions of the world. However, the genetic mechanism underlying muscle development and meat quality of goats is limited. Therefore, the aim of this study was to identify known and novel genes regulating muscle development and meat quality of goats using second- and third-generation sequencing technologies. To achieve this, the meat quality and transcriptomes of longissimus dorsi (LD) and biceps femoris (BF) muscle tissues of Lingqiu Greyback goats were examined and compared. Differentially expressed genes (DEGs) and isoforms (DEIs) were functionally annotated. Results showed that 45,574 full-length transcripts covering 18,491 loci were characterized, and 12,566 genes were co-expressed in all samples. Differential expression analysis identified 231 DEGs, including 45 novel genes in the LD and BF muscles of the goats. Additionally, 1173 DEIs were found, in which 642 novel isoforms were identified in this study. Functional annotation and pathway analysis of the DEGs and DEIs revealed that some of them were associated with muscle growth and lipid metabolism. Overall, the findings of this study contribute to the understanding of the transcriptomic diversity underlying meat quality and muscle development of goat.

Overdominance at the Gene Expression Level Plays a Critical Role in the Hybrid Root Growth of Brassica napus

Despite heterosis contributing to genetic improvements in crops, root growth heterosis in rapeseed plants is poorly understood at the molecular level. The current study was performed to discover key differentially expressed genes (DEGs) related to heterosis in two hybrids with contrasting root growth performance (FO; high hybrid and FV; low hybrid) based on analysis of the root heterosis effect. Based on comparative transcriptomic analysis, we believe that the overdominance at the gene expression level plays a critical role in hybrid roots’ early biomass heterosis. Our findings imply that a considerable increase in up-regulation of gene expression underpins heterosis. In the FO hybrid, high expression of DEGs overdominant in the starch/sucrose and galactose metabolic pathways revealed a link between hybrid vigor and root growth. DEGs linked to auxin, cytokinin, brassinosteroids, ethylene, and abscisic acid were also specified, showing that these hormones may enhance mechanisms of root growth and the development in the FO hybrid. Moreover, transcription factors such as MYB, ERF, bHLH, NAC, bZIP, and WRKY are thought to control downstream genes involved in root growth. Overall, this is the first study to provide a better understanding related to the regulation of the molecular mechanism of heterosis, which assists in rapeseed growth and yield improvement.

Long-read transcriptome and other genomic resources for the angiosperm Silene noctiflora

The angiosperm genus Silene is a model system for several traits of ecological and evolutionary significance in plants, including breeding system and sex chromosome evolution, host-pathogen interactions, invasive species biology, heavy metal tolerance, and cytonuclear interactions. Despite its importance, genomic resources for this large genus of approximately 850 species are scarce, with only one published whole-genome sequence (from the dioecious species Silene latifolia). Here, we provide genomic and transcriptomic resources for a hermaphroditic representative of this genus (S. noctiflora), including a PacBio Iso-Seq transcriptome, which uses long-read, single-molecule sequencing technology to analyze full-length mRNA transcripts. Using these data, we have assembled and annotated high-quality full-length cDNA sequences for approximately 14,126 S. noctiflora genes and 25,317 isoforms. We demonstrated the utility of these data to distinguish between recent and highly similar gene duplicates by identifying novel paralogous genes in an essential protease complex. Furthermore, we provide a draft assembly for the approximately 2.7-Gb genome of this species, which is near the upper range of genome-size values reported for diploids in this genus and threefold larger than the 0.9-Gb genome of Silene conica, another species in the same subgenus. Karyotyping confirmed that S. noctiflora is a diploid, indicating that its large genome size is not due to polyploidization. These resources should facilitate further study and development of this genus as a model in plant ecology and evolution.

Long-read RNA sequencing reveals widespread sex-specific alternative splicing in threespine stickleback fish

Alternate isoforms are important contributors to phenotypic diversity across eukaryotes. Although short-read RNA-sequencing has increased our understanding of isoform diversity, it is challenging to accurately detect full-length transcripts, preventing the identification of many alternate isoforms. Long-read sequencing technologies have made it possible to sequence full-length alternative transcripts, accurately characterizing alternative splicing events, alternate transcription start and end sites, and differences in UTR regions. Here, we use Pacific Biosciences (PacBio) long-read RNA-sequencing (Iso-Seq) to examine the transcriptomes of five organs in threespine stickleback fish (Gasterosteus aculeatus), a widely used genetic model species. The threespine stickleback fish has a refined genome assembly in which gene annotations are based on short-read RNA sequencing and predictions from coding sequence of other species. This suggests some of the existing annotations may be inaccurate or alternative transcripts may not be fully characterized. Using Iso-Seq we detected thousands of novel isoforms, indicating many isoforms are absent in the current Ensembl gene annotations. In addition, we refined many of the existing annotations within the genome. We noted many improperly positioned transcription start sites that were refined with long-read sequencing. The Iso-Seq-predicted transcription start sites were more accurate and verified through ATAC-seq. We also detected many alternative splicing events between sexes and across organs. We found a substantial number of genes in both somatic and gonadal samples that had sex-specific isoforms. Our study highlights the power of long-read sequencing to study the complexity of transcriptomes, greatly improving genomic resources for the threespine stickleback fish.

Gene Expression and Isoform Identification of PacBio Full-Length cDNA Sequences for Berberine Biosynthesis in Berberis koreana

Berberis koreana is a medicinal plant containing berberine, which is a bioactive compound of the benzylisoquinoline alkaloid (BIA) class. BIA is widely used in the food and drug industry for its health benefits. To investigate the berberine biosynthesis pathway, gene expression analysis was performed in leaves, flowers, and fruits at different stages of growth. This was followed by full-length cDNA sequencing analysis using the PacBio sequencer platform to determine the number of isoforms of those expressed genes. We identified 23,246 full-length unigenes, among which 8479 had more than one isoform. The number of isoforms ranged between two to thirty-one among all genes. Complete isoform analysis was carried out on the unigenes encoding BIA synthesis. Thirteen of the sixteen genes encoding enzymes for berberine synthesis were present in more than one copy. This demonstrates that gene duplication and translation into isoforms may contribute to the functional specificity of the duplicated genes and isoforms in plant alkaloid synthesis. Our study also demonstrated the streamlining of berberine biosynthesis via the absence of genes for enzymes of other BIAs, but the presence of all the genes for berberine biosynthesize in B. koreana. In addition to genes encoding enzymes for the berberine biosynthesis pathway, the genes encoding enzymes for other BIAs were not present in our dataset except for those encoding corytuberine synthase (CTS) and berbamunine synthase (BS). Therefore, this explains how B. koreana produces berberine by blocking the pathways leading to other BIAs, effectively only allowing the pathway to lead to berberine synthesis.

Characterization of full-length transcriptome in Saccharum officinarum and molecular insights into tiller development

BACKGROUND

Although extensive breeding efforts are ongoing in sugarcane (Saccharum officinarum L.), the average yield is far below the theoretical potential. Tillering is an important component of sugarcane yield, however, the molecular mechanism underlying tiller development is still elusive. The limited genomic data in sugarcane, particularly due to its complex and large genome, has hindered in-depth molecular studies.

RESULTS

Herein, we generated full-length (FL) transcriptome from developing leaf and tiller bud samples based on PacBio Iso-Seq. In addition, we performed RNA-seq from tiller bud samples at three developmental stages (T0, T1 and T2) to uncover key genes and biological pathways involved in sugarcane tiller development. In total, 30,360 and 20,088 high-quality non-redundant isoforms were identified in leaf and tiller bud samples, respectively, representing 41,109 unique isoforms in sugarcane. Likewise, we identified 1063 and 1037 alternative splicing events identified in leaf and tiller bud samples, respectively. We predicted the presence of coding sequence for 40,343 isoforms, 98% of which was successfully annotated. Comparison with previous FL transcriptomes in sugarcane revealed 2963 unreported isoforms. In addition, we characterized 14,946 SSRs from 11,700 transcripts and 310 lncRNAs. By integrating RNA-seq with the FL transcriptome, 468 and 57 differentially expressed genes (DEG) were identified in T1vsT0 and T2vsT0, respectively. Strong up-regulation of several pyruvate phosphate dikinase and phosphoenolpyruvate carboxylase genes suggests enhanced carbon fixation and protein synthesis to facilitate tiller growth. Similarly, up-regulation of linoleate 9S-lipoxygenase and lipoxygenase genes in the linoleic acid metabolism pathway suggests high synthesis of key oxylipins involved in tiller growth and development.

CONCLUSIONS

Collectively, we have enriched the genomic data available in sugarcane and provided candidate genes for manipulating tiller formation and development, towards productivity enhancement in sugarcane.

Aspergillus derived mycotoxins in food and the environment: Prevalence, detection, and toxicity

Aspergillus species are the paramount ubiquitous fungi that contaminate various food substrates and produce biochemicals known as mycotoxins. Aflatoxins (AFTs), ochratoxin A (OTA), patulin (PAT), citrinin (CIT), aflatrem (AT), secalonic acids (SA), cyclopiazonic acid (CPA), terrein (TR), sterigmatocystin (ST) and gliotoxin (GT), and other toxins produced by species of Aspergillus plays a major role in food and human health. Mycotoxins exhibited wide range of toxicity to the humans and animal models even at nanomolar (nM) concentration. Consumption of detrimental mycotoxins adulterated foodstuffs affects human and animal health even trace amounts. Bioaerosols consisting of spores and hyphal fragments are active elicitors of bronchial irritation and allergy, and challenging to the public health. Aspergillus is the furthermost predominant environmental contaminant unswervingly defile lives with a 40-90 % mortality risk in patients with conceded immunity. Genomics, proteomics, transcriptomics, and metabolomics approaches useful for mycotoxins' detection which are expensive. Antibody based detection of toxins chemotypes may result in cross-reactivity and uncertainty. Aptamers (APT) are single stranded DNA (ssDNA/RNA), are specifically binds to the target molecules can be generated by systematic evolution of ligands through exponential enrichment (SELEX). APT are fast, sensitive, simple, in-expensive, and field-deployable rapid point of care (POC) detection of toxins, and a better alternative to antibodies.

High-throughput transcriptome sequencing and comparative analysis of Escherichia coli and Schizosaccharomyces pombe in respiratory and fermentative growth

In spite of increased complexity in eukaryotes compared to prokaryotes, several basic metabolic and regulatory processes are conserved. Here we explored analogies in the eubacteria Escherichia coli and the unicellular fission yeast Schizosaccharomyces pombe transcriptomes under two carbon sources: 2% glucose; or a mix of 2% glycerol and 0.2% sodium acetate using the same growth media and growth phase. Overall, twelve RNA-seq libraries were constructed. A total of 593 and 860 genes were detected as differentially expressed for E. coli and S. pombe, respectively, with a log2 of the Fold Change ≥ 1 and False Discovery Rate ≤ 0.05. In aerobic glycolysis, most of the expressed genes were associated with cell proliferation in both organisms, including amino acid metabolism and glycolysis. In contrast in glycerol/acetate condition, genes related to flagellar assembly and membrane proteins were differentially expressed such as the general transcription factors fliA, flhD, flhC, and flagellum assembly genes were detected in E. coli, whereas in S. pombe genes for hexose transporters, integral membrane proteins, galactose metabolism, and ncRNAs related to cellular stress were overexpressed. In general, our study shows that a conserved "foraging behavior" response is observed in these eukaryotic and eubacterial organisms in gluconeogenic carbon sources.

Data-driven biomarker analysis using computational omics approaches to assess neurodegenerative disease progression

The complexity of biological systems suggests that current definitions of molecular dysfunctions are essential distinctions of a complex phenotype. This is well seen in neurodegenerative diseases (ND), such as Alzheimer's disease (AD) and Parkinson's disease (PD), multi-factorial pathologies characterized by high heterogeneity. These challenges make it necessary to understand the effectiveness of candidate biomarkers for early diagnosis, as well as to obtain a comprehensive mapping of how selective treatment alters the progression of the disorder. A large number of computational methods have been developed to explain network-based approaches by integrating individual components for modeling a complex system. In this review, high-throughput omics methodologies are presented for the identification of potent biomarkers associated with AD and PD pathogenesis as well as for monitoring the response of dysfunctional molecular pathways incorporating multilevel clinical information. In addition, principles for efficient data analysis pipelines are being discussed that can help address current limitations during the experimental process by increasing the reproducibility of benchmarking studies.

Genomics Armed With Diversity Leads the Way in Brassica Improvement in a Changing Global Environment

Meeting the needs of a growing world population in the face of imminent climate change is a challenge; breeding of vegetable and oilseed Brassica crops is part of the race in meeting these demands. Available genetic diversity constituting the foundation of breeding is essential in plant improvement. Elite varieties, land races, and crop wild species are important resources of useful variation and are available from existing genepools or genebanks. Conservation of diversity in genepools, genebanks, and even the wild is crucial in preventing the loss of variation for future breeding efforts. In addition, the identification of suitable parental lines and alleles is critical in ensuring the development of resilient Brassica crops. During the past two decades, an increasing number of high-quality nuclear and organellar Brassica genomes have been assembled. Whole-genome re-sequencing and the development of pan-genomes are overcoming the limitations of the single reference genome and provide the basis for further exploration. Genomic and complementary omic tools such as microarrays, transcriptomics, epigenetics, and reverse genetics facilitate the study of crop evolution, breeding histories, and the discovery of loci associated with highly sought-after agronomic traits. Furthermore, in genomic selection, predicted breeding values based on phenotype and genome-wide marker scores allow the preselection of promising genotypes, enhancing genetic gains and substantially quickening the breeding cycle. It is clear that genomics, armed with diversity, is set to lead the way in Brassica improvement; however, a multidisciplinary plant breeding approach that includes phenotype = genotype × environment × management interaction will ultimately ensure the selection of resilient Brassica varieties ready for climate change.

BSR-Seq analysis provides insights into the cold stress response of Actinidia arguta F1 populations

Background

Freezing injury, which is an important abiotic stress in horticultural crops, influences the growth and development and the production area of kiwifruit (Actinidia Lind1). Among Actinidia species, Actinidia arguta has excellent cold resistance, but knowledge relevant to molecular mechanisms is still limited. Understanding the mechanism underlying cold resistance in kiwifruit is important for breeding cold resistance.

Results

In our study, a population resulting from the cross of A. arguta ‘Ruby-3’ × ‘Kuilv’ male was generated for kiwifruit hardiness study, and 20 cold-tolerant and 20 cold-sensitive populations were selected from 492 populations according to their LT50. Then, we performed bulked segregant RNA-seq combined with single-molecule real-time sequencing to identify differentially expressed genes that provide cold hardiness. We found that the content of soluble sucrose and the activity of β-amylase were higher in the cold-tolerant population than in the cold-sensitive population. Upon − 30 °C low-temperature treatment, 126 differentially expressed genes were identify; the expression of 59 genes was up-regulated and that of 67 genes was down-regulated between the tolerant and sensitive pools, respectively. KEGG pathway analysis showed that the DEGs were primarily related to starch and sucrose metabolism, amino sugar and nucleotide sugar metabolism. Ten major key enzyme-encoding genes and two regulatory genes were up-regulated in the tolerant pool, and regulatory genes of the CBF pathway were found to be differentially expressed. In particular, a 14–3-3 gene was down-regulated and an EBF gene was up-regulated. To validate the BSR-Seq results, 24 DEGs were assessed via qRT-PCR, and the results were consistent with those obtained by BSR-Seq.

Conclusion

Our research provides valuable insights into the mechanism related to cold resistance in Actinidia and identified potential genes that are important for cold resistance in kiwifruit.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07369-9.

Whole blood transcriptome analysis using RNA sequencing in individuals with insomnia disorder and good sleepers: a pilot study

BACKGROUND

Insomnia is a highly prevalent condition that is associated with negative health outcomes, yet little is known about the underlying molecular mechanisms.

METHOD

RNA sequencing was conducted using blood samples from 15 individuals with primary insomnia and 15 age- and gender-matched good sleeper controls. The RNA library was sequenced with 150 base pair paired-ends on the Illumina NovaSeq-6000 platform. Alignment was performed using human reference genome hg38. Differential gene expression analysis was performed using DESeq2 following alignment, using log fold change ±0.50, and had a false discovery rate p-value <0.05. Pathway analysis was performed using Ingenuity Pathway Analysis.

RESULTS

We found 288 differentially expressed genes in insomnia patients when compared to controls. Upregulated genes included LINC02224 (Long Intergenic Non-Protein Coding RNA 2224), DUX4L9 (Double Homeobox 4 Like 9), and TUSC3 (Tumor Suppressor Candidate 3) and down regulated genes included CTXN2 (Cortexin 2), CSMD1 (CUB And Sushi Multiple Domains 1), and SLC12A1 (Solute Carrier Family 12 Member 1). Ingenuity® Pathway Analysis (IPA) revealed 3 associated networks (score>40) with genes and hubs related to inflammation (nuclear factor-kB), oxidative stress (Mitochondrial complex 1) and ubiquitination.

CONCLUSION

Differentially expressed genes in this analysis are functionally associated with inflammation and immune response, mitochondrial and metabolic processes. Further research into the transcriptomic changes in insomnia is needed to understand related pathways to the disorder and provide new avenues for diagnostics and therapeutics.