Local selection signals in the genome of blue tits emphasize regulatory and neuronal evolution.

Understanding the genomic landscape of adaptation is central to understanding microevolution in wild populations. Genomic targets of selection and the underlying genomic mechanisms of adaptation can be elucidated by genome-wide scans for past selective sweeps or by scans for direct fitness associations. We sequenced and assembled 150 haplotypes of 75 blue tits (Cyanistes caeruleus) of a single Central European population by a linked-read technology. We used these genome data in combination with coalescent simulations (i) to estimate an historical effective population size of ~250,000, which recently declined to ~10,000, and (ii) to identify genome-wide distributed selective sweeps of beneficial variants probably originating from standing genetic variation (soft sweeps). The genes linked to these soft sweeps, but also those linked to hard sweeps based on new beneficial mutants, showed a significant enrichment for functions associated with gene expression and transcription regulation. This emphasizes the importance of regulatory evolution in the population's adaptive history. Soft sweeps were further enriched for genes related to axon and synapse development, indicating the significance of neuronal connectivity changes in the brain potentially linked to behavioural adaptations. A previous scan of heterozygosity-fitness correlations revealed a consistent negative effect on arrival date at the breeding site for a single microsatellite in the MDGA2 gene. Here, we used the haplotype structure around this microsatellite to explain the effect as a local and direct outbreeding effect of a gene involved in synapse development.

Linked-read sequencing enables haplotype-resolved resequencing at population scale.

The feasibility to sequence entire genomes of virtually any organism provides unprecedented insights into the evolutionary history of populations and species. Nevertheless, many population genomic inferences - including the quantification and dating of admixture, introgression and demographic events, and inference of selective sweeps - are still limited by the lack of high-quality haplotype information. The newest generation of sequencing technology now promises significant progress. To establish the feasibility of haplotype-resolved genome resequencing at population scale, we investigated properties of linked-read sequencing data of songbirds of the genus Oenanthe across a range of sequencing depths. Our results based on the comparison of downsampled (25×, 20×, 15×, 10×, 7×, and 5×) with high-coverage data (46-68×) of seven bird genomes mapped to a reference suggest that phasing contiguities and accuracies adequate for most population genomic analyses can be reached already with moderate sequencing effort. At 15× coverage, phased haplotypes span about 90% of the genome assembly, with 50% and 90% of phased sequences located in phase blocks longer than 1.25-4.6 Mb (N50) and 0.27-0.72 Mb (N90). Phasing accuracy reaches beyond 99% starting from 15× coverage. Higher coverages yielded higher contiguities (up to about 7 Mb/1 Mb [N50/N90] at 25× coverage), but only marginally improved phasing accuracy. Phase block contiguity improved with input DNA molecule length; thus, higher-quality DNA may help keeping sequencing costs at bay. In conclusion, even for organisms with gigabase-sized genomes like birds, linked-read sequencing at moderate depth opens an affordable avenue towards haplotype-resolved genome resequencing at population scale.

Automation of Spatial Transcriptomics library preparation to enable rapid and robust insights into spatial organization of tissues.

BACKGROUND: Interest in studying the spatial distribution of gene expression in tissues is rapidly increasing. Spatial Transcriptomics is a novel sequencing-based technology that generates high-throughput information on the distribution, heterogeneity and co-expression of cells in tissues. Unfortunately, manual preparation of high-quality sequencing libraries is time-consuming and subject to technical variability due to human error during manual pipetting, which results in sample swapping and the accidental introduction of batch effects. All these factors complicate the production and interpretation of biological datasets. RESULTS: We have integrated an Agilent Bravo Automated Liquid Handling Platform into the Spatial Transcriptomics workflow. Compared to the previously reported Magnatrix 8000+ automated protocol, this approach increases the number of samples processed per run, reduces sample preparation time by 35%, and minimizes batch effects between samples. The new approach is also shown to be highly accurate and almost completely free from technical variability between prepared samples. CONCLUSIONS: The new automated Spatial Transcriptomics protocol using the Agilent Bravo Automated Liquid Handling Platform rapidly generates high-quality Spatial Transcriptomics libraries. Given the wide use of the Agilent Bravo Automated Liquid Handling Platform in research laboratories and facilities, this will allow many researchers to quickly create robust Spatial Transcriptomics libraries.

Genome-wide evidence supports mitochondrial relationships and pervasive parallel phenotypic evolution in open-habitat chats.

