Telomere-to-telomere Schizosaccharomyces japonicus genome assembly reveals hitherto unknown genome features.

Schizosaccharomyces japonicus belongs to the single-genus class Schizosaccharomycetes, otherwise known as "fission yeasts." As part of a composite model system with its widely studied S. pombe sister species, S. japonicus has provided critical insights into the workings and the evolution of cell biological mechanisms. Furthermore, its divergent biology makes S. japonicus a valuable model organism in its own right. However, the currently available genome assembly contains gaps and has been unable to resolve centromeres and other repeat-rich chromosomal regions. Here we present a telomere-to-telomere long-read genome assembly of the S. japonicus genome. This includes the three megabase-length chromosomes, with centromeres hundreds of kilobases long, rich in 5S ribosomal RNA genes, transfer RNA genes, long terminal repeats, and short repeats. We identify a gene-sparse region on chromosome 2 that resembles a 331 kb centromeric duplication. We revise the genome size of S. japonicus to at least 16.6 Mb and possibly up to 18.12 Mb, at least 30% larger than previous estimates. Our whole genome assembly will support the growing S. japonicus research community and facilitate research in new directions, including centromere and DNA repeat evolution, and yeast comparative genomics.

Schizosaccharomyces versatilis represents a distinct evolutionary lineage of fission yeast.

The fission yeast species Schizosaccharomyces japonicus is currently divided into two varieties-S. japonicus var. japonicus and S. japonicus var. versatilis. Here we examine the var. versatilis isolate CBS5679. The CBS5679 genome shows 88% identity to the reference genome of S. japonicus var. japonicus at the coding sequence level, with phylogenetic analyses suggesting that it has split from the S. japonicus lineage 25 million years ago. The CBS5679 genome contains a reciprocal translocation between chromosomes 1 and 2, together with several large inversions. The products of genes linked to the major translocation are associated with 'metabolism' and 'cellular assembly' ontology terms. We further show that CBS5679 does not generate viable progeny with the reference strain of S. japonicus. Although CBS5679 shares closer similarity to the 'type' strain of var. versatilis as compared to S. japonicus, it is not identical to the type strain, suggesting population structure within var. versatilis. We recommend that the taxonomic status of S. japonicus var. versatilis is raised, with it being treated as a separate species, Schizosaccharomyces versatilis.

Panhematopoietic RNA barcoding enables kinetic measurements of nucleate and anucleate lineages and the activation of myeloid clones following acute platelet depletion.

BACKGROUND: Platelets and erythrocytes constitute over 95% of all hematopoietic stem cell output. However, the clonal dynamics of HSC contribution to these lineages remains largely unexplored. RESULTS: We use lentiviral genetic labeling of mouse hematopoietic stem cells to quantify output from all lineages, nucleate, and anucleate, simultaneously linking these with stem and progenitor cell transcriptomic phenotypes using single-cell RNA-sequencing. We observe dynamic shifts of clonal behaviors through time in same-animal peripheral blood and demonstrate that acute platelet depletion shifts the output of multipotent hematopoietic stem cells to the exclusive production of platelets. Additionally, we observe the emergence of new myeloid-biased clones, which support short- and long-term production of blood cells. CONCLUSIONS: Our approach enables kinetic studies of multi-lineage output in the peripheral blood and transcriptional heterogeneity of individual hematopoietic stem cells. Our results give a unique insight into hematopoietic stem cell reactivation upon platelet depletion and of clonal dynamics in both steady state and under stress.

Extensive genome introgression between domestic ferret and European polecat during population recovery in Great Britain.

