Alternative expression analysis by RNA sequencing

The Domestication of Wild Boar Could Result in a Relaxed Selection for Maintaining Olfactory Capacity

Domesticated animals are artificially selected to exhibit desirable traits, however not all traits of domesticated animals are the result of deliberate selection. Loss of olfactory capacity in the domesticated pig (Sus scrofa domesticus) is one example. We used whole transcriptome analysis (RNA-Seq) to compare patterns of gene expression in the olfactory mucosa of the pig and two subspecies of wild boar (Sus scrofa), and investigate candidate genes that could be responsible for the loss of olfactory capacity. We identified hundreds of genes with reductions in transcript abundance in pig relative to wild boar as well as differences between the two subspecies of wild boar. These differences were detected mainly in genes involved in the formation and motility of villi, cilia and microtubules, functions associated with olfaction. In addition, differences were found in the abundances of transcripts of genes related to immune defenses, with the highest levels in continental wild boar subspecies. Overall, the loss of olfactory capacity in pigs appears to have been accompanied by reductions in the expression of candidate genes for olfaction. These changes could have resulted from unintentional selection for reduced olfactory capacity, relaxed selection for maintaining olfactory capacity, pleiotropic effects of genes under selection, or other non-selective processes. Our findings could be a cornerstone for future researches on wild boars, pigs, feral populations, and their evolutionary trajectories, aimed to provide tools to better calibrate species management as well as guidelines for breeders.

Differences in Oxygen-Induced Retinopathy Susceptibility Between Two Sprague Dawley Rat Vendors: A Comparison of Retinal Transcriptomes

PURPOSE

To determine the retinal transcriptomic differences underlying the oxygen-induced retinopathy phenotypes between Sprague Dawley rat pups from two commonly used commercial vendors. This will allow us to discover genes and pathways that may be related to differences in disease severity in similarly aged premature babies and suggest possible new treatment approaches.

METHODS

We analyzed retinal vascular morphometry and transcriptomes from Sprague Dawley rat pups from Charles River Laboratories and Envigo (previously Harlan). Room air control and oxygen-induced retinopathy groups were compared. Oxygen-induced retinopathy was induced with the rat 50/10 model.

RESULTS

Pups from Charles River Laboratories developed a more severe oxygen-induced retinopathy phenotype, with 3.6-fold larger percent avascular area at P15 and twofold larger % neovascular area at P20 than pups from Envigo. Changes in retinal transcriptomes of rat pups from both vendors were substantial at baseline and in response to oxygen-induced retinopathy. Baseline differences centered on activated pathways of neuronal development in Charles River Laboratories pups. In response to oxygen-induced retinopathy, during the neovascular phase, retinas from Charles River Laboratories pups exhibited activation of pathways regulating necrosis, neuroinflammation, and interferon signaling, supporting the observed increase of neovascularization. Conversely, retinas from Envigo pups showed decreased necrosis and increased focal adhesion kinase signaling, supporting more normal vascular development. Comparing oxygen-induced retinopathy transcriptomes at P15 to those at P20, canonical pathways such as phosphate and tensin homolog, interferon, and coordinated lysosomal expression and regulation element signaling were identified, highlighting potential novel mechanistic targets for future research.

CONCLUSION

Transcriptomic profiles differ substantially between rat pup retinas from Charles River Laboratories and Envigo at baseline and in response to oxygen-induced retinopathy, providing insight into vascular morphologic differences. Comparing transcriptomes identified new pathways for further research in oxygen-induced retinopathy pathogenesis and increased scientific rigor of this model.

Metagenomics and metatranscriptomics suggest pathways of 3-chloroaniline degradation in wastewater reactors

Biological wastewater treatment systems are often affected by shifts in influent quality, including the input of toxic chemicals. Yet the mechanisms underlying the adaptation of activated sludge process performance are rarely studied in a controlled and replicated experimental setting, particularly when challenged with a sustained toxin input. Three replicate bench-scale bioreactors were subjected to a chemical disturbance in the form of 3-chloroaniline (3-CA) over 132 days, after an acclimation period of 58 days, while three control reactors received no 3-CA input. Ammonia oxidation was initially affected by 3-CA. Within three weeks of the experiment, microbial communities in all three treatment reactors adapted to biologically degrade 3-CA resulting in partial ammonia oxidation recovery. Combining process and microbial community data from amplicon sequencing with potential functions gleaned from assembled metagenomics and metatranscriptomics data, two putative degradation pathways for 3-CA were identified. The first pathway, determined from metagenomics data, involves a benzoate dioxygenase and subsequent meta-cleavage of the aromatic ring. The second, determined from intensive short-term sampling for gene expression data in tandem with 3-CA degradation, involves a phenol monooxygenase followed by ortho-cleavage of the aromatic ring. The relative abundances of amplicon sequence variants associated with the genera Gemmatimonas, OLB8, and Taibaiella correlated significantly with 3-CA degradation. Metagenome-assembled genome data also showed the genus OLB8 to be differentially enriched in treatment reactors, making it a strong candidate as 3-CA degrader. Using replicated reactors, this study has demonstrated the impact of a sustained stress on the activated sludge process. The unique and novel features of this study include the identification of putative pathways and potential degraders of 3-CA using long-term and short-term sampling in tandem with multiple methods in a controlled and replicated experiment.

Transcriptome-Powered Pluripotent Stem Cell Differentiation for Regenerative Medicine

Pluripotent stem cells are endless sources for in vitro engineering human tissues for regenerative medicine. Extensive studies have demonstrated that transcription factors are the key to stem cell lineage commitment and differentiation efficacy. As the transcription factor profile varies depending on the cell type, global transcriptome analysis through RNA sequencing (RNAseq) has been a powerful tool for measuring and characterizing the success of stem cell differentiation. RNAseq has been utilized to comprehend how gene expression changes as cells differentiate and provide a guide to inducing cellular differentiation based on promoting the expression of specific genes. It has also been utilized to determine the specific cell type. This review highlights RNAseq techniques, tools for RNAseq data interpretation, RNAseq data analytic methods and their utilities, and transcriptomics-enabled human stem cell differentiation. In addition, the review outlines the potential benefits of the transcriptomics-aided discovery of intrinsic factors influencing stem cell lineage commitment, transcriptomics applied to disease physiology studies using patients' induced pluripotent stem cell (iPSC)-derived cells for regenerative medicine, and the future outlook on the technology and its implementation.

Host-pathogen interaction involving cytoskeleton changes as well as non-coding regulation as primary mechanisms for SRS resistance in Atlantic salmon

The salmonid rickettsial syndrome (SRS) is a systemic bacterial infection caused by Piscirickettsia salmonis that generates significant economic losses in Atlantic salmon (Salmo salar) aquaculture. Despite this disease's relevance, the mechanisms involved in resistance against P. salmonis infection are not entirely understood. Thus, we aimed at studying the pathways explaining SRS resistance using different approaches. First, we determined the heritability using pedigree data from a challenge test. Secondly, a genome-wide association analysis was performed following a complete transcriptomic profile of fish from genetically susceptible and resistant families within the challenge infection with P. salmonis. We found differentially expressed transcripts related to immune response, pathogen recognition, and several new pathways related to extracellular matrix remodelling and intracellular invasion. The resistant background showed a constrained inflammatory response, mediated by the Arp2/3 complex actin cytoskeleton remodelling polymerization pathway, probably leading to bacterial clearance. A series of biomarkers of SRS resistance, such as the beta-enolase (ENO-β), Tubulin G1 (TUBG1), Plasmin (PLG) and ARP2/3 Complex Subunit 4 (ARPC4) genes showed consistent overexpression in resistant individuals, showing promise as biomarkers for SRS resistance. All these results together with the differential expression of several long non-coding RNAs show the complexity of the host-pathogen interaction of S. salar and P. salmonis. These results provide valuable information on new models describing host-pathogen interaction and its role in SRS resistance.

Characterizing nutritional phenotypes using experimental nutrigenomics: Is there nutrient-specificity to different types of dietary stress?

The ability to directly measure and monitor poor nutrition in individual animals and ecological communities is hampered by methodological limitations. In this study, we use nutrigenomics to identify nutritional biomarkers in a freshwater zooplankter, Daphnia pulex, a ubiquitous primary consumer in lakes and a sentinel of environmental change. We grew animals in six ecologically relevant nutritional treatments: nutrient replete, low carbon (food), low phosphorus, low nitrogen, low calcium and high Cyanobacteria. We extracted RNA for transcriptome sequencing to identify genes that were nutrient responsive and capable of predicting nutritional status with a high degree of accuracy. We selected a list of 125 candidate genes, which were subsequently pruned to 13 predictive potential biomarkers. Using a nearest-neighbour classification algorithm, we demonstrate that these potential biomarkers are capable of classifying our samples into the correct nutritional group with 100% accuracy. The functional annotation of the selected biomarkers revealed some specific nutritional pathways and supported our hypothesis that animal responses to poor nutrition are nutrient specific and not simply different presentations of slow growth or energy limitation. This is a key step in uncovering the causes and consequences of nutritional limitation in animal consumers and their responses to small- and large-scale changes in biogeochemical cycles.

