An electrochemical genomagnetic assay for detection of SARS-CoV-2 and Influenza A viruses in saliva.

Viral respiratory infections represent a major threat to the population's health globally. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes COVID-19 disease and in some cases the symptoms can be confused with Influenza disease caused by the Influenza A viruses. A simple, fast, and selective assay capable of identifying the etiological agent and differentiating the diseases is essential to provide the correct clinical management to the patient. Herein, we described the development of a genomagnetic assay for the selective capture of viral RNA from SARS-CoV-2 and Influenza A viruses in saliva samples and employing a simple disposable electrochemical device for gene detection and quantification. The proposed method showed excellent performance detecting RNA of SARS-CoV-2 and Influenza A viruses, with a limit of detection (LoD) and limit of quantification (LoQ) of 5.0 fmol L-1 and 8.6 fmol L-1 for SARS-CoV-2, and 1.0 fmol L-1 and 108.9 fmol L-1 for Influenza, respectively. The genomagnetic assay was employed to evaluate the presence of the viruses in 36 saliva samples and the results presented similar responses to those obtained by the real-time reverse transcription-polymerase chain reaction (RT-PCR), demonstrating the reliability and capability of a method as an alternative for the diagnosis of COVID-19 and Influenza with point-of-care capabilities.

Unveiling the role of MGMT and DAPK hypermethylation in response to anti-EGFR agents: Molecular insights for advancing HNSCC treatment.

BACKGROUND: Epidermal growth factor receptor (EGFR) is frequently activated in head and neck squamous cell carcinoma (HNSCC) and serves as a valuable target for therapy. Despite the availability of the EGFR inhibitors Cetuximab, Afatinib, and Allitinib, there are limited predictive markers for their response. Understanding molecular aberrations in HNSCC could facilitate the identification of new strategies for patient clinical and biological classification, offering novel therapeutic avenues. METHODS: We assessed CCNA1, DCC, MGMT, CDKN2A/p16, and DAPK methylation status in HNSCC cell lines and their association with anti-EGFR treatment response. RESULTS: MGMT methylation status displayed high sensitivity and specificity in distinguishing sensitive and resistant HNSCC cell lines to Afatinib (AUC = 0.955) and Allitinib (AUC = 0.935). Moreover, DAPK methylation status predicted response to Allitinib with high accuracy (AUC = 0.852), indicating their putative predictive biomarker roles. CONCLUSION: These findings hold promise for the development of more personalized and effective treatment approaches for HNSCC patients.

Biological and therapeutic implications of RKIP in Gastrointestinal Stromal Tumor (GIST): an integrated transcriptomic and proteomic analysis.

BACKGROUND: Gastrointestinal stromal tumors (GIST) represent a significant clinical challenge due to their metastatic potential and limited treatment options. Raf kinase inhibitor protein (RKIP), a suppressor of the MAPK signaling pathway, is downregulated in various cancers and acts as a metastasis suppressor. Our previous studies demonstrated low RKIP expression in GIST and its association with poor outcomes. This study aimed to expand on the previous findings and investigate the biological and therapeutic implications of RKIP loss on GIST. METHODS: To validate the RKIP prognostic significance, its expression was evaluated by immunohistochemistry in 142 bona fide GIST cases. The functional role of RKIP was evaluated in vitro, using the GIST-T1 cell line, which was knocked out for RKIP. The biological and therapeutic implications of RKIP were evaluated by invasion, migration, apoptosis, and 2D / 3D viability assays. Additionally, the transcriptome and proteome of RKIP knockout cells were determined by NanoString and mass spectrometry, respectively. RESULTS: Immunohistochemical analysis revealed the absence of RKIP in 25.3% of GIST cases, correlating with a tendency toward poor prognosis. Functional assays demonstrated that RKIP knockout increased GIST cells' invasion and migration potential by nearly 60%. Moreover, we found that RKIP knockout cells exhibited reduced responsiveness to Imatinib treatment and higher cellular viability in 2D and 3D in vitro models, as assessed by apoptosis-related protein expression. Through comprehensive genetic and proteomic profiling of RKIP knockout cells, we identified several putative RKIP-regulated proteins in GIST, such as COL3A1. CONCLUSIONS: Using a multidimensional integrative analysis, we identified, for the first time in GIST, molecules and pathways modulated by RKIP that may potentially drive metastasis and, consequently, poor prognosis in this disease.

