The origin and evolutionary history of HIV-1 subtype C in Senegal

Background

The classification of HIV-1 strains in subtypes and Circulating Recombinant Forms (CRFs) has helped in tracking the course of the HIV pandemic. In Senegal, which is located at the tip of West Africa, CRF02_AG predominates in the general population and Female Sex Workers (FSWs). In contrast, 40% of Men having Sex with Men (MSM) in Senegal are infected with subtype C. In this study we analyzed the geographical origins and introduction dates of HIV-1 C in Senegal in order to better understand the evolutionary history of this subtype, which predominates today in the MSM population

Methodology/Principal Findings

We used a combination of phylogenetic analyses and a Bayesian coalescent-based approach, to study the phylogenetic relationships in pol of 56 subtype C isolates from Senegal with 3,025 subtype C strains that were sampled worldwide. Our analysis shows a significantly well supported cluster which contains all subtype C strains that circulate among MSM in Senegal. The MSM cluster and other strains from Senegal are widely dispersed among the different subclusters of African HIV-1 C strains, suggesting multiple introductions of subtype C in Senegal from many different southern and east African countries. More detailed analyses show that HIV-1 C strains from MSM are more closely related to those from southern Africa. The estimated date of the MRCA of subtype C in the MSM population in Senegal is estimated to be in the early 80''s.

Conclusions/Significance

Our evolutionary reconstructions suggest that multiple subtype C viruses with a common ancestor originating in the early 1970s entered Senegal. There was only one efficient spread in the MSM population, which most likely resulted from a single introduction, underlining the importance of high-risk behavior in spread of viruses.     

Deep Sequencing of HIV-1 near Full-Length Proviral Genomes Identifies High Rates of BF1 Recombinants Including Two Novel Circulating Recombinant Forms (CRF) 70_BF1 and a Disseminating 71_BF1 among Blood Donors in Pernambuco,Brazil

Background

The findings of frequent circulation of HIV-1 subclade F1 viruses and the scarcity of BF1 recombinant viruses based on pol subgenomic fragment sequencing among blood donors in Pernambuco (PE), Northeast of Brazil, were reported recently. Here, we aimed to determine whether the classification of these strains (n = 26) extends to the whole genome sequences.

Methods

Five overlapping amplicons spanning the HIV near full-length genomes (NFLGs) were PCR amplified from peripheral blood mononuclear cells (PBMCs) of 26 blood donors. The amplicons were molecularly bar-coded, pooled, and sequenced by Illumina paired-end protocol. The prevalence of viral variants containing drug resistant mutations (DRMs) was compared between plasma and PBMCs.

Results

Of the 26 samples studied, 20 NFLGs and 4 partial fragments were de novo assembled into contiguous sequences and successfully subtyped. Two distinct BF1 recombinant profiles designated CRF70_BF1 and CRF71_BF1, with 4 samples in profile I and 11 in profile II were detected and thus constitute two novel recombinant forms circulating in PE. Evidence of dual infections was detected in four patients co-infected with distinct HIV-1 subtypes. According to our estimate, the new CRF71_BF1 accounts for 10% of the HIV-1 circulating strains among blood donors in PE. Discordant data between the plasma and PBMCs-virus were found in 15 of 24 donors. Six of these strains displayed major DRMs only in PBMCs and four of which had detectable DRMs changes at prevalence between 1-20% of the sequenced population.

Conclusions

The high percentage of the new RF71_BF1 and other BF1 recombinants found among blood donors in Pernambuco, coupled with high rates of transmitted DRMs and dual infections confirm the need for effective surveillance to monitor the prevalence and distribution of HIV variants in a variety of settings in Brazil.     

Phylogeographic study reveals the origin and evolutionary history of a Rhododendron species complex in Taiwan

