Temporal dynamics of woolly mammoth genome erosion prior to extinction.

A number of species have recently recovered from near-extinction. Although these species have avoided the immediate extinction threat, their long-term viability remains precarious due to the potential genetic consequences of population declines, which are poorly understood on a timescale beyond a few generations. Woolly mammoths (Mammuthus primigenius) became isolated on Wrangel Island around 10,000 years ago and persisted for over 200 generations before becoming extinct around 4,000 years ago. To study the evolutionary processes leading up to the mammoths' extinction, we analyzed 21 Siberian woolly mammoth genomes. Our results show that the population recovered quickly from a severe bottleneck and remained demographically stable during the ensuing six millennia. We find that mildly deleterious mutations gradually accumulated, whereas highly deleterious mutations were purged, suggesting ongoing inbreeding depression that lasted for hundreds of generations. The time-lag between demographic and genetic recovery has wide-ranging implications for conservation management of recently bottlenecked populations.

Identifying Epistasis in Cancer Genomes: A Delicate Affair.

Recent studies of the tumor genome seek to identify cancer pathways as groups of genes in which mutations are epistatic with one another or, specifically, "mutually exclusive." Here, we show that most mutations are mutually exclusive not due to pathway structure but to interactions with disease subtype and tumor mutation load. In particular, many cancer driver genes are mutated preferentially in tumors with few mutations overall, causing mutations in these cancer genes to appear mutually exclusive with numerous others. Researchers should view current epistasis maps with caution until we better understand the multiple cause-and-effect relationships among factors such as tumor subtype, positive selection for mutations, and gross tumor characteristics including mutational signatures and load.

UVB-Induced Tumor Heterogeneity Diminishes Immune Response in Melanoma.

Although clonal neo-antigen burden is associated with improved response to immune therapy, the functional basis for this remains unclear. Here we study this question in a novel controlled mouse melanoma model that enables us to explore the effects of intra-tumor heterogeneity (ITH) on tumor aggressiveness and immunity independent of tumor mutational burden. Induction of UVB-derived mutations yields highly aggressive tumors with decreased anti-tumor activity. However, single-cell-derived tumors with reduced ITH are swiftly rejected. Their rejection is accompanied by increased T cell reactivity and a less suppressive microenvironment. Using phylogenetic analyses and mixing experiments of single-cell clones, we dissect two characteristics of ITH: the number of clones forming the tumor and their clonal diversity. Our analysis of melanoma patient tumor data recapitulates our results in terms of overall survival and response to immune checkpoint therapy. These findings highlight the importance of clonal mutations in robust immune surveillance and the need to quantify patient ITH to determine the response to checkpoint blockade.

Tumor and Microenvironment Evolution during Immunotherapy with Nivolumab.

The mechanisms by which immune checkpoint blockade modulates tumor evolution during therapy are unclear. We assessed genomic changes in tumors from 68 patients with advanced melanoma, who progressed on ipilimumab or were ipilimumab-naive, before and after nivolumab initiation (CA209-038 study). Tumors were analyzed by whole-exome, transcriptome, and/or T cell receptor (TCR) sequencing. In responding patients, mutation and neoantigen load were reduced from baseline, and analysis of intratumoral heterogeneity during therapy demonstrated differential clonal evolution within tumors and putative selection against neoantigenic mutations on-therapy. Transcriptome analyses before and during nivolumab therapy revealed increases in distinct immune cell subsets, activation of specific transcriptional networks, and upregulation of immune checkpoint genes that were more pronounced in patients with response. Temporal changes in intratumoral TCR repertoire revealed expansion of T cell clones in the setting of neoantigen loss. Comprehensive genomic profiling data in this study provide insight into nivolumab's mechanism of action.

A Small RNA-Catalytic Argonaute Pathway Tunes Germline Transcript Levels to Ensure Embryonic Divisions.

Multiple division cycles without growth are a characteristic feature of early embryogenesis. The female germline loads proteins and RNAs into oocytes to support these divisions, which lack many quality control mechanisms operating in somatic cells undergoing growth. Here, we describe a small RNA-Argonaute pathway that ensures early embryonic divisions in C. elegans by employing catalytic slicing activity to broadly tune, instead of silence, germline gene expression. Misregulation of one target, a kinesin-13 microtubule depolymerase, underlies a major phenotype associated with pathway loss. Tuning of target transcript levels is guided by the density of homologous small RNAs, whose generation must ultimately be related to target sequence. Thus, the tuning action of a small RNA-catalytic Argonaute pathway generates oocytes capable of supporting embryogenesis. We speculate that the specialized nature of germline chromatin led to the emergence of small RNA-catalytic Argonaute pathways in the female germline as a post-transcriptional control layer to optimize oocyte composition.

