Functional annotation of a divergent genome using sequence and structure-based similarity.

BACKGROUND: Microsporidia are a large taxon of intracellular pathogens characterized by extraordinarily streamlined genomes with unusually high sequence divergence and many species-specific adaptations. These unique factors pose challenges for traditional genome annotation methods based on sequence similarity. As a result, many of the microsporidian genomes sequenced to date contain numerous genes of unknown function. Recent innovations in rapid and accurate structure prediction and comparison, together with the growing amount of data in structural databases, provide new opportunities to assist in the functional annotation of newly sequenced genomes. RESULTS: In this study, we established a workflow that combines sequence and structure-based functional gene annotation approaches employing a ChimeraX plugin named ANNOTEX (Annotation Extension for ChimeraX), allowing for visual inspection and manual curation. We employed this workflow on a high-quality telomere-to-telomere sequenced tetraploid genome of Vairimorpha necatrix. First, the 3080 predicted protein-coding DNA sequences, of which 89% were confirmed with RNA sequencing data, were used as input. Next, ColabFold was used to create protein structure predictions, followed by a Foldseek search for structural matching to the PDB and AlphaFold databases. The subsequent manual curation, using sequence and structure-based hits, increased the accuracy and quality of the functional genome annotation compared to results using only traditional annotation tools. Our workflow resulted in a comprehensive description of the V. necatrix genome, along with a structural summary of the most prevalent protein groups, such as the ricin B lectin family. In addition, and to test our tool, we identified the functions of several previously uncharacterized Encephalitozoon cuniculi genes. CONCLUSION: We provide a new functional annotation tool for divergent organisms and employ it on a newly sequenced, high-quality microsporidian genome to shed light on this uncharacterized intracellular pathogen of Lepidoptera. The addition of a structure-based annotation approach can serve as a valuable template for studying other microsporidian or similarly divergent species.

Differences in structure and hibernation mechanism highlight diversification of the microsporidian ribosome.

Assembling and powering ribosomes are energy-intensive processes requiring fine-tuned cellular control mechanisms. In organisms operating under strict nutrient limitations, such as pathogenic microsporidia, conservation of energy via ribosomal hibernation and recycling is critical. The mechanisms by which hibernation is achieved in microsporidia, however, remain poorly understood. Here, we present the cryo-electron microscopy structure of the ribosome from Paranosema locustae spores, bound by the conserved eukaryotic hibernation and recycling factor Lso2. The microsporidian Lso2 homolog adopts a V-shaped conformation to bridge the mRNA decoding site and the large subunit tRNA binding sites, providing a reversible ribosome inactivation mechanism. Although microsporidian ribosomes are highly compacted, the P. locustae ribosome retains several rRNA segments absent in other microsporidia, and represents an intermediate state of rRNA reduction. In one case, the near complete reduction of an expansion segment has resulted in a single bound nucleotide, which may act as an architectural co-factor to stabilize a protein-protein interface. The presented structure highlights the reductive evolution in these emerging pathogens and sheds light on a conserved mechanism for eukaryotic ribosome hibernation.

Molecular and morphological characterisation of a novel microsporidian species, Tubulinosema suzukii, infecting Drosophila suzukii (Diptera: Drosophilidae).

A microsporidium showing morphological characteristics typical of a Tubulinosema species was discovered in Drosophila suzukii. All developmental stages were diplokaryotic and grew in direct contact with the host cell cytoplasm. Spores from fresh preparations were ovoid to slightly pyriform and measured 4.29 × 2.47 µm in wet mount preparations. The spore wall consisted of a 125 nm thick endospore covered by a double layered exospore of 39 nm and 18 nm. The polar filament measured 67 µm in length, was slightly anisofilar and was arranged in ten coils in one or rarely two rows. The two posterior coils were 95 nm in diameter while the anterior coils were 115 nm in diameter. Early developmental stages were surrounded by electron-dense, 35.3 nm diameter, surface ornaments scattered over the membrane. Tubular elements with diameters of approximately 75 nm were seen attaching to the periphery of meronts and sporonts. Tissues infected included fat body, midgut and muscle. A 1915 bp rDNA fragment, covering the small subunit (SSU), the internal transcribed spacer (ITS) and the 5' end of the large subunit ribosomal DNA, was amplified by PCR and sequenced. Phylogenetic analyses of the SSU rDNA fragment revealed closest relationship to Tubulinosema pampeana (Host: Bombus atratus, South America) and Tubulinosema loxostegi (Host: Loxostege sticticalis, ubiquitous), but using the complete dataset of SSU-ITS-LSU rDNA genes revealed T. hippodamiae (Host: Hippodamiae convergens) as the most closely related species. Based on the morphological and genetic features a new species, Tubulinosema suzukii sp. nov., is proposed for this microsporidium isolated from D. suzukii.

