Cortical plasticity is associated with blood-brain barrier modulation.

Brain microvessels possess the unique properties of a blood-brain barrier (BBB), tightly regulating the passage of molecules from the blood to the brain neuropil and vice versa. In models of brain injury, BBB dysfunction and the associated leakage of serum albumin to the neuropil have been shown to induce pathological plasticity, neuronal hyper-excitability, and seizures. The effect of neuronal activity on BBB function and whether it plays a role in plasticity in the healthy brain remain unclear. Here we show that neuronal activity induces modulation of microvascular permeability in the healthy brain and that it has a role in local network reorganization. Combining simultaneous electrophysiological recording and vascular imaging with transcriptomic analysis in rats, and functional and BBB-mapping MRI in human subjects, we show that prolonged stimulation of the limb induces a focal increase in BBB permeability in the corresponding somatosensory cortex that is associated with long-term synaptic plasticity. We further show that the increased microvascular permeability depends on neuronal activity and involves caveolae-mediated transcytosis and transforming growth factor ß signaling. Our results reveal a role of BBB modulation in cortical plasticity in the healthy brain, highlighting the importance of neurovascular interactions for sensory experience and learning.

loco-pipe: an automated pipeline for population genomics with low-coverage whole-genome sequencing.

Summary: We developed loco-pipe, a Snakemake pipeline that seamlessly streamlines a set of essential population genomic analyses for low-coverage whole genome sequencing (lcWGS) data. loco-pipe is highly automated, easily customizable, massively parallelized, and thus is a valuable tool for both new and experienced users of lcWGS. Availability and implementation: loco-pipe is published under the GPLv3. It is freely available on GitHub (github.com/sudmantlab/loco-pipe) and archived on Zenodo (doi.org/10.5281/zenodo.10425920).

Insights into non-crossover recombination from long-read sperm sequencing.

Meiotic recombination is a fundamental process that generates genetic diversity by creating new combinations of existing alleles. Although human crossovers have been studied at the pedigree, population and single-cell level, the more frequent non-crossover events that lead to gene conversion are harder to study, particularly at the individual level. Here we show that single high-fidelity long sequencing reads from sperm can capture both crossovers and non-crossovers, allowing effectively arbitrary sample sizes for analysis from one male. Using fifteen sperm samples from thirteen donors we demonstrate variation between and within donors for the rates of different types of recombination. Intriguingly, we observe a tendency for non-crossover gene conversions to occur upstream of nearby PRDM9 binding sites, whereas crossover locations have a slight downstream bias. We further provide evidence for two distinct non-crossover processes. One gives rise to the vast majority of non-crossovers with mean conversion tract length under 50bp, which we suggest is an outcome of standard PRDM9-induced meiotic recombination. In contrast ~2% of non-crossovers have much longer mean tract length, and potentially originate from the same process as complex events with more than two haplotype switches, which is not associated with PRDM9 binding sites and is also seen in somatic cells.

Structural variation in humans and our primate kin in the era of telomere-to-telomere genomes and pangenomics.

Structural variants (SVs) account for the majority of base pair differences both within and between primate species. However, our understanding of inter- and intra-species SV has been historically hampered by the quality of draft primate genomes and the absence of genome resources for key taxa. Recently, advances in long-read sequencing and genome assembly have begun to radically reshape our understanding of SVs. Two landmark achievements include the publication of a human telomere-to-telomere (T2T) genome as well as the development of the first human pangenome reference. In this review, we first look back to the major works laying the foundation for these projects. We then examine the ways in which T2T genome assemblies and pangenomes are transforming our understanding of and approach to primate SV. Finally, we discuss what the future of primate SV research may look like in the era of T2T genomes and pangenomics.

Hypoxia exposure blunts angiogenic signaling and upregulates the antioxidant system in endothelial cells derived from elephant seals.