miR-99a-5p as a biomarker for lymph node metastasis prediction in oral squamous cell carcinoma patients.

BACKGROUND: Metastatic lymph node involvement influences therapy decisions and serves as a prognostic indicator in oral squamous cell carcinoma (OSCC). However, many early-stage patients with clinically negative lymph nodes exhibit no metastasis upon surgical staging. This study aimed to identify differentially expressed miRNAs capable of distinguishing pathologically positive (pN+) from negative (pN0) nodes in OSCC patients without clinical evidence of lymph node metastases (cN0). METHODS: Expression levels of 798 miRNAs were assessed in tumor samples from 10 pN+ and 10 pN0 patients using the Nanostring nCounter platform. Validation was performed in an independent cohort of 15 pN+ and 24 pN0 patients through RT-qPCR. RESULTS: Eight miRNAs exhibited differential expression between pN0 and pN+ patients. Notably, hsa-miR-99a-5p demonstrated high sensitivity and specificity in predicting patients at higher risk of positive lymph nodes. CONCLUSIONS: These findings highlight hsa-miR-99a-5p as a potential biomarker for detecting lymph node metastasis in primary OSCC tumors.

RAD50 Deficient in a Breast Cancer Model Predicts Sensitivity to PARP Inhibitors.

BACKGROUND: Breast and ovarian tumors with pathogenic variants in BRCA1 or BRCA2 genes are more sensitive to poly (ADP-ribose) polymerase inhibitors (PARPi) treatment than wildtype tumors. Pathogenic variants in non-BRCA1/2 homologous recombination repair genes (HRR) also concede sensitivity to PARPi treatment. RAD50 participates in the Mre11-RAD50-Nbn (MRN) complex of the HRR pathway and plays an important role in DNA repair. OBJECTIVE: The objective of this study is to evaluate whether RAD50 protein deficiency modulates the PARPi response in breast cancer cell lines. METHODS: T47D breast cancer cell line was modified using small interfering RNA and CRISPR/Cas9 technology, to knockout the RAD50 gene. PARPi response (niraparib, olaparib and rucaparib alone or in combination with carboplatin), in T47D and T47D-edited clones, was evaluated by cell viability, cell cycle, apoptosis and protein expression analyses. RESULTS: Treatment with niraparib and carboplatin exerted a synergistic effect on T47D-RAD50 deficient cells and an antagonistic effect on T47D cells parental. Cell cycle analysis demonstrated an increase in the G2/M population in cells treated with niraparib or rucaparib alone or in combination with carboplatin. T47D-RAD50 deficient cells treated with rucaparib and carboplatin exhibited twofold levels in late apoptosis, also showing differences in PARP activation. All T47D RAD50 deficient clones treated with niraparib or rucaparib combined with carboplatin, or rucaparib alone showed increased levels of H2AX phosphorylation. CONCLUSIONS: T47D RAD50 deficient cells treated with PARP inhibitors alone or in combination with carboplatin showed cell cycle arrest in the G2/M phase, leading to death by apoptosis. Thus, RAD50 deficiency may be a good biomarker for predicting PARPi response.

The Barretos Cancer Hospital Animal Facility: Implementation and Results of a Dedicated Platform for Preclinical Oncology Models.

