RNA-seq Libraries Research Articles

Affordable, accurate and unbiased RNA sequencing by manual library miniaturization: A case study in barley.

We present an easy-to-reproduce manual miniaturized full-length RNA sequencing (RNAseq) library preparation workflow that does not require the upfront investment in expensive lab equipment or long setup times. With minimal adjustments to an established commercial protocol, we were able to manually miniaturize the RNAseq library preparation by a factor of up to 1:8. This led to cost savings for miniaturized library preparation of up to 86.1% compared to the gold standard. The resulting data were the basis of a rigorous quality control analysis that inspected: sequencing quality metrics, gene body coverage, raw read duplications, alignment statistics, read pair duplications, detected transcripts and sequence variants. We also included a deep dive data analysis identifying rRNA contamination and suggested ways to circumvent these. In the end, we could not find any indication of biases or inaccuracies caused by the RNAseq library miniaturization. The variance in detected transcripts was minimal and not influenced by the miniaturization level. Our results suggest that the workflow is highly reproducible and the sequence data suitable for downstream analyses such as differential gene expression analysis or variant calling.

The lncRNA-mediated ceRNA network of Altica viridicyanea is involved in the regulation of the Toll/Imd signaling pathway under antibiotic treatment.

Long noncoding RNAs (lncRNAs) play significant roles in the regulation of mRNA expression or in shaping the competing endogenous RNA (ceRNA) network by targeting miRNA. The insect gut is one of the most important tissues due to direct contact with external pathogens and functions in the immune defense against pathogen infection through the innate immune system and symbionts, but there are limited observations on the role of the lncRNA-involved ceRNA network of the Toll/Imd pathway and correlation analysis between this network and bacterial microbiota in the Altica viridicyanea gut. In this research, we constructed and sequenced six RNA sequencing libraries using normal and antibiotic-reared samples, generating a total of 17,193 lncRNAs and 26,361 mRNAs from massive clean data by quality control and bioinformatic analysis. Furthermore, a set of 8,539 differentially expressed lncRNAs (DELs) and 13,263 differentially expressed mRNAs (DEMs), of which related to various immune signaling pathways, such as the Toll/Imd, JAK/STAT, NF-κB, and PI3K-Akt signaling pathways, were obtained between the two experimental groups in A. viridicyanea. In addition, numerous GO and KEGG enrichment analyses were used to annotate the DELs and their target genes. Moreover, six Toll family members and nineteen signal genes from the Toll/Imd signaling pathway were identified and characterized using online tools, and phylogenetic analyses of the above genes proved their classification. Next, a lncRNA-miRNA-mRNA network of the Toll/Imd pathway was built, and it contained different numbers of DEMs in this pathway and related DELs based on prediction and annotation. In addition, qRT-PCR validation and sequencing data were conducted to show the expression patterns of the above DELs and DEMs related to the Toll/Imd signaling pathway. Finally, the correlated investigations between DELs or DEMs of the Toll/Imd signaling pathway and most changes in the gut bacterial microbiota revealed significantly positive or negative relationships between them. The present findings provide essential evidence for innate immune ceRNAs in the beetle gut and uncover new potential relationships between innate immune pathways and the gut bacterial microbiota in insects.

Differential gene expression of Asian citrus psyllids infected with 'Ca. Liberibacter asiaticus' reveals hyper-susceptibility to invasion by instar fourth-fifth and teneral adult stages.

The bacterial pathogen Candidatus Liberibacter asiaticus (CLas) is the causal agent of citrus greening disease. This unusual plant pathogenic bacterium also infects its psyllid host, the Asian citrus psyllid (ACP). To investigate gene expression profiles with a focus on genes involved in infection and circulation within the psyllid host of CLas, RNA-seq libraries were constructed from CLas-infected and CLas-free ACP representing the five different developmental stages, namely, nymphal instars 1-2, 3, and 4-5, and teneral and mature adults. The Gbp paired-end reads (296) representing the transcriptional landscape of ACP across all life stages and the official gene set (OGSv3) were annotated based on the chromosomal-length v3 reference genome and used for de novo transcript discovery resulting in 25,410 genes with 124,177 isoforms. Differential expression analysis across all ACP developmental stages revealed instar-specific responses to CLas infection, with greater overall responses by nymphal instars, compared to mature adults. More genes were over-or under-expressed in the 4-5th nymphal instars and young (teneral) adults than in instars 1-3, or mature adults, indicating that late immature instars and young maturing adults were highly responsive to CLas infection. Genes identified with potential for direct or indirect involvement in the ACP-CLas circulative, propagative transmission pathway were predominantly responsive during early invasion and infection processes and included canonical cytoskeletal remodeling and endo-exocytosis pathway genes. Genes with predicted functions in defense, development, and immunity exhibited the greatest responsiveness to CLas infection. These results shed new light on ACP-CLas interactions essential for pathogenesis of the psyllid host, some that share striking similarities with effector protein-animal host mechanisms reported for other culturable and/or fastidious bacterial- or viral- host pathosystems.

