LncRNAs of Saccharomyces cerevisiae bypass the cell cycle arrest imposed by ethanol stress.

Ethanol alters many subsystems of Saccharomyces cerevisiae, including the cell cycle. Two ethanol-responsive lncRNAs in yeast interact with cell cycle proteins, and here, we investigated the role of these RNAs in cell cycle. Our network dynamic modeling showed that higher and lower ethanol-tolerant strains undergo cell cycle arrest in mitosis and G1 phases, respectively, during ethanol stress. The higher population rebound of the lower ethanol-tolerant phenotype after stress relief responds to the late phase arrest. We found that the lncRNA lnc9136 of SEY6210 (a lower ethanol-tolerant strain) induces cells to skip mitosis arrest. Simulating an overexpression of lnc9136 and analyzing CRISPR-Cas9 mutants lacking this lncRNA suggest that lnc9136 induces a regular cell cycle even under ethanol stress, indirectly regulating Swe1p and Clb1/2 by binding to Gin4p and Hsl1p. Notably, lnc10883 of BY4742 (a higher ethanol-tolerant strain) does not prevent G1 arrest in this strain under ethanol stress. However, lnc19883 circumvents DNA and spindle damage checkpoints, maintaining a functional cell cycle by interacting with Mec1p or Bub1p even in the presence of DNA/spindle damage. Overall, we present the first evidence of direct roles for lncRNAs in regulating yeast cell cycle proteins, the dynamics of this system in different ethanol-tolerant phenotypes, and a new yeast cell cycle model.

Integrative Analysis of the Ethanol Tolerance of Saccharomyces cerevisiae.

Ethanol (EtOH) alters many cellular processes in yeast. An integrated view of different EtOH-tolerant phenotypes and their long noncoding RNAs (lncRNAs) is not yet available. Here, large-scale data integration showed the core EtOH-responsive pathways, lncRNAs, and triggers of higher (HT) and lower (LT) EtOH-tolerant phenotypes. LncRNAs act in a strain-specific manner in the EtOH stress response. Network and omics analyses revealed that cells prepare for stress relief by favoring activation of life-essential systems. Therefore, longevity, peroxisomal, energy, lipid, and RNA/protein metabolisms are the core processes that drive EtOH tolerance. By integrating omics, network analysis, and several other experiments, we showed how the HT and LT phenotypes may arise: (1) the divergence occurs after cell signaling reaches the longevity and peroxisomal pathways, with CTA1 and ROS playing key roles; (2) signals reaching essential ribosomal and RNA pathways via SUI2 enhance the divergence; (3) specific lipid metabolism pathways also act on phenotype-specific profiles; (4) HTs take greater advantage of degradation and membraneless structures to cope with EtOH stress; and (5) our EtOH stress-buffering model suggests that diauxic shift drives EtOH buffering through an energy burst, mainly in HTs. Finally, critical genes, pathways, and the first models including lncRNAs to describe nuances of EtOH tolerance are reported here.

Uncovering patterns of the evolution of genomic sequence entropy and complexity.

The lack of consensus concerning the biological meaning of entropy and complexity of genomes and the different ways to assess these data hamper conclusions concerning what are the causes of genomic entropy variation among species. This study aims to evaluate the entropy and complexity of genomic sequences of several species without using homologies to assess relationships among these variables and non-molecular data (e.g., the number of individuals) to seek a trigger of interspecific genomic entropy variation. The results indicate a relationship among genomic entropy, genome size, genomic complexity, and the number of individuals: species with a small number of individuals harbors large genome, and hence, low entropy but a higher complexity. We defined that the complexity of a genome relies on the entropy of each DNA segment within genome. Then, the entropy and complexity of a genome reflects its organization solely. Exons of vertebrates harbor smaller entropies than non-exon regions (likely by the repeats that accumulated from duplications), whereas other taxonomic groups do not present this pattern. Our findings suggest that small initial population might have defined current genomic entropy and complexity: actual genomes are less complex than ancestral ones. Besides, our data disagree with the relationship between phenotype and genomic entropies previously established. Finally, by establishing the relationship between genomic entropy/complexity with the number of individuals and genome size, under an evolutive perspective, ideas concerning the genomic variability may emerge.

DNA transposon invasion and microsatellite accumulation guide W chromosome differentiation in a Neotropical fish genome.

