Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing.

Transfer RNAs (tRNAs) play a central role in protein translation. Studying them has been difficult in part because a simple method to simultaneously quantify their abundance and chemical modifications is lacking. Here we introduce Nano-tRNAseq, a nanopore-based approach to sequence native tRNA populations that provides quantitative estimates of both tRNA abundances and modification dynamics in a single experiment. We show that default nanopore sequencing settings discard the vast majority of tRNA reads, leading to poor sequencing yields and biased representations of tRNA abundances based on their transcript length. Re-processing of raw nanopore current intensity signals leads to a 12-fold increase in the number of recovered tRNA reads and enables recapitulation of accurate tRNA abundances. We then apply Nano-tRNAseq to Saccharomyces cerevisiae tRNA populations, revealing crosstalks and interdependencies between different tRNA modification types within the same molecule and changes in tRNA populations in response to oxidative stress.

Epigenetic Response of Yarrowia lipolytica to Stress: Tracking Methylation Level and Search for Methylation Patterns via Whole-Genome Sequencing.

DNA methylation is a common, but not universal, epigenetic modification that plays an important role in multiple cellular processes. While definitely settled for numerous plant, mammalian, and bacterial species, the genome methylation in different fungal species, including widely studied and industrially-relevant yeast species, Yarrowia lipolytica, is still a matter of debate. In this paper, we report a differential DNA methylation level in the genome of Y. lipolytica subjected to sequential subculturing and to heat stress conditions. To this end, we adopted repeated batch bioreactor cultivations of Y. lipolytica subjected to thermal stress in specific time intervals. To analyze the variation in DNA methylation between stressed and control cultures, we (a) quantified the global DNA methylation status using an immuno-assay, and (b) studied DNA methylation patterns through whole-genome sequencing. Primarily, we demonstrated that 5 mC modification can be detected using a commercial immuno-assay, and that the modifications are present in Y. lipolytica's genome at ~0.5% 5 mC frequency. On the other hand, we did not observe any changes in the epigenetic response of Y. lipolytica to heat shock (HS) treatment. Interestingly, we identified a general phenomenon of decreased 5 mC level in Y. lipolytica's genome in the stationary phase of growth, when compared to a late-exponential epigenome. While this study provides an insight into the subculturing stress response and adaptation to the stress at epigenetic level by Y. lipolytica, it also leaves an open question of inability to detect any genomic DNA methylation level (either in CpG context or context-less) through whole-genome sequencing. The results of ONT sequencing, suggesting that 5 mC modification is either rare or non-existent in Y. lipolytica genome, are contradicted with the results of the immunoassay.

ModPhred: an integrative toolkit for the analysis and storage of nanopore sequencing DNA and RNA modification data.

MOTIVATION: DNA and RNA modifications can now be identified using nanopore sequencing. However, we currently lack a flexible software to efficiently encode, store, analyze and visualize DNA and RNA modification data. RESULTS: Here, we present ModPhred, a versatile toolkit that facilitates DNA and RNA modification analysis from nanopore sequencing reads in a user-friendly manner. ModPhred integrates probabilistic DNA and RNA modification information within the FASTQ and BAM file formats, can be used to encode multiple types of modifications simultaneously, and its output can be easily coupled to genomic track viewers, facilitating the visualization and analysis of DNA and RNA modification information in individual reads in a simple and computationally efficient manner. AVAILABILITY AND IMPLEMENTATION: ModPhred is available at https://github.com/novoalab/modPhred, is implemented in Python3, and is released under an MIT license. Docker images with all dependencies preinstalled are also provided. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome analysis of Candida subhashii reveals its hybrid nature and dual mitochondrial genome conformations.

Candida subhashii belongs to the CUG-Ser clade, a group of phylogenetically closely related yeast species that includes some human opportunistic pathogens, such as Candida albicans. Despite being present in the environment, C. subhashii was initially described as the causative agent of a case of peritonitis. Considering the relevance of whole-genome sequencing and analysis for our understanding of genome evolution and pathogenicity, we sequenced, assembled and annotated the genome of C. subhashii type strain. Our results show that C. subhashii presents a highly heterozygous genome and other signatures that point to a hybrid ancestry. The presence of functional pathways for assimilation of hydroxyaromatic compounds goes in line with the affiliation of this yeast with soil microbial communities involved in lignin decomposition. Furthermore, we observed that different clones of this strain may present circular or linear mitochondrial DNA. Re-sequencing and comparison of strains with differential mitochondrial genome topology revealed five candidate genes potentially associated with this conformational change: MSK1, SSZ1, ALG5, MRPL9 and OYE32.