A formal redefinition of the genera Nosema and Vairimorpha (Microsporidia: Nosematidae) and reassignment of species based on molecular phylogenetics.

The microsporidian genera Nosema and Vairimorpha comprise a clade described from insects. Currently the genus Nosema is defined as having a dimorphic life cycle characterized by diplokaryotic stages and diplosporoblastic sporogony with two functionally and morphologically distinct spore types ("early" or "primary" and "environmental"). The Vairimorpha life cycle, in addition to a Nosema-type diplokaryotic sporogony, includes an octosporoblastic sporogony producing eight uninucleate spores (octospores) within a sporophorous vesicle. Molecular phylogeny, however, has clearly demonstrated that the genera Nosema and Vairimorpha, characterized by the absence or presence of uninucleate octospores, respectively, represent two polyphyletic taxa, and that octosporogony is turned on and off frequently within taxa, depending on environmental factors such as host species and rearing temperature. In addition, recent studies have shown that both branches of the Vairimorpha-Nosema clade contain species that are uninucleate throughout their life cycle. The SSU rRNA gene sequence data reveal two distinct clades, those closely related to Vairimorpha necatrix, the type species for the genus Vairimorpha, and those closely related to Nosema bombycis, the type species for the genus Nosema. Here, we redefine the two genera, giving priority to molecular character states over those observed at the developmental, structural or ultrastructural levels and present a list of revised species designations. Using this approach, a series of species are renamed (combination novum) and members of two genera, Rugispora and Oligosporidium, are reassigned to Vairimorpha because of their phylogenetic position. Moreover, the family Nosematidae is redefined and includes the genera Nosema and Vairimorpha comprising a monophyletic lineage of Microsporidia.

Divergence of a Tandem Duplication of Manganese Superoxide Dismutase in Nosema bombycis.

Manganese superoxide dismutase (MnSOD) is a key enzyme in the protection of cells from oxidative stress. A tandem duplication of the MnSOD gene (NbMnSOD1 and NbMnSOD2) in the genome of Nosema bombycis, a parasite of the silkworm Bombyx mori, was previously identified. Here, we compare the protein structures of NbMnSOD1 and NbMnSOD2 and characterize these two proteins in terms of cellular localization, timing of transcription, protein structure, and enzyme activity. Despite a similarity in the primary sequence of NbMnSOD1 and NbMnSOD2, the latter shows a remarkable degree of amino acid sequence difference on the protein's surface and in the active site, where there is a substitution of a phenylalanine for a histidine in NbMnSOD2. Immuno-electron microscopy demonstrates that NbMnSOD1 is present in the cytosol of mature spores, whereas NbMnSOD2 is localized on the polar tube and the spore wall. Immunofluorescence confirms the localization of NbMnSOD2 on the polar tube of the germinated spore. Quantitative measurement of gene expression (qRT-PCR) demonstrates production of both alleles during the first day of infection followed by a dramatic decrease during the second to fourth day of infection. From the fifth day onward, the two alleles show a complementary pattern of expression. The qRT-PCR of the host manganese superoxide dismutase (BmMnSOD) shows a notable increase in transcription upon infection, leading to a three-fold spike by the first day of infection, followed by a decrease in transcription. Measurement of overall MnSOD activity shows a similar peak at day 1 followed by a decrease to a constant rate of enzyme activity. The differences in cellular localization and pattern of gene expression of NbMnSOD2 compared to NbMnSOD1, as well as the differences in protein structure seen for NbMnSOD2 compared to other microsporidial MnSODs, strongly suggest a unique, recently evolved role for NbMnSOD2.

Envenomation by Trachelas tranquillus (Araneae: Corrinidae) in Connecticut.