BACKGROUND: Elephant seals exhibit extreme hypoxemic tolerance derived from repetitive hypoxia/reoxygenation episodes they experience during diving bouts. Real-time assessment of the molecular changes underlying protection against hypoxic injury in seals remains restricted by their at-sea inaccessibility. Hence, we developed a proliferative arterial endothelial cell culture model from elephant seals and used RNA-seq, functional assays, and confocal microscopy to assess the molecular response to prolonged hypoxia. RESULTS: Seal and human endothelial cells exposed to 1% O2 for up to 6 h respond differently to acute and prolonged hypoxia. Seal cells decouple stabilization of the hypoxia-sensitive transcriptional regulator HIF-1α from angiogenic signaling. Rapid upregulation of genes involved in glutathione (GSH) metabolism supports the maintenance of GSH pools, and intracellular succinate increases in seal but not human cells. High maximal and spare respiratory capacity in seal cells after hypoxia exposure occurs in concert with increasing mitochondrial branch length and independent from major changes in extracellular acidification rate, suggesting that seal cells recover oxidative metabolism without significant glycolytic dependency after hypoxia exposure. CONCLUSIONS: We found that the glutathione antioxidant system is upregulated in seal endothelial cells during hypoxia, while this system remains static in comparable human cells. Furthermore, we found that in contrast to human cells, hypoxia exposure rapidly activates HIF-1 in seal cells, but this response is decoupled from the canonical angiogenesis pathway. These results highlight the unique mechanisms that confer extraordinary tolerance to limited oxygen availability in a champion diving mammal.

Global diversity, recurrent evolution, and recent selection on amylase structural haplotypes in humans.

The adoption of agriculture, first documented ~12,000 years ago in the Fertile Crescent, triggered a rapid shift toward starch-rich diets in human populations. Amylase genes facilitate starch digestion and increased salivary amylase copy number has been observed in some modern human populations with high starch intake, though evidence of recent selection is lacking. Here, using 52 long-read diploid assemblies and short read data from ~5,600 contemporary and ancient humans, we resolve the diversity, evolutionary history, and selective impact of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in populations with agricultural subsistence compared to fishing, hunting, and pastoral groups. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each exhibit multiple duplications/deletions with mutation rates >10,000-fold the SNP mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome graph-based approach to infer structural haplotypes across thousands of humans, we identify extensively duplicated haplotypes present at higher frequencies in modern day populations with traditionally agricultural diets. Leveraging 533 ancient human genomes we find that duplication-containing haplotypes (i.e. haplotypes with more amylase gene copies than the ancestral haplotype) have increased in frequency more than seven-fold over the last 12,000 years providing evidence for recent selection in West Eurasians. Together, our study highlights the potential impacts of the agricultural revolution on human genomes and the importance of long-read sequencing in identifying signatures of selection at structurally complex loci.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout ageing.

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mitochondrial genome in different tissues throughout ageing. We used ultrasensitive duplex sequencing to profile ~2.5 million mitochondrial genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloguing ~1.2 million mitochondrial somatic and ultralow-frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the light strand origin of replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared with primates with a surfeit of reactive oxygen species-associated G > T/C > A mutations, and that somatic mutations in protein-coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that 're-align' mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

The selection landscape and genetic legacy of ancient Eurasians.

The Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

Reference genome of Townsend's big-eared bat, Corynorhinus townsendii.

Townsend's big-eared bat, Corynorhinus townsendii, is a cave- and mine-roosting species found largely in western North America. Considered a species of conservation concern throughout much of its range, protection efforts would greatly benefit from understanding patterns of population structure, genetic diversity, and local adaptation. To facilitate such research, we present the first de novo genome assembly of C. townsendii as part of the California Conservation Genomics Project (CCGP). Pacific Biosciences HiFi long reads and Omni-C chromatin-proximity sequencing technologies were used to produce a de novo genome assembly, consistent with the standard CCGP reference genome protocol. This assembly comprises 391 scaffolds spanning 2.1 Gb, represented by a scaffold N50 of 174.6 Mb, a contig N50 of 23.4 Mb, and a benchmarking universal single-copy ortholog (BUSCO) completeness score of 96.6%. This high-quality genome will be a key tool for informed conservation and management of this vulnerable species in California and across its range.

