Evidence for frequency‐dependent selection maintaining polymorphism in the Batesian mimic Papilio polytes in multiple islands in the Ryukyus,Japan

Batesian mimicry is a well‐studied adaptation for predation avoidance, in which a mimetic species resembles an unpalatable model species. Batesian mimicry can be under positive selection because of the protection gained against predators, due to resemblance to unpalatable model species. However, in some mimetic species, nonmimetic individuals are present in populations, despite the benefits of mimicry. The mechanism for evolution of such mimetic polymorphism remains an open question. Here, we address the hypothesis that the abundance of mimics is limited by that of the models, leading to mimetic polymorphism. In addition, other forces such as the effects of common ancestry and/or isolation by distance may explain this phenomenon. To investigate this question, we focused on the butterfly, Papilio polytes, that exhibits mimetic polymorphism on multiple islands of the Ryukyus, Japan, and performed field surveys and genetic analysis. We found that the mimic ratio of P. polytes was strongly correlated with the model abundance observed on each of the five islands, suggesting negative frequency‐dependent selection is driving the evolution of polymorphism in P. polytes populations. Molecular phylogenetic analysis indicated that the southern island populations are the major source of genetic diversity, and the middle and northern island populations arose by relatively recent migration. This view was also supported by mismatch distribution and Tajima's D analyses, suggesting a recent population expansion on the middle and northern islands, and stable population persistence on the southern islands. The frequency of the mimetic forms within P. polytes populations is thus explained by variations in the model abundance rather than by population structure. Thus, we propose that predation pressure, rather than neutral forces, have shaped the Batesian mimicry polymorphism in P. polytes observed in the Ryukyus.     

Low‐coverage reduced representation sequencing reveals subtle within‐island genetic structure in Aldabra giant tortoises

Aldabrachelys gigantea (Aldabra giant tortoise) is one of only two giant tortoise species left in the world and survives as a single wild population of over 100,000 individuals on Aldabra Atoll, Seychelles. Despite this large current population size, the species faces an uncertain future because of its extremely restricted distribution range and high vulnerability to the projected consequences of climate change. Captive‐bred A. gigantea are increasingly used in rewilding programs across the region, where they are introduced to replace extinct giant tortoises in an attempt to functionally resurrect degraded island ecosystems. However, there has been little consideration of the current levels of genetic variation and differentiation within and among the islands on Aldabra. As previous microsatellite studies were inconclusive, we combined low‐coverage and double‐digest restriction‐associated DNA (ddRAD) sequencing to analyze samples from 33 tortoises (11 from each main island). Using 5426 variant sites within the tortoise genome, we detected patterns of within‐island population structure, but no differentiation between the islands. These unexpected results highlight the importance of using genome‐wide genetic markers to capture higher‐resolution genetic structure to inform future management plans, even in a seemingly panmictic population. We show that low‐coverage ddRAD sequencing provides an affordable alternative approach to conservation genomic projects of non‐model species with large genomes.     

On‐site genetic analysis for species identification using lab‐on‐a‐chip

This paper presents a microfluidic device capable of performing genetic analysis on dung samples to identify White Rhinoceros (Ceratotherium simum). The development of a microfluidic device, which can be used in the field, offers a portable and cost‐effective solution for DNA analysis and species identification to aid conservation efforts. Optimization of the DNA extraction processes produced equivalent yields compared to conventional kit‐based methods within just 5 minutes. The use of a color‐changing loop‐mediated isothermal amplification reaction for simultaneous detection of the cytochrome B sequence of C. simum enabled positive results to be obtained within as little as 30 minutes. Field testing was performed at Knowsley Safari to demonstrate real‐world applicability of the microfluidic device for testing of biological samples.     

Concordant geographic and genetic structure revealed by genotyping‐by‐sequencing in a New Zealand marine isopod