Transcriptomic responses of females to consumption of nuptial food gifts as a potential mediator of sexual conflict in decorated crickets

Nuptial food gift provisioning by males to females at mating is a strategy in many insects that is thought to be shaped by sexual conflict or sexual selection, as it affords males access to a female's physiology. While males often attempt to use these gifts to influence female behaviour to their own advantage, females can evolve counter mechanisms. In decorated crickets, the male's nuptial gift comprises part of the spermatophore, the spermatophylax, the feeding on which deters the female from prematurely terminating sperm transfer. However, ingested compounds in the spermatophylax and attachment of the sperm-containing ampulla could further influence female physiology and behaviour. We investigated how mating per se and these two distinct routes of potential male-mediated manipulation influence the female transcriptomic response. We conducted an RNA sequencing experiment on gut and head tissues from females for whom nuptial food gift consumption and receipt of an ejaculation were independently manipulated. In the gut tissue, we found that females not permitted to feed during mating exhibited decreased overall gene expression, possibly caused by a reduced gut function, but this was countered by feeding on the spermatophylax or a sham gift. In the head tissue, we found only low numbers of differentially expressed genes, but a gene co-expression network analysis revealed that ampulla attachment and spermatophylax consumption independently induce distinct gene expression patterns. This study provides evidence that spermatophylax feeding alters the female post-mating transcriptomic response in decorated crickets, highlighting its potential to mediate sexual conflict in this system.

Comparative transcriptomics of growth metabolism and virulence reveal distinct morphogenic profiles of yeast-like cells and hyphae of the fungus Metarhizium rileyi

Metarhizium rileyiis an entomopathogenic fungus with a narrow host range which distinguishes it from other Metarhiziumspecies with broad host ranges. This species is also unique because the initial yeast-like growth on solid media is only observed in liquid culture in other Metharizium species. A lack of knowledge about the metabolism and genetic signatures of M. rileyiduring this yeast-like phase on solid and in liquid media is a bottleneck for its large-scale production as a commercial biocontrol agent.In this study wefound that M. rileyiyeast-like cells produced on solid medium infected and killed the important insect pest Spodoptera frugiperda with comparable efficiency as yeast-like cells grown in liquid medium. Secondly, we used comparative transcriptomic analysis to investigate theactive genes and genomic signatures of the M. rileyi yeast-like morphotypes produced on solid and in liquid media. Yeast-like cells grown in liquid medium had upregulated genes relating specifically to signal transduction andparticular membrane transporters. Thirdly, we compared the transcriptomic profiles of yeast-like phases of M. rileyi with those of M. anisopliae. The yeast-like phase of M. rileyi grown on solid medium upregulated unique genes not found in otherMetarhiziumspecies including specific membrane proteins and several virulence factors. Orthologous genes associated with heat shock protein, iron permease, membrane proteins and key virulence traits (e.g. collagen-like protein Mcl1) were upregulated in both species. Comparative transcriptome analyses of gene expression showed more differences than similarities between M. anisopliae and M. rileyi yeast-like cells.

Selection of Ideal Reference Genes for Gene Expression Analysis in COVID-19 and Mucormycosis

Selection of reference genes during real-time quantitative PCR (qRT-PCR) is critical to determine accurate and reliable mRNA expression. Nonetheless, not a single study has investigated the expression stability of candidate reference genes to determine their suitability as internal controls in SARS-CoV-2 infection or COVID-19-associated mucormycosis (CAM). Using qRT-PCR, we determined expression stability of the nine most commonly used housekeeping genes, namely, TATA-box binding protein (TBP), cyclophilin (CypA), β-2-microglobulin (B2M), 18S rRNA (18S), peroxisome proliferator-activated receptor gamma (PPARG) coactivator 1 alpha (PGC-1α), glucuronidase beta (GUSB), hypoxanthine phosphoribosyltransferase 1 (HPRT-1), β-ACTIN, and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) in patients with COVID-19 of various severities (asymptomatic, mild, moderate, and severe) and those with CAM. We used statistical algorithms (delta-CT [threshold cycle], NormFinder, BestKeeper, GeNorm, and RefFinder) to select the most appropriate reference gene and observed that clinical severity profoundly influences expression stability of reference genes. CypA demonstrated the most consistent expression irrespective of disease severity and emerged as the most suitable reference gene in COVID-19 and CAM. Incidentally, GAPDH, the most commonly used reference gene, showed the maximum variations in expression and emerged as the least suitable. Next, we determined expression of nuclear factor erythroid 2-related factor 2 (NRF2), interleukin-6 (IL-6), and IL-15 using CypA and GAPDH as internal controls and show that CypA-normalized expression matches well with the RNA sequencing-based expression of these genes. Further, IL-6 expression correlated well with the plasma levels of IL-6 and C-reactive protein, a marker of inflammation. In conclusion, GAPDH emerged as the least suitable and CypA as the most suitable reference gene in COVID-19 and CAM. The results highlight the expression variability of housekeeping genes due to disease severity and provide a strong rationale for identification of appropriate reference genes in other chronic conditions as well. IMPORTANCE Gene expression studies are critical to develop new diagnostics, therapeutics, and prognostic modalities. However, accurate determination of expression requires data normalization with a reference gene, whose expression does not vary across different disease stages. Misidentification of a reference gene can produce inaccurate results. Unfortunately, despite the global impact of COVID-19 and an urgent unmet need for better treatment, not a single study has investigated the expression stability of housekeeping genes across the disease spectrum to determine their suitability as internal controls. Our study identifies CypA and then TBP as the two most suitable reference genes for COVID-19 and CAM. Further, GAPDH, the most commonly used reference gene in COVID-19 studies, turned out to be the least suitable. This work fills an important gap in the field and promises to facilitate determination of an accurate expression of genes to catalyze development of novel molecular diagnostics and therapeutics for improved patient care.

Comprehensive and scalable quantification of splicing differences with MntJULiP

Tools for differential splicing detection have failed to provide a comprehensive and consistent view of splicing variation. We present MntJULiP, a novel method for comprehensive and accurate quantification of splicing differences between two or more conditions. MntJULiP detects both changes in intron splicing ratios and changes in absolute splicing levels with high accuracy, and can find classes of variation overlooked by other tools. MntJULiP identifies over 29,000 differentially spliced introns in 1398 GTEx brain samples, including 11,242 novel introns discovered in this dataset. Highly scalable, MntJULiP can process thousands of samples within hours to reveal splicing constituents of phenotypic differentiation.

Fracture hematoma micro-architecture influences transcriptional profile and plays a crucial role in determining bone healing outcomes

The hematoma that forms between broken fragments of bone serves as a natural fibrin scaffold, and its removal from the defect site delays bone healing. The hypothesis of this study is that the microarchitectural and mechanical properties of the initially formed hematoma has a significant effect on the regulation of the biological process, which ultimately determines the outcome of bone healing. To mimic three healing conditions in the rat femur (normal, delayed, and non-healing bone defects), three different defect sizes of 0.5, 1.5, and 5.0 mm, are respectively used. The analysis of 3-day-old hematomas demonstrates clear differences in fibrin clot micro-architecture in terms of fiber diameter, fiber density, and porosity of the formed fibrin network, which result in different mechanical properties (stiffness) of the hematoma in each model. Those differences directly affect the biological processes involved. Specifically, RNA-sequencing reveals almost 700 differentially expressed genes between normally healing and non-healing defects, including significantly up-regulated essential osteogenic genes in normally healing defects, also differences in immune cell populations, activated osteogenic transcriptional regulators as well as potential novel marker genes. Most importantly, this study demonstrates that the healing outcome has already been determined during the hematoma phase of bone healing, three days post-surgery.

Comparative transcriptome analysis reveals the resistance regulation mechanism and fungicidal activity of the fungicide phenamacril in Fusarium oxysporum

