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.

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.

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.

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.

Combination of mutations in genes controlling DNA repair and high mutational load plays a prognostic role in pancreatic ductal adenocarcinoma (PDAC): a retrospective real-life study in Sardinian population.

BACKGROUND: Patients with pancreatic ductal adenocarcinoma (PDCA) carrying impaired mismatch repair mechanisms seem to have an outcome advantage under treatment with conventional chemotherapy, whereas the role for the tumor mutation burden on prognosis is controversial. In this study, we evaluated the prognostic role of the mutated genes involved in genome damage repair in a real-life series of PDAC patients in a hospital-based manner from the main Institution deputed to surgically treat such a disease in North Sardinia. METHODS: A cohort of fifty-five consecutive PDAC patients with potentially resectable/border line resectable PDAC (stage IIB-III) or oligometastatic disease (stage IV) and tumor tissue availability underwent next-generation sequencing (NGS)-based analysis using a panel containing driver oncogenes and tumor suppressor genes as well as genes controlling DNA repair mechanisms. RESULTS: Genes involved in the both genome damage repair (DR) and DNA mismatch repair (MMR) were found mutated in 17 (31%) and 15 (27%) cases, respectively. One fourth of PDAC cases (14/55; 25.5%) carried tumors presenting a combination of mutations in repair genes (DR and MMR) and the highest mutation load rates (MLR-H). After correction for confounders (surgery, adjuvant therapy, stage T, and metastasis), multivariate Cox regression analysis indicated that mutations in DR genes (HR = 3.0126, 95% CI 1.0707 to 8.4764, p = 0.0367) and the MLR (HR = 1.0018, 95%CI 1.0005 to 1.0032, p = 0.009) were significantly related to worse survival. CONCLUSIONS: The combination of mutated repair genes and MLR-H, which is associated with a worse survival in our series of PDAC patients treated with conventional chemotherapy protocols, might become a predictive biomarker of response to immunotherapy in addition to its prognostic role in predicting survival.

Construction and validation of a bladder cancer risk model based on autophagy-related genes.

Autophagy has an important association with tumorigenesis, progression, and prognosis. However, the mechanism of autophagy-regulated genes on the risk prognosis of bladder cancer (BC) patients has not been fully elucidated yet. In this study, we created a prognostic model of BC risk based on autophagy-related genes, which further illustrates the value of genes associated with autophagy in the treatment of BC. We first downloaded human autophagy-associated genes and BC datasets from Human Autophagy Database and The Cancer Genome Atlas (TCGA) database, and finally obtained differential prognosis-associated genes for autophagy by univariate regression analysis and differential analysis of cancer versus normal tissues. Subsequently, we downloaded two datasets from Gene Expression Omnibus (GEO), GSE31684 and GSE15307, to expand the total number of samples. Based on these genes, we distinguished the molecular subtypes (C1, C2) and gene classes (A, B) of BC by consistent clustering analysis. Using the genes merged from TCGA and the two GEO datasets, we conducted least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis to obtain risk genes and construct autophagy-related risk prediction models. The accuracy of this risk prediction model was assessed by receiver operating characteristic (ROC) and calibration curves, and then nomograms were constructed to predict the survival of bladder cancer patients at 1, 3, and 5 years, respectively. According to the median value of the risk score, we divided BC samples into the high- and low-risk groups. Kaplan-Meier (K-M) survival analysis was performed to compare survival differences between subgroups. Then, we used single sample gene set enrichment analysis (ssGSEA) for immune cell infiltration abundance, immune checkpoint genes, immunotherapy response, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, and tumor mutation burden (TMB) analysis for different subgroups. We also applied quantitative real-time polymerase chain reaction (PCR) and immunohistochemistry (IHC) techniques to verify the expression of these six genes in the model. Finally, we chose the IMvigor210 dataset for external validation. Six risk genes associated with autophagy (SPOCD1, FKBP10, NAT8B, LDLR, STMN3, and ANXA2) were finally screened by LASSO regression algorithm and multivariate Cox regression analysis. ROC and calibration curves showed that the model established was accurate and reliable. Univariate and multivariate regression analyses were used to verify that the risk model was an independent predictor. K-M survival analysis indicated that patients in the high-risk group had significantly worse overall survival than those in the low-risk group. Analysis by algorithms such as correlation analysis, gene set variation analysis (GSVA), and ssGSEA showed that differences in immune microenvironment, enrichment of multiple biologically active pathways, TMB, immune checkpoint genes, and human leukocyte antigens (HLAs) were observed in the different risk groups. Then, we constructed nomograms that predicted the 1-, 3-, and 5-year survival rates of different BC patients. In addition, we screened nine sensitive chemotherapeutic drugs using the correlation between the obtained expression status of risk genes and drug sensitivity results. Finally, the external dataset IMvigor210 verified that the model is reliable and efficient. We established an autophagy-related risk prognostic model that is accurate and reliable, which lays the foundation for future personalized treatment of bladder cancer.