Population genetic structure in the marine environment can be influenced by life‐history traits such as developmental mode (biphasic, with distinct adult and larval morphology, and direct development, in which larvae resemble adults) or habitat specificity, as well as geography and selection. Developmental mode is thought to significantly influence dispersal, with direct developers expected to have much lower dispersal potential. However, this prediction can be complicated by the presence of geophysical barriers to dispersal. In this study, we use a panel of 8,020 SNPs to investigate population structure and biogeography over multiple spatial scales for a direct‐developing species, the New Zealand endemic marine isopod Isocladus armatus. Because our sampling range is intersected by two well‐known biogeographic barriers (the East Cape and the Cook Strait), our study provides an opportunity to understand how such barriers influence dispersal in direct developers. On a small spatial scale (20 km), gene flow between locations is extremely high, suggestive of an island model of migration. However, over larger spatial scales (600 km), populations exhibit a clear pattern of isolation‐by‐distance. Our results indicate that I. armatus exhibits significant migration across the hypothesized barriers and suggest that large‐scale ocean currents associated with these locations do not present a barrier to dispersal. Interestingly, we find evidence of a north‐south population genetic break occurring between Māhia and Wellington. While no known geophysical barrier is apparent in this area, it coincides with the location of a proposed border between bioregions. Analysis of loci under selection revealed that both isolation‐by‐distance and adaption may be contributing to the degree of population structure we have observed here. We conclude that developmental life history largely predicts dispersal in the intertidal isopod I. armatus. However, localized biogeographic processes can disrupt this expectation, and this may explain the potential meta‐population detected in the Auckland region.     

Attenuation of SARS‐CoV‐2 replication and associated inflammation by concomitant targeting of viral and host cap 2'‐O‐ribose methyltransferases

The SARS‐CoV‐2 infection cycle is a multistage process that relies on functional interactions between the host and the pathogen. Here, we repurposed antiviral drugs against both viral and host enzymes to pharmaceutically block methylation of the viral RNA 2''‐O‐ribose cap needed for viral immune escape. We find that the host cap 2''‐O‐ribose methyltransferase MTr1 can compensate for loss of viral NSP16 methyltransferase in facilitating virus replication. Concomitant inhibition of MTr1 and NSP16 efficiently suppresses SARS‐CoV‐2 replication. Using in silico target‐based drug screening, we identify a bispecific MTr1/NSP16 inhibitor with anti‐SARS‐CoV‐2 activity in vitro and in vivo but with unfavorable side effects. We further show antiviral activity of inhibitors that target independent stages of the host SAM cycle providing the methyltransferase co‐substrate. In particular, the adenosylhomocysteinase (AHCY) inhibitor DZNep is antiviral in in vitro, in ex vivo, and in a mouse infection model and synergizes with existing COVID‐19 treatments. Moreover, DZNep exhibits a strong immunomodulatory effect curbing infection‐induced hyperinflammation and reduces lung fibrosis markers ex vivo. Thus, multispecific and metabolic MTase inhibitors constitute yet unexplored treatment options against COVID‐19.     

Molecular basis of cross‐species ACE2 interactions with SARS‐CoV‐2‐like viruses of pangolin origin

Pangolins have been suggested as potential reservoir of zoonotic viruses, including SARS‐CoV‐2 causing the global COVID‐19 outbreak. Here, we study the binding of two SARS‐CoV‐2‐like viruses isolated from pangolins, GX/P2V/2017 and GD/1/2019, to human angiotensin‐converting enzyme 2 (hACE2), the receptor of SARS‐CoV‐2. We find that the spike protein receptor‐binding domain (RBD) of pangolin CoVs binds to hACE2 as efficiently as the SARS‐CoV‐2 RBD in vitro. Furthermore, incorporation of pangolin CoV RBDs allows entry of pseudotyped VSV particles into hACE2‐expressing cells. A screen for binding of pangolin CoV RBDs to ACE2 orthologs from various species suggests a broader host range than that of SARS‐CoV‐2. Additionally, cryo‐EM structures of GX/P2V/2017 and GD/1/2019 RBDs in complex with hACE2 show their molecular binding in modes similar to SARS‐CoV‐2 RBD. Introducing the Q498H substitution found in pangolin CoVs into the SARS‐CoV‐2 RBD expands its binding capacity to ACE2 homologs of mouse, rat, and European hedgehog. These findings suggest that these two pangolin CoVs may infect humans, highlighting the necessity of further surveillance of pangolin CoVs.     

Engineering defensin α‐helix to produce high‐affinity SARS‐CoV‐2 spike protein binding ligands

