The massive 340 megabase genome of Anisogramma anomala, a biotrophic ascomycete that causes eastern filbert blight of hazelnut.

BACKGROUND: The ascomycete fungus Anisogramma anomala causes Eastern Filbert Blight (EFB) on hazelnut (Corylus spp.) trees. It is a minor disease on its native host, the American hazelnut (C. americana), but is highly destructive on the commercially important European hazelnut (C. avellana). In North America, EFB has historically limited commercial production of hazelnut to west of the Rocky Mountains. A. anomala is an obligately biotrophic fungus that has not been grown in continuous culture, rendering its study challenging. There is a 15-month latency before symptoms appear on infected hazelnut trees, and only a sexual reproductive stage has been observed. Here we report the sequencing, annotation, and characterization of its genome. RESULTS: The genome of A. anomala was assembled into 108 scaffolds totaling 342,498,352 nt with a GC content of 34.46%. Scaffold N50 was 33.3 Mb and L50 was 5. Nineteen scaffolds with lengths over 1 Mb constituted 99% of the assembly. Telomere sequences were identified on both ends of two scaffolds and on one end of another 10 scaffolds. Flow cytometry estimated the genome size of A. anomala at 370 Mb. The genome exhibits two-speed evolution, with 93% of the assembly as AT-rich regions (32.9% GC) and the other 7% as GC-rich (57.1% GC). The AT-rich regions consist predominantly of repeats with low gene content, while 90% of predicted protein coding genes were identified in GC-rich regions. Copia-like retrotransposons accounted for more than half of the genome. Evidence of repeat-induced point mutation (RIP) was identified throughout the AT-rich regions, and two copies of the rid gene and one of dim-2, the key genes in the RIP mutation pathway, were identified in the genome. Consistent with its homothallic sexual reproduction cycle, both MAT1-1 and MAT1-2 idiomorphs were found. We identified a large suite of genes likely involved in pathogenicity, including 614 carbohydrate active enzymes, 762 secreted proteins and 165 effectors. CONCLUSIONS: This study reveals the genomic structure, composition, and putative gene function of the important pathogen A. anomala. It provides insight into the molecular basis of the pathogen's life cycle and a solid foundation for studying EFB.

First record of Ixodes keiransi (Acari: Ixodidae) in New Jersey, USA.

The hard tick, Ixodes keiransi Beati, Nava, Venzal, & Guglielmone, formerly the North American lineage of Ixodes affinis Neumann, is expanding its range northward along the US East Coast. In July 2023, we collected I. keiransi adult female and nymph in a single sampling event, suggesting its range now includes southern New Jersey. In this area, I. keiransi is sympatric with northern populations of Ixodes scapularis Say (Acari: Ixodidae), the primary vector of Lyme disease. Given its status as an enzootic vector of spirochaetes in the Borrelia burgdorferi sensu lato complex, proper differentiation of these 2 species will be critical for accurate estimates of entomological risk. Targeted surveillance should be implemented to monitor further I. keiransi expansion and to elucidate the phenology and enzootic role of this and other understudied Ixodes spp. in the northeastern United States.

Endogenous Viral Elements in Ixodid Tick Genomes.

The documentation of endogenous viral elements (EVEs; virus-derived genetic material integrated into the genome of a nonviral host) has offered insights into how arthropods respond to viral infection via RNA interference pathways. Small non-coding RNAs derived from EVE loci serve to direct RNAi pathways in limiting replication and infection from cognate viruses, thus benefiting the host's fitness and, potentially, vectorial capacity. Here we use informatic approaches to analyze nine available genome sequences of hard ticks (Acari: Ixodidae; Rhipicephalus sanguineus, R. microplus, R. annulatus, Ixodes ricinus, I. persulcatus, I. scapularis, Hyalomma asiaticum, Haemaphysalis longicornis, and Dermacentor silvarum) to identify endogenous viral elements and to illustrate the shared ancestry of all elements identified. Our results highlight a broad diversity of viral taxa as having given rise to 1234 identified EVEs in ticks, with Mononegavirales (specifically Rhabdoviridae) well-represented in this subset of hard ticks. Further investigation revealed extensive adintovirus integrations in several Ixodes species, the prevalence of Bunyavirales EVEs (notably not observed in mosquitoes), and the presence of several elements similar to known emerging human and veterinary pathogens. These results will inform subsequent work on current and past associations with tick species with regard to the viruses from which their "viral fossils" are derived and may serve as a reference for quality control of various tick-omics data that may suffer from misidentification of EVEs as viral genetic material.