Fusarium oxysporum (Fo) is an important species complex of soil-borne pathogenic fungi that cause vascular wilt diseases of agricultural crops and some opportunistic diseases of humans. The fungicide phenamacril has been extensively reported to have antifungal activity against Fusarium graminearum and Fusarium fujikuroi. In this study, we found that the amino acid substitutions (V151A and S418T) in Type I myosin FoMyo5 cause natural low resistance to phenamacril in the plant pathogenic Fo isolates. Therefore, we compared the transcriptomes of two phenamacril-resistant Fo isolates FoII5, Fo1st and one phenamacril-sensitive isolate Fo3_a after 1 μg/mL phenamacril treatment. Among the 2728 differentially expressed genes (DEGs), 14 DEGs involved in oxidation–reduction processes and MFS transporters, were significantly up-regulated in phenamacril-resistant isolates. On the other hand, 14 DEGs involved in ATP-dependent RNA helicase and ribosomal biogenesis related proteins, showed significantly down-regulated expression in both phenamacril-resistant and -sensitive isolates. These results indicated that phenamacril not only seriously affected the cytoskeletal protein binding and ATPase activity of sensitive isolate, but also suppressed ribosome biogenesis in all the isolates. Hence, this study helps us better understand resistance regulation mechanism and fungicidal activity of phenamacril and provide reference for the development of new fungicides to control Fo.

Single-Cell RNA-Seq Analysis of Cells from Degenerating and Non-Degenerating Intervertebral Discs from the Same Individual Reveals New Biomarkers for Intervertebral Disc Degeneration

In this study, we used single-cell transcriptomic analysis to identify new specific biomarkers for nucleus pulposus (NP) and inner annulus fibrosis (iAF) cells, and to define cell populations within non-degenerating (nD) and degenerating (D) human intervertebral discs (IVD) of the same individual. Cluster analysis based on differential gene expression delineated 14 cell clusters. Gene expression profiles at single-cell resolution revealed the potential functional differences linked to degeneration, and among NP and iAF subpopulations. GO and KEGG analyses discovered molecular functions, biological processes, and transcription factors linked to cell type and degeneration state. We propose two lists of biomarkers, one as specific cell type, including C2orf40, MGP, MSMP, CD44, EIF1, LGALS1, RGCC, EPYC, HILPDA, ACAN, MT1F, CHI3L1, ID1, ID3 and TMED2. The second list proposes predictive IVD degeneration genes, including MT1G, SPP1, HMGA1, FN1, FBXO2, SPARC, VIM, CTGF, MGST1, TAF1D, CAPS, SPTSSB, S100A1, CHI3L2, PLA2G2A, TNRSF11B, FGFBP2, MGP, SLPI, DCN, MT-ND2, MTCYB, ADIRF, FRZB, CLEC3A, UPP1, S100A2, PRG4, COL2A1, SOD2 and MT2A. Protein and mRNA expression of MGST1, vimentin, SOD2 and SYF2 (p29) genes validated our scRNA-seq findings. Our data provide new insights into disc cells phenotypes and biomarkers of IVD degeneration that could improve diagnostic and therapeutic options.

RNA-Seq is not required to determine stable reference genes for qPCR normalization

Assessment of differential gene expression by qPCR is heavily influenced by the choice of reference genes. Although numerous statistical approaches have been proposed to determine the best reference genes, they can give rise to conflicting results depending on experimental conditions. Hence, recent studies propose the use of RNA-Seq to identify stable genes followed by the application of different statistical approaches to determine the best set of reference genes for qPCR data normalization. In this study, however, we demonstrate that the statistical approach to determine the best reference genes from commonly used conventional candidates is more important than the preselection of ‘stable’ candidates from RNA-Seq data. Using a qPCR data normalization workflow that we have previously established; we show that qPCR data normalization using conventional reference genes render the same results as stable reference genes selected from RNA-Seq data. We validated these observations in two distinct cross-sectional experimental conditions involving human iPSC derived microglial cells and mouse sciatic nerves. These results taken together show that given a robust statistical approach for reference gene selection, stable genes selected from RNA-Seq data do not offer any significant advantage over commonly used reference genes for normalizing qPCR assays.

RTqPCR is a powerful technique that is widely used to quantify gene expression in research and diagnostics of different diseases. The technique involves making multiple copies (amplification) of a specific target DNA. The amplified target DNA binds to a molecule that emits fluorescence upon binding. The extent of fluorescence correlates to the amount of DNA present. To precisely quantify this fluorescence (and thus the quantities of target DNA), internal control genes also called as reference genes need to be determined. Such genes, in principle, do not have varied expression across samples and would exhibit the same fluorescence in all samples. They can thus be used to normalize the expression of the Target DNA. Unfortunately, choosing the right reference gene is very tricky and poor choice of reference genes results in unreliable data both in research and in diagnostics. In this study, we validate a statistical approach to find stably expressed reference genes for any experimental setting using a given set of candidates. We compare our approach to RNA sequencing which quantifies the expression of thousands of genes at the same time. We highlight the advantages of our approach which is cost effective and saves a lot of time when compared to sequencing.

Differential gene expression associated with a floral scent polymorphism in the evening primrose Oenothera harringtonii (Onagraceae)

Background

Plant volatiles play an important role in both plant-pollinator and plant-herbivore interactions. Intraspecific polymorphisms in volatile production are ubiquitous, but studies that explore underlying differential gene expression are rare. Oenothera harringtonii populations are polymorphic in floral emission of the monoterpene (R)-(−)-linalool; some plants emit (R)-(−)-linalool (linalool+ plants) while others do not (linalool- plants). However, the genes associated with differential production of this floral volatile in Oenothera are unknown. We used RNA-Seq to broadly characterize differential gene expression involved in (R)-(−)-linalool biosynthesis. To identify genes that may be associated with the polymorphism for this trait, we used RNA-Seq to compare gene expression in six different Oenothera harringtonii tissues from each of three linalool+ and linalool- plants.

Results

Three clusters of differentially expressed genes were enriched for terpene synthase activity: two were characterized by tissue-specific upregulation and one by upregulation only in plants with flowers that produce (R)-(−)-linalool. A molecular phylogeny of all terpene synthases identified two putative (R)-(−)-linalool synthase transcripts in Oenothera harringtonii, a single allele of which is found exclusively in linalool+ plants.

Conclusions

By using a naturally occurring polymorphism and comparing different tissues, we were able to identify candidate genes putatively involved in the biosynthesis of (R)-(−)-linalool. Expression of these genes in linalool- plants, while low, suggests a regulatory polymorphism, rather than a population-specific loss-of-function allele. Additional terpene biosynthesis-related genes that are up-regulated in plants that emit (R)-(−)-linalool may be associated with herbivore defense, suggesting a potential economy of scale between plant reproduction and defense.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08370-6.

Interspecies transcriptomics identify genes that underlie disproportionate foot growth in jerboas

Despite the great diversity of vertebrate limb proportion and our deep understanding of the genetic mechanisms that drive skeletal elongation, little is known about how individual bones reach different lengths in any species. Here, we directly compare the transcriptomes of homologous growth cartilages of the mouse (Mus musculus) and bipedal jerboa (Jaculus jaculus), the latter of which has "mouse-like" arms but extremely long metatarsals of the feet. Intersecting gene-expression differences in metatarsals and forearms of the two species revealed that about 10% of orthologous genes are associated with the disproportionately rapid elongation of neonatal jerboa feet. These include genes and enriched pathways not previously associated with endochondral elongation as well as those that might diversify skeletal proportion in addition to their known requirements for bone growth throughout the skeleton. We also identified transcription regulators that might act as "nodes" for sweeping differences in genome expression between species. Among these, Shox2, which is necessary for proximal limb elongation, has gained expression in jerboa metatarsals where it has not been detected in other vertebrates. We show that Shox2 is sufficient to increase mouse distal limb length, and a nearby putative cis-regulatory region is preferentially accessible in jerboa metatarsals. In addition to mechanisms that might directly promote growth, we found evidence that jerboa foot elongation may occur in part by de-repressing latent growth potential. The genes and pathways that we identified here provide a framework to understand the modular genetic control of skeletal growth and the remarkable malleability of vertebrate limb proportion.

The emerging roles of NGS in clinical oncology and personalized medicine

Next-generation sequencing (NGS) has been increasingly popular in genomics studies over the last decade, as new sequencing technology has been created and improved. Recently, NGS started to be used in clinical oncology to improve cancer therapy through diverse modalities ranging from finding novel and rare cancer mutations, discovering cancer mutation carriers to reaching specific therapeutic approaches known as personalized medicine (PM). PM has the potential to minimize medical expenses by shifting the current traditional medical approach of treating cancer and other diseases to an individualized preventive and predictive approach. Currently, NGS can speed up in the early diagnosis of diseases and discover pharmacogenetic markers that help in personalizing therapies. Despite the tremendous growth in our understanding of genetics, NGS holds the added advantage of providing more comprehensive picture of cancer landscape and uncovering cancer development pathways. In this review, we provided a complete overview of potential NGS applications in scientific and clinical oncology, with a particular emphasis on pharmacogenomics in the direction of precision medicine treatment options.

