During the current coronavirus disease 2019 (COVID-19) pandemic, a variety of mutations have accumulated in the viral genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and, at the time of writing, four variants of concern are considered to be potentially hazardous to human society¹. The recently emerged B.1.617.2/Delta variant of concern is closely associated with the COVID-19 surge that occurred in India in the spring of 2021 (ref. ²). However, the virological properties of B.1.617.2/Delta remain unclear. Here we show that the B.1.617.2/Delta variant is highly fusogenic and notably more pathogenic than prototypic SARS-CoV-2 in infected hamsters. The P681R mutation in the spike protein, which is highly conserved in this lineage, facilitates cleavage of the spike protein and enhances viral fusogenicity. Moreover, we demonstrate that the P681R-bearing virus exhibits higher pathogenicity compared with its parental virus. Our data suggest that the P681R mutation is a hallmark of the virological phenotype of the B.1.617.2/Delta variant and is associated with enhanced pathogenicity.

Many SARS-CoV-2 variants with naturally acquired mutations have emerged. These mutations can affect viral properties such as infectivity and immune resistance. Although the sensitivity of naturally occurring SARS-CoV-2 variants to humoral immunity has been investigated, sensitivity to human leukocyte antigen (HLA)-restricted cellular immunity remains largely unexplored. Here, we demonstrate that two recently emerging mutations in the receptor-binding domain of the SARS-CoV-2 spike protein, L452R (in B.1.427/429 and B.1.617) and Y453F (in B.1.1.298), confer escape from HLA-A24-restricted cellular immunity. These mutations reinforce affinity toward the host entry receptor ACE2. Notably, the L452R mutation increases spike stability, viral infectivity, viral fusogenicity, and thereby promotes viral replication. These data suggest that HLA-restricted cellular immunity potentially affects the evolution of viral phenotypes and that a further threat of the SARS-CoV-2 pandemic is escape from cellular immunity.

The evolution of transcriptional regulatory networks entails the expansion and diversification of transcription factor (TF) families. The forkhead family of TFs, defined by a highly conserved winged helix DNA-binding domain (DBD), has diverged into dozens of subfamilies in animals, fungi, and related protists. We have used a combination of maximum-likelihood phylogenetic inference and independent, comprehensive functional assays of DNA-binding capacity to explore the evolution of DNA-binding specificity within the forkhead family. We present converging evidence that similar alternative sequence preferences have arisen repeatedly and independently in the course of forkhead evolution. The vast majority of DNA-binding specificity changes we observed are not explained by alterations in the known DNA-contacting amino acid residues conferring specificity for canonical forkhead binding sites. Intriguingly, we have found forkhead DBDs that retain the ability to bind very specifically to two completely distinct DNA sequence motifs. We propose an alternate specificity-determining mechanism whereby conformational rearrangements of the DBD broaden the spectrum of sequence motifs that a TF can recognize. DNA-binding bispecificity suggests a previously undescribed source of modularity and flexibility in gene regulation and may play an important role in the evolution of transcriptional regulatory networks.

Understanding regulatory mechanisms of protein synthesis in eukaryotes is essential for the accurate annotation of genome sequences. Kozak reported that the nucleotide sequence GCCGCC(A/G)CCAUGG (AUG is the initiation codon) was frequently observed in vertebrate genes and that this ‘consensus’ sequence enhanced translation initiation. However, later studies using invertebrate, fungal and plant genes reported different ‘consensus’ sequences. In this study, we conducted extensive comparative analyses of nucleotide sequences around the initiation codon by using genomic data from 47 eukaryote species including animals, fungi, plants and protists. The analyses revealed that preferred nucleotide sequences are quite diverse among different species, but differences between patterns of nucleotide bias roughly reflect the evolutionary relationships of the species. We also found strong biases of A/G at position −3, A/C at position −2 and C at position +5 that were commonly observed in all species examined. Genes with higher expression levels showed stronger signals, suggesting that these nucleotides are responsible for the regulation of translation initiation. The diversity of preferred nucleotide sequences around the initiation codon might be explained by differences in relative contributions from two distinct patterns, GCCGCCAUG and AAAAAAAUG, which implies the presence of multiple molecular mechanisms for controlling translation initiation.

SARS-CoV-2 Lambda, a variant of interest, has spread in some South American countries; however, its virological features and evolutionary traits remain unclear. In this study, we use pseudoviruses and reveal that the spike protein of the Lambda variant is more infectious than that of other variants due to the T76I and L452Q mutations. The RSYLTPGD246-253N mutation, a unique 7-amino acid deletion in the N-terminal domain of the Lambda spike protein, is responsible for evasion from neutralizing antibodies and further augments antibody-mediated enhancement of infection. Although this mutation generates a nascent N-linked glycosylation site, the additional N-linked glycan is dispensable for the virological property conferred by this mutation. Since the Lambda variant has dominantly spread according to the increasing frequency of the isolates harboring the RSYLTPGD246-253N mutation, our data suggest that the RSYLTPGD246-253N mutation is closely associated with the substantial spread of the Lambda variant in South America.

Rapid identification of bacterial pathogens is crucial for appropriate and adequate antibiotic treatment, which significantly improves patient outcomes. 16S ribosomal RNA (rRNA) gene amplicon sequencing has proven to be a powerful strategy for diagnosing bacterial infections. We have recently established a sequencing method and bioinformatics pipeline for 16S rRNA gene analysis utilizing the Oxford Nanopore Technologies MinION™ sequencer. In combination with our taxonomy annotation analysis pipeline, the system enabled the molecular detection of bacterial DNA in a reasonable time frame for diagnostic purposes. However, purification of bacterial DNA from specimens remains a rate‐limiting step in the workflow. To further accelerate the process of sample preparation, we adopted a direct PCR strategy that amplifies 16S rRNA genes from bacterial cell suspensions without DNA purification. Our results indicate that differences in cell wall morphology significantly affect direct PCR efficiency and sequencing data. Notably, mechanical cell disruption preceding direct PCR was indispensable for obtaining an accurate representation of the specimen bacterial composition. Furthermore, 16S rRNA gene analysis of mock polymicrobial samples indicated that primer sequence optimization is required to avoid preferential detection of particular taxa and to cover a broad range of bacterial species. This study establishes a relatively simple workflow for rapid bacterial identification via MinION™ sequencing, which reduces the turnaround time from sample to result, and provides a reliable method that may be applicable to clinical settings.