Litter Commensal Bacteria Can Limit the Horizontal Gene Transfer of Antimicrobial Resistance to Salmonella in Chickens.

Fostering a "balanced" gut microbiome through the administration of beneficial microbes that can competitively exclude pathogens has gained a lot of attention and use in human and animal medicine. However, little is known about how microbes affect the horizontal gene transfer of antimicrobial resistance (AMR). To shed more light on this question, we challenged neonatal broiler chicks raised on reused broiler chicken litter-a complex environment made up of decomposing pine shavings, feces, uric acid, feathers, and feed-with Salmonella enterica serovar Heidelberg (S. Heidelberg), a model pathogen. Neonatal chicks challenged with S. Heidelberg and raised on reused litter were more resistant to S. Heidelberg cecal colonization than chicks grown on fresh litter. Furthermore, chicks grown on reused litter were at a lower risk of colonization with S. Heidelberg strains that encoded AMR on IncI1 plasmids. We used 16S rRNA gene sequencing and shotgun metagenomics to show that the major difference between chicks grown on fresh litter and those grown on reused litter was the microbiome harbored in the litter and ceca. The microbiome of reused litter samples was more uniform and enriched in functional pathways related to the biosynthesis of organic and antimicrobial molecules than that in fresh litter samples. We found that Escherichia coli was the main reservoir of plasmids encoding AMR and that the IncI1 plasmid was maintained at a significantly lower copy per cell in reused litter compared to fresh litter. These findings support the notion that commensal bacteria play an integral role in the horizontal transfer of plasmids encoding AMR to pathogens like Salmonella. IMPORTANCE Antimicrobial resistance spread is a worldwide health challenge, stemming in large part from the ability of microorganisms to share their genetic material through horizontal gene transfer. To address this issue, many countries and international organizations have adopted a One Health approach to curtail the proliferation of antimicrobial-resistant bacteria. This includes the removal and reduction of antibiotics used in food animal production and the development of alternatives to antibiotics. However, there is still a significant knowledge gap in our understanding of how resistance spreads in the absence of antibiotic selection and the role commensal bacteria play in reducing antibiotic resistance transfer. In this study, we show that commensal bacteria play a key role in reducing the horizontal gene transfer of antibiotic resistance to Salmonella, provide the identity of the bacterial species that potentially perform this function in broiler chickens, and also postulate the mechanism involved.

Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs.

BACKGROUND: The genetic mechanisms that underlie phenotypic differentiation in breeding animals have important implications in evolutionary biology and agriculture. However, the contribution of cis-regulatory variants to pig phenotypes is poorly understood. Therefore, our aim was to elucidate the molecular mechanisms by which non-coding variants cause phenotypic differences in pigs by combining evolutionary biology analyses and functional genomics. RESULTS: We obtained a high-resolution phased chromosome-scale reference genome with a contig N50 of 18.03 Mb for the Luchuan pig breed (a representative eastern breed) and profiled potential selective sweeps in eastern and western pigs by resequencing the genomes of 234 pigs. Multi-tissue transcriptome and chromatin accessibility analyses of these regions suggest that tissue-specific selection pressure is mediated by promoters and distal cis-regulatory elements. Promoter variants that are associated with increased expression of the lysozyme (LYZ) gene in the small intestine might enhance the immunity of the gastrointestinal tract and roughage tolerance in pigs. In skeletal muscle, an enhancer-modulating single-nucleotide polymorphism that is associated with up-regulation of the expression of the troponin C1, slow skeletal and cardiac type (TNNC1) gene might increase the proportion of slow muscle fibers and affect meat quality. CONCLUSIONS: Our work sheds light on the molecular mechanisms by which non-coding variants shape phenotypic differences in pigs and provides valuable resources and novel perspectives to dissect the role of gene regulatory evolution in animal domestication and breeding.

Repeat elements organise 3D genome structure and mediate transcription in the filamentous fungus Epichloë festucae.

