Fine-scale contemporary recombination variation and its fitness consequences in adaptively diverging stickleback fish.

Despite deep evolutionary conservation, recombination rates vary greatly across the genome and among individuals, sexes and populations. Yet the impact of this variation on adaptively diverging populations is not well understood. Here we characterized fine-scale recombination landscapes in an adaptively divergent pair of marine and freshwater populations of threespine stickleback from River Tyne, Scotland. Through whole-genome sequencing of large nuclear families, we identified the genomic locations of almost 50,000 crossovers and built recombination maps for marine, freshwater and hybrid individuals at a resolution of 3.8 kb. We used these maps to quantify the factors driving variation in recombination rates. We found strong heterochiasmy between sexes but also differences in recombination rates among ecotypes. Hybrids showed evidence of significant recombination suppression in overall map length and in individual loci. Recombination rates were lower not only within individual marine-freshwater-adaptive loci, but also between loci on the same chromosome, suggesting selection on linked gene 'cassettes'. Through temporal sampling along a natural hybrid zone, we found that recombinants showed traits associated with reduced fitness. Our results support predictions that divergence in cis-acting recombination modifiers, whose functions are disrupted in hybrids, may play an important role in maintaining differences among adaptively diverging populations.

Multi-omics analysis reveals the dynamic interplay between Vero host chromatin structure and function during vaccinia virus infection.

The genome folds into complex configurations and structures thought to profoundly impact its function. The intricacies of this dynamic structure-function relationship are not well understood particularly in the context of viral infection. To unravel this interplay, here we provide a comprehensive investigation of simultaneous host chromatin structural (via Hi-C and ATAC-seq) and functional changes (via RNA-seq) in response to vaccinia virus infection. Over time, infection significantly impacts global and local chromatin structure by increasing long-range intra-chromosomal interactions and B compartmentalization and by decreasing chromatin accessibility and inter-chromosomal interactions. Local accessibility changes are independent of broad-scale chromatin compartment exchange (~12% of the genome), underscoring potential independent mechanisms for global and local chromatin reorganization. While infection structurally condenses the host genome, there is nearly equal bidirectional differential gene expression. Despite global weakening of intra-TAD interactions, functional changes including downregulated immunity genes are associated with alterations in local accessibility and loop domain restructuring. Therefore, chromatin accessibility and local structure profiling provide impactful predictions for host responses and may improve development of efficacious anti-viral counter measures including the optimization of vaccine design.

Human Coronavirus Infection Reorganizes Spatial Genomic Architecture in Permissive Lung Cells.

Chromatin conformation capture followed by next-generation sequencing in combination with large-scale polymer simulations (4DHiC) produces detailed information on genomic loci interactions, allowing for the interrogation of 3D spatial genomic structures. Here, Hi-C data was acquired from the infection of fetal lung fibroblast (MRC5) cells with α-coronavirus 229E (CoV229E). Experimental Hi-C contact maps were used to determine viral-induced changes in genomic architecture over a 48-hour time period following viral infection, revealing substantial alterations in contacts within chromosomes and in contacts between different chromosomes. To gain further structural insight and quantify the underlying changes, we applied the 4DHiC polymer simulation method to reconstruct the 3D genomic structures and dynamics corresponding to the Hi-C maps. The models successfully reproduced experimental Hi-C data, including the changes in contacts induced by viral infection. Our 3D spatial simulations uncovered widespread chromatin restructuring, including increased chromosome compactness and A-B compartment mixing arising from infection. Our model also suggests increased spatial accessibility to regions containing interferon-stimulated genes upon infection with CoV229E, followed by chromatin restructuring at later time points, potentially inducing the migration of chromatin into more compact regions. This is consistent with previously observed suppression of gene expression. Our spatial genomics study provides a mechanistic structural basis for changes in chromosome architecture induced by coronavirus infection in lung cells.

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information.

BACKGROUND: Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility. RESULTS: Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships. CONCLUSIONS: Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

Genome-wide recombination map construction from single individuals using linked-read sequencing.

Meiotic recombination rates vary across the genome, often involving localized crossover hotspots and coldspots. Studying the molecular basis and mechanisms underlying this variation has been challenging due to the high cost and effort required to construct individualized genome-wide maps of recombination crossovers. Here we introduce a new method, called ReMIX, to detect crossovers from gamete DNA of a single individual using Illumina sequencing of 10X Genomics linked-read libraries. ReMIX reconstructs haplotypes and identifies the valuable rare molecules spanning crossover breakpoints, allowing quantification of the genomic location and intensity of meiotic recombination. Using a single mouse and stickleback fish, we demonstrate how ReMIX faithfully recovers recombination hotspots and landscapes that have previously been built using hundreds of offspring. ReMIX provides a high-resolution, high-throughput, and low-cost approach to quantify recombination variation across the genome, providing an exciting opportunity to study recombination among multiple individuals in diverse organisms.