Longitudinal Multi-omics Analyses Identify Responses of Megakaryocytes, Erythroid Cells, and Plasmablasts as Hallmarks of Severe COVID-19.

Temporal resolution of cellular features associated with a severe COVID-19 disease trajectory is needed for understanding skewed immune responses and defining predictors of outcome. Here, we performed a longitudinal multi-omics study using a two-center cohort of 14 patients. We analyzed the bulk transcriptome, bulk DNA methylome, and single-cell transcriptome (>358,000 cells, including BCR profiles) of peripheral blood samples harvested from up to 5 time points. Validation was performed in two independent cohorts of COVID-19 patients. Severe COVID-19 was characterized by an increase of proliferating, metabolically hyperactive plasmablasts. Coinciding with critical illness, we also identified an expansion of interferon-activated circulating megakaryocytes and increased erythropoiesis with features of hypoxic signaling. Megakaryocyte- and erythroid-cell-derived co-expression modules were predictive of fatal disease outcome. The study demonstrates broad cellular effects of SARS-CoV-2 infection beyond adaptive immune cells and provides an entry point toward developing biomarkers and targeted treatments of patients with COVID-19.

Occasional paternal inheritance of the germline-restricted chromosome in songbirds.

Songbirds have one special accessory chromosome, the so-called germline-restricted chromosome (GRC), which is only present in germline cells and absent from all somatic tissues. Earlier work on the zebra finch (Taeniopygia guttata castanotis) showed that the GRC is inherited only through the female line-like the mitochondria-and is eliminated from the sperm during spermatogenesis. Here, we show that the GRC has the potential to be paternally inherited. Confocal microscopy using GRC-specific fluorescent in situ hybridization probes indicated that a considerable fraction of sperm heads (1 to 19%) in zebra finch ejaculates still contained the GRC. In line with these cytogenetic data, sequencing of ejaculates revealed that individual males from two families differed strongly and consistently in the number of GRCs in their ejaculates. Examining a captive-bred male hybrid of the two zebra finch subspecies (T. g. guttata and T. g. castanotis) revealed that the mitochondria originated from a castanotis mother, whereas the GRC came from a guttata father. Moreover, analyzing GRC haplotypes across nine castanotis matrilines, estimated to have diverged for up to 250,000 y, showed surprisingly little variability among GRCs. This suggests that a single GRC haplotype has spread relatively recently across all examined matrilines. A few diagnostic GRC mutations that arose since this inferred spreading suggest that the GRC has continued to jump across matriline boundaries. Our findings raise the possibility that certain GRC haplotypes could selfishly spread through the population via occasional paternal transmission, thereby outcompeting other GRC haplotypes that were limited to strict maternal inheritance, even if this was partly detrimental to organismal fitness.

An integrated molecular risk score early in life for subsequent childhood asthma risk.

BACKGROUND: Numerous children present with early wheeze symptoms, yet solely a subgroup develops childhood asthma. Early identification of children at risk is key for clinical monitoring, timely patient-tailored treatment, and preventing chronic, severe sequelae. For early prediction of childhood asthma, we aimed to define an integrated risk score combining established risk factors with genome-wide molecular markers at birth, complemented by subsequent clinical symptoms/diagnoses (wheezing, atopic dermatitis, food allergy). METHODS: Three longitudinal birth cohorts (PAULINA/PAULCHEN, n = 190 + 93 = 283, PASTURE, n = 1133) were used to predict childhood asthma (age 5-11) including epidemiological characteristics and molecular markers: genotype, DNA methylation and mRNA expression (RNASeq/NanoString). Apparent (ap) and optimism-corrected (oc) performance (AUC/R2) was assessed leveraging evidence from independent studies (Naïve-Bayes approach) combined with high-dimensional logistic regression models (LASSO). RESULTS: Asthma prediction with epidemiological characteristics at birth (maternal asthma, sex, farm environment) yielded an ocAUC = 0.65. Inclusion of molecular markers as predictors resulted in an improvement in apparent prediction performance, however, for optimism-corrected performance only a moderate increase was observed (upto ocAUC = 0.68). The greatest discriminate power was reached by adding the first symptoms/diagnosis (up to ocAUC = 0.76; increase of 0.08, p = .002). Longitudinal analysis of selected mRNA expression in PASTURE (cord blood, 1, 4.5, 6 years) showed that expression at age six had the strongest association with asthma and correlation of genes getting larger over time (r = .59, p < .001, 4.5-6 years). CONCLUSION: Applying epidemiological predictors alone showed moderate predictive abilities. Molecular markers from birth modestly improved prediction. Allergic symptoms/diagnoses enhanced the power of prediction, which is important for clinical practice and for the design of future studies with molecular markers.