We report a case of envenomation by Trachelas tranquillus (Hentz) in Connecticut in late September 2013. The bitten subject, a 50 yr-old-female Caucasian, reported a painful wasp-like sting and brushed the spider from her leg. An erythematous macule formed at the site of the bite. The macule was gone by the next day and there was no associated necrosis. The spider was collected and brought to our laboratory for identification. This is the second confirmed case of envenomation by T. tranquillus and the only case reported from Connecticut.

Searching for convergent evolution in manganese superoxidase dismutase using hydrophobic cluster analysis.

There are numerous examples of convergent evolution in nature. Major ecological adaptations such as flight, loss of limbs in vertebrates, pesticide resistance, adaptation to a parasitic way of life, etc., have all evolved more than once, as seen by their analogous functions in separate taxa. But what about protein evolution? Does the environment have a strong enough influence on intracellular processes that enzymes and other functional proteins play, to evolve similar functional roles separately in different organisms? Manganese Superoxide Dismutase (MnSOD) is a manganesedependant metallo-enzyme which plays a crucial role in protecting cells from anti-oxidative stress by eliminating reactive (superoxide) oxygen species. It is a ubiquitous housekeeping enzyme found in nearly all organisms. In this study we compare phylogenies based on MnSOD protein sequences to those based on scores from Hydrophobic Cluster Analysis (HCA). We calculated HCA similarity values for each pair of taxa to obtain a pair-wise distance matrix. A UPGMA tree based on the HCA distance matrix and a common tree based on the primary protein sequence for MnSOD was constructed. Differences between these two trees within animals, enterobacteriaceae, planctomycetes and cyanobacteria are presented and cited as possible examples of convergence. We note that several residue changes result in changes in hydrophobicity at positions which apparently are under the effect of positive selection.

Field Isolation and Laboratory Vector-Host Studies of Brazoran Virus (Peribunyaviridae: Orthobunyavirus) from Florida.

Brazoran virus was first isolated from Culex mosquitoes in Texas in 2012, yet little is known about this virus. We report the isolation of this virus from Culex erraticus from southern Florida during 2016. The Florida strain had a nucleotide identity of 96.3% (S segment), 99.1% (M segment), and 95.8% (L segment) to the Texas isolate. Culex quinquefasciatus and Aedes aegypti colonies were subsequently fed virus blood meals to determine their vector competence for Brazoran virus. Culex quinquefasciatus was susceptible to midgut infection, but few mosquitoes developed disseminated infections. Aedes aegypti supported disseminated infection, but virus transmission could not be demonstrated. Suckling mice became infected by intradermal inoculation without visible disease signs. The virus was detected in multiple mouse tissues but rarely infected the brain. This study documents the first isolation of Brazoran virus outside of Texas. Although this virus infected Ae. aegypti and Cx. quinquefasciatus in laboratory trials, their vector competence could not be demonstrated, suggesting they are unlikely vectors of Brazoran virus.

Morphological and molecular characterization of a microsporidian parasite, Takaokaspora nipponicus n. gen., n. sp. from the invasive rock pool mosquito, Ochlerotatus japonicus japonicus.

A new genus and species of Microsporidia, Takaokaspora nipponicus n. gen., n. sp. is described from Ochlerotatus japonicus japonicus (Theobald) and Ochlerotatus hatorii (Yamada) based on light microscope and ultrastructural morphology, developmental features, transmission cycles and comparative sequence analyses of the small subunit ribosomal DNA (SSU rDNA). The microsporidium is both vertically and horizontally transmitted, exhibits dimorphic development alternating between diplokaryotic and monokaryotic stages and produces two morphologically distinct spores, one in larvae and another in adult females. Horizontal transmission of infection to larval mosquitoes occurs via direct oral ingestion of uninucleate spores that are produced in vertically-infected larval hosts. Development in horizontally-infected hosts is diplokaryotic following karyokinesis of uninucleate schizonts and binary fission to produce small (4.3µm × 2.0µm) membrane free, ovoid, binucleate spores that are confined to adult female reproductive tissues (ovariole sheath and oviducts). Vertical transmission of the microsporidium from adult females to larval progeny takes place via surface contamination of the egg (transovum). Microsporidian development in vertically-infected larvae is haplophasic with unpaired nuclei throughout, producing rosette-shaped sporogonial plasmodia contained within a thin non-persistent sporophorous vesicle and culminating in the formation of membrane free, uninucleate, conical spores (7.0µm×2.8µm). Development is confined to host fat body tissue which appears as swollen white masses in the thorax and selected segments of the abdomen causing larvae to appear abnormally distorted and results in death during the third and fourth instar stages. The SSU rDNA sequences obtained from the two morphologically identical microsporidia isolated from Oc. j. japonicus and Oc. hatorii were nearly identical and unique when compared with GenBank entries of all other mosquito-parasitic species. Phylogenetic trees constructed by Maximum Parsimony, Maximum Likelihood and bootstrap analyses using the Neighbor Joining search parameter yielded similar typologies. In each case, the novel microsporidium was the sister group to the clade containing Parathelohania species from Anopheles mosquitoes and the monotypic Novothelohania ovalae from Ochlerotatus caspius showing approximately 10-13% sequence divergence to those two genera providing strong support for establishment as a separate genus.

