Next Article in Journal
Unusual Localization of AIDS-Related Kaposi’s Sarcoma in a Heterosexual Male during the COVID-19 Pandemic: A Case Report
Next Article in Special Issue
A Longitudinal Study in Tunisia to Assess the Anti-RBD IgG and IgA Responses Induced by Three Different COVID-19 Vaccine Platforms
Previous Article in Journal
Epidemiology of Dengue in Argentina during the 2010/11 to 2019/20 Seasons: A Contribution to the Burden of Disease
Previous Article in Special Issue
Mortality Trend of Severe COVID-19 in Under-Vaccinated Population Admitted to ICU in French Amazonia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic and Phylogenetic Characterisation of SARS-CoV-2 Genomes Isolated in Patients from Lambayeque Region, Peru

by
Sergio Luis Aguilar-Martinez
1,
Gustavo Adolfo Sandoval-Peña
1,
José Arturo Molina-Mora
2,
Pablo Tsukayama-Cisneros
3,
Cristian Díaz-Vélez
4,
Franklin Rómulo Aguilar-Gamboa
5,
D. Katterine Bonilla-Aldana
6,* and
Alfonso J. Rodriguez-Morales
7,8
1
Grupo de Investigación en Bioinformática y Biología Estructural, Facultad de Ciencias Biológicas, Universidad Nacional Mayor de San Marcos, Lima 15081, Peru
2
Centro de Investigación en Enfermedades Tropicales, Facultad de Microbiología, Universidad de Costa Rica, San Jose 11501-2060, Costa Rica
3
Facultad de Ciencias e Ingeniería, Universidad Peruana Cayetano Heredia, Lima 15102, Peru
4
Facultad de Medicina, Universidad Privada Antenor Orrego, Trujillo 13008, Peru
5
Grupo de Investigación en Inmunología y Virología del Norte, Laboratorio de Inmunología y Virología, Hospital Regional Lambayeque, Lambayeque 14011, Peru
6
Research Unit, Universidad Continental, Huancayo 12001, Peru
7
Faculty of Health Sciences, Universidad Científica del Sur, Lima 15307, Peru
8
Gilbert and Rose-Marie Chagoury School of Medicine, Lebanese American University, Beirut P.O. Box 13-5053, Lebanon
*
Author to whom correspondence should be addressed.
Trop. Med. Infect. Dis. 2024, 9(2), 46; https://doi.org/10.3390/tropicalmed9020046
Submission received: 22 November 2023 / Revised: 1 February 2024 / Accepted: 2 February 2024 / Published: 11 February 2024
(This article belongs to the Special Issue COVID-19 Variants, Vaccines and New Waves)

Abstract

:
Objective: this study aims to identify and characterise genomic and phylogenetically isolated SARS-CoV-2 viral isolates in patients from Lambayeque, Peru. Methods: Nasopharyngeal swabs were taken from patients from the Almanzor Aguinaga Asenjo Hospital, Chiclayo, Lambayeque, Peru, which had been considered mild, moderate, and severe cases of COVID-19. Patients had to have tested positive for COVID-19, using a positive RT-PCR for SARS-CoV-2. Subsequently, the SARS-CoV-2 complete viral genome sequencing was carried out using Illumina MiSeq®. The sequences obtained from the sequence were analysed in Nextclade V1.10.0 to assign the corresponding clades, identify mutations in the SARS-CoV-2 genes and perform quality control of the sequences obtained. All sequences were aligned using MAFFT v7.471. The SARS-CoV-2 isolate Wuhan NC 045512.2 was used as a reference sequence to analyse mutations at the amino acid level. The construction of the phylogenetic tree model was achieved with IQ-TREE v1.6.12. Results: It was determined that during the period from December 2020 to January 2021, the lineages s C.14, C.33, B.1.1.485, B.1.1, B.1.1.1, and B.1.111 circulated, with lineage C.14 being the most predominant at 76.7% (n = 23/30). These lineages were classified in clade 20D mainly and also within clades 20B and 20A. On the contrary, the variants found in the second batch of samples of the period from September to October 2021 were Delta (72.7%), Gamma (13.6%), Mu (4.6%), and Lambda (9.1%), distributed between clades 20J, 21G, 21H, 21J, and 21I. Conclusions: This study reveals updated information on the viral genomics of SARS-CoV-2 in the Lambayeque region, Peru, which is crucial to understanding the origins and dispersion of the virus and provides information on viral pathogenicity, transmission and epidemiology.

1. Introduction

COVID-19 is a respiratory disease caused by the SARS-CoV-2 virus (severe acute respiratory syndrome coronavirus 2), which was declared a pandemic by the World Health Organization (WHO) in early 2020. This disease has caused a health and economic emergency worldwide. Currently, research on SARS-CoV-2 is booming and great efforts are being made to characterise SARS-CoV-2 molecularly. The genomic and molecular variability of SARS-CoV-2 can be the basis for glimpses of etiological and pathological aspects of this virus, understanding that the virus can accumulate mutations of importance while expanding worldwide, as well as also be able to establish antiviral strategies designed and based on the molecular specificities of SARS-CoV-2.
One of the most striking aspects of COVID-19 is the marked difference in the evolution of the disease in patients. The spread and manifestations of COVID-19, an infectious disease, are influenced by multiple interrelated factors. These include the virus itself (SARS-CoV-2), the human host (comorbidities and genetics), and the environment (physical conditions, social interactions, containment measures). All of these play a role in determining the course of the disease and the pandemic [1]. By elucidating and obtaining these genomic data, it would be possible to reveal the evolutionary events of SARS-CoV-2, establish the types of circulating genomes, and determine in which parts of the genome these viral isolates differ [2,3].
In Peru, lineages of regional and global relevance variants have emerged; some researchers detected the circulation of SARS-CoV-2 strains with the D614G mutation in the Lambayeque region at the beginning of 2020. This mutation had already spread widely in Europe at that time. However, other uncommon mutations demonstrate the virus’s rapid evolutionary processes and adaptive capabilities [4]. Subsequent investigations corroborated the presence of a variant endemic to the region, which was designated the Lambda variant [5].
The genetic variability of SARS-CoV-2 requires continuous study to elucidate various aspects of its molecular biology. Due to this, various modifications or changes in the nucleotide sequence of the viral genome have been reported worldwide, causing the appearance of variants, which have been grouped into differentiated clades. Among the variants of interest of SARS-CoV-2 are Lambda and Mu, which were first identified in Peru and Colombia, respectively. Meanwhile, the variants of concern of SARS-CoV-2 identified and reported globally, in chronological order, are Alpha (British), Beta (South African), Gamma (Brazilian), Delta (Indian) and recently Omicron [6].
For this reason, sequencing of the SARS-CoV-2 viral genome in Peru is urgently required; this will provide information on the prevalence of viral clades belonging to SARS-CoV-2, which could lead to a better understanding of transmission patterns, outbreak monitoring and formulation of effective containment measures. Mutation data may also provide important clues for developing vaccines, antiviral drugs, and effective diagnostic assays.
This study aimed to investigate the genomic variation in SARS-CoV-2 in Peru through the whole genome sequencing of SARS-CoV-2 strains and compare their evolutionary trajectories with global strains through phylogenetic analysis.

2. Methods

2.1. Sample Collection and Complete Sequencing of the SARS-CoV-2 Genome

