A male-killing Wolbachia endosymbiont is concealed by another endosymbiont and a nuclear suppressor.

Bacteria that live inside the cells of insect hosts (endosymbionts) can alter the reproduction of their hosts, including the killing of male offspring (male killing, MK). MK has only been described in a few insects, but this may reflect challenges in detecting MK rather than its rarity. Here, we identify MK Wolbachia at a low frequency (around 4%) in natural populations of Drosophila pseudotakahashii. MK Wolbachia had a stable density and maternal transmission during laboratory culture, but the MK phenotype which manifested mainly at the larval stage was lost rapidly. MK Wolbachia occurred alongside a second Wolbachia strain expressing a different reproductive manipulation, cytoplasmic incompatibility (CI). A genomic analysis highlighted Wolbachia regions diverged between the 2 strains involving 17 genes, and homologs of the wmk and cif genes implicated in MK and CI were identified in the Wolbachia assembly. Doubly infected males induced CI with uninfected females but not females singly infected with CI-causing Wolbachia. A rapidly spreading dominant nuclear suppressor genetic element affecting MK was identified through backcrossing and subsequent analysis with ddRAD SNPs of the D. pseudotakahashii genome. These findings highlight the complexity of nuclear and microbial components affecting MK endosymbiont detection and dynamics in populations and the challenges of making connections between endosymbionts and the host phenotypes affected by them.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples.

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

A single mutation weakens symbiont-induced reproductive manipulation through reductions in deubiquitylation efficiency.

Animals interact with microbes that affect their performance and fitness, including endosymbionts that reside inside their cells. Maternally transmitted Wolbachia bacteria are the most common known endosymbionts, in large part because of their manipulation of host reproduction. For example, many Wolbachia cause cytoplasmic incompatibility (CI) that reduces host embryonic viability when Wolbachia-modified sperm fertilize uninfected eggs. Operons termed cifs control CI, and a single factor (cifA) rescues it, providing Wolbachia-infected females a fitness advantage. Despite CI's prevalence in nature, theory indicates that natural selection does not act to maintain CI, which varies widely in strength. Here, we investigate the genetic and functional basis of CI-strength variation observed among sister Wolbachia that infect Drosophila melanogaster subgroup hosts. We cloned, Sanger sequenced, and expressed cif repertoires from weak CI-causing wYak in Drosophila yakuba, revealing mutations suspected to weaken CI relative to model wMel in D. melanogaster A single valine-to-leucine mutation within the deubiquitylating (DUB) domain of the wYak cifB homolog (cidB) ablates a CI-like phenotype in yeast. The same mutation reduces both DUB efficiency in vitro and transgenic CI strength in the fly, each by about twofold. Our results map hypomorphic transgenic CI to reduced DUB activity and indicate that deubiquitylation is central to CI induction in cid systems. We also characterize effects of other genetic variation distinguishing wMel-like cifs Importantly, CI strength determines Wolbachia prevalence in natural systems and directly influences the efficacy of Wolbachia biocontrol strategies in transinfected mosquito systems. These approaches rely on strong CI to reduce human disease.

Patterns of Population Structure and Introgression Among Recently Differentiated Drosophila melanogaster Populations.

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

Baseline extracellular vesicle TGF-ß is a predictive biomarker for response to immune checkpoint inhibitors and survival in non-small cell lung cancer.

