Multi-platform assessment of transcriptome profiling using RNA-seq in the ABRF next-generation sequencing study

Transcription factor 3 is dysregulated in megakaryocytes in myelofibrosis

Transcription factor 3 (TCF3) is a DNA transcription factor that modulates megakaryocyte development. Although abnormal TCF3 expression has been identified in a range of hematological malignancies, to date, it has not been investigated in myelofibrosis (MF). MF is a Philadelphia-negative myeloproliferative neoplasm (MPN) that can arise de novo or progress from essential thrombocythemia [ET] and polycythemia vera [PV] and where dysfunctional megakaryocytes have a role in driving the fibrotic progression. We aimed to examine whether TCF3 is dysregulated in megakaryocytes in MPN, and specifically in MF. We first assessed TCF3 protein expression in megakaryocytes using an immunohistochemical approach analyses and showed that TCF3 was reduced in MF compared with ET and PV. Further, the TCF3-negative megakaryocytes were primarily located near trabecular bone and had the typical "MF-like" morphology as described by the WHO. Genomic analysis of isolated megakaryocytes showed three mutations, all predicted to result in a loss of function, in patients with MF; none were seen in megakaryocytes isolated from ET or PV marrow samples. We then progressed to transcriptomic sequencing of platelets which showed loss of TCF3 in MF. These proteomic, genomic and transcriptomic analyses appear to indicate that TCF3 is downregulated in megakaryocytes in MF. This infers aberrations in megakaryopoiesis occur in this progressive phase of MPN. Further exploration of this pathway could provide insights into TCF3 and the evolution of fibrosis and potentially lead to new preventative therapeutic targets.

Effect of Different Dietary Iron Contents on Liver Transcriptome Characteristics in Wujin Pigs

Iron is an important trace element that affects the growth and development of animals and regulates oxygen transport, hematopoiesis, and hypoxia adaptations. Wujin pig has unique hypoxic adaptability and iron homeostasis; however, the specific regulatory mechanisms have rarely been reported. This study randomly divided 18 healthy Wujin piglets into three groups: the control group, supplemented with 100 mg/kg iron (as iron glycinate); the low-iron group, no iron supplementation; and the high-iron group, supplemented with 200 mg/kg iron (as iron glycinate). The pre-feeding period was 5 days, and the formal period was 30 days. Serum was collected from empty stomachs before slaughter and at slaughter to detect changes in the serum iron metabolism parameters. Gene expression in the liver was analyzed via transcriptome analysis to determine the effects of low- and high-iron diets on transcriptome levels. Correlation analysis was performed for apparent serum parameters, and transcriptome sequencing was performed using weighted gene co-expression network analysis to reveal the key pathways underlying hypoxia regulation and iron metabolism. The main results are as follows. (1) Except for the hypoxia-inducible factor 1 (HIF-1) content (between the low- and high-iron groups), significant differences were not observed among the serum iron metabolic parameters. The serum HIF-1 content of the low-iron group was significantly higher than that of the high-iron group (p < 0.05). (2) Sequencing analysis of the liver transcriptome revealed 155 differentially expressed genes (DEGs) between the low-iron and control groups, 229 DEGs between the high-iron and control groups, and 279 DEGs between the low- and high-iron groups. Bioinformatics analysis showed that the HIF-1 and transforming growth factor-beta (TGF-β) signaling pathways were the key pathways for hypoxia regulation and iron metabolism. Four genes were selected for qPCR validation, and the results were consistent with the transcriptome sequencing data. In summary, the serum iron metabolism parameter results showed that under the influence of low- and high-iron diets, Wujin piglets maintain a steady state of physiological and biochemical indices via complex metabolic regulation of the body, which reflects their stress resistance and adaptability. The transcriptome results revealed the effects of low-iron and high-iron diets on the gene expression level in the liver and showed that the HIF-1 and TGF-β signaling pathways were key for regulating hypoxia adaptability and iron metabolism homeostasis under low-iron and high-iron diets. Moreover, HIF-1α and HEPC were the key genes. The findings provide a theoretical foundation for exploring the regulatory pathways and characteristics of iron metabolism in Wujin pigs.

A real-world multi-center RNA-seq benchmarking study using the Quartet and MAQC reference materials

Translating RNA-seq into clinical diagnostics requires ensuring the reliability and cross-laboratory consistency of detecting clinically relevant subtle differential expressions, such as those between different disease subtypes or stages. As part of the Quartet project, we present an RNA-seq benchmarking study across 45 laboratories using the Quartet and MAQC reference samples spiked with ERCC controls. Based on multiple types of 'ground truth', we systematically assess the real-world RNA-seq performance and investigate the influencing factors involved in 26 experimental processes and 140 bioinformatics pipelines. Here we show greater inter-laboratory variations in detecting subtle differential expressions among the Quartet samples. Experimental factors including mRNA enrichment and strandedness, and each bioinformatics step, emerge as primary sources of variations in gene expression. We underscore the profound influence of experimental execution, and provide best practice recommendations for experimental designs, strategies for filtering low-expression genes, and the optimal gene annotation and analysis pipelines. In summary, this study lays the foundation for developing and quality control of RNA-seq for clinical diagnostic purposes.

Quartet RNA reference materials improve the quality of transcriptomic data through ratio-based profiling

Certified RNA reference materials are indispensable for assessing the reliability of RNA sequencing to detect intrinsically small biological differences in clinical settings, such as molecular subtyping of diseases. As part of the Quartet Project for quality control and data integration of multi-omics profiling, we established four RNA reference materials derived from immortalized B-lymphoblastoid cell lines from four members of a monozygotic twin family. Additionally, we constructed ratio-based transcriptome-wide reference datasets between two samples, providing cross-platform and cross-laboratory 'ground truth'. Investigation of the intrinsically subtle biological differences among the Quartet samples enables sensitive assessment of cross-batch integration of transcriptomic measurements at the ratio level. The Quartet RNA reference materials, combined with the ratio-based reference datasets, can serve as unique resources for assessing and improving the quality of transcriptomic data in clinical and biological settings.

Reliable biological and multi-omics research through biometrology

Metrology is the science of measurement and its applications, whereas biometrology is the science of biological measurement and its applications. Biometrology aims to achieve accuracy and consistency of biological measurements by focusing on the development of metrological traceability, biological reference measurement procedures, and reference materials. Irreproducibility of biological and multi-omics research results from different laboratories, platforms, and analysis methods is hampering the translation of research into clinical uses and can often be attributed to the lack of biologists' attention to the general principles of metrology. In this paper, the progresses of biometrology including metrology on nucleic acid, protein, and cell measurements and its impacts on the improvement of reliability and comparability in biological research are reviewed. Challenges in obtaining more reliable biological and multi-omics measurements due to the lack of primary reference measurement procedures and new standards for biological reference materials faced by biometrology are discussed. In the future, in addition to establishing reliable reference measurement procedures, developing reference materials from single or multiple parameters to multi-omics scale should be emphasized. Thinking in way of biometrology is warranted for facilitating the translation of high-throughput omics research into clinical practices.

Precision and Accuracy of Single-Cell/Nuclei RNA Sequencing Data

Single-cell/nuclei RNA sequencing (sc/snRNA-Seq) is widely used for profiling cell-type gene expressions in biomedical research. An important but underappreciated issue is the quality of sc/snRNA-Seq data that would impact the reliability of downstream analyses. Here we evaluated the precision and accuracy in 18 sc/snRNA-Seq datasets. The precision was assessed on data from human brain studies with a total of 3,483,905 cells from 297 individuals, by utilizing technical replicates. The accuracy was evaluated with sample-matched scRNA-Seq and pooled-cell RNA-Seq data of cultured mononuclear phagocytes from four species. The results revealed low precision and accuracy at the single-cell level across all evaluated data. Cell number and RNA quality were highlighted as two key factors determining the expression precision, accuracy, and reproducibility of differential expression analysis in sc/snRNA-Seq. This study underscores the necessity of sequencing enough high-quality cells per cell type per individual, preferably in the hundreds, to mitigate noise in expression quantification.

Plasma-Free Blood as a Potential Alternative to Whole Blood for Transcriptomic Analysis