Increased homozygosity due to endogamy results in fitness consequences in a human population.

Recessive alleles have been shown to directly affect both human Mendelian disease phenotypes and complex traits. Pedigree studies also suggest that consanguinity results in increased childhood mortality and adverse health phenotypes, presumably through penetrance of recessive mutations. Here, we test whether the accumulation of homozygous, recessive alleles decreases reproductive success in a human population. We address this question among the Namibian Himba, an endogamous agro-pastoralist population, who until very recently practiced natural fertility. Using a sample of 681 individuals, we show that Himba exhibit elevated levels of "inbreeding," calculated as the fraction of the genome in runs of homozygosity (FROH). Many individuals contain multiple long segments of ROH in their genomes, indicating that their parents had high kinship coefficients. However, we do not find evidence that this is explained by first-cousin consanguinity, despite a reported social preference for cross-cousin marriages. Rather, we show that elevated haplotype sharing in the Himba is due to a bottleneck, likely in the past 60 generations. We test whether increased recessive mutation load results in observed fitness consequences by assessing the effect of FROH on completed fertility in a cohort of postreproductive women (n = 69). We find that higher FROH is significantly associated with lower fertility. Our data suggest a multilocus genetic effect on fitness driven by the expression of deleterious recessive alleles, especially those in long ROH. However, these effects are not the result of consanguinity but rather elevated background identity by descent.

Syllables and their beginnings have a special role in the mental lexicon.

The beginnings of words are, in some informal sense, special. This intuition is widely shared, for example, when playing word games. Less apparent is whether the intuition is substantiated empirically and what the underlying organizational principle(s) might be. Here, we answer this seemingly simple question in a quantitatively clear way. Based on arguments about the interplay between lexical storage and speech processing, we examine whether the distribution of information among different speech sounds of words is governed by a critical computational unit for online speech perception and production: syllables. By analyzing lexical databases of twelve languages, we demonstrate that there is a compelling asymmetry between syllable beginnings (onsets) versus ends (codas) in their involvement in distinguishing words stored in the lexicon. In particular, we show that the functional advantage of syllable onset reflects an asymmetrical distribution of lexical informativeness within the syllable unit but not an effect of a global decay of informativeness from the beginning to the end of a word. The converging finding across languages from a range of typological families supports the conjecture that the syllable unit, while being a critical primitive for both speech perception and production, is also a key organizational constraint for lexical storage.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake.

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Genome sequences and population genomics provide insights into the demographic history, inbreeding, and mutation load of two 'living fossil' tree species of Dipteronia.

'Living fossils', that is, ancient lineages of low taxonomic diversity, represent an exceptional evolutionary heritage, yet we know little about how demographic history and deleterious mutation load have affected their long-term survival and extinction risk. We performed whole-genome sequencing and population genomic analyses on Dipteronia sinensis and D. dyeriana, two East Asian Tertiary relict trees. We found large-scale genome reorganizations and identified species-specific genes under positive selection that are likely involved in adaptation. Our demographic analyses suggest that the wider-ranged D. sinensis repeatedly recovered from population bottlenecks over late Tertiary/Quaternary periods of adverse climate conditions, while the population size of the narrow-ranged D. dyeriana steadily decreased since the late Miocene, especially after the Last Glacial Maximum (LGM). We conclude that the efficient purging of deleterious mutations in D. sinensis facilitated its survival and repeated demographic recovery. By contrast, in D. dyeriana, increased genetic drift and reduced selection efficacy, due to recent severe population bottlenecks and a likely preponderance of vegetative propagation, resulted in fixation of strongly deleterious mutations, reduced fitness, and continuous population decline, with likely detrimental consequences for the species' future viability and adaptive potential. Overall, our findings highlight the significant impact of demographic history on levels of accumulation and purging of putatively deleterious mutations that likely determine the long-term survival and extinction risk of Tertiary relict trees.

Genomic signatures of inbreeding and mutation load in tree ferns.