BACKGROUND: Immune-checkpoint inhibitors (ICIs) are an effective therapeutic strategy, improving the survival of patients with lung cancer compared with conventional treatments. However, novel predictive biomarkers are needed to stratify which patients derive clinical benefit because the currently used and highly heterogenic histological PD-L1 has shown low accuracy. Liquid biopsy is the analysis of biomarkers in body fluids and represents a minimally invasive tool that can be used to monitor tumor evolution and treatment effects, potentially reducing biases associated with tumor heterogeneity associated with tissue biopsies. In this context, cytokines, such as transforming growth factor-ß (TGF-ß), can be found free in circulation in the blood and packaged into extracellular vesicles (EVs), which have a specific delivery tropism and can affect in tumor/immune system interaction. TGF-ß is an immunosuppressive cytokine that plays a crucial role in tumor immune escape, treatment resistance, and metastasis. Thus, we aimed to evaluate the predictive value of circulating and EV TGF-ß in patients with non-small-cell lung cancer receiving ICIs. METHODS: Plasma samples were collected in 33 patients with advanced non-small-cell lung cancer before and during treatment with ICIs. EV were isolated from plasma by serial ultracentrifugation methods and circulating and EV TGF-ß expression levels were evaluated by enzyme-linked immunosorbent assay. RESULTS: Baseline high expression of TGF-ß in EVs was associated with nonresponse to ICIs as well as shorter progression-free survival and overall survival, outperforming circulating TGF-ß levels and tissue PD-L1 as a predictive biomarker. CONCLUSION: If validated, EV TGF-ß could be used to improve patient stratification, increasing the effectiveness of treatment with ICIs and potentially informing combinatory treatments with TGF-ß blockade. PLAIN LANGUAGE SUMMARY: Treatment with immune-checkpoint inhibitors (ICIs) has improved the survival of some patients with lung cancer. However, the majority of patients do not benefit from this treatment, making it essential to develop more reliable biomarkers to identify patients most likely to benefit. In this pilot study, the expression of transforming growth factor-ß (TGF-ß) in blood circulation and in extracellular vesicles was analyzed. The levels of extracellular vesicle TGF-ß before treatment were able to determine which patients would benefit from treatment with ICIs and have a longer survival with higher accuracy than circulating TGF-ß and tissue PD-L1, which is the currently used biomarker in clinical practice.

Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian.

Simulations of laser-induced electron dynamics in a molecular system are performed using time-dependent (TD) equation-of-motion (EOM) coupled-cluster (CC) theory. The target system has been chosen to highlight potential shortcomings of truncated TD-EOM-CC methods [represented in this work by TD-EOM-CC with single and double excitations (TD-EOM-CCSD)], where unphysical spectroscopic features can emerge. Specifically, we explore driven resonant electronic excitations in magnesium fluoride in the proximity of an avoided crossing. Near the avoided crossing, the CCSD similarity-transformed Hamiltonian is defective, meaning that it has complex eigenvalues, and oscillator strengths may take on negative values. When an external field is applied to drive transitions to states exhibiting these traits, unphysical dynamics are observed. For example, the stationary states that make up the time-dependent state acquire populations that can be negative, exceed one, or even complex-valued.

Circulating tumor and invasive cell expression profiling predicts effective therapy in pancreatic cancer.

BACKGROUND: Pancreatic adenocarcinoma (PDAC) remains a refractory disease; however, modern cytotoxic chemotherapeutics can induce tumor regression and extend life. A blood-based, pharmacogenomic, chemosensitivity assay using gene expression profiling of circulating tumor and invasive cells (CTICs) to predict treatment response was previously developed. The combination regimen of 5-fluorouracil, leucovorin, irinotecan, and oxaliplatin (FOLFIRINOX) and gemcitabine/nab-paclitaxel (G/nab-P) are established frontline approaches for treating advanced PDAC; however, there are no validated biomarkers for treatment selection. A similar unmet need exists for choosing second-line therapy. METHODS: The chemosensitivity assay was evaluated in metastatic PDAC patients presenting for frontline treatment. A prospective study enrolled patients (n = 70) before receiving either FOLFIRINOX or G/nab-P at a 1:1 ratio. Six milliliters of peripheral blood was collected at baseline and at time of disease progression. CTICs were isolated, gene-expression profiling was performed, and the assay was used to predict effective and ineffective chemotherapeutic agents. Treating physicians were blinded to the assay prediction results. RESULTS: Patients receiving an effective regimen as predicted by the chemosensitivity assay experienced significantly longer median progression-free survival (mPFS; 7.8 months vs. 4.2 months; hazard ratio [HR], 0.35; p = .0002) and median overall survival (mOS; 21.0 months vs. 9.7 months; HR, 0.40; p = .005), compared with an ineffective regimen. Assay prediction for effective second-line therapy was explored. The entire study cohort experienced favorable outcomes compared with historical controls, 7.1-month mPFS and 12.3-month mOS. CONCLUSIONS: Chemosensitivity assay profiling is a promising tool for guiding therapy in advanced PDAC. Further prospective validation is under way (clinicaltrials.gov NCT03033927).