Intraspecific polymorphism of rDNA among five Nosema bombycis isolates from different geographic regions in China.

The microsporidian Nosema bombycis is the causative agent of pébrine, a highly infectious disease of the silkworm Bombyx mori. Three regions of the multicopy rDNA gene were examined in order to investigate the relationships among five Nosema isolates from various regions of China. Ribosomal DNA alleles are present on each of the 18 chromosomes of N. bombycis and show a high degree of variation. In this study the small subunit (SSU) rDNA, internal transcribed spacer (ITS) and intergenic spacer (IGS) regions for up to 10 different rDNA copies from each N. bombycis isolate were cloned and sequenced. As expected we see greater polymorphism in the ITS region (88 variable sites in 179 nucleotides) and IGS (200 variable sites in 279 nucleotides) than in the SSU rDNA (24 variable sites in 1232 nucleotides). Phylogenetic analysis shows greater differences between alleles within an isolate than between the same alleles from different isolates. The data reveal two very different groups, one from the Sichuan province and the other with a broad distribution including four provinces in southeast China and Japan. The Sichuan isolate does not have any rDNA alleles with sequences identical to those in the other isolates, implying that it is a separate, non-intermixing, population or perhaps a separate species from the other isolates. In light of the polymorphic nature of the rDNA alleles in N. bombycis and their presence on every chromosome, the rDNA gene may be useful for understanding the movement and ultimately the source of pébrine infections.

Ultrastructural characterization and comparative phylogenetic analysis of new microsporidia from Siberian mosquitoes: evidence for coevolution and host switching.

A survey of mosquito larvae infected with microsporidia was conducted from 2005 to 2008 in the Tomsk, Kemerovo and Novosibirsk regions of western Siberia, Russia. Twenty-one morphologically and genetically unique species of microsporidia were isolated from nine species of Anopheles, Aedes, Culex and Ochlerotatus mosquitoes including: (1) 14 proposed new species of Amblyospora (A. bakcharia, A. baritia, A. bogashovia, A. chulymia, A. hristinia, A. jurginia, A. kazankia, A. mavlukevia, A. mocrushinia, A. modestium, A. salairia, A. severinia, A. shegaria, and A. timirasia); (2) a newly proposed genus and species, Novothelohaniaovalae and; (3) six species of Amblyospora (A. flavescens, A. kolarovi, A. rugosa), Parathelohania (P. divulgata and P. tomski) and Trichoctosporea (T. pygopellita) from which gene sequences had not been previously obtained. Detailed ultrastructure of meiospores revealed unique cytological features associated with the length, arrangement and ratio of broad to narrow coils of the polar filament, comparative thickness of the exospore and endospore, and overall size of each species reaffirming their value in distinguishing taxonomic relationships. SSU rDNA sequences obtained from each species of microsporidia were unique when compared with GenBank entries. Phylogenetic trees constructed by Maximum Parsimony, Maximum Likelihood and Neighbor Joining analyses yielded similar topologies with a high degree of congruence between parasite and host at the generic level. Species that parasitize Aedes/Ochlerotatus and Culex mosquitoes segregate into distinct monophyletic groupings mirroring their host phylogeny, while species from Anopheles mosquitoes group as a sister clade basal to the entire group of mosquito-parasitic microsporidia as their Anopheles hosts cluster as a sister clade to the entire group of culicine mosquitoes. This provides strong evidence for host-parasite coevolution by descent at the generic level and limited host lineage switching between unrelated taxa. Among parasites of Aedes/Ochlerotatus and Anopheles mosquitoes, we found several instances where a single mosquito species serves as a host for two or more related species of microsporidia, an observation consistent with host switching and independent parasite speciation. Among the microsporidian parasites of Culex mosquitoes, we found only one parasite per host indicating a higher degree of host specificity and less host switching among parasites of this genus. Findings suggest a degree of host-parasite co-speciation with host switching occurring occasionally when the "normal" host is unavailable in the aquatic ecosystem. Frequency of host switching seems to be occurring in proportion to host relatedness and does not cross generic boundaries in this system.