Sex chromosome differentiation is subject to independent evolutionary processes among different lineages. The accumulation of repetitive DNAs and consequent crossing-over restriction guide the origin of the heteromorphic sex chromosome region. Several Neotropical fish species have emerged as interesting models for understanding evolution and genome diversity, although knowledge of their genomes is scarce. Here, we investigate the content of repetitive DNAs between males and females of Apareiodon sp. based on large-scale genomic data focusing on W sex chromosome differentiation. In Apareiodon, females are the heterogametic sex (ZW) and males are the homogametic sex (ZZ). The genome size estimate for Apareiodon was 1.2 Gb (with ~ 42× and ~ 47× coverage for males and females, respectively). In Apareiodon sp., approximately 36% of the genome was composed of repetitive DNAs and transposable elements (TEs) were the most abundant class. Read coverage analysis revealed different amounts of repetitive DNAs in males and females. The female-enriched clusters were located on the W sex chromosome and were mostly composed of microsatellite expansions and DNA transposons. Landscape analysis of TE contents demonstrated two major waves of invasions of TEs in the Apareiodon genome. Estimation of TE insertion times correlated with in situ locations permitted the inference that helitron, Tc1-mariner, and CMC EnSpm DNA transposons accumulated repeated copies during W chromosome differentiation between 20 and 12 million years ago. DNA transposons and microsatellite expansions appeared to be major players in W chromosome differentiation and to guide modifications in the genome content of the heteromorphic sex chromosomes.

Development and comparative analysis of yeast protein extraction protocols for mass spectrometry.

Mass spectrometry is the most used method for protein identification and quantification. Here we developed four protein extraction protocols precisely for mass spectrometry, and we compared with other ones already published. The best protocol developed by us consists on a simple extraction solution, a heat-shock step, and does not use protease inhibitor; moreover, it is the most efficient and uniform among replicates, besides to be safe, cheap and fast. That method also provided the highest number of proteins uniquely identified and allows finding a diversity of protein classes, which their absence is a problem to be avoided.

Centromeric enrichment of LINE-1 retrotransposons and its significance for the chromosome evolution of Phyllostomid bats.

Despite their ubiquitous incidence, little is known about the chromosomal distribution of long interspersed elements (LINEs) in mammalian genomes. Phyllostomid bats, characterized by lineages with distinct trends of chromosomal evolution coupled with remarkable ecological and taxonomic diversity, represent good models to understand how these repetitive sequences contribute to the evolution of genome architecture and its link to lineage diversification. To test the hypothesis that LINE-1 sequences were important modifiers of bat genome architecture, we characterized the distribution of LINE-1-derived sequences on genomes of 13 phyllostomid species within a phylogenetic framework. We found massive accumulation of LINE-1 elements in the centromeres of most species: a rare phenomenon on mammalian genomes. We hypothesize that expansion of these elements has occurred early in the radiation of phyllostomids and recurred episodically. LINE-1 expansions on centromeric heterochromatin probably spurred chromosomal change before the radiation of phyllostomids into the extant 11 subfamilies and contributed to the high degree of karyotypic variation observed among different lineages. Understanding centromere architecture in a variety of taxa promises to explain how lineage-specific changes on centromere structure can contribute to karyotypic diversity while not disrupting functional constraints for proper cell division.

Functional and evolutionary analyses of the miR156 and miR529 families in land plants.

BACKGROUND: MicroRNAs (miRNAs) are important regulatory elements of gene expression. Similarly to coding genes, miRNA genes follow a birth and death pattern of evolution likely reflecting functional relevance and divergence. For instance, miRNA529 is evolutionarily related to miRNA156 (a highly conserved miRNA in land plants), but it is lost in Arabidopsis thaliana. Interestingly, both miRNAs target sequences overlap in some members of the SQUAMOSA promoter-binding protein like (SPL) family, raising important questions regarding the diversification of the miR156/miR529-associated regulatory network in land plants. RESULTS: In this study, through phylogenic reconstruction of miR156/529 target sequences from several taxonomic groups, we have found that specific eudicot SPLs, despite miRNA529 loss, retained the corresponding target site. Detailed molecular evolutionary analyses of miR156/miR529-target sequence showed that loss of miR529 in core eudicots, such as Arabidopsis, is correlated with a more relaxed selection of the miRNA529 specific target element, while miRNA156-specific target sequence is under stronger selection, indicating that these two target sites might be under distinct evolutionary constraints. Importantly, over-expression in Arabidopsis of MIR529 precursor from a monocot, but not from a basal eudicot, demonstrates specific miR529 regulation of AtSPL9 and AtSPL15 genes, which contain conserved responsive elements for both miR156 and miR529. CONCLUSIONS: Our results suggest loss of functionality of MIR529 genes in the evolutionary history of eudicots and show that the miR529-responsive element present in some eudicot SPLs is still functional. Our data support the notion that particular miRNA156 family members might have compensated for the loss of miR529 regulation in eudicot species, which concomitantly may have favored diversification of eudicot SPLs.