Dynamic interactions within the host-associated microbiota cause tumor formation in the basal metazoan Hydra.

The extent to which disturbances in the resident microbiota can compromise an animal's health is poorly understood. Hydra is one of the evolutionary oldest animals with naturally occurring tumors. Here, we found a causal relationship between an environmental spirochete (Turneriella spec.) and tumorigenesis in Hydra. Unexpectedly, virulence of this pathogen requires the presence of Pseudomonas spec., a member of Hydra´s beneficial microbiome indicating that dynamic interactions between a resident bacterium and a pathogen cause tumor formation. The observation points to the crucial role of commensal bacteria in maintaining tissue homeostasis and adds support to the view that microbial community interactions are essential for disease. These findings in an organism that shares deep evolutionary connections with all animals have implications for our understanding of cancer.

MRI-Based Iron Phenotyping and Patient Selection for Next-Generation Sequencing of Non-Homeostatic Iron Regulator Hemochromatosis Genes.

BACKGROUND AND AIMS: High serum ferritin is frequent among patients with chronic liver disease and commonly associated with hepatic iron overload. Genetic causes of high liver iron include homozygosity for the p.Cys282Tyr variant in homeostatic iron regulator (HFE) and rare variants in non-HFE genes. The aims of the present study were to describe the landscape and frequency of mutations in hemochromatosis genes and determine whether patient selection by noninvasive hepatic iron quantification using MRI improves the diagnostic yield of next-generation sequencing (NGS) in patients with hyperferritinemia. APPROACH AND RESULTS: A cohort of 410 unselected liver clinic patients with high serum ferritin (defined as ≥200 µg/L for women and ≥300 µg/L for men) was investigated by HFE genotyping and abdominal MRI R2*. Forty-one (10%) patients were homozygous for the p.Cys282Tyr variant in HFE. Of the remaining 369 patients, 256 (69%) had high transferrin saturation (TSAT; ≥45%) and 199 (53%) had confirmed hepatic iron overload (liver R2* ≥70 s-1 ). NGS of hemochromatosis genes was carried out in 180 patients with hepatic iron overload, and likely pathogenic variants were identified in 68 of 180 (38%) patients, mainly in HFE (79%), ceruloplasmin (25%), and transferrin receptor 2 (19%). Low spleen iron (R2* <50 s-1 ), but not TSAT, was significantly associated with the presence of mutations. In 167 patients (93%), no monogenic cause of hepatic iron overload could be identified. CONCLUSIONS: In patients without homozygosity for p.Cys282Tyr, coincident pathogenic variants in HFE and non-HFE genes could explain hyperferritinemia with hepatic iron overload in a subset of patients. Unlike HFE hemochromatosis, this type of polygenic hepatic iron overload presents with variable TSAT. High ferritin in blood is an indicator of the iron storage disease, hemochromatosis. A simple genetic test establishes this diagnosis in the majority of patients affected. MRI of the abdomen can guide further genetic testing.

Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen.

The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).

Investigation of the host transcriptional response to intracellular bacterial infection using Dictyostelium discoideum as a host model.

BACKGROUND: During infection by intracellular pathogens, a highly complex interplay occurs between the infected cell trying to degrade the invader and the pathogen which actively manipulates the host cell to enable survival and proliferation. Many intracellular pathogens pose important threats to human health and major efforts have been undertaken to better understand the host-pathogen interactions that eventually determine the outcome of the infection. Over the last decades, the unicellular eukaryote Dictyostelium discoideum has become an established infection model, serving as a surrogate macrophage that can be infected with a wide range of intracellular pathogens. In this study, we use high-throughput RNA-sequencing to analyze the transcriptional response of D. discoideum when infected with Mycobacterium marinum and Legionella pneumophila. The results were compared to available data from human macrophages. RESULTS: The majority of the transcriptional regulation triggered by the two pathogens was found to be unique for each bacterial challenge. Hallmark transcriptional signatures were identified for each infection, e.g. induction of endosomal sorting complexes required for transport (ESCRT) and autophagy genes in response to M. marinum and inhibition of genes associated with the translation machinery and energy metabolism in response to L. pneumophila. However, a common response to the pathogenic bacteria was also identified, which was not induced by non-pathogenic food bacteria. Finally, comparison with available data sets of regulation in human monocyte derived macrophages shows that the elicited response in D. discoideum is in many aspects similar to what has been observed in human immune cells in response to Mycobacterium tuberculosis and L. pneumophila. CONCLUSIONS: Our study presents high-throughput characterization of D. discoideum transcriptional response to intracellular pathogens using RNA-seq. We demonstrate that the transcriptional response is in essence distinct to each pathogen and that in many cases, the corresponding regulation is recapitulated in human macrophages after infection by mycobacteria and L. pneumophila. This indicates that host-pathogen interactions are evolutionary conserved, derived from the early interactions between free-living phagocytic cells and bacteria. Taken together, our results strengthen the use of D. discoideum as a general infection model.

