Structural Modeling of the SARS-CoV-2 Spike/Human ACE2 Complex Interface can Identify High-Affinity Variants Associated with Increased Transmissibility.

The COVID-19 pandemic has triggered concerns about the emergence of more infectious and pathogenic viral strains. As a public health measure, efficient screening methods are needed to determine the functional effects of new sequence variants. Here we show that structural modeling of SARS-CoV-2 Spike protein binding to the human ACE2 receptor, the first step in host-cell entry, predicts many novel variant combinations with enhanced binding affinities. By focusing on natural variants at the Spike-hACE2 interface and assessing over 700 mutant complexes, our analysis reveals that high-affinity Spike mutations (including N440K, S443A, G476S, E484R, G502P) tend to cluster near known human ACE2 recognition sites (K31 and K353). These Spike regions are structurally flexible, allowing certain mutations to optimize interface interaction energies. Although most human ACE2 variants tend to weaken binding affinity, they can interact with Spike mutations to generate high-affinity double mutant complexes, suggesting variation in individual susceptibility to infection. Applying structural analysis to highly transmissible variants, we find that circulating point mutations S477N, E484K and N501Y form high-affinity complexes (~40% more than wild-type). By combining predicted affinities and available antibody escape data, we show that fast-spreading viral variants exploit combinatorial mutations possessing both enhanced affinity and antibody resistance, including S477N/E484K, E484K/N501Y and K417T/E484K/N501Y. Thus, three-dimensional modeling of the Spike/hACE2 complex predicts changes in structure and binding affinity that correlate with transmissibility and therefore can help inform future intervention strategies.

Nearly Complete Genome Sequence of Brugia malayi Strain FR3.

Lymphatic filariasis affects ∼120 million people and can result in elephantiasis and hydrocele. Here, we report the nearly complete genome sequence of the best-studied causative agent of lymphatic filariasis, Brugia malayi The assembly contains four autosomes, an X chromosome, and only eight gaps but lacks a contiguous sequence for the known Y chromosome.

Sex chromosome evolution in parasitic nematodes of humans.

Sex determination mechanisms often differ even between related species yet the evolution of sex chromosomes remains poorly understood in all but a few model organisms. Some nematodes such as Caenorhabditis elegans have an XO sex determination system while others, such as the filarial parasite Brugia malayi, have an XY mechanism. We present a complete B. malayi genome assembly and define Nigon elements shared with C. elegans, which we then map to the genomes of other filarial species and more distantly related nematodes. We find a remarkable plasticity in sex chromosome evolution with several distinct cases of neo-X and neo-Y formation, X-added regions, and conversion of autosomes to sex chromosomes from which we propose a model of chromosome evolution across different nematode clades. The phylum Nematoda offers a new and innovative system for gaining a deeper understanding of sex chromosome evolution.

Defining Brugia malayi and Wolbachia symbiosis by stage-specific dual RNA-seq.

BACKGROUND: Filarial nematodes currently infect up to 54 million people worldwide, with millions more at risk for infection, representing the leading cause of disability in the developing world. Brugia malayi is one of the causative agents of lymphatic filariasis and remains the only human filarial parasite that can be maintained in small laboratory animals. Many filarial nematode species, including B. malayi, carry an obligate endosymbiont, the alpha-proteobacteria Wolbachia, which can be eliminated through antibiotic treatment. Elimination of the endosymbiont interferes with development, reproduction, and survival of the worms within the mamalian host, a clear indicator that the Wolbachia are crucial for survival of the parasite. Little is understood about the mechanism underlying this symbiosis. METHODOLOGY/ PRINCIPLE FINDINGS: To better understand the molecular interplay between these two organisms we profiled the transcriptomes of B. malayi and Wolbachia by dual RNA-seq across the life cycle of the parasite. This helped identify functional pathways involved in this essential symbiotic relationship provided by the co-expression of nematode and bacterial genes. We have identified significant stage-specific and gender-specific differential expression in Wolbachia during the nematode's development. For example, during female worm development we find that Wolbachia upregulate genes involved in ATP production and purine biosynthesis, as well as genes involved in the oxidative stress response. CONCLUSIONS/ SIGNIFICANCE: This global transcriptional analysis has highlighted specific pathways to which both Wolbachia and B. malayi contribute concurrently over the life cycle of the parasite, paving the way for the development of novel intervention strategies.