Characterization of the Largest Secretory Protein Family, Ricin B Lectin-like Protein, in Nosema bombycis: Insights into Microsporidian Adaptation to Host.

Microsporidia are a group of obligate intracellular pathogens infecting nearly all animal phyla. The microsporidian Nosema bombycis has been isolated from several lepidopteran species, including the economy-important silkworms as well as several crop pests. Proteins secreted by parasites can be important virulent factors in modulating host pathways. Ricin is a two-chain lectin best known for its extreme vertebrate toxicity. Ricin B lectin-like proteins are widely distributed in microsporidia, especially in N. bombycis. In this study, we identify 52 Ricin B lectin-like proteins (RBLs) in N. bombycis. We show that the N. bombycis RBLs (NbRBLs) are classified into four subfamilies. The subfamily 1 was the most conserved, with all members having a Ricin B lectin domain and most members containing a signal peptide. The other three subfamilies were less conserved, and even lost the Ricin B lectin domain, suggesting that NbRBLs might be a multi-functional family. Our study here indicated that the NbRBL family had evolved by producing tandem duplications firstly and then expanded by segmental duplications, resulting in concentrated localizations mainly in three genomic regions. Moreover, based on RNA-seq data, we found that several Nbrbls were highly expressed during infection. Further, the results show that the NbRBL28 was secreted into host nucleus, where it promotes the expressions of genes involved in cell cycle progression. In summary, the great copy number, high divergence, and concentrated genome distribution of the NbRBLs demonstrated that these proteins might be adaptively evolved and played a vital role in the multi-host N. bombycis.

Structure of the reduced microsporidian proteasome bound by PI31-like peptides in dormant spores.

Proteasomes play an essential role in the life cycle of intracellular pathogens with extracellular stages by ensuring proteostasis in environments with limited resources. In microsporidia, divergent parasites with extraordinarily streamlined genomes, the proteasome complexity and structure are unknown, which limits our understanding of how these unique pathogens adapt and compact essential eukaryotic complexes. We present cryo-electron microscopy structures of the microsporidian 20S and 26S proteasome isolated from dormant or germinated Vairimorpha necatrix spores. The discovery of PI31-like peptides, known to inhibit proteasome activity, bound simultaneously to all six active sites within the central cavity of the dormant spore proteasome, suggests reduced activity in the environmental stage. In contrast, the absence of the PI31-like peptides and the existence of 26S particles post-germination in the presence of ATP indicates that proteasomes are reactivated in nutrient-rich conditions. Structural and phylogenetic analyses reveal that microsporidian proteasomes have undergone extensive reductive evolution, lost at least two regulatory proteins, and compacted nearly every subunit. The highly derived structure of the microsporidian proteasome, and the minimized version of PI31 presented here, reinforce the feasibility of the development of specific inhibitors and provide insight into the unique evolution and biology of these medically and economically important pathogens.

Envenomation by Steatoda borealis (Araneae: Theridiidae) in Connecticut, USA.

We report a case of envenomation by Steatoda borealis (Hentz, 1850) in Connecticut in late spring 2021. The bitten subject was a 60-yr-old female Caucasian. The bite occurred upon disturbing some books in her basement. She reported a painful stinging bite on her finger followed by the formation of an erythematous macule at the site of the bite which was gone the next day. There was no sign of necrosis. The spider was recovered directly at that time. The victim reported that later that night she experienced heart palpitations. This is the first report of envenomation by S. borealis, a common spider in Connecticut.