Comparative phylogenomics and phylotranscriptomics provide insights into the genetic complexity of nitrogen-fixing root-nodule symbiosis

Plant root-nodule symbiosis (RNS) with mutualistic nitrogen-fixing bacteria is restricted to a single clade of angiosperms, the Nitrogen-Fixing Nodulation Clade (NFNC), and is best understood in the legume family. Nodulating species share many commonalities, explained either by divergence from a common ancestor over 100 million years ago or by convergence following independent origins over that same time period. Regardless, comparative analyses of diverse nodulation syndromes can provide insights into constraints on nodulation—what must be acquired or cannot be lost for a functional symbiosis—and the latitude for variation in the symbiosis. However, much remains to be learned about nodulation, especially outside of legumes. Here, we employed a large-scale phylogenomic analysis across 88 species, complemented by 151 RNA-seq libraries, to elucidate the evolution of RNS. Our phylogenomic analyses further emphasize the uniqueness of the transcription factor NIN as a master regulator of nodulation and identify key mutations that affect its function across the NFNC. Comparative transcriptomic assessment revealed nodule-specific upregulated genes across diverse nodulating plants, while also identifying nodule-specific and nitrogen-response genes. Approximately 70% of symbiosis-related genes are highly conserved in the four representative species, whereas defense-related and host-range restriction genes tend to be lineage specific. Our study also identified over 900 000 conserved non-coding elements (CNEs), over 300 000 of which are unique to sampled NFNC species. NFNC-specific CNEs are enriched with the active H3K9ac mark and are correlated with accessible chromatin regions, thus representing a pool of candidate regulatory elements for genes involved in RNS. Collectively, our results provide novel insights into the evolution of nodulation and lay a foundation for engineering of RNS traits in agriculturally important crops.

Metabolic Pathway Reconstruction Indicates the Presence of Important Medicinal Compounds in Coffea Such as L-DOPA.

The use of transcriptomic data to make inferences about plant metabolomes is a useful tool to help the discovery of important compounds in the available biodiversity. To unveil previously undiscovered metabolites of Coffea, of phytotherapeutic and economic value, we employed 24 RNAseq libraries. These libraries were sequenced from leaves exposed to a diverse range of environmental conditions. Subsequently, the data were meticulously processed to create models of putative metabolic networks, which shed light on the production of potential natural compounds of significant interest. Then, we selected one of the predicted compounds, the L-3,4-dihydroxyphenylalanine (L-DOPA), to be analyzed by LC-MS/MS using three biological replicates of flowers, leaves, and fruits from Coffea arabica and Coffea canephora. We were able to identify metabolic pathways responsible for producing several compounds of economic importance. One of the identified pathways involved in isoquinoline alkaloid biosynthesis was found to be active and producing L-DOPA, which is a common product of POLYPHENOL OXIDASES (PPOs, EC 1.14.18.1 and EC 1.10.3.1). We show that coffee plants are a natural source of L-DOPA, a widely used medicine for treatment of the human neurodegenerative condition called Parkinson's disease. In addition, dozens of other compounds with medicinal significance were predicted as potential natural coffee products. By further refining analytical chemistry techniques, it will be possible to enhance the characterization of coffee metabolites, enabling a deeper understanding of their properties and potential applications in medicine.

The chromatin and single-cell transcriptional landscapes of CD4 T cells in inflammatory bowel disease link risk loci with a proinflammatory Th17 cell population.