Potential global cis and trans regulation of lncRNAs in Saccharomyces cerevisiae subjected to ethanol stress.

Long noncoding RNAs (lncRNAs) are regulatory RNAs. Saccharomyces cerevisiae strains transcribe hundreds of lncRNAs. LncRNAs can regulate the expression of adjacent genes (cis-regulation) or distant genes from lncRNAs (trans-regulation). Here, we analyzed the potential global cis and trans-regulation of lncRNAs of yeast subjected to ethanol stress. For potential cis regulation, for BMA641-A and S288C strains, we observed that most lncRNA-neighbor gene pairs increased the expression at a certain point followed by a decrease, and vice versa. Based on the transcriptome profile and triple helix prediction between lncRNAs and promoters of coding genes, we observed nine different ways of potential trans regulation that work in a strain-specific manner. Our data provide an initial landscape of potential cis and trans regulation in yeast, which seems to be strain-specific.

The joint action of yeast eisosomes and membraneless organelles in response to ethanol stress.

Elevated ethanol concentrations in yeast affect the plasma membrane. The plasma membrane in yeast has many lipid-protein complexes, such as Pma1 (MCP), Can1 (MCC), and the eisosome complex. We investigated the response of eisosomes, MCPs, and membraneless structures to ethanol stress. We found a correlation between ethanol stress and proton flux with quick acidification of the medium. Moreover, ethanol stress influences the symporter expression in stressed cells. We also suggest that acute stress from ethanol leads to increases in eisosome size and SG number: we hypothesized that eisosomes may protect APC symporters and accumulate an mRNA decay protein in ethanol-stressed cells. Our findings suggest that the joint action of these factors may provide a protective effect on cells under ethanol stress.

Saliva as a Biological Fluid in SARS-CoV-2 Detection.

BACKGROUND: The polymerase chain reaction of upper respiratory tract swab samples was established as the gold standard procedure for diagnosing SARS-CoV-2 during the COVID pandemic. However, saliva collection has attracted attention as an alternative diagnostic collection method. The goal of this study was to compare the use of saliva and nasopharyngeal swab (NPS) samples for the detection of SARS-CoV-2. METHODS: Ninety-nine paired samples were evaluated for the detection of SARS-CoV-2 by saliva and swab for a qualitative diagnosis and quantitative comparison of viral particles. Furthermore, the detection limits for each sample collection technique were determined. The cycle threshold (CT) values of the saliva samples, the vaccination status, and the financial costs associated with each collection technique were compared. RESULTS: The results showed qualitative equivalence in diagnosis (96.96%) comparing saliva and swab collection, although there was low quantitative agreement. Furthermore, the detection limit test demonstrated equivalence for both collection methods. We did not observe a statistically significant association between CT values and vaccination status, indicating that the vaccine had no influence on viral load at diagnosis. Finally, we observed that the use of saliva incurs lower financial costs and requires less use of plastic materials, making it more sustainable. CONCLUSIONS: These findings support the adoption of saliva collection as a feasible and sustainable alternative to the diagnosis of COVID-19.

Three topological features of regulatory networks control life-essential and specialized subsystems.