The European polecat (Mustela putorius) is a mammalian predator which occurs across much of Europe east to the Ural Mountains. In Great Britain, following years of persecution the range of the European polecat contracted and by the early 1900s was restricted to unmanaged forests of central Wales. The European polecat has recently undergone a population increase due to legal protection and its range now overlaps that of feral domestic ferrets (Mustela putorius furo). During this range expansion, European polecats hybridized with feral domestic ferrets producing viable offspring. Here, we carry out population-level whole-genome sequencing on 8 domestic ferrets, 19 British European polecats, and 15 European polecats from the European mainland. We used a range of population genomics methods to examine the data, including phylogenetics, phylogenetic graphs, model-based clustering, phylogenetic invariants, ABBA-BABA tests, topology weighting, and Fst. We found high degrees of genome introgression in British polecats outside their previous stronghold, even in those individuals phenotyped as "pure" polecats. These polecats ranged from presumed F1 hybrids (gamma = 0.53) to individuals that were much less introgressed (gamma = 0.2). We quantify this introgression and find introgressed genes containing Fst outliers associated with cognitive function and sight.

Whole genome resequencing data enables a targeted SNP panel for conservation and aquaculture of Oreochromis cichlid fishes.

Cichlid fish of the genus Oreochromis form the basis of the global tilapia aquaculture and fisheries industries. Broodstocks for aquaculture are often collected from wild populations, which in Africa may be from locations containing multiple Oreochromis species. However, many species are difficult to distinguish morphologically, hampering efforts to maintain good quality farmed strains. Additionally, non-native farmed tilapia populations are known to be widely distributed across Africa and to hybridize with native Oreochromis species, which themselves are important for capture fisheries. The morphological identification of these hybrids is particularly unreliable. Here, we describe the development of a single nucleotide polymorphism (SNP) genotyping panel from whole-genome resequencing data that enables targeted species identification in Tanzania. We demonstrate that an optimized panel of 96 genome-wide SNPs based on FST outliers performs comparably to whole genome resequencing in distinguishing species and identifying hybrids. We also show this panel outperforms microsatellite-based and phenotype-based classification methods. Case studies indicate several locations where introduced aquaculture species have become established in the wild, threatening native Oreochromis species. The novel SNP markers identified here represent an important resource for assessing broodstock purity in hatcheries and helping to conserve unique endemic biodiversity.

Describing variability in pig genes involved in coronavirus infections for a One Health perspective in conservation of animal genetic resources.

Coronaviruses silently circulate in human and animal populations, causing mild to severe diseases. Therefore, livestock are important components of a "One Health" perspective aimed to control these viral infections. However, at present there is no example that considers pig genetic resources in this context. In this study, we investigated the variability of four genes (ACE2, ANPEP and DPP4 encoding for host receptors of the viral spike proteins and TMPRSS2 encoding for a host proteinase) in 23 European (19 autochthonous and three commercial breeds and one wild boar population) and two Asian Sus scrofa populations. A total of 2229 variants were identified in the four candidate genes: 26% of them were not previously described; 29 variants affected the protein sequence and might potentially interact with the infection mechanisms. The results coming from this work are a first step towards a "One Health" perspective that should consider conservation programs of pig genetic resources with twofold objectives: (i) genetic resources could be reservoirs of host gene variability useful to design selection programs to increase resistance to coronaviruses; (ii) the described variability in genes involved in coronavirus infections across many different pig populations might be part of a risk assessment including pig genetic resources.

Single Nucleotide Polymorphisms in the Bovine TLR2 Extracellular Domain Contribute to Breed and Species-Specific Innate Immune Functionality.