The imbalance between Th17 and regulatory T cells in inflammatory bowel diseases (IBD) promotes intestinal epithelial cell damage. In this scenario, T helper cell lineage commitment is accompanied by dynamic changes to the chromatin that facilitate or repress gene expression. Here, we characterized the chromatin landscape and heterogeneity of intestinal and peripheral CD4 T cellsfrom IBD patients using in house ATAC-Seq and single cell RNA-Seq libraries. We show that chromatin accessibility profiles of CD4 T cells from inflamed intestinal biopsies relate to genes associated with a network of inflammatory processes. After integrating the chromatin profiles of tissue-derived CD4 T cells and in-vitro polarized CD4 T cell subpopulations, we found that the chromatin accessibility changes of CD4 T cells were associated with a higher predominance of pathogenic Th17 cells (pTh17 cells) in inflamed biopsies. In addition, IBD risk loci in CD4 T cells were colocalized with accessible chromatin changes near pTh17-related genes, as shown in intronic STAT3 and IL23R regions enriched in areas of active intestinal inflammation. Moreover, single cell RNA-Seq analysis revealed a population of pTh17 cells that co-expresses Th1 and cytotoxic transcriptional programs associated with IBD severity. Altogether, we show that cytotoxic pTh17 cells were specifically associated with IBD genetic variants and linked to intestinal inflammation of IBD patients.

Unraveling the regulatory role of MYC2 on ASMT gene expression in wheat: Implications for melatonin biosynthesis and drought tolerance.

Recognized for its multifaceted functions, melatonin is a hormone found in both animals and plants. In the plant kingdom, it plays diverse roles, regulating growth, development, and stress responses. Notably, melatonin demonstrates its significance by mitigating the effects of abiotic stresses like drought. However, understanding the precise regulatory mechanisms controlling melatonin biosynthesis genes, especially during monocots' response to stresses, requires further exploration. Seeking to understand the molecular basis of drought stress tolerance in wheat, we analyzed RNA-Seq libraries of wheat exposed to drought stress using bioinformatics methods. In light of our findings, we identified that the Myelocytomatosis oncogenes 2 (MYC2) transcription factor is a hub gene upstream of a main melatonin biosynthesis gene, N-acetylserotonin methyltransferase (ASMT), in the wheat drought response-gene network. Promoter analysis of the ASMT gene suggested that it might be a target gene of MYC2. We conducted a set of molecular and physiochemical assays along with robust machine learning approaches to elevate those findings further. MYC2 and ASMT were co-regulated under Jasmonate, drought, and a combination of them in the leaf tissues of wheat was detected. A meaningful correlation was observed among gene expression profiles, melatonin contents, photosynthetic activities, antioxidant enzyme activities, H2 O2 levels, and plasma membrane damage. The results indicated an evident relationship between jasmonic acid and the melatonin biosynthesis pathway. Moreover, it seems that the MYC2-ASMT module might contribute to wheat drought tolerance by regulating melatonin contents.

Transcriptomic and Functional Analyses of Two Cadmium Hyper-Enriched Duckweed Strains Reveal Putative Cadmium Tolerance Mechanisms.

Cadmium (Cd) is one of the most toxic metals in the environment and exerts deleterious effects on plant growth and production. Duckweed has been reported as a promising candidate for Cd phytoremediation. In this study, the growth, Cd enrichment, and antioxidant enzyme activity of duckweed were investigated. We found that both high-Cd-tolerance duckweed (HCD) and low-Cd-tolerance duckweed (LCD) strains exposed to Cd were hyper-enriched with Cd. To further explore the underlying molecular mechanisms, a genome-wide transcriptome analysis was performed. The results showed that the growth rate, chlorophyll content, and antioxidant enzyme activities of duckweed were significantly affected by Cd stress and differed between the two strains. In the genome-wide transcriptome analysis, the RNA-seq library generated 544,347,670 clean reads, and 1608 and 2045 differentially expressed genes were identified between HCD and LCD, respectively. The antioxidant system was significantly expressed during ribosomal biosynthesis in HCD but not in LCD. Fatty acid metabolism and ethanol production were significantly increased in LCD. Alpha-linolenic acid metabolism likely plays an important role in Cd detoxification in duckweed. These findings contribute to the understanding of Cd tolerance mechanisms in hyperaccumulator plants and lay the foundation for future phytoremediation studies.

Peonidin-3-O-glucoside and Resveratrol Increase the Viability of Cultured Human hFOB Osteoblasts and Alter the Expression of Genes Associated with Apoptosis, Osteoblast Differentiation and Osteoclastogenesis.

