{'fraction_non_alphanumeric': 0.03457953639860019, 'fraction_numerical': 0.019680537749904716, 'mean_word_length': 4.9632231404958675, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 1, 'https://': 7, 'lorem ipsum': 0, 'www.': 8, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Utility of in silico-identified-peptides in spike-S1 domain and nucleocapsid of SARS-CoV-2 for antibody detection in COVID-19 patients and antibody production

SARS-CoV-2 contains four structural proteins, two of which, the spike and nucleocapsid, are commonly used for the standardization of novel methods for antibody detection; however, some limitations in their use have been observed due to the homology of this virus with other phylogenetically-related viruses. We performed in silico analysis to search for novel immunogenic and antigenic peptides. A total of twenty-five peptides were preliminarily selected, located in the 3D structure of both proteins. Finally, eight peptides were selected: one located in the N protein and seven in the S1 domain of the spike protein. Additionally, the localization of selected peptides in 2D structures and possible changes in the sequences of these peptides in SARS-CoV-2 variants of concern were analyzed. All peptides were synthetized in MAP8 format, and recombinant S (trimer and RBD) and N proteins were used as antigens to search for antibodies in serum samples derived from COVID-19 patients, and for antibody response in New Zealand rabbits. Results showed high recognition of the serum derived from COVID-19 patients to all selected peptides; however, only the RBD3 peptide induced antibody production. In conclusion, this work provides evidence for a new strategy in peptide selection and its use for antibody detection or antibody production in animals.


Results
Identification of immunogenic and antigenic peptides using in silico analysis. Using the IEDB analysis resource, we searched for the presence of linear epitopes, sequential epitopes, surface accessibility and antigenicity; complementary analyses were performed in BepiPred-2.0 to predict antigenic peptides. The raw material of all peptides is shown in Supplementary Material S1, whereas in Table 1, a concentrated number of peptides found in each part of the S (S1 and S2) and N proteins is shown. During the first two rounds of selection, we observed that most of the peptides in the S protein were in the S2 domain. Furthermore, 59% were linear epitopes, 42% were sequential epitopes, 55% had surface accessibility and 52%-65% exhibited antigenicity, while the number of peptides found in the S1 NTD and S1 RBD was between 15 and 32% in each category. The S protein is highly antigenic; however, some peptides were shorter than eight amino acids, and during the second round of selection, they were eliminated: eleven in S1 NTD, five in S1 RBD and twenty-six in S2. With respect to the N protein, twelve linear epitopes and sequential epitopes were identified, ten with surface accessibility, fifteen with antigenicity and seven B cell epitopes, with no changes between the first and second rounds of selec-  (Table 1). For the third round of selection, only S1 (without subclassifications) and S2 domains were used considering immunogenicity, antigenicity, and the type of epitope (linear or conformational) for the selection of fifteen peptides in S1, four peptides in S2 and six peptides in N ( Table 2).

Localization of selected peptides in tertiary structures.
To analyze the localization of each peptide in the 3D structure, first, the trimeric structure of the glycosylated spike protein was constructed (Fig. 1A), and the localization of all preselected peptides was determined in the monomeric structure. Most of the peptides were exposed, and none of them were in glycosylation sites (Fig. 1B). The localization of each peptide and the color pattern used are shown in Supplementary Material S2. For the final selection step, we performed conservancy analysis (result of 100%) and protein BLAST using sequences derived from SARS-CoV-2 (Supplementary Material S3), and complementary conservancy analysis were performed considering other Coronavirus-related virus that infect human and animals (Supplementary materials S4 and S5), also, we considered the exposition in 3D structure. We selected seven peptides in the S1 domain, two in the S1 NTD from amino acids 13 to 20 (RTQLPPAYTNS) and 153 to 165 (MESEFRVYSSANN), while the selected peptides in the S1 RBD were located from amino acids 438 to 537. All selected peptides are marked in red in their monomeric structure (Fig. 1C). Additionally, all preselected and selected peptides were marked in the trimeric structure of the spike protein, and similar results were observed (Supplementary Material S6). Later, we analyzed the localization of preselected and selected peptides in the 3D structure of the N protein. Unlike the S protein, most of the peptides in the N protein were found internally, and only the selected peptide was fully exposed and located from amino acids 1 to 16 (MSDNGPQNQRNAPRIT) (Fig. 1D).
Characterization of selected peptides in the secondary structure and SARS-CoV-2 VOC. We analyzed the primary amino acid sequence of both proteins and localized specific domains as receptor binding domains or posttranscriptional modifications (PTMs) such as glycosylation in the spike protein and RNAbinding or dimerization domains in the nucleocapsid protein. Interestingly, analysis of the selected peptides in protein plots or snake diagrams provided information about the localization of selected peptides (extra, transmembrane, or intraviron) in the spike or N protein, as well as PTM, disulfide bonds and signal peptide. During the chemical synthesis of selected peptides, we obtained the recombinant proteins of trimeric spike (S trimer), recombinant RBD (rRBD) and recombinant N (rN), and once the purification procedure was completed, electrophoresis and Coomassie blue staining were performed. The S trimer and rRBD possessed approximate molecular weights (MWs) of 250 kDa and 25 kDa, respectively, whereas rN had an approximate MW of 50 kDa (Fig. 2   www.nature.com/scientificreports/ three amino acid changes; G496S, Q498R and Y505H in the RBD4 peptide sequence. The peptides RBD2 and RBD5 do not contain any amino acid changes or deletions, whereas, the sequence of N peptide presents a single amino acid change in the Alpha variant (D3L) ( Table 3). The eight selected peptides were synthetized in MAP8 format, provided that information regarding the synthesis included HPLC and MS. The length of the peptides was approximately 11-12 amino acids each with a molecular weight (MW) > 10 kDa due to the eight ramifications of the multiantigenic peptide. The predicted molecular weight correlated with the observed MW in MS, and HPLC results showed that MAP8 peptides had an estimated purity of 50-78% (Table 4). Crude data from the MS and HPLC are included in Supplementary Material S9. After synthesis was complete, MAP8 peptides were hydrated in sterile PBS to a final concentration of 1 mg/mL and used as antigens for antibody detection in COVID-19 patients and immune response induction in New Zealand rabbits.

Utility of MAP8 peptides and recombinant proteins for antibody detection in healthy and COVID-19 patients.
To evaluate the potential utility of the eight selected peptides and the full-length proteins for antibody detection by indirect ELISA, we used serum samples derived from healthy and COVID-19 patients (diagnosed by RT-PCR with at least 15 days since the onset of symptoms). The experiment included evaluation of IgM/IgG (whole molecule), IgM and IgG levels using recombinant proteins (S trimer, rRBD and rN), individual peptides (NT-1/2, RBD1-5 and N in MAP8 format) or a mixture of all the peptides found in protein S and all of them (S and N) at a final concentration of 0.1 µg/mL. We observed a high recognition of each individual peptide with at least four times the optical density of difference between healthy and COVID-19 patients in the detection of IgM/IgG. In relation to individual IgM and IgG detection, similar results were obtained, most of which exhibited statistically significant differences ( Immune response induction in rabbits using MAP8 peptides derived from S and N structural proteins. Finally, we selected four peptides (NT-1, RBD3, RBD4 and N) for the immunization of New Zealand rabbits to test the capacity of selected peptides for humoral immune response induction. Over thirty days, three immunizations were performed using 100 µg of each peptide, the mixture of the peptides found in S protein or all peptides. Blood samples were collected before immunizations (pre-immune serum) and 15 (first Table 3. Amino acid changes and deletions in SARS-CoV-2 variants of concern and their localization in the selected peptides. *Amino acid changes and deletions of each variants were collected from the database CoVariants (https:// covar iants. org/ varia nts). Δ: Deletion.  www.nature.com/scientificreports/ evaluation) and 35 (second evaluation) days after the first immunization. In rabbits immunized with the peptides NT-1, RBD4 and N, no antibody response was observed until day 35, and the only individual peptide that was able to induce antibody response was RBD3. Both, the mixture of peptides situated in the S protein and the mixture of all peptides induced an antibody response (Fig. 4).

Discussion
SARS-CoV-2 is responsible for COVID-19 infection; structurally, it is composed of four structural proteins called the spike (S), membrane (M), envelope (E) and nucleocapsid (N). Recombinant proteins of structural N and S (S1 and RBD subunit) proteins are commonly used as antigens in immunoassays for antibody detection, including lateral flow immunoassays (LFAs), enzyme-linked immunosorbent assays (ELISAs) and chemiluminescent immunoassays (CLIAs). However, some doubts regarding the utility of antibody tests have been voiced due to the antigenic similarity between SARS-CoV-2 and other circulating coronaviruses 21,22 . To avoid potential crossreactivity, in silico analysis can be performed for the selection of immunogenic peptides that could be used for the standardization of a novel ELISA for antibody detection or the production of monoclonal and polyclonal antibodies for its use in the standardization of novel methods for viral antigen detection, which was the primary objective of this work.  www.nature.com/scientificreports/ The spike glycoprotein mediates the entry of coronavirus into host cells through angiotensin-converting enzyme 2 (ACE2) recognition, which is an attractive antiviral target and the most common antigen for antibody detection. The S protein is composed of two functional subunits, S1 and S2, the S1 subunit is composed of an N-terminal domain and receptor binding domain (RBD) 5,22 . The nucleocapsid protein physically links the envelope to the + RNA genome and contains N-terminal RNA binding (N-NTD) and C-terminal dimerization domains 23 . Both proteins were considered primary targets for epitope identification in this work.
Several in silico approaches have been performed to identify antiviral peptides 24 , epitopes for the development of vaccine-based peptides or multiepitope-based vaccines or diagnostics 25,26 . In this study, a total of twenty-five peptides were reported, nineteen of which were in the S protein and six of in the N protein. We selected eight peptides, two of which were located in the S1 NTD, five in the RBD and one in N. Interestingly, previous studies reported that each monomer of the spike protein has twenty-two glycosylation sites 27 , and this post-translational modification plays an important role in immune response evasion, function and viral pathobiology 28 . This approach indicates that the selected peptides, in addition to all the characteristics associated with immune response induction, may be used as antigens in the development of immunoassays. Interestingly, the amino acids present in the RBD3 peptide are almost exclusively located in the RBD of SARS-CoV-2 but are not found in SARS-CoV according to previous alignments performed by other authors 29,30 .
When we searched for the presence of mutations in SARS-CoV-2 VOCs and their effect on the sequence of the selected peptides, we found that in six of the eight peptides at least one amino acid change was identified. The RBD4 peptide contains the N501Y amino acid change in the receptor-binding domain of the S protein, which may contribute to the 40-70% increase in transmission observed with the Alpha, Beta, Gamma and Omicron variants 6 . The Delta variant contains amino acid changes in the RBD1 (N440K and G446S) and RBD3 (S477N and T478K) peptides, whereas Omicron has 4 amino acid changes in the RBD4 peptide (G496S, Q498R, N501Y, Y505H). These changes are commonly associated with the prevalence, phylogenetic relation and the neutralizing potency of each variant, for which future studies should be performed 7,31 . The emergence of new SARS-CoV-2 variants is a matter of concern for diagnostic testing, however the use of peptides in MAP8 format could be a novel strategy to solve problems associated with cross reactivity of recombinant proteins. An important consideration is the possibility of obtaining a negative result if only one peptide is used, due to amino acid changes that could affect antibody recognition. Therefore, we propose the mixture of two or more peptides for serosurveillance.
The eight selected peptides were synthetized in MAP8. This format increases the size and structural complexity of the peptide, allowing the exposition of eight epitopes in a single molecule, which would favor the use of low concentrations of antigen during its use for immunoassays (0.1 µg/mL in this study). Another interesting characteristic is that this format stimulates antibody production 32 . Using the aforementioned peptides, their use as antigens for antibody detection in COVID-19 patients using indirect ELISA and humoral immune response induction in rabbits was assessed. The use of peptides found in the N protein is important due to IgG seroconversion appears an average 2 days sooner than the S protein 33 , so the use of this peptide could be helpful in the detection of early immune response products. To evaluate the presence of antibodies, in addition to selected peptides, we also obtained recombinant proteins derived from the S (trimer and RBD) and N proteins. In both antigens, the results were considerably higher in COVID-19-positive patients than in COVID-19-negative patients for IgG and total antibodies. In the case of IgM, no significant difference was observed, which may be due to the time that had passed from the onset of symptoms to serum collection. A complete standardization of this method is still necessary, but preliminary results are worthy of a more in-depth investigation. Previously, antibodies against N and S proteins were tested, and no significant difference between the two proteins was observed for the detection of total antibodies and IgG 34 , which could mean that all reported peptides could be mixed during the development of the immunoassay to increase the number of available epitopes or used for the standardization of ELISA methods for each individual protein.
In relation with the use of peptides as antigens for the standardization of methods for antibody detection, several evidence has been published. Woo et al. 35 , developed an ELISA test for IgM and IgG antibody detection against the Spike and Nucleocapsid protein of SARS-CoV. Whereas, during the COVID-19 pandemic the use of antigenic epitopes derived from the structural proteins of the SARS-CoV-2 (Spike, Nucleocapsid, Membrane and Envelope) have been employed as antigens for the standardization of novel ELISA methods with high sensitivity (91.4%) and specificity (83.7%) 36 . In another study performed by Polvere et al., reported ten predicted peptides derived from the S, M and N proteins of SARS-CoV-2, of which, seven were tested by ELISA the reactivity to a panel of serum samples from twenty-four COVID-19 convalescents patients. Results showed differences in the reactivity to selected peptides ranging from the recognition to only one peptide to all seven peptides 37 . The ELISA method used in this study has several similarities with the aforementioned examples, thus providing evidence about the potential clinical utility of the reported peptides in the standardization of an ELISA method.
Although all peptides could be used as antigens for the search for antibodies against SARS-CoV-2, not all of them could be used for humoral immune response induction. The RBD3 peptide was the only peptide that induced a humoral immune response by itself, and the mixture of peptides could induce an immune response due to the increased complexity and molecular weight of the combined peptides. The RBD3 peptide is highly conserved in VOCs Alpha, Beta, Gamma and Omicron, therefore the serum derived from the immunized rabbits could be used for analysis of neutralizing activity against these VOCs. A notable result was that the RBD4 peptide that contains N501Y (and three additional amino acid changes in the Omicron variant), did not induce humoral immune response in rabbits, but was recognized by the serum samples derived from patients.
Positive immune response induction in rabbits and the recognition of peptides by serum samples derived from COVID-19 patients provide evidence regarding the interaction of selected peptides with MHC molecules, T cell response, and B cell activation. The use of synthetic peptides for immune response induction has the advantage that polyclonal antiserum in theory will be as specific as monoclonal antibody but much easier to produce, and www.nature.com/scientificreports/ synthetic peptides are cheaper and less time-consuming to produce than the synthesis of other antigens, such as recombinant proteins 38 . Some experiments still need to be performed to assess the clinical application of peptides as antigens (for serosurveillance of VOCs) or in the production of antibodies for antigen detection; however, our preliminary results provide not only in silico approaches (localization in 2D and 3D structure, post-translational modifications and analysis of VOC) but also experimental data comparing peptides to the full-length protein.
In addition, we provide a list of seventeen other peptides that could be tested with several applications, such as immunoassay development, antibody production or vaccine candidates. Experimental procedures to be performed include the standardization of ELISA methods using reported peptides for antibody detection in vaccinated or infected patients (with circulating VOC or new variants) or the recognition of native protein using serum derived from immunized rabbits. In conclusion, this work provides evidence regarding the potential utility of MAP8 peptides selected based on in silico analysis of the primary amino acid structure of the spike glycoprotein and nucleocapsid protein of SARS-CoV-2 for the standardization of novel immunoassays for antibody detection in serum samples of COVID-19 patients and vaccinated individuals or for antibody production in animal models that can be used for the detection of antigens in infected individuals.

Materials and methods
In silico analysis. First, we downloaded the amino acid sequences of the S (access number QJX59884.1) and N (access number QJX60131.1) proteins of SARS-CoV-2 from the NCBI protein database (https:// www. ncbi. nlm. nih. gov/). The amino acid sequence of the S protein was divided into two segments: S1 (14-685) and S2 (686-1273). In S1, subdivision was performed as S1 NTD (N-terminal domain) from amino acids 14-305 and S1 receptor binding domain (RBD) from amino acids 319-685. Considering the importance of the S1 domain during SARS-CoV-2 infection, immune response, and diagnosis, most of our efforts were allocated to analyzing this specific region, although our analysis also included the S2 domain.
Using the aforementioned sequences, in silico analysis was performed using the B cell epitope prediction tools of the IEDB analysis resource, specifically the methods for predicting continuous antibody epitopes (http:// tools. iedb. org/ bcell/) from protein sequences as follows: Bepipred-1.0 (linear epitope prediction) and 2.0 (sequential B cell epitope predictor), Emini surface accessibility scale and Kolaskar and Tongaonkar antigenicity scale. Additionally, complementary analysis was performed using BepiPred-2.0 (http:// www. cbs. dtu. dk/ servi ces/ BepiP red/ index. php) to search for B cell epitopes and predicted antigenic peptides (http:// imed. med. ucm. es/ Tools/ antig enic. pl) to search for antigenic peptides. All predicted peptides are shown in Supplementary Material S1, and the number of obtained peptides is shown in Table 1. For the second selection, all peptides with a length < 8 amino acids were filtered and deleted, and we considered the presence of at least two different properties (accessibility, antigenicity, or immunogenicity) for preselection. Finally, we performed conservancy analysis using the IEBD analysis resource (http:// tools. iedb. org/ conse rvancy/), nucleotide BLAST (https:// blast. ncbi. nlm. nih. gov/ Blast. cgi) and considered the localization of each peptide in the primary, secondary and tertiary structures. PDB accession, 3D structure modeling and visualization of S and N proteins. The PDB of fulllength SARS-CoV-2 S protein in the open state and glycans was obtained from The Amaro Lab (https:// amaro lab. ucsd. edu/ covid 19. php), which was manually renumbered by eliminating the membrane, water molecules and ions. The 3D structure of the complete N protein was predicted for homology modeling using the Robetta program by Baker Lab, and PDB 7ACS and 6YUN (https:// robet ta. baker lab. org/) were used as templates. All 3D structures were visualized using the VMD 1.9.1 program (Visual Molecular Dynamics), which is a molecular visualization program for displaying, animating, and analyzing large biomolecular systems using 3D graphics with built-in scripting (http:// www. ks. uiuc. edu/ Resea rch/ vmd/).

Schematic representations of N and S protein.
The schematic representation (at scale) of the amino acid sequences of S and N proteins, including domains and localization, was illustrated using the Illustrator of Biological Sequences (IBS) 39 . The 2D schematic representation of both proteins was created in Protter (https:// wlab. ethz. ch/ prott er/ start/), an interactive and customizable web-based application that enables the integration and visualization of both annotated and predicted protein sequence features together with experimental proteomic evidence for peptides and PTMs, onto the transmembrane topology of a protein 40 . For both illustrations, we used UniProt protein accession P0DTC2 for the S protein and P0DTC9 for the N protein.
Purification of SARS-CoV-2 S, RBD and N proteins. The SARS-CoV-2 S HexaPro plasmid (Addgene Cat#154754) was transfected into HEK293F cells, a human embryonic kidney cell line, and the cells were grown in 293FreeStyle expression medium (Life Technologies Cat#12338002) and cultured at 37 °C with 5% CO 2 . The supernatant was recovered, centrifuged, filtered and purified by immobilized metal affinity chromatography (IMAC) under native conditions. The resulting protein had six prolines substituted, a substitution of "GSAS" (furin cleavage site) and the addition of a C-terminal foldon trimerization motif with an eight histidine tag that generates a prefusion-stabilized SARS-CoV-2 S protein. The SARS-CoV-2 receptor binding domain (RBD) was amplified from a template of cDNA derived from a clinical isolation of SARS-CoV-2 and cloned into pLATE-31 according to the protocol of the aLICator Ligation Independent Cloning and Expression System (Cat# K1261 Thermo Scientific). The cDNA of SARS-CoV-2 protein N was also obtained from clinical isolation and cloned into pLATE-51 (Cat# K1251 Thermo Scientific). pLATE-RBD with C-terminal hexahistadine (or 6xHIS) tags and pLATE-51-N with six N-terminal histidine tags were used for the transformation of the Rosetta strain of E. coli, and IMAC purification was performed under native conditions. Peptides synthesis. The multiantigenic peptide 8 (MAP8) format allows the synthesis of peptides with a length of < 15 residues. The synthesis was performed using PepMic (http:// www. pepmic. com/) and provided information including mass spectrometry (MS) and high-performance liquid chromatography (HPLC). We extracted the theoretical and observed molecular weight and purity from the technical sheet. Each peptide was dissolved to a final concentration of 1 mg/mL and used as antigen for antibody detection in serum samples derived from healthy and COVID-19 patients and for the immunization of New Zealand rabbits. Peptides were Sample selection and peptide evaluation using enzyme-linked immunosorbent assay (ELISA). Ten serum samples derived from patients with positive qRT-PCR for SARS-CoV-2 detection with at least fifteen days of infection and ten serum samples from patients with negative diagnoses for SARS-CoV-2 who reported no previous contact with an infected person were used. During all the performed experiments, potentially infectious samples were processed in a biosafety cabinet or inactivated at 56 °C during 1 h. For peptide evaluation, microtiter plates (Sigma Aldrich Cat # CLS3590) were coated with individual recombinant proteins (trimeric spike, recombinant RBD or recombinant N protein), individual peptides (NT1-2, RBD1-5, and N), a mixture of the peptides located in S or all peptides (at a final concentration of 0.1 μg/mL) by adding 100 μL of peptide solution to coating buffer (50 mM Na2CO3/NaCO3H pH 9.6) and incubated overnight at 4 °C. Plates were blocked for 40 min at room temperature (RT) in 200 μL 5% skimmed milk in PBS-Tween (0.05%) and then incubated with human serum (1:50 dilution) for 1 h at 37 °C. After washing, plates were incubated with 100 μL anti-human IgG (whole molecule dilution 1:4000; Sigma-Aldrich Cat# 401455), anti-human IgG (gamma-specific dilution 1:3000; Sigma-Aldrich Cat#A8419) and anti-human IgM (mu-specific dilution 1:4000; Sigma-Aldrich Cat#A0420) coupled to HRP for 1 h at 37 °C. The enzymatic reaction was developed using O-phenylenediamine dihydrochloride (Sigma Aldrich Cat # P1526) and was stopped by the addition of 2 N H 2 SO 4 . The optical density was measured at 492 nm using a microplate reader.

New Zealand rabbit immunization and immune response evaluation. Six New Zealand rabbits
were subcutaneously immunized with 100 µg of each selected peptide in the following order: NT-1 (rabbit #1), RBD3 (rabbit #2), RBD4 (rabbit #3), N (rabbit #4), mixture of all peptides in S protein (rabbit #5) and mixture of all four peptides (rabbit #6) using complete (Life Technologies Cat # P1503) (first immunization) and incomplete Freund's (Sigma Aldrich Cat # F5506) (second and third immunization) as adjuvants. A total of three immunizations were performed with an interval of ten days each. For immune response evaluation, rabbits were bled before immunizations (preimmune serum) and 15 (first evaluation) and 35 (second evaluation) days after the initial immunization. All collected blood samples were evaluated for immune response evaluation using indirect ELISA. Immune response evaluation was performed similarly to the abovementioned ELISA with the following specific differences: antigen fixation was performed using 1 µg/mL of each peptide incubated for 2 h at 37 °C, serial dilutions (1:1000-1:16,000) of the sera were performed and incubated overnight at 4 °C. Finally, an anti-rabbit IgG-HRP coupled (Santa Cruz Biotechnology Cat # sc-2357) was used as a secondary antibody and incubated for 2 h at 37 °C.
Analysis. Obtained data was analyzed using GraphPad Prism V.5.0. For statistical analysis, we performed ANOVA test with Bonferroni adjustment and a p value < 0.05 was considered statistically significant. Ethics declarations. All enrolled patients signed an informed consent form for serum collection considering the Helsinki Declaration, and this study was approved by the Ethics Committee of the Universidad Autónoma de Guerrero. Furthermore, the six rabbits were maintained in individual cages (with individual labels) and the number of animals was calculated considering the number of individual peptides and the potential combinations all procedures. All immunizations and bleeding were performed in the same order. For animal treatment, we consider the ARRIVE guidelines (https:// arriv eguid elines. org). Also, all procedures involving animal manipulation were approved by the ethical committee of the Universidad Autónoma de Guerrero with registry number CICUAL-01/2021, and experiments were performed in accordance with NOM-062-ZOO-1999.

Data availability
All data and supporting data associated with the conception and experiments is available in the supporting information.