Emergence of two distinct spatial folds in a pair of plant virus proteins encoded by nested genes.

Virus genomes may encode overlapping or nested open reading frames that increase their coding capacity. It is not known whether the constraints on spatial structures of the two encoded proteins limit the evolvability of nested genes. We examine the evolution of a pair of proteins, p22 and p19, encoded by nested genes in plant viruses from the genus Tombusvirus. The known structure of p19, a suppressor of RNA silencing, belongs to the RAGNYA fold from the alpha+beta class. The structure of p22, the cell-to-cell movement protein from the 30K family widespread in plant viruses, is predicted with the AlphaFold approach, suggesting a single jelly-roll fold core from the all-beta class, structurally similar to capsid proteins from plant and animal viruses. The nucleotide and codon preferences impose modest constraints on the types of secondary structures encoded in the alternative reading frames, nonetheless allowing for compact, well-ordered folds from different structural classes in two similarly-sized nested proteins. Tombusvirus p22 emerged through radiation of the widespread 30K family, which evolved by duplication of a virus capsid protein early in the evolution of plant viruses, whereas lineage-specific p19 may have emerged by a stepwise increase in the length of the overprinted gene and incremental acquisition of functionally active secondary structure elements by the protein product. This evolution of p19 toward the RAGNYA fold represents one of the first documented examples of protein structure convergence in naturally occurring proteins.

Amino acid sequence homology of monoclonal serum free light chain dimers and tissue deposited light chains in AL amyloidosis: a pilot study.

OBJECTIVES: Diagnosis of light chain amyloidosis (AL) requires demonstration of amyloid deposits in a tissue biopsy followed by appropriate typing. Previous studies demonstrated increased dimerization of monoclonal serum free light chains (FLCs) as a pathological feature of AL. To further examine the pathogenicity of FLC, we aimed at testing amino acid sequence homology between circulating and deposited light chains (LCs). METHODS: Matched tissue biopsy and serum of 10 AL patients were subjected to tissue proteomic amyloid typing and nephelometric FLC assay, respectively. Serum FLC monomers (M) and dimers (D) were analyzed by Western blotting (WB) and mass spectrometry (MS). RESULTS: WB of serum FLCs showed predominance of either κ or λ type, in agreement with the nephelometric assay data. Abnormal FLC M-D patterns typical of AL amyloidosis were demonstrated in 8 AL-λ patients and in one of two AL-κ patients: increased levels of monoclonal FLC dimers, high D/M ratio values of involved FLCs, and high ratios of involved to uninvolved dimeric FLCs. MS of serum FLC dimers showed predominant constant domain sequences, in concordance with the tissue proteomic amyloid typing. Most importantly, variable domain sequence homology between circulating and deposited LC species was demonstrated, mainly in AL-λ cases. CONCLUSIONS: This is the first study to demonstrate homology between circulating FLCs and tissue-deposited LCs in AL-λ amyloidosis. The applied methodology can facilitate studying the pathogenicity of circulating FLC dimers in AL amyloidosis. The study also highlights the potential of FLC monomer and dimer analysis as a non-invasive screening tool for this disease.

Crystal structure of adenylosuccinate lyase from the thermophilic bacterium Thermus thermophilus HB8.

Adenylosuccinate lyase (PurB) catalyzes two distinct reactions in the purine nucleotide biosynthetic pathway using the same active site. The ability to recognize two different sets of substrates is of structural and evolutionary interest. In the present study, the crystal structure of PurB from the thermophilic bacterium Thermus thermophilus HB8 (TtPurB) was determined at a resolution of 2.38 Šby molecular replacement using a structure predicted by AlphaFold2 as a template. The asymmetric unit of the TtPurB crystal contained two TtPurB molecules, and some regions were disordered in the crystal structure. The disordered regions were the substrate-binding site and domain 3. TtPurB forms a homotetramer and the monomer is composed of three domains (domains 1, 2 and 3), which is a typical structure for the aspartase/fumarase superfamily. Molecular dynamics simulations with and without substrate/product were performed using a full-length model of TtPurB which was obtained before deletion of the disordered regions. The substrates and products were bound to the model structures during the MD simulations. The fluctuations of amino-acid residues were greater in the disordered regions and became smaller upon the binding of substrate or product. These results demonstrate that the full-length model obtained using AlphaFold2 can be used to generate the coordinates of disordered regions within the crystal structure.