A genome assembly of the Yuma myotis bat, Myotis yumanensis.

The Yuma myotis bat (Myotis yumanensis) is a small vespertilionid bat and one of 52 species of new world Myotis bats in the subgenus Pizonyx. While M. yumanensis populations currently appear relatively stable, it is one of 12 bat species known or suspected to be susceptible to white-nose syndrome, the fungal disease causing declines in bat populations across North America. Only two of these 12 species have genome resources available, which limits the ability of resource managers to use genomic techniques to track the responses of bat populations to white-nose syndrome generally. Here we present the first de novo genome assembly for Yuma myotis, generated as a part of the California Conservation Genomics Project. The M. yumanensis genome was generated using a combination of PacBio HiFi long reads and Omni-C chromatin-proximity sequencing technology. This high-quality genome is one of the most complete bat assemblies available, with a contig N50 of 28.03 Mb, scaffold N50 of 99.14 Mb, and BUSCO completeness score of 93.7%. The Yuma myotis genome provides a high-quality resource that will aid in comparative genomic and evolutionary studies, as well as inform conservation management related to white-nose syndrome.

The evolution of aging and lifespan.

Aging is a nearly inescapable trait among organisms yet lifespan varies tremendously across different species and spans several orders of magnitude in vertebrates alone. This vast phenotypic diversity is driven by distinct evolutionary trajectories and tradeoffs that are reflected in patterns of diversification and constraint in organismal genomes. Age-specific impacts of selection also shape allele frequencies in populations, thus impacting disease susceptibility and environment-specific mortality risk. Further, the mutational processes that spawn this genetic diversity in both germline and somatic cells are strongly influenced by age and life history. We discuss recent advances in our understanding of the evolution of aging and lifespan at organismal, population, and cellular scales, and highlight outstanding questions that remain unanswered.

Rapid turnover and evolution of sex-determining regions in Sebastes rockfishes.

Nature has evolved a wealth of sex determination (SD) mechanisms, driven by both genetic and environmental factors. Recent studies of SD in fishes have shown that not all taxa fit the classic paradigm of sex chromosome evolution and diverse SD methods can be found even among closely related species. Here, we apply a suite of genomic approaches to investigate sex-biased genomic variation in eight species of Sebastes rockfish found in the northeast Pacific Ocean. Using recently assembled chromosome-level rockfish genomes, we leverage published sequence data to identify disparate sex chromosomes and sex-biased loci in five species. We identify two putative male sex chromosomes in S. diaconus, a single putative sex chromosome in the sibling species S. carnatus and S. chrysomelas, and an unplaced sex determining contig in the sibling species S. miniatus and S. crocotulus. Our study provides evidence for disparate means of sex determination within a recently diverged set of species and sheds light on the diverse origins of sex determination mechanisms present in the animal kingdom.

Hypoxia blunts angiogenic signaling and upregulates the antioxidant system in elephant seal endothelial cells.

Elephant seals experience extreme hypoxemia during diving bouts. Similar depletions in oxygen availability characterize pathologies including myocardial infarction and ischemic stroke in humans, but seals manage these repeated episodes without injury. However, the real-time assessment of the molecular changes underlying protection against hypoxic injury in seals remains restricted by their at-sea inaccessibility. Hence, we developed a proliferative arterial endothelial cell culture system to assess the molecular response to prolonged hypoxia. Seal and human cells exposed to 1% O 2 for up to 6 h demonstrated differential responses to both acute and prolonged hypoxia. Seal cells decouple stabilization of the hypoxia-sensitive transcriptional regulator HIF-1α from angiogenic signaling at both the transcriptional and cellular level. Rapid upregulation of genes involved in the glutathione (GSH) metabolism pathway supported maintenance of GSH pools and increases in intracellular succinate in seal but not human cells during hypoxia exposure. High maximal and spare respiratory capacity in seal cells after hypoxia exposure occurred in concert with increasing mitochondrial branch length and independent from major changes in extracellular acidification rate, suggesting seal cells recover oxidative metabolism without significant glycolytic dependency after hypoxia exposure. In sum, our studies show that in contrast to human cells, seal cells adapt to hypoxia exposure by dampening angiogenic signaling, increasing antioxidant protection, and maintaining mitochondrial morphological integrity and function.