Structural features of genomes, including the three-dimensional arrangement of DNA in the nucleus, are increasingly seen as key contributors to the regulation of gene expression. However, studies on how genome structure and nuclear organisation influence transcription have so far been limited to a handful of model species. This narrow focus limits our ability to draw general conclusions about the ways in which three-dimensional structures are encoded, and to integrate information from three-dimensional data to address a broader gamut of biological questions. Here, we generate a complete and gapless genome sequence for the filamentous fungus, Epichloë festucae. We use Hi-C data to examine the three-dimensional organisation of the genome, and RNA-seq data to investigate how Epichloë genome structure contributes to the suite of transcriptional changes needed to maintain symbiotic relationships with the grass host. Our results reveal a genome in which very repeat-rich blocks of DNA with discrete boundaries are interspersed by gene-rich sequences that are almost repeat-free. In contrast to other species reported to date, the three-dimensional structure of the genome is anchored by these repeat blocks, which act to isolate transcription in neighbouring gene-rich regions. Genes that are differentially expressed in planta are enriched near the boundaries of these repeat-rich blocks, suggesting that their three-dimensional orientation partly encodes and regulates the symbiotic relationship formed by this organism.

Differential paralog divergence modulates genome evolution across yeast species.

Evolutionary outcomes depend not only on the selective forces acting upon a species, but also on the genetic background. However, large timescales and uncertain historical selection pressures can make it difficult to discern such important background differences between species. Experimental evolution is one tool to compare evolutionary potential of known genotypes in a controlled environment. Here we utilized a highly reproducible evolutionary adaptation in Saccharomyces cerevisiae to investigate whether experimental evolution of other yeast species would select for similar adaptive mutations. We evolved populations of S. cerevisiae, S. paradoxus, S. mikatae, S. uvarum, and interspecific hybrids between S. uvarum and S. cerevisiae for ~200-500 generations in sulfate-limited continuous culture. Wild-type S. cerevisiae cultures invariably amplify the high affinity sulfate transporter gene, SUL1. However, while amplification of the SUL1 locus was detected in S. paradoxus and S. mikatae populations, S. uvarum cultures instead selected for amplification of the paralog, SUL2. We measured the relative fitness of strains bearing deletions and amplifications of both SUL genes from different species, confirming that, converse to S. cerevisiae, S. uvarum SUL2 contributes more to fitness in sulfate limitation than S. uvarum SUL1. By measuring the fitness and gene expression of chimeric promoter-ORF constructs, we were able to delineate the cause of this differential fitness effect primarily to the promoter of S. uvarum SUL1. Our data show evidence of differential sub-functionalization among the sulfate transporters across Saccharomyces species through recent changes in noncoding sequence. Furthermore, these results show a clear example of how such background differences due to paralog divergence can drive changes in genome evolution.

Neocentromeres Provide Chromosome Segregation Accuracy and Centromere Clustering to Multiple Loci along a Candida albicans Chromosome.

Assembly of kinetochore complexes, involving greater than one hundred proteins, is essential for chromosome segregation and genome stability. Neocentromeres, or new centromeres, occur when kinetochores assemble de novo, at DNA loci not previously associated with kinetochore proteins, and they restore chromosome segregation to chromosomes lacking a functional centromere. Neocentromeres have been observed in a number of diseases and may play an evolutionary role in adaptation or speciation. However, the consequences of neocentromere formation on chromosome missegregation rates, gene expression, and three-dimensional (3D) nuclear structure are not well understood. Here, we used Candida albicans, an organism with small, epigenetically-inherited centromeres, as a model system to study the functions of twenty different neocentromere loci along a single chromosome, chromosome 5. Comparison of neocentromere properties relative to native centromere functions revealed that all twenty neocentromeres mediated chromosome segregation, albeit to different degrees. Some neocentromeres also caused reduced levels of transcription from genes found within the neocentromere region. Furthermore, like native centromeres, neocentromeres clustered in 3D with active/functional centromeres, indicating that formation of a new centromere mediates the reorganization of 3D nuclear architecture. This demonstrates that centromere clustering depends on epigenetically defined function and not on the primary DNA sequence, and that neocentromere function is independent of its distance from the native centromere position. Together, the results show that a neocentromere can form at many loci along a chromosome and can support the assembly of a functional kinetochore that exhibits native centromere functions including chromosome segregation accuracy and centromere clustering within the nucleus.

