Search

©2014 Macmillan Publishers Limited. All rights reserved.We sequenced the genomes of a ∼7,000-year-old farmer from Germany and eight ∼8,000-year-old hunter-gatherers from Luxembourg and Sweden. We analysed these and other ancient genomes1-4 with

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder, a pipeline for constructing pangenome graphs without bias or exclusion. The PanGenome Graph Builder uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events and infer phylogenetic relationships., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

The rich species diversity of bats encompasses extraordinary adaptations, including extreme longevity and tolerance to infectious disease. While traditional approaches using genetic screens in model organisms have uncovered some fundamental processes underlying these traits, model organisms do not possess the variation required to understand the evolution of traits with complex genetic architectures. In contrast, the advent of genomics at tree-of-life scales enables us to study the genetic interactions underlying these processes by leveraging millions of years of evolutionary trial-and-error. Here, we use the rich species diversity of the genus Myotis - one of the longest-living clades of mammals - to study the evolution of longevity-associated traits and infectious disease using functional evolutionary genomics. We generated reference genome assemblies and cell lines for 8 closely-related (~11 MYA) species of Myotis rich in phenotypic and life history diversity. Using genome-wide screens of positive selection, analysis of structural variation and copy number variation, and functional experiments in primary cell lines, we identify new patterns of adaptation in longevity, cancer resistance, and viral interactions both within Myotis and across bats. We find that the rapid evolution of lifespan in Myotis has some of the most significant variations in cancer risk across mammals, and demonstrate a unique DNA damage response in the long-lived M. lucifugus using primary cell culture models. Furthermore, we find evidence of abundant adaptation in response to DNA viruses, but not RNA viruses, in Myotis and other bats. This is in contrast to these patterns of adaptation in humans, which might contribute to the importance of bats as a reservoir of zoonotic viruses. Together, our results demonstrate the utility of leveraging natural variation to understand the genomics of traits with implications for human health and suggest important pleiotropic relationships between infectious disease tolerance and cancer resistance., Competing Interests: Declaration of Interests: The authors declare no competing interests.

We present haplotype-resolved reference genomes and comparative analyses of six ape species, namely: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. We achieve chromosome-level contiguity with unparalleled sequence accuracy (<1 error in 500,000 base pairs), completely sequencing 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, providing more in-depth evolutionary insights. Comparative analyses, including human, allow us to investigate the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference. This includes newly minted gene families within lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes, and subterminal heterochromatin. This resource should serve as a definitive baseline for all future evolutionary studies of humans and our closest living ape relatives., Competing Interests: COMPETING INTERESTS E.E.E. is a scientific advisory board (SAB) member of Variant Bio, Inc. C.T.W. is a co-founder/CSO of Clareo Biosciences, Inc. W.L. is a co-founder/CIO of Clareo Biosciences, Inc. The other authors declare no competing interests.

The adoption of agriculture triggered a rapid shift towards starch-rich diets in human populations 1 . Amylase genes facilitate starch digestion, and increased amylase copy number has been observed in some modern human populations with high-starch intake 2 , although evidence of recent selection is lacking 3,4 . Here, using 94 long-read haplotype-resolved assemblies and short-read data from approximately 5,600 contemporary and ancient humans, we resolve the diversity and evolutionary history of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in agricultural populations than in fishing, hunting and pastoral populations. 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 underwent multiple duplication/deletion events with mutation rates up to more than 10,000-fold the single-nucleotide polymorphism mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome-based approach, we infer structural haplotypes across thousands of humans identifying extensively duplicated haplotypes at higher frequency in modern agricultural populations. Leveraging 533 ancient human genomes, we find that duplication-containing haplotypes (with more gene copies than the ancestral haplotype) have rapidly increased in frequency over the past 12,000 years in West Eurasians, suggestive of positive selection. Together, our study highlights the potential effects of the agricultural revolution on human genomes and the importance of structural variation in human adaptation., (© 2024. The Author(s).)

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., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interest., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

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., Competing Interests: ES, UM, LY, LS, LK, SM, IB, SU, RP, PS, OP, DK, AF No competing interests declared, (© 2023, Swissa et al.)

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)., Competing Interests: None declared., (© The Author(s) 2024. Published by Oxford University Press.)

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., Competing Interests: Competing interests S.Lee is a shareholder and an employee at Pacific Biosciences. P.J.C. is a co-founder, shareholder, and consultant for Quotient Therapeutics. R.S is a consultant to MyHeritage Ltd.

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.

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility 1 . The X chromosome is vital for reproduction and cognition 2 . Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

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% O 2 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., (© 2024. The Author(s).)

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., (© The American Genetic Association. 2023.)

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 American Genetic Association. 2023.)

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 genomes 1 , 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., (© 2024. The Author(s).)