North African fox genomes show signatures of repeated introgression and adaptation to life in deserts.

Elucidating the evolutionary process of animal adaptation to deserts is key to understanding adaptive responses to climate change. Here we generated 82 individual whole genomes of four fox species (genus Vulpes) inhabiting the Sahara Desert at different evolutionary times. We show that adaptation of new colonizing species to a hot arid environment has probably been facilitated by introgression and trans-species polymorphisms shared with older desert resident species, including a putatively adaptive 25 Mb genomic region. Scans for signatures of selection implicated genes affecting temperature perception, non-renal water loss and heat production in the recent adaptation of North African red foxes (Vulpes vulpes), after divergence from Eurasian populations approximately 78 thousand years ago. In the extreme desert specialists, Rueppell's fox (V. rueppellii) and fennec (V. zerda), we identified repeated signatures of selection in genes affecting renal water homeostasis supported by gene expression and physiological differences. Our study provides insights into the mechanisms and genetic underpinnings of a natural experiment of repeated adaptation to extreme conditions.

Multiscale analysis of pangenomes enables improved representation of genomic diversity for repetitive and clinically relevant genes.

Advancements in sequencing technologies and assembly methods enable the regular production of high-quality genome assemblies characterizing complex regions. However, challenges remain in efficiently interpreting variation at various scales, from smaller tandem repeats to megabase rearrangements, across many human genomes. We present a PanGenome Research Tool Kit (PGR-TK) enabling analyses of complex pangenome structural and haplotype variation at multiple scales. We apply the graph decomposition methods in PGR-TK to the class II major histocompatibility complex demonstrating the importance of the human pangenome for analyzing complicated regions. Moreover, we investigate the Y-chromosome genes, DAZ1/DAZ2/DAZ3/DAZ4, of which structural variants have been linked to male infertility, and X-chromosome genes OPN1LW and OPN1MW linked to eye disorders. We further showcase PGR-TK across 395 complex repetitive medically important genes. This highlights the power of PGR-TK to resolve complex variation in regions of the genome that were previously too complex to analyze.

Highly-multiplexed and efficient long-amplicon PacBio and Nanopore sequencing of hundreds of full mitochondrial genomes.

BACKGROUND: Mitochondrial genome sequences have become critical to the study of biodiversity. Genome skimming and other short-read based methods are the most common approaches, but they are not well-suited to scale up to multiplexing hundreds of samples. Here, we report on a new approach to sequence hundreds to thousands of complete mitochondrial genomes in parallel using long-amplicon sequencing. We amplified the mitochondrial genome of 677 specimens in two partially overlapping amplicons and implemented an asymmetric PCR-based indexing approach to multiplex 1,159 long amplicons together on a single PacBio SMRT Sequel II cell. We also tested this method on Oxford Nanopore Technologies (ONT) MinION R9.4 to assess if this method could be applied to other long-read technologies. We implemented several optimizations that make this method significantly more efficient than alternative mitochondrial genome sequencing methods. RESULTS: With the PacBio sequencing data we recovered at least one of the two fragments for 96% of samples (~ 80-90%) with mean coverage ~ 1,500x. The ONT data recovered less than 50% of input fragments likely due to low throughput and the design of the Barcoded Universal Primers which were optimized for PacBio sequencing. We compared a single mitochondrial gene alignment to half and full mitochondrial genomes and found, as expected, increased tree support with longer alignments, though whole mitochondrial genomes were not significantly better than half mitochondrial genomes. CONCLUSIONS: This method can effectively capture thousands of long amplicons in a single run and be used to build more robust phylogenies quickly and effectively. We provide several recommendations for future users depending on the evolutionary scale of their system. A natural extension of this method is to collect multi-locus datasets consisting of mitochondrial genomes and several long nuclear loci at once.