Identification of a novel interspecific hybrid yeast from a metagenomic spontaneously inoculated beer sample using Hi-C.

Interspecific hybridization is a common mechanism enabling genetic diversification and adaptation; however, the detection of hybrid species has been quite difficult. The identification of microbial hybrids is made even more complicated, as most environmental microbes are resistant to culturing and must be studied in their native mixed communities. We have previously adapted the chromosome conformation capture method Hi-C to the assembly of genomes from mixed populations. Here, we show the method's application in assembling genomes directly from an uncultured, mixed population from a spontaneously inoculated beer sample. Our assembly method has enabled us to de-convolute four bacterial and four yeast genomes from this sample, including a putative yeast hybrid. Downstream isolation and analysis of this hybrid confirmed its genome to consist of Pichia membranifaciens and that of another related, but undescribed, yeast. Our work shows that Hi-C-based metagenomic methods can overcome the limitation of traditional sequencing methods in studying complex mixtures of genomes. Copyright © 2017 John Wiley & Sons, Ltd.

The fitness consequences of aneuploidy are driven by condition-dependent gene effects.

Aneuploidy is a hallmark of tumor cells, and yet the precise relationship between aneuploidy and a cell's proliferative ability, or cellular fitness, has remained elusive. In this study, we have combined a detailed analysis of aneuploid clones isolated from laboratory-evolved populations of Saccharomyces cerevisiae with a systematic, genome-wide screen for the fitness effects of telomeric amplifications to address the relationship between aneuploidy and cellular fitness. We found that aneuploid clones rise to high population frequencies in nutrient-limited evolution experiments and show increased fitness relative to wild type. Direct competition experiments confirmed that three out of four aneuploid events isolated from evolved populations were themselves sufficient to improve fitness. To expand the scope beyond this small number of exemplars, we created a genome-wide collection of >1,800 diploid yeast strains, each containing a different telomeric amplicon (Tamp), ranging in size from 0.4 to 1,000 kb. Using pooled competition experiments in nutrient-limited chemostats followed by high-throughput sequencing of strain-identifying barcodes, we determined the fitness effects of these >1,800 Tamps under three different conditions. Our data revealed that the fitness landscape explored by telomeric amplifications is much broader than that explored by single-gene amplifications. As also observed in the evolved clones, we found the fitness effects of most Tamps to be condition specific, with a minority showing common effects in all three conditions. By integrating our data with previous work that examined the fitness effects of single-gene amplifications genome-wide, we found that a small number of genes within each Tamp are centrally responsible for each Tamp's fitness effects. Our genome-wide Tamp screen confirmed that telomeric amplifications identified in laboratory-evolved populations generally increased fitness. Our results show that Tamps are mutations that produce large, typically condition-dependent changes in fitness that are important drivers of increased fitness in asexually evolving populations.

An underlying mechanism for the increased mutagenesis of lagging-strand genes in Bacillus subtilis.

