Search

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.

Transcriptome and genome data from twenty stony coral species and a selection of reference bilaterians were studied to elucidate coral evolutionary history. We identified genes that encode the proteins responsible for the precipitation and aggregation of the aragonite skeleton on which the organisms live, and revealed a network of environmental sensors that coordinate responses of the host animals to temperature, light, and pH. Furthermore, we describe a variety of stress-related pathways, including apoptotic pathways that allow the host animals to detoxify reactive oxygen and nitrogen species that are generated by their intracellular photosynthetic symbionts, and determine the fate of corals under environmental stress. Some of these genes arose through horizontal gene transfer and comprise at least 0.2% of the animal gene inventory. Our analysis elucidates the evolutionary strategies that have allowed symbiotic corals to adapt and thrive for hundreds of millions of years.

The tribe Aedini (Family Culicidae) contains approximately one-quarter of the known species of mosquitoes, including vectors of deadly or debilitating disease agents. This tribe contains the genus Aedes, which is one of the three most familiar genera of mosquitoes. During the past decade, Aedini has been the focus of a series of extensive morphology-based phylogenetic studies published by Reinert, Harbach, and Kitching (RH&K). Those authors created 74 new, elevated or resurrected genera from what had been the single genus Aedes, almost tripling the number of genera in the entire family Culicidae. The proposed classification is based on subjective assessments of the "number and nature of the characters that support the branches" subtending particular monophyletic groups in the results of cladistic analyses of a large set of morphological characters of representative species. To gauge the stability of RH&K's generic groupings we reanalyzed their data with unweighted parsimony jackknife and maximum-parsimony analyses, with and without ordering 14 of the characters as in RH&K. We found that their phylogeny was largely weakly supported and their taxonomic rankings failed priority and other useful taxon-naming criteria. Consequently, we propose simplified aedine generic designations that 1) restore a classification system that is useful for the operational community; 2) enhance the ability of taxonomists to accurately place new species into genera; 3) maintain the progress toward a natural classification based on monophyletic groups of species; and 4) correct the current classification system that is subject to instability as new species are described and existing species more thoroughly defined. We do not challenge the phylogenetic hypotheses generated by the above-mentioned series of morphological studies. However, we reduce the ranks of the genera and subgenera of RH&K to subgenera or informal species groups, respectively, to preserve stability as new data become available.

Improving feedstock is critical to facilitate the commercial utilization of algae, in particular in open pond systems where, due to the presence of competitors and pests, high algal growth rates and stress tolerance are beneficial. Here we raised laboratory cultures of the model alga Chlamydomonas reinhardtii under serial dilution to explore the potential of crop improvement using natural selection. The alga was evolved for 1,880 generations in liquid medium under continuous light (EL population). At the end of the experiment, EL cells had a growth rate that was 35% greater than the progenitor population (PL). The removal of acetate from the medium demonstrated that EL growth enhancement largely relied on efficient usage of this organic carbon source. Genome re-sequencing uncovered 1,937 polymorphic DNA regions in the EL population with 149 single nucleotide polymorphisms resulting in amino acid substitutions. Transcriptome analysis showed, in the EL population, significant up regulation of genes involved in protein synthesis, the cell cycle and cellular respiration, whereas the DNA repair pathway and photosynthesis were down regulated. Like other algae, EL cells accumulated neutral lipids under nitrogen depletion. Our work demonstrates transcriptome and genome-wide impacts of natural selection on algal cells and points to a useful strategy for strain improvement.

Abstract 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.

RNAi (RNA interference) relies on the production of small RNAs (sRNAs) from double-stranded RNA and comprises a major pathway in eukaryotes to restrict the propagation of selfish genetic elements. Amplification of the initial RNAi signal by generation of multiple secondary sRNAs from a targeted mRNA is catalyzed by RNA-dependent RNA polymerases (RdRPs). This phenomenon is known as transitivity and is particularly important in plants to limit the spread of viruses. Here we describe, using a genome-wide approach, the distribution of sRNAs in the glaucophyte alga Cyanophora paradoxa. C. paradoxa is a member of the supergroup Plantae (also known as Archaeplastida) that includes red algae, green algae, and plants. The ancient (>1 billion years ago) split of glaucophytes within Plantae suggests that C. paradoxa may be a useful model to learn about the early evolution of RNAi in the supergroup that ultimately gave rise to plants. Using next-generation sequencing and bioinformatic analyses we find that sRNAs in C. paradoxa are preferentially associated with mRNAs, including a large number of transcripts that encode proteins arising from different functional categories. This pattern of exonic sRNAs appears to be a general trend that affects a large fraction of mRNAs in the cell. In several cases we observe that sRNAs have a bias for a specific strand of the mRNA, including many instances of antisense predominance. The genome of C. paradoxa encodes four sequences that are homologous to RdRPs in Arabidopsis thaliana. We discuss the possibility that exonic sRNAs in the glaucophyte may be secondarily derived from mRNAs by the action of RdRPs. If this hypothesis is confirmed, then transitivity may have had an ancient origin in Plantae.