Detection of multiple tick-borne pathogens in Ixodes scapularis from Hunterdon County, NJ, USA.

Several human pathogens vectored by the blacklegged tick (Ixodes scapularis Say; Acari: Ixodidae) are endemic in the state of New Jersey. Disease incidence data suggest that these conditions occur disproportionately in the northwestern portion of the state, including in the county of Hunterdon. We conducted active surveillance at three forested sites in Hunterdon County during 2020 and 2021, collecting 662 nymphal and adult I. scapularis. Ticks were tested for five pathogens by qPCR/qRT-PCR: Anaplasma phagocytophilum, Babesia microti, Borrelia burgdorferi, Borrelia miyamotoi, and Powassan virus (POWV) lineage 2. Over 2 years, 25.4% of nymphs and 58.4% of adults were found infected with at least one pathogen, with 10.6% of all ticks infected with more than one pathogen. We report substantial spatial and temporal variability of A. phagocytophilum and B. burgdorferi, with high relative abundance of the human-infective A. phagocytophilum variant Ap-ha. Notably, POWV was detected for the first time in Hunterdon, a county where human cases have not been reported. Based on comparisons with active surveillance initiatives in nearby counties, further investigation of non-entomological factors potentially influencing rates of tick-borne illness in Hunterdon is recommended.

Lone star ticks (Acari: Ixodidae) infected with Bourbon virus in New Jersey, USA.

Lone star ticks (Amblyomma americanum L.) are expanding within the northeast United States, a region historically focused on Ixodes scapularis-transmitted diseases. In Monmouth County, NJ, the shift has been dramatic, and lone star ticks now vastly outnumber blacklegged ticks. As a result, there is an enhanced need to focus on the potential health risks of A. americanum-transmitted pathogens, such as the emerging Heartland (HRTV) and Bourbon (BRBV) viruses. We screened 1,205 nymphal lone star ticks for HRTV and BRBV using RT-qPCR assays and detected BRBV in 3 ticks collected in Monmouth County, NJ, in 2021. Additionally, we sequenced a complete BRBV genome from a single infected specimen, finding 99.4% identity with human pathogenic isolates from the eastern-central United States. Our results have important public health implications for a region only recently becoming aware of public health risks posed by lone star ticks. Of note, we report successful detection of viral RNA in samples that were stored and intended for DNA preservation, for example, kept in ethanol at room temperature, which may reduce barriers for public health agencies seeking to expand their tick testing to include viruses.

Rickettsia and relapsing fever Borrelia in Alectorobius kelleyi (Ixodida: Argasidae) from peri domestic bats in the northeastern United States.

The soft ticks (Argasidae) are known vectors of human and animal pathogens around the globe and are relatively understudied. Our aim was to assess the presence of Rickettsia and Borrelia bacteria in Alectorobius kelleyi (Argasidae) parasitizing synanthropic bats in the highly urbanized northeastern United States. By collaborating with parasitologists, bat scientists and wildlife rehabilitators we were successful in obtaining A. kelleyi from five states. Since Argasid larvae will attach to their hosts for many days, most A. kelleyi examined (92%) were larvae collected from sick or injured big brown bats, Eptesicus fuscus, undergoing care at rehabilitation centers. In addition, we obtained adult A. kelleyi captured in residential living areas and trapped in attics. An in-depth analysis of a A. kelleyi found to be infected with a spotted fever group Rickettsia (SFGR) revealed a dual infection with a R. belli-like taxon (ancestral group) as well as an SFGR closely related to R. peacockii, likely the same previously found in A. kelleyi from Iowa and Kansas. We found that 36% of the A. kelleyi tested carried the SFGR. Furthermore, we detected a relapsing fever spirochete, likely Candidatus Borrelia johnsonii, in 25% of the A. kelleyi from Pennsylvania. While it is unclear if these bacteria constitute a health risk to either bats or humans, our study indicates that human exposure to ectoparasites infesting peridomestic wildlife should be considered in the epidemiology of tick-borne diseases.