We previously reported that lagging-strand genes accumulate mutations faster than those encoded on the leading strand in Bacillus subtilis. Although we proposed that orientation-specific encounters between replication and transcription underlie this phenomenon, the mechanism leading to the increased mutagenesis of lagging-strand genes remained unknown. Here, we report that the transcription-dependent and orientation-specific differences in mutation rates of genes require the B. subtilis Y-family polymerase, PolY1 (yqjH). We find that without PolY1, association of the replicative helicase, DnaC, and the recombination protein, RecA, with lagging-strand genes increases in a transcription-dependent manner. These data suggest that PolY1 promotes efficient replisome progression through lagging-strand genes, thereby reducing potentially detrimental breaks and single-stranded DNA at these loci. Y-family polymerases can alleviate potential obstacles to replisome progression by facilitating DNA lesion bypass, extension of D-loops, or excision repair. We find that the nucleotide excision repair (NER) proteins UvrA, UvrB, and UvrC, but not RecA, are required for transcription-dependent asymmetry in mutation rates of genes in the two orientations. Furthermore, we find that the transcription-coupling repair factor Mfd functions in the same pathway as PolY1 and is also required for increased mutagenesis of lagging-strand genes. Experimental and SNP analyses of B. subtilis genomes show mutational footprints consistent with these findings. We propose that the interplay between replication and transcription increases lesion susceptibility of, specifically, lagging-strand genes, activating an Mfd-dependent error-prone NER mechanism. We propose that this process, at least partially, underlies the accelerated evolution of lagging-strand genes.

Improvement of the Threespine Stickleback Genome Using a Hi-C-Based Proximity-Guided Assembly.

Scaffolding genomes into complete chromosome assemblies remains challenging even with the rapidly increasing sequence coverage generated by current next-generation sequence technologies. Even with scaffolding information, many genome assemblies remain incomplete. The genome of the threespine stickleback (Gasterosteus aculeatus), a fish model system in evolutionary genetics and genomics, is not completely assembled despite scaffolding with high-density linkage maps. Here, we first test the ability of a Hi-C based proximity-guided assembly (PGA) to perform a de novo genome assembly from relatively short contigs. Using Hi-C based PGA, we generated complete chromosome assemblies from a distribution of short contigs (20-100 kb). We found that 96.40% of contigs were correctly assigned to linkage groups (LGs), with ordering nearly identical to the previous genome assembly. Using available bacterial artificial chromosome (BAC) end sequences, we provide evidence that some of the few discrepancies between the Hi-C assembly and the existing assembly are due to structural variation between the populations used for the 2 assemblies or errors in the existing assembly. This Hi-C assembly also allowed us to improve the existing assembly, assigning over 60% (13.35 Mb) of the previously unassigned (~21.7 Mb) contigs to LGs. Together, our results highlight the potential of the Hi-C based PGA method to be used in combination with short read data to perform relatively inexpensive de novo genome assemblies. This approach will be particularly useful in organisms in which it is difficult to perform linkage mapping or to obtain high molecular weight DNA required for other scaffolding methods.

Accurate identification of centromere locations in yeast genomes using Hi-C.

Centromeres are essential for proper chromosome segregation. Despite extensive research, centromere locations in yeast genomes remain difficult to infer, and in most species they are still unknown. Recently, the chromatin conformation capture assay, Hi-C, has been re-purposed for diverse applications, including de novo genome assembly, deconvolution of metagenomic samples and inference of centromere locations. We describe a method, Centurion, that jointly infers the locations of all centromeres in a single genome from Hi-C data by exploiting the centromeres' tendency to cluster in three-dimensional space. We first demonstrate the accuracy of Centurion in identifying known centromere locations from high coverage Hi-C data of budding yeast and a human malaria parasite. We then use Centurion to infer centromere locations in 14 yeast species. Across all microbes that we consider, Centurion predicts 89% of centromeres within 5 kb of their known locations. We also demonstrate the robustness of the approach in datasets with low sequencing depth. Finally, we predict centromere coordinates for six yeast species that currently lack centromere annotations. These results show that Centurion can be used for centromere identification for diverse species of yeast and possibly other microorganisms.

GC-rich DNA elements enable replication origin activity in the methylotrophic yeast Pichia pastoris.