The Barretos Cancer Hospital Animal Facility (BCHAF) is a unique facility in Brazil exclusively dedicated to working with animal models for cancer research. In this article, we briefly present our modern facility and the main experiments performed, focusing on mutant strains of mice (PTCH-knockout and ApcMin mice), xenograft models, and patient-derived xenografts (PDXs). Our results show the progress and challenges in establishing these models and the need for having an appropriate representation of our cancer population to better understand tumor biology and to identify cancer biomarkers, which could be putatively targeted, allowing for personalized therapy.

Whole-Exome Sequencing Identifies Pathogenic Germline Variants in Patients with Lynch-Like Syndrome.

Lynch syndrome (LS) is the most common hereditary colorectal cancer (CRC) syndrome, characterized by germline pathogenic variants in mismatch repair (MMR)-related genes that lead to microsatellite instability. Patients who meet the clinical criteria for LS and MMR deficiency and without any identified germline pathogenic variants are frequently considered to have Lynch-like syndrome (LLS). These patients have a higher risk of CRC and extracolonic tumors, and little is known about their underlying genetic causes. We investigated the germline spectrum of LLS patients through whole-exome sequencing (WES). A total of 20 unrelated patients with MMR deficiency who met the clinical criteria for LS and had no germline variant were subjected to germline WES. Variant classification was performed according to the American College of Medical Genetics and Genomics (ACMG) criteria. Pathogenic/likely pathogenic variants were identified in 35% of patients in known cancer genes such as MUTYH and ATM. Besides this, rare and potentially pathogenic variants were identified in the DNA repair gene POLN and other cancer-related genes such as PPARG, CTC1, DCC and ALPK1. Our study demonstrates the germline mutational status of LLS patients, a population at high risk of colorectal cancer.

Efficacy of Combined Use of Everolimus and Second-Generation Pan-EGRF Inhibitors in KRAS Mutant Non-Small Cell Lung Cancer Cell Lines.

Background: EGFR mutations are present in approximately 15−50% of non-small cell lung cancer (NSCLC), which are predictive of anti-EGFR therapies. At variance, NSCLC patients harboring KRAS mutations are resistant to those anti-EGFR approaches. Afatinib and allitinib are second-generation pan-EGFR drugs, yet no predictive biomarkers are known in the NSCLC context. In the present study, we evaluated the efficacy of pan-EGFR inhibitors in a panel of 15 lung cancer cell lines associated with the KRAS mutations phenotype. Methods: KRAS wild-type sensitive NCI-H292 cell line was further transfected with KRAS mutations (p.G12D and p.G12S). The pan-EGFR inhibitors' activity and biologic effect of KRAS mutations were evaluated by cytotoxicity, MAPK phospho-protein array, colony formation, migration, invasion, and adhesion. In addition, in vivo chicken chorioallantoic membrane assay was performed in KRAS mutant cell lines. The gene expression profile was evaluated by NanoString. Lastly, everolimus and pan-EGFR combinations were performed to determine the combination index. Results: The GI50 score classified two cell lines treated with afatinib and seven treated with allitinib as high-sensitive phenotypes. All KRAS mutant cell lines demonstrated a resistant profile for both therapies (GI50 < 30%). The protein array of KRAS edited cells indicated a significant increase in AKT, CREB, HSP27, JNK, and, importantly, mTOR protein levels compared with KRAS wild-type cells. The colony formation, migration, invasion, adhesion, tumor perimeter, and mesenchymal phenotype were increased in the H292 KRAS mutated cells. Gene expression analysis showed 18 dysregulated genes associated with the focal adhesion-PI3K-Akt-mTOR-signaling correlated in KRAS mutant cell lines. Moreover, mTOR overexpression in KRAS mutant H292 cells was inhibited after everolimus exposure, and sensitivity to afatinib and allitinib was restored. Conclusions: Our results indicate that allitinib was more effective than afatinib in NSCLC cell lines. KRAS mutations increased aggressive behavior through upregulation of the focal adhesion-PI3K-Akt-mTOR-signaling in NSCLC cells. Significantly, everolimus restored sensibility and improved cytotoxicity of EGFR inhibitors in the KRAS mutant NSCLC cell lines.