Identification of a covert evolutionary pathway between two protein folds.

Although homologous protein sequences are expected to adopt similar structures, some amino acid substitutions can interconvert α-helices and ß-sheets. Such fold switching may have occurred over evolutionary history, but supporting evidence has been limited by the: (1) abundance and diversity of sequenced genes, (2) quantity of experimentally determined protein structures, and (3) assumptions underlying the statistical methods used to infer homology. Here, we overcome these barriers by applying multiple statistical methods to a family of ~600,000 bacterial response regulator proteins. We find that their homologous DNA-binding subunits assume divergent structures: helix-turn-helix versus α-helix + ß-sheet (winged helix). Phylogenetic analyses, ancestral sequence reconstruction, and AlphaFold2 models indicate that amino acid substitutions facilitated a switch from helix-turn-helix into winged helix. This structural transformation likely expanded DNA-binding specificity. Our approach uncovers an evolutionary pathway between two protein folds and provides a methodology to identify secondary structure switching in other protein families.

KinFams: De-Novo Classification of Protein Kinases Using CATH Functional Units.

Protein kinases are important targets for treating human disorders, and they are the second most targeted families after G-protein coupled receptors. Several resources provide classification of kinases into evolutionary families (based on sequence homology); however, very few systematically classify functional families (FunFams) comprising evolutionary relatives that share similar functional properties. We have developed the FunFam-MARC (Multidomain ARchitecture-based Clustering) protocol, which uses multi-domain architectures of protein kinases and specificity-determining residues for functional family classification. FunFam-MARC predicts 2210 kinase functional families (KinFams), which have increased functional coherence, in terms of EC annotations, compared to the widely used KinBase classification. Our protocol provides a comprehensive classification for kinase sequences from >10,000 organisms. We associate human KinFams with diseases and drugs and identify 28 druggable human KinFams, i.e., enriched in clinically approved drugs. Since relatives in the same druggable KinFam tend to be structurally conserved, including the drug-binding site, these KinFams may be valuable for shortlisting therapeutic targets. Information on the human KinFams and associated 3D structures from AlphaFold2 are provided via our CATH FTP website and Zenodo. This gives the domain structure representative of each KinFam together with information on any drug compounds available. For 32% of the KinFams, we provide information on highly conserved residue sites that may be associated with specificity.

CATHe: detection of remote homologues for CATH superfamilies using embeddings from protein language models.

MOTIVATION: CATH is a protein domain classification resource that exploits an automated workflow of structure and sequence comparison alongside expert manual curation to construct a hierarchical classification of evolutionary and structural relationships. The aim of this study was to develop algorithms for detecting remote homologues missed by state-of-the-art hidden Markov model (HMM)-based approaches. The method developed (CATHe) combines a neural network with sequence representations obtained from protein language models. It was assessed using a dataset of remote homologues having less than 20% sequence identity to any domain in the training set. RESULTS: The CATHe models trained on 1773 largest and 50 largest CATH superfamilies had an accuracy of 85.6 ± 0.4% and 98.2 ± 0.3%, respectively. As a further test of the power of CATHe to detect more remote homologues missed by HMMs derived from CATH domains, we used a dataset consisting of protein domains that had annotations in Pfam, but not in CATH. By using highly reliable CATHe predictions (expected error rate <0.5%), we were able to provide CATH annotations for 4.62 million Pfam domains. For a subset of these domains from Homo sapiens, we structurally validated 90.86% of the predictions by comparing their corresponding AlphaFold2 structures with structures from the CATH superfamilies to which they were assigned. AVAILABILITY AND IMPLEMENTATION: The code for the developed models is available on https://github.com/vam-sin/CATHe, and the datasets developed in this study can be accessed on https://zenodo.org/record/6327572. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Transcriptome Sequencing and Comparison of Venom Glands Revealed Intraspecific Differentiation and Expression Characteristics of Toxin and Defensin Genes in Mesobuthus martensii Populations.