Gene regulatory networks (GRNs) play key roles in development, phenotype plasticity, and evolution. Although graph theory has been used to explore GRNs, associations amongst topological features, transcription factors (TFs), and systems essentiality are poorly understood. Here we sought the relationship amongst the main GRN topological features that influence the control of essential and specific subsystems. We found that the Knn, page rank, and degree are the most relevant GRN features: the ones are conserved along the evolution and are also relevant in pluripotent cells. Interestingly, life-essential subsystems are governed mainly by TFs with intermediary Knn and high page rank or degree, whereas specialized subsystems are mainly regulated by TFs with low Knn. Hence, we suggest that the high probability of TFs be toured by a random signal, and the high probability of the signal propagation to target genes ensures the life-essential subsystems' robustness. Gene/genome duplication is the main evolutionary process to rise Knn as the most relevant feature. Herein, we shed light on unexplored topological GRN features to assess how they are related to subsystems and how the duplications shaped the regulatory systems along the evolution. The classification model generated can be found here: https://github.com/ivanrwolf/NoC/ .

Predicting metabolic pathways of plant enzymes without using sequence similarity: Models from machine learning.

Most of the bioinformatics tools for enzyme annotation focus on enzymatic function assignments. Sequence similarity to well-characterized enzymes is often used for functional annotation and to assign metabolic pathways. However, these approaches are not feasible for all sequences leading to inaccurate annotations or lack of metabolic pathway information. Here we present the mApLe (metabolic pathway predictor of plant enzymes), a high-performance machine learning-based tool with models to label the metabolic pathway of enzymes rather than specifying enzymes' reactions. The mApLe uses molecular descriptors of the enzyme sequences to perform predictions without considering sequence similarities with reference sequences. Hence, mApLe can classify a diversity of enzymes, even the ones without any homolog or with incomplete EC numbers. This tool can be used to improve the quality of genomic annotation of plants or to narrow down the number of candidate genes for metabolic engineering researches. The mApLe tool is available online, and the GUI can be locally installed.

A New Method for Next-Generation Sequencing of the Full Hepatitis B Virus Genome from A Clinical Specimen: Impact for Virus Genotyping.

Hepatitis B virus (HBV) is an enveloped virus that induces chronic liver disease. HBV has been classified into eight genotypes (A-H) according to its genome sequence by using Sanger sequencing or reverse hybridization. Sanger sequencing is often restricted to analyzing the S gene and is inaccurate for detecting minority genetic variants, whereas reverse hybridization detects only known mutations. Next-generation sequencing (NGS) is a robust tool for clinical virology with different protocols available. The objective of this study was to develop a new method for the study of viral genetic polymorphisms or more accurate genotyping using genome amplification followed by NGS. Plasma obtained from five chronically infected HBV individuals was used for viral DNA isolation. HBV full-genome PCR amplification was the enrichment method for NGS. Primers were used to amplify all HBV genotypes in three overlapping amplicons, following a tagmentation step and Illumina NGS. For phylogenetic analysis, sequences were extracted from the HBVdb database. We were able to amplify a full HBV genome; further, NGS was shown to be a robust method and allowed better genotyping, mainly in patients carrying mixed genotypes, classified according to other techniques. This new method may be significant for whole genome analyses, including other viruses.

Transcriptional profile of a bioethanol production contaminant Candida tropicalis.

The fermentation process is widely used in the industry for bioethanol production. Even though it is widely used, microbial contamination is unpredictable and difficult to control. The problem of reduced productivity is directly linked to competition for nutrients during contamination. Yeasts representing the Candida species are frequently isolated contaminants. Elucidating the behavior of a contaminant during the fermentation cycle is essential for combatting the contamination. Consequently, the aim of the current study was to better understand the functional and transcriptional behavior of a contaminating yeast Candida tropicalis. We used a global RNA sequencing approach (RNA-seq/MiSeq) to analyze gene expression. Genes with significantly repressed or induced expression, and related to the fermentations process, such as sugar transport, pyruvate decarboxylase, amino acid metabolism, membrane, tolerance to high concentrations of ethanol and temperatures, nutrient suppression), and transcription-linked processes, were identified. The expression pattern suggested that the functional and transcriptional behavior of the contaminating yeast during fermentation for bioethanol production is similar to that of the standard yeast Saccharomyces cerevisiae. In addition, the analysis confirmed that C. tropicalis is an important contaminant of the alcoholic fermentation process, generating bioethanol and viability through its tolerance to all the adversities of a fermentation process essential for the production of bioethanol. According on the gene expression profile, many of these mechanisms are similar to those of S. cerevisiae strains currently used for bioethanol production. These mechanisms can inform studies on antimicrobials, to combat yeast contamination during industrial bioethanol production.

