Interactions between Closely Related Bacterial Strains Are Revealed by Deep Transcriptome Sequencing

Comparative genomics, metagenomics, and single-cell technologies have shown that populations of microbial species encompass assemblages of closely related strains. This raises the question of whether individual bacterial lineages respond to the presence of their close relatives by modifying their gene expression or, instead, whether assemblages simply act as the arithmetic addition of their individual components. Here, we took advantage of transcriptome sequencing to address this question. For this, we analyzed the transcriptomes of two closely related strains of the extremely halophilic bacterium Salinibacter ruber grown axenically and in coculture. These organisms dominate bacterial assemblages in hypersaline environments worldwide. The strains used here cooccurred in the natural environment and are 100% identical in their 16S rRNA genes, and each strain harbors an accessory genome representing 10% of its complete genome. Overall, transcriptomic patterns from pure cultures were very similar for both strains. Expression was detected along practically the whole genome albeit with some genes at low levels. A subset of genes was very highly expressed in both strains, including genes coding for the light-driven proton pump xanthorhodopsin, genes involved in the stress response, and genes coding for transcriptional regulators. Expression differences between pure cultures affected mainly genes involved in environmental sensing. When the strains were grown in coculture, there was a modest but significant change in their individual transcription patterns compared to those in pure culture. Each strain sensed the presence of the other and responded in a specific manner, which points to fine intraspecific transcriptomic modulation.The group of J.A. is funded by grant CGL2012-39627-C03-01 from the Spanish Ministry of Economy and Competitiveness (MINECO), which is cofinanced with FEDER support from the European Union. P.G.-T. was an FPI-MINECO fellow. Research by the group of T.G. is funded in part by a grant from the Spanish Ministry of Economy and Competitiveness (BIO2012-37161), a grant from the Qatar National Research Fund (NPRP 5-298-3-086), and a grant from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC (grant agreement no. ERC-2012-StG-310325). L.P.P. was funded through the La Caixa-CRG international fellowship program

Genome analysis of five recently described species of the CUG-Ser clade uncovers Candida theae as a new hybrid lineage with pathogenic potential in the Candida parapsilosis species complex