Quantifying influenza virus diversity and transmission in humans.

Influenza A virus is characterized by high genetic diversity. However, most of what is known about influenza evolution has come from consensus sequences sampled at the epidemiological scale that only represent the dominant virus lineage within each infected host. Less is known about the extent of within-host virus diversity and what proportion of this diversity is transmitted between individuals. To characterize virus variants that achieve sustainable transmission in new hosts, we examined within-host virus genetic diversity in household donor-recipient pairs from the first wave of the 2009 H1N1 pandemic when seasonal H3N2 was co-circulating. Although the same variants were found in multiple members of the community, the relative frequencies of variants fluctuated, with patterns of genetic variation more similar within than between households. We estimated the effective population size of influenza A virus across donor-recipient pairs to be approximately 100-200 contributing members, which enabled the transmission of multiple lineages, including antigenic variants.

Intrahost dynamics of antiviral resistance in influenza A virus reflect complex patterns of segment linkage, reassortment, and natural selection.

UNLABELLED: Resistance following antiviral therapy is commonly observed in human influenza viruses. Although this evolutionary process is initiated within individual hosts, little is known about the pattern, dynamics, and drivers of antiviral resistance at this scale, including the role played by reassortment. In addition, the short duration of human influenza virus infections limits the available time window in which to examine intrahost evolution. Using single-molecule sequencing, we mapped, in detail, the mutational spectrum of an H3N2 influenza A virus population sampled from an immunocompromised patient who shed virus over a 21-month period. In this unique natural experiment, we were able to document the complex dynamics underlying the evolution of antiviral resistance. Individual resistance mutations appeared weeks before they became dominant, evolved independently on cocirculating lineages, led to a genome-wide reduction in genetic diversity through a selective sweep, and were placed into new combinations by reassortment. Notably, despite frequent reassortment, phylogenetic analysis also provided evidence for specific patterns of segment linkage, with a strong association between the hemagglutinin (HA)- and matrix (M)-encoding segments that matches that previously observed at the epidemiological scale. In sum, we were able to reveal, for the first time, the complex interaction between multiple evolutionary processes as they occur within an individual host. IMPORTANCE: Understanding the evolutionary forces that shape the genetic diversity of influenza virus is crucial for predicting the emergence of drug-resistant strains but remains challenging because multiple processes occur concurrently. We characterized the evolution of antiviral resistance in a single persistent influenza virus infection, representing the first case in which reassortment and the complex patterns of drug resistance emergence and evolution have been determined within an individual host. Deep-sequence data from multiple time points revealed that the evolution of antiviral resistance reflects a combination of frequent mutation, natural selection, and a complex pattern of segment linkage and reassortment. In sum, these data show how immunocompromised hosts may help reveal the drivers of strain emergence.

Direct regulation of diurnal Drd3 expression and cocaine reward by NPAS2.

