Genetic insights into the microevolutionary dynamics and early introductions of human monkeypox virus in Mexico.

The recent global outbreak of mpox, caused by monkeypox virus (MPV) emerged in Europe in 2022 and rapidly spread to over 40 countries. The Americas are currently facing the highest impact, reporting over 50,000 cases by early 2023. In this study, we analyzed 880 MPV isolates worldwide to gain insights into the evolutionary patterns and initial introduction events of the virus in Mexico. We found that MPV entered Mexico on multiple occasions, from the United Kingdom, Portugal, and Canada, and subsequently spread locally in different regions of Mexico. Additionally, we show that MPV has an open pangenome, highlighting the role of gene turnover in shaping its genomic diversity, rather than single-nucleotide polymorphisms (SNPs), which do not contribute significantly to genome diversity. Although the genome contains multiple SNPs in coding regions, these remain under purifying selection, suggesting their evolutionary conservation. One notable exception is amino acid position 63 of the protein encoded by the Cop-A4L gene, which is intricately related to viral maturity, which we found to be under strong positive selection. Ancestral state reconstruction indicated that the ancestral state at position 63 corresponds to the amino acid valine, which is present only in isolates of clade I. However, the isolates from the current outbreak contained threonine at position 63. Our findings contribute new information about the evolution of monkeypox virus.

Genome Evolution and Early Introductions of the SARS-CoV-2 Omicron Variant in Mexico.

A new variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), named Omicron (Pango lineage designation B.1.1.529), was first reported to the World Health Organization by South African health authorities on 24 November 2021. The Omicron variant possesses numerous mutations associated with increased transmissibility and immune escape properties. In November 2021, Mexican authorities reported Omicron's presence in the country. In this study, we infer the first introductory events of Omicron and the impact that human mobility has had on the spread of the virus. We also evaluated the adaptive evolutionary processes in Mexican SARS-CoV-2 genomes during the first month of the circulation of Omicron. We inferred 160 introduction events of Omicron in Mexico since its first detection in South Africa; subsequently, after the first introductions there was an evident increase in the prevalence of SARS-CoV-2 during January. This higher prevalence of the novel variant resulted in a peak of reported cases; on average 6 weeks after, a higher mobility trend was reported. During the peak of cases in the country from January to February 2022, the Omicron BA.1.1 sub-lineage dominated, followed by the BA.1 and BA.15 sub-lineages. Additionally, we identified the presence of diversifying natural selection in the genomes of Omicron and found six non-synonymous mutations in the receptor binding domain of the spike protein, all of them related to evasion of the immune response. In contrast, the other proteins in the genome are highly conserved; however, we identified homoplasic mutations in non-structural proteins, indicating a parallel evolution.

The Alpha Variant (B.1.1.7) of SARS-CoV-2 Failed to Become Dominant in Mexico.

During the coronavirus disease 2019 (COVID-19) pandemic, the emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in the United Kingdom in September 2020, was well documented in different areas of the world and became a global public health concern because of its increased transmissibility. The B.1.1.7 lineage was first detected in Mexico during December 2020, showing a slow progressive increase in its circulation frequency, which reached its maximum in May 2021 but never became predominant. In this work, we analyzed the patterns of diversity and distribution of this lineage in Mexico using phylogenetic and haplotype network analyses. Despite the reported increase in transmissibility of the B.1.1.7 lineage, in most Mexican states, it did not displace cocirculating lineages, such as B.1.1.519, which dominated the country from February to May 2021. Our results show that the states with the highest prevalence of B.1.1.7 were those at the Mexico-U.S. border. An apparent pattern of dispersion of this lineage from the northern states of Mexico toward the center or the southeast was observed in the largest transmission chains, indicating possible independent introduction events from the United States. However, other entry points cannot be excluded, as shown by multiple introduction events. Local transmission led to a few successful haplotypes with a localized distribution and specific mutations indicating sustained community transmission. IMPORTANCE The emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) throughout the world were due to its increased transmissibility. However, it did not displace cocirculating lineages in most of Mexico, particularly B.1.1.519, which dominated the country from February to May 2021. In this work, we analyzed the distribution of B.1.1.7 in Mexico using phylogenetic and haplotype network analyses. Our results show that the states with the highest prevalence of B.1.1.7 (around 30%) were those at the Mexico-U.S. border, which also exhibited the highest lineage diversity, indicating possible introduction events from the United States. Also, several haplotypes were identified with a localized distribution and specific mutations, indicating that sustained community transmission occurred in the country.