A tandem duplication of manganese superoxide dismutase in Nosema bombycis and its evolutionary origins.

Microsporidia are a group of obligate intracellular eukaryotic parasites with small genomes. They infect animals from a wide variety of phyla, including humans. Two manganese superoxide dismutase (MnSOD) genes, designated NbMnSOD1 and NbMnSOD2, were found to be organized in a tandem array within the Nosema bombycis genome. The genes, both 678 bp in length, were found to be more similar to each other than they are to homologous genes of other Microsporidia, suggesting that the tandem duplication occurred subsequent to the development of this lineage. Reverse transcript PCR shows that mRNA for both genes is present in the spores. Analysis of the primary structure, hydrophobic cluster analysis, target signal analysis, and phylogenetic analysis all indicate that NbMnSOD1 is dimeric and targeted to the cytosol. NbMnSOD2 seems to have changed more rapidly and is under less evolutionary constraint than NbMnSOD1 suggesting that NbMnSOD2 may function under different conditions or in different tissues of its host rather than simply resulting in an increase in expression. A phylogenetic analysis of MnSOD sequences from eukaryotes, Archaea, and bacteria shows the microsporidial MnSODs to be grouped with the bacteria suggesting a possible horizontal gene transfer.

Phylogenetic position of Octosporea muscaedomesticae (Microsporidia) and its relationship to Octosporea bayeri based on small subunit rDNA analysis.

Comparative phylogenetic analysis of the small subunit rDNA sequence of Octosporea muscaedomesticae (Flu, 1911) (type species) (Microsporidia) isolated from the blowfly Phormia regina (Diptera:Calliphoridae) is presented. Neighbor Joining bootstrap, Maximum Parsimony and Maximum Likelihood analyses with 38 microsporidian taxa representing five major clades of Microsporidia placed O. muscaedomesticae on a separate branch within a clade containing parasites of freshwater hosts. O. muscaedomesticae differed from Octosporea bayeri, a parasite of the microcrustacean, Daphnia magna (Cladocera:Daphniidae) by 29% demonstrating that the latter microsporidium is not closely related to the type species at the generic level, and should not be placed within the genus Octosporea, a conclusion that is further supported by morphological and developmental differences. Considering the number of disparately related hosts from which Octosporea species have been previously described based mostly on developmental and morphological characters it is likely that many will not fit the current definition of the genus, and it is possible that molecular analysis of these species will show that this genus as defined represents a polyphyletic grouping of unrelated taxa.

Occurrence of Nosema oryzaephili in Cryptolestes ferrugineus and transfer to the genus Paranosema.

A microsporidium that closely resembles Paranosema species at the level of the light microscope was isolated from the rusty grain beetle, Cryptolestes ferrugineus. It's identity as Nosema oryzaephili (originally described from Oryzaephilus surinamensis) was confirmed by comparison with a known isolate of N. oryzaephili based on spore size, small subunit rDNA sequence, and relative infectivity to O. surinamensis, Tribolium castaneum, and Ephestia kuehniella. Phylogenetic analysis of the small subunit rDNA indicates clearly that this species belongs in the genus Paranosema and thus the designation Paranosema oryzaephili (Burges, Canning and Hurst) is proposed. In spite of the abundance, economic importance, and world-wide distribution of C. ferrugineus, this is the first report of a microsporidial infection in this species. This is also the first report of P. oryzaephili in the new world.

Prevalence of Infection and Co-Infection and Presence of Rickettsial Endosymbionts in Ixodes Scapularis (Acari: Ixodidae) in Connecticut, USA.