This study infers a single origin and a once-widespread distribution of the Rhododendron pseudochrysanthum species complex in Taiwan based on chloroplast DNA (cpDNA) variation. In total, 124 individuals from five endemic Rhododendron species were used for amplifications of two chloroplast intergenic spacers: trnL-trnF and atpB-rbcL. The haplotype and nucleotide diversities were much lower for the R. pseudochrysanthum complex, comprised of the species R. pseudochrysanthum, R. morii, Rhododendron rubropunctatum, and Rhododendron hyperythrum, than for Rhododendron formosanum. Two measures of pairwise population differentiation, N(ST) and F(ST), consistently revealed mostly non-significant levels of genetic divergence between populations of the R. pseudochrysanthum complex. No genetic difference was found among the four species of the R. pseudochrysanthum complex by analysis of molecular variance (AMOVA), which is concordant with the parsimonious topology of cpDNA haplotypes for the complex. Nested clade analysis (NCA) of the cpDNA haplotypes indicated that restricted gene flow with isolation-by-distance characterized the recolonization pattern of the R. pseudochrysanthum complex. In contrast, the NCA analysis indicated a contiguous range expansion for cpDNA haplotypes of R. formosanum. This research suggests a once-widespread distribution of the R. pseudochrysanthum complex probably via north-to-south colonization of mid-elevations during low-temperature periods of the Pleistocene. Population fragmentation followed the warmer climate which began in the Holocene and resulted in the present-day range contraction into high elevations.     

Origin and epidemiological history of HIV-1 CRF14_BG

Background

CRF14_BG isolates, originally found in Spain, are characterized by CXCR4 tropism and rapid disease progression. This study aimed to identify the origin of CRF14_BG and reconstruct its epidemiological history based on new isolates from Portugal.

Methodology/Principal Findings

C2V3C3 env gene sequences were obtained from 62 samples collected in 1993–1998 from Portuguese HIV-1 patients. Full-length genomic sequences were obtained from three patients. Viral subtypes, diversity, divergence rate and positive selection were investigated by phylogenetic analysis. The molecular structure of the genomes was determined by bootscanning. A relaxed molecular clock model was used to date the origin of CRF14_BG. Geno2pheno was used to predict viral tropism. Subtype B was the most prevalent subtype (45 sequences; 73%) followed by CRF14_BG (8; 13%), G (4; 6%), F1 (2; 3%), C (2; 3%) and CRF02_AG (1; 2%). Three CRF14_BG sequences were derived from 1993 samples. Near full-length genomic sequences were strongly related to the CRF14_BG isolates from Spain. Genetic diversity of the Portuguese isolates was significantly higher than the Spanish isolates (0.044 vs 0.014, P<0.0001). The mean date of origin of the CRF14_BG cluster was estimated to be 1992 (range, 1989 and 1996) based on the subtype G genomic region and 1989 (range, 1984–1993) based on the subtype B genomic region. Most CRF14_BG strains (78.9%) were predicted to be CXCR4. Finally, up to five amino acids were under selective pressure in subtype B V3 loop whereas only one was found in the CRF14_BG cluster.

Conclusions

CRF14_BG emerged in Portugal in the early 1990 s soon after the beginning of the HIV-1 epidemics, spread to Spain in late 1990 s as a consequence of IVDUs migration and then to the rest of Europe. CXCR4 tropism is a general characteristic of this CRF that may have been selected for by escape from neutralizing antibody response.     

The evolution of HIV-1 and the origin of AIDS

The major cause of acquired immune deficiency syndrome (AIDS) is human immunodeficiency virus type 1 (HIV-1). We have been using evolutionary comparisons to trace (i) the origin(s) of HIV-1 and (ii) the origin(s) of AIDS. The closest relatives of HIV-1 are simian immunodeficiency viruses (SIVs) infecting wild-living chimpanzees (Pan troglodytes troglodytes) and gorillas (Gorilla gorilla gorilla) in west central Africa. Phylogenetic analyses have revealed the origins of HIV-1: chimpanzees were the original hosts of this clade of viruses; four lineages of HIV-1 have arisen by independent cross-species transmissions to humans and one or two of those transmissions may have been via gorillas. However, SIVs are primarily monkey viruses: more than 40 species of African monkeys are infected with their own, species-specific, SIV and in at least some host species, the infection seems non-pathogenic. Chimpanzees acquired from monkeys two distinct forms of SIVs that recombined to produce a virus with a unique genome structure. We have found that SIV infection causes CD4⁺ T-cell depletion and increases mortality in wild chimpanzees, and so the origin of AIDS is more ancient than the origin of HIV-1. Tracing the genetic changes that occurred as monkey viruses adapted to infect first chimpanzees and then humans may provide insights into the causes of the pathogenicity of these viruses.     

Analysis of the evolutionary relationships of HIV-1 and SIVcpz sequences using bayesian inference: implications for the origin of HIV-1