Evolutionary and functional relationships in the ribosome biogenesis SBDS and EFL1 protein families.

Nascent ribosomal 60S subunits undergo the last maturation steps in the cytoplasm. The last one involves removing the anti-association factor eIF6 from the 60S ribosomal surface by the joint action of the Elongation Factor-like 1 (EFL1) GTPase and the SBDS protein. Herein, we studied the evolutionary relationship of the EFL1 and EF-2 protein families and the functional conservation within EFL1 orthologues. Phylogenetic analysis demonstrated that the EFL1 proteins are exclusive of eukaryotes and share an evolutionary origin with the EF-2 and EF-G protein families. EFL1 proteins originated by gene duplication from the EF-2 proteins and specialized in ribosome maturation while the latter retained their function in translation. Some organisms have more than one EFL1 protein resulting from alternative splicing, while others are encoded in different genes originated by gene duplication. However, the function of these alternative EFL1 proteins is still unknown. We performed GTPase activity and complementation assays to study the functional conservation of EFL1 homologs alone and together with their SBDS counterparts. None of the orthologues or cross-species combinations could replace the function of the corresponding yeast EFL1•SBDS binomial. The complementation of SBDS interspecies chimeras indicates that domain 2 is vital for its function together with EFL1 and the 60S subunit. The results suggest a functional species-specificity and possible co-evolution between EFL1, SBDS, and the 60S ribosomal subunit. These findings set the basis for further studies directed to understand the molecular evolution of these proteins and their impact on ribosome biogenesis and disease.

Identification of the allosteric site for neutral amino acids in the maize C4 isozyme of phosphoenolpyruvate carboxylase: The critical role of Ser-100.

The isozymes of photosynthetic phosphoenolpyruvate carboxylase from C4 plants (PEPC-C4) play a critical role in their atmospheric CO2 assimilation and productivity. They are allosterically activated by phosphorylated trioses or hexoses, such as d-glucose 6-phosphate, and inhibited by l-malate or l-aspartate. Additionally, PEPC-C4 isozymes from grasses are activated by glycine, serine, or alanine, but the allosteric site for these compounds remains unknown. Here, we report a new crystal structure of the isozyme from Zea mays (ZmPEPC-C4) with glycine bound at the monomer-monomer interfaces of the two dimers of the tetramer, making interactions with residues of both monomers. This binding site is close to, but different from, the one proposed to bind glucose 6-phosphate. Docking experiments indicated that d/l-serine or d/l-alanine could also bind to this site, which does not exist in the PEPC-C4 isozyme from the eudicot plant Flaveria, mainly because of a lysyl residue at the equivalent position of Ser-100 in ZmPEPC-C4 Accordingly, the ZmPEPC-C4 S100K mutant is not activated by glycine, serine, or alanine. Amino acid sequence alignments showed that PEPC-C4 isozymes from the monocot family Poaceae have either serine or glycine at this position, whereas those from Cyperaceae and eudicot families have lysine. The size and charge of the residue equivalent to Ser-100 are not only crucial for the activation of PEPC-C4 isozymes by neutral amino acids but also affect their affinity for the substrate phosphoenolpyruvate and their allosteric regulation by glucose 6-phosphate and malate, accounting for the reported kinetic differences between PEPC-C4 isozymes from monocot and eudicot plants.

A novel motif in the NaTrxh N-terminus promotes its secretion, whereas the C-terminus participates in its interaction with S-RNase in vitro.