Mechanochemical Syntheses of Ln(hfac)3(H2O)x (Ln = La-Sm, Tb): Isolation of 10-, 9-, and 8-Coordinate Ln(hfac)n Complexes.

Volatile lanthanide coordination complexes are critical to the generation of new optical and magnetic materials. One of the most common precursors for preparing volatile lanthanide complexes is the hydrate with the general formula Ln(hfac)3(H2O)x (x = 3 for La-Nd, x = 2 for Sm) (hfac = 1,1,1,5,5,5-hexafluoroacetylacetonato). We have investigated the synthesis of Ln(hfac)3(H2O)x using more environmentally sustainable mechanochemical approaches. Characterization of the products using Fourier transform infrared spectroscopy, nuclear magnetic resonance spectroscopy, elemental analysis, and powder X-ray diffraction shows substantial differences in product distribution between methods. The mechanochemical synthesis of the hydrate complexes leads to a variety of coordination compounds including the expected hydrate product, the known retro-Claisen impurity, and hydrated protonated Hhfac ligand depending on the technique employed. Surprisingly, 10-coordinate complexes of the form Na2Ln(hfac)5·3H2O for Ln = La-Nd were also isolated from reactions using a mortar and pestle. The electrostatic bonding of lanthanide coordination complexes is a challenge for obtaining reproducible reactions and clean products. The reproducibility issues are most acute for the large, early lanthanides whereas for the mid to late lanthanides, reproducibility in terms of product distribution and yield is less of an issue because of their smaller size and greater charge to radius ratio. Ball milling increases reproducibility in terms of generating the desired Ln(hfac)3(H2O)x along with hydrated Hhfac (tetraol) and free Hhfac products. The results illustrate the dynamic behavior of lanthanide complexes in solution and the solid state as well as the structural diversity available to the early lanthanides.

Recurrent Collection of Drosophila melanogaster from Wild African Environments and Genomic Insights into Species History.

A long-standing enigma concerns the geographic and ecological origins of the intensively studied vinegar fly, Drosophila melanogaster. This globally distributed human commensal is thought to originate from sub-Saharan Africa, yet until recently, it had never been reported from undisturbed wilderness environments that could reflect its precommensal niche. Here, we document the collection of 288 D. melanogaster individuals from multiple African wilderness areas in Zambia, Zimbabwe, and Namibia. The presence of D. melanogaster in these remote woodland environments is consistent with an ancestral range in southern-central Africa, as opposed to equatorial regions. After sequencing the genomes of 17 wilderness-collected flies collected from Kafue National Park in Zambia, we found reduced genetic diversity relative to town populations, elevated chromosomal inversion frequencies, and strong differences at specific genes including known insecticide targets. Combining these genomes with existing data, we probed the history of this species' geographic expansion. Demographic estimates indicated that expansion from southern-central Africa began ∼10,000 years ago, with a Saharan crossing soon after, but expansion from the Middle East into Europe did not begin until roughly 1,400 years ago. This improved model of demographic history will provide an important resource for future evolutionary and genomic studies of this key model organism. Our findings add context to the history of D. melanogaster, while opening the door for future studies on the biological basis of adaptation to human environments.

A phylogeny for the Drosophila montium species group: A model clade for comparative analyses.

