Strong and consistent effects of waterbird composition on HPAI H5 occurrences across Europe.

Since 2014, highly pathogenic avian influenza (HPAI) H5 viruses of clade 2.3.4.4 have been dominating the outbreaks across Europe, causing massive deaths among poultry and wild birds. However, the factors shaping these broad-scale outbreak patterns, especially those related to waterbird community composition, remain unclear. In particular, we do not know whether these risk factors differ from those of other H5 clades. Addressing this knowledge gap is important for predicting and preventing future HPAI outbreaks. Using extensive waterbird survey datasets from about 6883 sites, we here explored the effect of waterbird community composition on HPAI H5Nx (clade 2.3.4.4) spatial patterns in the 2016/2017 and 2020/2021 epidemics in Europe, and compared it with the 2005/2006 HPAI H5N1 (clade 2.2) epidemic. We showed that HPAI H5 occurrences in wild birds in the three epidemics were strongly associated with very similar waterbird community attributes, which suggested that, in nature, similar interspecific transmission processes operate between the HPAI H5 subtypes or clades. Importantly, community phylogenetic diversity consistently showed a negative association with H5 occurrence in all three epidemics, suggesting a dilution effect of phylogenetic diversity. In contrast, waterbird community variables showed much weaker associations with HPAI H5Nx occurrence in poultry. Our results demonstrate that models based on previous epidemics can predict future HPAI H5 patterns in wild birds, implying that it is important to include waterbird community factors in future HPAI studies to predict outbreaks and improve surveillance activities.

Chromosome-level genome assembly of navel orange cv. Gannanzao (Citrus sinensis Osbeck cv. Gannanzao).

Navel orange cv. Gannanzao is a variant of the navel orange cv. Newhall (Citrus sinensis Osbeck cv. Newhall) that exhibits an earlier maturation, making it commercially valuable. However, the mechanisms underlying its early maturation remain obscure. To address this question, we conducted genome sequencing and de novo assembly of navel orange cv. Gannanzao. The assembled genome sequence is 334.57 Mb in length with a GC content of 31.48%. It comprises 318 contigs (N50 = 3.23 Mb) and 187 scaffolds (N50 = 31.86 Mb). The Benchmarking Universal Single-Copy Orthologs test demonstrates 94.6% completeness. The annotation revealed 23,037 gene models, 164.95 Mb of repetitive sequences, and 2,554 noncoding RNAs. A comparative analysis identified 323 fruit ripening-related genes in navel orange cv. Gannanzao genome, while navel orange cv. Newhall genome contained 345 such genes. These genes were organized into 320 orthologous gene families, with 30.3% of them exhibiting differences in gene copy numbers between the 2 genomes. Additionally, we identified 15 fruit ripening-related genes that have undergone adaptive evolution, suggesting their potential role in advancing fruit maturation in navel orange cv. Gannanzao. Whole-genome sequencing and annotation of navel orange cv. Gannanzao provides a valuable resource to unravel the early maturation mechanism of citrus and enriches the genomic resources for citrus research.

Plants acquired mitochondrial linear plasmids horizontally from fungi likely during the conquest of land.

Mitochondrial linear plasmids have been sporadically reported in fungi and plants. Yet, much remains obscure about the diversity, distribution, and evolution of mitochondrial linear plasmids. Here, through phylogenomic analyses across 7,163 cellular organisms (including 991 plants), we find that mitochondrial linear plasmids are widely present in land plants and fungi. Phylogenetic analyses indicate that plants are likely to have acquired mitochondrial linear plasmids horizontally from fungi before or during the conquest of terrestrial environments by plants. Gene content analyses show that mitochondrial linear plasmids harbor a highly dynamic and promiscuous repertoire of genes. Our study refines the understanding of the origin and evolution of mitochondrial linear plasmids.

Functional traits explain waterbirds' host status, subtype richness, and community-level infection risk for avian influenza.