Mesobuthus martensii, a famous and important Traditional Chinese Medicine has a long medical history and unique functions. It is the first scorpion species whose whole genome was sequenced worldwide. In addition, it is the most widespread and infamous poisonous animal in northern China with complex habitats. It possesses several kinds of toxins that can regulate different ion channels and serve as crucial natural drug resources. Extensive and in-depth studies have been performed on the structures and functions of toxins of M. martensii. In this research, we compared the morphology of M. martensii populations from different localities and calculated the COI genetic distance to determine intraspecific variations. Transcriptome sequencing by RNA-sequencing of the venom glands of M. martensii from ten localities and M. eupeus from one locality was analyzed. The results revealed intraspecific variation in the expression of sodium channel toxin genes, potassium channel toxin genes, calcium channel toxin genes, chloride channel toxin genes, and defensin genes that could be related to the habitats in which these populations are distributed, except the genetic relationships. However, it is not the same in different toxin families. M. martensii and M. eupeus exhibit sexual dimorphism under the expression of toxin genes, which also vary in different toxin families. The following order was recorded in the difference of expression of sodium channel toxin genes: interspecific difference; differences among different populations of the same species; differences between sexes in the same population, whereas the order in the difference of expression of potassium channel toxin genes was interspecific difference; differences between both sexes of same populations; differences among the same sex in different populations of the same species. In addition, there existed fewer expressed genes of calcium channel toxins, chloride channel toxins, and defensins (no more than four members in each family), and their expression differences were not distinct. Interestingly, the expression of two calcium channel toxin genes showed a preference for males and certain populations. We found a difference in the expression of sodium channel toxin genes, potassium channel toxin genes, and chloride channel toxin genes between M. martensii and M. eupeus. In most cases, the expression of one member of the toxin gene clusters distributed in series on the genome were close in different populations and genders, and the members of most clusters expressed in same population and gender tended to be the different. Twenty-one toxin genes were found with the MS/MS identification evidence of M. martensii venom. Since scorpions were not subjected to electrical stimulation or other special treatments before conducting the transcriptome extraction experiment, the results suggested the presence of intraspecific variation and sexual dimorphism of toxin components which revealed the expression characteristics of toxin and defensin genes in M. martensii. We believe this study will promote further in-depth research and use of scorpions and their toxin resources, which in turn will be helpful in standardizing the identification and medical applications of Quanxie in traditional Chinese medicine.

DeepRHD: An efficient hybrid feature extraction technique for protein remote homology detection using deep learning strategies.

In computational biology, the Protein Remote homology Detection technique (PRHD) has got undeniable significance. It is mostly important for structure and function identification of a protein sequence. The previous years have seen a challenge that lacks postulating a correlation among the sequences. However, the sequences are of variable length. Thereby, it inhibits the proper derivation of evolutionary information among the sequences. The challenges are the usage of physico-chemical properties as a source to get the evolutionary information and the number of sequences generated every day. This however facilitates a new technique to integrate huge amount of data with a massive feature set. In this article, a new and efficient technique is proposed to predict homology for distantly located sequences of proteins. Deep neural network(CNN-GRU model) is used for the classification of the protein sequences. This is based on different protein families and methods of feature extraction.The efficiency of the proposed model DeepRHD is tested on average 8000 sequences per superfamily taken from SCOP benchmark dataset and the results shows that the proposed model is better than other state of art methods. This model is useful in detecting diseases like sickle cell anemia and influenza and developing a drug thereafter.