The Drosophila montium species group is a clade of 94 named species, closely related to the model species D. melanogaster. The montium species group is distributed over a broad geographic range throughout Asia, Africa, and Australasia. Species of this group possess a wide range of morphologies, mating behaviors, and endosymbiont associations, making this clade useful for comparative analyses. We use genomic data from 42 available species to estimate the phylogeny and relative divergence times within the montium species group, and its relative divergence time from D. melanogaster. To assess the robustness of our phylogenetic inferences, we use 3 non-overlapping sets of 20 single-copy coding sequences and analyze all 60 genes with both Bayesian and maximum likelihood methods. Our analyses support monophyly of the group. Apart from the uncertain placement of a single species, D. baimaii, our analyses also support the monophyly of all seven subgroups proposed within the montium group. Our phylograms and relative chronograms provide a highly resolved species tree, with discordance restricted to estimates of relatively short branches deep in the tree. In contrast, age estimates for the montium crown group, relative to its divergence from D. melanogaster, depend critically on prior assumptions concerning variation in rates of molecular evolution across branches, and hence have not been reliably determined. We discuss methodological issues that limit phylogenetic resolution - even when complete genome sequences are available - as well as the utility of the current phylogeny for understanding the evolutionary and biogeographic history of this clade.

Pervasive effects of Wolbachia on host activity.

Heritable symbionts have diverse effects on the physiology, reproduction and fitness of their hosts. Maternally transmitted Wolbachia are one of the most common endosymbionts in nature, infecting about half of all insect species. We test the hypothesis that Wolbachia alter host behaviour by assessing the effects of 14 different Wolbachia strains on the locomotor activity of nine Drosophila host species. We find that Wolbachia alter the activity of six different host genotypes, including all hosts in our assay infected with wRi-like Wolbachia strains (wRi, wSuz and wAur), which have rapidly spread among Drosophila species in about the last 14 000 years. While Wolbachia effects on host activity were common, the direction of these effects varied unpredictably and sometimes depended on host sex. We hypothesize that the prominent effects of wRi-like Wolbachia may be explained by patterns of Wolbachia titre and localization within host somatic tissues, particularly in the central nervous system. Our findings support the view that Wolbachia have wide-ranging effects on host behaviour. The fitness consequences of these behavioural modifications are important for understanding the evolution of host-symbiont interactions, including how Wolbachia spread within host populations.

Short Iterative Lanczos Integration in Time-Dependent Equation-of-Motion Coupled-Cluster Theory.

A time-dependent (TD) formulation of equation-of-motion coupled-cluster (EOM-CC) theory can provide excited-state information over an arbitrarily wide energy window with a reduced memory footprint relative to conventional, frequency-domain EOM-CC theory. However, the floating-point costs of the time-integration required by TD-EOM-CC are generally far larger than those of the frequency-domain form of the approach. This work considers the potential of the short iterative Lanczos (SIL) integration scheme [J. Chem. Phys. 1986, 85, 5870-5876] to reduce the floating-point costs of TD-EOM-CC simulations. Low-energy and K-edge absorption features for small molecules are evaluated using TD-EOM-CC with single and double excitations, with the time-integrations carried out via SIL and fourth-order Runge-Kutta (RK4) schemes. Spectra derived from SIL- and RK4-driven simulations are nearly indistinguishable, and with an appropriately chosen subspace dimension, the SIL requires far fewer floating-point operations than are required by RK4. For K-edge spectra, SIL is the more efficient scheme by an average factor of 7.2.

Regulation of cisplatin-resistant head and neck squamous cell carcinoma by the SRC/ETS-1 signaling pathway.