Microhabitat modeling of the invasive Asian longhorned tick (Haemaphysalis longicornis) in New Jersey, USA.

The Asian longhorned tick (Haemaphysalis longicornis) is a vector of multiple arboviral and bacterial pathogens in its native East Asia and expanded distribution in Australasia. This species has both bisexual and parthenogenetic populations that can reach high population densities under favorable conditions. Established populations of parthenogenetic H. longicornis were detected in the eastern United States in 2017 and the possible range of this species at the continental level (North America) based on climatic conditions has been modeled. However, little is known about factors influencing the distribution of H. longicornis at geographic scales relevant to local surveillance and control. To examine the importance of local physiogeographic conditions such as geology, soil characteristics, and land cover on the distribution of H. longicornis we employed ecological niche modeling using three machine learning algorithms - Maxent, Random Forest (RF), and Generalized Boosting Method (GBM) to estimate probability of finding H. longicornis in a particular location in New Jersey (USA), based on environmental predictors. The presence of H. longicornis in New Jersey was positively associated with Piedmont physiogeographic province and two soil types - Alfisols and Inceptisols. Soil hydraulic conductivity was the most important predictor explaining H. longicornis habitat suitability, with more permeable sandy soils with higher hydraulic conductivity being less suitable than clay or loam soils. The models were projected over the state of New Jersey creating a probabilistic map of H. longicornis habitat suitability at a high spatial resolution of 90×90 meters. The model's sensitivity was 87% for locations sampled in 2017-2019 adding to the growing evidence of the importance of soil characteristics to the survival of ticks. For the 2020-2022 dataset the model fit was 57%, suggestive of spillover to less optimal habitats or, alternatively, heterogeneity in soil characteristics at the edges of broad physiographic zones. Further modeling should incorporate abundance and life-stage information as well as detailed characterization of the soil at collection sites. Once critical parameters that drive the survival and abundance of H. longicornis are identified they can be used to guide surveillance and control strategies for this invasive species.

Rapid Spread of Severe Fever with Thrombocytopenia Syndrome Virus by Parthenogenetic Asian Longhorned Ticks.

Severe fever with thrombocytopenia syndrome virus (SFTSV) is spreading rapidly in Asia. This virus is transmitted by the Asian longhorned tick (Haemaphysalis longicornis), which has parthenogenetically and sexually reproducing populations. Parthenogenetic populations were found in ≥15 provinces in China and strongly correlated with the distribution of severe fever with thrombocytopenia syndrome cases. However, distribution of these cases was poorly correlated with the distribution of populations of bisexual ticks. Phylogeographic analysis suggested that the parthenogenetic population spread much faster than bisexual population because colonization is independent of sexual reproduction. A higher proportion of parthenogenetic ticks was collected from migratory birds captured at an SFTSV-endemic area, implicating the contribution to the long-range movement of these ticks in China. The SFTSV susceptibility of parthenogenetic females was similar to that of bisexual females under laboratory conditions. These results suggest that parthenogenetic Asian longhorned ticks, probably transported by migratory birds, play a major role in the rapid spread of SFTSV.

Comparative hologenomics of two Ixodes scapularis tick populations in New Jersey.

Tick-borne diseases, such as those transmitted by the blacklegged tick Ixodes scapularis, are a significant and growing public health problem in the US. There is mounting evidence that co-occurring non-pathogenic microbes can also impact tick-borne disease transmission. Shotgun metagenome sequencing enables sampling of the complete tick hologenome-the collective genomes of the tick and all of the microbial species contained therein, whether pathogenic, commensal or symbiotic. This approach simultaneously uncovers taxonomic composition and allows the detection of intraspecific genetic variation, making it a useful tool to compare spatial differences across tick populations. We evaluated this approach by comparing hologenome data from two tick samples (N = 6 ticks per location) collected at a relatively fine spatial scale, approximately 23 km apart, within a single US county. Several intriguing variants in the data between the two sites were detected, including polymorphisms in both in the tick's own mitochondrial DNA and that of a rickettsial endosymbiont. The two samples were broadly similar in terms of the microbial species present, including multiple known tick-borne pathogens (Borrelia burgdorferi, Babesia microti, and Anaplasma phagocytophilum), filarial nematodes, and Wolbachia and Babesia species. We assembled the complete genome of the rickettsial endosymbiont (most likely Rickettsia buchneri) from both populations. Our results provide further evidence for the use of shotgun metagenome sequencing as a tool to compare tick hologenomes and differentiate tick populations across localized spatial scales.