RNA sequencing (RNAseq) technology has become increasingly important in precision medicine and clinical diagnostics, and emerged as a powerful tool for identifying protein-coding genes, performing differential gene analysis, and inferring immune cell composition. Human peripheral blood samples are widely used for RNAseq, providing valuable insights into individual biomolecular information. Blood samples can be classified as whole blood (WB), plasma, serum, and remaining sediment samples, including plasma-free blood (PFB) and serum-free blood (SFB) samples that are generally considered less useful byproducts during the processes of plasma and serum separation, respectively. However, the feasibility of using PFB and SFB samples for transcriptome analysis remains unclear. In this study, we aimed to assess the suitability of employing PFB or SFB samples as an alternative RNA source in transcriptomic analysis. We performed a comparative analysis of WB, PFB, and SFB samples for different applications. Our results revealed that PFB samples exhibit greater similarity to WB samples than SFB samples in terms of protein-coding gene expression patterns, detection of differentially expressed genes, and immunological characterizations, suggesting that PFB can serve as a viable alternative to WB for transcriptomic analysis. Our study contributes to the optimization of blood sample utilization and the advancement of precision medicine research.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00121-1.

Single cell RNA-seq: a novel tool to unravel virus-host interplay

Over the last decade, single cell RNA sequencing (scRNA-seq) technology has caught the momentum of being a vital revolutionary tool to unfold cellular heterogeneity by high resolution assessment. It evades the inadequacies of conventional sequencing technology which was able to detect only average expression level among cell populations. In the era of twenty-first century, several epidemic and pandemic viruses have emerged. Being an intracellular entity, viruses totally rely on host. Complex virus-host dynamics result when the virus tend to obtain factors from host cell required for its replication and establishment of infection. As a prevailing tool, scRNA-seq is able to understand virus-host interplay by comprehensive transcriptome profiling. Because of technological and methodological advancement, this technology is capable to recognize viral genome and host cell response heterogeneity. Further development in analytical methods with multiomics approach and increased availability of accessible scRNA-seq datasets will improve the understanding of viral pathogenesis that can be helpful for development of novel antiviral therapeutic strategies.

Standardizing digital biobanks: integrating imaging, genomic, and clinical data for precision medicine

Advancements in data acquisition and computational methods are generating a large amount of heterogeneous biomedical data from diagnostic domains such as clinical imaging, pathology, and next-generation sequencing (NGS), which help characterize individual differences in patients. However, this information needs to be available and suitable to promote and support scientific research and technological development, supporting the effective adoption of the precision medicine approach in clinical practice. Digital biobanks can catalyze this process, facilitating the sharing of curated and standardized imaging data, clinical, pathological and molecular data, crucial to enable the development of a comprehensive and personalized data-driven diagnostic approach in disease management and fostering the development of computational predictive models. This work aims to frame this perspective, first by evaluating the state of standardization of individual diagnostic domains and then by identifying challenges and proposing a possible solution towards an integrative approach that can guarantee the suitability of information that can be shared through a digital biobank. Our analysis of the state of the art shows the presence and use of reference standards in biobanks and, generally, digital repositories for each specific domain. Despite this, standardization to guarantee the integration and reproducibility of the numerical descriptors generated by each domain, e.g. radiomic, pathomic and -omic features, is still an open challenge. Based on specific use cases and scenarios, an integration model, based on the JSON format, is proposed that can help address this problem. Ultimately, this work shows how, with specific standardization and promotion efforts, the digital biobank model can become an enabling technology for the comprehensive study of diseases and the effective development of data-driven technologies at the service of precision medicine.

Transcriptome sequencing and screening of genes related to the MADS-box gene family in Clematis courtoisii

The MADS-box gene family controls plant flowering and floral organ development; therefore, it is particularly important in ornamental plants. To investigate the genes associated with the MADS-box family in Clematis courtoisii, we performed full-length transcriptome sequencing on C. courtoisii using the PacBio Sequel third-generation sequencing platform, as no reference genome data was available. A total of 12.38 Gb of data, containing 9,476,585 subreads and 50,439 Unigenes were obtained. According to functional annotation, a total of 37,923 Unigenes (75.18% of the total) were assigned with functional annotations, and 50 Unigenes were identified as MADS-box related genes. Subsequently, we employed hmmerscan to perform protein sequence similarity search for the translated Unigene sequences and successfully identified 19 Unigenes associated with the MADS-box gene family, including MIKC*(1) and MIKCC (18) genes. Furthermore, within the MIKCC group, six subclasses can be further distinguished.

A uniform data processing pipeline enables harmonized nanoparticle protein corona analysis across proteomics core facilities

Protein corona, a layer of biomolecules primarily comprising proteins, forms dynamically on nanoparticles in biological fluids and is crucial for predicting nanomedicine safety and efficacy. The protein composition of the corona layer is typically analyzed using liquid chromatography-mass spectrometry (LC-MS/MS). Our recent study, involving identical samples analyzed by 17 proteomics facilities, highlighted significant data variability, with only 1.8% of proteins consistently identified across these centers. Here, we implement an aggregated database search unifying parameters such as variable modifications, enzyme specificity, number of allowed missed cleavages and a stringent 1% false discovery rate at the protein and peptide levels. Such uniform search dramatically harmonizes the proteomics data, increasing the reproducibility and the percentage of consistency-identified unique proteins across distinct cores. Specifically, out of the 717 quantified proteins, 253 (35.3%) are shared among the top 5 facilities (and 16.2% among top 11 facilities). Furthermore, we note that reduction and alkylation are important steps in protein corona sample processing and as expected, omitting these steps reduces the number of total quantified peptides by around 20%. These findings underscore the need for standardized procedures in protein corona analysis, which is vital for advancing clinical applications of nanoscale biotechnologies.

miRNA and mRNA Signatures in Human Acute Kidney Injury Tissue

Acute kidney injury (AKI) is an important contributor to the development of chronic kidney disease (CKD). There is a need to understand molecular mediators that drive either recovery or progression to CKD. In particular, the role of miRNA and its regulatory role in AKI is poorly understood. We performed miRNA and mRNA sequencing on biobanked human kidney tissues obtained in the routine clinical care of patients with the diagnoses of AKI and minimal change disease (MCD), in addition to nephrectomized (Ref) tissue from individuals without known kidney disease. Transcriptomic analysis of mRNA revealed that Ref tissues exhibited a similar injury signature to AKI, not identified in MCD samples. The transcriptomic signature of human AKI was enriched with genes in pathways involved in cell adhesion and epithelial-to-mesenchymal transition (e.g., CDH6, ITGB6, CDKN1A ). miRNA DE analysis revealed upregulation of miRNA associated with immune cell recruitment and inflammation (e.g., miR-146a, miR-155, miR-142, miR-122). These miRNA (i.e., miR-122, miR-146) are also associated with downregulation of mRNA such as DDR2 and IGFBP6 , respectively. These findings suggest integrated interactions between miRNAs and target mRNAs in AKI-related processes such as inflammation, immune cell activation and epithelial-to-mesenchymal transition. These data contribute several novel findings when describing the epigenetic regulation of AKI by miRNA, and also underscores the importance of utilizing an appropriate reference control tissue to understand canonical pathway alterations in AKI.

Transcriptomics for Clinical and Experimental Biology Research: Hang on a Seq

Sequencing the human genome empowers translational medicine, facilitating transcriptome-wide molecular diagnosis, pathway biology, and drug repositioning. Initially, microarrays are used to study the bulk transcriptome; but now short-read RNA sequencing (RNA-seq) predominates. Positioned as a superior technology, that makes the discovery of novel transcripts routine, most RNA-seq analyses are in fact modeled on the known transcriptome. Limitations of the RNA-seq methodology have emerged, while the design of, and the analysis strategies applied to, arrays have matured. An equitable comparison between these technologies is provided, highlighting advantages that modern arrays hold over RNA-seq. Array protocols more accurately quantify constitutively expressed protein coding genes across tissue replicates, and are more reliable for studying lower expressed genes. Arrays reveal long noncoding RNAs (lncRNA) are neither sparsely nor lower expressed than protein coding genes. Heterogeneous coverage of constitutively expressed genes observed with RNA-seq, undermines the validity and reproducibility of pathway analyses. The factors driving these observations, many of which are relevant to long-read or single-cell sequencing are discussed. As proposed herein, a reappreciation of bulk transcriptomic methods is required, including wider use of the modern high-density array data-to urgently revise existing anatomical RNA reference atlases and assist with more accurate study of lncRNAs.

