Utility of CYP2D6 copy number variants as prognostic biomarker in localized anal squamous cell carcinoma.

BACKGROUND: Anal squamous cell carcinoma (ASCC) is an infrequent tumor whose treatment has not changed since the 1970s. The aim of this study is the identification of biomarkers allowing personalized treatments and improvement of therapeutic outcomes. METHODS: Forty-six paraffin tumor samples from ASCC patients were analyzed by whole-exome sequencing. Copy number variants (CNVs) were identified and their relation to disease-free survival (DFS) was studied and validated in an independent retrospective cohort of 101 ASCC patients from the Multidisciplinary Spanish Digestive Cancer Group (GEMCAD). GEMCAD cohort proteomics allowed assessing the biological features of these tumors. RESULTS: On the discovery cohort, the median age was 61 years old, 50% were males, stages I/II/III: 3 (7%)/16 (35%)/27 (58%), respectively, median DFS was 33 months, and overall survival was 45 months. Twenty-nine genes whose duplication was related to DFS were identified. The most representative was duplications of the CYP2D locus, including CYP2D6, CYP2D7P, and CYP2D8P genes. Patients with CYP2D6 CNV had worse DFS at 5 years than those with two CYP2D6 copies (21% vs. 84%; p < .0002, hazard ratio [HR], 5.8; 95% confidence interval [CI], 2.7-24.9). In the GEMCAD validation cohort, patients with CYP2D6 CNV also had worse DFS at 5 years (56% vs. 87%; p = .02, HR = 3.6; 95% CI, 1.1-5.7). Mitochondria and mitochondrial cell-cycle proteins were overexpressed in patients with CYP2D6 CNV. CONCLUSIONS: Tumor CYP2D6 CNV identified patients with a significantly worse DFS at 5 years among localized ASCC patients treated with 5-fluorouracil, mitomycin C, and radiotherapy. Proteomics pointed out mitochondria and mitochondrial cell-cycle genes as possible therapeutic targets for these high-risk patients. PLAIN LANGUAGE SUMMARY: Anal squamous cell carcinoma is an infrequent tumor whose treatment has not been changed since the 1970s. However, disease-free survival in late staged tumors is between 40% and 70%. The presence of an alteration in the number of copies of CYP2D6 gene is a biomarker of worse disease-free survival. The analysis of the proteins in these high-risk patients pointed out mitochondria and mitochondrial cell-cycle genes as possible therapeutic targets. Therefore, the determination of the number of copies of CYP2D6 allows the identification of anal squamous carcinoma patients with a high-risk of relapse that could be redirected to a clinical trial. Additionally, this study may be useful to suggest new treatment strategies to increase current therapy efficacy.

SeqEditor: an application for primer design and sequence analysis with or without GTF/GFF files.

MOTIVATION: Sequence analyses oriented to investigate specific features, patterns and functions of protein and DNA/RNA sequences usually require tools based on graphic interfaces whose main characteristic is their intuitiveness and interactivity with the user's expertise, especially when curation or primer design tasks are required. However, interface-based tools usually pose certain computational limitations when managing large sequences or complex datasets, such as genome and transcriptome assemblies. Having these requirments in mind we have developed SeqEditor an interactive software tool for nucleotide and protein sequences' analysis. RESULT: SeqEditor is a cross-platform desktop application for the analysis of nucleotide and protein sequences. It is managed through a Graphical User Interface and can work either as a graphical sequence browser or as a fasta task manager for multi-fasta files. SeqEditor has been optimized for the management of large sequences, such as contigs, scaffolds or even chromosomes, and includes a GTF/GFF viewer to visualize and manage annotation files. In turn, this allows for content mining from reference genomes and transcriptomes with similar efficiency to that of command line tools. SeqEditor also incorporates a set of tools for singleplex and multiplex PCR primer design and pooling that uses a newly optimized and validated search strategy for target and species-specific primers. All these features make SeqEditor a flexible application that can be used to analyses complex sequences, design primers in PCR assays oriented for diagnosis, and/or manage, edit and personalize reference sequence datasets. AVAILABILITYAND IMPLEMENTATION: SeqEditor was developed in Java using Eclipse Rich Client Platform and is publicly available at https://gpro.biotechvana.com/download/SeqEditor as binaries for Windows, Linux and Mac OS. The user manual and tutorials are available online at https://gpro.biotechvana.com/tool/seqeditor/manual. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genetic Profile and Functional Proteomics of Anal Squamous Cell Carcinoma: Proposal for a Molecular Classification.

