The Incompetence of Mosquitoes-Can Zika Virus Be Adapted To Infect Culex tarsalis Cells?

The molecular evolutionary mechanisms underpinning virus-host interactions are increasingly recognized as key drivers of virus emergence, host specificity, and the likelihood that viruses can undergo a host shift that alters epidemiology and transmission biology. Zika virus (ZIKV) is mainly transmitted between humans by Aedes aegypti mosquitoes. However, the 2015 to 2017 outbreak stimulated discussion regarding the role of Culex spp. mosquitoes in transmission. Reports of ZIKV-infected Culex mosquitoes, in nature and under laboratory conditions, resulted in public and scientific confusion. We previously found that Puerto Rican ZIKV does not infect colonized Culex quinquefasciatus, Culex pipiens, or Culex tarsalis, but some studies suggest they may be competent ZIKV vectors. Therefore, we attempted to adapt ZIKV to Cx. tarsalis by serially passaging virus on cocultured Ae. aegypti (Aag2) and Cx. tarsalis (CT) cells to identify viral determinants of species specificity. Increasing fractions of CT cells resulted in decreased overall virus titer and no enhancement of Culex cell or mosquito infection. Next-generation sequencing of cocultured virus passages revealed synonymous and nonsynonymous variants throughout the genome that arose as CT cell fractions increased. We generated nine recombinant ZIKVs containing combinations of the variants of interest. None of these viruses showed increased infection of Culex cells or mosquitoes, demonstrating that variants associated with passaging were not specific to increased Culex infection. These results reveal the challenge of a virus adapting to a new host, even when pushed to adapt artificially. Importantly, they also demonstrate that while ZIKV may occasionally infect Culex mosquitoes, Aedes mosquitoes likely drive transmission and human risk. IMPORTANCE ZIKV is mainly transmitted between humans by Aedes mosquitoes. In nature, ZIKV-infected Culex mosquitoes have been found, and ZIKV infrequently infects Culex mosquitoes under laboratory conditions. Yet, most studies show that Culex mosquitoes are not competent vectors for ZIKV. We attempted to adapt ZIKV to Culex cells to identify viral determinants of species specificity. We sequenced ZIKV after it was passaged on a mixture of Aedes and Culex cells and found that it acquired many variants. We generated recombinant viruses containing combinations of the variants of interest to determine if any of these changes enhance infection in Culex cells or mosquitoes. Recombinant viruses did not show increased infection in Culex cells or mosquitoes, but some variants increased infection in Aedes cells, suggesting adaptation to those cells instead. These results reveal that arbovirus species specificity is complex, and that virus adaptation to a new genus of mosquito vectors likely requires multiple genetic changes.

An evolutionary trade-off between host immunity and metabolism drives fatty liver in male mice.

Adaptations to infectious and dietary pressures shape mammalian physiology and disease risk. How such adaptations affect sex-biased diseases remains insufficiently studied. In this study, we show that sex-dependent hepatic gene programs confer a robust (~300%) survival advantage for male mice during lethal bacterial infection. The transcription factor B cell lymphoma 6 (BCL6), which masculinizes hepatic gene expression at puberty, is essential for this advantage. However, protection by BCL6 protein comes at a cost during conditions of dietary excess, which result in overt fatty liver and glucose intolerance in males. Deleting hepatic BCL6 reverses these phenotypes but markedly lowers male survival during infection, thus establishing a sex-dependent trade-off between host defense and metabolic systems. Our findings offer strong evidence that some current sex-biased diseases are rooted in ancient evolutionary trade-offs between immunity and metabolism.

Host adaptation of codon usage in SARS-CoV-2 from mammals indicates potential natural selection and viral fitness.

SARS-CoV-2 infection, which is the cause of the COVID-19 pandemic, has expanded across various animal hosts, and the virus can be transmitted particularly efficiently in minks. It is still not clear how SARS-CoV-2 is selected and evolves in its hosts, or how mutations affect viral fitness. In this report, sequences of SARS-CoV-2 isolated from human and animal hosts were analyzed, and the binding energy and capacity of the spike protein to bind human ACE2 and the mink receptor were compared. Codon adaptation index (CAI) analysis indicated the optimization of viral codons in some animals such as bats and minks, and a neutrality plot demonstrated that natural selection had a greater influence on some SARS-CoV-2 sequences than mutational pressure. Molecular dynamics simulation results showed that the mutations Y453F and N501T in mink SARS-CoV-2 could enhance the binding of the viral spike to the mink receptor, indicating the involvement of these mutations in natural selection and viral fitness. Receptor binding analysis revealed that the mink SARS-CoV-2 spike interacted more strongly with the mink receptor than the human receptor. Tracking the variations and codon bias of SARS-CoV-2 is helpful for understanding the fitness of the virus in virus transmission, pathogenesis, and immune evasion.