New insights on familial colorectal cancer type X syndrome.

Familial colorectal cancer type X (FCCTX) is a heterogeneous colorectal cancer predisposition syndrome that, although displays a cancer pattern similar to Lynch syndrome, is mismatch repair proficient and does not exhibit microsatellite instability. Besides, its genetic etiology remains to be elucidated. In this study we performed germline exome sequencing of 39 cancer-affected patients from 34 families at risk for FCCTX. Variant classification followed the American College of Medical Genetics and Genomics (ACMG) guidelines. Pathogenic/likely pathogenic variants were identified in 17.65% of the families. Rare and potentially pathogenic alterations were identified in known hereditary cancer genes (CHEK2), in putative FCCTX candidate genes (OGG1 and FAN1) and in other cancer-related genes such as ATR, ASXL1, PARK2, SLX4 and TREX1. This study provides novel important clues that can contribute to the understanding of FCCTX genetic basis.

One-carbon metabolism and global DNA methylation in mothers of individuals with Down syndrome.

Down syndrome (DS) is the most common chromosomal disorder, resulting from the failure of normal chromosome 21 segregation. Studies have suggested that impairments within the one-carbon metabolic pathway can be of relevance for the global genome instability observed in mothers of individuals with DS. Based on the association between global DNA hypomethylation, genome instability, and impairments within the one-carbon metabolic pathway, the present study aimed to identify possible predictors, within the one-carbon metabolism, of global DNA methylation, measured by methylation patterns of LINE-1 and Alu repetitive sequences, in mothers of individuals with DS and mothers of individuals without the syndrome. In addition, we investigated one-carbon genetic polymorphisms and metabolites as maternal predisposing factors for the occurrence of trisomy 21 in children. Eighty-three samples of mothers of children with DS with karyotypically confirmed free trisomy 21 (case group) and 84 of mothers who had at least one child without DS or any other aneuploidy were included in the study. Pyrosequencing assays were performed to access global methylation. The results showed that group affiliation (case or control), betaine-homocysteine methyltransferase (BHMT) G742A and transcobalamin 2 (TCN2) C776G polymorphisms, and folate concentration were identified as predictors of global Alu DNA methylation values. In addition, thymidylate synthase (TYMS) 28-bp repeats 2R/3R or 3R/3R genotypes are independent maternal predisposing factors for having a child with DS. This study adds evidence that supports the association of impairments in the one-carbon metabolism, global DNA methylation, and the possibility of having a child with DS.

Evaluation of acetylation and methylation in oral rinse of patients with head and neck cancer history exposed to valproic acid.

Evaluate the biological action of valproic acid in the acetylation of histones and in the methylation of tumor suppressor genes via oral rinse in patients with a previous history of head and neck squamous cell carcinoma (HNSCC). Forty-two active or former smokers were included in this randomized, double-blind, placebo-controlled trial. Oral rinse samples were collected prior to treatment with valproic acid or placebo and after 90 days of treatment. The methylation status of five tumor suppressor genes and histone acetylation were evaluated by pyrosequencing and ELISA techniques, respectively. Differences between the 90-day and baseline oral rinse acetylation and methylation results were analyzed by comparing groups. Thirty-four patients were considered for analysis. The mean percentage adherence in the valproic and placebo groups was 93.4 and 93.0, respectively (p = 0.718). There was no statistically significant difference between groups when comparing the medians of the histone acetylation ratio and the methylation ratio for most of the studied genes. A significant reduction in the DCC methylation pattern was observed in the valproic group (p = 0.023). The use of valproic acid was safe and accompanied by good therapeutic adherence. DCC methylation was lower in the valproic acid group than in the placebo group.

DNA Methylation Markers from Negative Surgical Margins Can Predict Recurrence of Oral Squamous Cell Carcinoma.