Recent evidence suggests that several cattle breeds may be more resistant to infection with the zoonotic pathogen Mycobacterium bovis. Our data presented here suggests that the response to mycobacterial antigens varies in macrophages generated from Brown Swiss (BS) and Holstein Friesian (HF) cattle, two breeds belonging to the Bos taurus family. Whole genome sequencing of the Brown Swiss genome identified several potential candidate genes, in particular Toll-like Receptor-2 (TLR2), a pattern recognition receptor (PRR) that has previously been described to be involved in mycobacterial recognition. Further investigation revealed single nucleotide polymorphisms (SNP) in TLR2 that were identified between DNA isolated from cells of BS and HF cows. Interestingly, one specific SNP, H326Q, showed a different genotype frequency in two cattle subspecies, Bos (B.) taurus and Bos indicus. Cloning of the TLR2 gene and subsequent gene-reporter and chemokine assays revealed that this SNP, present in BS and Bos indicus breeds, resulted in a significantly higher response to mycobacterial antigens as well as tri-acylated lipopeptide ligands in general. Comparing wild-type and H326Q containing TLR2 responses, wild-type bovine TLR2 response showed clear, diminished mycobacterial antigen responses compared to human TLR2, however bovine TLR2 responses containing H326Q were found to be partially recovered compared to human TLR2. The creation of human:bovine TLR2 chimeras increased the response to mycobacterial antigens compared to the full-length bovine TLR2, but significantly reduced the response compared to the full-length human TLR2. Thus, our data, not only present evidence that TLR2 is a major PRR in the mammalian species-specific response to mycobacterial antigens, but furthermore, that there are clear differences between the response seen in different cattle breeds, which may contribute to their enhanced or reduced susceptibility to mycobacterial infection.

Newly discovered cichlid fish biodiversity threatened by hybridization with non-native species.

Invasive freshwater fishes are known to readily hybridize with indigenous congeneric species, driving loss of unique and irreplaceable genetic resources. Here we reveal that newly discovered (2013-2016) evolutionarily significant populations of Korogwe tilapia (Oreochromis korogwe) from southern Tanzania are threatened by hybridization with the larger invasive Nile tilapia (Oreochromis niloticus). We use a combination of morphology, microsatellite allele frequencies and whole genome sequences to show that O. korogwe from southern lakes (Nambawala, Rutamba and Mitupa) are distinct from geographically disjunct populations in northern Tanzania (Zigi River and Mlingano Dam). We also provide genetic evidence of O. korogwe × niloticus hybrids in three southern lakes and demonstrate heterogeneity in the extent of admixture across the genome. Finally, using the least admixed genomic regions we estimate that the northern and southern O. korogwe populations most plausibly diverged ~140,000 years ago, suggesting that the geographical separation of the northern and southern groups is not a result of a recent translocation, and instead these populations represent independent evolutionarily significant units. We conclude that these newly discovered and phenotypically unique cichlid populations are already threatened by hybridization with an invasive species, and propose that these irreplaceable genetic resources would benefit from conservation interventions.

A single-cell RNA-sequencing training and analysis suite using the Galaxy framework.

BACKGROUND: The vast ecosystem of single-cell RNA-sequencing tools has until recently been plagued by an excess of diverging analysis strategies, inconsistent file formats, and compatibility issues between different software suites. The uptake of 10x Genomics datasets has begun to calm this diversity, and the bioinformatics community leans once more towards the large computing requirements and the statistically driven methods needed to process and understand these ever-growing datasets. RESULTS: Here we outline several Galaxy workflows and learning resources for single-cell RNA-sequencing, with the aim of providing a comprehensive analysis environment paired with a thorough user learning experience that bridges the knowledge gap between the computational methods and the underlying cell biology. The Galaxy reproducible bioinformatics framework provides tools, workflows, and trainings that not only enable users to perform 1-click 10x preprocessing but also empower them to demultiplex raw sequencing from custom tagged and full-length sequencing protocols. The downstream analysis supports a range of high-quality interoperable suites separated into common stages of analysis: inspection, filtering, normalization, confounder removal, and clustering. The teaching resources cover concepts from computer science to cell biology. Access to all resources is provided at the singlecell.usegalaxy.eu portal. CONCLUSIONS: The reproducible and training-oriented Galaxy framework provides a sustainable high-performance computing environment for users to run flexible analyses on both 10x and alternative platforms. The tutorials from the Galaxy Training Network along with the frequent training workshops hosted by the Galaxy community provide a means for users to learn, publish, and teach single-cell RNA-sequencing analysis.

Whole-genome sequencing of European autochthonous and commercial pig breeds allows the detection of signatures of selection for adaptation of genetic resources to different breeding and production systems.