Parental care results in a greater mutation load, for which it is also a phenotypic antidote.

Benevolent social behaviours, such as parental care, are thought to enable mildly deleterious mutations to persist. We tested this prediction experimentally using the burying beetle Nicrophorus vespilloides, an insect with biparental care. For 20 generations, we allowed replicate experimental burying beetle populations to evolve either with post-hatching care ('Full Care' populations) or without it ('No Care' populations). We then established new lineages, seeded from these experimental populations, which we inbred to assess their mutation load. Outbred lineages served as controls. We also tested whether the deleterious effects of a greater mutation load could be concealed by parental care by allowing half the lineages to receive post-hatching care, while half did not. We found that inbred lineages from the Full Care populations went extinct more quickly than inbred lineages from the No Care populations-but only when offspring received no post-hatching care. We infer that Full Care lineages carried a greater mutation load, but that the associated deleterious effects on fitness could be overcome if larvae received parental care. We suggest that the increased mutation load caused by parental care increases a population's dependence upon care. This could explain why care is seldom lost once it has evolved.

Sex differences in deleterious genetic variants in a haplodiploid social insect.

Deleterious variants are selected against but can linger in populations at low frequencies for long periods of time, decreasing fitness and contributing to disease burden in humans and other species. Deleterious variants occur at low frequency but distinguishing deleterious variants from low-frequency neutral variation is challenging based on population genomics data alone. As a result, we have little sense of the number and identity of deleterious variants in wild populations. For haplodiploid species, it has been hypothesised that deleterious alleles will be directly exposed to selection in haploid males, but selection can be masked in diploid females when deleterious variants are recessive, resulting in more efficient purging of deleterious mutations in males. Therefore, comparisons of the differences between haploid and diploid genomes from the same population may be a useful method for inferring rare deleterious variants. This study provides the first formal test of this hypothesis. Using wild populations of Northern paper wasps (Polistes fuscatus), we find that males have fewer missense and nonsense variants per generation than females from the same population. Allele frequency differences are especially pronounced for rare missense and nonsense variants and these differences lead to a lower mutational load in males than females. Based on these data we infer that many highly deleterious mutations are segregating in the paper wasp population. Stronger selection against deleterious alleles in haploid males may have implications for adaptation in other haplodiploid insects and provides evidence that wild populations harbour abundant deleterious variants.

Surveillance of cfDNA Hot Spot Mutations in NSCLC Patients during Disease Progression.

Non-small cell cancer (NSCLC) has been identified with a great variation of mutations that can be surveyed during disease progression. The aim of the study was to identify and monitor lung cancer-specific mutations incidence in cell-free DNA as well as overall plasma cell-free DNA load by means of targeted next-generation sequencing. Sequencing libraries were prepared from cell-free DNA (cfDNA) isolated from 72 plasma samples of 41 patients using the Oncomine Lung cfDNA panel covering hot spot regions of 11 genes. Sequencing was performed with the Ion Torrent™ Ion S5™ system. Four genes were detected with highest mutation incidence: KRAS (43.9% of all cases), followed by ALK (36.6%), TP53 (31.7%), and PIK3CA (29.3%). Seven patients had co-occurring KRAS + TP53 (6/41, 14.6%) or KRAS + PIK3CA (7/41, 17.1%) mutations. Moreover, the mutational status of TP53 as well an overall cell-free DNA load were confirmed to be predictors of poor progression-free survival (HR = 2.5 [0.8-7.7]; p = 0.029 and HR = 2.3 [0.9-5.5]; p = 0.029, respectively) in NSCLC patients. In addition, TP53 mutation status significantly predicts shorter overall survival (HR = 3.4 [1.2-9.7]; p < 0.001). We demonstrated that TP53 mutation incidence as well as a cell-free DNA load can be used as biomarkers for NSCLC monitoring and can help to detect the disease progression prior to radiological confirmation of the status.