The well-studied DNA replication origins of the model budding and fission yeasts are A/T-rich elements. However, unlike their yeast counterparts, both plant and metazoan origins are G/C-rich and are associated with transcription start sites. Here we show that an industrially important methylotrophic budding yeast, Pichia pastoris, simultaneously employs at least two types of replication origins--a G/C-rich type associated with transcription start sites and an A/T-rich type more reminiscent of typical budding and fission yeast origins. We used a suite of massively parallel sequencing tools to map and dissect P. pastoris origins comprehensively, to measure their replication dynamics, and to assay the global positioning of nucleosomes across the genome. Our results suggest that some functional overlap exists between promoter sequences and G/C-rich replication origins in P. pastoris and imply an evolutionary bifurcation of the modes of replication initiation.

High-resolution mapping, characterization, and optimization of autonomously replicating sequences in yeast.

DNA replication origins are necessary for the duplication of genomes. In addition, plasmid-based expression systems require DNA replication origins to maintain plasmids efficiently. The yeast autonomously replicating sequence (ARS) assay has been a valuable tool in dissecting replication origin structure and function. However, the dearth of information on origins in diverse yeasts limits the availability of efficient replication origin modules to only a handful of species and restricts our understanding of origin function and evolution. To enable rapid study of origins, we have developed a sequencing-based suite of methods for comprehensively mapping and characterizing ARSs within a yeast genome. Our approach finely maps genomic inserts capable of supporting plasmid replication and uses massively parallel deep mutational scanning to define molecular determinants of ARS function with single-nucleotide resolution. In addition to providing unprecedented detail into origin structure, our data have allowed us to design short, synthetic DNA sequences that retain maximal ARS function. These methods can be readily applied to understand and modulate ARS function in diverse systems.

Characterization of a panARS-based episomal vector in the methylotrophic yeast Pichia pastoris for recombinant protein production and synthetic biology applications.

BACKGROUND: Recombinant protein production in the methylotrophic yeast Pichia pastoris largely relies on integrative vectors. Although the stability of integrated expression cassettes is well appreciated for most applications, the availability of reliable episomal vectors for this host would represent a useful tool to expedite cloning and high-throughput screening, ameliorating also the relatively high clonal variability reported in transformants from integrative vectors caused by off-target integration in the P. pastoris genome. Recently, heterologous and endogenous autonomously replicating sequences (ARS) were identified in P. pastoris by genome mining, opening the possibility of expanding the available toolbox to include efficient episomal plasmids. The aim of this technical report is to validate a 452-bp sequence ("panARS") in context of P. pastoris expression vectors, and to compare their performance to classical integrative plasmids. Moreover, we aimed to test if such episomal vectors would be suitable to sustain in vivo recombination, using fragments for transformation, directly in P. pastoris cells. RESULTS: A panARS-based episomal vector was evaluated using blue fluorescent protein (BFP) as a reporter gene. Normalized fluorescence from colonies carrying panARS-BFP outperformed the level of signal obtained from integrative controls by several-fold, whereas endogenous sequences, identified from the P. pastoris genome, were not as efficient in terms of protein production. At the single cell level, panARS-BFP clones showed lower interclonal variability but higher intraclonal variation compared to their integrative counterparts, supporting the idea that heterologous protein production could benefit from episomal plasmids. Finally, efficiency of 2-fragment and 3-fragment in vivo recombination was tested using varying lengths of overlapping regions and molar ratios between fragments. Upon optimization, minimal background was obtained for in vivo assembled vectors, suggesting this could be a quick and efficient method to generate of episomal plasmids of interest. CONCLUSIONS: An expression vector based on the panARS sequence was shown to outperform its integrative counterparts in terms of protein productivity and interclonal variability, facilitating recombinant protein expression and screening. Using optimized fragment lengths and ratios, it was possible to perform reliable in vivo recombination of fragments in P. pastoris. Taken together, these results support the applicability of panARS episomal vectors for synthetic biology approaches.