The most plausible origin of HIV-1 group M is an SIV lineage currently represented by SIVcpz isolated from the chimpanzee subspecies Pan troglodytes troglodytes. The origin of HIV-1 group O is less clear. Putative recombination between any of the HIV-1 and SIVcpz sequences was tested using bootscanning and Bayesian-scanning plots, as well as a new method using a Bayesian multiple change-point (BMCP) model to infer parental sequences and crossing-over points. We found that in the case of highly divergent sequences, such as HIV-1/SIVcpz, Bayesian scanning and BMCP methods are more appropriate than bootscanning analysis to investigate spatial phylogenetic variation, including estimating the boundaries of the regions with discordant evolutionary relationships and the levels of support of the phylogenetic clusters under study. According to the Bayesian scanning plots and BMCP method, there was strong evidence for discordant phylogenetic clustering throughout the genome: (1) HIV-1 group O clustered with SIVcpzANT/TAN in middle pol, and partial vif/env; (2) SIVcpzGab1 clustered with SIVcpzANT/TAN in 3'pol/vif, and middle env; (3) HIV-1 group O grouped with SIVcpzCamUS and SIVcpzGab1 in p17/p24; (4) HIV-1 group M was more closely related to SIVcpzCamUS in 3'gag/pol and in middle pol, whereas in partial gp120 group M clustered with group O. Conditionally independent phylogenetic analysis inferred by maximum likelihood (ML) and Bayesian methods further confirmed these findings. The discordant phylogenetic relationships between the HIV-1/SIVcpz sequences may have been caused by ancient recombination events, but they are also due, at least in part, to altered rates of evolution between parental SIVcpz lineages.     

High-resolution molecular epidemiology and evolutionary history of HIV-1 subtypes in Albania

Background

HIV-1 epidemic in Western Europe is largely due to subtype B. Little is known about the HIV-1 in Eastern Europe, but a few studies have shown that non-B subtypes are quite common. In Albania, where a recent study estimated a ten-fold increase of AIDS incidence during the last six years, subtype A and B account for 90% of the know infections.

Methodology/Principal Findings

We investigated the demographic history of HIV-1 subtype A and B in Albania by using a statistical framework based on coalescent theory and phylogeography. High-resolution phylogenetic and molecular clock analysis showed a limited introduction to the Balkan country of subtype A during the late 1980s followed by an epidemic outburst in the early 1990s. In contrast, subtype B was apparently introduced multiple times between the mid-1970s and mid-1980s. Both subtypes are growing exponentially, although the HIV-1A epidemic displays a faster growth rate, and a significantly higher basic reproductive number R₀. HIV-1A gene flow occurs primarily from the capital Tirane, in the center of the country, to the periphery, while HIV-1B flow is characterized by a balanced exchange between center and periphery. Finally, we calculated that the actual number of infections in Albania is at least two orders of magnitude higher than previously thought.

Conclusions/Significance

Our analysis demonstrates the power of recently developed computational tools to investigate molecular epidemiology of pathogens, and emphasize the complex factors involved in the establishment of HIV-1 epidemics. We suggest that a significant correlation exists between HIV-1 exponential spread and the socio-political changes occurred during the Balkan wars. The fast growth of a relatively new non-B epidemic in the Balkans may have significant consequences for the evolution of HIV-1 epidemiology in neighboring countries in Eastern and Western Europe.     

Recent Transmission Clustering of HIV-1 C and CRF17_BF Strains Characterized by NNRTI-Related Mutations among Newly Diagnosed Men in Central Italy

Background

Increased evidence of relevant HIV-1 epidemic transmission in European countries is being reported, with an increased circulation of non-B-subtypes. Here, we present two recent HIV-1 non-B transmission clusters characterized by NNRTI-related amino-acidic mutations among newly diagnosed HIV-1 infected men, living in Rome (Central-Italy).

Methods

Pol and V3 sequences were available at the time of diagnosis for all individuals. Maximum-Likelihood and Bayesian phylogenetic-trees with bootstrap and Bayesian-probability supports defined transmission-clusters. HIV-1 drug-resistance and V3-tropism were also evaluated.

Results

Among 534 new HIV-1 non-B cases, diagnosed from 2011 to 2014, in Central-Italy, 35 carried virus gathering in two distinct clusters, including 27 HIV-1 C and 8 CRF17_BF subtypes, respectively. Both clusters were centralized in Rome, and their origin was estimated to have been after 2007. All individuals within both clusters were males and 37.1% of them had been recently-infected. While C-cluster was entirely composed by Italian men-who-have-sex-with-men, with a median-age of 34 years (IQR:30–39), individuals in CRF17_BF-cluster were older, with a median-age of 51 years (IQR:48–59) and almost all reported sexual-contacts with men and women. All carried R5-tropic viruses, with evidence of atypical or resistance amino-acidic mutations related to NNRTI-drugs (K103Q in C-cluster, and K101E+E138K in CRF17_BF-cluster).

