Background: Comprehensive annotation and quantification of transcriptomes are outstanding problems in functional genomics. While high throughput mRNA sequencing (RNA-Seq) has emerged as a powerful tool for addressing these problems, its success is dependent upon the availability and quality of reference genome sequences, thus limiting the organisms to which it can be applied. Results: Here, we describe Rnnotator, an automated software pipeline that generates transcript models by de novo assembly of RNA-Seq data without the need for a reference genome. We have applied the Rnnotator assembly pipeline to two yeast transcriptomes and compared the results to the reference gene catalogs of these organisms. The contigs produced by Rnnotator are highly accurate (95percent) and reconstruct full-length genes for the majority of the existing gene models (54.3percent). Furthermore, our analyses revealed many novel transcribed regions that are absent from well annotated genomes, suggesting Rnnotator serves as a complementary approach to analysis based on a reference genome for comprehensive transcriptomics. Conclusions: These results demonstrate that the Rnnotator pipeline is able to reconstruct full-length transcripts in the absence of a complete reference genome.

Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads

Ocimum L. of family Lamiaceae is a well known genus for its ethnobotanical, medicinal and aromatic properties, which are attributed to innumerable phenylpropanoid and terpenoid compounds produced by the plant. To enrich genomic resources for understanding various pathways, de novo transcriptome sequencing of two important species, O. sanctum and O. basilicum, was carried out by Illumina paired-end sequencing. The sequence assembly resulted in 69117 and 130043 transcripts with an average length of 1646 ± 1210.1 bp and 1363 ± 1139.3 bp for O. sanctum and O. basilicum, respectively. Out of the total transcripts, 59648 (86.30%) and 105470 (81.10%) from O. sanctum and O. basilicum, and respectively were annotated by uniprot blastx against Arabidopsis, rice and lamiaceae. KEGG analysis identified 501 and 952 transcripts from O. sanctum and O. basilicum, respectively, related to secondary metabolism with higher percentage of transcripts for biosynthesis of terpenoids in O. sanctum and phenylpropanoids in O. basilicum. Higher digital gene expression in O. basilicum was validated through qPCR and correlated to higher essential oil content and chromosome number (O. sanctum, 2n = 16; and O. basilicum, 2n = 48). Several CYP450 (26) and TF (40) families were identified having probable roles in primary and secondary metabolism. Also SSR and SNP markers were identified in the transcriptomes of both species with many SSRs linked to phenylpropanoid and terpenoid pathway genes. This is the first report of a comparative transcriptome analysis of Ocimum species and can be utilized to characterize genes related to secondary metabolism, their regulation, and breeding special chemotypes with unique essential oil composition in Ocimum.

/pdf/de-novo-sequencing-and-comparative-analysis-of-holy-and-2hbyukbcta.pdf

De novo sequencing and comparative analysis of holy and sweet basil transcriptomes

In plants, UGTs (UDP-glycosyltransferases) glycosylate various phytohormones and metabolites in response to biotic and abiotic stresses. Little is known about stress-responsive glycosyltransferases in plants. Therefore, it is important to understand the genomic and transcriptomic portfolio of plants with regard to biotic and abiotic stresses. Here, we identified 140, 154, and 251 putative UGTs in Brassica rapa, Brassica oleracea, and Brassica napus, respectively, and clustered them into 14 major phylogenetic groups (A–N). Fourteen major KEGG pathways and 24 biological processes were associated with the UGTs, highlighting them as unique modulators against environmental stimuli. Putative UGTs from B. rapa and B. oleracea showed a negative selection pressure and biased gene fractionation pattern during their evolution. Polyploidization increased the intron proportion and number of UGT-containing introns among Brassica. The putative UGTs were preferentially expressed in developing tissues and at the senescence stage. Differential expression of up- and down-regulated UGTs in response to phytohormone treatments, pathogen responsiveness and abiotic stresses, inferred from microarray and RNA-Seq data in Arabidopsis and Brassica broaden the glycosylation impact at the molecular level. This study identifies unique candidate UGTs for the manipulation of biotic and abiotic stress pathways in Brassica and Arabidopsis.

/pdf/comparative-genomic-and-transcriptomic-analyses-of-family-1-1i5vpoujj3.pdf

Comparative genomic and transcriptomic analyses of Family-1 UDP glycosyltransferase in three Brassica species and Arabidopsis indicates stress-responsive regulation

/pdf/the-p132h-mutation-in-the-main-protease-of-omicron-sars-cov-1qhx5oup.pdf

The P132H mutation in the main protease of Omicron SARS-CoV-2 decreases thermal stability without compromising catalysis or small-molecule drug inhibition

A large number of terpenoid compounds have been extracted from different tissues of S. officinalis. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 6.6 Gb of raw data were generated from the transcriptome of S. officinalis leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 21 million, comprising 98 million of high-quality nucleotide bases. 48,671 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 135 unigenes that are putatively involved in terpenoid metabolism, including 70 mevalonate and methyl-erythritol phosphate pathways, terpenoid backbone biosynthesis genes, and 65 terpene synthase genes. Moreover, five terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic tobacco; most transgenic tobacco plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. officinalis.

/pdf/transcriptome-and-metabolite-analyses-reveal-the-complex-15m0ciorr6.pdf

Transcriptome and metabolite analyses reveal the complex metabolic genes involved in volatile terpenoid biosynthesis in garden sage ( Salvia officinalis )

Background 
The cabbage, Brassica oleracea var. capitata L., has a distinguishable phenotype within the genus Brassica. Despite the economic and genetic importance of cabbage, there is little genomic data for cabbage, and most studies of Brassica are focused on other species or other B. oleracea subspecies. The lack of genomic data for cabbage, a non-model organism, hinders research on its molecular biology. Hence, the construction of reliable transcriptomic data based on high-throughput sequencing technologies is needed to enhance our understanding of cabbage and provide genomic information for future work.