Species functional traits can influence pathogen transmission processes, and consequently affect species' host status, pathogen diversity, and community-level infection risk. We here investigated, for 143 European waterbird species, effects of functional traits on host status and pathogen diversity (subtype richness) for avian influenza virus at species level. We then explored the association between functional diversity and HPAI H5Nx occurrence at the community level for 2016/17 and 2021/22 epidemics in Europe. We found that both host status and subtype richness were shaped by several traits, such as diet guild and dispersal ability, and that the community-weighted means of these traits were also correlated with community-level risk of H5Nx occurrence. Moreover, functional divergence was negatively associated with H5Nx occurrence, indicating that functional diversity can reduce infection risk. Our findings highlight the value of integrating trait-based ecology into the framework of diversity-disease relationship, and provide new insights for HPAI prediction and prevention.

Genome mining shows that retroviruses are pervasively invading vertebrate genomes.

Endogenous retroviruses (ERVs) record past retroviral infections, providing molecular archives for interrogating the evolution of retroviruses and retrovirus-host interaction. However, the vast majority of ERVs are not active anymore due to various disruptive mutations, and ongoing retroviral invasion of vertebrate genomes has been rarely documented. Here we analyze genomics data from 2004 vertebrates for mining invading ERVs (ERVi). We find that at least 412 ERVi elements representing 217 viral operational taxonomic units are invading the genomes of 123 vertebrates, 18 of which have been assessed to be threatened species. Our results reveal an unexpected prevalence of ongoing retroviral invasion in vertebrates and expand the diversity of retroviruses recently circulating in the wild. We characterize the pattern and nature of ERVi in the historical and biogeographical context of their hosts, for instance, the generation of model organisms, sympatric speciation, and domestication. We suspect that these ERVi are relevant to conservation of threatened species, zoonoses in the wild, and emerging infectious diseases in humans.

Origin of the OAS-RNase L innate immune pathway before the rise of jawed vertebrates via molecular tinkering.

Discriminating self from nonself is fundamental to immunity. Yet, it remains largely elusive how the mechanisms of self and nonself discrimination originated. Sensing double-stranded RNA as nonself, the 2',5'-oligoadenylate synthetase (OAS)-ribonuclease L (RNase L) pathway represents a crucial component of innate immunity. Here, we combine phylogenomic and functional analyses to show that the functional OAS-RNase L pathway likely originated through tinkering with preexisting proteins before the rise of jawed vertebrates during or before the Silurian period (444 to 419 Mya). Multiple concerted losses of OAS and RNase L occurred during the evolution of jawed vertebrates, further supporting the ancient coupling between OAS and RNase L. Moreover, both OAS and RNase L genes evolved under episodic positive selection across jawed vertebrates, suggesting a long-running evolutionary arms race between the OAS-RNase L pathway and microbes. Our findings illuminate how an innate immune pathway originated via molecular tinkering.

Jawless vertebrates do not escape retrovirus infection.

On occasion, retroviruses infect the genome of germline cell, forming endogenous retroviruses (ERVs), which provide molecular fossils for studying the deep evolution of retroviruses. While ERVs have been extensively characterized in the genomes of jawed vertebrates, much remains contentious and unexplored about the diversity and evolution of ERVs within jawless vertebrates. Here, we report the discovery of a novel ERV lineage, designated EbuERVs, in the genome of a hagfish Eptatretus burgeri. Phylogenetic analyses show that EbuERVs pertain to epsilon-retroviruses and might have derived from cross-species transmission from jawed vertebrates. EbuERVs are estimated to have invaded in the hagfish genome at least tens of millions of years ago. Evolutionary dynamics analyses indicate that EbuERVs might have experienced one proliferation peak and have been not active in transposition anymore. However, some EbuERVs can transcribe in embryo and might serve as lncRNA. Overall, these findings expand the distribution of retroviruses from jawed vertebrates to jawless vertebrates.

A Long-Running Arms Race between APOBEC1 Genes and Retroviruses in Tetrapods.