Anal squamous cell carcinoma is a rare tumor. Chemo-radiotherapy yields a 50% 3-year relapse-free survival rate in advanced anal cancer, so improved predictive markers and therapeutic options are needed. High-throughput proteomics and whole-exome sequencing were performed in 46 paraffin samples from anal squamous cell carcinoma patients. Hierarchical clustering was used to establish groups de novo Then, probabilistic graphical models were used to study the differences between groups of patients at the biological process level. A molecular classification into two groups of patients was established, one group with increased expression of proteins related to adhesion, T lymphocytes and glycolysis; and the other group with increased expression of proteins related to translation and ribosomes. The functional analysis by the probabilistic graphical model showed that these two groups presented differences in metabolism, mitochondria, translation, splicing and adhesion processes. Additionally, these groups showed different frequencies of genetic variants in some genes, such as ATM, SLFN11 and DST Finally, genetic and proteomic characteristics of these groups suggested the use of some possible targeted therapies, such as PARP inhibitors or immunotherapy.

ICTV Virus Taxonomy Profile: Belpaoviridae 2021.

The family Belpaoviridae comprises metazoan-infecting reverse-transcribing viruses with long terminal repeats, commonly known as Bel/Pao LTR retrotransposons. These viruses share evolutionary history and genes involved in genome replication and virion formation with reverse-transcribing viruses of the families Metaviridae, Pseudoviridae, Retroviridae and Caulimoviridae. These five families form the order Ortervirales. This is a summary of the ICTV Report on the family Belpaoviridae, which is available at ictv.global/report/belpaoviridae.

ICTV Virus Taxonomy Profile: Pseudoviridae.

Pseudoviridae is a family of reverse-transcribing viruses with long terminal repeats (LTRs) belonging to the order Ortervirales. Pseudoviruses are commonly found integrated in the genomes of diverse plants, fungi and animals and are broadly known as Ty1/Copia LTR retrotransposons. Inside the cell, they form icosahedral virus particles, but unlike most other viruses, do not have an extracellular phase. This is a summary of the ICTV Report on the family Pseudoviridae, which is available at ictv.global/report/pseudoviridae.

CROSSMAPPER: estimating cross-mapping rates and optimizing experimental design in multi-species sequencing studies.

MOTIVATION: Numerous sequencing studies, including transcriptomics of host-pathogen systems, sequencing of hybrid genomes, xenografts, mixed species systems, metagenomics and meta-transcriptomics, involve samples containing genetic material from divergent organisms. A crucial step in these studies is identifying from which organism each sequencing read originated, and the experimental design should be directed to minimize biases caused by cross-mapping of reads to incorrect source genomes. Additionally, pooling of sufficiently different genetic material into a single sequencing library could significantly reduce experimental costs but requires careful planning and assessment of the impact of cross-mapping. Having these applications in mind we designed Crossmapper, the first to our knowledge tool able to assess cross-mapping prior to sequencing, therefore allowing optimization of experimental design. RESULTS: Using any combination of reference genomes, Crossmapper performs read simulation and back-mapping of those reads to the pool of references, quantifies and reports the cross-mapping rates for each organism. Crossmapper performs these analyses with numerous user-specified parameters, including, among others, read length, read layout, coverage, mapping parameters, genomic or transcriptomic data. Additionally, it outputs the results in highly interactive and publication-ready reports. This allows the user to perform multiple comparisons at once and choose the experimental setup minimizing cross-mapping rates. Moreover, Crossmapper can be used for resource optimization in sequencing facilities by pooling different samples into one sequencing library. AVAILABILITY AND IMPLEMENTATION: Crossmapper is a command line tool implemented in Python 3.6 and available as a conda package, allowing effortless installation. The source code, detailed information and a step-by-step tutorial is available at our GitHub page https://github.com/Gabaldonlab/crossmapper. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Immune expression profile identification in a group of proliferative verrucous leukoplakia patients: a pre-cancer niche for oral squamous cell carcinoma development.