WheatCENet: A Database for Comparative Co-expression Networks Analysis of Allohexaploid Wheat and Its Progenitors

Genetic and epigenetic changes after polyploidization events could result in variable gene expression and modified regulatory networks. Here, using large-scale transcriptome data, we constructed co-expression networks for diploid, tetraploid, and hexaploid wheat species, and built a platform for comparing co-expression networks of allohexaploid wheat and its progenitors, named WheatCENet. WheatCENet is a platform for searching and comparing specific functional co-expression networks, as well as identifying the related functions of the genes clustered therein. Functional annotations like pathways, gene families, protein-protein interactions, microRNAs (miRNAs), and several lines of epigenome data are integrated into this platform, and Gene Ontology (GO) annotation, gene set enrichment analysis (GSEA), motif identification, and other useful tools are also included. Using WheatCENet, we found that the network of WHEAT ABERRANT PANICLE ORGANIZATION 1 (WAPO1) has more co-expressed genes related to spike development in hexaploid wheat than its progenitors. We also found a novel motif of CCWWWWWWGG (CArG) specifically in the promoter region of WAPO-A1, suggesting that neofunctionalization of the WAPO-A1 gene affects spikelet development in hexaploid wheat. WheatCENet is useful for investigating co-expression networks and conducting other analyses, and thus facilitates comparative and functional genomic studies in wheat. WheatCENet is freely available at http://bioinformatics.cpolar.cn/WheatCENet and http://bioinformatics.cau.edu.cn/WheatCENet.

A comprehensive genomic and transcriptomic dataset of triple-negative breast cancers

Molecular subtyping of triple-negative breast cancer (TNBC) is essential for understanding the mechanisms and discovering actionable targets of this highly heterogeneous type of breast cancer. We previously performed a large single-center and multiomics study consisting of genomics, transcriptomics, and clinical information from 465 patients with primary TNBC. To facilitate reusing this unique dataset, we provided a detailed description of the dataset with special attention to data quality in this study. The multiomics data were generally of high quality, but a few sequencing data had quality issues and should be noted in subsequent data reuse. Furthermore, we reconduct data analyses with updated pipelines and the updated version of the human reference genome from hg19 to hg38. The updated profiles were in good concordance with those previously published in terms of gene quantification, variant calling, and copy number alteration. Additionally, we developed a user-friendly web-based database for convenient access and interactive exploration of the dataset. Our work will facilitate reusing the dataset, maximize the values of data and further accelerate cancer research.

Measurement(s)	RNA expression profiling • whole-exome sequencing (WES) • somatic mutations • copy number alterations (CNAs)
Technology Type(s)	RNA sequencing • DNA sequencing • OncoScan CNV assay
Factor Type(s)	Intervention or procedure
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Location	China

The role of antigen expression in shaping the repertoire of HLA presented ligands

Human leukocyte antigen (HLA) presentation of peptides is a prerequisite of T cell immune activation. The understanding of the rules defining this event has large implications for our knowledge of basic immunology and for the rational design of immuno-therapeutics and vaccines. Historically, most of the available prediction methods have been solely focused on the information related to antigen processing and presentation. Recent work has, however, demonstrated that method performance can be boosted by integrating information related to antigen abundance. Here we expand on these later findings and develop an extended version of NetMHCpan, called NetMHCpanExp, integrating information on antigen abundance from RNA-Seq experiments. In line with earlier works, the model demonstrates improved performance for both HLA ligand and cancer neoantigen epitope prediction. Optimal results are obtained by use of sample-specific abundance information but also reference datasets can be applied with a limited performance drop. The developed tool is available at https://services.healthtech.dtu.dk/service.php?NetMHCpanExp-1.0.

Advancements in Genomic and Behavioral Neuroscience Analysis for the Study of Normal and Pathological Brain Function

Psychiatric and neurological disorders are influenced by an undetermined number of genes and molecular pathways that may differ among afflicted individuals. Functionally testing and characterizing biological systems is essential to discovering the interrelationship among candidate genes and understanding the neurobiology of behavior. Recent advancements in genetic, genomic, and behavioral approaches are revolutionizing modern neuroscience. Although these tools are often used separately for independent experiments, combining these areas of research will provide a viable avenue for multidimensional studies on the brain. Herein we will briefly review some of the available tools that have been developed for characterizing novel cellular and animal models of human disease. A major challenge will be openly sharing resources and datasets to effectively integrate seemingly disparate types of information and how these systems impact human disorders. However, as these emerging technologies continue to be developed and adopted by the scientific community, they will bring about unprecedented opportunities in our understanding of molecular neuroscience and behavior.

Whole Transcriptome Profiling of the Effects of Cadmium on the Liver of the Xiangxi Yellow Heifer

Cadmium (Cd) is a major heavy metal toxicant found in industrial zones. Humans and animals are exposed to it through their diet, which results in various physiological problems. In the current study, the toxic effects of Cd on the liver were investigated by whole-transcriptome sequencing (RNA-seq) of the livers of Xiangxi heifers fed a diet with excess Cd. We randomly divided six healthy heifers into two groups. The first group received a control diet, whereas the second group received Cd-exceeding diets for 100 days. After 100 days, the livers were collected. A total of 551 differentially expressed mRNAs, 24 differentially expressed miRNAs, and 169 differentially expressed lncRNAs were identified (p < 0.05, |log2FC| >1). Differentially expressed genes (DEGs) were analyzed by gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses. We found that under Cd exposure, DEGs were enriched in the adenosine 5'-monophosphate-activated protein kinase pathway, which is involved in autophagy regulation, and the peroxisome proliferator-activated receptor pathway, which is involved in lipid metabolism. In addition, the apolipoprotein A4 gene, which has anti-inflammatory and antioxidant effects, the anti-apoptotic gene ATPase H+/K+ transporting the nongastric alpha2 subunit, and the cholesterol metabolism-associated gene endothelial lipase gene were significantly downregulated. C-X-C motif chemokine ligand 3, cholesterol 7α-hydroxylase, and stearoyl-CoA desaturase, which are involved in the development of fatty liver, were significantly upregulated. These genes revealed the main effects of Cd on the liver of Xiangxi yellow heifers. The current study provides insightful information regarding the DEGs involved in autophagy regulation, apoptosis, lipid metabolism, anti-inflammation, and antioxidant enzyme activity. These may serve as useful biomarkers for predicting and treating Cd-related diseases in the future.

Sequencing of individual barcoded cDNAs using Pacific Biosciences and Oxford Nanopore Technologies reveals platform-specific error patterns

Long-read transcriptomics require understanding error sources inherent to technologies. Current approaches cannot compare methods for an individual RNA molecule. Here, we present a novel platform-comparison method that combines barcoding strategies and long-read sequencing to sequence cDNA copies representing an individual RNA molecule on both Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT). We compare these long-read pairs in terms of sequence content and isoform patterns. Although individual read pairs show high similarity, we find differences in (1) aligned length, (2) transcription start site (TSS), (3) polyadenylation site (poly(A)-site) assignment, and (4) exon-intron structures. Overall, 25% of read pairs disagree on either TSS, poly(A)-site, or splice site. Intron-chain disagreement typically arises from alignment errors of microexons and complicated splice sites. Our single-molecule technology comparison reveals that inconsistencies are often caused by sequencing error-induced inaccurate ONT alignments, especially to downstream GUNNGU donor motifs. However, annotation-disagreeing upstream shifts in NAGNAG acceptors in ONT are often confirmed by PacBio and are thus likely real. In both barcoded and nonbarcoded ONT reads, we find that intron number and proximity of GU/AGs better predict inconsistencies with the annotation than read quality alone. We summarize these findings in an annotation-based algorithm for spliced alignment correction that improves subsequent transcript construction with ONT reads.