Activation-induced cytidine deaminase/apolipoprotein B mRNA editing catalytic polypeptide-like (AID/APOBEC) proteins are cytosine deaminases implicated in diverse biological functions. APOBEC1 (A1) proteins have long been thought to regulate lipid metabolism, whereas the evolutionary significance of A1 proteins in antiviral defense remains largely obscure. Endogenous retroviruses (ERVs) document past retroviral infections and are ubiquitous within the vertebrate genomes. Here, we identify the A1 gene repertoire, characterize the A1-mediated mutation footprints in ERVs, and interrogate the evolutionary arms race between A1 genes and ERVs across vertebrate species. We find that A1 genes are widely present in tetrapods, recurrently amplified and lost in certain lineages, suggesting that A1 genes might have originated during the early evolution of tetrapods. A1-mediated mutation footprints can be detected in ERVs across tetrapods. Moreover, A1 genes appear to have experienced episodic positive selection in many tetrapod lineages. Taken together, we propose that a long-running arms race between A1 genes and retroviruses might have persisted throughout the evolutionary course of tetrapods. IMPORTANCE APOBEC3 (A3) genes have been thought to function in defense against retroviruses, whereas the evolutionary significance of A1 proteins in antiviral defense remains largely obscure. In this study, we identify the A1 gene repertoire, characterize the A1-mediated mutation footprints in endogenous retroviruses (ERVs), and explore the evolutionary arms race between A1 genes and ERVs across vertebrate species. We found A1 proteins originated during the early evolution of tetrapods, and detected the footprints of A1-induced hypermutations in retroviral fossils. A1 genes appear to have experienced pervasive positive selection in tetrapods. Our study indicates a long-running arms race between A1 genes and retroviruses taking place throughout the evolutionary course of tetrapods.

Two plant NLR proteins confer strain-specific resistance conditioned by an effector from Pseudomonas syringae pv. actinidiae.

Pseudomonas syringae pv. actinidiae (Psa) causes bacterial canker, a devastating disease threatening the Actinidia fruit industry. In a search for non-host resistance genes against Psa, we find that the nucleotide-binding leucine-rich repeat receptor (NLR) protein ZAR1 from both Arabidopsis and Nicotiana benthamiana (Nb) recognizes HopZ5 and triggers cell death. The recognition requires ZED1 in Arabidopsis and JIM2 in Nb plants, which are members of the ZRK pseudokinases and known components of the ZAR1 resistosome. Surprisingly, Arabidopsis ZAR1 and RPM1, another NLR known to recognize HopZ5, confer disease resistance to HopZ5 in a strain-specific manner. Thus, ZAR1, but not RPM1, is solely required for resistance to P. s. maculicola ES4326 (Psm) carrying hopZ5, whereas RPM1 is primarily required for resistance to P. s. tomato DC3000 (Pst) carrying hopZ5. Furthermore, the ZAR1-mediated resistance to Psm hopZ5 in Arabidopsis is insensitive to SOBER1, which encodes a deacetylase known to suppress the RPM1-mediated resistance to Pst hopZ5. In addition, hopZ5 enhances P. syringae virulence in the absence of ZAR1 or RPM1 and that SOBER1 abolishes such virulence function. Together the study suggests that ZAR1 may be used for improving Psa resistance in Actinidia and uncovers previously unknown complexity of effector-triggered immunity and effector-triggered virulence.

The genome of an apodid holothuroid (Chiridota heheva) provides insights into its adaptation to a deep-sea reducing environment.

Cold seeps and hydrothermal vents are deep-sea reducing environments that are characterized by lacking oxygen and photosynthesis-derived nutrients. Most animals acquire nutrition in cold seeps or hydrothermal vents by maintaining epi- or endosymbiotic relationship with chemoautotrophic microorganisms. Although several seep- and vent-dwelling animals hosting symbiotic microbes have been well-studied, the genomic basis of adaptation to deep-sea reducing environment in nonsymbiotic animals is still lacking. Here, we report a high-quality genome of Chiridota heheva Pawson & Vance, 2004, which thrives by extracting organic components from sediment detritus and suspended material, as a reference for nonsymbiotic animal's adaptation to deep-sea reducing environments. The expansion of the aerolysin-like protein family in C. heheva compared with other echinoderms might be involved in the disintegration of microbes during digestion. Moreover, several hypoxia-related genes (Pyruvate Kinase M2, PKM2; Phospholysine Phosphohistidine Inorganic Pyrophosphate Phosphatase, LHPP; Poly(A)-specific Ribonuclease Subunit PAN2, PAN2; and Ribosomal RNA Processing 9, RRP9) were subject to positive selection in the genome of C. heheva, which contributes to their adaptation to hypoxic environments.

The diversity and evolution of retroviruses: Perspectives from viral "fossils".