Mistranslation Reduces Mutation Load in Evolving Proteins through Negative Epistasis with DNA Mutations.

Translational errors during protein synthesis cause phenotypic mutations that are several orders of magnitude more frequent than DNA mutations. Such phenotypic mutations may affect adaptive evolution through their interactions with DNA mutations. To study how mistranslation may affect the adaptive evolution of evolving proteins, we evolved populations of green fluorescent protein (GFP) in either high-mistranslation or low-mistranslation Escherichia coli hosts. In both hosts, we first evolved GFP under purifying selection for the ancestral phenotype green fluorescence, and then under directional selection toward the new phenotype yellow fluorescence. High-mistranslation populations evolved modestly higher yellow fluorescence during each generation of evolution than low-mistranslation populations. We demonstrate by high-throughput sequencing that elevated mistranslation reduced the accumulation of deleterious DNA mutations under both purifying and directional selection. It did so by amplifying the fitness effects of deleterious DNA mutations through negative epistasis with phenotypic mutations. In contrast, mistranslation did not affect the incidence of beneficial mutations. Our findings show that phenotypic mutations interact epistatically with DNA mutations. By reducing a population's mutation load, mistranslation can affect an important determinant of evolvability.

Severe inbreeding, increased mutation load and gene loss-of-function in the critically endangered Devils Hole pupfish.

Small populations with limited range are often threatened by inbreeding and reduced genetic diversity, which can reduce fitness and exacerbate population decline. One of the most extreme natural examples is the Devils Hole pupfish (Cyprinodon diabolis), an iconic and critically endangered species with the smallest known range of any vertebrate. This species has experienced severe declines in population size over the last 30 years and suffered major bottlenecks in 2007 and 2013, when the population shrunk to 38 and 35 individuals, respectively. Here, we analysed 30 resequenced genomes of desert pupfishes from Death Valley, Ash Meadows and surrounding areas to examine the genomic consequences of small population size. We found extremely high levels of inbreeding (FROH = 0.34-0.81) and an increased amount of potentially deleterious genetic variation in the Devils Hole pupfish as compared to other species, including unique, fixed loss-of-function alleles and deletions in genes associated with sperm motility and hypoxia. Additionally, we successfully resequenced a formalin-fixed museum specimen from 1980 and found that the population was already highly inbred prior to recent known bottlenecks. We thus document severe inbreeding and increased mutation load in the Devils Hole pupfish and identify candidate deleterious variants to inform management of this conservation icon.

Not out of the woods yet: Signatures of the prolonged negative genetic consequences of a population bottleneck in a rapidly re-expanding wader, the black-faced spoonbill Platalea minor.

The long-term persistence of a population which has suffered a bottleneck partly depends on how historical demographic dynamics impacted its genetic diversity and the accumulation of deleterious mutations. Here we provide genomic evidence for the genetic effect of a recent population bottleneck in the endangered black-faced spoonbill (Platalea minor) after its rapid population recovery. Our data suggest that the bird's effective population size, Ne , had been relatively stable (7500-9000) since 22,000 years ago; however, a recent brief yet severe bottleneck (Ne  = 20) which we here estimated to occur around the 1940s wiped out >99% of its historical Ne in roughly three generations. Despite a >15-fold population recovery since 1988, we found that black-faced spoonbill population has higher levels of inbreeding (7.4 times more runs of homozygosity) than its sister species, the royal spoonbill (P. regia), which is not thought to have undergone a marked population contraction. Although the two spoonbills have similar levels of genome-wide genetic diversity, our results suggest that selection on more genes was relaxed in the black-faced spoonbill; moreover individual black-faced spoonbills carry more putatively deleterious mutations (Grantham's score > 50), and may therefore express more deleterious phenotypic effects than royal spoonbills. Here we demonstrate the value of using genomic indices to monitor levels of genetic erosion, inbreeding and mutation load in species with conservation concerns. To mitigate the prolonged negative genetic effect of a population bottleneck, we recommend that all possible measures should be employed to maintain population growth of a threatened species.