Quantitative profiling of pseudouridylation dynamics in native RNAs with nanopore sequencing.

Nanopore RNA sequencing shows promise as a method for discriminating and identifying different RNA modifications in native RNA. Expanding on the ability of nanopore sequencing to detect N6-methyladenosine, we show that other modifications, in particular pseudouridine (Ψ) and 2'-O-methylation (Nm), also result in characteristic base-calling 'error' signatures in the nanopore data. Focusing on Ψ modification sites, we detected known and uncovered previously unreported Ψ sites in mRNAs, non-coding RNAs and rRNAs, including a Pus4-dependent Ψ modification in yeast mitochondrial rRNA. To explore the dynamics of pseudouridylation, we treated yeast cells with oxidative, cold and heat stresses and detected heat-sensitive Ψ-modified sites in small nuclear RNAs, small nucleolar RNAs and mRNAs. Finally, we developed a software, nanoRMS, that estimates per-site modification stoichiometries by identifying single-molecule reads with altered current intensity and trace profiles. This work demonstrates that Nm and Ψ RNA modifications can be detected in cellular RNAs and that their modification stoichiometry can be quantified by nanopore sequencing of native RNA.

MetaPhOrs 2.0: integrative, phylogeny-based inference of orthology and paralogy across the tree of life.

Inferring homology relationships across genes in different species is a central task in comparative genomics. Therefore, a large number of resources and methods have been developed over the years. Some public databases include phylogenetic trees of homologous gene families which can be used to further differentiate homology relationships into orthology and paralogy. MetaPhOrs is a web server that integrates phylogenetic information from different sources to provide orthology and paralogy relationships based on a common phylogeny-based predictive algorithm and associated with a consistency-based confidence score. Here we describe the latest version of the web server which includes major new implementations and provides orthology and paralogy relationships derived from ∼8.2 million gene family trees-from 13 different source repositories across ∼4000 species with sequenced genomes. MetaPhOrs server is freely available, without registration, at http://orthology.phylomedb.org/.

MasterOfPores: A Workflow for the Analysis of Oxford Nanopore Direct RNA Sequencing Datasets.