An Evaluation of Function of Multicopy Noncoding RNAs in Mammals Using ENCODE/FANTOM Data and Comparative Genomics.

Mammalian diversification has coincided with a rapid proliferation of various types of noncoding RNAs, including members of both snRNAs and snoRNAs. The significance of this expansion however remains obscure. While some ncRNA copy-number expansions have been linked to functionally tractable effects, such events may equally likely be neutral, perhaps as a result of random retrotransposition. Hindering progress in our understanding of such observations is the difficulty in establishing function for the diverse features that have been identified in our own genome. Projects such as ENCODE and FANTOM have revealed a hidden world of genomic expression patterns, as well as a host of other potential indicators of biological function. However, such projects have been criticized, particularly from practitioners in the field of molecular evolution, where many suspect these data provide limited insight into biological function. The molecular evolution community has largely taken a skeptical view, thus it is important to establish tests of function. We use a range of data, including data drawn from ENCODE and FANTOM, to examine the case for function for the recent copy number expansion in mammals of six evolutionarily ancient RNA families involved in splicing and rRNA maturation. We use several criteria to assess evidence for function: conservation of sequence and structure, genomic synteny, evidence for transposition, and evidence for species-specific expression. Applying these criteria, we find that only a minority of loci show strong evidence for function and that, for the majority, we cannot reject the null hypothesis of no function.

Rates of Mutation and Recombination in Siphoviridae Phage Genome Evolution over Three Decades.

The evolution of asexual organisms is driven not only by the inheritance of genetic modification but also by the acquisition of foreign DNA. The contribution of vertical and horizontal processes to genome evolution depends on their rates per year and is quantified by the ratio of recombination to mutation. These rates have been estimated for bacteria; however, no estimates have been reported for phages. Here, we delineate the contribution of mutation and recombination to dsDNA phage genome evolution. We analyzed 34 isolates of the 936 group of Siphoviridae phages using a Lactococcus lactis strain from a single dairy over 29 years. We estimate a constant substitution rate of 1.9 × 10-4 substitutions per site per year due to mutation that is within the range of estimates for eukaryotic RNA and DNA viruses. The reconstruction of recombination events reveals a constant rate of five recombination events per year and 4.5 × 10-3 nucleotide alterations due to recombination per site per year. Thus, the recombination rate exceeds the substitution rate, resulting in a relative effect of recombination to mutation (r/m) of ∼24 that is homogenous over time. Especially in the early transcriptional region, we detect frequent gene loss and regain due to recombination with phages of the 936 group, demonstrating the role of the 936 group pangenome as a reservoir of genetic variation. The observed substitution rate homogeneity conforms to the neutral theory of evolution; hence, the neutral theory can be applied to phage genome evolution and also to genetic variation brought about by recombination.

A comprehensive, cell specific microRNA catalogue of human peripheral blood.

With this study, we provide a comprehensive reference dataset of detailed miRNA expression profiles from seven types of human peripheral blood cells (NK cells, B lymphocytes, cytotoxic T lymphocytes, T helper cells, monocytes, neutrophils and erythrocytes), serum, exosomes and whole blood. The peripheral blood cells from buffy coats were typed and sorted using FACS/MACS. The overall dataset was generated from 450 small RNA libraries using high-throughput sequencing. By employing a comprehensive bioinformatics and statistical analysis, we show that 3' trimming modifications as well as composition of 3' added non-templated nucleotides are distributed in a lineage-specific manner-the closer the hematopoietic progenitors are, the higher their similarities in sequence variation of the 3' end. Furthermore, we define the blood cell-specific miRNA and isomiR expression patterns and identify novel cell type specific miRNA candidates. The study provides the most comprehensive contribution to date towards a complete miRNA catalogue of human peripheral blood, which can be used as a reference for future studies. The dataset has been deposited in GEO and also can be explored interactively following this link: http://134.245.63.235/ikmb-tools/bloodmiRs.