OBJECTIVES: To explore the pathophysiology of proliferative verrucous leukoplakia, a rare oral disorder that exhibits high rates of recurrence and malignant transformation, through a RNAseq case-control study. MATERIAL AND METHODS: We obtained oral biopsies from 10 patients with verrucous leukoplakia lesions and from the mucosa of 5 healthy individuals for sequencing using RNAseq technology. Using bioinformatic methods, we investigated gene expression and enrichment differences between patients both with and without the disorder. We applied network biology methods to investigate functional relations among those genes that were differentially deregulated. RESULTS: We detected 140 differentially expressed genes with distinct roles in immune surveillance, tissue and organ morphogenesis, development, and organization. Of these 140 genes, 111 have been previously described as cancer expression biomarkers, being oral squamous cell carcinoma the most represented type of cancer among them. Of these 140 genes, 26 were prioritized for further investigation as biomarkers using larger sample sizes. CONCLUSIONS: The gene expression patterns of healthy and unhealthy patients differed in 140 genes whose deregulation has a functional impact on normal functioning of the immune system. This immune expression profile provides a plausible hypothesis to explain the transformation to oral squamous cell carcinoma observed in 6 of the 10 assayed cases. CLINICAL RELEVANCE: By determining the molecular bases of the proliferative verrucous leukoplakia disorder and identifying early biomarkers of malignancy, this can allow us to develop new treatment strategies.

Simulating the Influence of Conjugative-Plasmid Kinetic Values on the Multilevel Dynamics of Antimicrobial Resistance in a Membrane Computing Model.

Bacterial plasmids harboring antibiotic resistance genes are critical in the spread of antibiotic resistance. It is known that plasmids differ in their kinetic values, i.e., conjugation rate, segregation rate by copy number incompatibility with related plasmids, and rate of stochastic loss during replication. They also differ in cost to the cell in terms of reducing fitness and in the frequency of compensatory mutations compensating plasmid cost. However, we do not know how variation in these values influences the success of a plasmid and its resistance genes in complex ecosystems, such as the microbiota. Genes are in plasmids, plasmids are in cells, and cells are in bacterial populations and microbiotas, which are inside hosts, and hosts are in human communities at the hospital or the community under various levels of cross-colonization and antibiotic exposure. Differences in plasmid kinetics might have consequences on the global spread of antibiotic resistance. New membrane computing methods help to predict these consequences. In our simulation, conjugation frequency of at least 10-3 influences the dominance of a strain with a resistance plasmid. Coexistence of different antibiotic resistances occurs if host strains can maintain two copies of similar plasmids. Plasmid loss rates of 10-4 or 10-5 or plasmid fitness costs of ≥0.06 favor plasmids located in the most abundant species. The beneficial effect of compensatory mutations for plasmid fitness cost is proportional to this cost at high mutation frequencies (10-3 to 10-5). The results of this computational model clearly show how changes in plasmid kinetics can modify the entire population ecology of antibiotic resistance in the hospital setting.

ICTV Virus Taxonomy Profile: Metaviridae.

Metaviridae is a family of retrotransposons and reverse-transcribing viruses with long terminal repeats belonging to the order Ortervirales. Members of the genera Errantivirus and Metavirus include, respectively, Saccharomyces cerevisiae Ty3 virus and its Gypsy-like relatives in drosophilids. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Metaviridae, which is available at ictv.global/report/metaviridae.