BACKGROUND: NaTrxh, a thioredoxin type h, shows differential expression between self-incompatible and self-compatible Nicotiana species. NaTrxh interacts in vitro with S-RNase and co-localizes with it in the extracellular matrix of the stylar transmitting tissue. NaTrxh contains N- and C-terminal extensions, a feature shared by thioredoxin h proteins of subgroup 2. To ascertain the function of these extensions in NaTrxh secretion and protein-protein interaction, we performed a deletion analysis on NaTrxh and fused the resulting variants to GFP. RESULTS: We found an internal domain in the N-terminal extension, called Nß, that is essential for NaTrxh secretion but is not hydrophobic, a canonical feature of a signal peptide. The lack of hydrophobicity as well as the location of the secretion signal within the NaTrxh primary structure, suggest an unorthodox secretion route for NaTrxh. Notably, we found that the fusion protein NaTrxh-GFP(KDEL) is retained in the endoplasmic reticulum and that treatment of NaTrxh-GFP-expressing cells with Brefeldin A leads to its retention in the Golgi, which indicates that NaTrxh uses, to some extent, the endoplasmic reticulum and Golgi apparatus for secretion. Furthermore, we found that Nß contributes to NaTrxh tertiary structure stabilization and that the C-terminus functions in the protein-protein interaction with S-RNase. CONCLUSIONS: The extensions contained in NaTrxh sequence have specific functions on the protein. While the C-terminus directly participates in protein-protein interaction, particularly on its interaction with S-RNase in vitro; the N-terminal extension contains two structurally different motifs: Nα and Nß. Nß, the inner domain (Ala-17 to Pro-27), is essential and enough to target NaTrxh towards the apoplast. Interestingly, when it was fused to GFP, this protein was also found in the cell wall of the onion cells. Although the biochemical features of the N-terminus suggested a non-classical secretion pathway, our results provided evidence that NaTrxh at least uses the endoplasmic reticulum, Golgi apparatus and also vesicles for secretion. Therefore, the Nß domain sequence is suggested to be a novel signal peptide.

Metabolic analysis of Chlorobium chlorochromatii CaD3 reveals clues of the symbiosis in 'Chlorochromatium aggregatum'.

A symbiotic association occurs in 'Chlorochromatium aggregatum', a phototrophic consortium integrated by two species of phylogenetically distant bacteria composed by the green-sulfur Chlorobium chlorochromatii CaD3 epibiont that surrounds a central ß-proteobacterium. The non-motile chlorobia can perform nitrogen and carbon fixation, using sulfide as electron donors for anoxygenic photosynthesis. The consortium can move due to the flagella present in the central ß-protobacterium. Although Chl. chlorochromatii CaD3 is never found as free-living bacteria in nature, previous transcriptomic and proteomic studies have revealed that there are differential transcription patterns between the symbiotic and free-living status of Chl. chlorocromatii CaD3 when grown in laboratory conditions. The differences occur mainly in genes encoding the enzymatic reactions involved in nitrogen and amino acid metabolism. We performed a metabolic reconstruction of Chl. chlorochromatii CaD3 and an in silico analysis of its amino acid metabolism using an elementary flux modes approach (EFM). Our study suggests that in symbiosis, Chl. chlorochromatii CaD3 is under limited nitrogen conditions where the GS/GOGAT (glutamine synthetase/glutamate synthetase) pathway is actively assimilating ammonia obtained via N2 fixation. In contrast, when free-living, Chl. chlorochromatii CaD3 is in a condition of nitrogen excess and ammonia is assimilated by the alanine dehydrogenase (AlaDH) pathway. We postulate that 'Chlorochromatium aggregatum' originated from a parasitic interaction where the N2 fixation capacity of the chlorobia would be enhanced by injection of 2-oxoglutarate from the ß-proteobacterium via the periplasm. This consortium would have the advantage of motility, which is fundamental to a phototrophic bacterium, and the syntrophy of nitrogen and carbon sources.

Amino acid residues critical for the specificity for betaine aldehyde of the plant ALDH10 isoenzyme involved in the synthesis of glycine betaine.