A viable allele of Mcm4 causes chromosome instability and mammary adenocarcinomas in mice.

Mcm4 (minichromosome maintenance-deficient 4 homolog) encodes a subunit of the MCM2-7 complex (also known as MCM2-MCM7), the replication licensing factor and presumptive replicative helicase. Here, we report that the mouse chromosome instability mutation Chaos3 (chromosome aberrations occurring spontaneously 3), isolated in a forward genetic screen, is a viable allele of Mcm4. Mcm4(Chaos3) encodes a change in an evolutionarily invariant amino acid (F345I), producing an apparently destabilized MCM4. Saccharomyces cerevisiae strains that we engineered to contain a corresponding allele (resulting in an F391I change) showed a classical minichromosome loss phenotype. Whereas homozygosity for a disrupted Mcm4 allele (Mcm4(-)) caused preimplantation lethality, Mcm(Chaos3/-) embryos died late in gestation, indicating that Mcm4(Chaos3) is hypomorphic. Mutant embryonic fibroblasts were highly susceptible to chromosome breaks induced by the DNA replication inhibitor aphidicolin. Most notably, >80% of Mcm4(Chaos3/Chaos3) females succumbed to mammary adenocarcinomas with a mean latency of 12 months. These findings suggest that hypomorphic alleles of the genes encoding the subunits of the MCM2-7 complex may increase breast cancer risk.

An autonomously replicating sequence for use in a wide range of budding yeasts.

The initiation of DNA replication at replication origins is essential for the duplication of genomes. In yeast, the autonomously replicating sequence (ARS) property of replication origins is necessary for the stable maintenance of episomal plasmids. However, because the sequence determinants of ARS function differ among yeast species, current ARS modules are limited for use to a subset of yeasts. Here, we describe a short ARS sequence that functions in at least 10 diverse species of budding yeast. These include, but are not limited to members of the Saccharomyces, Lachancea, Kluyveromyces, and Pichia (Komagataella) genera spanning over 500 million years of evolution. In addition to its wide species range, this ARS and an optimized derivative confer improved plasmid stability relative to other currently used ARS modules.

Analysis and culturing of the prototypic crAssphage reveals a phage-plasmid lifestyle.

The prototypic crAssphage (Carjivirus communis) is one of the most abundant, prevalent, and persistent gut bacteriophages, yet it remains uncultured and its lifestyle uncharacterized. For the last decade, crAssphage has escaped plaque-dependent culturing efforts, leading us to investigate alternative lifestyles that might explain its widespread success. Through genomic analyses and culturing, we find that crAssphage uses a phage-plasmid lifestyle to persist extrachromosomally. Plasmid-related genes are more highly expressed than those implicated in phage maintenance. Leveraging this finding, we use a plaque-free culturing approach to measure crAssphage replication in culture with Phocaeicola vulgatus, Phocaeicola dorei, and Bacteroides stercoris, revealing a broad host range. We demonstrate that crAssphage persists with its hosts in culture without causing major cell lysis events or integrating into host chromosomes. The ability to switch between phage and plasmid lifestyles within a wide range of hosts contributes to the prolific nature of crAssphage in the human gut microbiome.

Proxi-RIMS-seq2 applied to native microbiomes uncovers hundreds of known and novel m5C methyltransferase specificities.

Methylation patterns in bacteria can be used to study Restriction-Modification (RM) or other defense systems with novel properties. While m4C and m6A methylation is well characterized mainly through PacBio sequencing, the landscape of m5C methylation is under-characterized. To bridge this gap, we performed RIMS-seq2 on microbiomes composed of resolved assemblies of distinct genomes through proximity ligation. This high-throughput approach enables the identification of m5C methylated motifs and links them to cognate methyltransferases directly on native microbiomes without the need to isolate bacterial strains. Methylation patterns can also be identified on viral DNA and compared to host DNA, strengthening evidence for virus-host interaction. Applied to three different microbiomes, the method unveils over 1900 motifs that were deposited in REBASE. The motifs include a novel 8-base recognition site (CATm5CGATG) that was experimentally validated by characterizing its cognate methyltransferase. Our findings suggest that microbiomes harbor arrays of untapped m5C methyltransferase specificities, providing insights to bacterial biology and biotechnological applications.