In wheatears and related species ('open-habitat chats'), molecular phylogenetics has led to a comprehensively revised understanding of species relationships and species diversity. Phylogenetic analyses have suggested that, in many cases, phenotypic similarities do not reflect species' relationships, revealing traditionally defined genera as non-monophyletic. This led to the suggestion of pervasive parallel evolution of open-habitat chats' plumage coloration and ecological phenotypes. However, to date, the molecular evidence for the phylogenetic relationships among open-habitat chats is mainly limited to mitochondrial DNA. Here, we assessed whether the mitochondrial relationships are supported by genome-wide data. To this end, we reconstructed the species tree among 14 open-habitat chat taxa using multi-species coalescent analyses based on ~1'300 SNPs. Our results confirm previous ones based chiefly on mitochondrial DNA; notably the paraphyly of the Oenanthe lugens complex and the clustering of individual species formerly placed in the genera Cercomela and Myrmecocichla within Oenanthe. Since several variable morphological and ecological characteristics occur in multiple places across the open-habitat chat phylogeny, our study consolidates the evidence for pervasive parallel evolution in the plumage coloration and ecology of open-habitat chats.

Comprehensive structural variation genome map of individuals carrying complex chromosomal rearrangements.

Complex chromosomal rearrangements (CCRs) are rearrangements involving more than two chromosomes or more than two breakpoints. Whole genome sequencing (WGS) allows for outstanding high resolution characterization on the nucleotide level in unique sequences of such rearrangements, but problems remain for mapping breakpoints in repetitive regions of the genome, which are known to be prone to rearrangements. Hence, multiple complementary WGS experiments are sometimes needed to solve the structures of CCRs. We have studied three individuals with CCRs: Case 1 and Case 2 presented with de novo karyotypically balanced, complex interchromosomal rearrangements (46,XX,t(2;8;15)(q35;q24.1;q22) and 46,XY,t(1;10;5)(q32;p12;q31)), and Case 3 presented with a de novo, extremely complex intrachromosomal rearrangement on chromosome 1. Molecular cytogenetic investigation revealed cryptic deletions in the breakpoints of chromosome 2 and 8 in Case 1, and on chromosome 10 in Case 2, explaining their clinical symptoms. In Case 3, 26 breakpoints were identified using WGS, disrupting five known disease genes. All rearrangements were subsequently analyzed using optical maps, linked-read WGS, and short-read WGS. In conclusion, we present a case series of three unique de novo CCRs where we by combining the results from the different technologies fully solved the structure of each rearrangement. The power in combining short-read WGS with long-molecule sequencing or optical mapping in these unique de novo CCRs in a clinical setting is demonstrated.

Transcriptomics and methylomics of CD4-positive T cells in arsenic-exposed women.

Arsenic, a carcinogen with immunotoxic effects, is a common contaminant of drinking water and certain food worldwide. We hypothesized that chronic arsenic exposure alters gene expression, potentially by altering DNA methylation of genes encoding central components of the immune system. We therefore analyzed the transcriptomes (by RNA sequencing) and methylomes (by target-enrichment next-generation sequencing) of primary CD4-positive T cells from matched groups of four women each in the Argentinean Andes, with fivefold differences in urinary arsenic concentrations (median concentrations of urinary arsenic in the lower- and high-arsenic groups: 65 and 276 µg/l, respectively). Arsenic exposure was associated with genome-wide alterations of gene expression; principal component analysis indicated that the exposure explained 53% of the variance in gene expression among the top variable genes and 19% of 28,351 genes were differentially expressed (false discovery rate <0.05) between the exposure groups. Key genes regulating the immune system, such as tumor necrosis factor alpha and interferon gamma, as well as genes related to the NF-kappa-beta complex, were significantly downregulated in the high-arsenic group. Arsenic exposure was associated with genome-wide DNA methylation; the high-arsenic group had 3% points higher genome-wide full methylation (>80% methylation) than the lower-arsenic group. Differentially methylated regions that were hyper-methylated in the high-arsenic group showed enrichment for immune-related gene ontologies that constitute the basic functions of CD4-positive T cells, such as isotype switching and lymphocyte activation and differentiation. In conclusion, chronic arsenic exposure from drinking water was related to changes in the transcriptome and methylome of CD4-positive T cells, both genome wide and in specific genes, supporting the hypothesis that arsenic causes immunotoxicity by interfering with gene expression and regulation.

Hierarchical molecular tagging to resolve long continuous sequences by massively parallel sequencing.

Here we demonstrate the use of short-read massive sequencing systems to in effect achieve longer read lengths through hierarchical molecular tagging. We show how indexed and PCR-amplified targeted libraries are degraded, sub-sampled and arrested at timed intervals to achieve pools of differing average length, each of which is indexed with a new tag. By this process, indices of sample origin, molecular origin, and degree of degradation is incorporated in order to achieve a nested hierarchical structure, later to be utilized in the data processing to order the reads over a longer distance than the sequencing system originally allows. With this protocol we show how continuous regions beyond 3000 bp can be decoded by an Illumina sequencing system, and we illustrate the potential applications by calling variants of the lambda genome, analysing TP53 in cancer cell lines, and targeting a variable canine mitochondrial region.