nextPARS: parallel probing of RNA structures in Illumina.

RNA molecules play important roles in virtually every cellular process. These functions are often mediated through the adoption of specific structures that enable RNAs to interact with other molecules. Thus, determining the secondary structures of RNAs is central to understanding their function and evolution. In recent years several sequencing-based approaches have been developed that allow probing structural features of thousands of RNA molecules present in a sample. Here, we describe nextPARS, a novel Illumina-based implementation of in vitro parallel probing of RNA structures. Our approach achieves comparable accuracy to previous implementations, while enabling higher throughput and sample multiplexing.

Genome of wild olive and the evolution of oil biosynthesis.

Here we present the genome sequence and annotation of the wild olive tree (Olea europaea var. sylvestris), called oleaster, which is considered an ancestor of cultivated olive trees. More than 50,000 protein-coding genes were predicted, a majority of which could be anchored to 23 pseudochromosomes obtained through a newly constructed genetic map. The oleaster genome contains signatures of two Oleaceae lineage-specific paleopolyploidy events, dated at ∼28 and ∼59 Mya. These events contributed to the expansion and neofunctionalization of genes and gene families that play important roles in oil biosynthesis. The functional divergence of oil biosynthesis pathway genes, such as FAD2, SACPD, EAR, and ACPTE, following duplication, has been responsible for the differential accumulation of oleic and linoleic acids produced in olive compared with sesame, a closely related oil crop. Duplicated oleaster FAD2 genes are regulated by an siRNA derived from a transposable element-rich region, leading to suppressed levels of FAD2 gene expression. Additionally, neofunctionalization of members of the SACPD gene family has led to increased expression of SACPD2, 3, 5, and 7, consequently resulting in an increased desaturation of steric acid. Taken together, decreased FAD2 expression and increased SACPD expression likely explain the accumulation of exceptionally high levels of oleic acid in olive. The oleaster genome thus provides important insights into the evolution of oil biosynthesis and will be a valuable resource for oil crop genomics.

Increased miR-7641 Levels in Peritoneal Hyalinizing Vasculopathy in Long-Term Peritoneal Dialysis Patients.

Peritoneal hyalinizing vasculopathy (PHV) represents the cornerstone of long-term peritoneal dialysis (PD), and especially characterizes patients associated with encapsulating peritoneal sclerosis. However, the mechanisms of PHV development remain unknown. A cross sectional study was performed in 100 non-selected peritoneal biopsies of PD patients. Clinical data were collected and lesions were evaluated by immunohistochemistry. In selected biopsies a microRNA (miRNA)-sequencing analysis was performed. Only fifteen patients (15%) showed PHV at different degrees. PHV prevalence was significantly lower among patients using PD fluids containing low glucose degradation products (GDP) (5.9% vs. 24.5%), angiotensin converting enzyme inhibitors (ACEIs) (7.5% vs. 23.4%), statins (6.5% vs. 22.6%) or presenting residual renal function, suggesting the existence of several PHV protective factors. Peritoneal biopsies from PHV samples showed loss of endothelial markers and induction of mesenchymal proteins, associated with collagen IV accumulation and wide reduplication of the basement membrane. Moreover, co-expression of endothelial and mesenchymal markers, as well as TGF-ß1/Smad3 signaling activation were found in PHV biopsies. These findings suggest that an endothelial-to-mesenchymal transition (EndMT) process was taking place. Additionally, significantly higher levels of miR-7641 were observed in severe PHV compared to non-PHV peritoneal biopsies. Peritoneal damage by GDPs induce miRNA deregulation and an EndMT process in submesothelial vessels, which could contribute to collagen IV accumulation and PHV.

A metagenomic study of patients with alveolar osteitis after tooth extraction. A preliminary case-control study.