Western Eurasia witnessed several large-scale human migrations during the Holocene 1-5 . Here, to investigate the cross-continental effects of these migrations, we shotgun-sequenced 317 genomes-mainly from the Mesolithic and Neolithic periods-from across northern and western Eurasia. These were imputed alongside published data to obtain diploid genotypes from more than 1,600 ancient humans. Our analyses revealed a 'great divide' genomic boundary extending from the Black Sea to the Baltic. Mesolithic hunter-gatherers were highly genetically differentiated east and west of this zone, and the effect of the neolithization was equally disparate. Large-scale ancestry shifts occurred in the west as farming was introduced, including near-total replacement of hunter-gatherers in many areas, whereas no substantial ancestry shifts happened east of the zone during the same period. Similarly, relatedness decreased in the west from the Neolithic transition onwards, whereas, east of the Urals, relatedness remained high until around 4,000 BP, consistent with the persistence of localized groups of hunter-gatherers. The boundary dissolved when Yamnaya-related ancestry spread across western Eurasia around 5,000 BP, resulting in a second major turnover that reached most parts of Europe within a 1,000-year span. The genetic origin and fate of the Yamnaya have remained elusive, but we show that hunter-gatherers from the Middle Don region contributed ancestry to them. Yamnaya groups later admixed with individuals associated with the Globular Amphora culture before expanding into Europe. Similar turnovers occurred in western Siberia, where we report new genomic data from a 'Neolithic steppe' cline spanning the Siberian forest steppe to Lake Baikal. These prehistoric migrations had profound and lasting effects on the genetic diversity of Eurasian populations., (© 2024. The Author(s).)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

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., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

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., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

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., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Elsevier Inc. All rights reserved.)

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., (© 2023. The Author(s).)