Distinct evolutionary trajectories of SARS-CoV-2-interacting proteins in bats and primates identify important host determinants of COVID-19.

The coronavirus disease 19 (COVID-19) pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a coronavirus that spilled over from the bat reservoir. Despite numerous clinical trials and vaccines, the burden remains immense, and the host determinants of SARS-CoV-2 susceptibility and COVID-19 severity remain largely unknown. Signatures of positive selection detected by comparative functional genetic analyses in primate and bat genomes can uncover important and specific adaptations that occurred at virus-host interfaces. We performed high-throughput evolutionary analyses of 334 SARS-CoV-2-interacting proteins to identify SARS-CoV adaptive loci and uncover functional differences between modern humans, primates, and bats. Using DGINN (Detection of Genetic INNovation), we identified 38 bat and 81 primate proteins with marks of positive selection. Seventeen genes, including the ACE2 receptor, present adaptive marks in both mammalian orders, suggesting common virus-host interfaces and past epidemics of coronaviruses shaping their genomes. Yet, 84 genes presented distinct adaptations in bats and primates. Notably, residues involved in ubiquitination and phosphorylation of the inflammatory RIPK1 have rapidly evolved in bats but not primates, suggesting different inflammation regulation versus humans. Furthermore, we discovered residues with typical virus-host arms race marks in primates, such as in the entry factor TMPRSS2 or the autophagy adaptor FYCO1, pointing to host-specific in vivo interfaces that may be drug targets. Finally, we found that FYCO1 sites under adaptation in primates are those associated with severe COVID-19, supporting their importance in pathogenesis and replication. Overall, we identified adaptations involved in SARS-CoV-2 infection in bats and primates, enlightening modern genetic determinants of virus susceptibility and severity.

Persisting uropathogenic Escherichia coli lineages show signatures of niche-specific within-host adaptation mediated by mobile genetic elements.

Large-scale genomic studies have identified within-host adaptation as a hallmark of bacterial infections. However, the impact of physiological, metabolic, and immunological differences between distinct niches on the pathoadaptation of opportunistic pathogens remains elusive. Here, we profile the within-host adaptation and evolutionary trajectories of 976 isolates representing 119 lineages of uropathogenic Escherichia coli (UPEC) sampled longitudinally from both the gastrointestinal and urinary tracts of 123 patients with urinary tract infections. We show that lineages persisting in both niches within a patient exhibit increased allelic diversity. Habitat-specific selection results in niche-specific adaptive mutations and genes, putatively mediating fitness in either environment. Within-lineage inter-habitat genomic plasticity mediated by mobile genetic elements (MGEs) provides the opportunistic pathogen with a mechanism to adapt to the physiological conditions of either habitat, and reduced MGE richness is associated with recurrence in gut-adapted UPEC lineages. Collectively, our results establish niche-specific adaptation as a driver of UPEC within-host evolution.

Malaria protection due to sickle haemoglobin depends on parasite genotype.

Host genetic factors can confer resistance against malaria1, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member2-4 PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

The FUR-like regulators PerRA and PerRB integrate a complex regulatory network that promotes mammalian host-adaptation and virulence of Leptospira interrogans.