Nasopharyngeal swabs were obtained from positive cases of COVID-19, and the first 30 sequenced samples were obtained at the Almanzor Aguinaga Asenjo Hospital belonging to Essalud, Chiclayo, Lambayeque, Peru, during the period from December 2020 to January 2021. The samples with a CT ≤28 were selected for genome sequencing; these samples (n = 30) were then sent to the Microbial Genomics laboratory of the Universidad Peruana Cayetano Heredia for subsequent analysis and sequencing. Another 44 samples of nasopharyngeal swabs from positive cases of COVID-19 confirmed by the reference laboratory of the Regional Health Management Lambayeque were included; these biological samples belong to the period from September to October 2021, thus completing 74 genomes sequenced by the authors of the present research.
Whole genome sequencing of SARS-CoV-2 isolates was performed with a MiSeq (Illumina, San Diego, CA, USA) using the COVIDseq Illumina kit. In detail, the samples were processed using the Sansure brand RNA kit and a Sansure Natch 48 automated device. Subsequently, genomic libraries were prepared with the Illumina COVIDSeq kit. The library sequencing was conducted using an instrument of the Illumina MiSeq model and the 300-cycle v2 kit to achieve an average coverage of 1500X. Once the sequences were obtained, they were processed and assembled with Illumina’s DRAGEN pipeline through its Base Space application. The Illumina MiSeq sequencing procedure was performed at the Microbial Genomics Laboratory, Department of Cellular and Molecular Sciences Faculty of Sciences and Philosophy (Universidad Peruana Cayetano Heredia, Lima, Peru).

2.2. Bioinformatics Analysis

All bioinformatic analyses used in the development of this research were based on the protocols proposed and designed for the genomic surveillance of SARS-CoV-2 in Costa Rica [7].

2.3. SARS-CoV-2 Sequences

Until 28 April 2022, all SARS-CoV-2 sequences available from the Lambayeque region of Peru were recovered and retrieved from the GISAID (Global Initiative on Sharing All Influenza Data, www.gisaid.org) database. A total of 714 sequences were recovered from the Lambayeque region, of which 74 were sequenced and uploaded to GISAID by the authors of this article.

2.4. Multi-Sequence Alignment

All sequences were aligned using MAFFT v7.471 [8]. The SARS-CoV-2 isolate Wuhan NC 045512.2 was used as a reference sequence to analyse mutations at the amino acid level.

2.5. Phylogenetic Analyses

The construction of the phylogenetic tree model was achieved with IQ-TREE v1.6.12 [9], including ModelFinder [10], to select the best nucleotide substitution model (using the Bayesian Information Criterion BIC, the best model was TN+F+I). The visualisation was performed using the iTOL v4 tool [11].

2.6. Selection and Identification of Mutations of Epidemiological Relevance of SARS-CoV-2 and Their Geographical Association

