A straightforward assay to evaluate DNA integrity and optimize next-generation sequencing for clinical diagnosis in oncology.

Next generation sequencing (NGS) has become an informative tool to guide cancer treatment and conduce a personalized approach in oncology. The biopsy collected for pathologic analysis is usually stored as formalin-fixed paraffin-embedded (FFPE) blocks and then availed for molecular diagnostic, resulting in DNA molecules that are invariably fragmented and chemically modified. In an attempt to improve NGS based diagnostics in oncology we developed a straightforward DNA integrity assessment assay based on qPCR, defining clear parameters to whether NGS sequencing results is accurate or when it should be analyzed with caution. We performed DNA extraction from 12 tumor samples from diverse tissues and accessed DNA integrity by straightforward qPCR assays. In order to perform a cancer panel NGS sequencing, DNA library preparation was performed using RNA capture baits. Reads were aligned to the reference human genome and mutation calls were further validated by Sanger sequencing. Results obtained by the DNA integrity assays correlated to the efficiency of the pre-capture library preparation in up to 0.94 (Pearson's test). Moreover, sequencing results showed that poor integrity DNA leads to high rates of false positive mutation calls, specially C:G>T:A and C:G>A:T. Poor quality FFPE DNA samples are prone to generating false positive mutation calls. These are especially perilous in cases in which subclonal populations are expected, such as in advance disease, since it could lead clinicians to erroneous conclusions and equivocated conduct.

The Stress-Induced Soybean NAC Transcription Factor GmNAC81 Plays a Positive Role in Developmentally Programmed Leaf Senescence.

The onset of leaf senescence is a highly regulated developmental change that is controlled by both genetics and the environment. Senescence is triggered by massive transcriptional reprogramming, but functional information about its underlying regulatory mechanisms is limited. In the current investigation, we performed a functional analysis of the soybean (Glycine max) osmotic stress- and endoplasmic reticulum (ER) stress-induced NAC transcription factor GmNAC81 during natural leaf senescence using overexpression studies and reverse genetics. GmNAC81-overexpressing lines displayed accelerated flowering and leaf senescence but otherwise developed normally. The precocious leaf senescence of GmNAC81-overexpressing lines was associated with greater Chl loss, faster photosynthetic decay and higher expression of hydrolytic enzyme-encoding GmNAC81 target genes, including the vacuolar processing enzyme (VPE), an executioner of vacuole-triggered programmed cell death (PCD). Conversely, virus-induced gene silencing-mediated silencing of GmNAC81 delayed leaf senescence and was associated with reductions in Chl loss, lipid peroxidation and the expression of GmNAC81 direct targets. Promoter-reporter studies revealed that the expression pattern of GmNAC81 was associated with senescence in soybean leaves. Our data indicate that GmNAC81 is a positive regulator of age-dependent senescence and may integrate osmotic stress- and ER stress-induced PCD responses with natural leaf senescence through the GmNAC81/VPE regulatory circuit.

Functional and regulatory conservation of the soybean ER stress-induced DCD/NRP-mediated cell death signaling in plants.

BACKGROUND: The developmental and cell death domain (DCD)-containing asparagine-rich proteins (NRPs) were first identified in soybean (Glycine max) as transducers of a cell death signal derived from prolonged endoplasmic reticulum (ER) stress, osmotic stress, drought or developmentally-programmed leaf senescence via the GmNAC81/GmNAC30/GmVPE signaling module. In spite of the relevance of the DCD/NRP-mediated signaling as a versatile adaptive response to multiple stresses, mechanistic knowledge of the pathway is lacking and the extent to which this pathway may operate in the plant kingdom has not been investigated. RESULTS: Here, we demonstrated that the DCD/NRP-mediated signaling also propagates a stress-induced cell death signal in other plant species with features of a programmed cell death (PCD) response. In silico analysis revealed that several plant genomes harbor conserved sequences of the pathway components, which share functional analogy with their soybean counterparts. We showed that GmNRPs, GmNAC81and VPE orthologs from Arabidopsis, designated as AtNRP-1, AtNRP-2, ANAC036 and gVPE, respectively, induced cell death when transiently expressed in N. benthamiana leaves. In addition, loss of AtNRP1 and AtNRP2 function attenuated ER stress-induced cell death in Arabidopsis, which was in marked contrast with the enhanced cell death phenotype displayed by overexpressing lines as compared to Col-0. Furthermore, atnrp-1 knockout mutants displayed enhanced sensitivity to PEG-induced osmotic stress, a phenotype that could be complemented with ectopic expression of either GmNRP-A or GmNRP-B. In addition, AtNRPs, ANAC036 and gVPE were induced by osmotic and ER stress to an extent that was modulated by the ER-resident molecular chaperone binding protein (BiP) similarly as in soybean. Finally, as putative downstream components of the NRP-mediated cell death signaling, the stress induction of AtNRP2, ANAC036 and gVPE was dependent on the AtNRP1 function. BiP overexpression also conferred tolerance to water stress in Arabidopsis, most likely due to modulation of the drought-induced NRP-mediated cell death response. CONCLUSION: Our results indicated that the NRP-mediated cell death signaling operates in the plant kingdom with conserved regulatory mechanisms and hence may be target for engineering stress tolerance and adaptation in crops.

ICRmax: an optimized approach to detect tumor-specific interchromosomal rearrangements for clinical application.

Somatically acquired chromosomal rearrangements occur at early stages during tumorigenesis and can be used to indirectly detect tumor cells, serving as highly sensitive and tumor-specific biomarkers. Advances in high-throughput sequencing have allowed the genome-wide identification of patient-specific chromosomal rearrangements to be used as personalized biomarkers to efficiently assess response to treatment, detect residual disease and monitor disease recurrence. However, sequencing and data processing costs still represent major obstacles for the widespread application of personalized biomarkers in oncology. We developed a computational pipeline (ICRmax) for the cost-effective identification of a minimal set of tumor-specific interchromosomal rearrangements (ICRs). We examined ICRmax performance on sequencing data from rectal tumors and simulated data achieving an average accuracy of 68% for ICR identification. ICRmax identifies ICRs from low-coverage sequenced tumors, eliminates the need to sequence a matched normal tissue and significantly reduces the costs that limit the utilization of personalized biomarkers in the clinical setting.

Exploring the versatility of sesquiterpene biosynthesis in guava plants: a comparative genome-wide analysis of two cultivars.

Psidium guajava L., a fruit crop belonging to the Myrtaceae family, is highly valued for its nutritional and medicinal properties. The family exhibits a diverse chemical profile of essential oils and serves as a valuable resource due to its ecological interactions, adaptability, and dispersal capacity. The Myrtaceae family has been extensively studied for its terpenoids. Genetic studies have focused on foliar terpene yield in species from the Eucalypteae and Melaleucaceae tribes. To understand the evolutionary trends in guava breeding, this study predicted terpene synthase genes (TPS) from different cultivars. Through this analysis, 43 full-length TPS genes were identified, and approximately 77% of them exhibited relative expression in at least one of the five investigated plant tissues (root, leaf, bud, flower, and fruit) of two guava cultivars. We identified intra-species variation in the terpene profile and single nucleotide polymorphisms (SNPs) in twelve TPS genes, resulting in the clustering of 62 genotypes according to their essential oil chemotypes. The high concentration of sesquiterpenes is supported by the higher number of TPS-a genes and their expression. The expansion for TPS sub-families in P. guajava occurred after the expansion of other rosids species. Providing insight into the origin of structural diversification and expansion in each clade of the TPS gene family within Myrtaceae. This study can provide insights into the diversity of genes for specialized metabolites such as terpenes, and their regulation, which can lead to a diverse chemotype of essential oil in different tissues and genotypes. This suggests a mode of enzymatic evolution that could lead to high sesquiterpene production, act as a chemical defense and contribute to the adaptive capacity of this species to different habitats.

Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers.

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Abiotic Stresses in Plants and Their Markers: A Practice View of Plant Stress Responses and Programmed Cell Death Mechanisms.

Understanding how plants cope with stress and the intricate mechanisms thereby used to adapt and survive environmental imbalances comprise one of the most powerful tools for modern agriculture. Interdisciplinary studies suggest that knowledge in how plants perceive, transduce and respond to abiotic stresses are a meaningful way to design engineered crops since the manipulation of basic characteristics leads to physiological remodeling for plant adaption to different environments. Herein, we discussed the main pathways involved in stress-sensing, signal transduction and plant adaption, highlighting biochemical, physiological and genetic events involved in abiotic stress responses. Finally, we have proposed a list of practice markers for studying plant responses to multiple stresses, highlighting how plant molecular biology, phenotyping and genetic engineering interconnect for creating superior crops.

Efficient method for isolation of high-quality RNA from Psidium guajava L. tissues.

Acquiring high-quality RNA in sufficient amounts is crucial in plant molecular biology and genetic studies. Several methods for RNA extraction from plants are available in the literature, mainly due to the great biochemical diversity present in each species and tissue, which can complicate or prevent the extraction. Psidium guajava (Myrtaceae family) is a perennial fruit tree of medicinal and economic value; nevertheless, only a few molecular studies are available for the species. One reason is the difficulty in obtaining RNA due to the content of the samples, which are rich in polyphenols, polysaccharides, and secondary metabolites. Furthermore, there are few studies available for the isolation of RNA from guava or Psidium samples, which hampers advances in the study of the genus. Here, quality and yields of RNA isolates were compared using six extraction protocols: two protocols based on the application of cetyltrimethylammonium bromide (CTAB) lysis buffer, one protocol which uses the TRIzol reagent, one which applies guanidine thiocyanate lysis buffer followed by organic phase extraction, and two commercial kits (PureLink RNA Mini Kit and RNeasy Plant Mini Kit). The CTAB-based method provided the highest RNA yields and quality for five different tissues (flower bud, immature leaf, young leaf, mature leaf, and root), genotypes, and stress conditions. For the most efficient protocol, the average yield of RNA from guava leaves was 203.06 µg/g of tissue, and the A260/A280 and A260/A230 ratios were 2.1 and 2.2, respectively. RT-qPCR analysis demonstrated that the purity of the samples was sufficient for molecular biology experiments. CTAB-based methods for RNA isolation were found to be the most efficient, providing the highest RNA yields and quality for tissues from P. guajava. Additionally, they were compatible for downstream RNA-based applications, besides being simple and cost-effective.

Remodeling of the cell wall as a drought-tolerance mechanism of a soybean genotype revealed by global gene expression analysis.

Drought stress is major abiotic stress that affects soybean production. Therefore, it is widely desirable that soybean becomes more tolerant to stress. To provide insights into regulatory mechanisms of the stress response, we compared the global gene expression profiles from leaves of two soybean genotypes that display different responses to water-deficit (BR 16 and Embrapa 48, drought-sensitive and drought-tolerant, respectively). After the RNA-seq analysis, a total of 5335 down-regulated and 3170 up-regulated genes were identified in the BR16. On the other hand, the number of genes differentially expressed was markedly lower in the Embrapa 48, 355 up-regulated and 471 down-regulated genes. However, induction and expression of protein kinases and transcription factors indicated signaling cascades involved in the drought tolerance. Overall, the results suggest that the metabolism of pectin is differently modulated in response to drought stress and may play a role in the soybean defense mechanism against drought. This occurs via an increase of the cell wall plasticity and crosslink, which contributed to a higher hydraulic conductance (K f) and relative water content (RWC%). The drought-tolerance mechanism of the Embrapa 48 genotype involves remodeling of the cell wall and increase of the hydraulic conductance to the maintenance of cell turgor and metabolic processes, resulting in the highest leaf RWC, photosynthetic rate (A), transpiration (E) and carboxylation (A/C i). Thus, we concluded that the cell wall adjustment under drought is important for a more efficient water use which promoted a more active photosynthetic metabolism, maintaining higher plant growth under drought stress. Supplementary Information: The online version contains supplementary material available at 10.1007/s42994-021-00043-4.

GmNAC81 Inversely Modulates Leaf Senescence and Drought Tolerance.

Glycine max NAC81 (GmNAC81) is a downstream effector of the DCD/NRP-mediated cell death signaling, which interacts with GmNAC30 to fully induce the caspase 1-like vacuolar processing enzyme (VPE) expression, the executioner of the cell death program. GmNAC81 has been previously shown to positively modulate leaf senescence via the NRP/GmNAC81/VPE signaling module. Here, we examined the transcriptome induced by GmNAC81 overexpression and leaf senescence and showed that GmNAC81 further modulates leaf senescence by regulating an extensive repertoire of functionally characterized senescence-associated genes (SAGs). Because the NRP/GmNAC81/VPE signaling circuit also relays stress-induced cell death signals, we examined the effect of GmNAC81 overexpression in drought responses. Enhanced GmNAC81 expression in the transgenic lines increased sensitivity to water deprivation. Under progressive drought, the GmNAC81-overexpressing lines displayed severe leaf wilting, a larger and faster decline in leaf Ψw, relative water content (RWC), photosynthesis rate, stomatal conductance, and transpiration rate, in addition to higher Ci/Ca and lower Fm/Fv ratios compared to the BR16 control line. Collectively, these results indicate that the photosynthetic activity and apparatus were more affected by drought in the transgenic lines. Consistent with hypersensitivity to drought, chlorophyll loss, and lipid peroxidation were higher in the GmNAC81-overexpressing lines than in BR16 under dehydration. In addition to inducing VPE expression, GmNAC81 overexpression uncovered the regulation of typical drought-responsive genes. In particular, key regulators and effectors of ABA signaling were suppressed by GmNAC81 overexpression. These results suggest that GmNAC81 may negatively control drought tolerance not only via VPE activation but also via suppression of ABA signaling.

The receptor-like kinase NIK1 targets FLS2/BAK1 immune complex and inversely modulates antiviral and antibacterial immunity.

Plants deploy various immune receptors to recognize pathogens and defend themselves. Crosstalk may happen among receptor-mediated signal transduction pathways in the same host during simultaneous infection of different pathogens. However, the related function of the receptor-like kinases (RLKs) in thwarting different pathogens remains elusive. Here, we report that NIK1, which positively regulates plant antiviral immunity, acts as an important negative regulator of antibacterial immunity. nik1 plants exhibit dwarfed morphology, enhanced disease resistance to bacteria and increased PAMP-triggered immunity (PTI) responses, which are restored by NIK1 reintroduction. Additionally, NIK1 negatively regulates the formation of the FLS2/BAK1 complex. The interaction between NIK1 and FLS2/BAK1 is enhanced upon flg22 perception, revealing a novel PTI regulatory mechanism by an RLK. Furthermore, flg22 perception induces NIK1 and RPL10A phosphorylation in vivo, activating antiviral signalling. The NIK1-mediated inverse modulation of antiviral and antibacterial immunity may allow bacteria and viruses to activate host immune responses against each other.

Revisiting the Soybean GmNAC Superfamily.

The NAC (NAM, ATAF, and CUC) genes encode transcription factors involved with the control of plant morph-physiology and stress responses. The release of the last soybean (Glycine max) genome assembly (Wm82.a2.v1) raised the possibility that new NAC genes would be present in the soybean genome. Here, we interrogated the last version of the soybean genome against a conserved NAC domain structure. Our analysis identified 32 putative novel NAC genes, updating the superfamily to 180 gene members. We also organized the genes in 15 phylogenetic subfamilies, which showed a perfect correlation among sequence conservation, expression profile, and function of orthologous Arabidopsis thaliana genes and NAC soybean genes. To validate our in silico analyses, we monitored the stress-mediated gene expression profiles of eight new NAC-genes by qRT-PCR and monitored the GmNAC senescence-associated genes by RNA-seq. Among ER stress, osmotic stress and salicylic acid treatment, all the novel tested GmNAC genes responded to at least one type of stress, displaying a complex expression profile under different kinetics and extension of the response. Furthermore, we showed that 40% of the GmNACs were differentially regulated by natural leaf senescence, including eight (8) newly identified GmNACs. The developmental and stress-responsive expression profiles of the novel NAC genes fitted perfectly with their phylogenetic subfamily. Finally, we examined two uncharacterized senescence-associated proteins, GmNAC065 and GmNAC085, and a novel, previously unidentified, NAC protein, GmNAC177, and showed that they are nuclear localized, and except for GmNAC065, they display transactivation activity in yeast. Consistent with a role in leaf senescence, transient expression of GmNAC065 and GmNAC085 induces the appearance of hallmarks of leaf senescence, including chlorophyll loss, leaf yellowing, lipid peroxidation and accumulation of H2O2. GmNAC177 was clustered to an uncharacterized subfamily but in close proximity to the TIP subfamily. Accordingly, it was rapidly induced by ER stress and by salicylic acid under late kinetic response and promoted cell death in planta. Collectively, our data further substantiated the notion that the GmNAC genes display functional and expression profiles consistent with their phylogenetic relatedness and established a complete framework of the soybean NAC superfamily as a foundation for future analyses.

E2F1 somatic mutation within miRNA target site impairs gene regulation in colorectal cancer.

BACKGROUND: Genetic studies have largely concentrated on the impact of somatic mutations found in coding regions, and have neglected mutations outside of these. However, 3' untranslated regions (3' UTR) mutations can also disrupt or create miRNA target sites, and trigger oncogene activation or tumor suppressor inactivation. METHODS: We used next-generation sequencing to widely screen for genetic alterations within predicted miRNA target sites of oncogenes associated with colorectal cancer, and evaluated the functional impact of a new somatic mutation. Target sequencing of 47 genes was performed for 29 primary colorectal tumor samples. For 71 independent samples, Sanger methodology was used to screen for E2F1 mutations in miRNA predicted target sites, and the functional impact of these mutations was evaluated by luciferase reporter assays. RESULTS: We identified germline and somatic alterations in E2F1. Of the 100 samples evaluated, 3 had germline alterations at the MIR205-5p target site, while one had a somatic mutation at MIR136-5p target site. E2F1 gene expression was similar between normal and tumor tissues bearing the germline alteration; however, expression was increased 4-fold in tumor tissue that harbored a somatic mutation compared to that in normal tissue. Luciferase reporter assays revealed both germline and somatic alterations increased E2F1 activity relative to wild-type E2F1. CONCLUSIONS: We demonstrated that somatic mutation within E2F1:MIR136-5p target site impairs miRNA-mediated regulation and leads to increased gene activity. We conclude that somatic mutations that disrupt miRNA target sites have the potential to impact gene regulation, highlighting an important mechanism of oncogene activation.

