Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation.

Long-read sequencing technologies substantially overcome the limitations of short-reads but have not been considered as a feasible replacement for population-scale projects, being a combination of too expensive, not scalable enough or too error-prone. Here we develop an efficient and scalable wet lab and computational protocol, Napu, for Oxford Nanopore Technologies long-read sequencing that seeks to address those limitations. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the National Institutes of Health Center for Alzheimer's and Related Dementias. Using a single PromethION flow cell, we can detect single nucleotide polymorphisms with F1-score comparable to Illumina short-read sequencing. Small indel calling remains difficult within homopolymers and tandem repeats, but achieves good concordance to Illumina indel calls elsewhere. Further, we can discover structural variants with F1-score on par with state-of-the-art de novo assembly methods. Our protocol phases small and structural variants at megabase scales and produces highly accurate, haplotype-specific methylation calls.

Diverse Molecular Mechanisms Contribute to Differential Expression of Human Duplicated Genes.

Emerging evidence links genes within human-specific segmental duplications (HSDs) to traits and diseases unique to our species. Strikingly, despite being nearly identical by sequence (>98.5%), paralogous HSD genes are differentially expressed across human cell and tissue types, though the underlying mechanisms have not been examined. We compared cross-tissue mRNA levels of 75 HSD genes from 30 families between humans and chimpanzees and found expression patterns consistent with relaxed selection on or neofunctionalization of derived paralogs. In general, ancestral paralogs exhibited greatest expression conservation with chimpanzee orthologs, though exceptions suggest certain derived paralogs may retain or supplant ancestral functions. Concordantly, analysis of long-read isoform sequencing data sets from diverse human tissues and cell lines found that about half of derived paralogs exhibited globally lower expression. To understand mechanisms underlying these differences, we leveraged data from human lymphoblastoid cell lines (LCLs) and found no relationship between paralogous expression divergence and post-transcriptional regulation, sequence divergence, or copy-number variation. Considering cis-regulation, we reanalyzed ENCODE data and recovered hundreds of previously unidentified candidate CREs in HSDs. We also generated large-insert ChIP-sequencing data for active chromatin features in an LCL to better distinguish paralogous regions. Some duplicated CREs were sufficient to drive differential reporter activity, suggesting they may contribute to divergent cis-regulation of paralogous genes. This work provides evidence that cis-regulatory divergence contributes to novel expression patterns of recent gene duplicates in humans.

Complete de novo assembly of Wolbachia endosymbiont of Drosophila willistoni using long-read genome sequencing.

Wolbachia is an obligate intracellular 𝛼-proteobacterium which commonly infects arthropods and filarial nematodes. Different strains of Wolbachia are capable of a wide range of regulatory manipulations in many hosts and modulate host cellular differentiation to influence host reproduction. The genetic basis for the majority of these phenotypes is unknown. The w Wil strain from the neotropical fruit fly, Drosophila willistoni , exhibits a remarkably high affinity for host germline-derived cells relative to the soma. This trait could be leveraged for understanding how Wolbachia influences the host germline and for controlling host populations in the field. To further the use of this strain in biological and biomedical research, we sequenced the genome of the w Wil strain isolated from host cell culture cells. Here, we present the first high quality nanopore assembly of w Wil, the Wolbachia endosymbiont of D. willistoni . Our assembly resulted in a circular genome of 1.27 Mb with a BUSCO completeness score of 99.7%. Consistent with other insect-associated Wolbachia strains, comparative genomic analysis revealed that wWil has a highly mosaic genome relative to the closely related wMel strain from Drosophila melanogaster .

Persistence of Matrilocal Postmarital Residence Across Multiple Generations in Southern Africa.