BACKGROUND: We investigated the role of the ETS-1 transcription factor in Head and Neck Squamous Cell Carcinoma (HNSCC) in multiple cisplatin-resistant HNSCC cell lines. METHODS: We examined its molecular link with SRC and MEK/ERK pathways and determined the efficacy of either MEK/ERK inhibitor PD0325901 or SRC inhibitor Dasatinib on cisplatin-resistant HNSCC inhibition. RESULTS: We found that ETS-1 protein expression levels in a majority of cisplatin-resistant HNSCC cell types were higher than those in their parental cisplatin sensitive partners. High ETS-1 expression was also found in patient-derived, cisplatin-resistant HNSCC cells. While ETS-1 knockdown inhibited cell proliferation, migration, and invasion, it could still re-sensitize cells to cisplatin treatment. Interestingly, previous studies have shown that MER/ERK pathways could regulate ETS-1 through its phosphorylation at threonine 38 (T38). Although almost all cisplatin-resistant HNSCC cells we tested showed higher ETS-1 phosphorylation levels at T38, we found that inhibition of MEK/ERK pathways with the MEK inhibitor PD0325901 did not block this phosphorylation. In addition, treatment of cisplatin-resistant HNSCC cells with the MEK inhibitor completely blocked ERK phosphorylation but did not re-sensitize cells to cisplatin treatment. Furthermore, we found that, consistent with ETS-1 increase, SRC phosphorylation dramatically increased in cisplatin-resistant HNSCC, and treatment of cells with the SRC inhibitor, Dasatinib, blocked SRC phosphorylation and decreased ETS-1 expression. Importantly, we showed that Dasatinib, as a single agent, significantly suppressed cell proliferation, migration, and invasion, in addition to survival. CONCLUSIONS: Our results demonstrate that the SRC/ETS-1 pathway plays a crucial role and could be a key therapeutic target in cisplatin-resistant HNSCC treatment.

Revisiting classic clines in Drosophila melanogaster in the age of genomics.

Adaptation to spatially varying environments has been studied for decades, but advances in sequencing technology are now enabling researchers to investigate the landscape of genetic variation underlying this adaptation genome wide. In this review we highlight some of the decades-long research on local adaptation in Drosophila melanogaster from well-studied clines in North America and Australia. We explore the evidence for parallel adaptation and identify commonalities in the genes responding to clinal selection across continents as well as discussing instances where patterns differ among clines. We also investigate recent studies utilizing whole-genome data to identify clines in D. melanogaster and several other systems. Although connecting segregating genomic variation to variation in phenotypes and fitness remains challenging, clinal genomics is poised to increase our understanding of local adaptation and the selective pressures that drive the extensive phenotypic diversity observed in nature.

Synthesis, characterization and antineoplastic activity of bis-aziridinyl dimeric naphthoquinone - A novel class of compounds with potent activity against acute myeloid leukemia cells.

The synthesis, characterization and antileukemic activity of rationally designed amino dimeric naphthoquinone (BiQ) possessing aziridine as alkylating moiety is described. Bis-aziridinyl BiQ decreased proliferation of acute myeloid leukemia (AML) cell lines and primary cells from patients, and exhibited potent (nanomolar) inhibition of colony formation and overall cell survival in AML cells. Effective production of reactive oxygen species (ROS) and double stranded DNA breaks (DSB) induced by bis-aziridinyl BiQ is reported. Bis-dimethylamine BiQ, as the isostere of bis-aziridinyl BiQ but without the alkylating moiety did not show as potent anti-AML activity. Systemic administration of bis-aziridinyl BiQ was well tolerated in NSG mice.

Hyperglycemia and redox status regulate RUNX2 DNA-binding and an angiogenic phenotype in endothelial cells.

Angiogenesis is regulated by hyperglycemic conditions, which can induce cellular stress responses, reactive oxygen species (ROS), and anti-oxidant defenses that modulate intracellular signaling to prevent oxidative damage. The RUNX2 DNA-binding transcription factor is activated by a glucose-mediated intracellular pathway, plays an important role in endothelial cell (EC) function and angiogenesis, and is a target of oxidative stress. RUNX2 DNA-binding and EC differentiation in response to glucose were conserved in ECs from different tissues and inhibited by hyperglycemia, which stimulated ROS production through the aldose reductase glucose-utilization pathway. Furthermore, the redox status of cysteine and methionine residues regulated RUNX2 DNA-binding and reversal of oxidative inhibition was consistent with an endogenous Methionine sulfoxide reductase-A (MsrA) activity. Low molecular weight MsrA substrates and sulfoxide scavengers were potent inhibitors of RUNX2 DNA binding in the absence of oxidative stress, but acted as antioxidants to increase DNA binding in the presence of oxidants. MsrA was associated with RUNX2:DNA complexes, as measured by a sensitive, quantitative DNA-binding ELISA. The related RUNX2 protein family member, RUNX1, which contains an identical DNA-binding domain, was a catalytic substrate of recombinant MsrA. These findings define novel redox pathways involving aldose reductase and MsrA that regulate RUNX2 transcription factor activity and biological function in ECs. Targeting of these pathways could result in more effective strategies to alleviate the vascular dysfunction associated with diabetes or cancer.