Advancing remote homology detection: A step toward understanding and accurately predicting protein function.

Identifying homologous proteins with divergent amino acid sequences can add to our understanding of protein evolution, structure, and function. A new study reports the development of a deep-network-based method to identify 6.8 million new Pfam members, a dramatic singular increase that exceeds a decade of accumulation using traditional approaches.

Genome wide identification of MADS box gene family in Musa balbisiana and their divergence during evolution.

MADS box gene family is transcription factor gene family that is involved in growth and development of eukaryotes. In plants the MADS box gene family is mainly associated with floral meristem identity and flower development, apart from being involved in nearly all the phases of plant growth. The MADS box gene family has also been shown to be involved during fruit development and ripening. In this study the MADS box gene family from Musa balbisiana was identified and the divergence of this gene family between Musa balbisiana and Musa acuminata studied. A total of 97 MADS box genes were identified from the genome of Musa balbisiana. Phylogenetic analysis showed that the MbMADS box genes were categorised into type I (α and γ; the ß group was not distinguishable) and type II groups (MIKCc and MIKC* and MIKCc was further divided into 13 subfamilies). The typeII group has the largest number of genes and also showed the most expansion which could be correlated with the whole genome duplications. There were significant differences in the MADS box genes from Musa acuminata and Musa balbisiana during evolution that can be correlated with different floral phenotype and fruit ripening pattern. The divergence of the MADS RIN genes in Musa balbisiana as compared to Musa acuminata might play an important role in the slow ripening of Musa balbisiana fruits.

On the reliability and the limits of inference of amino acid sequence alignments.

MOTIVATION: Alignments are correspondences between sequences. How reliable are alignments of amino acid sequences of proteins, and what inferences about protein relationships can be drawn? Using techniques not previously applied to these questions, by weighting every possible sequence alignment by its posterior probability we derive a formal mathematical expectation, and develop an efficient algorithm for computation of the distance between alternative alignments allowing quantitative comparisons of sequence-based alignments with corresponding reference structure alignments. RESULTS: By analyzing the sequences and structures of 1 million protein domain pairs, we report the variation of the expected distance between sequence-based and structure-based alignments, as a function of (Markov time of) sequence divergence. Our results clearly demarcate the 'daylight', 'twilight' and 'midnight' zones for interpreting residue-residue correspondences from sequence information alone. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cross-Recognition of SARS-CoV-2 B-Cell Epitopes with Other Betacoronavirus Nucleoproteins.

The B and T lymphocytes of the adaptive immune system are important for the control of most viral infections, including COVID-19. Identification of epitopes recognized by these cells is fundamental for understanding how the immune system detects and removes pathogens, and for antiviral vaccine design. Intriguingly, several cross-reactive T lymphocyte epitopes from SARS-CoV-2 with other betacoronaviruses responsible for the common cold have been identified. In addition, antibodies that cross-recognize the spike protein, but not the nucleoprotein (N protein), from different betacoronavirus have also been reported. Using a consensus of eight bioinformatic methods for predicting B-cell epitopes and the collection of experimentally detected epitopes for SARS-CoV and SARS-CoV-2, we identified four surface-exposed, conserved, and hypothetical antigenic regions that are exclusive of the N protein. These regions were analyzed using ELISA assays with two cohorts: SARS-CoV-2 infected patients and pre-COVID-19 samples. Here we describe four epitopes from SARS-CoV-2 N protein that are recognized by the humoral response from multiple individuals infected with COVID-19, and are conserved in other human coronaviruses. Three of these linear surface-exposed sequences and their peptide homologs in SARS-CoV-2 and HCoV-OC43 were also recognized by antibodies from pre-COVID-19 serum samples, indicating cross-reactivity of antibodies against coronavirus N proteins. Different conserved human coronaviruses (HCoVs) cross-reactive B epitopes against SARS-CoV-2 N protein are detected in a significant fraction of individuals not exposed to this pandemic virus. These results have potential clinical implications.