Evaluation of Acute Myeloid Leukemia Genomes using Genomic Proximity Mapping.

Background: Cytogenetic analysis encompasses a suite of standard-of-care diagnostic testing methods that is routinely applied in cases of acute myeloid leukemia (AML) to assess chromosomal changes that are clinically relevant for risk classification and treatment decisions. Objective: In this study, we assess the use of Genomic Proximity Mapping (GPM) for cytogenomic analysis of AML diagnostic specimens for detection of cytogenetic risk variants included in the European Leukemia Network (ELN) risk stratification guidelines. Methods: Archival patient samples (N=48) from the Fred Hutchinson Cancer Center leukemia bank with historical clinical cytogenetic data were processed for GPM and analyzed with the CytoTerra® cloud-based analysis platform. Results: GPM showed 100% concordance for all specific variants that have associated impacts on risk stratification as defined by ELN 2022 criteria, and a 72% concordance rate when considering all variants reported by the FH cytogenetic lab. GPM identified 39 additional variants, including variants of known clinical impact, not observed by cytogenetics. Conclusions: GPM is an effective solution for the evaluation of known AML-associated risk variants and a source for biomarker discovery.

A comprehensive genome-wide map of autonomously replicating sequences in a naive genome.

Eukaryotic chromosomes initiate DNA synthesis from multiple replication origins. The machinery that initiates DNA synthesis is highly conserved, but the sites where the replication initiation proteins bind have diverged significantly. Functional comparative genomics is an obvious approach to study the evolution of replication origins. However, to date, the Saccharomyces cerevisiae replication origin map is the only genome map available. Using an iterative approach that combines computational prediction and functional validation, we have generated a high-resolution genome-wide map of DNA replication origins in Kluyveromyces lactis. Unlike other yeasts or metazoans, K. lactis autonomously replicating sequences (KlARSs) contain a 50 bp consensus motif suggestive of a dimeric structure. This motif is necessary and largely sufficient for initiation and was used to dependably identify 145 of the up to 156 non-repetitive intergenic ARSs projected for the K. lactis genome. Though similar in genome sizes, K. lactis has half as many ARSs as its distant relative S. cerevisiae. Comparative genomic analysis shows that ARSs in K. lactis and S. cerevisiae preferentially localize to non-syntenic intergenic regions, linking ARSs with loci of accelerated evolutionary change.

Integrating chromatin conformation information in a self-supervised learning model improves metagenome binning.

Metagenome binning is a key step, downstream of metagenome assembly, to group scaffolds by their genome of origin. Although accurate binning has been achieved on datasets containing multiple samples from the same community, the completeness of binning is often low in datasets with a small number of samples due to a lack of robust species co-abundance information. In this study, we exploited the chromatin conformation information obtained from Hi-C sequencing and developed a new reference-independent algorithm, Metagenome Binning with Abundance and Tetra-nucleotide frequencies-Long Range (metaBAT-LR), to improve the binning completeness of these datasets. This self-supervised algorithm builds a model from a set of high-quality genome bins to predict scaffold pairs that are likely to be derived from the same genome. Then, it applies these predictions to merge incomplete genome bins, as well as recruit unbinned scaffolds. We validated metaBAT-LR's ability to bin-merge and recruit scaffolds on both synthetic and real-world metagenome datasets of varying complexity. Benchmarking against similar software tools suggests that metaBAT-LR uncovers unique bins that were missed by all other methods. MetaBAT-LR is open-source and is available at https://bitbucket.org/project-metabat/metabat-lr.