BACKGROUND: Natural and artificial directional selection in cosmopolitan and autochthonous pig breeds and wild boars have shaped their genomes and resulted in a reservoir of animal genetic diversity. Signatures of selection are the result of these selection events that have contributed to the adaptation of breeds to different environments and production systems. In this study, we analysed the genome variability of 19 European autochthonous pig breeds (Alentejana, Bísara, Majorcan Black, Basque, Gascon, Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano, Sarda, Krskopolje pig, Black Slavonian, Turopolje, Moravka, Swallow-Bellied Mangalitsa, Schwäbisch-Hällisches Schwein, Lithuanian indigenous wattle and Lithuanian White old type) from nine countries, three European commercial breeds (Italian Large White, Italian Landrace and Italian Duroc), and European wild boars, by mining whole-genome sequencing data obtained by using a DNA-pool sequencing approach. Signatures of selection were identified by using a single-breed approach with two statistics [within-breed pooled heterozygosity (HP) and fixation index (FST)] and group-based FST approaches, which compare groups of breeds defined according to external traits and use/specialization/type. RESULTS: We detected more than 22 million single nucleotide polymorphisms (SNPs) across the 23 compared populations and identified 359 chromosome regions showing signatures of selection. These regions harbour genes that are already known or new genes that are under selection and relevant for the domestication process in this species, and that affect several morphological and physiological traits (e.g. coat colours and patterns, body size, number of vertebrae and teats, ear size and conformation, reproductive traits, growth and fat deposition traits). Wild boar related signatures of selection were detected across all the genome of several autochthonous breeds, which suggests that crossbreeding (accidental or deliberate) occurred with wild boars. CONCLUSIONS: Our findings provide a catalogue of genetic variants of many European pig populations and identify genome regions that can explain, at least in part, the phenotypic diversity of these genetic resources.

Sequencing smart: De novo sequencing and assembly approaches for a non-model mammal.

BACKGROUND: Whilst much sequencing effort has focused on key mammalian model organisms such as mouse and human, little is known about the relationship between genome sequencing techniques for non-model mammals and genome assembly quality. This is especially relevant to non-model mammals, where the samples to be sequenced are often degraded and of low quality. A key aspect when planning a genome project is the choice of sequencing data to generate. This decision is driven by several factors, including the biological questions being asked, the quality of DNA available, and the availability of funds. Cutting-edge sequencing technologies now make it possible to achieve highly contiguous, chromosome-level genome assemblies, but rely on high-quality high molecular weight DNA. However, funding is often insufficient for many independent research groups to use these techniques. Here we use a range of different genomic technologies generated from a roadkill European polecat (Mustela putorius) to assess various assembly techniques on this low-quality sample. We evaluated different approaches for de novo assemblies and discuss their value in relation to biological analyses. RESULTS: Generally, assemblies containing more data types achieved better scores in our ranking system. However, when accounting for misassemblies, this was not always the case for Bionano and low-coverage 10x Genomics (for scaffolding only). We also find that the extra cost associated with combining multiple data types is not necessarily associated with better genome assemblies. CONCLUSIONS: The high degree of variability between each de novo assembly method (assessed from the 7 key metrics) highlights the importance of carefully devising the sequencing strategy to be able to carry out the desired analysis. Adding more data to genome assemblies does not always result in better assemblies, so it is important to understand the nuances of genomic data integration explained here, in order to obtain cost-effective value for money when sequencing genomes.

Genes Encoding Recognition of the Cladosporium fulvum Effector Protein Ecp5 Are Encoded at Several Loci in the Tomato Genome.