The binding of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) spike protein to the angiotensin‐converting enzyme 2 (ACE2) receptor expressed on the host cells is a critical initial step for viral infection. This interaction is blocked through competitive inhibition by soluble ACE2 protein. Therefore, developing high‐affinity and cost‐effective ACE2 mimetic ligands that disrupt this protein–protein interaction is a promising strategy for viral diagnostics and therapy. We employed human and plant defensins, a class of small (2–5 kDa) and highly stable proteins containing solvent‐exposed alpha‐helix, conformationally constrained by two disulfide bonds. Therefore, we engineered the amino acid residues on the constrained alpha‐helix of defensins to mimic the critical residues on the ACE2 helix 1 that interact with the SARS‐CoV‐2 spike protein. The engineered proteins (h‐deface2, p‐deface2, and p‐deface2‐MUT) were soluble and purified to homogeneity with a high yield from a bacterial expression system. The proteins demonstrated exceptional thermostability (Tm 70.7°C), high‐affinity binding to the spike protein with apparent K _d values of 54.4 ± 11.3, 33.5 ± 8.2, and 14.4 ± 3.5 nM for h‐deface2, p‐deface2, and p‐deface2‐MUT, respectively, and were used in a diagnostic assay that detected SARS‐CoV‐2 neutralizing antibodies. This work addresses the challenge of developing helical ACE2 mimetics by demonstrating that defensins provide promising scaffolds to engineer alpha‐helices in a constrained form for designing of high‐affinity ligands.     

Seed dispersal by wind decreases when plants are water‐stressed,potentially counteracting species coexistence and niche evolution

Hydrology is a major environmental factor determining plant fitness, and hydrological niche segregation (HNS) has been widely used to explain species coexistence. Nevertheless, the distribution of plant species along hydrological gradients does not only depend on their hydrological niches but also depend on their seed dispersal, with dispersal either weakening or reinforcing the effects of HNS on coexistence. However, it is poorly understood how seed dispersal responds to hydrological conditions. To close this gap, we conducted a common‐garden experiment exposing five wind‐dispersed plant species (Bellis perennis, Chenopodium album, Crepis sancta, Hypochaeris glabra, and Hypochaeris radicata) to different hydrological conditions. We quantified the effects of hydrological conditions on seed production and dispersal traits, and simulated seed dispersal distances with a mechanistic dispersal model. We found species‐specific responses of seed production, seed dispersal traits, and predicted dispersal distances to hydrological conditions. Despite these species‐specific responses, there was a general positive relationship between seed production and dispersal distance: Plants growing in favorable hydrological conditions not only produce more seeds but also disperse them over longer distances. This arises mostly because plants growing in favorable environments grow taller and thus disperse their seeds over longer distances. We postulate that the positive relationship between seed production and dispersal may reduce the concentration of each species to the environments favorable for it, thus counteracting species coexistence. Moreover, the resulting asymmetrical gene flow from favorable to stressful habitats may slow down the microevolution of hydrological niches, causing evolutionary niche conservatism. Accounting for context‐dependent seed dispersal should thus improve ecological and evolutionary models for the spatial dynamics of plant populations and communities.     

CellRegMap: a statistical framework for mapping context‐specific regulatory variants using scRNA‐seq

Single‐cell RNA sequencing (scRNA‐seq) enables characterizing the cellular heterogeneity in human tissues. Recent technological advances have enabled the first population‐scale scRNA‐seq studies in hundreds of individuals, allowing to assay genetic effects with single‐cell resolution. However, existing strategies to analyze these data remain based on principles established for the genetic analysis of bulk RNA‐seq. In particular, current methods depend on a priori definitions of discrete cell types, and hence cannot assess allelic effects across subtle cell types and cell states. To address this, we propose the Cell Regulatory Map (CellRegMap), a statistical framework to test for and quantify genetic effects on gene expression in individual cells. CellRegMap provides a principled approach to identify and characterize genotype–context interactions of known eQTL variants using scRNA‐seq data. This model‐based approach resolves allelic effects across cellular contexts of different granularity, including genetic effects specific to cell subtypes and continuous cell transitions. We validate CellRegMap using simulated data and apply it to previously identified eQTL from two recent studies of differentiating iPSCs, where we uncover hundreds of eQTL displaying heterogeneity of genetic effects across cellular contexts. Finally, we identify fine‐grained genetic regulation in neuronal subtypes for eQTL that are colocalized with human disease variants.     