Kinase-dead mutation: A novel strategy for improving soybean resistance to soybean cyst nematode Heterodera glycines

Protein kinases phosphorylate proteins for functional changes and are involved in nearly all cellular processes, thereby regulating almost all aspects of plant growth and development, and responses to biotic and abiotic stresses. We generated two independent co-expression networks of soybean genes using control and stress response gene expression data and identified 392 differentially highly interconnected kinase hub genes among the two networks. Of these 392 kinases, 90 genes were identified as "syncytium highly connected hubs", potentially essential for activating kinase signalling pathways in the nematode feeding site. Overexpression of wild-type coding sequences of five syncytium highly connected kinase hub genes using transgenic soybean hairy roots enhanced plant susceptibility to soybean cyst nematode (SCN; Heterodera glycines) Hg Type 0 (race 3). In contrast, overexpression of kinase-dead variants of these five syncytium kinase hub genes significantly enhanced soybean resistance to SCN. Additionally, three of the five tested kinase hub genes enhanced soybean resistance to SCN Hg Type 1.2.5.7 (race 2), highlighting the potential of the kinase-dead approach to generate effective and durable resistance against a wide range of SCN Hg types. Subcellular localization analysis revealed that kinase-dead mutations do not alter protein cellular localization, confirming the structure-function of the kinase-inactive variants in producing loss-of-function phenotypes causing significant decrease in nematode susceptibility. Because many protein kinases are highly conserved and are involved in plant responses to various biotic and abiotic stresses, our approach of identifying kinase hub genes and their inactivation using kinase-dead mutation could be translated for biotic and abiotic stress tolerance.

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.

Deciphering the Genome-Wide Transcriptomic Changes during Interactions of Resistant and Susceptible Genotypes of American Elm with Ophiostoma novo-ulmi

Dutch elm disease (DED), caused by Ophiostoma novo-ulmi (Onu), is a destructive disease of American elm (Ulmus americana L.). The molecular mechanisms of resistance and susceptibility against DED in American elm are still largely uncharacterized. In the present study, we performed a de novo transcriptome (RNA-sequencing; RNA-Seq) assembly of U. americana and compared the gene expression in a resistant genotype, 'Valley Forge', and a susceptible (S) elm genotype at 0 and 96 h post-inoculation of Onu. A total of 85,863 non-redundant unigenes were identified. Compared to the previously characterized U. minor transcriptome, U. americana has 35,290 similar and 55,499 unique genes. The transcriptomic variations between 'Valley Forge' and 'S' were found primarily in the photosynthesis and primary metabolism, which were highly upregulated in the susceptible genotype irrespective of the Onu inoculation. The resistance to DED was associated with the activation of RPM1-mediated effector-triggered immunity that was demonstrated by the upregulation of genes involved in the phenylpropanoids biosynthesis and PR genes. The most significantly enriched gene ontology (GO) terms in response to Onu were response to stimulus (GO:0006950), response to stress (GO:0050896), and secondary metabolic process (GO:0008152) in both genotypes. However, only in the resistant genotype, the defense response (GO:0006952) was among the topmost significantly enriched GO terms. Our findings revealed the molecular regulations of DED resistance and susceptibility and provide a platform for marker-assisted breeding of resistant American elm genotypes.

A method for simultaneous detection of small and long RNA biotypes by ribodepleted RNA-Seq

RNA sequencing offers unprecedented access to the transcriptome. Key to this is the identification and quantification of many different species of RNA from the same sample at the same time. In this study we describe a novel protocol for simultaneous detection of coding and non-coding transcripts using modifications to the Ion Total RNA-Seq kit v2 protocol, with integration of QIASeq FastSelect rRNA removal kit. We report highly consistent sequencing libraries can be produced from both frozen high integrity mouse hippocampal tissue and the more challenging post-mortem human tissue. Removal of rRNA using FastSelect was extremely efficient, resulting in less than 1.5% rRNA content in the final library. We identified > 30,000 unique transcripts from all samples, including protein-coding genes and many species of non-coding RNA, in biologically-relevant proportions. Furthermore, the normalized sequencing read count for select genes significantly negatively correlated with Ct values from qRT-PCR analysis from the same samples. These results indicate that this protocol accurately and consistently identifies and quantifies a wide variety of transcripts simultaneously. The highly efficient rRNA depletion, coupled with minimized sample handling and without complicated and high-loss size selection protocols, makes this protocol useful to researchers wishing to investigate whole transcriptomes.

The Vacc-SeqQC project: Benchmarking RNA-Seq for clinical vaccine studies

Introduction

Over the last decade, the field of systems vaccinology has emerged, in which high throughput transcriptomics and other omics assays are used to probe changes of the innate and adaptive immune system in response to vaccination. The goal of this study was to benchmark key technical and analytical parameters of RNA sequencing (RNA-seq) in the context of a multi-site, double-blind randomized vaccine clinical trial.

Methods

We collected longitudinal peripheral blood mononuclear cell (PBMC) samples from 10 subjects before and after vaccination with a live attenuated Francisella tularensis vaccine and performed RNA-Seq at two different sites using aliquots from the same sample to generate two replicate datasets (5 time points for 50 samples each). We evaluated the impact of (i) filtering lowly-expressed genes, (ii) using external RNA controls, (iii) fold change and false discovery rate (FDR) filtering, (iv) read length, and (v) sequencing depth on differential expressed genes (DEGs) concordance between replicate datasets. Using synthetic mRNA spike-ins, we developed a method for empirically establishing minimal read-count thresholds for maintaining fold change accuracy on a per-experiment basis. We defined a reference PBMC transcriptome by pooling sequence data and established the impact of sequencing depth and gene filtering on transcriptome representation. Lastly, we modeled statistical power to detect DEGs for a range of sample sizes, effect sizes, and sequencing depths.

Results and Discussion

Our results showed that (i) filtering lowly-expressed genes is recommended to improve fold-change accuracy and inter-site agreement, if possible guided by mRNA spike-ins (ii) read length did not have a major impact on DEG detection, (iii) applying fold-change cutoffs for DEG detection reduced inter-set agreement and should be used with caution, if at all, (iv) reduction in sequencing depth had a minimal impact on statistical power but reduced the identifiable fraction of the PBMC transcriptome, (v) after sample size, effect size (i.e. the magnitude of fold change) was the most important driver of statistical power to detect DEG. The results from this study provide RNA sequencing benchmarks and guidelines for planning future similar vaccine studies.

RNA isolation from micro-quantity of articular cartilage for quantitative gene expression by microarray analysis

Isolation of quality RNA from articular cartilage has been challenging due to low cellularity and the high abundance of extracellular matrix and proteoglycan proteins. Recently developed methods for isolation of high quality RNA from cartilage are more applicable to larger cartilage specimens typically weighing at least 25 mg. While these methods generate RNA suitable for analysis, they are less successful with smaller tissue inputs. For the study of small focal defect cartilage specimens an improved RNA extraction method is needed. Here we report a protocol for direct RNA isolation from less than 3 mg of wet weight rabbit articular cartilage for quantitative microarray gene profiling. This protocol is useful for identifying differentially expressed genes in chondrocytes following focal cartilage repair and can potentially be adopted for gene expression analysis of cartilage biopsy specimens from human joints.

Flimma: a federated and privacy-aware tool for differential gene expression analysis