The molecular interactions between tomato and Cladosporium fulvum have been an important model for molecular plant pathology. Complex genetic loci on tomato chromosomes 1 and 6 harbor genes for resistance to Cladosporium fulvum, encoding receptor like-proteins that perceive distinct Cladosporium fulvum effectors and trigger plant defenses. Here, we report classical mapping strategies for loci in tomato accessions that respond to Cladosporium fulvum effector Ecp5, which is very sequence-monomorphic. We screened 139 wild tomato accessions for an Ecp5-induced hypersensitive response, and in five accessions, the Ecp5-induced hypersensitive response segregated as a monogenic trait, mapping to distinct loci in the tomato genome. We identified at least three loci on chromosomes 1, 7 and 12 that harbor distinct Cf-Ecp5 genes in four different accessions. Our mapping showed that the Cf-Ecp5 in Solanum pimpinellifolium G1.1161 is located at the Milky Way locus. The Cf-Ecp5 in Solanum pimpinellifolium LA0722 was mapped to the bottom arm of chromosome 7, while the Cf-Ecp5 genes in Solanum lycopersicum Ontario 7522 and Solanum pimpinellifolium LA2852 were mapped to the same locus on the top arm of chromosome 12. Bi-parental crosses between accessions carrying distinct Cf-Ecp5 genes revealed putative genetically unlinked suppressors of the Ecp5-induced hypersensitive response. Our mapping also showed that Cf-11 is located on chromosome 11, close to the Cf-3 locus. The Ecp5-induced hypersensitive response is widely distributed within tomato species and is variable in strength. This novel example of convergent evolution could be used for choosing different functional Cf-Ecp5 genes according to individual plant breeding needs.

A Galaxy-based training resource for single-cell RNA-sequencing quality control and analyses.

BACKGROUND: It is not a trivial step to move from single-cell RNA-sequencing (scRNA-seq) data production to data analysis. There is a lack of intuitive training materials and easy-to-use analysis tools, and researchers can find it difficult to master the basics of scRNA-seq quality control and the later analysis. RESULTS: We have developed a range of practical scripts, together with their corresponding Galaxy wrappers, that make scRNA-seq training and quality control accessible to researchers previously daunted by the prospect of scRNA-seq analysis. We implement a "visualize-filter-visualize" paradigm through simple command line tools that use the Loom format to exchange data between the tools. The point-and-click nature of Galaxy makes it easy to assess, visualize, and filter scRNA-seq data from short-read sequencing data. CONCLUSION: We have developed a suite of scRNA-seq tools that can be used for both training and more in-depth analyses.

Adaptation and conservation insights from the koala genome.

The koala, the only extant species of the marsupial family Phascolarctidae, is classified as 'vulnerable' due to habitat loss and widespread disease. We sequenced the koala genome, producing a complete and contiguous marsupial reference genome, including centromeres. We reveal that the koala's ability to detoxify eucalypt foliage may be due to expansions within a cytochrome P450 gene family, and its ability to smell, taste and moderate ingestion of plant secondary metabolites may be due to expansions in the vomeronasal and taste receptors. We characterized novel lactation proteins that protect young in the pouch and annotated immune genes important for response to chlamydial disease. Historical demography showed a substantial population crash coincident with the decline of Australian megafauna, while contemporary populations had biogeographic boundaries and increased inbreeding in populations affected by historic translocations. We identified genetically diverse populations that require habitat corridors and instituting of translocation programs to aid the koala's survival in the wild.

Arabidopsis late blight: infection of a nonhost plant by Albugo laibachii enables full colonization by Phytophthora infestans.

The oomycete pathogen Phytophthora infestans causes potato late blight, and as a potato and tomato specialist pathogen, is seemingly poorly adapted to infect plants outside the Solanaceae. Here, we report the unexpected finding that P. infestans can infect Arabidopsis thaliana when another oomycete pathogen, Albugo laibachii, has colonized the host plant. The behaviour and speed of P. infestans infection in Arabidopsis pre-infected with A. laibachii resemble P. infestans infection of susceptible potato plants. Transcriptional profiling of P. infestans genes during infection revealed a significant overlap in the sets of secreted-protein genes that are induced in P. infestans upon colonization of potato and susceptible Arabidopsis, suggesting major similarities in P. infestans gene expression dynamics on the two plant species. Furthermore, we found haustoria of A. laibachii and P. infestans within the same Arabidopsis cells. This Arabidopsis-A. laibachii-P. infestans tripartite interaction opens up various possibilities to dissect the molecular mechanisms of P. infestans infection and the processes occurring in co-infected Arabidopsis cells.