HIV Reverse Transcriptase and Protease Genes Variability Can Be a Biomarker Associated with HIV and Hepatitis B or C Coinfection.

Variability of the HIV reverse transcriptase (RT) and protease (PR) genes has been used as indicators of drug resistance and as a mean to evaluate phylogenetic relationships among circulating virus. However, these studies have been carried in HIV mono-infected populations. The goal of this study was to evaluate, for the first time, the HIV PR and RT sequences from HIV/HBV and HIV/HCV co-infected patients. HIV PR and RT genes were amplificated and sequenced to resistance analysis. The bioinformatics analysis was performed to infer about sequences clustering and molecular evolution. The results showed that the most frequent amino acid substitutions in RT were L214F (67.6%), I135T (55.9%), and in PR was V15I (41.2%). The molecular clock analysis showed that the HIV circulating in co-infected patients were separated in two clusters in the years 1999-2000. Some patients included as HIV mono-infected according patients' medical records and inside the co-infected cluster were, in fact, co-infected by PCR analysis. Analysis of the decision trees showed susceptibility to lamivudine and emtricitabine were important attribute to characterize co-infected patients. In conclusion, the results obtained in this study suggest, for the first time, that HIV RT and PR genes variability could be a genetic biomarker to coinfection.

Dimerization and Transactivation Domains as Candidates for Functional Modulation and Diversity of Sox9.

Sox9 plays an important role in a large variety of developmental pathways in vertebrates. It is composed of three domains: high-mobility group box (HMG box), dimerization (DIM) and transactivation (TAD). One of the main processes for regulation and variability of the pathways involving Sox9 is the self-gene expression regulation of Sox9. However, the subsequent roles of the Sox9 domains can also generate regulatory modulations. Studies have shown that TADs can bind to different types of proteins and its function seems to be influenced by DIM. Therefore, we hypothesized that both domains are directly associated and can be responsible for the functional variability of Sox9. We applied a method based on a broad phylogenetic context, using sequences of the HMG box domain, to ensure the homology of all the Sox9 copies used herein. The data obtained included 4,921 sequences relative to 657 metazoan species. Based on coevolutionary and selective pressure analyses of the Sox9 sequences, we observed coevolutions involving DIM and TADs. These data, along with the experimental data from literature, indicate a functional relationship between these domains. Moreover, DIM and TADs may be responsible for the functional plasticity of Sox9 because they are more tolerant for molecular changes (higher Ka/Ks ratio than the HMG box domain). This tolerance could allow a differential regulation of target genes or promote novel targets during transcriptional activation. In conclusion, we suggest that DIM and TADs functional association may regulate differentially the target genes or even promote novel targets during transcription activation mediated by Sox9 paralogs, contributing to the subfunctionalization of Sox9a and Sox9b in teleosts.

Sequence analyses and chromosomal distribution of the Tc1/Mariner element in Parodontidae fish (Teleostei: Characiformes).

Transposable elements are able to move along eukaryotic genomes. They are divided into two classes according to their transposition intermediate: RNA (class I or retrotransposons) or DNA (class II or DNA transposons). Most of these sequences are inactive or non-autonomous in eukaryotic genomes. Inactivate transposons can accumulate mutations at neutral rates until losing their molecular identity. They may either be eliminated from the genome or take on different molecular functions. Transposable elements may also participate in the differentiation of sex chromosomes. Therefore, the structural variations and nucleotide similarity of Tc1/Mariner sequences were analyzed along with their potential participation in the differentiation processes of sex chromosomes in the genomes of Parodontidae fish. All Parodontidae species presented non-autonomous copies of Tc1/Mariner with structural variation, different levels of deterioration (genetic distance), and variations in insertion and deletion patterns. The physical mapping of Tc1/Mariner on chromosomes revealed dispersed signals in euchromatins, with small accumulations in terminal regions and in the sex chromosomes. The gene dosage ratios indicated copy number variations of Tc1/Mariner among the genomes and high transposase open reading frame deterioration in Parodon hilarii and Parodon pongoensis genomes. This transposon presented transcriptional activity in gonads, but there was no significant difference between sexes. This may indicate non-functional protein expression or may correspond to DNA binding proteins derived from Tc1/Mariner. Thus, our results show Tc1/Mariner inactivation along with a diversity in Parodontidae genomes and its participation in the differentiation of the W sex chromosome.