Ferns (Pteridophyta), as the second largest group of vascular plants, play important roles in ecosystem functioning. Homosporous ferns exhibit a remarkable range of mating systems, from extreme inbreeding to obligate outcrossing, which may have significant evolutionary and ecological implications. Despite their significance, the impact of genome-wide inbreeding on genetic diversity and mutation load within the fern lineage remain largely unexplored. In this study, we utilized whole-genome sequencing to investigate the genomic signatures of inbreeding and genetic load in three Alsophila tree fern species. Our analysis revealed extremely high inbreeding in A. spinulosa, in contrast to the predominantly outcrossing observed in A. costularis and A. latebrosa. This difference likely reflects divergent mating systems and demographic histories. Consistent with its extreme inbreeding propensity, A. spinulosa exhibits reduced genetic diversity and a pronounced decline in effective population size. Comparison of genetic load revealed an overall reduction in deleterious mutations in the highly inbred A. spinulosa, highlighting that long-term inbreeding may have contributed to the purging of strongly deleterious mutations, thereby prolonging the survival of A. spinulosa. Despite this, however, A. spinulosa carries a substantive realized genetic load that may potentially instigate future fitness decline. Our findings illuminate the complex evolutionary interplay between inbreeding and mutation load in homosporous ferns, yielding insights with important implications for the conservation and management of these species.

Whole Genomes Inform Genetic Rescue Strategy for Montane Red Foxes in North America.

A few iconic examples have proven the value of facilitated gene flow for counteracting inbreeding depression and staving off extinction; yet, the practice is often not implemented for fear of causing outbreeding depression. Using genomic sequencing, climatic niche modeling, and demographic reconstruction, we sought to assess the risks and benefits of using translocations as a tool for recovery of endangered montane red fox (Vulpes vulpes) populations in the western United States. We demonstrated elevated inbreeding and homozygosity of deleterious alleles across all populations, but especially those isolated in the Cascade and Sierra Nevada ranges. Consequently, translocations would be expected to increase population growth by masking deleterious recessive alleles. Demographic reconstructions further indicated shallow divergences of less than a few thousand years among montane populations, suggesting low risk of outbreeding depression. These genomic-guided findings set the stage for future management, the documentation of which will provide a roadmap for recovery of other data-deficient taxa.

Genomic Reconstruction of the Successful Establishment of a Feralized Bovine Population on the Subantarctic Island of Amsterdam.

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19th century and revealed contributions from European taurine and Indian Ocean Zebu in the founder ancestry. Comparative analysis of whole-genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. Unlike some populations that have experienced prolonged reductions in effective population size, we did not observe any significant purging of highly deleterious variants. Interestingly, the population's success in the harsh environment can be attributed to preadaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. The Amsterdam Island cattle offers valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It also sheds light on the unique genetic legacies of feral populations, raising ethical questions according to conservation efforts.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa.

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

Measuring the Efficiency of Purging by non-random Mating in Human Populations.

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.

A joint Bayesian hierarchical model for estimating SARS-CoV-2 genomic and subgenomic RNA viral dynamics and seroconversion.

Understanding the viral dynamics of and natural immunity to the severe acute respiratory syndrome coronavirus 2 is crucial for devising better therapeutic and prevention strategies for coronavirus disease 2019 (COVID-19). Here, we present a Bayesian hierarchical model that jointly estimates the genomic RNA viral load, the subgenomic RNA (sgRNA) viral load (correlated to active viral replication), and the rate and timing of seroconversion (correlated to presence of antibodies). Our proposed method accounts for the dynamical relationship and correlation structure between the two types of viral load, allows for borrowing of information between viral load and antibody data, and identifies potential correlates of viral load characteristics and propensity for seroconversion. We demonstrate the features of the joint model through application to the COVID-19 post-exposure prophylaxis study and conduct a cross-validation exercise to illustrate the model's ability to impute the sgRNA viral trajectories for people who only had genomic RNA viral load data.

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression.

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Phonological properties of logographic words modulate brain activation in bilinguals: a comparative study of Chinese characters and Japanese Kanji.