Presence of isocitrate dehydrogenase mutations may predict clinical response to hypomethylating agents in patients with acute myeloid leukemia.

Mutations in IDH1 and IDH2 occur in 15-20% of AML cases, resulting in the production of 2-hydroxyglutarate, which promotes aberrant hypermethylation of DNA in leukemic cells. Although these mutations have been shown to have prognostic implications for patients with AML, optimal treatment strategies have yet to be defined. We retrospectively identified forty-two patients with AML treated with DNA methyltransferase inhibitors (DNMTIs) decitabine (n = 36) or azacitidine (n = 6) and performed analysis of stored samples for the presence of IDH1 and IDH2 mutations. Of the forty-two samples analyzed, seven (16.7%) had IDH mutations. Thirteen patients (31%) achieved remission [(complete remission (CR)/complete remission with incomplete count recovery (CRi)/partial response (PR)] after treatment with a DNMTI, five of seven (71.4%) with IDH mutations and eight of thirty-five (22.9%) without IDH mutations (P = 0.01). When adjusted for age at diagnosis, sex, bone marrow blast percentage and cytogenetic, the odds of achieving response after administration of a DNMTI among patients with an IDH mutation was 14.2 when compared to patients without an IDH mutation (95%CI: 1.3-150.4). IDH1 and IDH2 mutations may predict a favorable response to DNMTI in patients with AML.

Describing endosymbiont-host interactions within the parasitism-mutualism continuum.

Endosymbionts are widespread in arthropods, living in host cells with effects that extend from parasitic to mutualistic. Newly acquired endosymbionts tend to be parasitic, but vertical transmission favors coevolution toward mutualism, with hosts sometimes developing dependency. Endosymbionts negatively affecting host fitness may still spread by impacting host reproductive traits, referred to as reproductive "manipulation," although costs for hosts are often assumed rather than demonstrated. For cytoplasmic incompatibility (CI) that involves endosymbiont-mediated embryo death, theory predicts directional shifts away from "manipulation" toward reduced CI strength; moreover, CI-causing endosymbionts need to increase host fitness to initially spread. In nature, endosymbiont-host interactions and dynamics are complex, often depending on environmental conditions and evolutionary history. We advocate for capturing this complexity through appropriate datasets, rather than relying on terms like "manipulation." Such imprecision can lead to the misclassification of endosymbionts along the parasitism-mutualism continuum.

Aedes albopictus is present in the lowlands of southern Zambia.

Identifying the current geographic range of disease vectors is a critical first step towards determining effective mechanisms for controlling and potentially eradicating them. This is particularly true given that historical vector ranges may expand due to changing climates and human activity. The Aedes subgenus Stegomyia contains over 100 species, and among them, Ae. aegypti and Ae. albopictus mosquitoes represent the largest concern for public health, spreading dengue, chikungunya, and zika viruses. While Ae. aegypti has been observed in the country of Zambia for decades, Ae. albopictus has not. In 2015 we sampled four urban and three rural areas in Zambia for Aedes species. Using DNA barcoding, we confirmed the presence of immature and adult Ae. albopictus at two sites: Siavonga and Livingstone. These genotypes seem most closely related to specimens previously collected in Mozambique based on mtDNA barcoding. We resampled Siavonga and Livingstone sites in 2019, again observing immature and adult Ae. albopictus at both sites. Relative Ae. albopictus frequencies were similar between sites, with the exception of immature life stages, which were higher in Siavonga than in Livingstone in 2019. While Ae. albopictus frequencies did not vary through time in Livingstone, both immature and adult frequencies increased through time in Siavonga. This report serves to document the presence of Ae. albopictus in Zambia, which will contribute to understanding the potential public health implications of this disease vector in southern Africa.