Genome-wide analysis of Claviceps paspali: insights into the secretome of the main species causing ergot disease in Paspalum spp

BACKGROUND

The phytopatogen Claviceps paspali is the causal agent of Ergot disease in Paspalum spp., which includes highly productive forage grasses such as P. dilatatum. This disease impacts dairy and beef production by affecting seed quality and producing mycotoxins that can affect performance in feeding animals. The molecular basis of pathogenicity of C. paspali remains unknown, which makes it more difficult to find solutions for this problem. Secreted proteins are related to fungi virulence and can manipulate plant immunity acting on different subcellular localizations. Therefore, identifying and characterizing secreted proteins in phytopathogenic fungi will provide a better understanding of how they overcome host defense and cause disease. The aim of this work is to analyze the whole genome sequences of three C. paspali isolates to obtain a comparative genome characterization based on possible secreted proteins and pathogenicity factors present in their genome. In planta RNA-seq analysis at an early stage of the interaction of C. paspali with P. dilatatum stigmas was also conducted in order to determine possible secreted proteins expressed in the infection process.

RESULTS

C. paspali isolates had compact genomes and secretome which accounted for 4.6-4.9% of the predicted proteomes. More than 50% of the predicted secretome had no homology to known proteins. RNA-Seq revealed that three protein-coding genes predicted as secreted have mayor expression changes during 1 dpi vs 4 dpi. Also, three of the first 10 highly expressed genes in both time points were predicted as effector-like. CAZyme-like proteins were found in the predicted secretome and the most abundant family could be associated to pectine degradation. Based on this, pectine could be a main component affected by the cell wall degrading enzymes of C. paspali.

CONCLUSIONS

Based on predictions from DNA sequence and RNA-seq, unique probable secreted proteins and probable pathogenicity factors were identified in C. paspali isolates. This information opens new avenues in the study of the biology of this fungus and how it modulates the interaction with its host. Knowledge of the diversity of the secretome and putative pathogenicity genes should facilitate future research in disease management of Claviceps spp.

Alternative RNA splicing in stem cells and cancer stem cells: Importance of transcript-based expression analysis

Alternative ribonucleic acid (RNA) splicing can lead to the assembly of different protein isoforms with distinctive functions. The outcome of alternative splicing (AS) can result in a complete loss of function or the acquisition of new functions. There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells (CSCs), and their prospective contribution in cancer progression. AS directly regulates the self-renewal features of stem cells (SCs) and stem-like cancer cells. Notably, octamer-binding transcription factor 4A spliced variant of octamer-binding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs. The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness. The alternative spliced variants of CSCs markers, including cluster of differentiation 44, aldehyde dehydrogenase, and doublecortin-like kinase, α6β1 integrin, have pivotal roles in increasing self-renewal properties and maintaining the pluripotency of CSCs. Various splicing analysis tools are considered in this study. LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events. Additionally, LeafCutter can be used for efficient mapping splicing quantitative trait loci. Altogether, the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

Global Transcriptome Characterization and Assembly of the Thermophilic Ascomycete Chaetomium thermophilum

A correct genome annotation is fundamental for research in the field of molecular and structural biology. The annotation of the reference genome of Chaetomium thermophilum has been reported previously, but it is essentially limited to open reading frames (ORFs) of protein coding genes and contains only a few noncoding transcripts. In this study, we identified and annotated full-length transcripts of C. thermophilum by deep RNA sequencing. We annotated 7044 coding genes and 4567 noncoding genes. Astonishingly, 23% of the coding genes are alternatively spliced. We identified 679 novel coding genes as well as 2878 novel noncoding genes and corrected the structural organization of more than 50% of the previously annotated genes. Furthermore, we substantially extended the Gene Ontology (GO) and Enzyme Commission (EC) lists, which provide comprehensive search tools for potential industrial applications and basic research. The identified novel transcripts and improved annotation will help to understand the gene regulatory landscape in C. thermophilum. The analysis pipeline developed here can be used to build transcriptome assemblies and identify coding and noncoding RNAs of other species.

Predicting Agents That Can Overcome 5-FU Resistance in Colorectal Cancers via Pharmacogenomic Analysis

5-Fluorouracil (5-FU) is one of several chemotherapeutic agents in clinical use as a standard of care to treat colorectal cancers (CRCs). As an antimetabolite, 5-FU inhibits thymidylate synthase to disrupt the synthesis and repair of DNA and RNA. However, only a small proportion of patients benefit from 5-FU treatment due to the development of drug resistance. This study applied pharmacogenomic analysis using two public resources, the Genomics of Drug Sensitivity in Cancer (GDSC) and the Connectivity Map, to predict agents overcoming 5-FU resistance in CRC cells based on their genetic background or gene expression profile. Based on the genetic status of adenomatous polyposis coli (APC), the most frequent mutated gene found in CRC, we found that combining a MEK inhibitor with 5-FU exhibited synergism effects on CRC cells with APC truncations. While considering the gene expression in 5-FU resistant cells, we demonstrated that targeting ROCK is a potential avenue to restore 5-FU response to resistant cells with wild-type APC background. Our results reveal MEK signaling plays a pivotal role in loss-of-function, APC-mediated 5-FU resistance, and ROCK activation serves as a signature in APC-independent 5-FU resistance. Through the use of these available database resources, we highlight possible approaches to predict potential drugs for combinatorial therapy for patients developing resistance to 5-FU treatment.

Gene and protein expression in human megakaryocytes derived from induced pluripotent stem cells

BACKGROUND

There is interest in deriving megakaryocytes (MKs) from pluripotent stem cells (iPSC) for biological studies. We previously found that genomic structural integrity and genotype concordance is maintained in iPSC-derived MKs.

OBJECTIVE

To establish a comprehensive dataset of genes and proteins expressed in iPSC-derived MKs.

METHODS

iPSCs were reprogrammed from peripheral blood mononuclear cells (MNCs) and MKs were derived from the iPSCs in 194 healthy European American and African American subjects. mRNA was isolated and gene expression measured by RNA sequencing. Protein expression was measured in 62 of the subjects using mass spectrometry.

RESULTS AND CONCLUSIONS

MKs expressed genes and proteins known to be important in MK and platelet function and demonstrated good agreement with previous studies in human MKs derived from CD34+ progenitor cells. The percent of cells expressing the MK markers CD41 and CD42a was consistent in biological replicates, but variable across subjects, suggesting that unidentified subject-specific factors determine differentiation of MKs from iPSCs. Gene and protein sets important in platelet function were associated with increasing expression of CD41/42a, while those related to more basic cellular functions were associated with lower CD41/42a expression. There was differential gene expression by the sex and race (but not age) of the subject. Numerous genes and proteins were highly expressed in MKs but not known to play a role in MK or platelet function; these represent excellent candidates for future study of hematopoiesis, platelet formation, and/or platelet function.

Lamin A/C Is Dispensable to Mechanical Repression of Adipogenesis

Mesenchymal stem cells (MSCs) maintain the musculoskeletal system by differentiating into multiple lineages, including osteoblasts and adipocytes. Mechanical signals, including strain and low-intensity vibration (LIV), are important regulators of MSC differentiation via control exerted through the cell structure. Lamin A/C is a protein vital to the nuclear architecture that supports chromatin organization and differentiation and contributes to the mechanical integrity of the nucleus. We investigated whether lamin A/C and mechanoresponsiveness are functionally coupled during adipogenesis in MSCs. siRNA depletion of lamin A/C increased the nuclear area, height, and volume and decreased the circularity and stiffness. Lamin A/C depletion significantly decreased markers of adipogenesis (adiponectin, cellular lipid content) as did LIV treatment despite depletion of lamin A/C. Phosphorylation of focal adhesions in response to mechanical challenge was also preserved during loss of lamin A/C. RNA-seq showed no major adipogenic transcriptome changes resulting from LIV treatment, suggesting that LIV regulation of adipogenesis may not occur at the transcriptional level. We observed that during both lamin A/C depletion and LIV, interferon signaling was downregulated, suggesting potentially shared regulatory mechanism elements that could regulate protein translation. We conclude that the mechanoregulation of adipogenesis and the mechanical activation of focal adhesions function independently from those of lamin A/C.

Alternative splicing acts as an independent prognosticator in ovarian carcinoma