High-throughput RNA-sequencing can determine the impact of nutrients and their combinations on gene transcription levels in osteocytes, and clarify the biological pathways associated with their impact on bone tissues. Previously, we reported that resveratrol (RES) and peonidin-3-O-glucoside (POG) increased osteoblastogenesis, as well as reduced osteoclastogenesis in transgenic teleost fish models. Here, we perform whole-genome transcriptomic profiling of osteoblasts treated with POG or RES to provide a comprehensive understanding of alterations in gene expression and the molecular mechanisms involved. Cultured human fetal osteoblastic hFOB 1.19 cells were treated with the test compounds, and then RNA was used to prepare RNA-seq libraries, that were sequenced using a NovaSeq 6000. Treatment with POG or RES increased osteoblast proliferation and reduced apoptosis. Transcriptomic profiling showed that of the 29,762 genes investigated, 3177 were differentially expressed (1481 upregulated, 1696 downregulated, FDR ≤ 0.05) in POG-treated osteoblasts. In the RES-treated osteoblasts, 2288 genes were differentially expressed (DGEs, 1068 upregulated, 1220 downregulated, FDR ≤ 0.05). Ingenuity® Pathway Analysis (IPA) of DGEs from RES or POG-treated osteoblasts revealed significant downregulation of the apoptosis, osteoarthritis and HIF1α canonical pathways, and a significant reduction in Rankl mRNA expression. The data suggest that RES and POG have both anabolic and anticlastogenic effects.

Fixation and Laser Capture Microdissection of Plant Tissue for RNA Extraction and RNASeq Library Preparation.

To study the transcriptome of individual plant cells at specific points in time, we developed protocols for fixation, embedding, and sectioning of plant tissue followed by laser capture microdissection (LCM) and processing for RNA recovery. LCM allows the isolation of individual cell types from heterogeneous tissue sections and is particularly suited to plant processing because it does not require the breakdown of cell walls. This approach allows accurate separation of a small volume of cells that can be used to study gene expression profiles in different tissues or cell layers. The technique requires neither separation of cells by enzymatic digestion of any kind nor cell-specific reporter genes, and it allows storage of fixed and embedded tissue for months before capture. The methods for fixation, embedding, sectioning, and capturing of plant cells that we describe yield high-quality RNA suitable for making libraries for RNASeq. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Tissue Preparation for Laser Capture Microdissection Basic Protocol 2: Tissue Sectioning Basic Protocol 3: Laser Capture Microdissection of Embedded Tissue Basic Protocol 4: RNA Extraction from Laser Capture Microdissection Samples.

Meta-analysis of Diurnal Transcriptomics in Mouse Liver Reveals Low Repeatability of Rhythm Analyses.

To assess the consistency of biological rhythms across studies, 57 public mouse liver tissue timeseries totaling 1096 RNA-seq samples were obtained and analyzed. Only the control groups of each study were included, to create comparable data. Technical factors in RNA-seq library preparation were the largest contributors to transcriptome-level differences, beyond biological or experiment-specific factors such as lighting conditions. Core clock genes were remarkably consistent in phase across all studies. Overlap of genes identified as rhythmic across studies was generally low, with no pair of studies having over 60% overlap. Distributions of phases of significant genes were remarkably inconsistent across studies, but the genes that consistently identified as rhythmic had acrophase clustering near ZT0 and ZT12. Despite the discrepancies between single-study analyses, cross-study analyses found substantial consistency. Running compareRhythms on each pair of studies identified a median of only 11% of the identified rhythmic genes as rhythmic in only 1 of the 2 studies. Data were integrated across studies in a joint and individual variance estimate (JIVE) analysis, which showed that the top 2 components of joint within-study variation are determined by time of day. A shape-invariant model with random effects was fit to the genes to identify the underlying shape of the rhythms, consistent across all studies, including identifying 72 genes with consistently multiple peaks.

Neural circuit-wide analysis of changes to gene expression during deafening-induced birdsong destabilization.

Sensory feedback is required for the stable execution of learned motor skills, and its loss can severely disrupt motor performance. The neural mechanisms that mediate sensorimotor stability have been extensively studied at systems and physiological levels, yet relatively little is known about how disruptions to sensory input alter the molecular properties of associated motor systems. Songbird courtship song, a model for skilled behavior, is a learned and highly structured vocalization that is destabilized following deafening. Here, we sought to determine how the loss of auditory feedback modifies gene expression and its coordination across the birdsong sensorimotor circuit. To facilitate this system-wide analysis of transcriptional responses, we developed a gene expression profiling approach that enables the construction of hundreds of spatially-defined RNA-sequencing libraries. Using this method, we found that deafening preferentially alters gene expression across birdsong neural circuitry relative to surrounding areas, particularly in premotor and striatal regions. Genes with altered expression are associated with synaptic transmission, neuronal spines, and neuromodulation and show a bias toward expression in glutamatergic neurons and Pvalb/Sst-class GABAergic interneurons. We also found that connected song regions exhibit correlations in gene expression that were reduced in deafened birds relative to hearing birds, suggesting that song destabilization alters the inter-region coordination of transcriptional states. Finally, lesioning LMAN, a forebrain afferent of RA required for deafening-induced song plasticity, had the largest effect on groups of genes that were also most affected by deafening. Combined, this integrated transcriptomics analysis demonstrates that the loss of peripheral sensory input drives a distributed gene expression response throughout associated sensorimotor neural circuitry and identifies specific candidate molecular and cellular mechanisms that support the stability and plasticity of learned motor skills.