Factors such as subsistence turnover, warfare, or interaction between different groups can be major sources of cultural change in human populations. Global demographic shifts such as the transition to agriculture during the Neolithic and more recently the urbanization and globalization of the twentieth century have been major catalysts for cultural change. Here, we test whether cultural traits such as patri/matrilocality and postmarital migration persist in the face of social upheaval and gene flow during the past 150 years in postcolonial South Africa. The recent history of South Africa has seen major demographic shifts that resulted in the displacement and forced sedentism of indigenous Khoekhoe and San populations. During the expansion of the colonial frontier, the Khoe-San admixed with European colonists and enslaved individuals from West/Central Africa, Indonesia, and South Asia, introducing novel cultural norms. We conducted demographic interviews among Nama and Cederberg communities representing nearly 3,000 individuals across three generations. Despite the history of colonial expansion, and the subsequent incorporation of Khoe-San and Khoe-San-descendant communities into a colonial society with strong patrilocal norms, patrilocality is the least common postmarital residence pattern in our study populations today. Our results suggest that more recent forces of integration into the market economy are likely the primary drivers of change in the cultural traits examined in our study. Birthplace had a strong effect on an individual's odds of migration, distance moved, and postmarital residence form. These effects are at least partially explained by the population size of the birthplace. Our results suggest that market factors local to birthplaces are important drivers of residence decisions, although the frequency of matrilocal residence and a geographic and temporal cline in migration and residence patterns also indicate the persistence of some historic Khoe-San cultural traits in contemporary groups.

Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation.

Long-read sequencing technologies substantially overcome the limitations of short-reads but to date have not been considered as feasible replacement at scale due to a combination of being too expensive, not scalable enough, or too error-prone. Here, we develop an efficient and scalable wet lab and computational protocol for Oxford Nanopore Technologies (ONT) long-read sequencing that seeks to provide a genuine alternative to short-reads for large-scale genomics projects. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the NIH Center for Alzheimer's and Related Dementias (CARD). Using a single PromethION flow cell, we can detect SNPs with F1-score better than Illumina short-read sequencing. Small indel calling remains to be difficult inside homopolymers and tandem repeats, but is comparable to Illumina calls elsewhere. Further, we can discover structural variants with F1-score comparable to state-of the-art methods involving Pacific Biosciences HiFi sequencing and trio information (but at a lower cost and greater throughput). Using ONT based phasing, we can then combine and phase small and structural variants at megabase scales. Our protocol also produces highly accurate, haplotype-specific methylation calls. Overall, this makes large-scale long-read sequencing projects feasible; the protocol is currently being used to sequence thousands of brain-based genomes as a part of the NIH CARD initiative. We provide the protocol and software as open-source integrated pipelines for generating phased variant calls and assemblies.

The recombination landscape of the Khoe-San likely represents the upper limits of recombination divergence in humans.

BACKGROUND: Recombination maps are important resources for epidemiological and evolutionary analyses; however, there are currently no recombination maps representing any African population outside of those with West African ancestry. We infer the demographic history for the Nama, an indigenous Khoe-San population of southern Africa, and derive a novel, population-specific recombination map from the whole genome sequencing of 54 Nama individuals. We hypothesise that there are no publicly available recombination maps representative of the Nama, considering the deep population divergence and subsequent isolation of the Khoe-San from other African groups. RESULTS: We show that the recombination landscape of the Nama does not cluster with any continental groups with publicly available representative recombination maps. Finally, we use selection scans as an example of how fine-scale differences between the Nama recombination map and the combined Phase II HapMap recombination map can impact the outcome of selection scans. CONCLUSIONS: Fine-scale differences in recombination can meaningfully alter the results of a selection scan. The recombination map we infer likely represents an upper bound on the extent of divergence we expect to see for a recombination map in humans and would be of interest to any researcher that wants to test the sensitivity of population genetic or GWAS analysis to recombination map input.

Effects of pairing on color change and central gene expression in lined seahorses.