Alternative splicing (AS) events associated with oncogenic processes present anomalous perturbations in many cancers, including ovarian carcinoma. There are no reliable features to predict survival outcomes for ovarian cancer patients. In this study, comprehensive profiling of AS events was conducted by integrating AS data and clinical information of ovarian serous cystadenocarcinoma (OV). Survival-related AS events were identified by Univariate Cox regression analysis. Then, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis were used to construct the prognostic signatures within each AS type. Furthermore, we established a splicing-related network to reveal the potential regulatory mechanisms between splicing factors and candidate AS events. A total of 730 AS events were identified as survival-associated splicing events, and the final prognostic signature based on all seven types of AS events could serve as an independent prognostic indicator and had powerful efficiency in distinguishing patient outcomes. In addition, survival-related AS events might be involved in tumor-related pathways including base excision repair and pyrimidine metabolism pathways, and some splicing factors might be correlated with prognosis-related AS events, including SPEN, SF3B5, RNPC3, LUC7L3, SRSF11 and PRPF38B. Our study constructs an independent prognostic signature for predicting ovarian cancer patients’ survival outcome and contributes to elucidating the underlying mechanism of AS in tumor development.

Effects of purposeful soccer heading on circulating small extracellular vesicle concentration and cargo

BACKGROUND

Considering the potential cumulative effects of repetitive head impact (HI) exposure, we need sensitive biomarkers to track short- and long-term effects. Circulating small extracellular vesicles (sEVs) (<200 nm) traffic biological molecules throughout the body and may have diagnostic value as biomarkers for disease. The purpose of this study was to identify the microRNA (miRNA) profile in circulating sEVs derived from human plasma following repetitive HI exposure.

METHODS