Abstract PR001: An endogenous molecular brake preventing APOBEC-driven tumor mutational burden, heterogeneity and AR therapy resistance

Abstract Tumor mutational burden and heterogeneity fuel resistance to many targeted therapies, which inevitably limits the clinical outcome of patients. The cytosine deaminase APOBEC proteins are major drivers of mutagenesis in human cancer, with over 70% of tumors exhibiting a mutational signature that is impacted by APOBEC activity. Nevertheless, the mechanisms through which tumor cells hijack the powerful APOBEC mutagenesis machinery to promote heterogeneity and therapy resistance remain largely unknown. Through a multi-disciplinary approach integrating bulk and single cell RNA-Seq (scRNA-Seq), whole-genome exome-sequencing (WES), and CRISPR library screening, we identified a long sought-after cell-intrinsic mechanism that prevents APOBEC-driven mutagenesis in non-malignant cells. The loss of this “molecular brake” unleashes APOBEC3B-driven mutagenesis in malignant cells, which then becomes a key mutator and represents the long sought-after molecular source of driver mutations in some frequently mutated genes in cancers, including FOXA1, EP300, and AR. Functional screening identified eight crucial drivers for Androgen Receptor (AR)-targeted therapy resistance in prostate cancer that are mutated by APOBEC3B: BRD7, CBX8, EP300, FOXA1, HDAC5, HSF4, STAT3, and AR. Finally, contrary to conventional understanding, our study reveals that driver mutations in FOXA1 induced by APOBEC3B, not mutations in AR, evolutionarily outcompete other driver mutations and eventually dominate the resistant tumors. Collectively, these results uncover a cell-intrinsic mechanism that unleashes APOBEC-driven mutagenesis, which plays a significant role in conferring targeted therapy resistance and could be the potential therapeutic targets to overcome resistance. Citation Format: Ping Mu, Xiaoling Li, Yunguan Wang. An endogenous molecular brake preventing APOBEC-driven tumor mutational burden, heterogeneity and AR therapy resistance [abstract]. In: Proceedings of the AACR Special Conference: Advances in Prostate Cancer Research; 2023 Mar 15-18; Denver, Colorado. Philadelphia (PA): AACR; Cancer Res 2023;83(11 Suppl):Abstract nr PR001.

Abstract B021: Examining the role of HER2 signaling in promoting prostate cancer proliferation in an androgen receptor-low environment