Leptospira interrogans, the causative agent of most cases of human leptospirosis, must respond to myriad environmental signals during its free-living and pathogenic lifestyles. Previously, we compared L. interrogans cultivated in vitro and in vivo using a dialysis membrane chamber (DMC) peritoneal implant model. From these studies emerged the importance of genes encoding the Peroxide responsive regulators PerRA and PerRB. First described in in Bacillus subtilis, PerRs are widespread in Gram-negative and -positive bacteria, where regulate the expression of gene products involved in detoxification of reactive oxygen species and virulence. Using perRA and perRB single and double mutants, we establish that L. interrogans requires at least one functional PerR for infectivity and renal colonization in a reservoir host. Our finding that the perRA/B double mutant survives at wild-type levels in DMCs is noteworthy as it demonstrates that the loss of virulence is not due to a metabolic lesion (i.e., metal starvation) but instead reflects dysregulation of virulence-related gene products. Comparative RNA-Seq analyses of perRA, perRB and perRA/B mutants cultivated within DMCs identified 106 genes that are dysregulated in the double mutant, including ligA, ligB and lvrA/B sensory histidine kinases. Decreased expression of LigA and LigB in the perRA/B mutant was not due to loss of LvrAB signaling. The majority of genes in the perRA and perRB single and double mutant DMC regulons were differentially expressed only in vivo, highlighting the importance of host signals for regulating gene expression in L. interrogans. Importantly, the PerRA, PerRB and PerRA/B DMC regulons each contain multiple genes related to environmental sensing and/or transcriptional regulation. Collectively, our data suggest that PerRA and PerRB are part of a complex regulatory network that promotes host adaptation by L. interrogans within mammals.

Behavioral and Physiological Plasticity Provides Insights into Molecular Based Adaptation Mechanism to Strain Shift in Spodoptera frugiperda.

How herbivorous insects adapt to host plants is a key question in ecological and evolutionary biology. The fall armyworm, (FAW) Spodoptera frugiperda (J.E. Smith), although polyphagous and a major pest on various crops, has been reported to have a rice and corn (maize) feeding strain in its native range in the Americas. The species is highly invasive and has recently established in China. We compared behavioral changes in larvae and adults of a corn population (Corn) when selected on rice (Rice) and the molecular basis of these adaptational changes in midgut and antennae based on a comparative transcriptome analysis. Larvae of S. frugiperda reared on rice plants continuously for 20 generations exhibited strong feeding preference for with higher larval performance and pupal weight on rice than on maize plants. Similarly, females from the rice selected population laid significantly more eggs on rice as compared to females from maize population. The most highly expressed DEGs were shown in the midgut of Rice vs. Corn. A total of 6430 DEGs were identified between the populations mostly in genes related to digestion and detoxification. These results suggest that potential adaptations for feeding on rice crops, may contribute to the current rapid spread of fall armyworm on rice crops in China and potentially elsewhere. Consistently, highly expressed DEGs were also shown in antennae; a total of 5125 differentially expressed genes (DEGs) s were identified related to the expansions of major chemosensory genes family in Rice compared to the Corn feeding population. These results not only provide valuable insight into the molecular mechanisms in host plants adaptation of S. frugiperda but may provide new gene targets for the management of this pest.

The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages.

The independent emergence late in 2020 of the B.1.1.7, B.1.351, and P.1 lineages of SARS-CoV-2 prompted renewed concerns about the evolutionary capacity of this virus to overcome public health interventions and rising population immunity. Here, by examining patterns of synonymous and non-synonymous mutations that have accumulated in SARS-CoV-2 genomes since the pandemic began, we find that the emergence of these three "501Y lineages" coincided with a major global shift in the selective forces acting on various SARS-CoV-2 genes. Following their emergence, the adaptive evolution of 501Y lineage viruses has involved repeated selectively favored convergent mutations at 35 genome sites, mutations we refer to as the 501Y meta-signature. The ongoing convergence of viruses in many other lineages on this meta-signature suggests that it includes multiple mutation combinations capable of promoting the persistence of diverse SARS-CoV-2 lineages in the face of mounting host immune recognition.

Host adaptation in gut Firmicutes is associated with sporulation loss and altered transmission cycle.