OBJECTIVE: To identify the microbiome in sockets with alveolar osteitis and compare it with a control group using metagenomic techniques. MATERIALS AND METHODS: A case-control study was conducted in subjects that had undergone a tooth extraction. Microbiological samples were taken from the sockets of 10 patients with dry socket after tooth extraction (AO group) and 10 patients in whom exodontia resulted in no postoperative complications (control group). Bacterial DNA was isolated, and the 16S rRNA gene was amplified and sequenced. Multiplexed tag-encoded sequencing of DNA from the samples was performed, and the reads were processed by Metagenomic Rapid Annotation. RESULTS: A total of 151 different species were found: 55 bacteria were only found in the AO group, 51 were specific to the control group, and 45 were common to both groups. The most frequently found genera in both groups were Prevotella. Prevotella nanceiensis, Actinomyces odontolyticus, Treponema maltophilum, Veillonella dispar, Tannerella forsythia, and Leuconostoc mesenteroides were found in several patients with alveolar osteitis, with an abundance greater than 0.5%, and were absent in all the control group samples. CONCLUSIONS: Patients who develop alveolar osteitis after dental extractions might have a different microbiota from that of patients without postoperative complications. Since this is a preliminary report, further research is needed to assess whether bacteria play an important role in the etiology of dry socket. CLINICAL RELEVANCE: This study seems to indicate that bacteria may play an important role in the alveolar osteitis etiology. Thus, new prevention and treatment strategies should be considered.

Differential distribution and enrichment of non-coding RNAs in exosomes from normal and Cancer-associated fibroblasts in colorectal cancer.

Exosome production from cancer-associated fibroblasts seems to be an important driver of tumor progression. We report the first in-depth biotype characterization of ncRNAs, analyzed by Next Generation Sequencing and Bioinformatics, expressed in established primary human normal and cancer-associated fibroblasts (CAFs) from cancer and normal mucosa tissues from 9 colorectal cancer patients, and/or packaged in their derived exosomes. Differential representation and enrichment analyses based on these ncRNAs revealed a significant number of differences between the ncRNA content of exosomes and the expression patterns of the normal and cancer-associated fibroblast cells. ncRNA regulatory elements are specifically packaged in CAF-derived exosomes, supporting a specific cross-talk between CAFs and colon cancer cells and/or other stromal cells, mediated by exosomes. These sncRNAs are potential biomarkers present in cancer-associated fibroblast-derived exosomes, which should thereby contribute to developing new non-invasive diagnostic, prognostic and predictive methods for clinical applications in management of cancer patients.

Analysis of the Transcriptome of the Red Seaweed Grateloupia imbricata with Emphasis on Reproductive Potential.

Grateloupia imbricata is an intertidal marine seaweed and candidate model organism for both industry and academic research, owing to its ability to produce raw materials such as carrageenan. Here we report on the transcriptome of G. imbricata with the aim of providing new insights into the metabolic pathways and other functional pathways related to the reproduction of Grateloupia species. Next-generation sequencing was carried out with subsequent de novo assembly and annotation using state-of-the-art bioinformatic protocols. The results show the presence of transcripts required for the uptake of glycerol, which is a specific carbon source for in vitro culture of G. imbricata and nucleotide sequences that are involved in polyamine-based biosynthesis, polyamine degradation, and metabolism of jasmonates and ethylene. Polyamines, ethylene and methyl jasmonate are plant growth regulators that elicit the development and maturation of cystocarps and the release of spores from seaweeds. Our results will inform studies of the mechanisms that control polysaccharide accumulation, cystocarp formation and spore release. Moreover, our transcriptome information clarifies aspects of red seaweed carposporogenesis with potential benefits for enhancing reproduction.

Acoustic Parametric Signal Generation for Underwater Communication.

This paper presents a study of different types of parametric signals with application to underwater acoustic communications. In all the signals, the carrier frequency is 200 kHz, which corresponds to the resonance frequency of the transducer under study and different modulations are presented and compared. In this sense, we study modulations with parametric sine sweeps (4 to 40 kHz) that represent binary codes (zeros and ones), getting closer to the application in acoustic communications. The different properties of the transmitting signals in terms of bit rate reconstruction, directivity, efficiency, and power needed are discussed as well.