The mitochondrial unfolded protein response regulates hippocampal neural stem cell aging.

Aging results in a decline in neural stem cells (NSCs), neurogenesis, and cognitive function, and evidence is emerging to demonstrate disrupted adult neurogenesis in the hippocampus of patients with several neurodegenerative disorders. Here, single-cell RNA sequencing of the dentate gyrus of young and old mice shows that the mitochondrial protein folding stress is prominent in activated NSCs/neural progenitors (NPCs) among the neurogenic niche, and it increases with aging accompanying dysregulated cell cycle and mitochondrial activity in activated NSCs/NPCs in the dentate gyrus. Increasing mitochondrial protein folding stress results in compromised NSC maintenance and reduced neurogenesis in the dentate gyrus, neural hyperactivity, and impaired cognitive function. Reducing mitochondrial protein folding stress in the dentate gyrus of old mice improves neurogenesis and cognitive function. These results establish the mitochondrial protein folding stress as a driver of NSC aging and suggest approaches to improve aging-associated cognitive decline.

Building pangenome graphs.

Pangenome graphs can represent all variation between multiple genomes, but existing methods for constructing them are biased due to reference-guided approaches. In response, we have developed PanGenome Graph Builder (PGGB), a reference-free pipeline for constructing unbi-ased pangenome graphs. PGGB uses all-to-all whole-genome alignments and learned graph embeddings to build and iteratively refine a model in which we can identify variation, measure conservation, detect recombination events, and infer phylogenetic relationships.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout aging.

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here, we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mt-genome in different tissues throughout aging. We used ultra-sensitive Duplex Sequencing to profile ~2.5 million mt-genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloging ~1.2 million mitochondrial somatic and ultra low frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the Light Strand Origin of Replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared to primates with a surfeit of reactive oxygen species-associated G>T/C>A mutations, and that somatic mutations in protein coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that "re-align" mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Transcriptomic profiles of stress susceptibility and resilience in the amygdala and hippocampus.

A single, severe episode of stress can bring about myriad responses amongst individuals, ranging from cognitive enhancement to debilitating and persistent anxiety; however, the biological mechanisms that contribute to resilience versus susceptibility to stress are poorly understood. The dentate gyrus (DG) of the hippocampus and the basolateral nucleus of the amygdala (BLA) are key limbic regions that are susceptible to the neural and hormonal effects of stress. Previous work has also shown that these regions contribute to individual variability in stress responses; however, the molecular mechanisms underlying the role of these regions in susceptibility and resilience are unknown. In this study, we profiled the transcriptomic signatures of the DG and BLA of rats with divergent behavioral outcomes after a single, severe stressor. We subjected rats to three hours of immobilization with exposure to fox urine and conducted a behavioral battery one week after stress to identify animals that showed persistent, high anxiety-like behavior. We then conducted bulk RNA sequencing of the DG and BLA from susceptible, resilient, and unexposed control rats. Differential gene expression analyses revealed that the molecular signatures separating each of the three groups were distinct and non-overlapping between the DG and BLA. In the amygdala, key genes associated with insulin and hormonal signaling corresponded with vulnerability. Specifically, Inhbb, Rab31 , and Ncoa3 were upregulated in the amygdala of stress-susceptible animals compared to resilient animals. In the hippocampus, increased expression of Cartpt - which encodes a key neuropeptide involved in reward, reinforcement, and stress responses - was strongly correlated with vulnerability to anxiety-like behavior. However, few other genes distinguished stress-susceptible animals from control animals, while a larger number of genes separated stress-resilient animals from control and stress-susceptible animals. Of these, Rnf112, Tbx19 , and UBALD1 distinguished resilient animals from both control and susceptible animals and were downregulated in resilience, suggesting that an active molecular response in the hippocampus facilitates protection from the long-term consequences of severe stress. These results provide novel insight into the mechanisms that bring about individual variability in the behavioral responses to stress and provide new targets for the advancement of therapies for stress-induced neuropsychiatric disorders.