Global characterization of the Dicer-like protein DrnB roles in miRNA biogenesis in the social amoeba Dictyostelium discoideum.

Micro (mi)RNAs regulate gene expression in many eukaryotic organisms where they control diverse biological processes. Their biogenesis, from primary transcripts to mature miRNAs, have been extensively characterized in animals and plants, showing distinct differences between these phylogenetically distant groups of organisms. However, comparably little is known about miRNA biogenesis in organisms whose evolutionary position is placed in between plants and animals and/or in unicellular organisms. Here, we investigate miRNA maturation in the unicellular amoeba Dictyostelium discoideum, belonging to Amoebozoa, which branched out after plants but before animals. High-throughput sequencing of small RNAs and poly(A)-selected RNAs demonstrated that the Dicer-like protein DrnB is required, and essentially specific, for global miRNA maturation in D. discoideum. Our RNA-seq data also showed that longer miRNA transcripts, generally preceded by a T-rich putative promoter motif, accumulate in a drnB knock-out strain. For two model miRNAs we defined the transcriptional start sites (TSSs) of primary (pri)-miRNAs and showed that they carry the RNA polymerase II specific m7G-cap. The generation of the 3'-ends of these pri-miRNAs differs, with pri-mir-1177 reading into the downstream gene, and pri-mir-1176 displaying a distinct end. This 3´-end is processed to shorter intermediates, stabilized in DrnB-depleted cells, of which some carry a short oligo(A)-tail. Furthermore, we identified 10 new miRNAs, all DrnB dependent and developmentally regulated. Thus, the miRNA machinery in D. discoideum shares features with both plants and animals, which is in agreement with its evolutionary position and perhaps also an adaptation to its complex lifestyle: unicellular growth and multicellular development.

Opportunities and challenges of whole-genome and -exome sequencing.

Recent advances in the development of sequencing technologies provide researchers with unprecedented possibilities for genetic analyses. In this review, we will discuss the history of genetic studies and the progress driven by next-generation sequencing (NGS), using complex inflammatory bowel diseases as an example. We focus on the opportunities, but also challenges that researchers are facing when working with NGS data to unravel the genetic causes underlying diseases.

Does the Ribosome Challenge our Understanding of the RNA World?

In a recent article published in these pages, Bowman and colleagues propose that the ribosome represents a challenge to the RNA world model, a long-standing framework to explain the origin of DNA and genetically encoded proteins from a hypothetical RNA-based system. Specifically, they outline a scenario for the emergence and subsequent coevolution of the peptidyl transferase centre (PTC) of the ribosome with non-templated peptide products of this RNA through chemical evolution. They also propose that the PTC would have predated the emergence of enzymatic RNA replication, and that this in turn indicates that the RNA world never existed. We and others have previously incorporated non-templated peptide production as an early stage in the evolution of protein synthesis, which we would count as a chemical process, in agreement with Bowman and colleagues' model. However, their model raises an important question: to what extent could early protein synthesis and its products have evolved in the absence of Darwinian processes? We argue that evolution of the early ribosome requires Darwinian evolution, and that, while chemical evolution could give rise to peptidyl transferase activity, it is insufficient for subsequent improvement of a proto-PTC, or for ongoing coevolution of the proto-PTC with its early non-templated peptide products. We conclude that it is difficult to preclude the involvement of replicative processes, themselves subject to Darwinian evolution, from the evolution of the PTC. Finally, Bowman et al. call into question current models for the RNA to protein transition. We show that the difficulty that Bowman et al. have with this scenario is down to a misreading of our previous work.

Linked genetic variants on chromosome 10 control ear morphology and body mass among dog breeds.