BACKGROUND: Human-to-human transmission of symbiotic, anaerobic bacteria is a fundamental evolutionary adaptation essential for membership of the human gut microbiota. However, despite its importance, the genomic and biological adaptations underpinning symbiont transmission remain poorly understood. The Firmicutes are a dominant phylum within the intestinal microbiota that are capable of producing resistant endospores that maintain viability within the environment and germinate within the intestine to facilitate transmission. However, the impact of host transmission on the evolutionary and adaptive processes within the intestinal microbiota remains unknown. RESULTS: We analyze 1358 genomes of Firmicutes bacteria derived from host and environment-associated habitats. Characterization of genomes as spore-forming based on the presence of sporulation-predictive genes reveals multiple losses of sporulation in many distinct lineages. Loss of sporulation in gut Firmicutes is associated with features of host-adaptation such as genome reduction and specialized metabolic capabilities. Consistent with these data, analysis of 9966 gut metagenomes from adults around the world demonstrates that bacteria now incapable of sporulation are more abundant within individuals but less prevalent in the human population compared to spore-forming bacteria. CONCLUSIONS: Our results suggest host adaptation in gut Firmicutes is an evolutionary trade-off between transmission range and colonization abundance. We reveal host transmission as an underappreciated process that shapes the evolution, assembly, and functions of gut Firmicutes.

The dinucleotide composition of the Zika virus genome is shaped by conflicting evolutionary pressures in mammalian hosts and mosquito vectors.

Most vertebrate RNA viruses show pervasive suppression of CpG and UpA dinucleotides, closely resembling the dinucleotide composition of host cell transcriptomes. In contrast, CpG suppression is absent in both invertebrate mRNA and RNA viruses that exclusively infect arthropods. Arthropod-borne (arbo) viruses are transmitted between vertebrate hosts by invertebrate vectors and thus encounter potentially conflicting evolutionary pressures in the different cytoplasmic environments. Using a newly developed Zika virus (ZIKV) model, we have investigated how demands for CpG suppression in vertebrate cells can be reconciled with potentially quite different compositional requirements in invertebrates and how this affects ZIKV replication and transmission. Mutant viruses with synonymously elevated CpG or UpA dinucleotide frequencies showed attenuated replication in vertebrate cell lines, which was rescued by knockout of the zinc-finger antiviral protein (ZAP). Conversely, in mosquito cells, ZIKV mutants with elevated CpG dinucleotide frequencies showed substantially enhanced replication compared to wild type. Host-driven effects on virus replication attenuation and enhancement were even more apparent in mouse and mosquito models. Infections with CpG- or UpA-high ZIKV mutants in mice did not cause typical ZIKV-induced tissue damage and completely protected mice during subsequent challenge with wild-type virus, which demonstrates their potential as live-attenuated vaccines. In contrast, the CpG-high mutants displayed enhanced replication in Aedes aegypti mosquitoes and a larger proportion of mosquitoes carried infectious virus in their saliva. These findings show that mosquito cells are also capable of discriminating RNA based on dinucleotide composition. However, the evolutionary pressure on the CpG dinucleotides of viral genomes in arthropod vectors directly opposes the pressure present in vertebrate host cells, which provides evidence that an adaptive compromise is required for arbovirus transmission. This suggests that the genome composition of arbo flaviviruses is crucial to maintain the balance between high-level replication in the vertebrate host and persistent replication in the mosquito vector.

Genetic variation analyses indicate conserved SARS-CoV-2-host interaction and varied genetic adaptation in immune response factors in modern human evolution.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a pandemic as of early 2020. Upon infection, SARS-CoV-2 attaches to its receptor, that is, angiotensin-converting enzyme 2 (ACE2), on the surface of host cells and is then internalized into host cells via enzymatic machineries. This subsequently stimulates immune response factors. Since the host immune response and severity of COVID-19 vary among individuals, genetic risk factors for severe COVID-19 cases have been investigated. Our research group recently conducted a survey of genetic variants among SARS-CoV-2-interacting molecules across populations, noting near absence of difference in allele frequency spectrum between populations in these genes. Recent genome-wide association studies have identified genetic risk factors for severe COVID-19 cases in a segment of chromosome 3 that involves six genes encoding three immune-regulatory chemokine receptors and another three molecules. The risk haplotype seemed to be inherited from Neanderthals, suggesting genetic adaptation against pathogens in modern human evolution. Therefore, SARS-CoV-2 uses highly conserved molecules as its virion interaction, whereas its immune response appears to be genetically biased in individuals to some extent. We herein review the molecular process of SARS-CoV-2 infection as well as our further survey of genetic variants of its related immune effectors. We also discuss aspects of modern human evolution.