BACKGROUND: Circadian gene disruptions are associated with the development of psychiatric disorders, including addiction. However, the mechanisms by which circadian genes regulate reward remain poorly understood. METHODS: We used mice with a mutation in Npas2 and adeno-associated virus-short hairpin RNA mediated knockdown of Npas2 and Clock in the nucleus accumbens (NAc). We performed conditioned place preference assays. We utilized cell sorting quantitative real-time polymerase chain reaction, and chromatin immunoprecipitation followed by deep sequencing. RESULTS: Npas2 mutants exhibit decreased sensitivity to cocaine reward, which is recapitulated with a knockdown of neuronal PAS domain protein 2 (NPAS2) specifically in the NAc, demonstrating the importance of NPAS2 in this region. Interestingly, reducing circadian locomotor output cycles kaput (CLOCK) (a homologue of NPAS2) in the NAc had no effect, suggesting an important distinction in NPAS2 and CLOCK function. Furthermore, we found that NPAS2 expression is restricted to Drd1 expressing neurons while CLOCK is ubiquitous. Moreover, NPAS2 and CLOCK have distinct temporal patterns of DNA binding, and we identified novel and unique binding sites for each protein. We identified the Drd3 dopamine receptor as a direct transcriptional target of NPAS2 and found that NPAS2 knockdown in the NAc disrupts its diurnal rhythm in expression. Chronic cocaine treatment likewise disrupts the normal rhythm in Npas2 and Drd3 expression in the NAc, which may underlie behavioral plasticity in response to cocaine. CONCLUSIONS: Together, these findings identify an important role for the circadian protein, NPAS2, in the NAc in the regulation of dopamine receptor expression and drug reward.

Genome scale evolution of myxoma virus reveals host-pathogen adaptation and rapid geographic spread.

The evolutionary interplay between myxoma virus (MYXV) and the European rabbit (Oryctolagus cuniculus) following release of the virus in Australia in 1950 as a biological control is a classic example of host-pathogen coevolution. We present a detailed genomic and phylogeographic analysis of 30 strains of MYXV, including the Australian progenitor strain Standard Laboratory Strain (SLS), 24 Australian viruses isolated from 1951 to 1999, and three isolates from the early radiation in Britain from 1954 and 1955. We show that in Australia MYXV has spread rapidly on a spatial scale, with multiple lineages cocirculating within individual localities, and that both highly virulent and attenuated viruses were still present in the field through the 1990s. In addition, the detection of closely related virus lineages at sites 1,000 km apart suggests that MYXV moves freely in geographic space, with mosquitoes, fleas, and rabbit migration all providing means of transport. Strikingly, despite multiple introductions, all modern viruses appear to be ultimately derived from the original introductions of SLS. The rapidity of MYXV evolution was also apparent at the genomic scale, with gene duplications documented in a number of viruses. Duplication of potential virulence genes may be important in increasing the expression of virulence proteins and provides the basis for the evolution of novel functions. Mutations leading to loss of open reading frames were surprisingly frequent and in some cases may explain attenuation, but no common mutations that correlated with virulence or attenuation were identified.

Mammalian adaptation in the PB2 gene of avian H5N1 influenza virus.

The substitution of glutamic acid (E) for lysine (K) at position 627 of the PB2 protein of avian H5N1 viruses has been identified as a virulence and host range determinant for infection of mammals. Here, we report that the E-to-K host-adaptive mutation in the PB2 gene appeared from day 4 and 5 along the respiratory tracts of mice and was complete by day 6 postinoculation. This mutation correlated with efficient replication of the virus in mice.

Sequence analysis of in vivo defective interfering-like RNA of influenza A H1N1 pandemic virus.

Influenza virus defective interfering (DI) particles are naturally occurring noninfectious virions typically generated during in vitro serial passages in cell culture of the virus at a high multiplicity of infection. DI particles are recognized for the role they play in inhibiting viral replication and for the impact they have on the production of infectious virions. To date, influenza virus DI particles have been reported primarily as a phenomenon of cell culture and in experimentally infected embryonated chicken eggs. They have also been isolated from a respiratory infection of chickens. Using a sequencing approach, we characterize several subgenomic viral RNAs from human nasopharyngeal specimens infected with the influenza A(H1N1)pdm09 virus. The distribution of these in vivo-derived DI-like RNAs was similar to that of in vitro DIs, with the majority of the defective RNAs generated from the PB2 (segment 1) of the polymerase complex, followed by PB1 and PA. The lengths of the in vivo-derived DI-like segments also are similar to those of known in vitro DIs, and the in vivo-derived DI-like segments share internal deletions of the same segments. The presence of identical DI-like RNAs in patients linked by direct contact is compatible with transmission between them. The functional role of DI-like RNAs in natural infections remains to be established.