Abstract Background: Intense neoadjuvant androgen deprivation therapy (inADT) reduces prostate cancer recurrence in approximately 40% of patients with intermediate/high risk prostate cancer. Using samples from our recent phase 2 study (NCT02430480) from patients undergoing 6 months of inADT prior to radical prostatectomy, we sought to examine molecular features present at baseline that correlate with pathologic response following surgery. Methods: Tumors from 37 patients with locally advanced prostate cancer were subjected to immunohistochemistry (IHC) and laser capture microdissection to obtain 143 histologically distinct foci. RNA was isolated from each focus and used to generate RNA-seq libraries. Differentially-expressed genes were determined per-unit of residual tumor volume using linear mixed-effect models. Histologic features and staining intensities were obtained from IHC of posttreatment specimens. Ingenuity Pathway Analysis (IPA) identified upstream regulators at baseline that interacted with pathologic outcome. Cancer cell lines were treated with a panel of EGFR family inhibitors and androgen receptor (AR) axis inhibitors. Protein and RNA expression were measured to identify the impact of dual inhibition. Results: In pre-treatment targeted biopsies from men with locally advanced prostate cancer, we observed that tumor samples that exhibited worse response to inADT had lower AR activity at baseline, presumably due to a decreased dependency on AR, potentially related to increased dependency on alternative growth pathways. These incomplete/nonresponding tumors displayed a significant increase in expression of EFGR family signaling (EGFR, HER2, and HER3) using IPA, which have previously been shown to serve as bypass mechanisms to AR signaling to drive downstream PI3K activity. HER2 was the most abundantly expressed protein by IHC in posttreatment tumors. Given the multitude of receptor tyrosine kinases (RTKs) regulated by EGFR family signaling and RTK inhibitors targeting this activity, we examined the impact of combining RTK inhibitors with ADT to overcome resistance to ADT. Afatinib, lapatinib, and neratinib showed the greatest anti-tumor efficacy in vitro. Treatment with these agents rapidly and dynamically increased AR expression and activity. We observed synergistic effects on cell viability and proliferation by treating cell lines with combination therapy of antiandrogen therapy and these inhibitors. Conclusions: Our data suggest that elevated HER2 activity defines an AR-low subtype of prostate cancer that is resistant to inADT. Disrupting the HER2-mediated bypass of AR blockade re-sensitized resistant PC cells to AR-targeted therapy. These effects were observed in an untreated setting and represent a divergent path in tumorigenesis where prostate tumors develop with intrinsically lower AR activity. These findings will be validated using genetic abrogation of HER2 in additional clinical samples, organoids, and patient-derived xenograft models, and represent new opportunities for combination therapy in newly-diagnosed high-risk disease. Citation Format: Scott Wilkinson, Anson Ku, Isaiah King, Anna Baj, Daniel Low, Huihui Ye, Stephanie A. Harmon, Nicholas T. Terrigino, John R. Bright, Rayann Atway, Shana Y. Trostel, Rosina T. Lis, Baris Turkbey, William L. Dahut, Fatima Karzai, Adam G. Sowalsky. Examining the role of HER2 signaling in promoting prostate cancer proliferation in an androgen receptor-low environment [abstract]. In: Proceedings of the AACR Special Conference: Advances in Prostate Cancer Research; 2023 Mar 15-18; Denver, Colorado. Philadelphia (PA): AACR; Cancer Res 2023;83(11 Suppl):Abstract nr B021.

Establishing clinical performance characteristics for Duoseq, a new assay for DNA and RNA sequencing.

e15059 Background: Despite advances in genomics and next generation sequencing (NGS), clinical implementation of NGS remains challenging. This challenge is critical in hematologic malignancies where NGS has become an essential part of the clinical workup. Many clinical labs resort to send-out testing that can delay appropriate treatment and drive up the cost of care. To address these issues, we have developed Duoseq as a unified approach for DNA and RNA NGS. Duoseq assesses mutations (single nucleotide variants [SNVs], insertions/deletions [indels] and splice site variants [SSVs]), structural variants, and gene expression using a panel of over 400 genes. It is supplied as a kit that enables preparation of DNA and RNA sequencing libraries in a single workflow performed at the point-of-care. It includes complete bioinformatics support that enables a clinical laboratory with an NGS instrument to perform the assay with a turn-around time of 2-3 days. Duoseq is optimized to work with formalin-fixed, paraffin embedded (FFPE) samples that are the standard fixation practice. Methods: We evaluated Duoseq on four dimensions critical to clinical performance: reproducibility, precision, analytical sensitivity (limit of detection, LoD), and analytical specificity. Identical experiments were run at clinical laboratories of two different institutions. Analytical validation samples included 4 FFPE cell lines and 2 synthetic DNA samples designed to harbor variants of interest at known allele frequencies. Results: We determined the LoD to be 5% mutant allele fraction (MAF) for SNVs, 10% MAF for indels/SSVs, and 10% tumor purity for SVs. Duoseq achieved 100% precision and reproducibility for each variant class for events at or above the LoD. Analytical specificity was 100%, indicating no instances of cross-contamination of samples within or across runs. We further evaluated Duoseq performance in 54 clinical cases of myeloid and lymphoid malignancies compared against clinical results provided by FoundationOne Heme NGS panel. Collectively, we assessed 119 SNVs and 49 indels/SSVs. Duoseq performed with 100% specificity for both SNVs and indels/SSVs. Duoseq achieved 98% sensitivity for SNVs and 96% for indels/SSVs. For example, Duoseq detected TP53 variants in 6 cases of myeloid neoplasms, which has relevance for the 2022 WHO and ICC classifications. Duoseq also identified therapeutically relevant variants in IDH1/2, EZH2, FLT3, KRAS and TET2. Conclusions: This study demonstrates performance characteristics that suggest Duoseq would perform well as a clinical-grade assay. Local laboratory factors that include staffing, case volumes, and reimbursement affect adoption and turnaround times in the individual clinical laboratory. However, we anticipate that Duoseq implementation will accelerate the diagnosis of hematologic malignancies and facilitate rapid, more cost-effective care to improve patient outcomes.