Nonalcoholic fatty liver disease (NAFLD) can be ameliorated by calorie restriction, which leads to the suppressed somatotroph axis. Paradoxically, the suppressed somatotroph axis is associated with patients with NAFLD and is correlated with the severity of fibrosis. How the somatotroph axis becomes dysregulated and whether the repressed somatotroph axis impacts liver damage during the progression of NAFLD are unclear. Here, we identify a regulatory branch of the hepatic integrated stress response (ISR), which represses the somatotroph axis in hepatocytes through ATF3, resulting in enhanced cell survival and reduced cell proliferation. In mouse models of NAFLD, the ISR represses the somatotroph axis, leading to reduced apoptosis and inflammation but decreased hepatocyte proliferation and exacerbated fibrosis in the liver. NAD + repletion reduces the ISR, rescues the dysregulated somatotroph axis, and alleviates NAFLD. These results establish that the hepatic ISR suppresses the somatotroph axis to control cell fate decisions and liver damage in NAFLD., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Age is the primary risk factor for many common human diseases. Here, we quantify the relative contributions of genetics and aging to gene expression patterns across 27 tissues from 948 humans. We show that the predictive power of expression quantitative trait loci is impacted by age in many tissues. Jointly modelling the contributions of age and genetics to transcript level variation we find expression heritability (h 2 ) is consistent among tissues while the contribution of aging varies by >20-fold with [Formula: see text] in 5 tissues. We find that while the force of purifying selection is stronger on genes expressed early versus late in life (Medawar's hypothesis), several highly proliferative tissues exhibit the opposite pattern. These non-Medawarian tissues exhibit high rates of cancer and age-of-expression-associated somatic mutations. In contrast, genes under genetic control are under relaxed constraint. Together, we demonstrate the distinct roles of aging and genetics on expression phenotypes., (© 2022. The Author(s).)

The variegated toad-headed agama, Phrynocephalus versicolor, lives in the arid landscape of the Chinese Gobi Desert. We analyzed populations from three different locations which vary in substrate color and altitude: Heishankou (HSK), Guazhou County (GZ), and Ejin Banner (EJN). The substrate color is either light-yellow (GZ-y), yellow (EJN-y), or black (HSK-b); the corresponding lizard population colors largely match their substrate in the degree of melanism. We assembled the P. versicolor genome and sequenced over 90 individuals from the three different populations. Genetic divergence between populations corresponds to their geographic distribution. We inferred the genetic relationships among these populations and used selection scans and differential expression to identify genes that show signatures of selection. Slc2a11 and akap12, among other genes, are highly differentiated and may be responsible for pigment adaptation to substrate color in P. versicolor., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Gametogenesis is an evolutionarily conserved developmental program whereby a diploid progenitor cell undergoes meiosis and cellular remodeling to differentiate into haploid gametes, the precursors for sexual reproduction. Even in the simple eukaryotic organism Saccharomyces cerevisiae, the meiotic transcriptome is very rich and complex, thereby necessitating new tools for functional studies. Here, we report the construction of 5 stage-specific, inducible complementary DNA libraries from meiotic cells that represent over 84% of the genes found in the budding yeast genome. We employed computational strategies to detect endogenous meiotic transcript isoforms as well as library-specific gene truncations. Furthermore, we developed a robust screening pipeline to test the effect of each complementary DNA on competitive fitness. Our multiday proof-of-principle time course revealed 877 complementary DNAs that were detrimental for competitive fitness when overexpressed. The list included mitochondrial proteins that cause dose-dependent disruption of cellular respiration as well as library-specific gene truncations that expose a dominant negative effect on competitive growth. Together, these high-quality complementary DNA libraries provide an important tool for systematically identifying meiotic genes, transcript isoforms, and protein domains that are important for a specific biological function., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background: Cellular senescence is a complex stress response that impacts cellular function and organismal health. Multiple developmental and environmental factors, such as intrinsic cellular cues, radiation, oxidative stress, oncogenes, and protein accumulation, activate genes and pathways that can lead to senescence. Enormous efforts have been made to identify and characterize senescence genes (SnGs) in stress and disease systems. However, the prevalence of senescent cells in healthy human tissues and the global SnG expression signature in different cell types are poorly understood., Methods: This study performed an integrative gene network analysis of bulk and single-cell RNA-seq data in non-diseased human tissues to investigate SnG co-expression signatures and their cell-type specificity., Results: Through a comprehensive transcriptomic network analysis of 50 human tissues in the Genotype-Tissue Expression Project (GTEx) cohort, we identified SnG-enriched gene modules, characterized SnG co-expression patterns, and constructed aggregated SnG networks across primary tissues of the human body. Our network approaches identified 51 SnGs highly conserved across the human tissues, including CDKN1A (p21)-centered regulators that control cell cycle progression and the senescence-associated secretory phenotype (SASP). The SnG-enriched modules showed remarkable cell-type specificity, especially in fibroblasts, endothelial cells, and immune cells. Further analyses of single-cell RNA-seq and spatial transcriptomic data independently validated the cell-type specific SnG signatures predicted by the network analysis., Conclusions: This study systematically revealed the co-regulated organizations and cell type specificity of SnGs in major human tissues, which can serve as a blueprint for future studies to map senescent cells and their cellular interactions in human tissues., (© 2021. The Author(s).)

Pacific Ocean rockfishes (genus Sebastes ) exhibit extreme variation in life span, with some species being among the most long-lived extant vertebrates. We de novo assembled the genomes of 88 rockfish species and from these identified repeated signatures of positive selection in DNA repair pathways in long-lived taxa and 137 longevity-associated genes with direct effects on life span through insulin signaling and with pleiotropic effects through size and environmental adaptations. A genome-wide screen of structural variation reveals copy number expansions in the immune modulatory butyrophilin gene family in long-lived species. The evolution of different rockfish life histories is coupled to genetic diversity and reshapes the mutational spectrum driving segregating CpG→TpG variants in long-lived species. These analyses highlight the genetic innovations that underlie life history trait adaptations and, in turn, how they shape genomic diversity.

Elephant seals experience natural periods of prolonged food deprivation while breeding, molting, and undergoing postnatal development. Prolonged food deprivation in elephant seals increases circulating glucocorticoids without inducing muscle atrophy, but the cellular mechanisms that allow elephant seals to cope with such conditions remain elusive. We generated a cellular model and conducted transcriptomic, metabolic, and morphological analyses to study how seal cells adapt to sustained glucocorticoid exposure. Seal muscle progenitor cells differentiate into contractile myotubes with a distinctive morphology, gene expression profile, and metabolic phenotype. Exposure to dexamethasone at three ascending concentrations for 48 h modulated the expression of six clusters of genes related to structural constituents of muscle and pathways associated with energy metabolism and cell survival. Knockdown of the glucocorticoid receptor (GR) and downstream expression analyses corroborated that GR mediates the observed effects. Dexamethasone also decreased cellular respiration, shifted the metabolic phenotype toward glycolysis, and induced mitochondrial fission and dissociation of mitochondria-endoplasmic reticulum (ER) interactions without decreasing cell viability. Knockdown of DNA damage-inducible transcript 4 (DDIT4), a GR target involved in the dissociation of mitochondria-ER membranes, recovered respiration and modulated antioxidant gene expression in myotubes treated with dexamethasone. These results show that adaptation to sustained glucocorticoid exposure in elephant seal myotubes involves a metabolic shift toward glycolysis, which is supported by alterations in mitochondrial morphology and a reduction in mitochondria-ER interactions, resulting in decreased respiration without compromising cell survival.

A progressive loss of protein homeostasis is characteristic of aging and a driver of neurodegeneration. To investigate this process quantitatively, we characterized proteome dynamics during brain aging in the short-lived vertebrate Nothobranchius furzeri combining transcriptomics and proteomics. We detected a progressive reduction in the correlation between protein and mRNA, mainly due to post-transcriptional mechanisms that account for over 40% of the age-regulated proteins. These changes cause a progressive loss of stoichiometry in several protein complexes, including ribosomes, which show impaired assembly/disassembly and are enriched in protein aggregates in old brains. Mechanistically, we show that reduction of proteasome activity is an early event during brain aging and is sufficient to induce proteomic signatures of aging and loss of stoichiometry in vivo. Using longitudinal transcriptomic data, we show that the magnitude of early life decline in proteasome levels is a major risk factor for mortality. Our work defines causative events in the aging process that can be targeted to prevent loss of protein homeostasis and delay the onset of age-related neurodegeneration., (©2020 The Authors. Published under the terms of the CC BY 4.0 license.)

Aging involves a decline in neural function that contributes to cognitive impairment and disease. However, the mechanisms underlying the transition from a young-and-healthy to aged-and-dysfunctional brain are not well understood. Here, we report breakdown of the vascular blood-brain barrier (BBB) in aging humans and rodents, which begins as early as middle age and progresses to the end of the life span. Gain-of-function and loss-of-function manipulations show that this BBB dysfunction triggers hyperactivation of transforming growth factor-β (TGFβ) signaling in astrocytes, which is necessary and sufficient to cause neural dysfunction and age-related pathology in rodents. Specifically, infusion of the serum protein albumin into the young rodent brain (mimicking BBB leakiness) induced astrocytic TGFβ signaling and an aged brain phenotype including aberrant electrocorticographic activity, vulnerability to seizures, and cognitive impairment. Furthermore, conditional genetic knockdown of astrocytic TGFβ receptors or pharmacological inhibition of TGFβ signaling reversed these symptomatic outcomes in aged mice. Last, we found that this same signaling pathway is activated in aging human subjects with BBB dysfunction. Our study identifies dysfunction in the neurovascular unit as one of the earliest triggers of neurological aging and demonstrates that the aging brain may retain considerable latent capacity, which can be revitalized by therapeutic inhibition of TGFβ signaling., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Particular brain regions and cell populations exhibit increased susceptibility to aging-related stresses. Here, we describe the age-specific and brain-region-specific accumulation of ribosome-associated 3' UTR RNAs that lack the 5' UTR and open reading frame. Our study reveals that this phenomenon impacts hundreds of genes in aged D1 spiny projection neurons of the mouse striatum and also occurs in the aging human brain. Isolated 3' UTR accumulation is tightly correlated with mitochondrial gene expression and oxidative stress, with full-length mRNA expression that is reduced but not eliminated, and with production of short 3' UTR-encoded peptides. Depletion of the oxidation-sensitive Fe-S cluster ribosome recycling factor ABCE1 induces the accumulation of 3' UTRs, consistent with a model in which ribosome stalling and mRNA cleavage by No-Go decay yields isolated 3' UTR RNAs protected by ribosomes. Isolated 3' UTR accumulation is a hallmark of brain aging, likely reflecting regional differences in metabolism and oxidative stress., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Many RNA binding proteins (RBPs) bind specific RNA sequence motifs, but only a small fraction (∼15%-40%) of RBP motif occurrences are occupied in vivo. To determine which contextual features discriminate between bound and unbound motifs, we performed an in vitro binding assay using 12,000 mouse RNA sequences with the RBPs MBNL1 and RBFOX2. Surprisingly, the strength of binding to motif occurrences in vitro was significantly correlated with in vivo binding, developmental regulation, and evolutionary age of alternative splicing. Multiple lines of evidence indicate that the primary context effect that affects binding in vitro and in vivo is RNA secondary structure. Large-scale combinatorial mutagenesis of unfavorable sequence contexts revealed a consistent pattern whereby mutations that increased motif accessibility improved protein binding and regulatory activity. Our results indicate widespread inhibition of motif binding by local RNA secondary structure and suggest that mutations that alter sequence context commonly affect RBP binding and regulation., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Genetic differences that specify unique aspects of human evolution have typically been identified by comparative analyses between the genomes of humans and closely related primates, including more recently the genomes of archaic hominins. Not all regions of the genome, however, are equally amenable to such study. Recurrent copy number variation (CNV) at chromosome 16p11.2 accounts for approximately 1% of cases of autism and is mediated by a complex set of segmental duplications, many of which arose recently during human evolution. Here we reconstruct the evolutionary history of the locus and identify bolA family member 2 (BOLA2) as a gene duplicated exclusively in Homo sapiens. We estimate that a 95-kilobase-pair segment containing BOLA2 duplicated across the critical region approximately 282 thousand years ago (ka), one of the latest among a series of genomic changes that dramatically restructured the locus during hominid evolution. All humans examined carried one or more copies of the duplication, which nearly fixed early in the human lineage--a pattern unlikely to have arisen so rapidly in the absence of selection (P < 0.0097). We show that the duplication of BOLA2 led to a novel, human-specific in-frame fusion transcript and that BOLA2 copy number correlates with both RNA expression (r = 0.36) and protein level (r = 0.65), with the greatest expression difference between human and chimpanzee in experimentally derived stem cells. Analyses of 152 patients carrying a chromosome 16p11. rearrangement show that more than 96% of breakpoints occur within the H. sapiens-specific duplication. In summary, the duplicative transposition of BOLA2 at the root of the H. sapiens lineage about 282 ka simultaneously increased copy number of a gene associated with iron homeostasis and predisposed our species to recurrent rearrangements associated with disease.

Background: Differences in gene expression drive phenotypic differences between species, yet major organs and tissues generally have conserved gene expression programs. Several comparative transcriptomic studies have observed greater similarity in gene expression between homologous tissues from different vertebrate species than between diverse tissues of the same species. However, a recent study by Lin and colleagues reached the opposite conclusion. These studies differed in the species and tissues analyzed, and in technical details of library preparation, sequencing, read mapping, normalization, gene sets, and clustering methods., Results: To better understand gene expression evolution we reanalyzed data from four studies, including that of Lin, encompassing 6-13 tissues each from 11 vertebrate species using standardized mapping, normalization, and clustering methods. An analysis of independent data showed that the set of tissues chosen by Lin et al. were more similar to each other than those analyzed by previous studies. Comparing expression in five common tissues from the four studies, we observed that samples clustered exclusively by tissue rather than by species or study, supporting conservation of organ physiology in mammals. Furthermore, inter-study distances between homologous tissues were generally less than intra-study distances among different tissues, enabling informative meta-analyses. Notably, when comparing expression divergence of tissues over time to expression variation across 51 human GTEx tissues, we could accurately predict the clustering of expression for arbitrary pairs of tissues and species., Conclusions: These results provide a framework for the design of future evolutionary studies of gene expression and demonstrate the utility of comparing RNA-seq data across studies.

Structural variants are implicated in numerous diseases and make up the majority of varying nucleotides among human genomes. Here we describe an integrated set of eight structural variant classes comprising both balanced and unbalanced variants, which we constructed using short-read DNA sequencing data and statistically phased onto haplotype blocks in 26 human populations. Analysing this set, we identify numerous gene-intersecting structural variants exhibiting population stratification and describe naturally occurring homozygous gene knockouts that suggest the dispensability of a variety of human genes. We demonstrate that structural variants are enriched on haplotypes identified by genome-wide association studies and exhibit enrichment for expression quantitative trait loci. Additionally, we uncover appreciable levels of structural variant complexity at different scales, including genic loci subject to clusters of repeated rearrangement and complex structural variants with multiple breakpoints likely to have formed through individual mutational events. Our catalogue will enhance future studies into structural variant demography, functional impact and disease association.

Cross-species transmissions of viruses from animals to humans are at the origin of major human pathogenic viruses. While the role of ecological and epidemiological factors in the emergence of new pathogens is well documented, the importance of host factors is often unknown. Chimpanzees are the closest relatives of humans and the animal reservoir at the origin of the human AIDS pandemic. However, despite being regularly exposed to monkey lentiviruses through hunting, chimpanzees are naturally infected by only a single simian immunodeficiency virus, SIVcpz. Here, we asked why chimpanzees appear to be protected against the successful emergence of other SIVs. In particular, we investigated the role of the chimpanzee APOBEC3 genes in providing a barrier to infection by most monkey lentiviruses. We found that most SIV Vifs, including Vif from SIVwrc infecting western-red colobus, the chimpanzee's main monkey prey in West Africa, could not antagonize chimpanzee APOBEC3G. Moreover, chimpanzee APOBEC3D, as well as APOBEC3F and APOBEC3H, provided additional protection against SIV Vif antagonism. Consequently, lentiviral replication in primary chimpanzee CD4(+) T cells was dependent on the presence of a lentiviral vif gene that could antagonize chimpanzee APOBEC3s. Finally, by identifying and functionally characterizing several APOBEC3 gene polymorphisms in both common chimpanzees and bonobos, we found that these ape populations encode APOBEC3 proteins that are uniformly resistant to antagonism by monkey lentiviruses.

In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide-variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load., (Copyright © 2015, American Association for the Advancement of Science.)

Short tandem repeats (STRs) are highly mutable genetic elements that often reside in regulatory and coding DNA. The cumulative evidence of genetic studies on individual STRs suggests that STR variation profoundly affects phenotype and contributes to trait heritability. Despite recent advances in sequencing technology, STR variation has remained largely inaccessible across many individuals compared to single nucleotide variation or copy number variation. STR genotyping with short-read sequence data is confounded by (1) the difficulty of uniquely mapping short, low-complexity reads; and (2) the high rate of STR amplification stutter. Here, we present MIPSTR, a robust, scalable, and affordable method that addresses these challenges. MIPSTR uses targeted capture of STR loci by single-molecule Molecular Inversion Probes (smMIPs) and a unique mapping strategy. Targeted capture and our mapping strategy resolve the first challenge; the use of single molecule information resolves the second challenge. Unlike previous methods, MIPSTR is capable of distinguishing technical error due to amplification stutter from somatic STR mutations. In proof-of-principle experiments, we use MIPSTR to determine germline STR genotypes for 102 STR loci with high accuracy across diverse populations of the plant A. thaliana. We show that putatively functional STRs may be identified by deviation from predicted STR variation and by association with quantitative phenotypes. Using DNA mixing experiments and a mutant deficient in DNA repair, we demonstrate that MIPSTR can detect low-frequency somatic STR variants. MIPSTR is applicable to any organism with a high-quality reference genome and is scalable to genotyping many thousands of STR loci in thousands of individuals., (© 2015 Carlson et al.; Published by Cold Spring Harbor Laboratory Press.)

Mountain gorillas are an endangered great ape subspecies and a prominent focus for conservation, yet we know little about their genomic diversity and evolutionary past. We sequenced whole genomes from multiple wild individuals and compared the genomes of all four Gorilla subspecies. We found that the two eastern subspecies have experienced a prolonged population decline over the past 100,000 years, resulting in very low genetic diversity and an increased overall burden of deleterious variation. A further recent decline in the mountain gorilla population has led to extensive inbreeding, such that individuals are typically homozygous at 34% of their sequence, leading to the purging of severely deleterious recessive mutations from the population. We discuss the causes of their decline and the consequences for their future survival., (Copyright © 2015, American Association for the Advancement of Science.)

The human genome is arguably the most complete mammalian reference assembly, yet more than 160 euchromatic gaps remain and aspects of its structural variation remain poorly understood ten years after its completion. To identify missing sequence and genetic variation, here we sequence and analyse a haploid human genome (CHM1) using single-molecule, real-time DNA sequencing. We close or extend 55% of the remaining interstitial gaps in the human GRCh37 reference genome--78% of which carried long runs of degenerate short tandem repeats, often several kilobases in length, embedded within (G+C)-rich genomic regions. We resolve the complete sequence of 26,079 euchromatic structural variants at the base-pair level, including inversions, complex insertions and long tracts of tandem repeats. Most have not been previously reported, with the greatest increases in sensitivity occurring for events less than 5 kilobases in size. Compared to the human reference, we find a significant insertional bias (3:1) in regions corresponding to complex insertions and long short tandem repeats. Our results suggest a greater complexity of the human genome in the form of variation of longer and more complex repetitive DNA that can now be largely resolved with the application of this longer-read sequencing technology.

Recurrent deletions of chromosome 15q13.3 associate with intellectual disability, schizophrenia, autism and epilepsy. To gain insight into the instability of this region, we sequenced it in affected individuals, normal individuals and nonhuman primates. We discovered five structural configurations of the human chromosome 15q13.3 region ranging in size from 2 to 3 Mb. These configurations arose recently (∼0.5-0.9 million years ago) as a result of human-specific expansions of segmental duplications and two independent inversion events. All inversion breakpoints map near GOLGA8 core duplicons-a ∼14-kb primate-specific chromosome 15 repeat that became organized into larger palindromic structures. GOLGA8-flanked palindromes also demarcate the breakpoints of recurrent 15q13.3 microdeletions, the expansion of chromosome 15 segmental duplications in the human lineage and independent structural changes in apes. The significant clustering (P = 0.002) of breakpoints provides mechanistic evidence for the role of this core duplicon and its palindromic architecture in promoting the evolutionary and disease-related instability of chromosome 15.