Aggregating transcriptomics data across hospitals can increase sensitivity and robustness of differential expression analyses, yielding deeper clinical insights. As data exchange is often restricted by privacy legislation, meta-analyses are frequently employed to pool local results. However, the accuracy might drop if class labels are inhomogeneously distributed among cohorts. Flimma ( https://exbio.wzw.tum.de/flimma/ ) addresses this issue by implementing the state-of-the-art workflow limma voom in a federated manner, i.e., patient data never leaves its source site. Flimma results are identical to those generated by limma voom on aggregated datasets even in imbalanced scenarios where meta-analysis approaches fail.

Revealing protein-protein interactions at the transcriptome scale by sequencing

We describe PROPER-seq (protein-protein interaction sequencing) to map protein-protein interactions (PPIs) en masse. PROPER-seq first converts transcriptomes of input cells into RNA-barcoded protein libraries, in which all interacting protein pairs are captured through nucleotide barcode ligation, recorded as chimeric DNA sequences, and decoded at once by sequencing and mapping. We applied PROPER-seq to human embryonic kidney cells, T lymphocytes, and endothelial cells and identified 210,518 human PPIs (collected in the PROPER v.1.0 database). Among these, 1,365 and 2,480 PPIs are supported by published co-immunoprecipitation (coIP) and affinity purification-mass spectrometry (AP-MS) data, 17,638 PPIs are predicted by the prePPI algorithm without previous experimental validation, and 100 PPIs overlap human synthetic lethal gene pairs. In addition, four previously uncharacterized interaction partners with poly(ADP-ribose) polymerase 1 (PARP1) (a critical protein in DNA repair) known as XPO1, MATR3, IPO5, and LEO1 are validated in vivo. PROPER-seq presents a time-effective technology to map PPIs at the transcriptome scale, and PROPER v.1.0 provides a rich resource for studying PPIs.

Performance assessment of DNA sequencing platforms in the ABRF Next-Generation Sequencing Study

Assessing the reproducibility, accuracy and utility of massively parallel DNA sequencing platforms remains an ongoing challenge. Here the Association of Biomolecular Resource Facilities (ABRF) Next-Generation Sequencing Study benchmarks the performance of a set of sequencing instruments (HiSeq/NovaSeq/paired-end 2 × 250-bp chemistry, Ion S5/Proton, PacBio circular consensus sequencing (CCS), Oxford Nanopore Technologies PromethION/MinION, BGISEQ-500/MGISEQ-2000 and GS111) on human and bacterial reference DNA samples. Among short-read instruments, HiSeq 4000 and X10 provided the most consistent, highest genome coverage, while BGI/MGISEQ provided the lowest sequencing error rates. The long-read instrument PacBio CCS had the highest reference-based mapping rate and lowest non-mapping rate. The two long-read platforms PacBio CCS and PromethION/MinION showed the best sequence mapping in repeat-rich areas and across homopolymers. NovaSeq 6000 using 2 × 250-bp read chemistry was the most robust instrument for capturing known insertion/deletion events. This study serves as a benchmark for current genomics technologies, as well as a resource to inform experimental design and next-generation sequencing variant calling.

Loop-Mediated Isothermal Amplification Detection of SARS-CoV-2 and Myriad Other Applications

As the second year of the COVID-19 pandemic begins, it remains clear that a massive increase in the ability to test for SARS-CoV-2 infections in a myriad of settings is critical to controlling the pandemic and to preparing for future outbreaks. The current gold standard for molecular diagnostics is the polymerase chain reaction (PCR), but the extraordinary and unmet demand for testing in a variety of environments means that both complementary and supplementary testing solutions are still needed. This review highlights the role that loop-mediated isothermal amplification (LAMP) has had in filling this global testing need, providing a faster and easier means of testing, and what it can do for future applications, pathogens, and the preparation for future outbreaks. This review describes the current state of the art for research of LAMP-based SARS-CoV-2 testing, as well as its implications for other pathogens and testing. The authors represent the global LAMP (gLAMP) Consortium, an international research collective, which has regularly met to share their experiences on LAMP deployment and best practices; sections are devoted to all aspects of LAMP testing, including preanalytic sample processing, target amplification, and amplicon detection, then the hardware and software required for deployment are discussed, and finally, a summary of the current regulatory landscape is provided. Included as well are a series of first-person accounts of LAMP method development and deployment. The final discussion section provides the reader with a distillation of the most validated testing methods and their paths to implementation. This review also aims to provide practical information and insight for a range of audiences: for a research audience, to help accelerate research through sharing of best practices; for an implementation audience, to help get testing up and running quickly; and for a public health, clinical, and policy audience, to help convey the breadth of the effect that LAMP methods have to offer.

Novel Facets of the Liver Transcriptome Are Associated with the Susceptibility and Resistance to Lipid-Related Metabolic Disorders in Periparturient Holstein Cows

Simple Summary

Energy and nutrient demands of the early lactation period can result in the development of metabolic disorders, such as ketosis and fatty liver, in dairy cows. Variability in the incidence of these disorders suggests that some cows have an ability to adapt. The objective of this study was to discover differences in liver gene expression that are associated with a cow’s susceptibility (disposition to disorder during typical conditions) or resistance (disposition to disorder onset and severity when presented a challenge) to metabolic disorders. Cows in a control treatment and a ketosis induction protocol treatment were retrospectively grouped into susceptibility and resistance groups, respectively, by a machine learning algorithm using lipid biomarker concentrations. Whole-transcriptome RNA sequencing was performed on liver samples from these cows. More susceptible cows had lower expression of glutathione metabolism genes, while less resistant cows had greater expression of eicosanoid-metabolism-related genes. Additionally, differentially expressed genes suggested a role for immune-response-related genes in conferring susceptibility and resistance to metabolic disorders. The overall inferred metabolism suggests that responses to oxidative stress may determine susceptibility and resistance to metabolic disorders, with novel implications for immunometabolism.

Abstract

Lipid-related metabolic disorders (LRMD) are prevalent in early lactation dairy cows, and have detrimental effects on productivity and health. Our objectives were to identify cows resistant or susceptible to LRMD using a ketosis induction protocol (KIP) to discover differentially expressed liver genes and metabolic pathways associated with disposition. Clustering cows based on postpartum lipid metabolite concentrations within dietary treatments identified cows more or less susceptible (MS vs. LS) to LRMD within the control treatment, and more or less resistant (MR vs. LR) within the KIP treatment. Whole-transcriptome RNA sequencing was performed on liver samples (−28, +1, and +14 days relative to calving) to assess differential gene and pathway expression (LS vs. MS; MR vs. LR; n = 3 cows per cluster). Cows within the MS and LR clusters had evidence of greater blood serum β-hydroxybutyrate concentration and liver triglyceride content than the LS and MR clusters, respectively. The inferred metabolism of differentially expressed genes suggested a role of immune response (i.e., interferon-inducible proteins and major histocompatibility complex molecules). Additionally, unique roles for glutathione metabolism and eicosanoid metabolism in modulating susceptibility and resistance, respectively, were implicated. Overall, this research provides novel insight into the role of immunometabolism in LRMD pathology, and suggests the potential for unique control points for LRMD progression and severity.

Development of omics biomarkers for estrogen exposure using mRNA, miRNA and piRNAs

The number of chemicals requiring risk evaluation exceeds our capacity to provide the underlying data using traditional methodology. This has led to an increased focus on the development of novel approach methodologies. This work aimed to expand the panel of gene expression-based biomarkers to include responses to estrogens, to identify training strategies that maximize the range of applicable concentrations, and to evaluate the potential for two classes of small non-coding RNAs (sncRNAs), microRNA (miRNA) and piwi-interacting RNA (piRNA), as biomarkers. To this end larval Pimephales promelas (96 hpf +/- 1h) were exposed to 5 concentrations of 17α- ethinylestradiol (0.12, 1.25, 2.5, 5.0, 10.0 ng/L) for 48 h. For mRNA-based biomarker development, RNA-seq was conducted across all concentrations. For sncRNA biomarkers, small RNA libraries were prepared only for the control and 10.0 ng/L EE2 treatment. In order to develop an mRNA classifier that remained accurate over the range of exposure concentrations, three different training strategies were employed that focused on 10 ng/L, 2.5 ng/L or a combination of both. Classifiers were tested against an independent test set of individuals exposed to the same concentrations used in training and subsequently against concentrations not included in model training. Both random forest (RF) and logistic regression with elastic net regularizations (glmnet) models trained on 10 ng/L EE2 performed poorly when applied to lower concentrations. RF models trained with either the 2.5 ng/L or combination (2.5 + 10 ng/L) treatments were able to accurately discriminate exposed vs. non-exposed across all but the lowest concentrations. glmnet models were unable to accurately classify below 5 ng/L. With the exception of the 10 ng/L treatment, few mRNA differentially expressed genes (DEG) were observed, however, there was marked overlap of DEGs across treatments. Overlapping DEGs have well established linkages to estrogen and several of the 81 DEGs identified in the 10 ng/L treatment have been previously utilized as estrogenic biomarkers (vitellogenin, estrogen receptor-β). Following multiple test correction, no sncRNAs were found to be differentially expressed, however, both miRNA and piRNA classifiers were able to accurately discriminate control and 10 ng/L exposed organisms with AUCs of 0.83 and 1.0 respectively. We have developed a highly discriminative estrogenic mRNA biomarker that is accurate over a range of concentrations likely to be found in real-world exposures. We have demonstrated that both miRNA and piRNA are responsive to estrogenic exposure, suggesting the need to further investigate their regulatory roles in the estrogenic response.

Constrained Standardization of Count Data from Massive Parallel Sequencing

In high-throughput omics disciplines like transcriptomics, researchers face a need to assess the quality of an experiment prior to an in-depth statistical analysis. To efficiently analyze such voluminous collections of data, researchers need triage methods that are both quick and easy to use. Such a normalization method for relative quantitation, CONSTANd, was recently introduced for isobarically-labeled mass spectra in proteomics. It transforms the data matrix of abundances through an iterative, convergent process enforcing three constraints: (I) identical column sums; (II) each row sum is fixed (across matrices) and (III) identical to all other row sums. In this study, we investigate whether CONSTANd is suitable for count data from massively parallel sequencing, by qualitatively comparing its results to those of DESeq2. Further, we propose an adjustment of the method so that it may be applied to identically balanced but differently sized experiments for joint analysis. We find that CONSTANd can process large data sets at well over 1 million count records per second whilst mitigating unwanted systematic bias and thus quickly uncovering the underlying biological structure when combined with a PCA plot or hierarchical clustering. Moreover, it allows joint analysis of data sets obtained from different batches, with different protocols and from different labs but without exploiting information from the experimental setup other than the delineation of samples into identically processed sets (IPSs). CONSTANd's simplicity and applicability to proteomics as well as transcriptomics data make it an interesting candidate for integration in multi-omics workflows.

Gene expression analysis of combined RNA-seq experiments using a receiver operating characteristic calibrated procedure

Because of rapid advancements in sequencing technology, the experimental platforms of RNA-seq are updated frequently. It is quite common to combine data sets from several experimental platforms for analysis in order to increase the sample size and achieve more powerful tests for detecting the presence of differential gene expression. The data sets combined from different experimental platforms will have a complex data distribution, which causes a major problem in statistical modeling as well as in multiple testing. Although plenty of research have studied this problem by modeling the batch effects, there are no general and robust data-driven procedures for RNA-seq analysis. In this paper we propose a new robust procedure which combines the use of popular methods (packages) with a data-driven simulation (DDS). We construct the average receiver operating characteristic curve through the DDS to provide the calibrated levels of significance for multiple testing. Instead of further modifying the adjusted p-values, we calibrated the levels of significance for each specific method and mean effect model. The procedure was demonstrated with several popular RNA-seq analysis methods (edgeR, DEseq2, limma+voom). The proposed procedure relaxes the stringent assumptions of data distributions for RNA-seq analysis methods and is illustrated using colorectal cancer studies from The Cancer Genome Atlas database.

Bias in RNA-seq Library Preparation: Current Challenges and Solutions

Although RNA sequencing (RNA-seq) has become the most advanced technology for transcriptome analysis, it also confronts various challenges. As we all know, the workflow of RNA-seq is extremely complicated and it is easy to produce bias. This may damage the quality of RNA-seq dataset and lead to an incorrect interpretation for sequencing result. Thus, our detailed understanding of the source and nature of these biases is essential for the interpretation of RNA-seq data, finding methods to improve the quality of RNA-seq experimental, or development bioinformatics tools to compensate for these biases. Here, we discuss the sources of experimental bias in RNA-seq. And for each type of bias, we discussed the method for improvement, in order to provide some useful suggestions for researcher in RNA-seq experimental.

Performance comparison: exome sequencing as a single test replacing Sanger sequencing

Next generation sequencing tests are used routinely as first-choice tests in the clinic. However, systematic performance comparing the results of exome sequencing as a single test replacing Sanger sequencing of targeted gene(s) is still lacking. Performance comparison data are critically important for clinical case management. In this study, we compared Sanger-sequencing results of 258 genes to those obtained from next generation sequencing (NGS) using two exome-sequencing enrichment kits: Agilent-SureSelectQXT and Illumina-Nextera. Sequencing was performed on leukocytes and buccal-derived DNA from a single individual, and all 258 genes were sequenced a total of 11 times (using different sequencing methods and DNA sources). Sanger sequencing was completed for all exons, including flanking ± 8 bp regions. For the 258 genes, NGS mean coverage was > 20 × for > 98 and > 91% of the regions targeted by SureSelect and Nextera, respectively. Overall, 449 variants were identified in at least one experiment, and 407/449 (90.6%) were detected by all. Of the 42 discordant variants, 23 were determined as true calls, summing-up to a truth set of 430 variants. Sensitivity of true-variant detection was 99% for Sanger sequencing and 97-100% for the NGS experiments. Mean false-positive rates were 3.7E-6 for Sanger sequencing, 2.5E-6 for SureSelect-NGS and 5.2E-6 for Nextera-NGS. Our findings suggest a high overall concordance between Sanger sequencing and NGS performances. Both methods demonstrated false-positive and false-negative calls. High clinical suspicion for a specific diagnosis should, therefore, override negative results of either Sanger sequencing or NGS.

ReSeq simulates realistic Illumina high-throughput sequencing data

In high-throughput sequencing data, performance comparisons between computational tools are essential for making informed decisions at each step of a project. Simulations are a critical part of method comparisons, but for standard Illumina sequencing of genomic DNA, they are often oversimplified, which leads to optimistic results for most tools. ReSeq improves the authenticity of synthetic data by extracting and reproducing key components from real data. Major advancements are the inclusion of systematic errors, a fragment-based coverage model and sampling-matrix estimates based on two-dimensional margins. These improvements lead to more faithful performance evaluations. ReSeq is available at https://github.com/schmeing/ReSeq.

Targeted RNAseq assay incorporating unique molecular identifiers for improved quantification of gene expression signatures and transcribed mutation fraction in fixed tumor samples

BACKGROUND

Our objective was to assess whether modifications to a customized targeted RNA sequencing (RNAseq) assay to include unique molecular identifiers (UMIs) that collapse read counts to their source mRNA counts would improve quantification of transcripts from formalin-fixed paraffin-embedded (FFPE) tumor tissue samples. The assay (SET4) includes signatures that measure hormone receptor and PI3-kinase related transcriptional activity (SETER/PR and PI3Kges), and measures expression of selected activating point mutations and key breast cancer genes.

METHODS

Modifications included steps to introduce eight nucleotides-long UMIs during reverse transcription (RT) in bulk solution, followed by polymerase chain reaction (PCR) of labeled cDNA in droplets, with optimization of the polymerase enzyme and reaction conditions. We used Lin's concordance correlation coefficient (CCC) to measure concordance, including precision (Rho) and accuracy (Bias), and nonparametric tests (Wilcoxon, Levene's) to compare the modified (NEW) SET4 assay to the original (OLD) SET4 assay and to whole transcriptome RNAseq using RNA from matched fresh frozen (FF) and FFPE samples from 12 primary breast cancers.