Artificial intelligence analysis of CCNG1 expression in tumors: A novel biomarker in development for anti-CCNG1 targeted therapy.

e15083 Background: Metastatic cancer is associated with an invariably fatal outcome. Therefore, innovative therapies are urgently needed. Although expanded access for DeltaRex-G, a tumor targeted retrovector encoding a CCNG1 inhibitor gene, is on-going for an intermediate-size (n = up to 40) population of advanced sarcoma and pancreatic cancer, more data is needed to identify patients who are likely to benefit from DeltaRex-G gene therapy. In this study, we retrospectively analyzed CCNG1 expression in archived tumors of patients who were previously treated with DeltaRex-G and who are active candidates for DeltaRex-G therapy. Methods: Archived formalin-fixed paraffin embedded (FFPE) tumor specimens (n = 58) from patients with solid malignancies who are actively followed at the Cancer Center of Southern California were collected, processed, and subjected to RNA sequencing. Briefly, RNA-seq libraries were sequenced to generate 50 million reads that were aligned using Kallisto v0.42.4 to GENCODE v23 transcripts with default parameters. Only protein-coding, IGH/K/L- and TCR-related transcripts were retained for downstream processing, resulting in 20,062 protein coding genes. Gene expression was quantified as transcripts per million (TPM) and log2-transformed (Goldman et al., 2020). A gene expression level is presented as low/medium/high depending on the expression level of a such gene in patients of the reference cohort. Low = < 17%; Medium = 17%-83%; High = > 83%. Results: Thirty-two male and 26 female subjects, ages ranging from 16 to 86 years were studied. Forty-nine (84.4%) patients had sarcoma, 3 (5.2%) had urothelial carcinoma, 2 (3.5%) had breast carcinoma, 2 (3.5%) had pancreatic cancer, 1 (1.7%) had Sertoli cell tumor, 1 (1.7%) had adenocarcinoma of appendix. Eleven (19%) tumors showed high CCNG1 expression, 44 (76%) tumors had medium expression, and 3 (5.2%) tumors had low CCNG1 expression. Of note, the tumor of one 14-year survivor with metastatic pancreatic adenocarcinoma in sustained remission had 24 % CCNG1 expression, one 3-year survivor with metastatic chondrosarcoma metastatic to lung with stable disease had 74% CCNG1 expression, one 2-year survivor with early stage HR+ HER2+ breast cancer in remission had 23% CCNG1 expression, and one 2-year survivor with early stage triple negative breast cancer in remission had 74% CCNG1 expression. Conclusions: Taken together, these data indicate that (1) Medium to high CCNG1 expression was found in 95% of tumors studied, and (2) Patients with medium CCNG1 expression who received DeltaRex-G had clinical benefit and are alive in sustained remission or with stable disease 2-14 years from DeltaRex-G treatment initiation, and (3) Prospective studies are warranted to correlate CCNG1 expression level and response to DeltaRex-G therapy.

Chromosome-level genome assemblies from two sandalwood species provide insights into the evolution of the Santalales

Sandalwood is one of the most expensive woods in the world and is well known for its long-lasting and distinctive aroma. In our study, chromosome-level genome assemblies for two sandalwood species (Santalum album and Santalumyasi) were constructed by integrating NGS short reads, RNA-seq, and Hi-C libraries with PacBio HiFi long reads. The S. album and S. yasi genomes were both assembled into 10 pseudochromosomes with a length of 229.59 Mb and 232.64 Mb, containing 21,673 and 22,816 predicted genes and a repeat content of 28.93% and 29.54% of the total genomes, respectively. Further analyses resolved a Santalum-specific whole-genome triplication event after divergence from ancestors of the Santalales lineage Malania, yet due to dramatic differences in transposon content, the Santalum genomes were only one-sixth the size of the Malania oleifera genome. Examination of RNA-seq data revealed a suite of genes that are differentially expressed in haustoria and might be involved in host hemiparasite interactions. The two genomes presented here not only provide an important comparative dataset for studying genome evolution in early diverging eudicots and hemiparasitic plants but will also hasten the application of conservation genomics for a lineage of trees recovering from decades of overexploitation.

The draft genome sequence of the Japanese rhinoceros beetle Trypoxylus dichotomus septentrionalis towards an understanding of horn formation