Retroviruses exclusively infect vertebrates, causing a variety of diseases. The replication of retroviruses requires reverse transcription and integration into host genomes. When infecting germline cells, retroviruses become inherited vertically, forming endogenous retroviruses (ERVs). ERVs document past viral infections, providing molecular fossils for studying the evolutionary history of retroviruses. In this review, we summarize the recent advances in understanding the diversity and evolution of retroviruses from the perspectives of viral fossils, and discuss the effects of ERVs on the evolution of host biology.

A Missing Link between Retrotransposons and Retroviruses.

The origin and deep evolution of retroviruses remain largely unclear. It has been proposed that retroviruses might have originated from a Ty3/Gypsy retrotransposon, but all known Ty3/Gypsy retrotransposons are only distantly related to retroviruses. Retroviruses and some plant Athila/Tat elements (within Ty3/Gypsy retrotransposons) independently evolved a dual RNase H domain and an env/env-like gene. Here, we reported the discovery of a novel lineage of retrotransposons, designated Odin retrotransposons, in the genomes of eight sea anemones (order Actinaria) within the Cnidaria phylum. Odin retrotransposons exhibited unique genome features, encoding a dual RNase H domain (like retroviruses) but no env gene (like most Ty3/Gypsy retrotransposons). Phylogenetic analyses based on reverse transcriptase showed that Odin retrotransposons formed a sister group to lokiretroviruses, and lokiretroviruses and Odin retrotransposons together were sister to canonical retroviruses. Moreover, phylogenetic analyses based on RNase H and integrase also supported the hypothesis that Odin retrotransposons were sisters to lokiretroviruses. Lokiretroviruses and canonical retroviruses did not form a monophyletic group, indicating that lokiretroviruses and canonical retroviruses might represent two distinct virus families. Taken together, the discovery of Odin retrotransposons narrowed down the evolutionary gaps between retrotransposons and canonical retroviruses and lokiretroviruses. IMPORTANCE The origin of retroviruses remains largely unclear. In this study, we discovered a novel retrotransposon lineage, Odin retrotransposons, within the genomes of sea anemones (order Actinaria). In contrast to retroviruses and most retrotransposons, Odin retrotransposons encode a dual RNase H domain but no env gene. Phylogenetic analyses showed that Odin retrotransposons were sisters to lokiretroviruses, and lokiretroviruses and Odin retrotransposons were sisters to retroviruses, establishing an evolutionary framework to decipher the origin of retroviruses (canonical retroviruses and lokiretroviruses). Our results provided insights into the diversity and deep evolution of LTR retrotransposons closely related to retroviruses.

The origin and evolution of a plant resistosome.

The nucleotide-binding, leucine-rich receptor (NLR) protein HOPZ-ACTIVATED RESISTANCE 1 (ZAR1), an immune receptor, interacts with HOPZ-ETI-DEFICIENT 1 (ZED1)-related kinases (ZRKs) and AVRPPHB SUSCEPTIBLE 1-like proteins to form a pentameric resistosome, triggering immune responses. Here, we show that ZAR1 emerged through gene duplication and that ZRKs were derived from the cell surface immune receptors wall-associated protein kinases (WAKs) through the loss of the extracellular domain before the split of eudicots and monocots during the Jurassic period. Many angiosperm ZAR1 orthologs, but not ZAR1 paralogs, are capable of oligomerization in the presence of AtZRKs and triggering cell death, suggesting that the functional ZAR1 resistosome might have originated during the early evolution of angiosperms. Surprisingly, inter-specific pairing of ZAR1 and AtZRKs sometimes results in the formation of a resistosome in the absence of pathogen stimulation, suggesting within-species compatibility between ZAR1 and ZRKs as a result of co-evolution. Numerous concerted losses of ZAR1 and ZRKs occurred in angiosperms, further supporting the ancient co-evolution between ZAR1 and ZRKs. Our findings provide insights into the origin of new plant immune surveillance networks.

Segregating Complete Tf2 Elements Are Largely Neutral in Fission Yeast.