Conclusions

These two epidemiological clusters provided evidence of a strong and recent circulation of C and CRF17_BF strains in central Italy, characterized by NNRTI-related mutations among men engaging in high-risk behaviours. These findings underline the role of molecular epidemiology in identifying groups at increased risk of HIV-1 transmission, and in enhancing additional prevention efforts.     

Analysis of introduction history of invasive plants in Brazil reveals patterns of association between biogeographical origin and reason for introduction

Knowledge of historical factors associated with biological invasions in a region can help identify source regions, vectors and pathways more likely to originate potential invaders as well as prioritize resource allocation for selective prevention and early detection strategies. In Brazil, little is known about the introduction history of many invasive plant species, and analyses of historical factors associated with invasions are lacking. To fill this gap, I used a dataset of 117 invasive alien plants across 13 habitats in Brazil to identify potential patterns of continent of origin, reason for introduction, and to test the hypotheses that (i) more Eurasian species are invasive in Brazil than species native from other continents, that (ii) more horticultural species are invasive in Brazil than species introduced for other reasons, and that (iii) continent of origin and reason for introduction are associated. I found that significantly more invasive plant species in Brazil are native to Africa and Asia, were introduced for horticulture and forage, and are part of the families Poaceae, Fabaceae, and Pinaceae. I also found a significant association between continent of origin and reason for introduction, with more invasive species than average being African forage grasses and Asian agroforestry or ornamental plants. In conclusion, the results suggest that the current invasive alien flora of Brazil results from the combination of patterns of recent human migration waves and deliberate species introductions for technological and commercial reasons. These results can help prioritize invasive species and vectors in prevention, early detection, and control strategies.     

Unequal evolutionary rates in the human immunodeficiency virus type 1 (HIV-1) pandemic: the evolutionary rate of HIV-1 slows down when the epidemic rate increases

HIV-1 sequences in intravenous drug user (IDU) networks are highly homogenous even after several years, while this is not observed in most sexual epidemics. To address this disparity, we examined the human immunodeficiency virus type 1 (HIV-1) evolutionary rate on the population level for IDU and heterosexual transmissions. All available HIV-1 env V3 sequences from IDU outbreaks and heterosexual epidemics with known sampling dates were collected from the Los Alamos HIV sequence database. Evolutionary rates were calculated using phylogenetic trees with a t test root optimization of dated samples. The evolutionary rate of HIV-1 subtype A1 was found to be 8.4 times lower in fast spread among IDUs in the former Soviet Union (FSU) than in slow spread among heterosexual individuals in Africa. Mixed epidemics (IDU and heterosexual) showed intermediate evolutionary rates, indicating a combination of fast- and slow-spread patterns. Hence, if transmissions occur repeatedly during the initial stage of host infection, before selective pressures of the immune system have much impact, the rate of HIV-1 evolution on the population level will decrease. Conversely, in slow spread, where HIV-1 evolves under the pressure of the immune system before a donor infects a recipient, the virus evolution at the population level will increase. Epidemiological modeling confirmed that the evolutionary rate of HIV-1 depends on the rate of spread and predicted that the HIV-1 evolutionary rate in a fast-spreading epidemic, e.g., for IDUs in the FSU, will increase as the population becomes saturated with infections and the virus starts to spread to other risk groups.     

Long-term experimental evolution of HIV-1 reveals effects of environment and mutational history

An often-returning question for not only HIV-1, but also other organisms, is how predictable evolutionary paths are. The environment, mutational history, and random processes can all impact the exact evolutionary paths, but to which extent these factors contribute to the evolutionary dynamics of a particular system is an open question. Especially in a virus like HIV-1, with a large mutation rate and large population sizes, evolution is expected to be highly predictable if the impact of environment and history is low, and evolution is not neutral. We investigated the effect of environment and mutational history by analyzing sequences from a long-term evolution experiment, in which HIV-1 was passaged on 2 different cell types in 8 independent evolutionary lines and 8 derived lines, 4 of which involved a switch of the environment. The experiments lasted for 240–300 passages, corresponding to approximately 400–600 generations or almost 3 years. The sequences show signs of extensive parallel evolution—the majority of mutations that are shared between independent lines appear in both cell types, but we also find that both environment and mutational history significantly impact the evolutionary paths. We conclude that HIV-1 evolution is robust to small changes in the environment, similar to a transmission event in the absence of an immune response or drug pressure. We also find that the fitness landscape of HIV-1 is largely smooth, although we find some evidence for both positive and negative epistatic interactions between mutations.