Review of the Effects of Offshore Seismic Surveys in Cetaceans: Are Mass Strandings a Possibility?

Displacement of cetaceans is commonly reported during offshore seismic surveys. Speculation concerning possible links between seismic survey noise and cetacean strandings is available for a dozen events but without convincing causal evidence. This lack of evidence should not be considered conclusive but rather as reflecting the absence of a comprehensive analysis of the circumstances. Current mitigation guidelines are inadequate for long-range effects such as displacements and the potential for strandings. This review presents the available information for ten documented strandings that were possibly linked to seismic surveys and recommends initial measures and actions to further evaluate this potential link.

Transducer Development and Characterization for Underwater Acoustic Neutrino Detection Calibration.

A short bipolar pressure pulse with "pancake" directivity is produced and propagated when an Ultra-High Energy (UHE) neutrino interacts with a nucleus in water. Nowadays, acoustic sensor networks are being deployed in deep seas to detect this phenomenon as a first step toward building a neutrino telescope. In order to study the feasibility of the method, it is critical to have a calibrator that is able to mimic the neutrino signature. In previous works the possibility of using the acoustic parametric technique for this aim was proven. In this study, the array is operated at a high frequency and, by means of the parametric effect, the emission of the low-frequency acoustic bipolar pulse is generated mimicking the UHE neutrino acoustic pulse. To this end, the development of the transducer to be used in the parametric array is described in all its phases. The transducer design process, the characterization tests for the bare piezoelectric ceramic, and the addition of backing and matching layers are presented. The efficiencies and directivity patterns obtained for both primary and parametric beams confirm that the design of the proposed calibrator meets all the requirements for the emitter.

Membrane computing simulation of sexually transmitted bacterial infections in hotspots of individuals with various risk behaviors.

The epidemiology of sexually transmitted infections (STIs) is complex due to the coexistence of various pathogens, the variety of transmission modes derived from sexual orientations and behaviors at different ages and genders, and sexual contact hotspots resulting in network transmission. There is also a growing proportion of recreational drug users engaged in high-risk sexual activities, as well as pharmacological self-protection routines fostering non-condom practices. The frequency of asymptomatic patients makes it difficult to develop a comprehensive approach to STI epidemiology. Modeling approaches are required to deal with such complexity. Membrane computing is a natural computing methodology for the virtual reproduction of epidemics under the influence of deterministic and stochastic events with an unprecedented level of granularity. The application of the LOIMOS program to STI epidemiology illustrates the possibility of using it to shape appropriate interventions. Under the conditions of our basic landscape, including sexual hotspots of individuals with various risk behaviors, an increase in condom use reduces STIs in a larger proportion of heterosexuals than in same-gender sexual contacts and is much more efficient for reducing Neisseria gonorrhoeae than Chlamydia and lymphogranuloma venereum infections. Amelioration from diagnostic STI screening could be instrumental in reducing N. gonorrhoeae infections, particularly in men having sex with men (MSM), and Chlamydia trachomatis infections in the heterosexual population; however, screening was less effective in decreasing lymphogranuloma venereum infections in MSM. The influence of STI epidemiology of sexual contacts between different age groups (<35 and ≥35 years) and in bisexual populations was also submitted for simulation.IMPORTANCEThe epidemiology of sexually transmitted infections (STIs) is complex and significantly influences sexual and reproductive health worldwide. Gender, age, sexual orientation, sexual behavior (including recreational drug use and physical and pharmacological protection practices), the structure of sexual contact networks, and the limited application or efficiency of diagnostic screening procedures create variable landscapes in different countries. Modeling techniques are required to deal with such complexity. We propose the use of a simulation technology based on membrane computing, mimicking in silico STI epidemics under various local conditions with an unprecedented level of detail. This approach allows us to evaluate the relative weight of the various epidemic drivers in various populations at risk and the possible outcomes of interventions in particular epidemiological landscapes.