Hunger for sex: Abundant, heterogeneous resources select for sexual reproduction in the field.

Major hypotheses on sex evolution predict that resource abundance and heterogeneity should either select for or against sexual reproduction. However, seldom have these predictions been explicitly tested in the field. Here, we investigated this question using soil oribatid mites, a diverse and abundant group of soil arthropods whose local communities can be dominated by either sexual or asexual species. First, we refined theoretical predictions by addressing how the effects of resource abundance, heterogeneity and abiotic conditions could modify each other. Then, we estimated the strength of selection for sexual species in local communities while controlling for phylogeny and neutral processes (ecological drift and dispersal), and tested its relation to resource and abiotic gradients. We show that sexual species tended to be favoured with increasing litter amount, a measure of basal resource abundance. Further, there was some evidence that this response occurred mainly under higher tree species richness, a measure of basal resource heterogeneity. This response to resources is unlikely to reflect niche partitioning between reproductive modes, as sexual and asexual species overlapped in trophic niche according to a comparative analysis using literature data on stable isotope ratios. Rather, these findings are consistent with the hypothesis that sex facilitates adaptation by breaking unfavourable genetic associations, an advantage that should increase with effective population size when many loci are under selection and, thus, with resource abundance.

Preparation of Duplex Sequencing Libraries for Archival Paraffin-Embedded Tissue Samples Using Single-Strand-Specific Nuclease P1.

DNA from formalin-fixed paraffin-embedded (FFPE) tissues, which are frequently utilized in cancer research, is significantly affected by chemical degradation. It was suggested that approaches that are based on duplex sequencing can significantly improve the accuracy of mutation detection in FFPE-derived DNA. However, the original duplex sequencing method cannot be utilized for the analysis of formalin-fixed paraffin-embedded (FFPE) tissues, as FFPE DNA contains an excessive number of damaged bases, and these lesions are converted to false double-strand nucleotide substitutions during polymerase-driven DNA end repair process. To resolve this drawback, we replaced DNA polymerase by a single strand-specific nuclease P1. Nuclease P1 was shown to efficiently remove RNA from DNA preparations, to fragment the FFPE-derived DNA and to remove 5'/3'-overhangs. To assess the performance of duplex sequencing-based methods in FFPE-derived DNA, we constructed the Bottleneck Sequencing System (BotSeqS) libraries from five colorectal carcinomas (CRCs) using either DNA polymerase or nuclease P1. As expected, the number of identified mutations was approximately an order of magnitude higher in libraries prepared with DNA polymerase vs. nuclease P1 (626 ± 167/Mb vs. 75 ± 37/Mb, paired t-test p-value 0.003). Furthermore, the use of nuclease P1 but not polymerase-driven DNA end repair allowed a reliable discrimination between CRC tumors with and without hypermutator phenotypes. The utility of newly developed modification was validated in the collection of 17 CRCs and 5 adjacent normal tissues. Nuclease P1 can be recommended for the use in duplex sequencing library preparation from FFPE-derived DNA.

Purifying Selection in Corvids Is Less Efficient on Islands.

Theory predicts that deleterious mutations accumulate more readily in small populations. As a consequence, mutation load is expected to be elevated in species where life-history strategies and geographic or historical contingencies reduce the number of reproducing individuals. Yet, few studies have empirically tested this prediction using genome-wide data in a comparative framework. We collected whole-genome sequencing data for 147 individuals across seven crow species (Corvus spp.). For each species, we estimated the distribution of fitness effects of deleterious mutations and compared it with proxies of the effective population size Ne. Island species with comparatively smaller geographic range sizes had a significantly increased mutation load. These results support the view that small populations have an elevated risk of mutational meltdown, which may contribute to the higher extinction rates observed in island species.