The direct RNA sequencing platform offered by Oxford Nanopore Technologies allows for direct measurement of RNA molecules without the need of conversion to complementary DNA, fragmentation or amplification. As such, it is virtually capable of detecting any given RNA modification present in the molecule that is being sequenced, as well as provide polyA tail length estimations at the level of individual RNA molecules. Although this technology has been publicly available since 2017, the complexity of the raw Nanopore data, together with the lack of systematic and reproducible pipelines, have greatly hindered the access of this technology to the general user. Here we address this problem by providing a fully benchmarked workflow for the analysis of direct RNA sequencing reads, termed MasterOfPores. The pipeline starts with a pre-processing module, which converts raw current intensities into multiple types of processed data including FASTQ and BAM, providing metrics of the quality of the run, quality-filtering, demultiplexing, base-calling and mapping. In a second step, the pipeline performs downstream analyses of the mapped reads, including prediction of RNA modifications and estimation of polyA tail lengths. Four direct RNA MinION sequencing runs can be fully processed and analyzed in 10 h on 100 CPUs. The pipeline can also be executed in GPU locally or in the cloud, decreasing the run time fourfold. The software is written using the NextFlow framework for parallelization and portability, and relies on Linux containers such as Docker and Singularity for achieving better reproducibility. The MasterOfPores workflow can be executed on any Unix-compatible OS on a computer, cluster or cloud without the need of installing any additional software or dependencies, and is freely available in Github (https://github.com/biocorecrg/master_of_pores). This workflow simplifies direct RNA sequencing data analyses, facilitating the study of the (epi)transcriptome at single molecule resolution.

Genome Sequence of the Yeast Saprochaete ingens CBS 517.90.

Chromosome-scale genome assembly of the yeast Saprochaete ingens CBS 517.90 was determined by a combination of technologies producing short (HiSeq X; Illumina) and long (MinION; Oxford Nanopore Technologies) reads. The 21.2-Mbp genome sequence has a GC content of 36.9% and codes for 6,475 predicted proteins.

Transcriptome analyses of cells carrying the Type II Csp231I restriction-modification system reveal cross-talk between two unrelated transcription factors: C protein and the Rac prophage repressor.

Restriction-modification (R-M) systems represent an effective mechanism of defence against invading bacteriophages, and are widely spread among bacteria and archaea. In acquiring a Type II R-M system via horizontal gene transfer, the new hosts become more resistant to phage infection, through the action of a restriction endonuclease (REase), which recognizes and cleaves specific target DNAs. To protect the host cell's DNA, there is also a methyltransferase (MTase), which prevents DNA cleavage by the cognate REase. In some R-M systems, the host also accepts a cis-acting transcription factor (C protein), which regulates the counteracting activities of REase and MTase to avoid host self-restriction. Our study characterized the unexpected phenotype of Escherichia coli cells, which manifested as extensive cell filamentation triggered by acquiring the Csp231I R-M system from Citrobacter sp. Surprisingly, we found that the cell morphology defect was solely dependent on the C regulator. Our transcriptome analysis supported by in vivo and in vitro assays showed that C protein directly silenced the expression of the RacR repressor to affect the Rac prophage-related genes. The rac locus ydaST genes, when derepressed, exerted a toxicity indicated by cell filamentation through an unknown mechanism. These results provide an apparent example of transcription factor cross-talk, which can have significant consequences for the host, and may represent a constraint on lateral gene transfer.

Genome Sequence of Flavor-Producing Yeast Saprochaete suaveolens NRRL Y-17571.

Saprochaete suaveolens is an ascomycetous yeast that produces a range of fruity flavors and fragrances. Here, we report the high-contiguity genome sequence of the ex-holotype strain, NRRL Y-17571 (CBS 152.25). The nuclear genome sequence contains 24.4 Mbp and codes for 8,119 predicted proteins.

Interplay of Chimeric Mating-Type Loci Impairs Fertility Rescue and Accounts for Intra-Strain Variability in Zygosaccharomyces rouxii Interspecies Hybrid ATCC42981.

The pre-whole genome duplication (WGD) Zygosaccharomyces clade comprises several allodiploid strain/species with industrially interesting traits. The salt-tolerant yeast ATCC42981 is a sterile and allodiploid strain which contains two subgenomes, one of them resembling the haploid parental species Z. rouxii. Recently, different mating-type-like (MTL) loci repertoires were reported for ATCC42981 and the Japanese strain JCM22060, which are considered two stocks of the same strain. MTL reconstruction by direct sequencing approach is challenging due to gene redundancy, structure complexities, and allodiploid nature of ATCC42981. Here, DBG2OLC and MaSuRCA hybrid de novo assemblies of ONT and Illumina reads were combined with in vitro long PCR to definitively solve these incongruences. ATCC42981 exhibits several chimeric MTL loci resulting from reciprocal translocation between parental haplotypes and retains two MATa/MATα expression loci, in contrast to MATα in JCM22060. Consistently to these reconstructions, JCM22060, but not ATCC42981, undergoes mating and meiosis. To ascertain whether the damage of one allele at the MAT locus regains the complete sexual cycle in ATCC42981, we removed the MATα expressed locus by gene deletion. The resulting MATa/- hemizygous mutants did not show any evidence of sporulation, as well as of self- and out-crossing fertility, probably because incomplete silencing at the chimeric HMLα cassette masks the loss of heterozygosity at the MAT locus. We also found that MATα deletion switched off a2 transcription, an activator of a-specific genes in pre-WGD species. These findings suggest that regulatory scheme of cell identity needs to be further investigated in Z. rouxii protoploid yeast.

Draft Genome Sequences of the Highly Halotolerant Strain Zygosaccharomyces rouxii ATCC 42981 and the Novel Allodiploid Strain Zygosaccharomyces sapae ATB301T Obtained Using the MinION Platform.

Here, we report draft genome sequences of the halotolerant and allodiploid strains Zygosaccharomyces rouxii ATCC 42981 and Zygosaccharomyces sapae ABT301T. Illumina and Oxford Nanopore MinION sequencing revealed genome sizes of 20.9 and 24.7 Mb, respectively. This information will be useful for deciphering the genetics of hybrid adaptation to high salt and sugar concentrations in nonconventional yeasts.

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.

Distinct roles of the polarity factors Boi1 and Boi2 in the control of exocytosis and abscission in budding yeast.

Boi1 and Boi2 (Boi1/2) are budding yeast plasma membrane proteins that function in polarized growth, and in cytokinesis inhibition in response to chromosome bridges via the NoCut abscission checkpoint. How Boi1/2 act in these two distinct processes is not understood. We demonstrate that Boi1/2 are required for a late step in the fusion of secretory vesicles with the plasma membrane of the growing bud. Cells lacking Boi1/2 accumulate secretory vesicles and are defective in bud growth. In contrast, Boi2 is specifically required for abscission inhibition in cells with chromatin bridges. The SH3 domain of Boi2, which is dispensable for bud growth and targets Boi2 to the site of abscission, is necessary and sufficient for abscission inhibition. Gain of function of the exocyst, a conserved protein complex involved in tethering of exocytic vesicles to the plasma membrane, rescued secretion and bud growth defects in boi mutant cells, and abrogated NoCut checkpoint function. Thus Boi2 functions redundantly with Boi1 to promote the fusion of secretory vesicles with the plasma membrane at sites of polarized growth, and acts as an abscission inhibitor during cytokinesis in response to chromatin bridges.

Eukaryotic transporters for hydroxyderivatives of benzoic acid.

Several yeast species catabolize hydroxyderivatives of benzoic acid. However, the nature of carriers responsible for transport of these compounds across the plasma membrane is currently unknown. In this study, we analyzed a family of genes coding for permeases belonging to the major facilitator superfamily (MFS) in the pathogenic yeast Candida parapsilosis. Our results revealed that these transporters are functionally equivalent to bacterial aromatic acid: H+ symporters (AAHS) such as GenK, MhbT and PcaK. We demonstrate that the genes HBT1 and HBT2 encoding putative transporters are highly upregulated in C. parapsilosis cells assimilating hydroxybenzoate substrates and the corresponding proteins reside in the plasma membrane. Phenotypic analyses of knockout mutants and hydroxybenzoate uptake assays provide compelling evidence that the permeases Hbt1 and Hbt2 transport the substrates that are metabolized via the gentisate (3-hydroxybenzoate, gentisate) and 3-oxoadipate pathway (4-hydroxybenzoate, 2,4-dihydroxybenzoate and protocatechuate), respectively. Our data support the hypothesis that the carriers belong to the AAHS family of MFS transporters. Phylogenetic analyses revealed that the orthologs of Hbt permeases are widespread in the subphylum Pezizomycotina, but have a sparse distribution among Saccharomycotina lineages. Moreover, these analyses shed additional light on the evolution of biochemical pathways involved in the catabolic degradation of hydroxyaromatic compounds.

Mitochondrial Carriers Link the Catabolism of Hydroxyaromatic Compounds to the Central Metabolism in Candida parapsilosis.

The pathogenic yeast Candida parapsilosis metabolizes hydroxyderivatives of benzene and benzoic acid to compounds channeled into central metabolism, including the mitochondrially localized tricarboxylic acid cycle, via the 3-oxoadipate and gentisate pathways. The orchestration of both catabolic pathways with mitochondrial metabolism as well as their evolutionary origin is not fully understood. Our results show that the enzymes involved in these two pathways operate in the cytoplasm with the exception of the mitochondrially targeted 3-oxoadipate CoA-transferase (Osc1p) and 3-oxoadipyl-CoA thiolase (Oct1p) catalyzing the last two reactions of the 3-oxoadipate pathway. The cellular localization of the enzymes indicates that degradation of hydroxyaromatic compounds requires a shuttling of intermediates, cofactors, and products of the corresponding biochemical reactions between cytosol and mitochondria. Indeed, we found that yeast cells assimilating hydroxybenzoates increase the expression of genes SFC1, LEU5, YHM2, and MPC1 coding for succinate/fumarate carrier, coenzyme A carrier, oxoglutarate/citrate carrier, and the subunit of pyruvate carrier, respectively. A phylogenetic analysis uncovered distinct evolutionary trajectories for sparsely distributed gene clusters coding for enzymes of both pathways. Whereas the 3-oxoadipate pathway appears to have evolved by vertical descent combined with multiple losses, the gentisate pathway shows a striking pattern suggestive of horizontal gene transfer to the evolutionarily distant Mucorales.

Correction: The Genomic Aftermath of Hybridization in the Opportunistic Pathogen Candida metapsilosis.

[This corrects the article DOI: 10.1371/journal.pgen.1005626.].

Redundans: an assembly pipeline for highly heterozygous genomes.

Many genomes display high levels of heterozygosity (i.e. presence of different alleles at the same loci in homologous chromosomes), being those of hybrid organisms an extreme such case. The assembly of highly heterozygous genomes from short sequencing reads is a challenging task because it is difficult to accurately recover the different haplotypes. When confronted with highly heterozygous genomes, the standard assembly process tends to collapse homozygous regions and reports heterozygous regions in alternative contigs. The boundaries between homozygous and heterozygous regions result in multiple assembly paths that are hard to resolve, which leads to highly fragmented assemblies with a total size larger than expected. This, in turn, causes numerous problems in downstream analyses such as fragmented gene models, wrong gene copy number, or broken synteny. To circumvent these caveats we have developed a pipeline that specifically deals with the assembly of heterozygous genomes by introducing a step to recognise and selectively remove alternative heterozygous contigs. We tested our pipeline on simulated and naturally-occurring heterozygous genomes and compared its accuracy to other existing tools. Our method is freely available at https://github.com/Gabaldonlab/redundans.

Standardized benchmarking in the quest for orthologs.

Achieving high accuracy in orthology inference is essential for many comparative, evolutionary and functional genomic analyses, yet the true evolutionary history of genes is generally unknown and orthologs are used for very different applications across phyla, requiring different precision-recall trade-offs. As a result, it is difficult to assess the performance of orthology inference methods. Here, we present a community effort to establish standards and an automated web-based service to facilitate orthology benchmarking. Using this service, we characterize 15 well-established inference methods and resources on a battery of 20 different benchmarks. Standardized benchmarking provides a way for users to identify the most effective methods for the problem at hand, sets a minimum requirement for new tools and resources, and guides the development of more accurate orthology inference methods.

The Genomic Aftermath of Hybridization in the Opportunistic Pathogen Candida metapsilosis.

Candida metapsilosis is a rarely-isolated, opportunistic pathogen that belongs to a clade of pathogenic yeasts known as the C. parapsilosis sensu lato species complex. To gain insight into the recent evolution of C. metapsilosis and the genetic basis of its virulence, we sequenced the genome of 11 clinical isolates from various locations, which we compared to each other and to the available genomes of the two remaining members of the complex: C. orthopsilosis and C. parapsilosis. Unexpectedly, we found compelling genomic evidence that C. metapsilosis is a highly heterozygous hybrid species, with all sequenced clinical strains resulting from the same past hybridization event involving two parental lineages that were approximately 4.5% divergent in sequence. This result indicates that the parental species are non-pathogenic, but that hybridization between them formed a new opportunistic pathogen, C. metapsilosis, that has achieved a worldwide distribution. We show that these hybrids are diploid and we identified strains carrying loci for both alternative mating types, which supports mating as the initial mechanism for hybrid formation. We trace the aftermath of this hybridization at the genomic level, and reconstruct the evolutionary relationships among the different strains. Recombination and introgression -resulting in loss of heterozygosis- between the two subgenomes have been rampant, and includes the partial overwriting of the MTLa mating locus in all strains. Collectively, our results shed light on the recent genomic evolution within the C. parapsilosis sensu lato complex, and argue for a re-definition of species within this clade, with at least five distinct homozygous lineages, some of which having the ability to form hybrids.