Candida parapsilosis species complex comprises three important pathogenic species: Candida parapsilosis sensu stricto, Candida orthopsilosis and Candida metapsilosis. The majority of C. orthopsilosis and all C. metapsilosis isolates sequenced thus far are hybrids, and most of the parental lineages remain unidentified. This led to the hypothesis that hybrids with pathogenic potential were formed by the hybridization of non-pathogenic lineages that thrive in the environment. In a search for the missing hybrid parentals, and aiming to get a better understanding of the evolution of the species complex, we sequenced, assembled and analysed the genome of five close relatives isolated from the environment: Candida jiufengensis, Candida pseudojiufengensis, Candida oxycetoniae, Candida margitis and Candida theae. We found that the linear conformation of mitochondrial genomes in Candida species emerged multiple times independently. Furthermore, our analyses discarded the possible involvement of these species in the mentioned hybridizations, but identified C. theae as an additional hybrid in the species complex. Importantly, C. theae was recently associated with a case of infection, and we also uncovered the hybrid nature of this clinical isolate. Altogether, our results reinforce the hypothesis that hybridization is widespread among Candida species, and potentially contributes to the emergence of lineages with opportunistic pathogenic behaviour.This work received funding from the European Union’s Horizon 2020 Research and Innovation Program under the Marie Skłodowska-Curie Grant Agreement No. H2020-MSCA-ITN-2014-642095. T.G. group also acknowledges support from the Spanish Ministry of Science and Innovation for grant PGC2018-099921-B-I00, cofounded by European Regional Development Fund (ERDF); from the Catalan Research Agency (AGAUR) SGR423; from the European Union’s Horizon 2020 research and innovation programme (ERC-2016-724173); from the Gordon and Betty Moore Foundation (Grant GBMF9742) and from the Instituto de Salud Carlos III (INB Grant PT17/0009/0023—ISCIII-SGEFI/ERDF). J.N. group was supported by the Slovak Research and Development Agency (APVV-18-0239, APVV-20-0166) and the Scientific Grant Agency of the Ministry of Education, Science, Research and Sport of the Slovak Republic (VEGA 1/0027/19).Peer ReviewedPostprint (published version

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\u3b1 expression loci, in contrast to MAT\u3b1 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\u3b1 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\u3b1 cassette masks the loss of heterozygosity at the MAT locus. We also found that MAT\u3b1 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\u2009Mb, respectively. This information will be useful for deciphering the genetics of hybrid adaptation to high salt and sugar concentrations in nonconventional yeasts

Transcriptome Analyses of Mosaic (MSC) Mitochondrial Mutants of Cucumber in a Highly Inbred Nuclear Background.

Cucumber (Cucumis sativus L.) has a large, paternally transmitted mitochondrial genome. Cucumber plants regenerated from cell cultures occasionally show paternally transmitted mosaic (MSC) phenotypes, characterized by slower growth, chlorotic patterns on the leaves and fruit, lower fertility, and rearrangements in their mitochondrial DNAs (mtDNAs). MSC lines 3, 12, and 16 originated from different cell cultures all established using the highly inbred, wild-type line B. These MSC lines possess different rearrangements and under-represented regions in their mtDNAs. We completed RNA-seq on normalized and non-normalized cDNA libraries from MSC3, MSC12, and MSC16 to study their nuclear gene-expression profiles relative to inbred B. Results from both libraries indicated that gene expression in MSC12 and MSC16 were more similar to each other than MSC3. Forty-one differentially expressed genes (DEGs) were upregulated and one downregulated in the MSC lines relative to B. Gene functional classifications revealed that more than half of these DEGs are associated with stress-response pathways. Consistent with this observation, we detected elevated levels of hydrogen peroxide throughout leaf tissue in all MSC lines compared to wild-type line B. These results demonstrate that independently produced MSC lines with different mitochondrial polymorphisms show unique and shared nuclear responses. This study revealed genes associated with stress response that could become selection targets to develop cucumber cultivars with increased stress tolerance, and further support of cucumber as a model plant to study nuclear-mitochondrial interactions

MetaPhOrs: orthology and paralogy predictions from multiple phylogenetic evidence using a consistency-based confidence score

Reliable prediction of orthology is central to comparative genomics. Approaches based on phylogenetic analyses closely resemble the original definition of orthology and paralogy and are known to be highly accurate. However, the large computational cost associated to these analyses is a limiting factor that often prevents its use at genomic scales. Recently, several projects have addressed the reconstruction of large collections of high-quality phylogenetic trees from which orthology and paralogy relationships can be inferred. This provides us with the opportunity to infer the evolutionary relationships of genes from multiple, independent, phylogenetic trees. Using such strategy, we combine phylogenetic information derived from different databases, to predict orthology and paralogy relationships for 4.1 million proteins in 829 fully sequenced genomes. We show that the number of independent sources from which a prediction is made, as well as the level of consistency across predictions, can be used as reliable confidence scores. A webserver has been developed to easily access these data (http://orthology.phylomedb.org), which provides users with a global repository of phylogeny-based orthology and paralogy predictions

Ancient dispersal of the human fungal pathogen Cryptococcus gattii from the Amazon rainforest.

Over the past two decades, several fungal outbreaks have occurred, including the high-profile 'Vancouver Island' and 'Pacific Northwest' outbreaks, caused by Cryptococcus gattii, which has affected hundreds of otherwise healthy humans and animals. Over the same time period, C. gattii was the cause of several additional case clusters at localities outside of the tropical and subtropical climate zones where the species normally occurs. In every case, the causative agent belongs to a previously rare genotype of C. gattii called AFLP6/VGII, but the origin of the outbreak clades remains enigmatic. Here we used phylogenetic and recombination analyses, based on AFLP and multiple MLST datasets, and coalescence gene genealogy to demonstrate that these outbreaks have arisen from a highly-recombining C. gattii population in the native rainforest of Northern Brazil. Thus the modern virulent C. gattii AFLP6/VGII outbreak lineages derived from mating events in South America and then dispersed to temperate regions where they cause serious infections in humans and animals

Standardised Benchmarking in the Quest for Orthologs

The identification of evolutionarily related genes across different species—orthologs in particular—forms the backbone of many comparative, evolutionary, and functional genomic analyses. Achieving high accuracy in orthology inference is thus essential. Yet the true evolutionary history of genes, required to ascertain orthology, is generally unknown. Furthermore, orthologs are used for very different applications across different phyla, with different requirements in terms of the precision-recall trade-off. As a result, assessing the performance of orthology inference methods remains difficult for both users and method developers. Here, we present a community effort to establish standards in orthology benchmarking and facilitate orthology benchmarking through an automated web-based service (http://orthology.benchmarkservice.org). Using this new service, we characterise the performance of 15 well-established orthology inference methods and resources on a battery of 20 different benchmarks. Standardised benchmarking provides a way for users to identify the most effective methods for the problem at hand, sets a minimal requirement for new tools and resources, and guides the development of more accurate orthology inference methods

Comparative genomics to unravel virulence mechanisms in fungal human pathogens

Las Candidas forman uno de los grupos con mayor n´umero de hongos pat´ogenos en humanos. La relaci ´on que tienen desde un punto de vista filogen´etico demuestra que la capacidad de infectar humanos ha surgido varias veces de forma independiente en este clado. El complejo de Candida parapsilosis es ideal para investigar la aparici´on de virulencia puesto que contiene tres especies cercanas que muestran diferentes grados de virulencia y de importancia: C. parapsilosis, C. orthopsilosis y C. metapsilosis. En esta tesis presento la secuenciaci´on y posterior an´alisis de quince cepas de origen cl´ınico y medioambiental que forman parte de este clado. Cabe remarcar que los genomas de C. orthopsilosis tipo 1 y de C. metapsilosis no se hab´ıan secuenciado previamente. Nuestros resultados muestran evidencia gen´omica de la existencia de recombinaci´on, apareamiento e hibridaci´on en este clado, que previamente se habia considerado asexual. Proponemos que las cepas cl´ınicas emergieron de forma independiente en linajes de cepas encontradas en el medioambiente y una posible relaci ´on entre la hibridaci´on y la adquisici´on de caracter´ısticas relevantes para la virulencia. Finalmente, para estudiar la aparici´on de virulencia en este clado, hemos comparado, por primera vez, los genomas de las tres especies que se encuentran dentro del complejo de Candida parapsilosis. Hemos encontrado expansiones de familias g´enicas previamente implicadas en virulencia, como es el caso de las adhesinas, transportadores de membrana y enzimas extracelulares. Tambi´en hemos observado expansiones de familias de genes que hasta el momento no han estado relacionadas con virulencia. En resumen, nuestros resultados aportan una base para poder elaborar numerosas hip´ otesis sobre la aparici´on de virulencia en las especies de Candida y para que estas hip´ otesis puedan ser comprobadas posteriormente en el laboratorio.Candida species constitute one of the most prevalent groups of fungal pathogens. From their phylogenetic relationships it is clear that virulence to humans has emerged in this clade several, independent times. The Candida parapsilosis complex is particularly suitable to investigate the emergence of virulence, with three closely-related species of varying degree of pathogenicity and of growing relevance: C. parapsilosis, C. orthopsilosis and C. metapsilosis. In this thesis I present the genome sequencing and analysis of fifteen strains from this clade, sampled from clinical and environmental sources. Notably, genomes of C. orthopsilosis Type 1 and C. metapsilosis were sequenced for the first time. Our results show for the first time the genomic evidence for the existence of recombination, mating and hybridization in this clade, previously considered asexual. We propose the independent emergence of clinical isolates from environmental lineages and a possible role of hybridization in the acquisition of relevant traits for pathogenesis. Finally, in order to gain insight into the emergence of virulence in this clade, we have compared the genomes of all three species of C. parapsilosis complex. We have found expansions in gene families known to be involved in virulence, like adhesins, membrane transporters and extracellular enzymes, as well as expansions in gene families not implicated in virulence so far. Altogether, our findings provide the grounds for numerous hypotheses about the emergence of virulence in Candida spp. and for their future experimental testing.

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.Spanish Ministry of Economy and Competitiveness grants, ‘Centro de Excelencia Severo Ochoa [2013–2017’ SEV-2012-0208, BFU2015-67107 to TG group] cofounded by European Regional Development Fund (ERDF); European Union and ERC Seventh Framework Programme [FP7/2007-2013] under grant agreements [FP7-PEOPLE-2013-ITN-606786 and ERC-2012-StG-310325]; European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie [H2020-MSCA-ITN-2014-642095]; La Caixa-CRG International Fellowship Program (to L.P.P.). Funding for open access charge: ERC Seventh Framework Programme [FP7/2007-2013] under grant agreements [FP7-PEOPLE-2013-ITN-606786 and ERC-2012-StG-310325]; European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie [H2020-MSCA-ITN-2014-642095]; Spanish Ministry of Economy and Competitiveness grants, ‘Centro de Excelencia Severo Ochoa [2013–2017’ SEV-2012-0208, BFU2015-67107] cofounded by ERDF; Catalan Research Agency (AGAUR) [SGR857]