RESULTS

The modified (NEW) SET4 assay measured single transcripts (p< 0.001) and SETER/PR (p=0.002) more reproducibly in technical replicates from FFPE samples. The modified SET4 assay was more precise for measuring single transcripts (Rho 0.966 vs 0.888, p< 0.01) but not multigene expression signatures SETER/PR (Rho 0.985 vs 0.968) or PI3Kges (Rho 0.985 vs 0.946) in FFPE, compared to FF samples. It was also more precise than wtRNAseq of FFPE for measuring transcripts (Rho 0.986 vs 0.934, p< 0.001) and SETER/PR (Rho 0.993 vs 0.915, p=0.004), but not PI3Kges (Rho 0.988 vs 0.945, p=0.051). Accuracy (Bias) was comparable between protocols. Two samples carried a PIK3CA mutation, and measurements of transcribed mutant allele fraction was similar in FF and FFPE samples and appeared more precise with the modified SET4 assay. Amplification efficiency (reads per UMI) was consistent in FF and FFPE samples, and close to the theoretically expected value, when the library size exceeded 400,000 aligned reads.

CONCLUSIONS

Modifications to the targeted RNAseq protocol for SET4 assay significantly increased the precision of UMI-based and reads-based measurements of individual transcripts, multi-gene signatures, and mutant transcript fraction, particularly with FFPE samples.