Plant Aldehyde Dehydrogenase10 (ALDH10) enzymes catalyze the oxidation of ω-primary or ω-quaternary aminoaldehydes, but, intriguingly, only some of them, such as the spinach (Spinacia oleracea) betaine aldehyde dehydrogenase (SoBADH), efficiently oxidize betaine aldehyde (BAL) forming the osmoprotectant glycine betaine (GB), which confers tolerance to osmotic stress. The crystal structure of SoBADH reported here shows tyrosine (Tyr)-160, tryptophan (Trp)-167, Trp-285, and Trp-456 in an arrangement suitable for cation-π interactions with the trimethylammonium group of BAL. Mutation of these residues to alanine (Ala) resulted in significant K(m)(BAL) increases and V(max)/K(m)(BAL) decreases, particularly in the Y160A mutant. Tyr-160 and Trp-456, strictly conserved in plant ALDH10s, form a pocket where the bulky trimethylammonium group binds. This space is reduced in ALDH10s with low BADH activity, because an isoleucine (Ile) pushes the Trp against the Tyr. Those with high BADH activity instead have Ala (Ala-441 in SoBADH) or cysteine, which allow enough room for binding of BAL. Accordingly, the mutation A441I decreased the V(max)/K(m)(BAL) of SoBADH approximately 200 times, while the mutation A441C had no effect. The kinetics with other ω-aminoaldehydes were not affected in the A441I or A441C mutant, demonstrating that the existence of an Ile in the second sphere of interaction of the aldehyde is critical for discriminating against BAL in some plant ALDH10s. A survey of the known sequences indicates that plants have two ALDH10 isoenzymes: those known to be GB accumulators have a high-BAL-affinity isoenzyme with Ala or cysteine in this critical position, while non GB accumulators have low-BAL-affinity isoenzymes containing Ile. Therefore, BADH activity appears to restrict GB synthesis in non-GB-accumulator plants.

The family of maize D-type cyclins: genomic organization, phylogeny and expression patterns.

Cyclin proteins, associated to cyclin-dependent kinases (CDKs), play fundamental roles in cell cycle control as they constitute a very important driving force to allow cell cycle progression. D-type cyclins (CycDs) are important both for interpreting external mitogenic signals and in the control of the G1 phase. The maize (Zea mays) genome appears to contain at least 17 different CycD genes, and they fall into the subgroups previously described for other plants. Maize CycDs have been named according to identity percentages of the corresponding orthologs in rice and Arabidopsis. In silico analysis confirmed the presence of characteristic cyclin domains in each maize CycD gene and showed that their genomic organization is similar to their orthologs in rice and Arabidopsis. The expression of maize CycD genes was followed in seeds, during germination in the presence/absence of exogenously added hormones, and also in different plantlet tissues (mesocotyl, root tips and first leaf). Most cyclins were expressed in germinating seeds and at least in one of the plantlet tissues tested; almost all of the detected cyclins show an accumulating pattern of mRNA along germination (0-24 h) and higher levels in root tissue. Interestingly, some cyclins show high levels in non-proliferating tissues as leaf. Addition of auxins or cytokinins does not seem to importantly modify transcript levels; on the other hand, addition of abscisic acid repressed the expression of several cyclins. The role of each CycD during germination and plant growth and its interaction with other cell cycle proteins becomes a topic of the highest interest.

A score of the ability of a three-dimensional protein model to retrieve its own sequence as a quantitative measure of its quality and appropriateness.

BACKGROUND: Despite the remarkable progress of bioinformatics, how the primary structure of a protein leads to a three-dimensional fold, and in turn determines its function remains an elusive question. Alignments of sequences with known function can be used to identify proteins with the same or similar function with high success. However, identification of function-related and structure-related amino acid positions is only possible after a detailed study of every protein. Folding pattern diversity seems to be much narrower than sequence diversity, and the amino acid sequences of natural proteins have evolved under a selective pressure comprising structural and functional requirements acting in parallel. PRINCIPAL FINDINGS: The approach described in this work begins by generating a large number of amino acid sequences using ROSETTA [Dantas G et al. (2003) J Mol Biol 332:449-460], a program with notable robustness in the assignment of amino acids to a known three-dimensional structure. The resulting sequence-sets showed no conservation of amino acids at active sites, or protein-protein interfaces. Hidden Markov models built from the resulting sequence sets were used to search sequence databases. Surprisingly, the models retrieved from the database sequences belonged to proteins with the same or a very similar function. Given an appropriate cutoff, the rate of false positives was zero. According to our results, this protocol, here referred to as Rd.HMM, detects fine structural details on the folding patterns, that seem to be tightly linked to the fitness of a structural framework for a specific biological function. CONCLUSION: Because the sequence of the native protein used to create the Rd.HMM model was always amongst the top hits, the procedure is a reliable tool to score, very accurately, the quality and appropriateness of computer-modeled 3D-structures, without the need for spectroscopy data. However, Rd.HMM is very sensitive to the conformational features of the models' backbone.