Scientific highlights

Africa's black (Diceros bicornis) and white (Ceratotherium simum) rhinoceros are closely related sister-taxa that evolved highly divergent obligate browsing and grazing feeding strategies. Although their precursor species D. praecox and C. mauritanicum appear in the fossil record ∼5.2 million years ago (Ma), by 4 Ma both were still mixed feeders, and were even spatio-temporally sympatric at several Pliocene sites in what is today Africa's Rift Valley. Here, we ask whether or not D. praecox and C. mauritanicum were reproductively isolated when they came into Pliocene secondary contact. We sequenced and de novo assembled the first annotated black rhinoceros reference genome, and compared it with available genomes of other black and white rhinoceros. We show that ancestral gene flow between D. praecox and C. mauritanicum ceased sometime between 3.3 and 4.1 Ma, despite conventional methods for the detection of gene flow from whole genome data returning false positive signatures of recent interspecific migration due to incomplete lineage sorting. We propose that ongoing Pliocene genetic exchange, for up to 2 million years after initial divergence, could have potentially hindered the development of obligate feeding strategies until both species were fully reproductively isolated, but that the more severe and shifting palaeoclimate of the early Pleistocene was likely the ultimate driver of ecological specialisation in African rhinoceros.

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 (25x, 20x, 15x, 10x, 7x, and 5x) with high-coverage data (46-68x) 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 15x coverage, phased haplotypes span about 90% of the genome assembly, with 50 and 90 percent 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 15x coverage. Higher coverages yielded higher contiguities (up to about 7 Mb/1Mb (N50/N90) at 25x 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.

Pathogen-mediated balancing selection is commonly considered to play an important role in the maintenance of genetic diversity, in particular in immune genes. However, the factors that may influence which immune genes are the targets of such selection are largely unknown. To address this, here we focus on Pattern Recognition Receptor (PRR) signalling pathways, which play a key role in innate immunity. We used whole-genome resequencing data from a population of bank voles (Myodes glareolus) to test for associations between balancing selection, pleiotropy and gene function in a set of 123 PRR signalling pathway genes. To investigate the effect of gene function, we compared genes encoding (a) receptors for microbial ligands versus downstream signalling proteins, and (b) receptors recognizing components of microbial cell walls, flagella and capsids versus receptors recognizing features of microbial nucleic acids. Analyses based on the nucleotide diversity of full coding sequences showed that balancing selection primarily targeted receptor genes with a low degree of pleiotropy. Moreover, genes encoding receptors recognizing components of microbial cell walls etc. were more important targets of balancing selection than receptors recognizing nucleic acids. Tests for localized signatures of balancing selection in coding and noncoding sequences showed that such signatures were mostly located in introns, and more evenly distributed among different functional categories of PRR pathway genes. The finding that signatures of balancing selection in full coding sequences primarily occur in receptor genes, in particular those encoding receptors for components of microbial cell walls etc., is consistent with the idea that coevolution between hosts and pathogens is an important cause of balancing selection on immune genes.

The standardization, portability and reproducibility of analysis pipelines are key issues within the bioinformatics community. Most bioinformatics pipelines are designed for use on-premises; as a result, the associated software dependencies and execution logic are likely to be tightly coupled with proprietary computing environments. This can make it difficult or even impossible for others to reproduce the ensuing results, which is a fundamental requirement for the validation of scientific findings. Here, we introduce the nf-core framework as a means for the development of collaborative, peer-reviewed, best-practice analysis pipelines (Fig. 1). All nf-core pipelines are written in Nextflow and so inherit the ability to be executed on most computational infrastructures, as well as having native support for container technologies such as Docker and Singularity. The nf-core community (Supplementary Fig. 1) has developed a suite of tools that automate pipeline creation, testing, deployment and synchronization. Our goal is to provide a framework for high-quality bioinformatics pipelines that can be used across all institutions and research facilities.

As the only endemic neotropical parrot to have recently lived in the northern hemisphere, the Carolina parakeet (Conuropsis carolinensis) was an iconic North American bird. The last surviving specimen died in the Cincinnati Zoo in 1918 [1]. The cause of its extinction remains contentious: besides excessive mortality associated to habitat destruction and active hunting, their survival could have been negatively affected by its range having become increasingly patchy [2] or by the exposure to poultry pathogens [3, 4]. In addition, the Carolina parakeet showed a predilection for cockleburs, an herbaceous plant that contains a powerful toxin, carboxyatractyloside, or CAT [5], which did not seem to affect them but made the birds notoriously toxic to most predators [3]. To explore the demographic history of this bird, we generated the complete genomic sequence of a preserved specimen held in a private collection in Espinelves (Girona, Spain), as well as of a close extant relative, Aratinga solstitialis. We identified two non-synonymous genetic changes in two highly conserved proteins known to interact with CAT that could underlie a specific dietary adaptation to this toxin. Our genomic analyses did not reveal evidence of a dramatic past demographic decline in the Carolina parakeet; also, its genome did not exhibit the long runs of homozygosity that are signals of recent inbreeding and are typically found in endangered species. As such, our results suggest its extinction was an abrupt process and thus likely solely attributable to human causes.

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

Long-read sequencing can resolve regions of the genome that are inaccessible to short reads, and therefore are ideal for genome-gap closure, solving structural rearrangements and sequencing through repetitive elements. Here we introduce the Xdrop technology: a novel microfluidic-based system that allows for targeted enrichment of long DNA molecules starting from only a few nanograms of DNA. Xdrop is based on isolation of long DNA fragments in millions of droplets, where the droplets containing a target sequence of interest are fluorescently labeled and sorted using flow cytometry. The final product from the Xdrop procedure is an enriched population of long DNA molecules that can be investigated by sequencing. To demonstrate the capability of Xdrop, we performed enrichment of the human papilloma virus (HPV) 18 integrated in the genome of human HeLa cells. Analysis of the sequencing reads resolved three HPV18-chr8 integrations at base pair resolution, and the captured fragments extended up to 30 kb into the human genome at the integration sites. Further, we enriched the complete TP53 locus in a leukemia cell line and could successfully phase co-existing mutations using PacBio sequencing. In summary, our results show that Xdrop is an efficient enrichment technology for studying complex genomic regions. This article is protected by copyright. All rights reserved.

Population health research is increasingly focused on the genetic determinants of healthy ageing, but there is no public resource of whole genome sequences and phenotype data from healthy elderly individuals. Here we describe the first release of the Medical Genome Reference Bank (MGRB), comprising whole genome sequence and phenotype of 2570 elderly Australians depleted for cancer, cardiovascular disease, and dementia. We analyse the MGRB for single-nucleotide, indel and structural variation in the nuclear and mitochondrial genomes. MGRB individuals have fewer disease-associated common and rare germline variants, relative to both cancer cases and the gnomAD and UK Biobank cohorts, consistent with risk depletion. Age-related somatic changes are correlated with grip strength in men, suggesting blood-derived whole genomes may also provide a biologic measure of age-related functional deterioration. The MGRB provides a broadly applicable reference cohort for clinical genetics and genomic association studies, and for understanding the genetics of healthy ageing.

Genomic rearrangements associated with speciation often result in chromosome number variation among closely related species. Malassezia species show variable karyotypes ranging between 6 and 9 chromosomes. Here, we experimentally identified all 8 centromeres in M. sympodialis as 3 to 5 kb long kinetochore-bound regions spanning an AT-rich core and depleted of the canonical histone H3. Centromeres of similar sequence features were identified as CENP-A-rich regions in Malassezia furfur with 7 chromosomes, and histone H3 depleted regions in Malassezia slooffiae and Malassezia globosa with 9 chromosomes each. Analysis of synteny conservation across centromeres with newly generated chromosome-level genome assemblies suggests two distinct mechanisms of chromosome number reduction from an inferred 9-chromosome ancestral state: (a) chromosome breakage followed by loss of centromere DNA and (b) centromere inactivation accompanied by changes in DNA sequence following chromosome-chromosome fusion. We propose AT-rich centromeres drive karyotype diversity in the Malassezia species complex through breakage and inactivation.

We previously observed an unexpected fivefold (35 vs. 200 days) difference in the survival of respiratory chain complex III (CIII) deficient Bcs1l p.S78G mice between two congenic backgrounds. Here, we identify a spontaneous homoplasmic mtDNA variant (m.G14904A, mt-Cybp.D254N), affecting the CIII subunit cytochrome b (MT-CYB), in the background with short survival. We utilize maternal inheritance of mtDNA to confirm this as the causative variant and show that it further decreases the low CIII activity in Bcs1lp.S78G tissues to below survival threshold by 35 days of age. Molecular dynamics simulations predict D254N to restrict the flexibility of MT-CYB ef loop, potentially affecting RISP dynamics. In Rhodobacter cytochrome bc1 complex the equivalent substitution causes a kinetics defect with longer occupancy of RISP head domain towards the quinol oxidation site. These findings represent a unique case of spontaneous mitonuclear epistasis and highlight the role of mtDNA variation as modifier of mitochondrial disease phenotypes.

To better understand the health benefits of lifelong exercise in humans, we conduct global skeletal muscle transcriptomic analyses of long-term endurance- (9 men, 9 women) and strength-trained (7 men) humans compared with age-matched untrained controls (7 men, 8 women). Transcriptomic analysis, Gene Ontology, and genome-scale metabolic modeling demonstrate changes in pathways related to the prevention of metabolic diseases, particularly with endurance training. Our data also show prominent sex differences between controls and that these differences are reduced with endurance training. Additionally, we compare our data with studies examining muscle gene expression before and after a months-long training period in individuals with metabolic diseases. This analysis reveals that training shifts gene expression in individuals with impaired metabolism to become more similar to our endurance-trained group. Overall, our data provide an extensive examination of the accumulated transcriptional changes that occur with decades-long training and identify important "exercise-responsive" genes that could attenuate metabolic disease.

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

Cardiac progenitor formation is one of the earliest committed steps of human cardiogenesis and requires the cooperation of multiple gene sets governed by developmental signaling cascades. To determine the key regulators for cardiac progenitor formation, we have developed a two-stage genome-wide CRISPR-knockout screen. We mimicked the progenitor formation process by differentiating human pluripotent stem cells (hPSCs) into cardiomyocytes, monitored by two distinct stage markers of early cardiac mesodermal formation and commitment to a multipotent heart progenitor cell fate: MESP1 and ISL1, respectively. From the screen output, we compiled a list of 15 candidate genes. After validating seven of them, we identified ZIC2 as an essential gene for cardiac progenitor formation. ZIC2 is known as a master regulator of neurogenesis. hPSCs with ZIC2 mutated still express pluripotency markers. However, their ability to differentiate into cardiomyocytes was greatly attenuated. RNA-Seq profiling of the ZIC2-mutant cells revealed that the mutants switched their cell fate alternatively to the noncardiac cell lineage. Further, single cell RNA-seq analysis showed the ZIC2 mutants affected the apelin receptor-related signaling pathway during mesoderm formation. Our results provide a new link between ZIC2 and human cardiogenesis and document the potential power of a genome-wide unbiased CRISPR-knockout screen to identify the key steps in human mesoderm precursor cell- and heart progenitor cell-fate determination during in vitro hPSC cardiogenesis.

Biodegradation by microorganisms is a useful tool that helps alleviating hydrocarbon pollution in nature. Microbes are more efficient in degradation under aerobic than anaerobic conditions, but the majority of sediment by volume is generally anoxic. Incubation experiments were conducted to study the biodegradation potential of naphthalene-a common polycyclic aromatic hydrocarbon (PAH)-and the diversity of microbial communities in presence/absence of activated carbon (AC) under aerobic/anaerobic conditions. Radio-respirometry experiments with endogenous microorganisms indicated that degradation of naphthalene was strongly stimulated (96%) by the AC addition under anaerobic conditions. In aerobic conditions, however, AC had no effects on naphthalene biodegradation. Bioaugmentation tests with cultured microbial populations grown on naphthalene showed that AC further stimulated (92%) naphthalene degradation in anoxia. Analysis of the 16S rRNA gene sequences implied that sediment amendment with AC increased microbial community diversity and changed community structure. Moreover, the relative abundance of Geobacter, Thiobacillus, Sulfuricurvum, and methanogenic archaea increased sharply after amendment with AC under anaerobic conditions. These results may be explained by the fact that AC particles promoted direct interspecies electron transfer (DIET) between microorganisms involved in PAH degradation pathways. We suggest that important ecosystem functions mediated by microbes-such as hydrocarbon degradation-can be induced and that AC enrichment strategies can be exploited for facilitating bioremediation of anoxic oil-contaminated sediments and soils.

Single-cell RNA sequencing (scRNA-seq) is the leading technique for characterizing the transcriptomes of individual cells in a sample. The latest protocols are scalable to thousands of cells and are being used to compile cell atlases of tissues, organs and organisms. However, the protocols differ substantially with respect to their RNA capture efficiency, bias, scale and costs, and their relative advantages for different applications are unclear. In the present study, we generated benchmark datasets to systematically evaluate protocols in terms of their power to comprehensively describe cell types and states. We performed a multicenter study comparing 13 commonly used scRNA-seq and single-nucleus RNA-seq protocols applied to a heterogeneous reference sample resource. Comparative analysis revealed marked differences in protocol performance. The protocols differed in library complexity and their ability to detect cell-type markers, impacting their predictive value and suitability for integration into reference cell atlases. These results provide guidance both for individual researchers and for consortium projects such as the Human Cell Atlas.

The bowerbirds in New Guinea and Australia include species that build the largest and perhaps most elaborately decorated constructions outside of humans. The males use these courtship bowers, along with their displays, to attract females. In these species, the mating system is polygynous and the females alone incubate and feed the nestlings. The bowerbirds also include ten species of the socially monogamous catbirds in which the male participates in most aspects of raising the young. How the bower-building behavior evolved has remained poorly understood, as no comprehensive phylogeny exists for the family. It has been assumed that the monogamous catbird clade is sister to all polygynous species. We here test this hypothesis using a newly developed pipeline for obtaining homologous alignments of thousands of exonic and intronic regions from genomic data to build a phylogeny. Our well-supported species tree shows that the polygynous, bower-building species are not monophyletic. The result suggests either that bower-building behavior is an ancestral condition in the family that was secondarily lost in the catbirds, or that it has arisen in parallel in two lineages of bowerbirds. We favor the latter hypothesis based on an ancestral character reconstruction showing that polygyny but not bower-building is ancestral in bowerbirds, and on the observation that Scenopoeetes dentirostris, the sister species to one of the bower-building clades, does not build a proper bower but constructs a court for male display. This species is also sexually monomorphic in plumage despite having a polygynous mating system. We argue that the relatively stable tropical and subtropical forest environment in combination with low predator pressure and rich food access (mostly fruit) facilitated the evolution of these unique life-history traits.

Endomyocardial biopsy is a valuable tool in cardiac diagnostics but is limited by low diagnostic yield and significant complication risks. Meanwhile, recent developments in transcriptomic and proteomic technologies promise a wealth of biological data from minimal tissue samples. To take advantage of the minimal tissue amount needed for molecular analyses, we have developed a sub-millimeter endovascular biopsy device, considerably smaller than current clinical equipment, and devised a low-input RNA-sequencing protocol for analyzing small tissue samples. In in vivo evaluation in swine, 81% of biopsy attempts (n = 157) were successful. High quality RNA-sequencing data was generated from 91% of the sequenced cardiac micro-biopsy samples (n = 32). Gene expression signatures of samples taken with the novel device were comparable with a conventional device. No major complications were detected either during procedures or during 7 days' follow-up, despite acquiring a relatively large number of biopsies (median 30) in each animal. In conclusion, the novel device coupled with RNA-sequencing provides a feasible method to obtain molecular data from the myocardium. The method is less traumatic and has a higher flexibility compared to conventional methods, enabling safer and more targeted sampling from different parts of the myocardium.

Litter decomposing Agaricales play key role in terrestrial carbon cycling, but little is known about their decomposition mechanisms. We assembled datasets of 42 gene families involved in plant-cell-wall decomposition from seven newly sequenced litter decomposers and 35 other Agaricomycotina members, mostly white-rot and brown-rot species. Using sequence similarity and phylogenetics, we split the families into phylogroups and compared their gene composition across nutritional strategies. Subsequently, we used Raman spectroscopy to examine the ability of litter decomposers, white-rot fungi, and brown-rot fungi to decompose crystalline cellulose. Both litter decomposers and white-rot fungi share the enzymatic cellulose decomposition, whereas brown-rot fungi possess a distinct mechanism that disrupts cellulose crystallinity. However, litter decomposers and white-rot fungi differ with respect to hemicellulose and lignin degradation phylogroups, suggesting adaptation of the former group to the litter environment. Litter decomposers show high phylogroup diversity, which is indicative of high functional versatility within the group, whereas a set of white-rot species shows adaptation to bulk-wood decomposition. In both groups, we detected species that have unique characteristics associated with hitherto unknown adaptations to diverse wood and litter substrates. Our results suggest that the terms white-rot fungi and litter decomposers mask a much larger functional diversity.

Large-scale sequencing of RNA from individual cells can reveal patterns of gene, isoform and allelic expression across cell types and states 1. However, current short-read single-cell RNA-sequencing methods have limited ability to count RNAs at allele and isoform resolution, and long-read sequencing techniques lack the depth required for large-scale applications across cells2,3. Here we introduce Smart-seq3, which combines full-length transcriptome coverage with a 5' unique molecular identifier RNA counting strategy that enables in silico reconstruction of thousands of RNA molecules per cell. Of the counted and reconstructed molecules, 60% could be directly assigned to allelic origin and 30-50% to specific isoforms, and we identified substantial differences in isoform usage in different mouse strains and human cell types. Smart-seq3 greatly increased sensitivity compared to Smart-seq2, typically detecting thousands more transcripts per cell. We expect that Smart-seq3 will enable large-scale characterization of cell types and states across tissues and organisms.