LOXL1‐AS1 communicating with TIAR modulates vasculogenic mimicry in glioma via regulation of the miR‐374b‐5p/MMP14 axis

At present, growing evidence indicates that long non‐coding RNAs (lncRNAs) participate in the progression of glioma. The function of LOXL1‐AS1 in vasculogenic mimicry (VM) in glioma remains unclear. First, the expressions of TIAR, the lncRNA LOXL1‐AS1, miR‐374b‐5p and MMP14 were examined by qRT‐PCR and Western blot in both, glioma tissues and glioma cell lines. Proliferation, migration, invasion and tube formation assays were conducted to evaluate the roles of TIAR, LOXL1‐AS1, miR‐374b‐5p and MMP14 in malignant cellular behaviours in glioma cells. A nude mouse xenograft model and dual staining for CD34 and PAS were used to assess whether VM was affected by TIAR, LOXL1‐AS1 or miR‐374b‐5p in vivo. In this study, low levels of TIAR and high levels of LOXL1‐AS1 were found in glioma cells and tissues. TIAR downregulated the expression of LOXL1‐AS1 by destabilizing it. LOXL1‐AS1 acted like a miRNA sponge towards miR‐374b‐5p so that downregulation of the former greatly inhibited cell proliferation, migration, invasion and VM. Additionally, miR‐374b‐5p overexpression repressed malignant biological behaviours and VM in glioma by modifying MMP14. In summary, we demonstrated that TIAR combined with LOXL1‐AS1 modulates VM in glioma via the miR‐374b‐5p/MMP14 axis, revealing novel targets for glioma therapy.     

Genome‐scale metabolic modeling reveals SARS‐CoV‐2‐induced metabolic changes and antiviral targets

Tremendous progress has been made to control the COVID‐19 pandemic caused by the SARS‐CoV‐2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS‐CoV‐2 infection using genome‐scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS‐CoV‐2 infection. We next applied the GEM‐based metabolic transformation algorithm to predict anti‐SARS‐CoV‐2 targets that counteract the virus‐induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco‐2 cells. Further generating and analyzing RNA‐sequencing data of remdesivir‐treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti‐SARS‐CoV‐2 drug. Our study provides clinical data‐supported candidate anti‐SARS‐CoV‐2 targets for future evaluation, demonstrating host metabolism targeting as a promising antiviral strategy.     

Genome‐wide selective detection of Mile red‐bone goat using next‐generation sequencing technology

The ecotype population of goats (Capra hircus) was created by long‐term artificial selection and natural adaptation. Mile red‐bone goat is an indigenous breed with visible red bones, and its special bone structure has received extensive attention. This study aimed to identify genetic variants and candidate genes associated with specific bone phenotypes using next‐generation sequencing technology (NGS). The results revealed that 31,828,206 single nucleotide polymorphisms (SNPs) were obtained from 72 goats (20 Mile red‐bone goats and 52 common goats) by NGS. A total of 100 candidate genes were identified on the basis top 1% window interaction from nucleotide diversity (π), π ratio (π _A/π _B), and pairwise fixation index (F _ST). Exactly 77 known signaling pathways were enriched. Specifically, three coding genes (NMNAT2, LOC102172983, and PNLIP) were annotated in the vitamin metabolism signaling pathways, and NCF2 was annotated to the osteoclast (OC) differentiation pathway. Furthermore, 5862 reliable copy number variations (CNVs) were obtained, and 14 and 24 genes were annotated with the top 1‰ CNV based on F _ST (>0.490) and V _ST (>0.527), respectively. Several pathways related to bone development and metabolism of exogenous substances in vivo, including calcium signaling pathway, OC differentiation, and glycerophospholipid metabolism, were annotated. Specifically, six genes from 19 candidate CNVs, which were obtained by interaction of the top 1‰ CNVs with F _ST and V _ST, were annotated to mucin‐type O‐glycan biosynthesis and metabolic pathways. Briefly, the results implied that pseudopurpurin and specific genetic variants work together to contribute to the red‐bone color and specific bone structure of Mile red‐bone goat. This study is helpful to understanding the genetic basis of the unique bone phenotype of Mile red‐bone goats.     

How rapidly do self‐compatible populations evolve selfing? Mating system estimation within recently evolved self‐compatible populations of Azorean Tolpis succulenta (Asteraceae)