The identification of molecular markers in negative surgical margins of oral squamous cell carcinoma (OSCC) might help in identifying residual molecular aberrations, and potentially improve the prediction of prognosis. We performed an Infinium MethylationEPIC BeadChip array on 32 negative surgical margins stratified based on the status of tumor recurrence in order to identify recurrence-specific aberrant DNA methylation (DNAme) markers. We identified 2512 recurrence-associated Differentially Methylated Positions (DMPs) and 392 Differentially Methylated Regions (DMRs) which were enriched in cell signaling and cancer-related pathways. A set of 14-CpG markers was able to discriminate recurrent and non-recurrent cases with high specificity and sensitivity rates (AUC 0.98, p = 3 × 10-6; CI: 0.95-1). A risk score based on the 14-CpG marker panel was applied, with cases classified within higher risk scores exhibiting poorer survival. The results were replicated using tumor-adjacent normal HNSCC samples from The Cancer Genome Atlas (TCGA). We identified residual DNAme aberrations in the negative surgical margins of OSCC patients, which could be informative for patient management by improving therapeutic intervention. This study proposes a novel DNAme-based 14-CpG marker panel as a promising predictor for tumor recurrence, which might contribute to improved decision-making for the personalized treatment of OSCC cases.

Electrochemical and optical detection and machine learning applied to images of genosensors for diagnosis of prostate cancer with the biomarker PCA3.

The development of simple detection methods aimed at widespread screening and testing is crucial for many infections and diseases, including prostate cancer where early diagnosis increases the chances of cure considerably. In this paper, we report on genosensors with different detection principles for a prostate cancer specific DNA sequence (PCA3). The genosensors were made with carbon printed electrodes or quartz coated with layer-by-layer (LbL) films containing gold nanoparticles and chondroitin sulfate and a layer of a complementary DNA sequence (PCA3 probe). The highest sensitivity was reached with electrochemical impedance spectroscopy with the detection limit of 83 pM in solutions of PCA3, while the limits of detection were 2000 pM and 900 pM for cyclic voltammetry and UV-vis spectroscopy, respectively. That detection could be performed with an optical method is encouraging, as one may envisage extending it to colorimetric tests. Since the morphology of sensing units is known to be affected in detection experiments, we applied machine learning algorithms to classify scanning electron microscopy images of the genosensors and managed to distinguish those exposed to PCA3-containing solutions from control measurements with an accuracy of 99.9%. The performance in distinguishing each individual PCA3 concentration in a multiclass task was lower, with an accuracy of 88.3%, which means that further developments in image analysis are required for this innovative approach.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland.

BACKGROUND: Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design. METHODS: RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment. RESULTS: We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin. CONCLUSION: It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.

Feasibility of methylated ctDNA detection in plasma samples of oropharyngeal squamous cell carcinoma patients.

BACKGROUND: Oropharyngeal squamous cell carcinomas (OpSCCs) are commonly associated with high rates of treatment failure. OBJECTIVES: To evaluate methylation-based markers in plasma from OpSCC patients as emerging tools for accurate/noninvasive follow-up. METHODS: Pretreatment formalin-fixed paraffin-embedded (FFPE) biopsies (n = 52) and paired plasma (n = 15) were tested for the methylation of CCNA1, DAPK, CDH8, and TIMP3 by droplet digital PCR (ddPCR). RESULTS: Seventy-one percent (37/52) of the biopsies showed methylation of at least one of the evaluated genes and tumor CCNA1 methylation was associated with recurrence-free survival. Methylated circulating tumor DNA (meth-ctDNA) was detected in 11/15 (73.3%) plasma samples; conversely, plasma samples from healthy controls were all negative for DNA methylation (area under the curve = 0.867; 95% confidence interval = 0.720-1.000). Additionally, preliminary results on the detection of meth-ctDNA in plasma collected during follow-up closely matched patient outcome. CONCLUSIONS: The results suggest the feasibility of detecting meth-ctDNA in plasma using ddPCR and a possible application on routine setting after further validation.