BACKGROUND: The domestic dog is a rich resource for mapping the genetic components of phenotypic variation due to its unique population history involving strong artificial selection. Genome-wide association studies have revealed a number of chromosomal regions where genetic variation associates with morphological characters that typify dog breeds. A region on chromosome 10 is among those with the highest levels of genetic differentiation between dog breeds and is associated with body mass and ear morphology, a common motif of animal domestication. We characterised variation in this region to uncover haplotype structure and identify candidate functional variants. RESULTS: We first identified SNPs that strongly associate with body mass and ear type by comparing sequence variation in a 3 Mb region between 19 breeds with a variety of phenotypes. We next genotyped a subset of 123 candidate SNPs in 288 samples from 46 breeds to identify the variants most highly associated with phenotype and infer haplotype structure. A cluster of SNPs that associate strongly with the drop ear phenotype is located within a narrow interval downstream of the gene MSRB3, which is involved in human hearing. These SNPs are in strong genetic linkage with another set of variants that correlate with body mass within the gene HMGA2, which affects human height. In addition we find evidence that this region has been under selection during dog domestication, and identify a cluster of SNPs within MSRB3 that are highly differentiated between dogs and wolves. CONCLUSIONS: We characterise genetically linked variants that potentially influence ear type and body mass in dog breeds, both key traits that have been modified by selective breeding that may also be important for domestication. The finding that variants on long haplotypes have effects on more than one trait suggests that genetic linkage can be an important determinant of the phenotypic response to selection in domestic animals.

Unsupervised genome-wide recognition of local relationship patterns.

BACKGROUND: Phenomena such as incomplete lineage sorting, horizontal gene transfer, gene duplication and subsequent sub- and neo-functionalisation can result in distinct local phylogenetic relationships that are discordant with species phylogeny. In order to assess the possible biological roles for these subdivisions, they must first be identified and characterised, preferably on a large scale and in an automated fashion. RESULTS: We developed Saguaro, a combination of a Hidden Markov Model (HMM) and a Self Organising Map (SOM), to characterise local phylogenetic relationships among aligned sequences using cacti, matrices of pair-wise distance measures. While the HMM determines the genomic boundaries from aligned sequences, the SOM hypothesises new cacti in an unsupervised and iterative fashion based on the regions that were modelled least well by existing cacti. After testing the software on simulated data, we demonstrate the utility of Saguaro by testing two different data sets: (i) 181 Dengue virus strains, and (ii) 5 primate genomes. Saguaro identifies regions under lineage-specific constraint for the first set, and genomic segments that we attribute to incomplete lineage sorting in the second dataset. Intriguingly for the primate data, Saguaro also classified an additional ~3% of the genome as most incompatible with the expected species phylogeny. A substantial fraction of these regions was found to overlap genes associated with both the innate and adaptive immune systems. CONCLUSIONS: Saguaro detects distinct cacti describing local phylogenetic relationships without requiring any a priori hypotheses. We have successfully demonstrated Saguaro's utility with two contrasting data sets, one containing many members with short sequences (Dengue viral strains: n = 181, genome size = 10,700 nt), and the other with few members but complex genomes (related primate species: n = 5, genome size = 3 Gb), suggesting that the software is applicable to a wide variety of experimental populations. Saguaro is written in C++, runs on the Linux operating system, and can be downloaded from http://saguarogw.sourceforge.net/.

Comparative analysis of RNA families reveals distinct repertoires for each domain of life.

The RNA world hypothesis, that RNA genomes and catalysts preceded DNA genomes and genetically-encoded protein catalysts, has been central to models for the early evolution of life on Earth. A key part of such models is continuity between the earliest stages in the evolution of life and the RNA repertoires of extant lineages. Some assessments seem consistent with a diverse RNA world, yet direct continuity between modern RNAs and an RNA world has not been demonstrated for the majority of RNA families, and, anecdotally, many RNA functions appear restricted in their distribution. Despite much discussion of the possible antiquity of RNA families, no systematic analyses of RNA family distribution have been performed. To chart the broad evolutionary history of known RNA families, we performed comparative genomic analysis of over 3 million RNA annotations spanning 1446 families from the Rfam 10 database. We report that 99% of known RNA families are restricted to a single domain of life, revealing discrete repertoires for each domain. For the 1% of RNA families/clans present in more than one domain, over half show evidence of horizontal gene transfer (HGT), and the rest show a vertical trace, indicating the presence of a complex protein synthesis machinery in the Last Universal Common Ancestor (LUCA) and consistent with the evolutionary history of the most ancient protein-coding genes. However, with limited interdomain transfer and few RNA families exhibiting demonstrable antiquity as predicted under RNA world continuity, our results indicate that the majority of modern cellular RNA repertoires have primarily evolved in a domain-specific manner.

Host and microbiome jointly contribute to environmental adaptation.