Analysis of the longest evolutionary experiment with HIV-1 to-date reveals continuous viral adaptation over several years. The authors quantify the environment-specific mutations that arise and determine the fraction of mutations that co-occur with significantly different frequencies than expected by chance.     

A multi-gene approach reveals a complex evolutionary history in the Cyanistes species group

Quaternary climatic oscillations have been considered decisive in shaping much of the phylogeographic structure around the Mediterranean Basin. Within this paradigm, peripheral islands are usually considered as the endpoints of the colonization processes. Here, we use nuclear and mitochondrial markers to investigate the phylogeography of the blue tit complex (blue tit Cyanistes caeruleus, Canary blue tit C. teneriffae and azure tit C. cyanus), and assess the role of the Canary Islands for the geographic structuring of genetic variation. The Canary blue tit exhibits strong genetic differentiation within the Canary Islands and, in combination with other related continental species, provides an ideal model in which to examine recent differentiation within a closely related group of continental and oceanic island avian species. We analysed DNA sequences from 51 breeding populations and more than 400 individuals in the blue tit complex. Discrepancies in the nuclear and mitochondrial gene trees provided evidence of a complex evolutionary process around the Mediterranean Basin. Coalescent analyses revealed gene flow between C. caeruleus and C. teneriffae suggesting a dynamic process with multiple phases of colonization and geographic overlapping ranges. Microsatellite data indicated strong genetic differentiation among the Canary Islands and between the Canary archipelago and the close continental areas, indicating limited contemporary gene flow. Diversification of the blue tit complex is estimated to have started during the early Pliocene (≈ 5 Ma), coincident with the end of Messinian salinity crisis. Phylogenetic analyses indicated that the North African blue tit is derived from the Canary blue tits, a pattern is avian 'back colonization' that contrasts with more traditionally held views of islands being sinks rather than sources.     

IL-12 and GM-CSF in DNA/MVA immunizations against HIV-1 CRF12_BF Nef induced T-cell responses with an enhanced magnitude, breadth and quality

In Argentina, the HIV epidemic is characterized by the co-circulation of subtype B and BF recombinant viral variants. Nef is an HIV protein highly variable among subtypes, making it a good tool to study the impact of HIV variability in the vaccine design setting. We have previously reported a specific cellular response against NefBF with low cross-reactivity to NefB in mice. The aim of this work was to analyze whether the co-administration of IL-12 and GM-CSF, using DNA and MVA vaccine vectors, could improve the final cellular response induced. Mice received three DNA priming doses of a plasmid that express NefBF plus DNAs expressing IL-12 and/or GM-CSF. Afterwards, all the groups were boosted with a MVAnefBF dose. The highest increase in the magnitude of the NefBF response, compared to that induced in the control was found in the IL-12 group. Importantly, a response with higher breadth was detected in groups which received IL-12 or GM-CSF, evidenced as an increased frequency of recognition of homologous (BF) and heterologous (B) Nef peptides, as well as a higher number of other Nef peptide pools representing different viral subtypes. However, these improvements were lost when both DNA cytokines were simultaneously administered, as the response was focused against the immunodominant peptide with a detrimental response towards subdominant epitopes. The pattern of cytokines secreted and the specific-T-cell proliferative capacity were improved in IL-12 and IL-12+GM-CSF groups. Importantly IL-12 generated a significant higher T-cell avidity against a B heterologous peptide.This study indicates that the incorporation of DNA expressing IL-12 in DNA/MVA schemes produced the best results in terms of improvements of T-cell-response key properties such as breadth, cross-reactivity and quality (avidity and pattern of cytokines secreted). These relevant results contribute to the design of strategies aimed to induce T-cell responses against HIV antigens with higher quality.     

The evolutionary history of testicular externalization and the origin of the scrotum