Immunosensors containing solution blow spun fibers of poly(lactic acid) to detect p53 biomarker.

This paper reports on biosensors made with a matrix of polylactic acid (PLA) fibers, which are suitable for immobilization of the anti-p53 active layer for detection of p53 biomarker. The PLA fibers were produced with solution blow spinning, a method that is advantageous for its simplicity and possibility to tune the fiber properties. For the biosensors, the optimized time to deposit the fibers was 60 s, with which detection of p53 could be achieved with the limit of detection of 11 pg/mL using electrical impedance spectroscopy. This sensitivity is also sufficient to detect the p53 biomarker in patient samples, which was confirmed by distinguishing samples from cell lines with distinct p53 concentrations in a plot where the impedance spectra were visualized with the interactive document mapping (IDMAP) technique. The high sensitivity and selectivity of the biosensors may be attributed to the specific interaction between the active layer and p53 modeled with a Langmuir-Freundlich and Freundlich isotherms and inferred from the analysis of the vibrational bands at 1550, 1650 and 1757 cm-1 using polarization-modulated infrared reflection absorption spectroscopy (PM-IRRAS). The successful immobilization of the active layer is evidence that the approach based on solution blown spun fibers may be replicated to other types of biosensors.

Impact of genetic variants in clinical outcome of a cohort of patients with oropharyngeal squamous cell carcinoma.

Tobacco- or human papillomavirus- driven oropharyngeal squamous cell carcinomas (OpSCC) represent distinct clinical, biological and epidemiological entities. The aim of this study was to identify genetic variants based on somatic alterations in OpSCC samples from an admixed population, and to test for association with clinical features. The entire coding region of 15 OpSCC driver genes was sequenced by next-generation sequencing in 51 OpSCC FFPE samples. Thirty-five percent of the patients (18/51) were HPV-positive and current or past tobacco consumption was reported in 86.3% (44/51). The mutation profile identified an average of 2.67 variants per sample. Sixty-three percent of patients (32/51; 62.7%) were mutated for at least one of the genes tested and TP53 was the most frequently mutated gene. The presence of mutation in NOTCH1 and PTEN, significantly decreased patient's recurrence-free survival, but only NOTCH1 mutation remained significant after stepwise selection, with a risk of recurrence of 4.5 (HR 95% CI = 1.11-14.57; Cox Regression p = 0.034). These results show that Brazilian OpSCC patients exhibit a similar clinical and genetic profile in comparison to other populations. Molecular characterization is a promising tool for the definition of clinical subgroups, aiding in a more precise tailoring of treatment and prognostication.

Genosensor made with a self-assembled monolayer matrix to detect MGMT gene methylation in head and neck cancer cell lines.

DNA methylation is involved in the oncogenesis of head and neck squamous cell carcinoma and could be used for early detection of cancer to increase the chances of cure, but unfortunately diagnosis is usually made at late stages of the disease. In this work we developed genosensors to detect DNA methylation of the MGMT gene in head and neck cancer cell lines. The probe for MGMT promoter methylation was immobilized on gold electrodes modified with 11-mercaptoundecanoic acid (11-MUA) self-assembled monolayers (SAM). Detection was performed with electrochemical impedance spectroscopy, with clear distinction between methylated and non-methylated DNA from head and neck cell lines. The genosensor is sensitive with a low detection limit of 0.24 × 10-12 mol L-1. In addition, the cell lines FaDu, JHU28 and SCC25 for the MGMT gene, could be distinguished from the HN13 cell line which has a high degree of MGMT methylation (97%), thus confirming the selectivity. Samples with different percentages of MGMT DNA methylation could be separated in multidimensional projections using the visualization technique interactive document mapping (IDMAP). The genosensor matrix and the immobilization procedures are generic, and can be extended to other DNA methylation biomarkers.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design. Methods: RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment. Results: We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin. Conclusion: It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.(AU)

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Background: Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design. Methods: RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment. Results: We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin. Conclusion: It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.