Tissue-Associated Bacterial Alterations in Rectal Carcinoma Patients Revealed by 16S rRNA Community Profiling.

Sporadic and inflammatory forms of colorectal cancer (CRC) account for more than 80% of cases. Recent publications have shown mechanistic evidence for the involvement of gut bacteria in the development of both CRC-forms. Whereas, colon and rectal cancer have been routinely studied together as CRC, increasing evidence show these to be distinct diseases. Also, the common use of fecal samples to study microbial communities may reflect disease state but possibly not the tumor microenvironment. We performed this study to evaluate differences in bacterial communities found in tissue samples of 18 rectal-cancer subjects when compared to 18 non-cancer controls. Samples were collected during exploratory colonoscopy (non-cancer group) or during surgery for tumor excision (rectal-cancer group). High throughput 16S rRNA amplicon sequencing of the V4-V5 region was conducted on the Ion PGM platform, reads were filtered using Qiime and clustered using UPARSE. We observed significant increases in species richness and diversity in rectal cancer samples, evidenced by the total number of OTUs and the Shannon and Simpson indexes. Enterotyping analysis divided our cohort into two groups, with the majority of rectal cancer samples clustering into one enterotype, characterized by a greater abundance of Bacteroides and Dorea. At the phylum level, rectal-cancer samples had increased abundance of candidate phylum OD1 (also known as Parcubacteria) whilst non-cancer samples had increased abundance of Planctomycetes. At the genera level, rectal-cancer samples had higher abundances of Bacteroides, Phascolarctobacterium, Parabacteroides, Desulfovibrio, and Odoribacter whereas non-cancer samples had higher abundances of Pseudomonas, Escherichia, Acinetobacter, Lactobacillus, and Bacillus. Two Bacteroides fragilis OTUs were more abundant among rectal-cancer patients seen through 16S rRNA amplicon sequencing, whose presence was confirmed by immunohistochemistry and enrichment verified by digital droplet PCR. Our findings point to increased bacterial richness and diversity in rectal cancer, along with several differences in microbial community composition. Our work is the first to present evidence for a possible role of bacteria such as B. fragilis and the phylum Parcubacteria in rectal cancer, emphasizing the need to study tissue-associated bacteria and specific regions of the gastrointestinal tract in order to better understand the possible links between the microbiota and rectal cancer.