Summary: Temperature is increasing globally, and vector-borne diseases are particularly responsive to such increases. While it is known that temperature influences mosquito life history traits, transmission models have not historically considered population-specific effects of temperature. We assessed the interaction between Culex pipiens population and temperature in New York State (NYS) and utilized novel empirical data to inform predictive models of West Nile virus (WNV) transmission. Genetically and regionally distinct populations from NYS were reared at various temperatures, and life history traits were monitored and used to inform trait-based models. Variation in Cx. pipiens life history traits and population-dependent thermal responses account for a predicted 2.9°C difference in peak transmission that is reflected in regional differences in WNV prevalence. We additionally identified genetic signatures that may contribute to distinct thermal responses. Together, these data demonstrate how population variation contributes to significant geographic variability in arbovirus transmission with changing climates.

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.

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.

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.

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.

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.

Group II introns are closely linked to eukaryote evolution because nuclear spliceosomal introns and the small RNAs associated with the spliceosome are thought to trace their ancient origins to these mobile elements. Therefore, elucidating how group II introns move, and how they lose mobility can potentially shed light on fundamental aspects of eukaryote biology. To this end, we studied five strains of the unicellular red alga Porphyridium purpureum that surprisingly contain 42 group II introns in their plastid genomes. We focused on a subset of these introns that encode mobility-conferring intron-encoded proteins (IEPs) and found them to be distributed among the strains in a lineage-specific manner. The reverse transcriptase and maturase domains were present in all lineages but the DNA endonuclease domain was deleted in vertically inherited introns, demonstrating a key step in the loss of mobility. P. purpureum plastid intron RNAs had a classic group IIB secondary structure despite variability in the DIII and DVI domains. We report for the first time the presence of twintrons (introns-within-introns, derived from the same mobile element) in Rhodophyta. The P. purpureum IEPs and their mobile introns provide a valuable model for the study of mobile retroelements in eukaryotes and offer promise for biotechnological applications.

Culex pipiens, an invasive mosquito and vector of West Nile virus in the US, has two morphologically indistinguishable forms that differ dramatically in behavior and physiology. Cx. pipiens form pipiens is primarily a bird-feeding temperate mosquito, while the sub-tropical Cx. pipiens form molestus thrives in sewers and feeds on mammals. Because the feral form can diapause during the cold winters but the domestic form cannot, the two Cx. pipiens forms are allopatric in northern Europe and, although viable, hybrids are rare. Cx. pipiens form molestus has spread across all inhabited continents and hybrids of the two forms are common in the US. Here we elucidate the genes and gene families with the greatest divergence rates between these phenotypically diverged mosquito populations, and discuss them in light of their potential biological and ecological effects. After generating and assembling novel transcriptome data for each population, we performed pairwise tests for nonsynonymous divergence (Ka) of homologous coding sequences and examined gene ontology terms that were statistically over-represented in those sequences with the greatest divergence rates. We identified genes involved in digestion (serine endopeptidases), innate immunity (fibrinogens and α-macroglobulins), hemostasis (D7 salivary proteins), olfaction (odorant binding proteins) and chitin binding (peritrophic matrix proteins). By examining molecular divergence between closely related yet phenotypically divergent forms of the same species, our results provide insights into the identity of rapidly-evolving genes between incipient species. Additionally, we found that families of signal transducers, ATP synthases and transcription regulators remained identical at the amino acid level, thus constituting conserved components of the Cx. pipiens proteome. We provide a reference with which to gauge the divergence reported in this analysis by performing a comparison of transcriptome sequences from conspecific (yet allopatric) populations of another member of the Cx. pipiens complex, Cx. quinquefasciatus.