Host adaptation and microbial competition drive Ralstonia solanacearum phylotype I evolution in the Republic of Korea.

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) threatens the cultivation of important crops worldwide. We sequenced 30 RSSC phylotype I (R. pseudosolanacearum) strains isolated from pepper (Capsicum annuum) and tomato (Solanum lycopersicum) across the Republic of Korea. These isolates span the diversity of phylotype I, have extensive effector repertoires and are subject to frequent recombination. Recombination hotspots among South Korean phylotype I isolates include multiple predicted contact-dependent inhibition loci, suggesting that microbial competition plays a significant role in Ralstonia evolution. Rapid diversification of secreted effectors presents challenges for the development of disease-resistant plant varieties. We identified potential targets for disease resistance breeding by testing for allele-specific host recognition of T3Es present among South Korean phyloype I isolates. The integration of pathogen population genomics and molecular plant pathology contributes to the development of location-specific disease control and development of plant cultivars with durable resistance to relevant threats.

Evolutionary Selection of the Nuclear Localization Signal in the Viral Nucleoprotein Leads to Host Adaptation of the Genus Orthobornavirus.

Adaptation of the viral life cycle to host cells is necessary for efficient viral infection and replication. This evolutionary process has contributed to the mechanism for determining the host range of viruses. Orthobornaviruses, members of the family Bornaviridae, are non-segmented, negative-strand RNA viruses, and several genotypes have been isolated from different vertebrate species. Previous studies revealed that some genotypes isolated from avian species can replicate in mammalian cell lines, suggesting the zoonotic potential of avian orthobornaviruses. However, the mechanism by which the host specificity of orthobornaviruses is determined has not yet been identified. In this study, we found that the infectivity of orthobornaviruses is not determined at the viral entry step, mediated by the viral glycoprotein and matrix protein. Furthermore, we demonstrated that the nuclear localization signal (NLS) sequence in the viral nucleoprotein (N) has evolved under natural selection and determines the host-specific viral polymerase activity. A chimeric mammalian orthobornavirus, which has the NLS sequence of avian orthobornavirus N, exhibited a reduced propagation efficiency in mammalian cells. Our findings indicated that nuclear transport of the viral N is a determinant of the host range of orthobornaviruses, providing insights into the evolution and host adaptation of orthobornaviruses.

Analysis of the Codon Usage Pattern of HA and NA Genes of H7N9 Influenza A Virus.

Novel H7N9 influenza virus transmitted from birds to human and, since March 2013, it has caused five epidemic waves in China. Although the evolution of H7N9 viruses has been investigated, the evolutionary changes associated with codon usage are still unclear. Herein, the codon usage pattern of two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA), was studied to understand the evolutionary changes in relation to host, epidemic wave, and pathogenicity. Both genes displayed a low codon usage bias, with HA higher than NA. The codon usage was driven by mutation pressure and natural selection, although the main contributing factor was natural selection. Additionally, the codon adaptation index (CAI) and deoptimization (RCDI) illustrated the strong adaptability of H7N9 to Gallus gallus. Similarity index (SiD) analysis showed that Homo sapiens posed a stronger selection pressure than Gallus gallus. Thus, we assume that this may be related to the gradual adaptability of the virus to human. In addition, the host strong selection pressure was validated based on CpG dinucleotide content. In conclusion, this study analyzed the usage of codons of two genes of H7N9 and expanded our understanding of H7N9 host specificity. This aids into the development of control measures against H7N9 influenza virus.

Genome comparisons suggest an association between Ceratocystis host adaptations and effector clusters in unique transposable element families.

Ceratocystis fimbriata is a host specific fungal pathogen of sweet potato (Ipomoea batatas). The closely related species, C. manginecans, is an important pathogen of trees (e.g. Acacia mangium and Mangifera indica) but has never been isolated from tuber crops. The genetic factors that determine the host range and host specificity of these species have not been determined. The aim of this study was to compare the genomes of C. fimbriata and C. manginecans in order to identify species-specific genetic differences that could be associated with host specificity. This included whole-genome alignments as well as comparisons of gene content and transposable elements (TEs). The genomes of the two species were found to be very similar, sharing similar catalogues of CAZymes, peptidases and lipases. However, the genomes of the two species also varied, harbouring species-specific genes (e.g. small secreted effectors, nutrient processing proteins and stress response proteins). A portion of the TEs identified (17%) had a unique distribution in each species. Transposable elements appeared to have played a prominent role in the divergence of the two species because they were strongly associated with chromosomal translocations and inversions as well as with unique genomic regions containing species-specific genes. Two large effector clusters, with unique TEs in each species, were identified. These effectors displayed non-synonymous mutations and deletions, conserved within a species, and could serve as mutational hot-spots for the development of host specificity in the two species.