Mutations Beget More Mutations-Rapid Evolution of Mutation Rate in Response to the Risk of Runaway Accumulation.

The rapidity with which the mutation rate evolves could greatly impact evolutionary patterns. Nevertheless, most studies simply assume a constant rate in the time scale of interest (Kimura 1983; Drake 1991; Kumar 2005; Li 2007; Lynch 2010). In contrast, recent studies of somatic mutations suggest that the mutation rate may vary by several orders of magnitude within a lifetime (Kandoth et al. 2013; Lawrence et al. 2013). To resolve the discrepancy, we now propose a runaway model, applicable to both the germline and soma, whereby mutator mutations form a positive-feedback loop. In this loop, any mutator mutation would increase the rate of acquiring the next mutator, thus triggering a runaway escalation in mutation rate. The process can be initiated more readily if there are many weak mutators than a few strong ones. Interestingly, even a small increase in the mutation rate at birth could trigger the runaway process, resulting in unfit progeny. In slowly reproducing species, the need to minimize the risk of this uncontrolled accumulation would thus favor setting the mutation rate low. In comparison, species that starts and ends reproduction sooner do not face the risk and may set the baseline mutation rate higher. The mutation rate would evolve in response to the risk of runaway mutation, in particular, when the generation time changes. A rapidly evolving mutation rate may shed new lights on many evolutionary phenomena (Elango et al. 2006; Thomas et al. 2010, 2018; Langergraber et al. 2012; Besenbacher et al. 2019).

Genomic signatures of inbreeding and mutation load in a threatened rattlesnake.

Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long-term population bottlenecks. Empirical studies that examine this idea using genome-wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding, S. catenatus genomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared to S. tergeminus. At the species level, outbred S. tergeminus had higher genome-wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared to S. catenatus, presumably due to long-term purging of deleterious mutations in S. catenatus. In contrast, mutations that escaped species-level drift effects within S. catenatus populations were in general more frequent and more often found in homozygous genotypes than in S. tergeminus, suggesting a reduced efficiency of purifying selection in smaller S. catenatus populations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.

The roles of hybridization and habitat fragmentation in the evolution of Brazil's enigmatic longwing butterflies, Heliconius nattereri and H. hermathena.

BACKGROUND: Heliconius butterflies are widely distributed across the Neotropics and have evolved a stunning array of wing color patterns that mediate Müllerian mimicry and mating behavior. Their rapid radiation has been strongly influenced by hybridization, which has created new species and allowed sharing of color patterning alleles between mimetic species pairs. While these processes have frequently been observed in widespread species with contiguous distributions, many Heliconius species inhabit patchy or rare habitats that may strongly influence the origin and spread of species and color patterns. Here, we assess the effects of historical population fragmentation and unique biology on the origins, genetic health, and color pattern evolution of two rare and sparsely distributed Brazilian butterflies, Heliconius hermathena and Heliconius nattereri. RESULTS: We assembled genomes and re-sequenced whole genomes of eight H. nattereri and 71 H. hermathena individuals. These species harbor little genetic diversity, skewed site frequency spectra, and high deleterious mutation loads consistent with recent population bottlenecks. Heliconius hermathena consists of discrete, strongly isolated populations that likely arose from a single population that dispersed after the last glacial maximum. Despite having a unique color pattern combination that suggested a hybrid origin, we found no genome-wide evidence that H. hermathena is a hybrid species. However, H. hermathena mimicry evolved via introgression, from co-mimetic Heliconius erato, of a small genomic region upstream of the color patterning gene cortex. CONCLUSIONS: Heliconius hermathena and H. nattereri population fragmentation, potentially driven by historical climate change and recent deforestation, has significantly reduced the genetic health of these rare species. Our results contribute to a growing body of evidence that introgression of color patterning alleles between co-mimetic species appears to be a general feature of Heliconius evolution.