Ixodes scapularis is currently known to transmit 7 pathogens responsible for Lyme disease, anaplasmosis, babesiosis, tick-borne relapsing fever, ehrlichiosis, and Powassan encephalitis. Ixodes scapularis can also be colonized by endosymbiotic bacteria including those in the genus of Rickettsia. We screened 459 I. scapularis ticks submitted to the Connecticut Agricultural Experiment Station Tick Testing Laboratory with the objectives to (1) examine differences in infection prevalence of Borrelia burgdorferi, Anaplasma phagocytophilum, Babesia microti, and Borrelia miyamotoi, (2) evaluate whether prevalence of co-infections occur at the same frequency that would be expected based on single infection, and (3) determine the presence of rickettsial endosymbionts in I. scapularis. The prevalence of infection in I. scapularis was highest with Bo. burgdorferi sensu lato (nymph = 45.8%; female = 47.0%), followed by A. phagocytophilum (nymph = 4.0%; female = 6.9%), Ba. microti (nymph = 5.7%; female = 4.7%), and Bo. miyamotoi (nymph = 0%; female = 7.3%). We also identified rickettsial endosymbionts in 93.3% of I. scapularis. Nymphs were significantly more likely to be infected with Bo. burgdorferi if they were infected with Ba. microti, whereas adult females were significantly more likely to be infected with Bo. burgdorferi if they were infected with A. phagocytophilum. Our study suggests that the infection prevalence of Bo. burgdorferi is not independent of other co-circulating pathogens and that there is a substantially higher infection of Bo. miyamotoi in I. scapularis females compared with nymphs in this study. High prevalence of infection and co-infection with multiple pathogens in I. scapularis highlights the public health consequences in Connecticut, a state endemic for Lyme and other tick-borne diseases.

Genomic analysis of Asian honeybee populations in China reveals evolutionary relationships and adaptation to abiotic stress.

The geographic and biological diversity of China has resulted in the differential adaptation of the eastern honeybee, Apis cerana, to these varied habitats. A. cerana were collected from 14 locations in China. Their genomes were sequenced, and nucleotide polymorphisms were identified at more than 9 million sites. Both STRUCTURE and principal component analysis placed the bees into seven groups. Phylogenomic analysis groups the honeybees into many of the same clusters with high bootstrap values (91%-100%). Populations from Tibet and South Yunnan are sister taxa and together represent the earliest diverging lineage included in this study. We propose that the evolutionary origin of A. cerana in China was in the southern region of Yunnan Province and expanded from there into the southeastern regions and into the northeastern mountain regions. The Cold-Temperate West Sichuan Plateau and Tropical Diannan populations were compared to identify genes under adaptive selection in these two habitats. Pathway enrichment analysis showing genes under selection, including the Hippo signaling pathway, GABAergic pathway, and trehalose-phosphate synthase, indicates that most genes under selection pressure are involved in the process of signal transduction and energy metabolism. qRT-PCR analysis reveals that one gene under selection, the AcVIAAT gene, involved in the GABAergic pathway, is responding to cold temperature stress. Through homologous recombination, we show that the AcVIAAT gene is able to replace the CNAG_01904 gene in the fungus Cryptococcus neoformans and that it makes the fungus less sensitive to conditions of oxidative stress and variations in temperature. Our results contribute to our understanding of the evolutionary origin of A. cerana in China and the molecular basis of environmental adaptation.

Evolutionary compaction and adaptation visualized by the structure of the dormant microsporidian ribosome.

Microsporidia are eukaryotic parasites that infect essentially all animal species, including many of agricultural importance1-3, and are significant opportunistic parasites of humans4. They are characterized by having a specialized infection apparatus, an obligate intracellular lifestyle5, rudimentary mitochondria and the smallest known eukaryotic genomes5-7. Extreme genome compaction led to minimal gene sizes affecting even conserved ancient complexes such as the ribosome8-10. In the present study, the cryo-electron microscopy structure of the ribosome from the microsporidium Vairimorpha necatrix is presented, which illustrates how genome compaction has resulted in the smallest known eukaryotic cytoplasmic ribosome. Selection pressure led to the loss of two ribosomal proteins and removal of essentially all eukaryote-specific ribosomal RNA (rRNA) expansion segments, reducing the rRNA to a functionally conserved core. The structure highlights how one microsporidia-specific and several repurposed existing ribosomal proteins compensate for the extensive rRNA reduction. The microsporidian ribosome is kept in an inactive state by two previously uncharacterized dormancy factors that specifically target the functionally important E-site, P-site and polypeptide exit tunnel. The present study illustrates the distinct effects of evolutionary pressure on RNA and protein-coding genes, provides a mechanism for ribosome inhibition and can serve as a structural basis for the development of inhibitors against microsporidian parasites.