Clinical, histopathologic and molecular features of idiopathic and diabetic nodular mesangial sclerosis in humans

BACKGROUND

Idiopathic nodular mesangial sclerosis, also called idiopathic nodular glomerulosclerosis (ING), is a rare clinical entity with an unclear pathogenesis. The hallmark of this disease is the presence of nodular mesangial sclerosis on histology without clinical evidence of diabetes mellitus or other predisposing diagnoses. To achieve insights into its pathogenesis, we queried the clinical, histopathologic and transcriptomic features of ING and nodular diabetic nephropathy (DN).

METHODS

All renal biopsy reports accessioned at Indiana University Health from 2001 to 2016 were reviewed to identify 48 ING cases. Clinical and histopathologic features were compared between individuals with ING and DN (n = 751). Glomeruli of ING (n = 5), DN (n = 18) and reference (REF) nephrectomy (n = 9) samples were isolated by laser microdissection and RNA was sequenced. Immunohistochemistry of proline-rich 36 (PRR36) protein was performed.

RESULTS

ING subjects were frequently hypertensive (95.8%) with a smoking history (66.7%). ING subjects were older, had lower proteinuria and had less hyaline arteriolosclerosis than DN subjects. Butanoate metabolism was an enriched pathway in ING samples compared with either REF or DN samples. The top differentially expressed gene, PRR36, had increased expression in glomeruli 248-fold [false discovery rate (FDR) P = 5.93 × 10-6] compared with the REF and increased 109-fold (FDR P = 1.85 × 10-6) compared with DN samples. Immunohistochemistry revealed a reduced proportion of cells with perinuclear reaction in ING samples as compared to DN.

CONCLUSIONS

Despite similar clinical and histopathologic characteristics in ING and DN, the uncovered transcriptomic signature suggests that ING has distinct molecular features from nodular DN. Further study is warranted to understand these relationships.

Molecular characterization of the human kidney interstitium in health and disease

The gene expression signature of the human kidney interstitium is incompletely understood. The cortical interstitium (excluding tubules, glomeruli, and vessels) in reference nephrectomies (N = 9) and diabetic kidney biopsy specimens (N = 6) was laser microdissected (LMD) and sequenced. Samples underwent RNA sequencing. Gene signatures were deconvolved using single nuclear RNA sequencing (snRNAseq) data derived from overlapping specimens. Interstitial LMD transcriptomics uncovered previously unidentified markers including KISS1, validated with in situ hybridization. LMD transcriptomics and snRNAseq revealed strong correlation of gene expression within corresponding kidney regions. Relevant enriched interstitial pathways included G-protein coupled receptor. binding and collagen biosynthesis. The diabetic interstitium was enriched for extracellular matrix organization and small-molecule catabolism. Cell type markers with unchanged expression (NOTCH3, EGFR, and HEG1) and those down-regulated in diabetic nephropathy (MYH11, LUM, and CCDC3) were identified. LMD transcriptomics complements snRNAseq; together, they facilitate mapping of interstitial marker genes to aid interpretation of pathophysiology in precision medicine studies.

Multiple freeze-thaw cycles lead to a loss of consistency in poly(A)-enriched RNA sequencing

Background

Both RNA-Seq and sample freeze-thaw are ubiquitous. However, knowledge about the impact of freeze-thaw on downstream analyses is limited. The lack of common quality metrics that are sufficiently sensitive to freeze-thaw and RNA degradation, e.g. the RNA Integrity Score, makes such assessments challenging.

Results

Here we quantify the impact of repeated freeze-thaw cycles on the reliability of RNA-Seq by examining poly(A)-enriched and ribosomal RNA depleted RNA-seq from frozen leukocytes drawn from a toddler Autism cohort. To do so, we estimate the relative noise, or percentage of random counts, separating technical replicates. Using this approach we measured noise associated with RIN and freeze-thaw cycles. As expected, RIN does not fully capture sample degradation due to freeze-thaw. We further examined differential expression results and found that three freeze-thaws should extinguish the differential expression reproducibility of similar experiments. Freeze-thaw also resulted in a 3′ shift in the read coverage distribution along the gene body of poly(A)-enriched samples compared to ribosomal RNA depleted samples, suggesting that library preparation may exacerbate freeze-thaw-induced sample degradation.

Conclusion

The use of poly(A)-enrichment for RNA sequencing is pervasive in library preparation of frozen tissue, and thus, it is important during experimental design and data analysis to consider the impact of repeated freeze-thaw cycles on reproducibility.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07381-z.

Molecular intrinsic versus clinical subtyping in breast cancer: A comprehensive review

Breast cancer is a heterogeneous disease manifesting diversity at the molecular, histological and clinical level. The development of breast cancer classification was centered on informing clinical decisions. The current approach to the classification of breast cancer, which categorizes this disease into clinical subtypes based on the detection of estrogen receptor, progesterone receptor, human epidermal growth factor receptor 2, and proliferation marker Ki67, is not ideal. This is manifested as a heterogeneity of therapeutic responses and outcomes within the clinical subtypes. The newer classification model, based on gene expression profiling (intrinsic subtyping) informs about transcriptional responses downstream from IHC single markers, revealing deeper appreciation for the disease heterogeneity and capturing tumor biology in a more comprehensive way than an expression of a single protein or gene alone. While accumulating evidences suggest that intrinsic subtypes provide clinically relevant information beyond clinical surrogates, it is imperative to establish whether the current conventional immunohistochemistry-based clinical subtyping approach could be improved by gene expression profiling and if this approach has a potential to translate into clinical practice.

RNA sequencing: new technologies and applications in cancer research

Over the past few decades, RNA sequencing has significantly progressed, becoming a paramount approach for transcriptome profiling. The revolution from bulk RNA sequencing to single-molecular, single-cell and spatial transcriptome approaches has enabled increasingly accurate, individual cell resolution incorporated with spatial information. Cancer, a major malignant and heterogeneous lethal disease, remains an enormous challenge in medical research and clinical treatment. As a vital tool, RNA sequencing has been utilized in many aspects of cancer research and therapy, including biomarker discovery and characterization of cancer heterogeneity and evolution, drug resistance, cancer immune microenvironment and immunotherapy, cancer neoantigens and so on. In this review, the latest studies on RNA sequencing technology and their applications in cancer are summarized, and future challenges and opportunities for RNA sequencing technology in cancer applications are discussed.

RNAseq Analysis of Rodent Spaceflight Experiments Is Confounded by Sample Collection Techniques

To understand the physiological changes that occur in response to spaceflight, mice are transported to the International Space Station (ISS) and housed for variable periods of time before euthanasia on-orbit or return to Earth. Sample collection under such difficult conditions introduces confounding factors that need to be identified and addressed. We found large changes in the transcriptome of mouse tissues dissected and preserved on-orbit compared with tissues from mice euthanized on-orbit, preserved, and dissected after return to Earth. Changes due to preservation method eclipsed those between flight and ground samples, making it difficult to identify spaceflight-specific changes. Follow-on experiments to interrogate the roles of euthanasia methods, tissue and carcass preservation protocols, and library preparation methods suggested that differences due to preservation protocols are exacerbated when coupled with polyA selection. This has important implications for the interpretation of existing datasets and the design of future experiments.

•

Experimentation is necessary to understand how organisms respond to space

•

Specialized protocols are used for preserving biological samples on the ISS

•

RNAseq datasets are impacted by preservation protocols used on the ISS

•

Impacts can be alleviated with improved carcass preservation protocols

Utility of Extrapolating Human S1500+ Genes to the Whole Transcriptome: Tunicamycin Case Study