The brain networks for the first (L1) and second (L2) languages are dynamically formed in the bilingual brain. This study delves into the neural mechanisms associated with logographic-logographic bilingualism, where both languages employ visually complex and conceptually rich logographic scripts. Using functional Magnetic Resonance Imaging, we examined the brain activity of Chinese-Japanese bilinguals and Japanese-Chinese bilinguals as they engaged in rhyming tasks with Chinese characters and Japanese Kanji. Results showed that Japanese-Chinese bilinguals processed both languages using common brain areas, demonstrating an assimilation pattern, whereas Chinese-Japanese bilinguals recruited additional neural regions in the left lateral prefrontal cortex for processing Japanese Kanji, reflecting their accommodation to the higher phonological complexity of L2. In addition, Japanese speakers relied more on the phonological processing route, while Chinese speakers favored visual form analysis for both languages, indicating differing neural strategy preferences between the 2 bilingual groups. Moreover, multivariate pattern analysis demonstrated that, despite the considerable neural overlap, each bilingual group formed distinguishable neural representations for each language. These findings highlight the brain's capacity for neural adaptability and specificity when processing complex logographic languages, enriching our understanding of the neural underpinnings supporting bilingual language processing.

Increased copy-number variant load of associated risk genes in sporadic cases of amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is an age-related neurodegenerative disease characterized by selective loss of motor neurons in the brainstem and spinal cord. Several genetic factors have been associated to ALS, ranging from causal genes and potential risk factors to disease modifiers. The search for pathogenic variants in these genes has mostly focused on single nucleotide variants (SNVs) while relatively understudied and not fully elucidated is the contribution of structural variants, such as copy number variations (CNVs). Here, we applied an exon-centric aCGH method to investigate, in sporadic ALS patients, the load of CNVs in 131 genes previously associated to ALS. Our approach revealed that CNV load, defined as the total number of CNVs or their size, was significantly higher in ALS cases than controls. About 87% of patients harbored multiple CNVs in ALS-related genes, and 75% structural variants compromised genes directly implicated in ALS pathogenesis (C9orf72, CHCHD10, EPHA4, FUS, HNRNPA1, KIF5A, NEK1, OPTN, PFN1, SOD1, TARDBP, TBK1, UBQLN2, UNC13A, VAPB, VCP). CNV load was also associated to higher onset age and disease progression rate. Although the contribution of individual CNVs in ALS is still unknown, their extensive load in disease-related genes may have relevant implications for the diagnostic, prognostic and therapeutical management of this devastating disorder.

Measuring infectious SARS-CoV-2 in clinical samples reveals a higher viral titer:RNA ratio for Delta and Epsilon vs. Alpha variants.

Novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants pose a challenge to controlling the COVID-19 pandemic. Previous studies indicate that clinical samples collected from individuals infected with the Delta variant may contain higher levels of RNA than previous variants, but the relationship between levels of viral RNA and infectious virus for individual variants is unknown. We measured infectious viral titer (using a microfocus-forming assay) and total and subgenomic viral RNA levels (using RT-PCR) in a set of 162 clinical samples containing SARS-CoV-2 Alpha, Delta, and Epsilon variants that were collected in identical swab kits from outpatient test sites and processed soon after collection. We observed a high degree of variation in the relationship between viral titers and RNA levels. Despite this, the overall infectivity differed among the three variants. Both Delta and Epsilon had significantly higher infectivity than Alpha, as measured by the number of infectious units per quantity of viral E gene RNA (5.9- and 3.0-fold increase; P < 0.0001, P = 0.014, respectively) or subgenomic E RNA (14.3- and 6.9-fold increase; P < 0.0001, P = 0.004, respectively). In addition to higher viral RNA levels reported for the Delta variant, the infectivity (amount of replication competent virus per viral genome copy) may be increased compared to Alpha. Measuring the relationship between live virus and viral RNA is an important step in assessing the infectivity of novel SARS-CoV-2 variants. An increase in the infectivity for Delta may further explain increased spread, suggesting a need for increased measures to prevent viral transmission.

Purging of deleterious burden in the endangered Iberian lynx.

SignificanceThe dynamics of deleterious variation under contrasting demographic scenarios remain poorly understood in spite of their relevance in evolutionary and conservation terms. Here we apply a genomic approach to study differences in the burden of deleterious alleles between the endangered Iberian lynx (Lynx pardinus) and the widespread Eurasian lynx (Lynx lynx). Our analysis unveils a significantly lower deleterious burden in the former species that should be ascribed to genetic purging, that is, to the increased opportunities of selection against recessive homozygotes due to the inbreeding caused by its smaller population size, as illustrated by our analytical predictions. This research provides theoretical and empirical evidence on the evolutionary relevance of genetic purging under certain demographic conditions.