Most animals and plants have associated microorganisms, collectively referred to as their microbiomes, which can provide essential functions. Given their importance, host-associated microbiomes have the potential to contribute substantially to adaptation of the host-microbiome assemblage (the "metaorganism"). Microbiomes may be especially important for rapid adaptation to novel environments because microbiomes can change more rapidly than host genomes. However, it is not well understood how hosts and microbiomes jointly contribute to metaorganism adaptation. We developed a model system with which to disentangle the contributions of hosts and microbiomes to metaorganism adaptation. We established replicate mesocosms containing the nematode Caenorhabditis elegans co-cultured with microorganisms in a novel complex environment (laboratory compost). After approximately 30 nematode generations (100 days), we harvested worm populations and associated microbiomes, and subjected them to a common garden experiment designed to unravel the impacts of microbiome composition and host genetics on metaorganism adaptation. We observed that adaptation took different trajectories in different mesocosm lines, with some increasing in fitness and others decreasing, and that interactions between host and microbiome played an important role in these contrasting evolutionary paths. We chose two exemplary mesocosms (one with a fitness increase and one with a decrease) for detailed study. For each example, we identified specific changes in both microbiome composition (for both bacteria and fungi) and nematode gene expression associated with each change in fitness. Our study provides experimental evidence that adaptation to a novel environment can be jointly influenced by host and microbiome.

Comparative genomics of eukaryotic small nucleolar RNAs reveals deep evolutionary ancestry amidst ongoing intragenomic mobility.

BACKGROUND: Small nucleolar (sno)RNAs are required for posttranscriptional processing and modification of ribosomal, spliceosomal and messenger RNAs. Their presence in both eukaryotes and archaea indicates that snoRNAs are evolutionarily ancient. The location of some snoRNAs within the introns of ribosomal protein genes has been suggested to belie an RNA world origin, with the exons of the earliest protein-coding genes having evolved around snoRNAs after the advent of templated protein synthesis. Alternatively, this intronic location may reflect more recent selection for coexpression of snoRNAs and ribosomal components, ensuring rRNA modification by snoRNAs during ribosome synthesis. To gain insight into the evolutionary origins of this genetic organization, we examined the antiquity of snoRNA families and the stability of their genomic location across 44 eukaryote genomes. RESULTS: We report that dozens of snoRNA families are traceable to the Last Eukaryotic Common Ancestor (LECA), but find only weak similarities between the oldest eukaryotic snoRNAs and archaeal snoRNA-like genes. Moreover, many of these LECA snoRNAs are located within the introns of host genes independently traceable to the LECA. Comparative genomic analyses reveal the intronic location of LECA snoRNAs is not ancestral however, suggesting the pattern we observe is the result of ongoing intragenomic mobility. Analysis of human transcriptome data indicates that the primary requirement for hosting intronic snoRNAs is a broad expression profile. Consistent with ongoing mobility across broadly-expressed genes, we report a case of recent migration of a non-LECA snoRNA from the intron of a ubiquitously expressed non-LECA host gene into the introns of two LECA genes during the evolution of primates. CONCLUSIONS: Our analyses show that snoRNAs were a well-established family of RNAs at the time when eukaryotes began to diversify. While many are intronic, this association is not evolutionarily stable across the eukaryote tree; ongoing intragenomic mobility has erased signal of their ancestral gene organization, and neither introns-first nor evolved co-expression adequately explain our results. We therefore present a third model - constrained drift - whereby individual snoRNAs are intragenomically mobile and may occupy any genomic location from which expression satisfies phenotype.

High-throughput method for the hybridisation-based targeted enrichment of long genomic fragments for PacBio third-generation sequencing.

Hybridisation-based targeted enrichment is a widely used and well-established technique in high-throughput second-generation short-read sequencing. Despite the high potential to genetically resolve highly repetitive and variable genomic sequences by, for example PacBio third-generation sequencing, targeted enrichment for long fragments has not yet established the same high-throughput due to currently existing complex workflows and technological dependencies. We here describe a scalable targeted enrichment protocol for fragment sizes of >7 kb. For demonstration purposes we developed a custom blood group panel of challenging loci. Test results achieved > 65% on-target rate, good coverage (142.7×) and sufficient coverage evenness for both non-paralogous and paralogous targets, and sufficient non-duplicate read counts (83.5%) per sample for a highly multiplexed enrichment pool of 16 samples. We genotyped the blood groups of nine patients employing highly accurate phased assemblies at an allelic resolution that match reference blood group allele calls determined by SNP array and NGS genotyping. Seven Genome-in-a-Bottle reference samples achieved high recall (96%) and precision (99%) rates. Mendelian error rates were 0.04% and 0.13% for the included Ashkenazim and Han Chinese trios, respectively. In summary, we provide a protocol and first example for accurate targeted long-read sequencing that can be used in a high-throughput fashion.