Transposable elements (TEs) comprise a large proportion of the eukaryote genomes. Yet it remains poorly understood how TEs influence the fitness of the hosts carrying them. Here, we empirically test the impact of TEs on the host fitness in the fission yeast Schizosaccharomyces pombe. We find that two families of TEs (Tf1 and Tf2 elements), both of which belong to long terminal repeat retrotransposons, are highly polymorphic among individual S. pombe strains. Only 13 complete Tf2 elements are identified in S. pombe laboratory strain 972. These 13 Tf2 elements integrated into host genomes in very recent time and are segregating within the S. pombe population. Through knocking out each of the 13 Tf2 elements in S. pombe strain 972, we find Tf2 knockout does not affect the host fitness, and Tf2 elements do not alter the expression of nearby genes. Challenged by diverse forms of stress, the Tf2 knockout strains do not exhibit different growth rates from wild-type strain. Together, we conclude that segregating complete Tf2 elements insertions are largely neutral to host fitness in the fission yeast. Our study provides genome-wide empirical support for the selfish nature of TEs in fission yeast.

Ancient Adaptative Evolution of ACE2 in East Asians.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been posing an unprecedented challenge to global public health. SARS-CoV-2 and several other coronaviruses utilize angiotensin-converting enzyme 2 (ACE2) as their entry receptors. The ACE2 gene has been found to experience episodic positive selection across mammals. However, much remains unknown about how the ACE2 gene evolved in human populations. Here, we use population genetics approaches to investigate the evolution of the ACE2 gene in 26 human populations sampled globally. We find the ACE2 gene exhibits an extremely low nucleotide diversity in the East Asian populations. Strong signals of selective sweep are detected in the East Asian populations, but not in the other human populations. The selective sweep in ACE2 is estimated to begin in East Asian populations ∼23,600 years ago. Our study provides novel insights into the evolution of the ACE2 gene within human populations.

Viromics unveils extraordinary genetic diversity of the family Closteroviridae in wild citrus.

Our knowledge of citrus viruses is largely skewed toward virus pathology in cultivated orchards. Little is known about the virus diversity in wild citrus species. Here, we used a metatranscriptomics approach to characterize the virus diversity in a wild citrus habitat within the proposed center of the origin of citrus plants. We discovered a total of 44 virus isolates that could be classified into species Citrus tristeza virus and putative species citrus associated ampelovirus 1, citrus associated ampelovirus 2, and citrus virus B within the family Closteroviridae, providing important information to explore the factors facilitating outbreaks of citrus viruses and the evolutionary history of the family Closteroviridae. We found that frequent horizontal gene transfer, gene duplication, and alteration of expression strategy have shaped the genome complexity and diversification of the family Closteroviridae. Recombination frequently occurred among distinct Closteroviridae members, thereby facilitating the evolution of Closteroviridae. Given the potential emergence of similar wild-citrus-originated novel viruses as pathogens, the need for surveillance of their pathogenic and epidemiological characteristics is of utmost priority for global citrus production.

Molecular fossils illuminate the evolution of retroviruses following a macroevolutionary transition from land to water.

The ancestor of cetaceans underwent a macroevolutionary transition from land to water early in the Eocene Period >50 million years ago. However, little is known about how diverse retroviruses evolved during this shift from terrestrial to aquatic environments. Did retroviruses transition into water accompanying their hosts? Did retroviruses infect cetaceans through cross-species transmission after cetaceans invaded the aquatic environments? Endogenous retroviruses (ERVs) provide important molecular fossils for tracing the evolution of retroviruses during this macroevolutionary transition. Here, we use a phylogenomic approach to study the origin and evolution of ERVs in cetaceans. We identify a total of 8,724 ERVs within the genomes of 25 cetaceans, and phylogenetic analyses suggest these ERVs cluster into 315 independent lineages, each of which represents one or more independent endogenization events. We find that cetacean ERVs originated through two possible routes. 298 ERV lineages may derive from retrovirus endogenization that occurred before or during the transition from land to water of cetaceans, and most of these cetacean ERVs were reaching evolutionary dead-ends. 17 ERV lineages are likely to arise from independent retrovirus endogenization events that occurred after the split of mysticetes and odontocetes, indicating that diverse retroviruses infected cetaceans through cross-species transmission from non-cetacean mammals after the transition to aquatic life of cetaceans. Both integration time and synteny analyses support the recent or ongoing activity of multiple retroviral lineages in cetaceans, some of which proliferated into hundreds of copies within the host genomes. Although ERVs only recorded a proportion of past retroviral infections, our findings illuminate the complex evolution of retroviruses during one of the most marked macroevolutionary transitions in vertebrate history.