According to the table showing the top 5% mutations observed in the 714 genomes analysed in the Lambayeque region of Peru (Figure 1), a manual search was performed for the missense or non-synonymous mutations of each SARS-CoV-2 gene on the Nextstrain® (https://nextstrain.org/) online server, identifying the lineage to which this genomic sequence belonged. Subsequently, this lineage was verified in Outbreak.info® (https://outbreak.info/) to determine the geographical distribution of said mutations according to the GISAID database. Only those mutations predominant in Peru and mainly in the Lambayeque region were selected for detailed characterisation.

2.7. Ethical Considerations

Ethical approval for sample collection and analysis protocols was granted by the ethics committee of the Almanzor Aguinaga Asenjo Hospital, Chiclayo, Peru, through the ICIS-RPL code 066-DEC-2021. Participation in the study was voluntary, with the signing of an informed consent approved by the same ethics committee; in the case of patients, the consent was signed by their family member or proxy. All information obtained from participants was used only for this research. Therefore, such information will not be stored or used for further studies.

3. Results

Seventy-four nasopharyngeal swab samples were collected from SARS-CoV-2 positive patients (cycle threshold values [CT] obtained by qPCR, ≤28). The Illumina MiSeq sequencing procedure was performed in the Microbial Genomics Laboratory, Department of Cellular and Molecular Sciences Faculty of Sciences and Philosophy (Universidad Peruana Cayetano Heredia, Peru). The viral genomes were assembled by mapping the genome of the Wuhan Hu-1 strain deposited in GenBank as a reference (account number NC_045512).
Seventy-four genomic sequences were obtained from SARS-CoV-2 viral isolates from patients at the Almanzor Aguinaga Asenjo Hospital, Chiclayo, Peru. These sequences were divided into two batches of samples: the first 30 were obtained from COVID-19-positive nasopharyngeal swabs from December 2020 to January 2021. According to the PANGOLIN software v1, the 34 analysed sequences from the first batch were classified within the lineages C.14, C.33, B.1.1.485, B.1.1, B.1.1.1, and B.1.111. The C.14 lineage is the most predominant, at 76.7% (n = 23), and the other lineages range between 3.3–6.7%, respectively.
The second batch of samples consisted of 44 nasopharyngeal swabs positive for COVID-19 from September and October 2021. After analysing these sequences in Nextclade V1.10.0 (https://clades.nextstrain.org/) and PANGOLIN V3.1.16 (https://pangolin.cog-uk.io/), it was determined that the predominant variants in this sampling period were the Delta (72.7%), Gamma (13.6%), Mu (4.6%), and Lambda (9.1%) variants, which were distributed between clades 20J, 21G, 21H, 21J, and 21I.
Another variant of interest is C.37 (Lambda variant); the present investigation determined that 04 sequences (n = 4/44) of the second batch of samples analysed belonged to C.37, classified within clade 21G. C.37 is considered a variant native to Peru, also called the Andean variant; the first reports of C.37 began in Lima, Peru, in approximately August 2020. Subsequently, this variant has been predominant in the sequencing results detected in Peru since its first report, and it has spread to most countries in South America [12,13].
The sequences of the first and second batch of sequenced samples were analysed in the online programme Nextclade V1.10.0 (https://clades.nextstrain.org/) to assign the corresponding clades, identify the mutations in each SARS-CoV-2 gene, and also to perform quality control of these sequences.
Characteristic mutations that have been found in the C.14 lineage include the T1246I and G3278S mutations in the ORFIa gene; P314L in the ORFIb gene; D614G in spike protein; and R203K and G204R in gene N (Table 1).
In addition, in the C.14 lineage, it was possible to identify that the most frequent and relevant amino acid change in the SARS-CoV-2 Nucleocapsid gene was R203K and G204R (N gene, n = 34/35). Likewise, in the case of the changes in nearby nucleotides GGG>AAC at positions 28881–28883, the triplet was also present in n = 34/35 of the sequences analysed.
Using the Nextclade tool of Nextstrain, the Delta variant of SARS-CoV-2 (n = 32) was classified into clades 21J (n = 27) and 21I (n = 5). It was observed that within the Delta variant of SARS-CoV-2, there were several sublineages, among which were AY.26, AY.39.2, AY.100, AY.122, AY.43, AY.102, and B.1.617.2. The characteristic mutations found in the Delta variant and its sublineages are shown in Table 2.
As for the Gamma variant of SARS-CoV-2, according to the Nextclade tool of Nextstrain, the Gamma variant (n = 6) was classified in clade 20J. It was observed that within the Gamma variant, there were two sublineages, among which the following stand out: P.1 (n = 2) and P.1.12 (n = 4). The characteristic mutations found in the Gamma variant and its sublineages are shown in Table 3.
Concerning the Lambda variant of SARS-CoV-2, according to the Nextclade tool of Nextstrain, the Gamma variant (n = 4) was classified in clade 21G. The characteristic mutations found in the Lambda variant (C.37) are shown in Table 4.
A phylogenetic tree of the 714 genomes obtained on 28 April 2022 from the GISAID International Base (www.gisaid.org) of the Lambayeque region, Peru, was created. It can be seen in the phylogenetic tree that the sequences belong mainly to the variants Delta, Omicron, Mu, Lambda, and Gamma of SARS-CoV-2, as shown in Figure 2.
A potency law pattern was recognised in the analysed genomes, and the presence/absence of variants in the 714 sequences is evident. Few variants are widely distributed across genomes, and many are uniquely present in a single genome (Figure 1). This means that only a few variants are present in several genomes, and further analysis can focus on those variants.
The identification of variants of the SARS-CoV-2 genome observed in more than 5% of the genomes analysed in the Lambayeque region, Peru, was also carried out; a specific analysis of these variants is shown in Table 5. The identification of relevant mutations of SARS-CoV-2 and their geographical association was also carried out; for this, a manual review was conducted mutation by mutation in the Nextstrain online tool of those mutations with a frequency greater than 5% to determine where they had been reported according to the GISAID international database. A specific analysis was carried out if Peru emerged as one of the predominant countries for a particular mutation.
The K1817N mutations in the ORF1a gene, S74F in the ORF3a gene, and the P80R mutation in the N gene observed in the SARS-CoV-2 genomes from the Lambayeque region of Peru were selected for analysis because in Peru, a notable increase was found in these mutations during the period from July to August 2021, in contrast to the prevalence worldwide, which remains relatively constant and low.
For each of these mutations, its associated lineage or sublineage was determined; the prevalence statistics of that lineage or sublineage were also investigated worldwide, specifically in Peru. After carrying out all these searches, it was determined that the predominant SARS-CoV-2 sublineages in the 714 sequences analysed were AY.39.2 and P.1.12.1. It should be noted that both sublineages have been reported in various regions of Peru, with Lambayeque being the region with the highest prevalence (Figure 3 and Figure 4).

4. Discussion

Genomic surveillance of SARS-CoV-2 plays a critical role in understanding and responding to the pandemic. Tracking the emergence of mutations and variants through whole genome sequencing enables early detection of novel variants of concern and monitoring their spread, allowing public health officials to implement timely tailored containment measures. The study of the biological and pathogenic properties of new variants, including their transmissibility, virulence, and immune evasion ability, improves existing diagnostics and treatments to ensure they remain effective against new variants. Identifying key mutations correlating with concerning properties provides valuable insights into the virus’s adaptation. Therefore, genomic surveillance is essential to maintain pace with the evolution of SARS-CoV-2 and adapt strategies to combat current and future variants, thus improving our response to the pandemic [6].
The present study allowed us to identify and characterise the genome and lineage of 74 viral strains of the SARS-CoV-2 virus obtained from patients at the Almanzor Aguinaga Asenjo Hospital in Chiclayo, Peru, through next-generation sequencing (NGS) with the Illumina MiSeq system. NGS sequencing allows us to determine the molecular epidemiology of SARS-CoV-2 and learn about the virus’s evolution, transmission, virulence and pathology. Before the arrival of the COVID-19 pandemic, several researchers worldwide began to sequence the complete genome of SARS-CoV-2 to genetically understand this virus, try to elucidate its origin and find a molecular target that serves as a basis for the development of a biological product or vaccine against SARS-CoV-2. The implementation of bioinformatics and genomic tools allowed the active surveillance of SARS-CoV-2, as well as the identification of new lineages and the registration of new mutations in the viral genome, which will allow a better understanding of the evolution and replacement rates of the virus. Several countries in Latin America, Asia, Europe and Africa have published their sequencing results of the whole genome of SARS-CoV-2 in the GISAID international database [14,15].
The predominant lineage in this sampling period was lineage C.14 in 76.7% (n = 23) belonging to clade 20D, which also detected the circulation of lineages C.33, B.1.1.485, B.1.1, B.1.1.1, and B.1.111 in a percentage ranging between 3.3% and 6.7%, respectively, distributed between clades 20B and 20A. The results agreed with those reported by [4], who sequenced five genomes obtained from patients from the Lambayeque region at the end of April 2020, reporting the circulation of lineage B.1.1.1, which was classified according to Nextclade of Nextstrain in clade 20B. Also, our results agree with what was described by [3]. These authors indicate that SARS-CoV-2 isolates during the initial period of the pandemic in Peru belong or are grouped mainly in clade 20B; this clade is very characteristic of isolates obtained from patients with COVID-19 in the European continent. Likewise, these authors identified nine predominant lineages: A.1, A.2, A.5, B.1, B.1.1, B.1.1.1, B.1.5, B.1.8, and B.2, with the most predominant being B.1 and B.1.1.1.
According to the analyses carried out in GISAID and Pangolin [16] (B1.1.1, B1.5), the results highlight that most Peruvian SARS-CoV-2 sequences are classified within clade B.1 and within subclade B.1.1.1. The results described above differ from our results because the predominant lineage in our first sequences was C.14; this difference in results can be attributed to the sampling period in which the samples of nasopharyngeal swabs positive for COVID-19 were collected. Although there are no reports of the C.14 lineage, the GISAID-enabled outbreak.info mutation tracker (https://outbreak.info/situation-reports?pango=C.14) indicates that this lineage has been reported in the following countries: Peru (93.0%), United States (2.0%), Japan (2.0%), Democratic Republic of Congo (1.0%), and Brazil (1.0%) and was first reported on 20 March 2020.
The study reports that the principal variant detected in the second batch was B.1.617.2 (Delta), which has greater transmissibility virulence and can cause cases of reinfection and outbreaks due to the presence of a high number of mutations in the spike protein that allow more significant resistance to the action of antibodies or immune escape. The Delta variant has been reported in several countries worldwide and can replace other regional variants in circulation. The high infectivity of Delta is linked to its high viral load and the short incubation period until the appearance of symptoms. The study also reports that the Delta variant has been found to have immune evasion in patients who received doses of Pfizer® (Pfizer Inc., New York, NY, USA), Moderna® (Moderna, Inc. Cambridge, MA, USA), and Covax®, suggesting that the variant may require updated vaccines to provide better protection [17,18].
In addition, our sequences assigned or classified as Delta variant (21I, 21J) presented various mutations in the spike gene (S gene) such as L452R, T478K, D614G, and P618R. These mutations have been reported worldwide by various researchers, and indicate that they provide biological advantages, among which are: an increase in binding to the ACE-2 receptor, increased transmissibility, risk of hospitalisation, and immune escape or resistance to specific antibodies [19,20]. Some reports indicate that the Delta variant has undergone another mutation, K417N, T95I, and W258L, calling it the Delta Plus variant; however, in our results, we did not find this mutation in any of the analysed sequences belonging to the Delta variant. Some research indicates that the K417N, T95I, and W258L mutation of the spike protein increases the viral ability to achieve immune evasion; however, little is still known about the pathogenicity and virulence of this new variant of SARS-CoV-2 [21].
The study identified that the most frequent and relevant mutation in the SARS-CoV-2 spike gene of the C.14 lineage was D614G, which is associated with more significant pathogenesis and virulence, and evidence suggests that it can improve the transmission of the virus by increasing the amount of viral load in the upper respiratory tract. The sequences analysed and classified as C.37 contain a characteristic deletion in the gene S and mutations not synonymous in the gene spike. These could provide biological advantages such as increased transmissibility, virulence, viral invasion into host cells, and immune escape properties. The Gamma variant of SARS-CoV-2 was also detected in the study, which has lineage-defining mutations, including K417T, E484K, and N501Y, and mutations that allow this variant to increase ACE-2 receptor binding affinity, cause reinfection, increased transmissibility, higher viral load, and immune evasion [22]. The study also cites reports that patients infected with SARS-CoV-2 who carried the D614G mutation developed a moderate/severe COVID-19 condition, while patients infected with SARS-CoV-2 who did not carry this mutation developed mild symptoms and that the Gamma variant was the predominant lineage in the second wave of COVID-19 cases in Brazil.
During the third wave of the COVID-19 pandemic in Peru, the Delta and Gamma variants predominated until the emergence of the Omicron variant. Studies from other countries, such as one from Pakistan, report that the Delta, Beta, and Gamma variants had specific mutations that could provide various biological advantages to these variants. The simultaneous coexistence of highly transmissible SARS-CoV-2 variants could lead to evolutionary competition, where specific variants with mutations that improve their infectious capacity compete with others characterised by their immunological evasion capacity. In Peru, the Lambda variant (C.37) became the predominant variant in the coastal and Andean region, surpassing other circulating Variants of Concern (VOC) such as Gamma and Delta, despite Gamma having a higher frequency during the second wave in the Northwest region due to its proximity to Brazil (Vargas-Herrera et al., 2022a). The Lambda variant (identified as C.37) was initially identified in Peru in August 2020 where it became the predominant variant in early 2021. While it was detected to have a significant prevalence in Peru, its global spread was more limited than other variants, such as Delta and Omicron [23]. There are several possible explanations for Lambda’s initial prevalence in Peru and its lack of being widespread.
Local factors such as population density, mobility patterns, and social behaviours during restriction periods in Peru may have favoured Lambda variant transmission. It could also compete with other variants present in other regions with more extraordinary transmission and immune evasion capacity. Another factor was Peru’s effective implementation of public health measures that limited its expansion. Likewise, the Lambda variant probably needed the optimal characteristics to spread efficiently globally, given that its genetic properties and ability to infect cells and infect from person to person may have been less favourable than other variants. More research is required to fully understand the behaviour of this variant and why it spreads differentially in some regions.
The molecular mechanisms involved in the genetic variability of SARS-CoV-2 are multiple, and substitutions, insertions, deletions, and even genetic recombination mutations have been reported. It is known that the virus regions with the essential mutations reside in the spike protein (S) and the non-structural proteins (NSPs) such as ORF3b, ORF6, ORF7a, etc., which antagonise interferon signalling through different mechanisms and play a crucial role in altering the virus phenotype [24].
Some mutations, such as K417T, N501Y, and E484K, have been associated with higher affinity binding to the ACE2 receptor and resistance to neutralising antibodies, which could affect the infectivity and immune evasion of the virus. In addition, the error-prone replication process of the virus generates a broad spectrum of mutations, and dominant variants exhibit rapid propagation, contributing to the virus’s evolutionary dynamics and subsequent distribution [25]. The presence of specific mutations, such as Spike_D614G, has also been related to increased infectivity and transmissibility of the virus. These mutations can confer advantages by improving the ability of the virus to bind to host cells and evade immune responses, which could influence the appearance of mutations in virus variants [26].
It is known that the Beta (B.1.351) and Gamma (P.1) variants of SARS-CoV-2 present deletions in the spike protein region, such as the 242–244 deletion. These deletions can alter the virus’s interaction with cellular receptors and immune response. On the other hand, the Delta (B.1.617.2) variant shows the insertion of the L452R mutation in the spike protein, which seems to affect its interaction with ACE2 receptors and potentially increase its transmissibility. The Omicron (B.1.1.529) variant stands out for having many mutations, including insertions, deletions and multiple changes in the spike protein and other areas of the viral genome. An insertional genetic sequence in Omicron (insertion: ins214EPE) has even been detected that is not present in earlier versions of SARS-CoV-2. Still, it is common in other viruses, including those causing the common cold and in the human genome. By inserting this fragment into itself, Omicron could mimic the host and better evade the human immune system response, contributing to its transmissibility and the possibility of causing mild or asymptomatic disease [27].
In our results, we can also observe that the N gene of SARS-CoV-2 presents several amino acid mutations that confer various changes or biological advantages. Worldwide, several reports indicate that mutations in the N gene of SARS-CoV-2 reduce the sensitivity of molecular tests (RT-PCR) for the detection of SARS-CoV-2, thus causing the appearance of false negative results for this gene. The N gene of SARS-CoV-2 is of vital importance in the structure and viral cycle, as it is involved in viral assembly, replication and the immune response of the host; also, this SARS-CoV-2 gene is a gene not conserved due to its mutation rate. All these characteristics described above make this gene a target or target to update tests that allow viral diagnosis through and for the development of vaccines [28,29].
The natural evolution of SARS-CoV-2 has led to the emergence of multiple genetic variants with various biological properties, including increased transmission, immune escape, infectivity, and lethality. The initiation of mass vaccination could also be associated with an increase in selective pressure, leading to the appearance of escape mutants. Large-scale whole genome sequencing of SARS-CoV-2 is vital to track the spread of the virus, study local outbreaks, and identify critical mutations in SARS-CoV-2 genes. However, sharing sequencing results in the GISAID database is crucial for almost real-time genomic surveillance worldwide, providing a better understanding of the transmission and viral evolution dynamics of SARS-CoV-2 [30,31].
Finally, the importance of recognizing a power law pattern in the analysed genomes is fundamental as it allows us to identify critical mutations in the SARS-CoV-2 genes. This pattern shows the presence/absence of variants in the 714 genomes analysed. This pattern indicates that certain variants are widely distributed among the analysed genomes and many other variants are uniquely present in a single genome. This means that only a few variants are present in several genomes, and subsequent analysis can focus on those variants of epidemiological interest [7].
Although this study provides valuable information, it has some limitations related mainly to the sample size and the need for clinical data. Thus, only 74 viral sequences were analysed in two groups, which may limit the ability to detect some low-frequency circulating variants. The lack of clinical and epidemiological data associated with the cases analysed could be considered a limitation since the study focused on the genomic analysis of the samples but needed to report data on the severity of the cases, hospitalisations, contacts, etc. Although it was not the present study’s objective, this information is essential to determine the clinical and epidemiological impact of the detected variants. Lastly, it is necessary to consider a possible geographic bias since the study focused on a single city in the region; therefore, it may have captured only some of the diversity of variants circulating in other areas. Future genomic surveillance studies should ideally include representative samples from the entire region. Genomic surveillance in Peru, as in other countries of Latin America, has been vital in the understanding of the COVID-19 pandemic evolution during these almost four years [28,29,30].
Finally, regarding the Lambda variant, it is important to note that the global spread and dominance of specific SARS-CoV-2 variants can be influenced by various factors, including local epidemiological conditions, public health measures, vaccination rates, and the interplay between the virus and the host population. There are several potential reasons why the Lambda variant may have been predominant in Peru but not as widespread globally: (1) Local Factors: Regional differences in population density, healthcare infrastructure, and public health responses can impact the transmission and prevalence of a particular variant. Factors such as crowded living conditions, healthcare capacity, and adherence to preventive measures may contribute. (2) Competing Variants: other variants of concern, like Alpha, Beta, Delta, or others, may have outcompeted Lambda in other regions due to differences in transmissibility or immune escape properties during that time. (3) Vaccination Status: The spread of variants can also be influenced by the level of population immunity through vaccination. If other variants with a competitive advantage emerged in regions with higher vaccination coverage, they might have become more prevalent. (4) Human Mobility: International travel patterns and restrictions can impact the spread of variants. If a particular variant is prevalent in a country with limited international travel or strict entry protocols, its global dissemination may be limited. (5) Adaptation to Host Population: The success of a variant may also be influenced by its ability to adapt to the local population’s genetic and immune characteristics. It is essential to keep in mind that the situation with SARS-CoV-2 is dynamic, and ongoing research is continuously providing new insights. Monitoring the spread and impact of different variants is crucial for understanding the evolving nature of the virus and adapting public health strategies accordingly, including the importance of genomic surveillance and epidemiology [12,13,28,29,30], then continuing to make a molecular diagnosis of COVID-19 as a standard rule [31].
Genomic surveillance is a powerful tool for monitoring and understanding the dynamics of infectious diseases, enabling a more proactive and effective response to emerging threats, such as SARS-CoV-2 and future pandemic pathogens.
In conclusion, this study of SARS-CoV-2 genomes in Chiclayo, Peru, highlights the presence of multiple lineages and variants of concern circulating in the region. The emergence of Delta as the most common variant in later samples from 2021, along with other variants like Gamma, Mu, and Lambda, is particularly alarming due to their potential increased transmission and virulence and their possible ability to escape the immune response. The use of whole genome sequencing and data sharing in databases like GISAID is crucial for understanding the evolution and epidemiology of SARS-CoV-2, which can inform response measures and aid in detecting emerging strains. The findings underscore the importance of genomic surveillance in tracking the spread of SARS-CoV-2 variants and developing tailored public health strategies to limit their transmission. Continued monitoring and sequencing efforts are necessary to stay ahead of the virus’s evolution and ensure effective pandemic control.

Author Contributions

S.L.A.-M.: Conceptualisation, methodology, formal analysis, investigation, writing—original draft preparation, and writing—review and editing; G.A.S.-P.: Conceptualisation, formal analysis, investigation, writing—original draft preparation, and writing—review and editing; J.A.M.-M.: Conceptualisation, formal analysis, investigation, writing—original draft preparation, and writing—review and editing; P.T.-C.: Conceptualisation, investigation, writing—original draft preparation, and writing review and editing; C.D.-V.: Conceptualisation, methodology, investigation, writing—original draft preparation, and writing—review and editing; F.R.A.-G.: investigation, writing—original draft preparation, and writing—review and editing; D.K.B.-A.: investigation, writing—original draft preparation, and writing—review and editing; A.J.R.-M.: investigation, writing—original draft preparation, and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by CONCYTEC–PROCIENCIA Perú—Proyecto de investigación básica 2019–01 (Contrato N°390–2019-FONDECYT) and Vicerrectorado de Investigación y Posgrado Universidad Nacional Mayor de San Marcos (Código N°19101141). This report is part of a dissertation for an MSc degree in SLAM to Graduate School in Molecular Biology, Universidad Nacional Mayor de San Marcos. Universidad Continental, Huancayo, Peru, covered the APC in this article.

Institutional Review Board Statement

Ethical approval for sample collection and analysis protocols was granted by the ethics committee of the Almanzor Aguinaga Asenjo Hospital, Chiclayo, Peru, through the ICIS-RPL code 066-DEC-2021.

Informed Consent Statement

Participation in the study was voluntary, with the signing of an informed consent approved by the same ethics committee; in the case of patients, the consent was signed by their family member or proxy. All information obtained from participants will be used only for this research. Therefore, such information will not be stored or used for further studies.

Data Availability Statement

Available upon reasonable request.

Conflicts of Interest

A.J.R.-M. has been a consultant/speaker for AstraZeneca, Valneva, and Moderna regarding COVID-19 vaccines and long COVID. The rest of the authors declare no conflicts of interest.

References

  1. Molina-Mora, J.A.; Reales-González, J.; Camacho, E.; Duarte-Martínez, F.; Tsukayama, P.; Soto-Garita, C.; Brenes, H.; Cordero-Laurent, E.; dos Santos, A.R.; Salgado, C.G.; et al. Overview of the SARS-CoV-2 genotypes circulating in Latin America during 2021. Front. Public Health 2023, 11, 1095202. [Google Scholar] [CrossRef]
  2. Chong, Y.M.; Sam, I.C.; Ponnampalavanar, S.; Omar, S.F.S.; Kamarulzaman, A.; Munusamy, V.; Wong, C.K.; Jamaluddin, F.H.; Gan, H.M.; Chong, J.; et al. Complete Genome Sequences of SARS-CoV-2 Strains Detected in Malaysia. Microbiol. Resour. Announc. 2020, 9, e00383-20. [Google Scholar] [CrossRef]
  3. Juscamayta-López, E.; Tarazona, D.; Valdivia, F.; Rojas, N.; Carhuaricra, D.; Maturrano, L.; Gavilán, R. Phylogenomics reveals multiple introductions and early spread of SARS-CoV-2 into Peru. bioRxiv 2020. Available online: https://www.biorxiv.org/content/10.1101/2020.09.14.296814v2 (accessed on 3 October 2021).
  4. Aguilar-Gamboa, F.R.; Salcedo-Mejía, L.A.; Serquén-López, L.M.; Mechan-Llontop, M.E.; Tullume-Vergara, P.O.; Bonifacio-Briceño, J.J.; Salas-Asencios, R.; Díaz, H.S.; Cárdenas, J.P. Genomic Sequences and Analysis of Five SARS-CoV-2 Variants Obtained from Patients in Lambayeque, Peru. Microbiol. Resour. Announc. 2021, 10, e00789-21. [Google Scholar] [CrossRef]
  5. Romero, P.E.; Dávila-Barclay, A.; Salvatierra, G.; González, L.; Cuicapuza, D.; Solis, L.; Marcos-Carbajal, P.; Huancachoque, J.; Maturrano, L.; Tsukayama, P. The Emergence of SARS-CoV-2 Variant Lambda (C.37) in South America. medRxiv 2021. Available online: https://www.medrxiv.org/content/10.1101/2021.06.26.21259487v1 (accessed on 3 October 2021).
  6. Aguilar-Gamboa, F.R.; Suclupe-Campos, D.O.; Vega-Fernández, J.A.; Silva-Diaz, H. Diversidad genómica en SARS-CoV-2: Mutaciones y variantes. Rev. Cuerpo Méd. Hosp. Nac. Almanzor Aguinaga Asenjo 2021, 14, 572–582. Available online: http://cmhnaaa.org.pe/ojs/index.php/rcmhnaaa/article/view/1465/556 (accessed on 14 May 2023). [CrossRef]
  7. Molina-Mora, J.A.; Cordero-Laurent, E.; Godínez, A.; Calderón-Osorno, M.; Brenes, H.; Soto-Garita, C.; Pérez-Corrales, C.; Drexler, J.F.; Moreira-Soto, A.; Corrales-Aguilar, E.; et al. SARS-CoV-2 genomic surveillance in Costa Rica: Evidence of a divergent population and an increased detection of a spike T1117I mutation. Infect. Genet. Evol. 2021, 92, 104872. [Google Scholar] [CrossRef] [PubMed]
  8. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772. [Google Scholar] [CrossRef]
  9. Nguyen, L.T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  10. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates. Nat. Methods 2017, 14, 587. [Google Scholar] [CrossRef] [PubMed]
  11. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019, 47, W256. [Google Scholar] [CrossRef] [PubMed]
  12. Darvishi, M.; Rahimi, F.; Abadi, A.T.B. SARS-CoV-2 Lambda (C.37): An emerging variant of concern? Gene Rep. 2021, 25, 101378. [Google Scholar] [CrossRef] [PubMed]
  13. Padilla-Rojas, C.; Jimenez-Vasquez, V.; Hurtado, V.; Mestanza, O.; Molina, I.S.; Barcena, L.; Ruiz, S.M.; Acedo, S.; Lizarraga, W.; Bailon, H.; et al. Genomic analysis reveals a rapid spread and predominance of lambda (C.37) SARS-CoV-2 lineage in Peru despite circulation of variants of concern. J. Med. Virol. 2021, 93, 6845–6849. [Google Scholar] [CrossRef]
  14. Romero, P.E. Escasa información genómica en bases de datos públicas para investigar el SARS-CoV-2 en Latinoamérica. Rev. Peru. Med. Exp. Salud Publica 2020, 37, 374. [Google Scholar] [CrossRef]
  15. Mahmood, T.B.; Saha, A.; Hossan, M.I.; Mizan, S.; Arman, S.M.A.S.; Chowdhury, A.S. A next-generation sequencing (NGS) analysis to reveal genomic and proteomic mutation landscapes of SARS-CoV-2 in South Asia. Curr. Res. Microb. Sci. 2021, 2, 100065. [Google Scholar] [CrossRef]
  16. Padilla-Rojas, C.; Vega-Chozo, K.; Galarza-Perez, M.; Calderon, H.B.; Lope-Pari, P.; Balbuena-Torres, J.; Neyra, D.G.; Huaringa-Nuñez, M.; Rojas-Serrano, N.; Caceres-Rey, O.A. Genomic analysis reveals local transmission of SARS-CoV-2 in the early pandemic phase in Peru. bioRxiv 2020. [Google Scholar] [CrossRef]
  17. Tareq, A.M.; Emran, T.B.; Dhama, K.; Dhawan, M.; Tallei, T.E. Impact of SARS-CoV-2 delta variant (B.1.617.2) in surging second wave of COVID-19 and efficacy of vaccines in tackling the ongoing pandemic. Hum. Vaccines Immunother. 2021, 17, 4126–4127. [Google Scholar] [CrossRef]
  18. Luo, C.H.; Morris, C.P.; Sachithanandham, J.; Amadi, A.; Gaston, D.; Li, M.; Swanson, N.J.; Schwartz, M.; Klein, E.Y.; Pekosz, A.; et al. Infection with the SARS-CoV-2 Delta Variant is Associated with Higher Infectious Virus Loads Compared to the Alpha Variant in both Unvaccinated and Vaccinated Individuals. medRxiv 2021. [Google Scholar] [CrossRef]
  19. Liu, Y.; Liu, J.; Johnson, B.A.; Xia, H.; Ku, Z.; Schindewolf, C.; Widen, S.G.; An, Z.; Weaver, S.C.; Menachery, V.D. Delta spike P681R mutation enhances SARS-CoV-2 fitness over Alpha variant. bioRxiv 2021. [Google Scholar] [CrossRef] [PubMed]
  20. Cherian, S.; Potdar, V.; Jadhav, S.; Yadav, P.; Gupta, N.; Das, M.; Rakshit, P.; Singh, S.; Abraham, P.; Panda, S.; et al. SARS-CoV-2 Spike Mutations, L452R, T478K, E484Q and P681R, in the Second Wave of COVID-19 in Maharashtra, India. Microorganisms 2021, 9, 1542. [Google Scholar] [CrossRef] [PubMed]
  21. Rahimi, F.; Talebi Bezmin Abadi, A. Emergence of the Delta Plus variant of SARS-CoV-2 in Iran. Gene Rep. 2021, 25, 101341. [Google Scholar] [CrossRef] [PubMed]
  22. Nonaka, C.K.V.; Gräf, T.; de Lorenzo Barcia, C.A.; Costa, V.F.; de Oliveira, J.L.; da Hora Passos, R.; Bastos, I.N.; de Santana, M.C.B.; Santos, I.M.; de Sousa, K.A.F.; et al. SARS-CoV-2 variant of concern P.1 (Gamma) infection in young and middle-aged patients admitted to the intensive care units of a single hospital in Salvador, Northeast Brazil, February 2021. Int. J. Infect. Dis. 2021, 111, 47–54. [Google Scholar] [CrossRef] [PubMed]
  23. Vargas-Herrera, N.; Araujo-Castillo, R.V.; Mestanza, O.; Galarza, M.; Rojas-Serrano, N.; Solari-Zerpa, L. SARS-CoV-2 Lambda and Gamma variants competition in Peru, a country with high seroprevalence. Lancet Reg. Health Am. 2022, 6, 100112. [Google Scholar] [CrossRef] [PubMed]
  24. Leelawong, M.; Mitchell, S.L.; Fowler, R.C.; Gonzalez, E.; Hughes, S.; Griffith, M.P.; Marsh, J.W.; Harrison, L.H.; Rakeman, J.L. SARS-CoV-2 N gene mutations impact detection by clinical molecular diagnostics: Reports in two cities in the United States. Diagn. Microbiol. Infect. Dis. 2021, 101, 115468. [Google Scholar] [CrossRef]
  25. Lee, S.; Won, D.; Kim, C.K.; Ahn, J.; Lee, Y.; Na, H.; Kim, Y.T.; Lee, M.K.; Choi, J.R.; Lim, H.S.; et al. Novel indel mutation in the N gene of SARS-CoV-2 clinical samples that were diagnosed positive in a commercial RT-PCR assay. Virus Res. 2021, 297, 198398. [Google Scholar] [CrossRef] [PubMed]
  26. Chiara, M.; D’Erchia, A.M.; Gissi, C.; Manzari, C.; Parisi, A.; Resta, N.; Zambelli, F.; Picardi, E.; Pavesi, G.; Horner, D.S.; et al. Next generation sequencing of SARS-CoV-2 genomes: Challenges, applications and opportunities. Brief. Bioinform. 2021, 22, 616–630. [Google Scholar] [CrossRef]
  27. Álvarez-Díaz, D.A.; Laiton-Donato, K.; Franco-Muñoz, C.; Mercado-Reyes, M. Secuenciación del SARS-CoV-2: La iniciativa tecnológica para fortalecer los sistemas de alerta temprana ante emergencias de salud pública en Latinoamérica y el Caribe. Biomedica 2020, 40 (Suppl. S2), 188. [Google Scholar] [CrossRef]
  28. Rodriguez-Morales, A.J.; Balbin-Ramon, G.J.; Rabaan, A.A.; Sah, R.; Dhama, K.; Paniz-Mondolfi, A.; Pagliano, P.; Esposito, S. Genomic Epidemiology and its importance in the study of the COVID-19 pandemic. Infez. Med. 2020, 28, 139–142. [Google Scholar] [PubMed]
  29. Rabaan, A.A.; Al-Ahmed, S.H.; Sah, R.; Al-Tawfiq, J.A.; Haque, S.; Harapan, H.; Arteaga-Livias, K.; Aldana, D.K.B.; Kumar, P.; Dhama, K.; et al. Genomic Epidemiology and Recent Update on Nucleic Acid-Based Diagnostics for COVID-19. Curr. Trop. Med. Rep. 2020, 7, 113–119. [Google Scholar] [CrossRef]
  30. Rodriguez-Morales, A.J.; Rodriguez-Morales, A.G.; Méndez, C.A.; Hernández-Botero, S. Tracing New Clinical Manifestations in Patients with COVID-19 in Chile and Its Potential Relationship with the SARS-CoV-2 Divergence. Curr. Trop. Med. Rep. 2020, 7, 75–78. [Google Scholar] [CrossRef]
  31. Nawab, M.; Riaz, S.K.; Ismail, E.; Ahamed, A.; Tariq, A.; Malik, M.F.A.; Qusty, N.F.; Bantun, F.; Slama, P.; Umair, M.; et al. Integrated approach for detection of SARS-CoV-2 and its variant by utilizing LAMP and ARMS-PCR. Ann. Clin. Microbiol. Antimicrob. 2024, 23, 11. [Google Scholar] [CrossRef]
Figure 1. Law of potency of mutations observed in the analysis of 714 genomes from the Lambayeque region, Peru. The presence/absence of variants in the 714 genomes analysed is evident. Few variants are widely distributed among genomes, and many are uniquely present in a single genome.
Figure 1. Law of potency of mutations observed in the analysis of 714 genomes from the Lambayeque region, Peru. The presence/absence of variants in the 714 genomes analysed is evident. Few variants are widely distributed among genomes, and many are uniquely present in a single genome.
Tropicalmed 09 00046 g001
Figure 2. A phylogenetic tree was created in IQ-TREE v1.6.12 of the 714 genomes of the Lambayeque region, Peru (until 28 April 2022). The genomes are classified within the variants Mu, Delta, Gamma, Omicron, and Lambda.
Figure 2. A phylogenetic tree was created in IQ-TREE v1.6.12 of the 714 genomes of the Lambayeque region, Peru (until 28 April 2022). The genomes are classified within the variants Mu, Delta, Gamma, Omicron, and Lambda.
Tropicalmed 09 00046 g002
Figure 3. Epidemiological and genomic characteristics of the AY.39.2 lineage (A). The complete SARS-CoV-2 genome sequence, including the genes identified with different colours, is represented. The circles represent mutations of the AY.39.2 lineage. The lineage is characterised by various mutations distributed in the ORF1ab genes, S gene, ORF3a, and N gene. (B) The AY.39.2 lineage has been found in various regions of Peru, with Lambayeque being the region with the highest lineage prevalence. Image obtained from https://outbreak.info/situation-reports?xmin=2023-07-27&xmax=2024-01-27&pango=AY.39.2#geographic (accessed on 1 November 2023).
Figure 3. Epidemiological and genomic characteristics of the AY.39.2 lineage (A). The complete SARS-CoV-2 genome sequence, including the genes identified with different colours, is represented. The circles represent mutations of the AY.39.2 lineage. The lineage is characterised by various mutations distributed in the ORF1ab genes, S gene, ORF3a, and N gene. (B) The AY.39.2 lineage has been found in various regions of Peru, with Lambayeque being the region with the highest lineage prevalence. Image obtained from https://outbreak.info/situation-reports?xmin=2023-07-27&xmax=2024-01-27&pango=AY.39.2#geographic (accessed on 1 November 2023).
Tropicalmed 09 00046 g003
Figure 4. Epidemiological and genomic characteristics of the P.1.12.1 lineage (A). The complete SARS-CoV-2 genome sequence, including the genes identified with different colours, is represented. Table 1. 12.1 lineage. The lineage is characterised by various mutations distributed in the ORF1ab genes, S gene, ORF3a, and N gene. (B) The P.1.12.1 lineage has been found in various regions of Peru, with Lambayeque being the region with the highest lineage prevalence. Image obtained from https://outbreak.info/situation-reports?xmin=2023-07-27&xmax=2024-01-27&pango=P.1.12.1 (accessed on 1 November 2023).
Figure 4. Epidemiological and genomic characteristics of the P.1.12.1 lineage (A). The complete SARS-CoV-2 genome sequence, including the genes identified with different colours, is represented. Table 1. 12.1 lineage. The lineage is characterised by various mutations distributed in the ORF1ab genes, S gene, ORF3a, and N gene. (B) The P.1.12.1 lineage has been found in various regions of Peru, with Lambayeque being the region with the highest lineage prevalence. Image obtained from https://outbreak.info/situation-reports?xmin=2023-07-27&xmax=2024-01-27&pango=P.1.12.1 (accessed on 1 November 2023).
Tropicalmed 09 00046 g004
Table 1. Mutations found in the C.14 SARS-CoV-2 lineage of patients from the Lambayeque region, Peru.
Table 1. Mutations found in the C.14 SARS-CoV-2 lineage of patients from the Lambayeque region, Peru.
Lineage C.14Genes Affected by Mutations
ORF1aORF1bSORF3aORF9bN
C.14P2144L
T1246I
G3278S
P2685T
P314L
S638I
H1087Y
V2073L
A222V
D253E
D614G
L101F
L140F
S171L
V225F
T83IH145Y
R203K
G204K
Table 2. Mutations found in the sublineages of the Delta SARS-CoV-2 variant of patients from the Lambayeque region, Peru.
Table 2. Mutations found in the sublineages of the Delta SARS-CoV-2 variant of patients from the Lambayeque region, Peru.
Gene Sublineages of Delta Variant SARS-CoV-2
AY.26AY.39.2AY.122AY.100AY.43AY.102B.1.617.2
ORF1aP1640L
A3209V
V3718A
T3750I
E743D
A1306S
K1817N
P2046L
P2287S
V2930L
T3255I
T3646A
K261N
A1306S
P2046L
P2287S
V2930L
T3255I
T3646A
T403I
A1306S
P2046L
P2287S
V2930L
T3255I
T3646A
A1306S
P2046L
P2287S
V2930L
T3255I
T3646A
A1306S
P2046L
P2287S
V2930L
T3255I
T3646A
A1306S
T3255I
T3646A
ORF1bP314L
G662S
P1000L
P314L
G662S
P1000L
A1918V
Q2635H
P314L
G662S
P1000L
A1918V
P314L
G662S
P1000L
A1219S
A1918V
P314L
G662S
L829I
P1000L
A1918V
P314L
G662S
P1000L
A1918V
P314L
G662S
P1000L
A1918V
ST19R
R158G
Δ156/157
A222V
L452R
T478K
D614G
P681R
D950N
V1264L
T19R
R158G
Δ156/157
L452R
T478K
D614G
P681R
D950N
K1073N
T19R
R158G
Δ156/157
L452R
T478K
D614G
P681R
D950N
T19R
R158G
Δ156/157
L452R
T478K
D614G
P681R
G769V
D950N
T19R
R158G
Δ156/157
L452R
T478K
D614G
P681R
D950N
T19R
R158G
Δ156/157
L452R
T478K
D614G
P681R
D950N
T19R
R158G
Δ156/157
L452R
T478K
D614G
P681R
D950N
ORF3aS26LS26LS26LS26LS26L
T34A
S26LS26L
MI82TI82TI82TI82TI82TI82TI82T
ORF6K48N------------------------
ORF7aV82A
T120I
V71I
V82A
T120I
V82A
T120I
V82A
T120I
V82A
T120I
V82A
T120I
V82A
T120I
ORF7b----T40IT40IT40IT40IT40IT40I
ORF 8S84L
Δ119/120
S84L
Δ119/120
Δ119/120Δ119/120Δ119/120Δ119/120Δ119/120
ND63G
R203M
D377Y
D63G
R203M
G215C
D377Y
D63G
R203M
G215C
D377Y
D63G
R195K
R203M
G215C
D377Y
Q9L
D63G
R203M
G215C
D377Y
D63G
R203M
G215C
D377Y
D63G
R203M
G215C
D377Y
Table 3. Mutations found in the sublineages of the SARS-CoV-2 Gamma variant of patients from the Lambayeque region, Peru.
Table 3. Mutations found in the sublineages of the SARS-CoV-2 Gamma variant of patients from the Lambayeque region, Peru.
Gamma Variant SublineagesGenes Affected by Mutations
ORF1aORF1bSORF3aORF8N
P.1.12S1118L
K1795Q
Δ3675/3677
P314L
E1264D
L18F
T20N
P26S
D138Y
R190S
K417T
N501Y
D614G
H655Y
T1027I
V1176F
S253PE92KP80R
P.1
S1118L
K1795Q
Δ3675/3677

P314L
E1264D
L18F
T20N
P26S
D138Y
R190S
K417T
E484K
N501Y
D614G
H655Y
T1027I
V1176F
S253PE92KP80R
R203K
G204R
Table 4. Mutations found in the Lambda (C.37) SARS-CoV-2 variant from patients in the Lambayeque region, Peru.
Table 4. Mutations found in the Lambda (C.37) SARS-CoV-2 variant from patients in the Lambayeque region, Peru.
Lambda VariantGenes Affected by Mutations
ORF1aORF1bSORF3aORF9bMN
C.37T1246I
P1659T
P2287S
F2387V
P2483S
L3201P
T3255I
G3278S
A3620V
Δ3675/3677
S59F
P314L
T1137I
A1643V
Y1784C
K2385E
K2674R
L5F
G75V
Δ246–252
L452Q
A475V
E484K
P499R
N501T
D614G
H655Y
P681R
T859N
P240HP10SI82TP13L
R203K
G204R
G214C
Table 5. SARS-CoV-2 genome variants were observed in more than 5% of genomes analysed from the Lambayeque region, Peru.
Table 5. SARS-CoV-2 genome variants were observed in more than 5% of genomes analysed from the Lambayeque region, Peru.
Mutation NumberPOSREFALTTotal Number of Genomes with MutationClass of MutationEffectGeneTranscriptAASequence in TranscriptSequence ProteinPatterns in the World
533037CT712SynonymousLowORF1abc.2772C>Tp.Phe924Phe2772/21291924/7096N/A
49023403AG712MissenseModerateSc.1841A>Gp.Asp614Gly1841/3822614/1273Disseminated in the world
23610029CT518MissenseModerateORF1abc.9764C>Tp.Thr3255Ile9764/212913255/7096Disseminated in the world
37315451GA330MissenseModerateORF1abc.15187G>Ap.Gly5063Ser15187/212915063/7096Disseminated in the world
54125469CT328MissenseModerateORF3ac.77C>Tp.Ser26Leu77/82826/275Disseminated in the world
64428461AG328MissenseModerateNc.188A>Gp.Asp63Gly188/126063/419Disseminated in the world
2138986CT309SynonymousLowORF1abc.8721C>Tp.Asp2907Asp8721/212912907/7096N/A
2169053GT309MissenseModerateORF1abc.8788G>Tp.Val2930Leu8788/212912930/7096Disseminated in the world
24411332AG309SynonymousLowORF1abc.11067A>Gp.Val3689Val11067/212913689/7096N/A
974181GT308MissenseModerateORF1abc.3916G>Tp.Ala1306Ser3916/212911306/7096Disseminated in the world
62428311CT200MissenseModerateNc.38C>Tp.Pro13Leu38/126013/419Disseminated in the world
914002CT194MissenseModerateORF1abc.3737C>Tp.Thr1246Ile3737/212911246/7096Very little disseminated in the world
1575716GT121MissenseModerateORF1abc.5451G>Tp.Lys1817Asn5451/212911817/7096Very little disseminated in the world
2269867TC115MissenseModerateORF1abc.9602T>Cp.Leu3201Pro9602/212913201/7096Very little disseminated in the world
2259857CT111SynonymousLowORF1abc.9592C>Tp.Leu3198Leu9592/212913198/7096N/A
50825000CT87SynonymousLowSc.3438C>Tp.Asp1146Asp3438/38221146/1273N/A
56425584CT87SynonymousLowORF3ac.192C>Tp.Thr64Thr192/82864/275N/A
1375386TG86SynonymousLowORF1abc.5121T>Gp.Ala1707Ala5121/212911707/7096N/A
25911537AG86MissenseModerateORF1abc.11272A>Gp.Ile3758Val11272/212913758/7096Disseminated in the world
33813195TC86SynonymousLowORF1abc.12930T>Cp.Val4310Val12930/212914310/7096N/A
60426270CT86MissenseModerateEc.26C>Tp.Thr9Ile26/228Set-75Disseminated in the world
40617259GT72MissenseModerateORF1abc.16995G>Tp.Glu5665Asp16995/212915665/7096Very little disseminated in the world
1535648AC71MissenseModerateORF1abc.5383A>Cp.Lys1795Gln5383/212911795/7096Very little disseminated in the world
51425088GT71MissenseModerateSc.3526G>Tp.Val1176Phe3526/38221176/1273Very little disseminated in the world
14733TC70SynonymousLowORF1abc.468T>Cp.Asp156Asp468/21291156/7096N/A
31212778CT70SynonymousLowORF1abc.12513C>Tp.Tyr4171Tyr12513/212914171/7096N/A
34713860CT70SynonymousLowORF1abc.13596C>Tp.Asp4532Asp13596/212914532/7096N/A
64628512CG70MissenseModerateNc.239C>Gp.Pro80Arg239/126080/419Very little disseminated in the world
311048GT66MissenseModerateORF1abc.783G>Tp.Lys261Asn783/21291261/7096Disseminated in the world
47720937GT58SynonymousLowORF1abc.20673G>Tp.Thr6891Thr20673/212916891/7096N/A
59825844CT44MissenseModerateORF3ac.452C>Tp.Thr151Ile452/828151/275Very little disseminated in the world
1455515GT41SynonymousLowORF1abc.5250G>Tp.Val1750Val5250/212911750/7096N/A
56625613CT38MissenseModerateORF3ac.221C>Tp.Ser74Phe221/82874/275Very little disseminated in the world
Disseminated in the world = high proportion of disseminated variants. Very little disseminated in the world = low proportion of disseminated variants.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aguilar-Martinez, S.L.; Sandoval-Peña, G.A.; Molina-Mora, J.A.; Tsukayama-Cisneros, P.; Díaz-Vélez, C.; Aguilar-Gamboa, F.R.; Bonilla-Aldana, D.K.; Rodriguez-Morales, A.J. Genomic and Phylogenetic Characterisation of SARS-CoV-2 Genomes Isolated in Patients from Lambayeque Region, Peru. Trop. Med. Infect. Dis. 2024, 9, 46. https://doi.org/10.3390/tropicalmed9020046

AMA Style

Aguilar-Martinez SL, Sandoval-Peña GA, Molina-Mora JA, Tsukayama-Cisneros P, Díaz-Vélez C, Aguilar-Gamboa FR, Bonilla-Aldana DK, Rodriguez-Morales AJ. Genomic and Phylogenetic Characterisation of SARS-CoV-2 Genomes Isolated in Patients from Lambayeque Region, Peru. Tropical Medicine and Infectious Disease. 2024; 9(2):46. https://doi.org/10.3390/tropicalmed9020046

Chicago/Turabian Style

Aguilar-Martinez, Sergio Luis, Gustavo Adolfo Sandoval-Peña, José Arturo Molina-Mora, Pablo Tsukayama-Cisneros, Cristian Díaz-Vélez, Franklin Rómulo Aguilar-Gamboa, D. Katterine Bonilla-Aldana, and Alfonso J. Rodriguez-Morales. 2024. "Genomic and Phylogenetic Characterisation of SARS-CoV-2 Genomes Isolated in Patients from Lambayeque Region, Peru" Tropical Medicine and Infectious Disease 9, no. 2: 46. https://doi.org/10.3390/tropicalmed9020046

APA Style

Aguilar-Martinez, S. L., Sandoval-Peña, G. A., Molina-Mora, J. A., Tsukayama-Cisneros, P., Díaz-Vélez, C., Aguilar-Gamboa, F. R., Bonilla-Aldana, D. K., & Rodriguez-Morales, A. J. (2024). Genomic and Phylogenetic Characterisation of SARS-CoV-2 Genomes Isolated in Patients from Lambayeque Region, Peru. Tropical Medicine and Infectious Disease, 9(2), 46. https://doi.org/10.3390/tropicalmed9020046

Article Metrics

Back to TopTop