The TempO-Seq S1500+ platform(s), now available for human, mouse, rat, and zebrafish, measures a discrete number of genes that are representative of biological and pathway co-regulation across the entire genome in a given species. While measurement of these genes alone provides a direct assessment of gene expression activity, extrapolating expression values to the whole transcriptome (~26 000 genes in humans) can estimate measurements of non-measured genes of interest and increases the power of pathway analysis algorithms by using a larger background gene expression space. Here, we use data from primary hepatocytes of 54 donors that were treated with the endoplasmic reticulum (ER) stress inducer tunicamycin and then measured on the human S1500+ platform containing ~3000 representative genes. Measurements for the S1500+ genes were then used to extrapolate expression values for the remaining human transcriptome. As a case study of the improved downstream analysis achieved by extrapolation, the “measured only” and “whole transcriptome” (measured + extrapolated) gene sets were compared. Extrapolation increased the number of significant genes by 49%, bringing to the forefront many that are known to be associated with tunicamycin exposure. The extrapolation procedure also correctly identified established tunicamycin-related functional pathways reflected by coordinated changes in interrelated genes while maintaining the sample variability observed from the “measured only” genes. Extrapolation improved the gene- and pathway-level biological interpretations for a variety of downstream applications, including differential expression analysis, gene set enrichment pathway analysis, DAVID keyword analysis, Ingenuity Pathway Analysis, and NextBio correlated compound analysis. The extrapolated data highlight the role of metabolism/metabolic pathways, the ER, immune response, and the unfolded protein response, each of which are key activities associated with tunicamycin exposure that were unrepresented or underrepresented in one or more of the analyses of the original “measured only” dataset. Furthermore, the inclusion of the extrapolated genes raised “tunicamycin” from third to first upstream regulator in Ingenuity Pathway Analysis and from sixth to second most correlated compound in NextBio analysis. Therefore, our case study suggests an approach to extend and enhance data from the S1500+ platform for improved insight into biological mechanisms and functional outcomes of diseases, drugs, and other perturbations.

Physiological RNA dynamics in RNA-Seq analysis

Physiological RNA dynamics cause problems in transcriptome analysis. Physiological RNA accumulation affects the analysis of RNA quantification, and physiological RNA degradation affects the analysis of the RNA sequence length, feature site and quantification. In the present article, we review the effects of physiological degradation and accumulation of RNA on analysing RNA sequencing data. Physiological RNA accumulation and degradation probably led to such phenomena as incorrect estimations of transcription quantification, differential expressions, co-expressions, RNA decay rates, alternative splicing, boundaries of transcription, novel genes, new single-nucleotide polymorphisms, small RNAs and gene fusion. Thus, the transcriptomic data obtained up to date warrant further scrutiny. New and improved techniques and bioinformatics software are needed to produce accurate data in transcriptome research.

A Comprehensive Coexpression Network Analysis in Vibrio cholerae

Cholera is a devastating illness that kills tens of thousands of people annually. Vibrio cholerae, the causative agent of cholera, is an important model organism to investigate both bacterial pathogenesis and the impact of horizontal gene transfer on the emergence and dissemination of new virulent strains. Despite the importance of this pathogen, roughly one-third of V. cholerae genes are functionally unannotated, leaving large gaps in our understanding of this microbe. Through coexpression network analysis of existing RNA sequencing data, this work develops an approach to uncover novel gene-gene relationships and contextualize genes with no known function, which will advance our understanding of V. cholerae virulence and evolution.

Research into the evolution and pathogenesis of Vibrio cholerae has benefited greatly from the generation of high-throughput sequencing data to drive molecular analyses. The steady accumulation of these data sets now provides a unique opportunity for in silico hypothesis generation via coexpression analysis. Here, we leverage all published V. cholerae RNA sequencing data, in combination with select data from other platforms, to generate a gene coexpression network that validates known gene interactions and identifies novel genetic partners across the entire V. cholerae genome. This network provides direct insights into genes influencing pathogenicity, metabolism, and transcriptional regulation, further clarifies results from previous sequencing experiments in V. cholerae (e.g., transposon insertion sequencing [Tn-seq] and chromatin immunoprecipitation sequencing [ChIP-seq]), and expands upon microarray-based findings in related Gram-negative bacteria.

IMPORTANCE Cholera is a devastating illness that kills tens of thousands of people annually. Vibrio cholerae, the causative agent of cholera, is an important model organism to investigate both bacterial pathogenesis and the impact of horizontal gene transfer on the emergence and dissemination of new virulent strains. Despite the importance of this pathogen, roughly one-third of V. cholerae genes are functionally unannotated, leaving large gaps in our understanding of this microbe. Through coexpression network analysis of existing RNA sequencing data, this work develops an approach to uncover novel gene-gene relationships and contextualize genes with no known function, which will advance our understanding of V. cholerae virulence and evolution.

The NASA Twins Study: The Effect of One Year in Space on Long-Chain Fatty Acid Desaturases and Elongases

BACKGROUND

At present, there is no clear understanding of the effect of long-duration spaceflight on the major enzymes that govern the metabolism of omega-6 and omega-3 fatty acids. To address this gap in knowledge, we used data from the NASA Twins Study, which includes a multiscale omics investigation of the changes that occurred during a year-long (340 days) human spaceflight. Embedded within the NASA Twins data are specific analytes associated with fatty acid metabolism.

OBJECTIVES

To examine the long-chain fatty acid desaturases and elongases in a single human during 1 year in space.

METHOD

One male twin was on board the International Space Station (ISS) for 1 year, while his monozygotic twin served as a genetically matched ground control. Longitudinal assessments included the genome, epige-nome, transcriptome, proteome, metabolome, microbiome, and immunome during the mission, as well as 6 months before and after. The gene-specific fatty acid desaturase and elongase transcriptome data (FADS1, FADS2, ELOVL2, and ELOVL5) were extracted from untargeted RNA-seq measurements derived from white blood cell fractions.

RESULTS

Most data from the elongases and desaturases exhibited relatively similar expression profiles (R2 >0.6) over time for the CD8, CD19, and lymphocyte-depleted (LD) cell fractions, indicating overall conservation of function within and between the subjects. Both cell-type and temporal specificity was observed in some cases, and some differences were also apparent between the polyadenylated (polyA) fraction of processed RNAs versus the ribodepleted (ribo-) fraction. The flight subject showed a stronger enrichment of the fatty acid metabolic process pathway across almost all cell types (columns, CD4, CD8, CPT, and LD), most especially in the ribodepleted fraction of RNA, but also with the polyA+ fraction of RNA. Gene set enrichment analysis (GSEA) measures across three related fatty acid metabolism pathways showed a differential between the ground and the flight subject.

CONCLUSIONS

There appears to be no persistent alteration of desaturase and elongase gene expression associated with 1 year in space. However, these data provide evidence that cellular lipid metabolism can be responsive and dynamic to spaceflight, even though it appears cell-type and context specific, most notably in terms of the fraction of RNA measured and the collection protocols. These results also provide new evidence of mid-flight spikes in expression of selected genes, which may indicate transient responses to specific insults during spaceflight.

A Review of the Scientific Rigor, Reproducibility, and Transparency Studies Conducted by the ABRF Research Groups

Shared research resource facilities, also known as core laboratories (Cores), are responsible for generating a significant and growing portion of the research data in academic biomedical research institutions. Cores represent a central repository for institutional knowledge management, with deep expertise in the strengths and limitations of technology and its applications. They inherently support transparency and scientific reproducibility by protecting against cognitive bias in research design and data analysis, and they have institutional responsibility for the conduct of research (research ethics, regulatory compliance, and financial accountability) performed in their Cores. The Association of Biomolecular Resource Facilities (ABRF) is a FASEB-member scientific society whose members are scientists and administrators that manage or support Cores. The ABRF Research Groups (RGs), representing expertise for an array of cutting-edge and established technology platforms, perform multicenter research studies to determine and communicate best practices and community-based standards. This review provides a summary of the contributions of the ABRF RGs to promote scientific rigor and reproducibility in Cores from the published literature, ABRF meetings, and ABRF RGs communications.