Amoeba Genome Reveals Dominant Host Contribution to Plastid Endosymbiosis.

Eukaryotic photosynthetic organelles, plastids, are the powerhouses of many aquatic and terrestrial ecosystems. The canonical plastid in algae and plants originated >1 Ga and therefore offers limited insights into the initial stages of organelle evolution. To address this issue, we focus here on the photosynthetic amoeba Paulinella micropora strain KR01 (hereafter, KR01) that underwent a more recent (∼124 Ma) primary endosymbiosis, resulting in a photosynthetic organelle termed the chromatophore. Analysis of genomic and transcriptomic data resulted in a high-quality draft assembly of size 707 Mb and 32,361 predicted gene models. A total of 291 chromatophore-targeted proteins were predicted in silico, 208 of which comprise the ancestral organelle proteome in photosynthetic Paulinella species with functions, among others, in nucleotide metabolism and oxidative stress response. Gene coexpression analysis identified networks containing known high light stress response genes as well as a variety of genes of unknown function ("dark" genes). We characterized diurnally rhythmic genes in this species and found that over 49% are dark. It was recently hypothesized that large double-stranded DNA viruses may have driven gene transfer to the nucleus in Paulinella and facilitated endosymbiosis. Our analyses do not support this idea, but rather suggest that these viruses in the KR01 and closely related P. micropora MYN1 genomes resulted from a more recent invasion.

Leveraging the Expertise of the New Jersey Mosquito Control Community to Jump Start Standardized Tick Surveillance.

Despite the rising incidence of tick-borne diseases (TBD) in the northeastern United States (US), information and expertise needed to assess risk, inform the public and respond proactively is highly variable across states. Standardized and well-designed tick surveillance by trained personnel can facilitate the development of useful risk maps and help target resources, but requires nontrivial start-up costs. To address this challenge, we tested whether existing personnel in New Jersey's 21 county mosquito control agencies could be trained and interested to participate in a one-day collection of American dog ticks (Dermacentor variabilis), a presumably widespread species never before surveyed in this state. A workshop was held offering training in basic tick biology, identification, and standard operating procedures (SOPs) for surveillance, followed by a one-day simultaneous collection of D. variabilis across the state (the "NJ Tick Blitz"). In total, 498 D. variabilis were collected from 21 counties and follow-up participant surveys demonstrated an increase in knowledge and interest in ticks: 41.7% of respondents reported collecting ticks outside the Tick Blitz. We hope that the success of this initiative may provide a template for researchers and officials in other states with tick-borne disease concerns to obtain baseline tick surveillance data by training and partnering with existing personnel.

The genomes of polyextremophilic cyanidiales contain 1% horizontally transferred genes with diverse adaptive functions.

The role and extent of horizontal gene transfer (HGT) in eukaryotes are hotly disputed topics that impact our understanding of the origin of metabolic processes and the role of organelles in cellular evolution. We addressed this issue by analyzing 10 novel Cyanidiales genomes and determined that 1% of their gene inventory is HGT-derived. Numerous HGT candidates share a close phylogenetic relationship with prokaryotes that live in similar habitats as the Cyanidiales and encode functions related to polyextremophily. HGT candidates differ from native genes in GC-content, number of splice sites, and gene expression. HGT candidates are more prone to loss, which may explain the absence of a eukaryotic pan-genome. Therefore, the lack of a pan-genome and cumulative effects fail to provide substantive arguments against our hypothesis of recurring HGT followed by differential loss in eukaryotes. The maintenance of 1% HGTs, even under selection for genome reduction, underlines the importance of non-endosymbiosis related foreign gene acquisition.

Analysis of an improved Cyanophora paradoxa genome assembly.