The use of personalized biomarkers and liquid biopsies to monitor treatment response and disease recurrence in locally advanced rectal cancer after neoadjuvant chemoradiation.

Neoadjuvant chemoradiotherapy (nCRT) followed by surgery is the mainstay treatment for locally advanced rectal cancer. Variable degrees of tumor regression are observed after nCRT and alternative treatment strategies, including close surveillance without immediate surgery, have been investigated to spare patients with complete tumor regression from potentially adverse outcomes of radical surgery. However, clinical and radiological assessment of response does not allow accurate identification of patients with complete response. In addition, surveillance for recurrence is similarly important for these patients, as early detection of recurrence allows salvage resections and adjuvant interventions. We report the use of liquid biopsies and personalized biomarkers for monitoring treatment response to nCRT and detecting residual disease and recurrence in patients with rectal cancer. We sequenced the whole-genome of four rectal tumors to identify patient-specific chromosomal rearrangements that were used to monitor circulating tumor DNA (ctDNA) in liquid biopsies collected at diagnosis and during nCRT and follow-up. We compared ctDNA levels to clinical, radiological and pathological response to nCRT. Our results indicate that personalized biomarkers and liquid biopsies may not be sensitive for the detection of microscopic residual disease. However, it can be efficiently used to monitor treatment response to nCRT and detect disease recurrence, preceding increases in CEA levels and radiological diagnosis. Similar good results were observed when assessing tumor response to systemic therapy and disease progression. Our study supports the use of personalized biomarkers and liquid biopsies to tailor the management of rectal cancer patients, however, replication in a larger cohort is necessary to introduce this strategy into clinical practice.