Genome‐wide genotyping and Bayesian inference method (BORICE) were employed to estimate outcrossing rates and paternity in two small plant populations of Tolpis succulenta (Asteraceae) on Graciosa island in the Azores. These two known extant populations of T. succulenta on Graciosa have recently evolved self‐compatibility. Despite the expectation that selfing would occur at an appreciable rate (self‐incompatible populations of the same species show low but nonzero selfing), high outcrossing was found in progeny arrays from maternal plants in both populations. This is inconsistent with an immediate transition to high selfing following the breakdown of a genetic incompatibility system. This finding is surprising given the small population sizes and the recent colonization of an island from self‐incompatible colonists of T. succulenta from another island in the Azores, and a potential paucity of pollinators, all factors selecting for selfing through reproductive assurance. The self‐compatible lineage(s) likely have high inbreeding depression (ID) that effectively halts the evolution of increased selfing, but this remains to be determined. Like their progeny, all maternal plants in both populations are fully outbred, which is consistent with but not proof of high ID. High multiple paternity was found in both populations, which may be due in part to the abundant pollinators observed during the flowering season.     

Aging‐related cell type‐specific pathophysiologic immune responses that exacerbate disease severity in aged COVID‐19 patients

Coronavirus disease 2019 (COVID‐19) is especially severe in aged patients, defined as 65 years or older, for reasons that are currently unknown. To investigate the underlying basis for this vulnerability, we performed multimodal data analyses on immunity, inflammation, and COVID‐19 incidence and severity as a function of age. Our analysis leveraged age‐specific COVID‐19 mortality and laboratory testing from a large COVID‐19 registry, along with epidemiological data of ~3.4 million individuals, large‐scale deep immune cell profiling data, and single‐cell RNA‐sequencing data from aged COVID‐19 patients across diverse populations. We found that decreased lymphocyte count and elevated inflammatory markers (C‐reactive protein, D‐dimer, and neutrophil–lymphocyte ratio) are significantly associated with age‐specific COVID‐19 severities. We identified the reduced abundance of naïve CD8 T cells with decreased expression of antiviral defense genes (i.e., IFITM3 and TRIM22) in aged severe COVID‐19 patients. Older individuals with severe COVID‐19 displayed type I and II interferon deficiencies, which is correlated with SARS‐CoV‐2 viral load. Elevated expression of SARS‐CoV‐2 entry factors and reduced expression of antiviral defense genes (LY6E and IFNAR1) in the secretory cells are associated with critical COVID‐19 in aged individuals. Mechanistically, we identified strong TGF‐beta‐mediated immune–epithelial cell interactions (i.e., secretory‐non‐resident macrophages) in aged individuals with critical COVID‐19. Taken together, our findings point to immuno‐inflammatory factors that could be targeted therapeutically to reduce morbidity and mortality in aged COVID‐19 patients.     

Male‐biased dispersal in a fungus‐gardening ant symbiosis

For nearly all organisms, dispersal is a fundamental life‐history trait that can shape their ecology and evolution. Variation in dispersal capabilities within a species exists and can influence population genetic structure and ecological interactions. In fungus‐gardening (attine) ants, co‐dispersal of ants and mutualistic fungi is crucial to the success of this obligate symbiosis. Female‐biased dispersal (and gene flow) may be favored in attines because virgin queens carry the responsibility of dispersing the fungi, but a paucity of research has made this conclusion difficult. Here, we investigate dispersal of the fungus‐gardening ant Trachymyrmex septentrionalis using a combination of maternally (mitochondrial DNA) and biparentally inherited (microsatellites) markers. We found three distinct, spatially isolated mitochondrial DNA haplotypes; two were found in the Florida panhandle and the other in the Florida peninsula. In contrast, biparental markers illustrated significant gene flow across this region and minimal spatial structure. The differential patterns uncovered from mitochondrial DNA and microsatellite markers suggest that most long‐distance ant dispersal is male‐biased and that females (and concomitantly the fungus) have more limited dispersal capabilities. Consequently, the limited female dispersal is likely an important bottleneck for the fungal symbiont. This bottleneck could slow fungal genetic diversification, which has significant implications for both ant hosts and fungal symbionts regarding population genetics, species distributions, adaptive responses to environmental change, and coevolutionary patterns.     