Loss of a Conserved tRNA Anticodon Modification Perturbs Plant Immunity.

tRNA is the most highly modified class of RNA species, and modifications are found in tRNAs from all organisms that have been examined. Despite their vastly different chemical structures and their presence in different tRNAs, occurring in different locations in tRNA, the biosynthetic pathways of the majority of tRNA modifications include a methylation step(s). Recent discoveries have revealed unprecedented complexity in the modification patterns of tRNA, their regulation and function, suggesting that each modified nucleoside in tRNA may have its own specific function. However, in plants, our knowledge on the role of individual tRNA modifications and how they are regulated is very limited. In a genetic screen designed to identify factors regulating disease resistance and activation of defenses in Arabidopsis, we identified SUPPRESSOR OF CSB3 9 (SCS9). Our results reveal SCS9 encodes a tRNA methyltransferase that mediates the 2´-O-ribose methylation of selected tRNA species in the anticodon loop. These SCS9-mediated tRNA modifications enhance during the course of infection with the bacterial pathogen Pseudomonas syringae DC3000, and lack of such tRNA modification, as observed in scs9 mutants, severely compromise plant immunity against the same pathogen without affecting the salicylic acid (SA) signaling pathway which regulates plant immune responses. Our results support a model that gives importance to the control of certain tRNA modifications for mounting an effective immune response in Arabidopsis, and therefore expands the repertoire of molecular components essential for an efficient disease resistance response.

bio-samtools 2: a package for analysis and visualization of sequence and alignment data with SAMtools in Ruby.

MOTIVATION: bio-samtools is a Ruby language interface to SAMtools, the highly popular library that provides utilities for manipulating high-throughput sequence alignments in the Sequence Alignment/Map format. Advances in Ruby, now allow us to improve the analysis capabilities and increase bio-samtools utility, allowing users to accomplish a large amount of analysis using a very small amount of code. bio-samtools can also be easily developed to include additional SAMtools methods and hence stay current with the latest SAMtools releases. RESULTS: We have added new Ruby classes for the MPileup and Variant Call Format (VCF) data formats emitted by SAMtools and introduced more analysis methods for variant analysis, including alternative allele calculation and allele frequency calling for SNPs. Our new implementation of bio-samtools also ensures that all the functionality of the SAMtools library is now supported and that bio-samtools can be easily extended to include future changes in SAMtools. bio-samtools 2 also provides methods that allow the user to directly produce visualization of alignment data.

EXPRSS: an Illumina based high-throughput expression-profiling method to reveal transcriptional dynamics.

BACKGROUND: Next Generation Sequencing technologies have facilitated differential gene expression analysis through RNA-seq and Tag-seq methods. RNA-seq has biases associated with transcript lengths, lacks uniform coverage of regions in mRNA and requires 10-20 times more reads than a typical Tag-seq. Most existing Tag-seq methods either have biases or not high throughput due to use of restriction enzymes or enzymatic manipulation of 5' ends of mRNA or use of RNA ligations. RESULTS: We have developed EXpression Profiling through Randomly Sheared cDNA tag Sequencing (EXPRSS) that employs acoustic waves to randomly shear cDNA and generate sequence tags at a relatively defined position (~150-200 bp) from the 3' end of each mRNA. Implementation of the method was verified through comparative analysis of expression data generated from EXPRSS, NlaIII-DGE and Affymetrix microarray and through qPCR quantification of selected genes. EXPRSS is a strand specific and restriction enzyme independent tag sequencing method that does not require cDNA length-based data transformations. EXPRSS is highly reproducible, is high-throughput and it also reveals alternative polyadenylation and polyadenylated antisense transcripts. It is cost-effective using barcoded multiplexing, avoids the biases of existing SAGE and derivative methods and can reveal polyadenylation position from paired-end sequencing. CONCLUSIONS: EXPRSS Tag-seq provides sensitive and reliable gene expression data and enables high-throughput expression profiling with relatively simple downstream analysis.