Healthy adult (aged 18-35 years) soccer players were randomly assigned to one of 3 groups: the HI group performed 10 standing headers, the leg impact group performed 10 soccer ball trapping maneuvers over 10 min, and the control group did not participate in any soccer drills. Plasma was collected before testing and 24 h afterward, and sEVs were isolated and characterized via nanoparticle tracking analysis. Next-generation sequencing was utilized to identify candidate miRNAs isolated from sEVs, and candidate microRNAs were analyzed via quantitative polymerase chain reaction. In silico target prediction was performed using TargetScan (Version 7.0; targetscan.org) and miRWalk (http://mirwalk.umm.uni-heidelberg.de/) programs, and target validation was performed using luciferase reporter vectors with a miR-7844-5p mimic in human embryonic kidney (HEK) 293T/17 cells.

RESULTS

Plasma sEV concentration and size were not affected across time and group following repetitive HI exposure. After 24 h, the HI read count from next-generation sequencing showed a 4-fold or greater increase in miR-92b-5p, miR-423-5p, and miR-24-3p and a 3-fold or greater decrease in miR-7844-5p, miR-144-5p, miR-221-5p, and miR-22-3p. Analysis of quantitative polymerase chain reaction revealed that leg impact did not alter the candidate miRNA levels. To our knowledge, miR-7844-5p is a previously unknown miRNA. We identified 8 miR-7844-5p mRNA targets: protein phosphatase 1 regulatory inhibitor subunit 1B (PPP1R1B), LIM and senescent cell antigen-like domains 1 (LIMS1), autophagy-related 12 (ATG12), microtubule-associated protein 1 light chain 3 beta (MAP1LC3B), integrin subunit alpha-1 (ITGA1), mitogen-activated protein kinase 1 (MAPK1), glycogen synthase kinase 3β (GSK3β), and mitogen-activated protein kinase 8 (MAPK8).

CONCLUSION

Collectively, these data indicate repetitive HI exposure alters plasma sEV miRNA content, but not sEV size or number. Furthermore, for the first time we demonstrate that previously unknown miR-7844-5p targets mRNAs known to be involved in mitochondrial apoptosis, autophagy regulation, mood disorders, and neurodegenerative disease.

Decoding the Reproductive System of the Olive Fruit Fly, Bactrocera oleae

In most diploid organisms, mating is a prerequisite for reproduction and, thus, critical to the maintenance of their population and the perpetuation of the species. Besides the importance of understanding the fundamentals of reproduction, targeting the reproductive success of a pest insect is also a promising method for its control, as a possible manipulation of the reproductive system could affect its destructive activity. Here, we used an integrated approach for the elucidation of the reproductive system and mating procedures of the olive fruit fly, Bactrocera oleae. Initially, we performed a RNAseq analysis in reproductive tissues of virgin and mated insects. A comparison of the transcriptomes resulted in the identification of genes that are differentially expressed after mating. Functional annotation of the genes showed an alteration in the metabolic, catalytic, and cellular processes after mating. Moreover, a functional analysis through RNAi silencing of two differentially expressed genes, yellow-g and troponin C, resulted in a significantly reduced oviposition rate. This study provided a foundation for future investigations into the olive fruit fly's reproductive biology to the development of new exploitable tools for its control.

Characterization of lncRNA and mRNA profiles in rats with diabetic macroangiopathy

Diabetic macroangiopathy is part of the most common serious complications of diabetes. Previous studies indicate that lncRNAs involved in the process of diabetes and another vascular disease. However, their detailed mechanism of the lncRNAs involved in diabetic macroangiopathy has not been well characterized. In the present study, we generated rat models of diabetic macroangiopathy induced by High fat of 16weeks. A total of 15 GK rats were constructed as a test group, along with 15 Wistar rats set as control group, and thoracic aorta tissue from each group was collected. Whole genomic RNA sequencing was performed on thoracic aorta tissue; 3223 novel lncRNAs and 20367 annotated lncRNAs were indemnified in thoracic aorta samples, and 864 lncRNAs were expressed differently in the test and control groups. Gene ontology term enrichment showed the apparent enrichment of inflammatory response and cell apoptosis, which consistent with the results of H&E Staining, TUNEL Assay, and ELISA; Extensive literature reveals inflammatory response and cell apoptosis play an important role in the process of diabetic macroangiopathy. The results of the present study indicated that lncRNAs, especially Nrep. bSep08, Col5a1, aSep0, soygee.aSep08-unspliced, NONRATT013247.2, votar.aSep08-unspliced, etc, both participate in and mediate the process of inflammatory response, cell apoptosis. What’s more. Our research provides further insights into understanding of the basic molecular mechanisms underlying diabetic macroangiopathy.

Multi-Tissue Multi-Omics Nutrigenomics Indicates Context-Specific Effects of Docosahexaenoic Acid on Rat Brain

SCOPE

The influence of docosahexaenoic acid (DHA) on cardiometabolic and cognitive phenotypes, and multi-omic alterations in the brain under two metabolic conditions is explored to understand context-specific nutritional effects.

METHODS AND RESULTS

Rats are randomly assigned to a DHA-rich or a control chow diet while drinking water or high fructose solution, followed by profiling of metabolic and cognitive phenotypes and the transcriptome and DNA methylome of the hypothalamus and hippocampus. DHA reduces serum triglyceride and improves insulin resistance and memory exclusively in the fructose-consuming rats. In hippocampus, DHA affects genes related to synapse functions in the chow group but immune functions in the fructose group; in hypothalamus, DHA alters immune pathways in the chow group but metabolic pathways in the fructose group. Network modeling reveals context-specific regulators of DHA effects, including Klf4 and Dusp1 for chow condition and Lum, Fn1, and Col1a1 for fructose condition in hippocampus, as well as Cyr61, JunB, Ier2, and Pitx2 under chow condition and Hcar1, Cdh1, and Osr1 under fructose condition in hypothalamus.

CONCLUSION

DHA exhibits differential influence on epigenetic loci, genes, pathways, and metabolic and cognitive phenotypes under different dietary contexts, supporting population stratification in DHA studies to achieve precision nutrition.

Carbohydrate Accumulation and Differential Transcript Expression in Winter Wheat Lines with Different Levels of Snow Mold and Freezing Tolerance after Cold Treatment

Winter wheat (Triticum aestivum L.) undergoes a period of cold acclimation in order to survive the ensuing winter, which can bring freezing temperatures and snow mold infection. Tolerance of these stresses is conferred in part by accumulation of carbohydrates in the crown region. This study investigates the contributions of carbohydrate accumulation during a cold treatment among wheat lines that differ in their snow mold tolerance (SMT) or susceptibility (SMS) and freezing tolerance (FrT) or susceptibility (FrS). Two parent varieties and eight recombinant inbred lines (RILs) were analyzed. The selected RILs represent four combinations of tolerance: SMT/FrT, SMT/FrS, SMS/FrT, and SMS/FrS. It is hypothesized that carbohydrate accumulation and transcript expression will differ between sets of RILs. Liquid chromatography with a refractive index detector was used to quantify carbohydrate content at eight time points over the cold treatment period. Polysaccharide and sucrose content differed between SMT and SMS RILs at various time points, although there were no significant differences in glucose or fructose content. Glucose and fructose content differed between FrT and FrS RILs in this study, but no significant differences in polysaccharide or sucrose content. RNAseq was used to investigate differential transcript expression, followed by modular enrichment analysis, to reveal potential candidates for other mechanisms of tolerance, which included expected pathways such as oxidative stress, chitinase activity, and unexpected transcriptional pathways. These differences in carbohydrate accumulation and differential transcript expression begin to give insight into the differences of wheat lines when exposed to cold temperatures.

Identification and validation of critical alternative splicing events and splicing factors in gastric cancer progression

Gene expression and alternative splicing (AS) interact in complex ways to regulate biological process which is associated with cancer development. Here, by integrated analysis of gene expression and AS events, we aimed to identify the hub AS events and splicing factors relevant in gastric cancer development (GC). RNA‐seq data, clinical data and AS events of 348 GC samples were obtained from the TCGA and TCGASpliceSeq databases. Cox univariable and multivariable analyses, KEGG and GO pathway analyses were performed to identify hub AS events and splicing factor/spliceosome genes, which were further validated in 53 GCs. By bioinformatics methods, we found that gene AS event‐ and gene expression‐mediated GC progression shared the same mechanisms, such as PI3K/AKT pathway, but the involved genes were different. Though expression of 17 hub AS events were confirmed in 53 GC tissues, only 10 AS events in seven genes were identified as critical candidates related to GC progression, notably the AS events (Exon Skip) in CLSTN1 and SEC16A. Expression of these AS events in GC correlated with activation of the PI3K/AKT pathway. Genes with AS events associated with clinical parameters and prognosis were different from the genes whose mRNA levels were related to clinical parameters and prognosis. Besides, we further revealed that QKI and NOVA1 were the crucial splicing factors regulating expression of AS events in GC, but not spliceosome genes. Our integrated analysis revealed hub AS events in GC development, which might be the potential therapeutic targets for GC.

Genome-Wide Profiling of Alternative Splicing Signature Reveals Prognostic Predictor for Esophageal Carcinoma

Background

Alternative splicing (AS) is a molecular event that drives protein diversity through the generation of multiple mRNA isoforms. Growing evidence demonstrates that dysregulation of AS is associated with tumorigenesis. However, an integrated analysis in identifying the AS biomarkers attributed to esophageal carcinoma (ESCA) is largely unexplored.

Methods

AS percent-splice-in (PSI) data were obtained from the TCGA SpliceSeq database. Univariate and multivariate Cox regression analysis was successively performed to identify the overall survival (OS)-associated AS events, followed by the construction of AS predictor through different splicing patterns. Then, a nomogram that combines the final AS predictor and clinicopathological characteristics was established. Finally, a splicing regulatory network was created according to the correlation between the AS events and the splicing factors (SF).

Results

We identified a total of 2389 AS events with the potential to be used as prognostic markers that are associated with the OS of ESCA patients. Based on splicing patterns, we then built eight AS predictors that are highly capable in distinguishing high- and low-risk patients, and in predicting ESCA prognosis. Notably, the area under curve (AUC) value for the exon skip (ES) prognostic predictor was shown to reach a score of 0.885, indicating that ES has the highest prediction strength in predicting ESCA prognosis. In addition, a nomogram that comprises the pathological stage and risk group was shown to be highly efficient in predicting the survival possibility of ESCA patients. Lastly, the splicing correlation network analysis revealed the opposite roles of splicing factors (SFs) in ESCA.

Conclusion

In this study, the AS events may provide reliable biomarkers for the prognosis of ESCA. The splicing correlation networks could provide new insights in the identification of potential regulatory mechanisms during the ESCA development.

Comparative RNAseq Analysis of the Insect-Pathogenic Fungus Metarhizium anisopliae Reveals Specific Transcriptome Signatures of Filamentous and Yeast-Like Development

The fungus Metarhizium anisopliae is a facultative insect pathogen used as biological control agent of several agricultural pests worldwide. It is a dimorphic fungus that is able to display two growth morphologies, a filamentous phase with formation of hyphae and a yeast-like phase with formation of single-celled blastospores. Blastospores play an important role for M. anisopliae pathogenicity during disease development. They are formed solely in the hemolymph of infected insects as a fungal strategy to quickly multiply and colonize the insect's body. Here, we use comparative genome-wide transcriptome analyses to determine changes in gene expression between the filamentous and blastospore growth phases in vitro to characterize physiological changes and metabolic signatures associated with M. anisopliae dimorphism. Our results show a clear molecular distinction between the blastospore and mycelial phases. In total 6.4% (n = 696) out of 10,981 predicted genes in M. anisopliae were differentially expressed between the two phases with a fold-change > 4. The main physiological processes associated with up-regulated gene content in the single-celled yeast-like blastospores during liquid fermentation were oxidative stress, amino acid metabolism (catabolism and anabolism), respiration processes, transmembrane transport and production of secondary metabolites. In contrast, the up-regulated gene content in hyphae were associated with increased growth, metabolism and cell wall re-organization, which underlines the specific functions and altered growth morphology of M. anisopliae blastospores and hyphae, respectively. Our study revealed significant transcriptomic differences between the metabolism of blastospores and hyphae. These findings illustrate important aspects of fungal morphogenesis in M. anisopliae and highlight the main metabolic activities of each propagule under in vitro growth conditions.

Gene Biomarkers Derived from Clinical Data of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a common cancer of high mortality, mainly due to the difficulty in diagnosis during its clinical stage. Here we aim to find the gene biomarkers, which are of important significance for diagnosis and treatment. In this work, 3682 differentially expressed genes on HCC were firstly differentiated based on the Cancer Genome Atlas database (TCGA). Co-expression modules of these differentially expressed genes were then constructed based on the weighted correlation network algorithm. The correlation coefficient between the co-expression module and clinical data from the Broad GDAC Firehose was thereafter derived. Finally, the interactive network of genes was then constructed. Then, the hub genes were used to implement enrichment analysis and pathway analysis in the Database for Annotation, Visualization and Integrated Discovery (DAVID) database. Results revealed that the abnormally expressed genes in the module played an important role in the biological process including cell division, sister chromatid cohesion, DNA repair, and G1/S transition of mitotic cell cycle. Meanwhile, these genes also enriched in a few crucial pathways related to Cell cycle, Oocyte meiosis, and p53 signaling. Via investigating the closeness centrality of the interactive network, eight gene biomarkers including the CKAP2, TPX2, CDCA8, KIFC1, MELK, SGO1, RACGAP1, and KIAA1524 gene were discovered, whose functions had been indeed revealed to be correlated with HCC. This study, therefore, suggests that the abnormal expression of those eight genes may be taken as gene biomarkers of HCC.

Loss of lenalidomide-induced megakaryocytic differentiation leads to therapy resistance in del(5q) myelodysplastic syndrome

Interstitial deletion of the long arm of chromosome 5 (del(5q)) is the most common structural genomic variant in myelodysplastic syndromes (MDS)¹. Lenalidomide (LEN) is the treatment of choice for patients with del(5q) MDS, but half of the responding patients become resistant² within 2 years. TP53 mutations are detected in ~20% of LEN-resistant patients³. Here we show that patients who become resistant to LEN harbour recurrent variants of TP53 or RUNX1. LEN upregulated RUNX1 protein and function in a CRBN- and TP53-dependent manner in del(5q) cells, and mutation or downregulation of RUNX1 rendered cells resistant to LEN. LEN induced megakaryocytic differentiation of del(5q) cells followed by cell death that was dependent on calpain activation and CSNK1A1 degradation^4,5. We also identified GATA2 as a LEN-responsive gene that is required for LEN-induced megakaryocyte differentiation. Megakaryocytic gene-promoter analyses suggested that LEN-induced degradation of IKZF1 enables a RUNX1-GATA2 complex to drive megakaryocytic differentiation. Overexpression of GATA2 restored LEN sensitivity in the context of RUNX1 or TP53 mutations by enhancing LEN-induced megakaryocytic differentiation. Screening for mutations that block LEN-induced megakaryocytic differentiation should identify patients who are resistant to LEN.

Systematic evaluation of differential splicing tools for RNA-seq studies

Differential splicing (DS) is a post-transcriptional biological process with critical, wide-ranging effects on a plethora of cellular activities and disease processes. To date, a number of computational approaches have been developed to identify and quantify differentially spliced genes from RNA-seq data, but a comprehensive intercomparison and appraisal of these approaches is currently lacking. In this study, we systematically evaluated 10 DS analysis tools for consistency and reproducibility, precision, recall and false discovery rate, agreement upon reported differentially spliced genes and functional enrichment. The tools were selected to represent the three different methodological categories: exon-based (DEXSeq, edgeR, JunctionSeq, limma), isoform-based (cuffdiff2, DiffSplice) and event-based methods (dSpliceType, MAJIQ, rMATS, SUPPA). Overall, all the exon-based methods and two event-based methods (MAJIQ and rMATS) scored well on the selected measures. Of the 10 tools tested, the exon-based methods performed generally better than the isoform-based and event-based methods. However, overall, the different data analysis tools performed strikingly differently across different data sets or numbers of samples.

Transcriptomic analysis of Macrobrachium rosenbergii (giant fresh water prawn) post-larvae in response to M. rosenbergii nodavirus (MrNV) infection: de novo assembly and functional annotation

Background

Macrobrachium rosenbergii, is one of a major freshwater prawn species cultured in Southeast Asia. White tail disease (WTD), caused by Macrobrachium rosenbergii nodavirus (MrNV), is a serious problem in farm cultivation and is responsible for up to 100% mortality in the post larvae stage. Molecular data on how M. rosenbergii post-larvae launches an immune response to an infection with MrNV is not currently available. We therefore compared the whole transcriptomic sequence of M. rosenbergii post-larvae before and after MrNV infection.

Results

Transcriptome for M. rosenbergii post-larvae demonstrated high completeness (BUSCO Complete: 83.4%, fragmentation: 13%, missing:3.3%, duplication:16.2%; highest ExN50 value: 94%). The assembled transcriptome consists of 96,362 unigenes with N₅₀ of 1308 bp. The assembled transcriptome was successfully annotated against the NCBI non-redundant arthropod database (33.75%), UniProt database (26.73%), Gene Ontology (GO) (18.98%), Evolutionary Genealogy of Genes: Non-supervised Orthologous Groups (EggNOG) (20.88%), and Kyoto Encyclopedia of Genes and Genome pathway (KEGG) (20.46%). GO annotations included immune system process, signaling, response to stimulus, and antioxidant activity. Differential abundance analysis using EdgeR showed 2413 significantly up-regulated genes and 3125 significantly down-regulated genes during the infection of MrNV.

Conclusions

This study reported a highly complete transcriptome from the post-larvae stage of giant river prawn, M. rosenbergii. Differential abundant transcripts during MrNV infection were identified and validated by qPCR, many of these differentially abundant transcripts as key players in antiviral immunity. These include known members of the innate immune response with the largest expression change occurring in the M. rosenbergii post-larvae after MrNV infection such as antiviral protein, C-type lectin, prophenol oxidase, caspase, ADP ribosylation factors, and dicer.

Using Large Datasets to Understand Nanotechnology

Advances in sequencing technologies have made studying biological processes with genomics, transcriptomics, and proteomics commonplace. As a result, this suite of increasingly integrated techniques is well positioned to study drug delivery, a process that is influenced by many biomolecules working in concert. Omics-based approaches can be used to study the vast nanomaterial chemical space as well as the biological factors that affect the safety, toxicity, and efficacy of nanotechnologies. The generation and analysis of large datasets, methods to interpret them, and dataset applications to nanomaterials to date, are demonstrated here. Finally, new approaches for how sequencing-generated datasets can answer fundamental questions in nanotechnology based drug delivery are proposed.

Genome-wide profiling of prognosis-related alternative splicing signatures in sarcoma

BACKGROUND

Sarcomas (SARCs) are rare malignant tumors with poor prognosis. Increasing evidence has suggested that aberrant alternative splicing (AS) is strongly associated with tumor initiation and progression. We considered whether survival-related AS events might serve as prognosis predictors and underlying targeted molecules in SARC treatment.

METHODS

RNA-Seq data of the SARC cohort were downloaded from The Cancer Genome Atlas (TCGA) database. Survival-related AS events were selected by univariate and multivariate Cox regression analyses. Metascape was used for constructing a gene interaction network and performing functional enrichment analysis. Then, prognosis predictors were established based on statistically significant survival-related AS events and evaluated by receiver operator characteristic (ROC) curve analysis. Finally, the potential regulatory network was analyzed via Pearson's correlation between survival-related AS events and splicing factors (SFs).

RESULTS

A total of 3,610 AS events and 2,291 genes were found to be prognosis-related in 261 SARC samples. The focal adhesion pathway was identified as the most critical molecular mechanism corresponding to poor prognosis. Notably, several prognosis predictors based on survival-related AS events showed excellent performance in prognosis prediction. The area under the curve of the ROC of the risk score was 0.85 in the integrated predictor. The splicing network proved complicated regulation between prognosis-related SFs and AS events. Also, driver gene mutations were significantly associated with AS in SARC patients.

CONCLUSIONS

Survival-related AS events may become ideal indictors for the prognosis prediction of SARCs. Corresponding splicing regulatory mechanisms are worth further exploration.

Determining the prognostic significance of alternative splicing events in soft tissue sarcoma using data from The Cancer Genome Atlas

BACKGROUND

Surgery, adjuvant chemotherapy, and radiotherapy are the primary treatment options for soft tissue sarcomas (STSs). However, identifying ways to improve the prognosis of patients with STS remains a considerable challenge. Evidence shows that the dysregulation of alternative splicing (AS) events is involved in tumor pathogenesis and progression. The present study objective was to identify survival-associated AS events that could serve as prognostic biomarkers and potentially serve as tumor-selective STS drug targets.

METHODS

STS-specific 'percent spliced in' (PSI) values for splicing events in 206 STS samples were downloaded from The Cancer Genome Atlas SpliceSeq^® database. Prognostic analyses were performed on seven types of AS events to determine their prognostic value in STS patients, for which prediction models were constructed with the risk score formula [Formula: see text]. Prediction models were also constructed to determine the prognostic value of AS events, and Spearman's rank correlation coefficients were calculated to determine the degree of correlation between splicing factor expression and the PSI values.

RESULTS

A total 10,439 events were found to significantly correlate with patient survival rates. The area under the time-dependent receiver operating characteristic curve for the prognostic predictor of STS overall survival was 0.826. Notably, the splicing events of certain STS key genes were significantly associated with STS 2-year overall survival in the present study, including exon skip (ES) events in MDM2 and EWSR1, alternate terminator events in CDKN2A and HMGA2 for dedifferentiated liposarcoma, ES in MDM2 and alternate promoter events in CDKN2A for leiomyosarcoma, and ES in EWSR1 for undifferentiated pleomorphic sarcoma. Moreover, splicing correlation networks between AS events and splicing factors revealed that almost all of the AS events showed negatively correlations with the expression of splicing factors.

CONCLUSION

An in-depth analysis of alternative RNA splicing could provide new insights into the mechanisms of STS oncogenesis and the potential for novel approaches to this type of cancer therapy.

Reconstruction of the full-length transcriptome atlas using PacBio Iso-Seq provides insight into the alternative splicing in Gossypium australe

BACKGROUND

Gossypium australe F. Mueller (2n = 2x = 26, G₂ genome) possesses valuable characteristics. For example, the delayed gland morphogenesis trait causes cottonseed protein and oil to be edible while retaining resistance to biotic stress. However, the lack of gene sequences and their alternative splicing (AS) in G. australe remain unclear, hindering to explore species-specific biological morphogenesis.

RESULTS

Here, we report the first sequencing of the full-length transcriptome of the Australian wild cotton species, G. australe, using Pacific Biosciences single-molecule long-read isoform sequencing (Iso-Seq) from the pooled cDNA of ten tissues to identify transcript loci and splice isoforms. We reconstructed the G. australe full-length transcriptome and identified 25,246 genes, 86 pre-miRNAs and 1468 lncRNAs. Most genes (12,832, 50.83%) exhibited two or more isoforms, suggesting a high degree of transcriptome complexity in G. australe. A total of 31,448 AS events in five major types were found among the 9944 gene loci. Among these five major types, intron retention was the most frequent, accounting for 68.85% of AS events. 29,718 polyadenylation sites were detected from 14,536 genes, 7900 of which have alternative polyadenylation sites (APA). In addition, based on our AS events annotations, RNA-Seq short reads from germinating seeds showed that differential expression of these events occurred during seed germination. Ten AS events that were randomly selected were further confirmed by RT-PCR amplification in leaf and germinating seeds.

CONCLUSIONS

The reconstructed gene sequences and their AS in G. australe would provide information for exploring beneficial characteristics in G. australe.

Characterization of Long Noncoding RNA and mRNA Profiles in Sepsis-Induced Myocardial Depression

Septic shock with heart dysfunction is very common in intensive care units. However, whether long noncoding RNA (lncRNA) and mRNA profiles differ between patients with and without myocardial depression is unknown. We generated rat models of hypodynamic septic shock induced by lipopolysaccharide. A total of 12 rat models was constructed and heart tissue from each was collected. Whole genomic RNA sequencing was performed on left ventricular tissue; 6,508 novel lncRNAs and 432 annotated lncRNAs were identified in heart samples, and 74 lncRNAs were expressed differently in the sepsis and control groups. Gene ontology term enrichment indicated apoptosis and its related pathways showed obvious enrichment, which suggested cell apoptosis could play a critical role in the process of myocardial depression. Furthermore, we focused on one lncRNA from the Pvt1 gene. By silencing this lncRNA, we demonstrated knockdown of Pvt1 expression could induce cell apoptosis in lipopolysaccharide-induced heart cells, through increasing the expression of c-Myc, Bid, Bax, and caspase-3 and decreasing the expression of Myd88 and Bcl-2, thereby proving its functional role in myocardial depression. These results demonstrate that lncRNAs both participate in and mediate the pathological process of myocardial depression. Our study improves the understanding of the basic molecular mechanisms underlying myocardial depression.

Genome-wide identification and characterization of long non-coding RNAs expressed during sheep fetal and postnatal hair follicle development

Long non-coding RNAs (lncRNAs), >200 nt in length, are transcribed from mammalian genomes. They play important regulatory roles in various biological processes; However, the function and expression profile of lncRNAs involved in the development of hair follicles in the fetus, have been relatively under-explored area. To investigate the specific role of lncRNAs and mRNAs that regulate hair follicle development, we herein performed a comprehensive study on the lncRNA and mRNA expression profiles of sheep at multiple embryonic days (E65, E85, E105, and E135) and six lambs aged one week (D7) and one month (D30) using RNA-seq technology. The number of genes (471 lncRNAs and 12,812 mRNAs) differentially expressed and potential targets of differentially expressed lncRNAs were predicted. Differentially expressed lncRNAs were grouped into 10 clusters based on their expression pattern by K-means clustering. Moreover, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that some differentially expressed mRNAs, such as DKK1, DSG4, FOXE1, Hoxc13, SFRP1, SFRP2, and Wnt10A overlapped with lncRNAs targets, and enriched in important hair follicle developmental pathways, including Wnt, TNF, and MAPK signaling pathways. In addition, 9 differentially expressed lncRNAs and 4 differentially expressed mRNAs were validated using quantitative real-time PCR (qRT-PCR). This study helps enrich the Ovis lncRNA databases and provides a comprehensive lncRNA transcriptome profile of fetal and postnatal skin of sheep. Additionally, it provides a foundation for further experiments on the role of lncRNAs in the regulation of hair growth in sheep.

Long Noncoding RNA in Preeclampsia: Transcriptional Noise or Innovative Indicators?

Preeclampsia (PE) is termed as an obstetric issue that is characterized by hypertension (≧140/90 mm Hg), together with proteinuria following 20 weeks of pregnancy. Until today, PE still constitutes a severe threat to the lives of both the mothers and fetuses. In the past, long noncoding RNAs (lncRNAs) were considered as the transcriptional noise. However, some investigations have indicated that lncRNAs could be used as innovative indicators in PE. The current review aims to discuss the relationship between lncRNAs and PE in recent years. According to the retrieved data, we concluded that lncRNAs can exert an impact on both the occurrence and development of PE through the changes in the biological functions of trophoblasts, immune regulation, epigenetic regulation, decidualization, and energy metabolism. The mechanisms of lncRNAs in PE will help us to better understand the pathogenesis of PE and help us to find targets for predicting and diagnosing PE in the future.

Analysis of mRNA abundance for histone variants, histone- and DNA-modifiers in bovine in vivo and in vitro oocytes and embryos

Transcript abundance of histone variants, modifiers of histone and DNA in bovine in vivo oocytes and embryos were measured as mean transcripts per million (TPM). Six of 14 annotated histone variants, 8 of 52 histone methyl-transferases, 5 of 29 histone de-methylases, 5 of 20 acetyl-transferases, 5 of 19 de-acetylases, 1 of 4 DNA methyl-transferases and 0 of 3 DNA de-methylases were abundant (TPM >50) in at least one stage studied. Overall, oocytes and embryos contained more varieties of mRNAs for histone modification than for DNA. Three expression patterns were identified for histone modifiers: (1) transcription before embryonic genome activation (EGA) and down-regulated thereafter such as PRMT1; (2) low in oocytes but transiently increased for EGA such as EZH2; (3) high in oocytes but decreased by EGA such as SETD3. These expression patterns were altered by in vitro culture. Additionally, the presence of mRNAs for the TET enzymes throughout pre-implantation development suggests persistent de-methylation. Together, although DNA methylation changes are well-recognized, the first and second orders of significance in epigenetic changes by in vivo embryos may be histone variant replacements and modifications of histones.

Temporal dynamics in meta longitudinal RNA-Seq data

Identification of differentially expressed genes has been a high priority task of downstream analyses to further advances in biomedical research. Investigators have been faced with an array of issues in dealing with more complicated experiments and metadata, including batch effects, normalization, temporal dynamics (temporally differential expression), and isoform diversity (isoform-level quantification and differential splicing events). To date, there are currently no standard approaches to precisely and efficiently analyze these moderate or large-scale experimental designs, especially with combined metadata. In this report, we propose comprehensive analytical pipelines to precisely characterize temporal dynamics in differential expression of genes and other genomic features, i.e., the variability of transcripts, isoforms and exons, by controlling batch effects and other nuisance factors that could have significant confounding effects on the main effects of interest in comparative models and may result in misleading interpretations.

Varying dietary protein and fat elicits differential transcriptomic expression within stress response pathways in preweaned Holstein heifers

Increases in milk replacer dietary energy subsequently increase growth and weight in preweaned dairy heifers. However, the underlying effects of dietary component increases on key functional pathways have yet to be fully investigated. Elucidating these relationships may provide insights into the mechanisms through which protein and fat are partitioned for tissue growth and metabolism. We hypothesized that genes within key growth and metabolic pathways would be differentially expressed between calves fed a protein- and fat-restricted diet and calves fed a protein- and fat-enhanced diet. The objectives of this study were to (1) identify genes differentially expressed between dietary restricted calves and enhanced calves and (2) determine the key regulatory pathways influenced by these genes. Preweaned Holstein heifers (n = 12; 6 ± 0.02 d of age) were randomly assigned to 1 of 2 milk replacer diets: enhanced (28.9% crude protein, 26.2% fat; n = 6) or restricted (20.9% crude protein, 19.8% fat; n = 6). Growth measures included average daily gain and gain-to-feed ratio. After 56 d, calves were killed for tissue collection. Samples from longissimus dorsi, adipose, and liver tissues were collected and RNA was isolated for RNA sequencing analysis. The MIXED procedure of SAS (SAS Institute Inc., Cary, NC) was used to evaluate relationships of growth with dietary energy. Fixed effects included date of collection and time (day). Random effects included sire and birth weight. The RNA sequencing analysis was performed using CLC Genomics Workbench (Qiagen, Germantown, MD), and the Robinson and Smith exact test was used to identify differentially expressed genes between diets. The Protein Analysis Through Evolutionary Relationships (PANTHER) database was then used to identify functional categories of differentially expressed genes. Enhanced calves had increased growth rates and feed efficiency compared with restricted calves (average daily gain = 0.76 and 0.22, respectively; gain-to-feed ratio = 0.10 and 0.06, respectively). There were 238 differentially expressed genes in adipose, 227 in longissimus dorsi, and 40 in liver. We identified 10 genes concordant among tissues. As expected, functional analyses suggested that the majority of genes were associated with metabolic or cellular processes, predominantly cell communication and cell cycle. Overall, it appears that varying levels of dietary protein and fat influence calf growth and development through metabolic processes, including oxidative phosphorylation and glyceroneogenesis. However, protein- and fat-restricted calves appeared to experience metabolic stress at a cellular level, as evidenced by an upregulation in stress response pathways, including genes in the p53 pathway. Calves could be fed at a higher level of protein and fat to decrease the prevalence of metabolic stress at the cellular level, but evidence indicating the presence of inflammatory stress and adipose fibrosis in enhanced calves prompts further investigation of the effects of milk replacer component levels.

Simulation-based benchmarking of isoform quantification in single-cell RNA-seq

Single-cell RNA-seq has the potential to facilitate isoform quantification as the confounding factor of a mixed population of cells is eliminated. However, best practice for using existing quantification methods has not been established. We carry out a benchmark for five popular isoform quantification tools. Performance is generally good for simulated data based on SMARTer and SMART-seq2 data. The reduction in performance compared with bulk RNA-seq is small. An important biological insight comes from our analysis of real data which shows that genes that express two isoforms in bulk RNA-seq predominantly express one or neither isoform in individual cells.