Adaptation of H9N2 Influenza Viruses to Mammalian Hosts: A Review of Molecular Markers.

As the number of human infections with avian and swine influenza viruses continues to rise, the pandemic risk posed by zoonotic influenza viruses cannot be underestimated. Implementation of global pandemic preparedness efforts has largely focused on H5 and H7 avian influenza viruses; however, the pandemic threat posed by other subtypes of avian influenza viruses, especially the H9 subtype, should not be overlooked. In this review, we summarize the literature pertaining to the emergence, prevalence and risk assessment of H9N2 viruses, and add new molecular analyses of key mammalian adaptation markers in the hemagglutinin and polymerase proteins. Available evidence has demonstrated that H9N2 viruses within the Eurasian lineage continue to evolve, leading to the emergence of viruses with an enhanced receptor binding preference for human-like receptors and heightened polymerase activity in mammalian cells. Furthermore, the increased prevalence of certain mammalian adaptation markers and the enhanced transmissibility of selected viruses in mammalian animal models add to the pandemic risk posed by this virus subtype. Continued surveillance of zoonotic H9N2 influenza viruses, inclusive of close genetic monitoring and phenotypic characterization in animal models, should be included in our pandemic preparedness efforts.

H7N9 Influenza Virus Containing a Polybasic HA Cleavage Site Requires Minimal Host Adaptation to Obtain a Highly Pathogenic Disease Phenotype in Mice.

Low pathogenic avian influenza (LPAI) H7N9 viruses have recently evolved to gain a polybasic cleavage site in the hemagglutinin (HA) protein, resulting in variants with increased lethality in poultry that meet the criteria for highly pathogenic avian influenza (HPAI) viruses. Both LPAI and HPAI variants can cause severe disease in humans (case fatality rate of ~40%). Here, we investigated the virulence of HPAI H7N9 viruses containing a polybasic HA cleavage site (H7N9-PBC) in mice. Inoculation of mice with H7N9-PBC did not result in observable disease; however, mice inoculated with a mouse-adapted version of this virus, generated by a single passage in mice, caused uniformly lethal disease. In addition to the PBC site, we identified three other mutations that are important for host-adaptation and virulence in mice: HA (A452T), PA (D347G), and PB2 (M483K). Using reverse genetics, we confirmed that the HA mutation was the most critical for increased virulence in mice. Our study identifies additional disease determinants in a mammalian model for HPAI H7N9 virus. Furthermore, the ease displayed by the virus to adapt to a new host highlights the potential for H7N9-PBC viruses to rapidly acquire mutations that may enhance their risk to humans or other animal species.

Extensive chromosomal rearrangements and rapid evolution of novel effector superfamilies contribute to host adaptation and speciation in the basal ascomycetous fungi.

The basal ascomycetes in genus Taphrina have strict host specificity and coevolution with their host plants, making them appealing models for studying the genomic basis of ecological divergence and host adaption. We therefore performed genome sequencing and comparative genomics of different Taphrina species with distinct host ranges to reveal their evolution. We identified frequent chromosomal rearrangements and highly dynamic lineage-specific (LS) genomic regions in Taphrina genomes. The LS regions occur at the flanking regions of chromosomal breakpoints, and are greatly enriched for DNA repeats, non-core genes, and in planta up-regulated genes. Furthermore, we identified hundreds of candidate secreted effector proteins (CSEPs) that are commonly organized in gene clusters that form distinct AT-rich isochore-like regions. Nearly half of the CSEPs constitute two novel superfamilies with modular structures unique to Taphrina. These CSEPs are commonly up-regulated during infection, enriched in the LS regions, evolved faster, and underwent extensive gene gain and loss in different species. In addition to displaying signatures of positive selection, functional characterization of selected CSEP genes confirmed their roles in suppression of plant defence responses. Overall, our results showed that extensive chromosomal rearrangements and rapidly evolving CSEP superfamilies play important roles in speciation and host adaptation in the early-branching ascomycetous fungi.