The Japanese rhinoceros beetle Trypoxylus dichotomus is a giant beetle with distinctive exaggerated horns present on the head and prothoracic regions of the male. T. dichotomus has been used as a research model in various fields such as evolutionary developmental biology, ecology, ethology, biomimetics, and drug discovery. In this study, de novo assembly of 615 Mb, representing 80% of the genome estimated by flow cytometry, was obtained using the 10 × Chromium platform. The scaffold N50 length of the genome assembly was 8.02 Mb, with repetitive elements predicted to comprise 49.5% of the assembly. In total, 23,987 protein-coding genes were predicted in the genome. In addition, de novo assembly of the mitochondrial genome yielded a contig of 20,217 bp. We also analyzed the transcriptome by generating 16 RNA-seq libraries from a variety of tissues of both sexes and developmental stages, which allowed us to identify 13 co-expressed gene modules. We focused on the genes related to horn formation and obtained new insights into the evolution of the gene repertoire and sexual dimorphism as exemplified by the sex-specific splicing pattern of the doublesex gene. This genomic information will be an excellent resource for further functional and evolutionary analyses, including the evolutionary origin and genetic regulation of beetle horns and the molecular mechanisms underlying sexual dimorphism.

Pervasive mRNA uridylation in fission yeast is catalysed by both Cid1 and Cid16 terminal uridyltransferases.

Messenger RNA uridylation is pervasive and conserved among eukaryotes, but the consequences of this modification for mRNA fate are still under debate. Utilising a simple model organism to study uridylation may facilitate efforts to understand the cellular function of this process. Here we demonstrate that uridylation can be detected using simple bioinformatics approach. We utilise it to unravel widespread transcript uridylation in fission yeast and demonstrate the contribution of both Cid1 and Cid16, the only two annotated terminal uridyltransferases (TUT-ases) in this yeast. To detect uridylation in transcriptome data, we used a RNA-sequencing (RNA-seq) library preparation protocol involving initial linker ligation to fragmented RNA-an approach borrowed from small RNA sequencing that was commonly used in older RNA-seq protocols. We next explored the data to detect uridylation marks. Our analysis show that uridylation in yeast is pervasive, similarly to the one in multicellular organisms. Importantly, our results confirm the role of the cytoplasmic uridyltransferase Cid1 as the primary uridylation catalyst. However, we also observed an auxiliary role of the second uridyltransferase, Cid16. Thus both fission yeast uridyltransferases are involved in mRNA uridylation. Intriguingly, we found no physiological phenotype of the single and double deletion mutants of cid1 and cid16 and only minimal impact of uridylation on steady-state mRNA levels. Our work establishes fission yeast as a potent model to study uridylation in a simple eukaryote, and we demonstrate that it is possible to detect uridylation marks in RNA-seq data without the need for specific methodologies.

RANK+TLR2+ myeloid subpopulation converts autoimmune to joint destruction in rheumatoid arthritis.

Joint destruction is the major clinic burden in patients with rheumatoid arthritis (RA). It is unclear, though, how this autoimmune disease progresses to the point of deterioration of the joint. Here, we report that in a mouse model of RA the upregulation of TLR2 expression and its α(2,3) sialylation in RANK+ myeloid monocytes mediate the transition from autoimmunity to osteoclast fusion and bone resorption, resulting in joint destruction. The expression of α(2,3) sialyltransferases was significantly increased in RANK+TLR2+ myeloid monocytes, and their inhibition or treatment with a TLR2 inhibitor blocked osteoclast fusion. Notably, analysis of our single-cell RNA-sequencing (scRNA-seq) libraries generated from RA mice revealed a novel RANK+TLR2- a subset that negatively regulated osteoclast fusion. Importantly, the RANK+TLR2+ subset was significantly diminished with the treatments, whereas the RANK+TLR2- subset was expanded. Moreover, the RANK+TLR2- subset could differentiate into a TRAP+ osteoclast lineage, but the resulting cells did not fuse to form osteoclasts. Our scRNA-seq data showed that Maf is highly expressed in the RANK+TLR2- subset, and the α(2,3) sialyltransferase inhibitor-induced Maf expression in the RANK+TLR2+ subset. The identification of a RANK+TLR2- subset provides a potential explanation for TRAP+ mononuclear cells in bone and their anabolic activity. Further, TLR2 expression and its α(2,3) sialylation in the RANK+ myeloid monocytes could be effective targets to prevent autoimmune-mediated joint destruction.