Social monogamy is a reproductive strategy characterized by pair living and defense of a common territory. Pair bonding, sometimes displayed by monogamous species, is an affective construct that includes preference for a specific partner, distress upon separation, and the ability of the partner to buffer against stress. Many seahorse species show a monogamous social structure in the wild, but their pair bond has not been well studied. We examined the gene expression of lined seahorses (Hippocampus erectus) during and after the process of pairing in the laboratory as well as color change (luminance), a potential form of social communication and behavioral synchrony between pair mates. When a seahorse of either sex was interacting with its pair mate, their changes in luminance ("brightness") were correlated and larger than when interacting with an opposite-sex stranger. At the conclusion of testing, subjects were euthanized, RNA was extracted from whole brains and analyzed via RNA sequencing. Changes in gene expression in paired males versus those that were unpaired included processes governing metabolic activity, hormones and cilia. Perhaps most interesting is the overlap in gene expression change induced by pairing in both male seahorses and male prairie voles, including components of hormone systems regulating reproduction. Because of our limited sample size, we consider our results and interpretations to be preliminary, and prompts for further exploration. Future studies will expand upon these findings and investigate the neuroendocrine and genetic basis of these behaviors.

Hunter-gatherer genomes reveal diverse demographic trajectories during the rise of farming in Eastern Africa.

The fate of hunting and gathering populations following the rise of agriculture and pastoralism remains a topic of debate in the study of human prehistory. Studies of ancient and modern genomes have found that autochthonous groups were largely replaced by expanding farmer populations with varying levels of gene flow, a characterization that is influenced by the almost universal focus on the European Neolithic.1-5 We sought to understand the demographic impact of an ongoing cultural transition to farming in Southwest Ethiopia, one of the last regions in Africa to experience such shifts.6 Importantly, Southwest Ethiopia is home to several of the world's remaining hunter-gatherer groups, including the Chabu people, who are currently transitioning away from their traditional mode of subsistence.7 We generated genome-wide data from the Chabu and four neighboring populations, the Majang, Shekkacho, Bench, and Sheko, to characterize their genetic ancestry and estimate their effective population sizes over the last 60 generations. We show that the Chabu are a distinct population closely related to ancient people who occupied Southwest Ethiopia >4,500 years ago. Furthermore, the Chabu are undergoing a severe population bottleneck, which began approximately 1,400 years ago. By analyzing eleven Eastern African populations, we find evidence for divergent demographic trajectories among hunter-gatherer-descendant groups. Our results illustrate that although foragers respond to encroaching agriculture and pastoralism with multiple strategies, including cultural adoption of agropastoralism, gene flow, and economic specialization, they often face population decline.

Identification of Structural Variation in Chimpanzees Using Optical Mapping and Nanopore Sequencing.

Recent efforts to comprehensively characterize great ape genetic diversity using short-read sequencing and single-nucleotide variants have led to important discoveries related to selection within species, demographic history, and lineage-specific traits. Structural variants (SVs), including deletions and inversions, comprise a larger proportion of genetic differences between and within species, making them an important yet understudied source of trait divergence. Here, we used a combination of long-read and -range sequencing approaches to characterize the structural variant landscape of two additional Pan troglodytes verus individuals, one of whom carries 13% admixture from Pan troglodytes troglodytes. We performed optical mapping of both individuals followed by nanopore sequencing of one individual. Filtering for larger variants (>10 kbp) and combined with genotyping of SVs using short-read data from the Great Ape Genome Project, we identified 425 deletions and 59 inversions, of which 88 and 36, respectively, were novel. Compared with gene expression in humans, we found a significant enrichment of chimpanzee genes with differential expression in lymphoblastoid cell lines and induced pluripotent stem cells, both within deletions and near inversion breakpoints. We examined chromatin-conformation maps from human and chimpanzee using these same cell types and observed alterations in genomic interactions at SV breakpoints. Finally, we focused on 56 genes impacted by SVs in >90% of chimpanzees and absent in humans and gorillas, which may contribute to chimpanzee-specific features. Sequencing a greater set of individuals from diverse subspecies will be critical to establish the complete landscape of genetic variation in chimpanzees.