This paper re-examines the evolution of the scrotum and testicular descent in the context of the recent phylogeny of mammals. The adaptive significance of testicular descent and scrotality is briefly discussed. We mapped four character states reflecting the position of testes and presence of scrotum onto recent mammalian phylogeny. Our results are interpreted as follows: as to the presence of testicondy in Monotremata and most of Atlantogenata, which represent the basal group of all eutherians, we argue that primary testicondy represents a plesiomorphic condition for Eutheria as well as for all mammals. This is in opposition to the previous hypothesis of Werdelin and Nilsonne that the scrotum may have evolved before the origin of mammals and then repeatedly disappeared in many groups including monotremes. We suggest that the scrotum evolved at least twice during the evolutionary history of mammals, within Marsupialia and Boreoeutheria, and has subsequently been lost by many groups; this trend is especially strong in Laurasiatheria. We suggest that the recent diversity in testicular position within mammals is the result of multiple selection pressures stemming from the need to provide conditions suitable for sperm development and storage, or to protect the male gonads from excessive physical and physiological disturbance.     

The mid-depth method and HIV-1: a practical approach for testing hypotheses of viral epidemic history.

We introduce the mid-depth method, a practical approach for testing hypotheses of demographic history using genealogies reconstructed from sequence data. The relative positions of internal nodes within a genealogy contain information about past population dynamics. We explain how this information can be used to (1) test the null hypothesis of constant population size and (2) estimate the growth rate and current population size of an exponentially growing population. Simulation tests indicate that, as expected, estimates of exponential growth rates are sometimes biased. The mid-depth method is computationally rapid and does not require knowledge of the sample's mutation rate. However, it does assume that the reconstructed genealogy is correct and is therefore best suited to the analysis of variation-rich viral data sets. When applied to HIV-1 sequence data, the mid-depth method provides phylogenetic evidence of different exponential growth rates for subtypes A and B. We posit that this difference in growth rate reflects the different transmission routes and epidemiological histories of the two subtypes.     

Genome mining reveals the evolutionary origin and biosynthetic potential of basidiomycete polyketide synthases

Numerous polyketides are known from bacteria, plants, and fungi. However, only a few have been isolated from basidiomycetes. Large scale genome sequencing projects now help anticipate the capacity of basidiomycetes to synthesize polyketides. In this study, we identified and annotated 111 type I and three type III polyketide synthase (PKS) genes from 35 sequenced basidiomycete genomes. Phylogenetic analysis of PKS genes suggests that all main types of fungal iterative PKS had already evolved before the Ascomycota and Basidiomycota diverged. A comparison of genomic and metabolomic data shows that the number of polyketide genes exceeds the number of known polyketide structures by far. Exploiting these results to design degenerate PCR primers, we amplified and cloned the complete sequence of armB, a PKS gene from the melleolide producer Armillaria mellea. We expect this study will serve as a guide for future genomic mining projects to discover structurally diverse mushroom-derived polyketides.     

Phylogenomic analysis reveals dynamic evolutionary history of the Drosophila heterochromatin protein 1 (HP1) gene family

Heterochromatin is the gene-poor, satellite-rich eukaryotic genome compartment that supports many essential cellular processes. The functional diversity of proteins that bind and often epigenetically define heterochromatic DNA sequence reflects the diverse functions supported by this enigmatic genome compartment. Moreover, heterogeneous signatures of selection at chromosomal proteins often mirror the heterogeneity of evolutionary forces that act on heterochromatic DNA. To identify new such surrogates for dissecting heterochromatin function and evolution, we conducted a comprehensive phylogenomic analysis of the Heterochromatin Protein 1 gene family across 40 million years of Drosophila evolution. Our study expands this gene family from 5 genes to at least 26 genes, including several uncharacterized genes in Drosophila melanogaster. The 21 newly defined HP1s introduce unprecedented structural diversity, lineage-restriction, and germline-biased expression patterns into the HP1 family. We find little evidence of positive selection at these HP1 genes in both population genetic and molecular evolution analyses. Instead, we find that dynamic evolution occurs via prolific gene gains and losses. Despite this dynamic gene turnover, the number of HP1 genes is relatively constant across species. We propose that karyotype evolution drives at least some HP1 gene turnover. For example, the loss of the male germline-restricted HP1E in the obscura group coincides with one episode of dramatic karyotypic evolution, including the gain of a neo-Y in this lineage. This expanded compendium of ovary- and testis-restricted HP1 genes revealed by our study, together with correlated gain/loss dynamics and chromosome fission/fusion events, will guide functional analyses of novel roles supported by germline chromatin.