Methodology/Principal Findings 
We constructed cDNAs from total RNA isolated from the roots, leaves, flowers, seedlings, and calcium-limited seedling tissues of two cabbage genotypes: 102043 and 107140. We sequenced a total of six different samples using the Illumina HiSeq platform, producing 40.5 Gbp of sequence data comprising 401,454,986 short reads. We assembled 205,046 transcripts (≥ 200 bp) using the Velvet and Oases assembler and predicted 53,562 loci from the transcripts. We annotated 35,274 of the loci with 55,916 plant peptides in the Phytozome database. The average length of the annotated loci was 1,419 bp. We confirmed the reliability of the sequencing assembly using reverse-transcriptase PCR to identify tissue-specific gene candidates among the annotated loci.


Conclusion 
Our study provides valuable transcriptome sequence data for B. oleracea var. capitata L., offering a new resource for studying B. oleracea and closely related species. Our transcriptomic sequences will enhance the quality of gene annotation and functional analysis of the cabbage genome and serve as a material basis for future genomic research on cabbage. The sequencing data from this study can be used to develop molecular markers and to identify the extreme differences among the phenotypes of different species in the genus Brassica.

/pdf/high-throughput-sequencing-and-de-novo-assembly-of-brassica-27wizpw5g7.pdf

High-Throughput Sequencing and De Novo Assembly of Brassica oleracea var. Capitata L. for Transcriptome Analysis

Abstract The small GTPase Cdc42 is an integral component of the cytoskeleton, and its dysregulation leads to pathophysiological conditions, such as cancer. Binding of Cdc42 to the scaffold protein IQGAP1 stabilizes Cdc42 in its active form. The interaction between Cdc42 and IQGAP1 enhances migration and invasion of cancer cells. Disrupting this association could impair neoplastic progression and metastasis; however, no effective means to achieve this has been described. Here, we screened 78,500 compounds using a homogeneous time resolved fluorescence-based assay to identify small molecules that disrupt the binding of Cdc42 to IQGAP1. From the combined results of the validation assay and counter-screens, we selected 44 potent compounds for cell-based experiments. Immunoprecipitation and cell viability analysis rendered four lead compounds, namely NCGC00131308, NCGC00098561, MLS000332963 and NCGC00138812, three of which inhibited proliferation and migration of breast carcinoma cells. Microscale thermophoresis revealed that two compounds bind directly to Cdc42. One compound reduced the amount of active Cdc42 in cells and effectively impaired filopodia formation. Docking analysis provided plausible models of the compounds binding to the hydrophobic pocket adjacent to the GTP binding site of Cdc42. In conclusion, we identified small molecules that inhibit binding between Cdc42 and IQGAP1, which could potentially yield chemotherapeutic agents. 

https://www.nature.com/articles/s41598-022-21342-w.pdf

Discovery of small molecule inhibitors that effectively disrupt IQGAP1-Cdc42 interaction in breast cancer cells

The COVID-19 pandemic has had enormous health, economic, and social consequences. Vaccines have been successful in reducing rates of infection and hospitalization, but there is still a need for an acute treatment for the disease. We investigate whether compounds that bind the human ACE2 protein can interrupt SARS-CoV-2 replication without damaging ACE2’s natural enzymatic function. Initial compounds were screened for binding to ACE2 but little interruption of ACE2 enzymatic activity. This set of compounds was extended by application of quantitative structure-activity analysis, which resulted in 512 virtual hits for further confirmatory screening. A subsequent SARS-CoV-2 replication assay revealed that five of these compounds inhibit SARS-CoV-2 replication in human cells. Further effort is required to completely determine the antiviral mechanism of these compounds, but they serve as a strong starting point for both development of acute treatments for COVID-19 and research into the mechanism of infection. TOC Graphic: Overall study design.

/pdf/allosteric-binders-of-ace2-are-promising-anti-sars-cov-2-2z6di0er.pdf

Allosteric binders of ACE2 are promising anti-SARS-CoV-2 agents

The COVID-19 pandemic has had enormous health, economic, and social consequences. Vaccines have been successful in reducing rates of infection and hospitalization, but there is still a need for acute treatment of the disease. We investigate whether compounds that bind the human angiotensin-converting enzyme 2 (ACE2) protein can decrease SARS-CoV-2 replication without impacting ACE2’s natural enzymatic function. Initial screening of a diversity library resulted in hit compounds active in an ACE2-binding assay, which showed little inhibition of ACE2 enzymatic activity (116 actives, success rate ∼4%), suggesting they were allosteric binders. Subsequent application of in silico techniques boosted success rates to ∼14% and resulted in 73 novel confirmed ACE2 binders with Kd values as low as 6 nM. A subsequent SARS-CoV-2 assay revealed that five of these compounds inhibit the viral life cycle in human cells. Further effort is required to completely elucidate the antiviral mechanism of these ACE2-binders, but they present a valuable starting point for both the development of acute treatments for COVID-19 and research into the host-directed therapy.

https://pubs.acs.org/doi/pdf/10.1021/acsptsci.2c00049

Jun Kyoung Choe

Papers

High-Throughput Sequencing and De Novo Assembly of Brassica oleracea var. Capitata L. for Transcriptome Analysis

Discovery of small molecule inhibitors that effectively disrupt IQGAP1-Cdc42 interaction in breast cancer cells

Allosteric binders of ACE2 are promising anti-SARS-CoV-2 agents

Allosteric Binders of ACE2 Are Promising Anti-SARS-CoV-2 Agents