MicroRNA‐495‐3p diminishes doxorubicin‐induced cardiotoxicity through activating AKT

Doxorubicin (Dox) is a broad‐spectrum antitumour agent; however, its clinical application is impeded due to the cumulative cardiotoxicity. The present study aims to investigate the role and underlying mechanisms of microRNA‐495‐3p (miR‐495‐3p) in Dox‐induced cardiotoxicity. Herein, we found that cardiac miR‐495‐3p expression was significantly decreased in Dox‐treated hearts, and that the miR‐495‐3p agomir could prevent oxidative stress, cell apoptosis, cardiac mass loss, fibrosis and cardiac dysfunction upon Dox stimulation. In contrast, the miR‐495‐3p antagomir dramatically aggravated Dox‐induced cardiotoxicity in mice. Besides, we found that the miR‐495‐3p agomir attenuated, while the miR‐495‐3p antagomir exacerbated Dox‐induced oxidative stress and cellular injury in vitro. Mechanistically, we demonstrated that miR‐495‐3p directly bound to the 3′‐untranslational region of phosphate and tension homology deleted on chromosome ten (PTEN), downregulated PTEN expression and subsequently activated protein kinase B (PKB/AKT) pathway, and that PTEN overexpression or AKT inhibition completely abolished the cardioprotective effects of the miR‐495‐3p agomir. Our study for the first time identify miR‐495‐3p as an endogenous protectant against Dox‐induced cardiotoxicity through activating AKT pathway in vivo and in vitro.     

A genotyping‐in‐thousands by sequencing panel to inform invasive deer management using noninvasive fecal and hair samples

Studies in ecology, evolution, and conservation often rely on noninvasive samples, making it challenging to generate large amounts of high‐quality genetic data for many elusive and at‐risk species. We developed and optimized a Genotyping‐in‐Thousands by sequencing (GT‐seq) panel using noninvasive samples to inform the management of invasive Sitka black‐tailed deer (Odocoileus hemionus sitkensis) in Haida Gwaii (Canada). We validated our panel using paired high‐quality tissue and noninvasive fecal and hair samples to simultaneously distinguish individuals, identify sex, and reconstruct kinship among deer sampled across the archipelago, then provided a proof‐of‐concept application using field‐collected feces on SGang Gwaay, an island of high ecological and cultural value. Genotyping success across 244 loci was high (90.3%) and comparable to that of high‐quality tissue samples genotyped using restriction‐site associated DNA sequencing (92.4%), while genotyping discordance between paired high‐quality tissue and noninvasive samples was low (0.50%). The panel will be used to inform future invasive species operations in Haida Gwaii by providing individual and population information to inform management. More broadly, our GT‐seq workflow that includes quality control analyses for targeted SNP selection and a modified protocol may be of wider utility for other studies and systems where noninvasive genetic sampling is employed.     

Mito‐nuclear discordance at a mimicry color transition zone in bumble bee Bombus melanopygus

As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee, Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern “Rocky Mountain''’ color form with ferruginous mid‐abdominal segments (B. m. melanopygus) and a southern “Pacific''’ form with black mid‐abdominal segments (B. m. edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data for B. melanopygus from the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color‐controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito‐nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.     

SARS‐CoV‐2 infection remodels the host protein thermal stability landscape

The severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) is a global threat to human health and has compromised economic stability. In addition to the development of an effective vaccine, it is imperative to understand how SARS‐CoV‐2 hijacks host cellular machineries on a system‐wide scale so that potential host‐directed therapies can be developed. In situ proteome‐wide abundance and thermal stability measurements using thermal proteome profiling (TPP) can inform on global changes in protein activity. Here we adapted TPP to high biosafety conditions amenable to SARS‐CoV‐2 handling. We discovered pronounced temporal alterations in host protein thermostability during infection, which converged on cellular processes including cell cycle, microtubule and RNA splicing regulation. Pharmacological inhibition of host proteins displaying altered thermal stability or abundance during infection suppressed SARS‐CoV‐2 replication. Overall, this work serves as a framework for expanding TPP workflows to globally important human pathogens that require high biosafety containment and provides deeper resolution into the molecular changes induced by SARS‐CoV‐2 infection.