Crystal structures of TTHA1265 and TTHA1264/TTHA1265 complex reveal an intrinsic heterodimeric assembly.

Zinc peptidase M16 family members are widely distributed in most prokaryotic and eukaryotic organisms. M16 family has been divided into three subfamilies, M16A, M16B and M16C, based on sequence alignments and subunit connectivity. TTHA1264, an M16B protein found in Thermus thermophiles HB8, possesses an HXXEH motif essential for Zn2+ binding and catalytic activity. TTHA1265 is another member of M16B, which lacks the metal-binding motif but with a conserved active-site R/Y pair commonly found in the C-terminal half of M16 enzymes. Sequence analysis showed that two genes coding for TTHA1264 and TTHA1265 assemble into a single operon in the bacterial genome. Here, we report the crystal structure of TTHA1265 and TTHA1264/TTHA1265 complex from T. thermophilus HB8. Interestingly, when TTHA1264 and TTHA1265 are present alone, TTHA1264 forms a monomer, TTHA1265 forms a homodimer, respectively. However, TTHA264 and TTHA1265 assembled into a heterodimeric complex, indicating that they prefer to form heterodimer. Biochemical data further confirmed the heterodimeric assembly indicating intrinsic heterodimeric assembly of TTHA1264 and TTHA1265. This property of TTHA1264 and TTHA1265 is consistent with the characteristics of the M16B family.

The microprotein Nrs1 rewires the G1/S transcriptional machinery during nitrogen limitation in budding yeast.

Commitment to cell division at the end of G1 phase, termed Start in the budding yeast Saccharomyces cerevisiae, is strongly influenced by nutrient availability. To identify new dominant activators of Start that might operate under different nutrient conditions, we screened a genome-wide ORF overexpression library for genes that bypass a Start arrest caused by absence of the G1 cyclin Cln3 and the transcriptional activator Bck2. We recovered a hypothetical gene YLR053c, renamed NRS1 for Nitrogen-Responsive Start regulator 1, which encodes a poorly characterized 108 amino acid microprotein. Endogenous Nrs1 was nuclear-localized, restricted to poor nitrogen conditions, induced upon TORC1 inhibition, and cell cycle-regulated with a peak at Start. NRS1 interacted genetically with SWI4 and SWI6, which encode subunits of the main G1/S transcription factor complex SBF. Correspondingly, Nrs1 physically interacted with Swi4 and Swi6 and was localized to G1/S promoter DNA. Nrs1 exhibited inherent transactivation activity, and fusion of Nrs1 to the SBF inhibitor Whi5 was sufficient to suppress other Start defects. Nrs1 appears to be a recently evolved microprotein that rewires the G1/S transcriptional machinery under poor nitrogen conditions.

The ER membrane complex (EMC) can functionally replace the Oxa1 insertase in mitochondria.

Two multisubunit protein complexes for membrane protein insertion were recently identified in the endoplasmic reticulum (ER): the guided entry of tail anchor proteins (GET) complex and ER membrane complex (EMC). The structures of both of their hydrophobic core subunits, which are required for the insertion reaction, revealed an overall similarity to the YidC/Oxa1/Alb3 family members found in bacteria, mitochondria, and chloroplasts. This suggests that these membrane insertion machineries all share a common ancestry. To test whether these ER proteins can functionally replace Oxa1 in yeast mitochondria, we generated strains that express mitochondria-targeted Get2-Get1 and Emc6-Emc3 fusion proteins in Oxa1 deletion mutants. Interestingly, the Emc6-Emc3 fusion was able to complement an Δoxa1 mutant and restored its respiratory competence. The Emc6-Emc3 fusion promoted the insertion of the mitochondrially encoded protein Cox2, as well as of nuclear encoded inner membrane proteins, although was not able to facilitate the assembly of the Atp9 ring. Our observations indicate that protein insertion into the ER is functionally conserved to the insertion mechanism in bacteria and mitochondria and adheres to similar topological principles.