Glaucophyta are members of the Archaeplastida, the founding group of photosynthetic eukaryotes that also includes red algae (Rhodophyta), green algae, and plants (Viridiplantae). Here we present a high-quality assembly, built using long-read sequences, of the ca. 100 Mb nuclear genome of the model glaucophyte Cyanophora paradoxa. We also conducted a quick-freeze deep-etch electron microscopy (QFDEEM) analysis of C. paradoxa cells to investigate glaucophyte morphology in comparison to other organisms. Using the genome data, we generated a resolved 115-taxon eukaryotic tree of life that includes a well-supported, monophyletic Archaeplastida. Analysis of muroplast peptidoglycan (PG) ultrastructure using QFDEEM shows that PG is most dense at the cleavage-furrow. Analysis of the chlamydial contribution to glaucophytes and other Archaeplastida shows that these foreign sequences likely played a key role in anaerobic glycolysis in primordial algae to alleviate ATP starvation under night-time hypoxia. The robust genome assembly of C. paradoxa significantly advances knowledge about this model species and provides a reference for exploring the panoply of traits associated with the anciently diverged glaucophyte lineage.

Genome analysis of the rice coral Montipora capitata.

Corals comprise a biomineralizing cnidarian, dinoflagellate algal symbionts, and associated microbiome of prokaryotes and viruses. Ongoing efforts to conserve coral reefs by identifying the major stress response pathways and thereby laying the foundation to select resistant genotypes rely on a robust genomic foundation. Here we generated and analyzed a high quality long-read based ~886 Mbp nuclear genome assembly and transcriptome data from the dominant rice coral, Montipora capitata from Hawai'i. Our work provides insights into the architecture of coral genomes and shows how they differ in size and gene inventory, putatively due to population size variation. We describe a recent example of foreign gene acquisition via a bacterial gene transfer agent and illustrate the major pathways of stress response that can be used to predict regulatory components of the transcriptional networks in M. capitata. These genomic resources provide insights into the adaptive potential of these sessile, long-lived species in both natural and human influenced environments and facilitate functional and population genomic studies aimed at Hawaiian reef restoration and conservation.

Active Host Response to Algal Symbionts in the Sea Slug Elysia chlorotica.

Sacoglossan sea slugs offer fascinating systems to study the onset and persistence of algal-plastid symbioses. Elysia chlorotica is particularly noteworthy because it can survive for months, relying solely on energy produced by ingested plastids of the stramenopile alga Vaucheria litorea that are sequestered in cells lining its digestive diverticula. How this animal can maintain the actively photosynthesizing organelles without replenishment of proteins from the lost algal nucleus remains unknown. Here, we used RNA-Seq analysis to test the idea that plastid sequestration leaves a significant signature on host gene expression during E. chlorotica development. Our results support this hypothesis and show that upon exposure to and ingestion of V. litorea plastids, genes involved in microbe-associated molecular patterns and oxidative stress-response mechanisms are significantly up-regulated. Interestingly, our results with E. chlorotica mirror those found with corals that maintain dinoflagellates as intact cells in symbiosomes, suggesting parallels between these animal-algal symbiotic interactions.

Genome wide analysis of the transition to pathogenic lifestyles in Magnaporthales fungi.

The rice blast fungus Pyricularia oryzae (syn. Magnaporthe oryzae, Magnaporthe grisea), a member of the order Magnaporthales in the class Sordariomycetes, is an important plant pathogen and a model species for studying pathogen infection and plant-fungal interaction. In this study, we generated genome sequence data from five additional Magnaporthales fungi including non-pathogenic species, and performed comparative genome analysis of a total of 13 fungal species in the class Sordariomycetes to understand the evolutionary history of the Magnaporthales and of fungal pathogenesis. Our results suggest that the Magnaporthales diverged ca. 31 millon years ago from other Sordariomycetes, with the phytopathogenic blast clade diverging ca. 21 million years ago. Little evidence of inter-phylum horizontal gene transfer (HGT) was detected in Magnaporthales. In contrast, many genes underwent positive selection in this order and the majority of these sequences are clade-specific. The blast clade genomes contain more secretome and avirulence effector genes, which likely play key roles in the interaction between Pyricularia species and their plant hosts. Finally, analysis of transposable elements (TE) showed differing proportions of TE classes among Magnaporthales genomes, suggesting that species-specific patterns may hold clues to the history of host/environmental adaptation in these fungi.

Robust Dinoflagellata phylogeny inferred from public transcriptome databases.