SIRT2 regulates proliferation and chemotherapy response of MLL-ENL-driven acute myeloid leukemia.

Both MLL-AF9 and MLL-ENL leukemia fusion proteins drive oncogenic transformation of hematopoietic cells through their N-terminal DNA/histone binding mixed-lineage leukemia 1 domain and C-terminal fragment of AF9 or ENL containing an unstructured linker region and the ANC1 homology domain, which recruits transcription factors. Despite of their structural similarity, acute myeloid leukemia (AML) patients bearing MLL-ENL show more adverse outcomes compared to those with MLL-AF9. We recapitulated the clinical patterns of these two MLL-fusions driven AMLs using murine models and found that MLL-ENL AML cells showed slower cell cycle progression and more resistance to standard chemotherapy than MLL-AF9 cells. These phenotypes were primarily controlled by the linker regions of ENL and a highly conserved lysine residue K469 within. Substitution of K469 with an acetylated mimic glutamine abolished the ability of MLL-ENL to suppress proliferation and promote chemo-resistance. We showed that deacetylase Sirt2 might act as an upstream regulator of MLL-ENL. Deletion of Sirt2 promoted proliferation of AML cells with either MLL fusions. Importantly, loss of Sirt2 greatly enhanced the sensitivity of the MLL-ENL AML cells to chemo-treatment. Taken together, our study uncovered a unique regulatory role of Sirt2 in leukemogenesis and suggested targeting SIRT2 as a new way to sensitize MLL-ENL AML patience for chemotherapy.

Specialization of actin isoforms derived from the loss of key interactions with regulatory factors.

A paradox of eukaryotic cells is that while some species assemble a complex actin cytoskeleton from a single ortholog, other species utilize a greater diversity of actin isoforms. The physiological consequences of using different actin isoforms, and the molecular mechanisms by which highly conserved actin isoforms are segregated into distinct networks, are poorly known. Here, we sought to understand how a simple biological system, composed of a unique actin and a limited set of actin-binding proteins, reacts to a switch to heterologous actin expression. Using yeast as a model system and biomimetic assays, we show that such perturbation causes drastic reorganization of the actin cytoskeleton. Our results indicate that defective interaction of a heterologous actin for important regulators of actin assembly limits certain actin assembly pathways while reinforcing others. Expression of two heterologous actin variants, each specialized in assembling a different network, rescues cytoskeletal organization and confers resistance to external perturbation. Hence, while species using a unique actin have homeostatic actin networks, actin assembly pathways in species using several actin isoforms may act more independently.

Drosophila carrying epilepsy-associated variants in the vitamin B6 metabolism gene PNPO display allele- and diet-dependent phenotypes.

Pyridox(am)ine 5 ' -phosphate oxidase (PNPO) catalyzes the rate-limiting step in the synthesis of pyridoxal 5 ' -phosphate (PLP), the active form of vitamin B6 required for the synthesis of neurotransmitters gamma-aminobutyric acid (GABA) and the monoamines. Pathogenic variants in PNPO have been increasingly identified in patients with neonatal epileptic encephalopathy and early-onset epilepsy. These patients often exhibit different types of seizures and variable comorbidities. Recently, the PNPO gene has also been implicated in epilepsy in adults. It is unclear how these phenotypic variations are linked to specific PNPO alleles and to what degree diet can modify their expression. Using CRISPR-Cas9, we generated four knock-in Drosophila alleles, hWT , hR116Q , hD33V , and hR95H , in which the endogenous Drosophila PNPO was replaced by wild-type human PNPO complementary DNA (cDNA) and three epilepsy-associated variants. We found that these knock-in flies exhibited a wide range of phenotypes, including developmental impairments, abnormal locomotor activities, spontaneous seizures, and shortened life span. These phenotypes are allele dependent, varying with the known biochemical severity of these mutations and our characterized molecular defects. We also showed that diet treatments further diversified the phenotypes among alleles, and PLP supplementation at larval and adult stages prevented developmental impairments and seizures in adult flies, respectively. Furthermore, we found that hR95H had a significant dominant-negative effect, rendering heterozygous flies susceptible to seizures and premature death. Together, these results provide biological bases for the various phenotypes resulting from multifunction of PNPO, specific molecular and/or genetic properties of each PNPO variant, and differential allele-diet interactions.