Dinoflagellates are dominant members of the plankton and play key roles in ocean ecosystems as primary producers, predators, parasites, coral photobionts, and causative agents of algal blooms that produce toxins harmful to humans and commercial fisheries. These unicellular protists exhibit remarkable trophic and morphological diversity and include species with some of the largest reported nuclear genomes. Despite their high ecological and economic importance, comprehensive genome (or transcriptome) based dinoflagellate trees of life are few in number. To address this issue, we used recently generated public sequencing data, including from the Moore Microbial Eukaryote Transcriptome Sequencing Project, to identify dinoflagellate-specific ortholog groups. These orthologs were combined to create a broadly sampled and highly resolved phylogeny of dinoflagellates. Our results emphasize the scope and utility of public sequencing databases in creating broad and robust phylogenies for large and complex taxonomic lineages, while also providing unique insights into the evolution of thecate dinoflagellates.

Impact of light intensity and quality on chromatophore and nuclear gene expression in Paulinella chromatophora, an amoeba with nascent photosynthetic organelles.

Plastid evolution has been attributed to a single primary endosymbiotic event that occurred about 1.6 billion years ago (BYA) in which a cyanobacterium was engulfed and retained by a eukaryotic cell, although early steps in plastid integration are poorly understood. The photosynthetic amoeba Paulinella chromatophora represents a unique model for the study of plastid evolution because it contains cyanobacterium-derived photosynthetic organelles termed 'chromatophores' that originated relatively recently (0.09-0.14 BYA). The chromatophore genome is about a third the size of the genome of closely related cyanobacteria, but 10-fold larger than most plastid genomes. Several genes have been transferred from the chromatophore genome to the host nuclear genome through endosymbiotic gene transfer (EGT). Some EGT-derived proteins could be imported into chromatophores for function. Two photosynthesis-related genes (psaI and csos4A) are encoded by both the nuclear and chromatophore genomes, suggesting that EGT in Paulinella chromatophora is ongoing. Many EGT-derived genes encode proteins that function in photosynthesis and photoprotection, including an expanded family of high-light-inducible (ncHLI) proteins. Cyanobacterial hli genes are high-light induced and required for cell viability under excess light. We examined the impact of light on Paulinella chromatophora and found that this organism is light sensitive and lacks light-induced transcriptional regulation of chromatophore genes and most EGT-derived nuclear genes. However, several ncHLI genes have reestablished light-dependent regulation, which appears analogous to what is observed in cyanobacteria. We postulate that expansion of the ncHLI gene family and its regulation may reflect the light/oxidative stress experienced by Paulinella chromatophora as a consequence of the as yet incomplete integration of host and chromatophore metabolisms.

Gene transfers from diverse bacteria compensate for reductive genome evolution in the chromatophore of Paulinella chromatophora.

Plastids, the photosynthetic organelles, originated >1 billion y ago via the endosymbiosis of a cyanobacterium. The resulting proliferation of primary producers fundamentally changed global ecology. Endosymbiotic gene transfer (EGT) from the intracellular cyanobacterium to the nucleus is widely recognized as a critical factor in the evolution of photosynthetic eukaryotes. The contribution of horizontal gene transfers (HGTs) from other bacteria to plastid establishment remains more controversial. A novel perspective on this issue is provided by the amoeba Paulinella chromatophora, which contains photosynthetic organelles (chromatophores) that are only 60-200 million years old. Chromatophore genome reduction entailed the loss of many biosynthetic pathways including those for numerous amino acids and cofactors. How the host cell compensates for these losses remains unknown, because the presence of bacteria in all available P. chromatophora cultures excluded elucidation of the full metabolic capacity and occurrence of HGT in this species. Here we generated a high-quality transcriptome and draft genome assembly from the first bacteria-free P. chromatophora culture to deduce rules that govern organelle integration into cellular metabolism. Our analyses revealed that nuclear and chromatophore gene inventories provide highly complementary functions. At least 229 nuclear genes were acquired via HGT from various bacteria, of which only 25% putatively arose through EGT from the chromatophore genome. Many HGT-derived bacterial genes encode proteins that fill gaps in critical chromatophore pathways/processes. Our results demonstrate a dominant role for HGT in compensating for organelle genome reduction and suggest that phagotrophy may be a major driver of HGT.