Synthetic virions reveal fatty acid-coupled adaptive immunogenicity of SARS-CoV-2 spike glycoprotein.

SARS-CoV-2 infection is a major global public health concern with incompletely understood pathogenesis. The SARS-CoV-2 spike (S) glycoprotein comprises a highly conserved free fatty acid binding pocket (FABP) with unknown function and evolutionary selection advantage1,2. Deciphering FABP impact on COVID-19 progression is challenged by the heterogenous nature and large molecular variability of live virus. Here we create synthetic minimal virions (MiniVs) of wild-type and mutant SARS-CoV-2 with precise molecular composition and programmable complexity by bottom-up assembly. MiniV-based systematic assessment of S free fatty acid (FFA) binding reveals that FABP functions as an allosteric regulatory site enabling adaptation of SARS-CoV-2 immunogenicity to inflammation states via binding of pro-inflammatory FFAs. This is achieved by regulation of the S open-to-close equilibrium and the exposure of both, the receptor binding domain (RBD) and the SARS-CoV-2 RGD motif that is responsible for integrin co-receptor engagement. We find that the FDA-approved drugs vitamin K and dexamethasone modulate S-based cell binding in an FABP-like manner. In inflammatory FFA environments, neutralizing immunoglobulins from human convalescent COVID-19 donors lose neutralization activity. Empowered by our MiniV technology, we suggest a conserved mechanism by which SARS-CoV-2 dynamically couples its immunogenicity to the host immune response.

Mass Screening of SARS-CoV-2 Variants using Sanger Sequencing Strategy in Hiroshima, Japan.

This study aimed to develop the feasible and effective universal screening strategy of the notable SARS-CoV-2 variants by Sanger Sequencing Strategy and then practically applied it for mass screening in Hiroshima, Japan. A total of 734 samples from COVID-19 confirmed cases in Hiroshima were screened for the notable SARS-CoV-2 variants (B.1.1.7, B.1.351, P.1, B.1.617.2, B.1.617.1, C.37, B.1.1.529, etc.). The targeted spike region is amplified by nested RT-PCR using in-house designed primer set hCoV-Spike-A and standard amplification protocol. Additionally, randomly selected 96 samples were also amplified using primer sets hCoV-Spike-B and hCoV-Spike-C. The negative amplified samples were repeated for second attempt of amplification by volume-up protocol. Thereafter, the amplified products were assigned for Sanger sequencing using corresponding primers. The positive amplification rate of primer set hCoV-Spike-A, hCoV-Spike-B and hCoV-Spike-C were 87.3%, 83.3% and 93.8% respectively for standard protocol and increased to 99.6%, 95.8% and 96.9% after second attempt by volume-up protocol. The readiness of genome sequences was 96.9%, 100% and 100% respectively. Among 48 mutant isolates, 26 were B.1.1.7 (Alpha), 7 were E484K single mutation and the rest were other types of mutation. Moreover, 5 cluster cases with single mutation at N501S were firstly reported in Hiroshima. This study indicates the reliability and effectiveness of Sanger sequencing to screen large number of samples for the notable SARS-CoV-2 variants. Compared to the Next Generation Sequencing (NGS), our method introduces the feasible, universally applicable, and practically useful tool for identification of the emerging variants with less expensive and time consuming especially in those countries where the NGS is not practically available. Our method allows not only to identify the pre-existing variants but also to examine other rare type of mutation or newly emerged variants and is crucial for prevention and control of pandemic.