Search

Capsicum annuum var. glabriusculum is an economically important horticultural crop and is considered the wild genetic ancestor of chili peppers. The distribution range extends from southern North America, through Central America, to South America. Approximately 226 million 150 paired-end reads were generated from CHMX_Ch1 (a C. annuum from Chihuahua, Mexico). To compare with the CHMX_Ch1 genome, high-quality reads from QO (a C. annuum from Querétaro, Mexico) were downloaded from the NCBI database. A total of 210,324 variants were detected in CHMX_Ch1, whereas 169,718 variants were identified in QO, all compared to the domesticated C. annuum reference genome, UCD10Xv1.1. This comprised 203,990 SNPs and 6334 InDels in CHMX_Ch1 and 164,955 SNPs and 4763 InDels in QO. The variants with high and moderate impact were identified as missense, splice acceptor, splice donor, start lost, stop gain, stop lost, frameshift, insertion, and deletion effects. The candidate genes with the highest fold enrichment values among the SNPs were predominantly involved in gene regulation and metabolic processes. InDels were associated with nuclear and transcriptional regulator activity in both genomes. Overall, a greater number of variants were found in CHMX_Ch1 compared to QO. This study provides knowledge of the principal functions associated with high- and moderate-impact variants and supplies a resource for further investigations of the genetic characteristics of these chiltepin peppers.

Resistance to cisplatin is the main cause of treatment failure in lung adenocarcinoma. Drug-tolerant-persister (DTP) cells are responsible for intrinsic resistance, since they survive the initial cycles of treatment, representing a reservoir for the emergence of clones that display acquired resistance. Although the molecular mechanisms of DTP cells have been described, few studies have investigated the earliest molecular alterations of DTP cells in intrinsic resistance to cisplatin. In this work, we report a gene expression signature associated with the emergence of cisplatin-DTP cells in lung adenocarcinoma cell lines. After a single exposure to cisplatin, we sequenced the transcriptome of cisplatin-DTPs to identify differentially expressed genes. Bioinformatic analysis revealed that early cisplatin-DTP cells deregulate metabolic and proliferative pathways to survive the drug insult. Interaction network analysis identified three highly connected submodules in which SOCS1 had a significant participation in controlling the proliferation of cisplatin-DTP cells. Expression of the candidate genes and their corresponding protein was validated in lung adenocarcinoma cell lines. Importantly, the expression level of SOCS1 was different between CDDP-susceptible and CDDP-resistant lung adenocarcinoma cell lines. Moreover, knockdown of SOCS1 in the CDDP-resistant cell line partially promoted its susceptibility to CDDP. Finally, the clinical relevance of the candidate genes was analyzed in silico, according to the overall survival of cisplatin-treated patients from The Cancer Genome Atlas. Survival analysis showed that downregulation or upregulation of the selected genes was associated with overall survival. The results obtained indicate that these genes could be employed as predictive biomarkers or potential targets to improve the effectiveness of CDDP treatment in lung cancer patients.

Abstract Mexico City has one of the highest incidences of acute lymphoblastic leukemia (ALL) globally, with patients showing low survival, and high relapse rates. To gain more insight into the molecular features of B-ALL in Mexican children, we isolated CD10 + /CD19 + precursor B lymphoblasts from four bone marrow and nine peripheral blood samples of B-ALL patients using a fluorescence-activated cell sorting protocol. The global gene expression profile (BM vs PB) revealed 136 differentially expressed genes; 62 were upregulated (45.6%) and 74 were downregulated (54.4%). Pearson’s correlation coefficient was calculated to determine the similarity between pre-B lymphoblast populations. We selected 26 highly significant genes and validated 21 by RT-qPCR (CNN3, STON2, CALN1, RUNX2, GADD45A, CDC45, CDC20, PLK1, AIDA, HCK, LY86, GPR65, PIK3CG, LILRB2, IL7R, TCL1A, DOCK1, HIST1H3G, PTPN14, CD72, and NT5E). The gene set enrichment analysis of the total expression matrix and the ingenuity pathway analysis of the 136 differentially expressed genes showed that the cell cycle was altered in the bone marrow with four overexpressed genes (PLK1, CDC20, CDC45, and GADD45A) and a low expression of IL7R and PIK3CG, which are involved in B cell differentiation. A comparative bioinformatics analysis of 15 bone marrow and 10 peripheral blood samples from Hispanic B-ALL patients collected by the TARGET program, corroborated the genes observed, except for PIK3CG. We conclude the Mexican and the Hispanic B-ALL patients studied present common driver alterations and histotype-specific mutations that could facilitate risk stratification and diagnostic accuracy and serve as potential therapeutic targets.

Klebsiella pneumoniae is recognized as a common cause of nosocomial infections and outbreaks causing pneumonia, septicemia, and urinary tract infections. This opportunistic bacterium shows an increasing acquisition of antibiotic-resistance genes, which complicates treatment of infections. Hence, fast reliable strain typing methods are paramount for the study of this opportunistic pathogen’s multi-drug resistance genetic profiles. In this study, thirty-eight strains of K. pneumoniae isolated from the blood of pediatric patients were characterized by whole-genome sequencing and genomic clustering methods. Genes encoding β-lactamase were found in all the bacterial isolates, among which the blaSHV variant was the most prevalent (53%). Moreover, genes encoding virulence factors such as fimbriae, capsule, outer membrane proteins, T4SS and siderophores were investigated. Additionally, a multi-locus sequence typing (MLST) analysis revealed 24 distinct sequence types identified within the isolates, among which the most frequently represented were ST76 (16%) and ST70 (11%). Based on LPS structure, serotypes O1 and O3 were the most prevalent, accounting for approximately 63% of all infections. The virulence capsular types K10, K136, and K2 were present in 16, 13, and 8% of the isolates, respectively. Phylogenomic analysis based on virtual genome fingerprints correlated with the MLST data. The phylogenomic reconstruction also denoted association between strains with a higher abundance of virulence genes and virulent serotypes compared to strains that do not possess these traits. This study highlights the value of whole-genomic sequencing in the surveillance of virulence attributes among clinical K. pneumoniae strains.

[This corrects the article DOI: 10.1371/journal.pgen.1006546.].

For the last 500 years, the Americas have been a melting pot both for genetically diverse humans and for the pathogenic and commensal organisms associated with them. One such organism is the stomach-dwelling bacterium Helicobacter pylori, which is highly prevalent in Latin America where it is a major current public health challenge because of its strong association with gastric cancer. By analyzing the genome sequence of H. pylori isolated in North, Central and South America, we found evidence for admixture between H. pylori of European and African origin throughout the Americas, without substantial input from pre-Columbian (hspAmerind) bacteria. In the US, strains of African and European origin have remained genetically distinct, while in Colombia and Nicaragua, bottlenecks and rampant genetic exchange amongst isolates have led to the formation of national gene pools. We found three outer membrane proteins with atypical levels of Asian ancestry in American strains, as well as alleles that were nearly fixed specifically in South American isolates, suggesting a role for the ethnic makeup of hosts in the colonization of incoming strains. Our results show that new H. pylori subpopulations can rapidly arise, spread and adapt during times of demographic flux, and suggest that differences in transmission ecology between high and low prevalence areas may substantially affect the composition of bacterial populations.

Abstract Intrinsically disordered regions (IDRs) are highly enriched in the nucleolar proteome but their physiological role in ribosome assembly remains poorly understood. Our study reveals the functional plasticity of the extremely abundant lysine-rich IDRs of small nucleolar ribonucleoprotein particles (snoRNPs) from protists to mammalian cells. We show in Saccharomyces cerevisiae that the electrostatic properties of this lysine-rich IDR, the KKE/D domain, promote snoRNP accumulation in the vicinity of nascent rRNAs, facilitating their modification. Under stress conditions reducing the rate of ribosome assembly, they are essential for nucleolar compaction and sequestration of key early-acting ribosome biogenesis factors, including RNA polymerase I, owing to their self-interaction capacity in a latent, non-rRNA-associated state. We propose that such functional plasticity of these lysine-rich IDRs may represent an ancestral eukaryotic regulatory mechanism, explaining how nucleolar morphology is continuously adapted to rRNA production levels.

Plant-associated bacteria are an important tool for sustainable agriculture since the have been shown to play a significant role in the health and productivity of plants through various mechanisms. However, commercial inoculants often fail to produce desired effects due to colonization issues and lack of adaptation to local conditions. In this study, we aimed to evaluate the potential of plant growth-promoting bacteria (PGPB) isolated from Solanum hindsianum, a native plant from the Sonoran Desert, to colonize and promote growth in Solanaceae crops under different irrigation conditions. First, endophytic bacteria were screened for various plant growth-promoting traits, and selected strains were identified as Bacillus, Brevibacillus, Pseudomonas, Acinetobacter, and Siccibacter. Then, endophytic and rhizospheric Pseudomonas isolates, were tested for their tolerance to salinity and osmotic stress and for their ability to induce auxin-responsive gene expression and promote plant development in Arabidopsis. All strains induced strong expression of an auxin reporter gene, while significantly improving root and plant morphology, and showed tolerance to osmotic and salinity stress. In colonization experiments, these strains successfully colonized the root crown and rhizosphere of chili, tomatillo, and tomato plants. However, the endophytic strain showed better colonization ability, especially under low irrigation conditions. Overall, this study demonstrates the potential of native PGPR strains to colonize and promote growth in Solanaceae crops under specific environmental conditions.

Abstract Background Rhizosphere microbiomes are fundamental for plant health, development, and productivity, but can be altered by the incidence of soil-borne pathogens. The dysbiosis (disturbance of the microbiome structure of healthy plants) caused by these pathogens, combined with the recruitment of beneficial microorganisms by the diseased plant, may cause shifts in the rhizosphere microbiome during the infection process. These shifts are likely to be associated with changes in the rhizosphere metabolic profile, as the biochemical dialog, or crosstalk, between host plants and their microbiome is mostly mediated by root exudates. Our objective was to elucidate the shifts in the avocado rhizosphere microbiome and associated changes in the rhizosphere metabolome induced by the infection of the oomycete Phytophthora cinnamomi. We also evaluated the effect of inoculating a bacterial biological control agent (BCA) of P. cinnamomi on the avocado rhizosphere microbiome, in the presence and absence of the pathogen, and on morphological and physiological plant variables, to confirm the potential of the BCA to alleviate the stress induced by the disease. Dataset presentation Here, we present a novel dataset collected from a time-course experiment with four treatments: (1) control trees; (2) trees infected with P. cinnamomi; (3) trees inoculated with the BCA; (4) trees infected with P. cinnamomi and inoculated with the BCA. During the infection process, we measured plant morphological and physiological variables and collected rhizosphere soil samples for bacterial and fungal amplicon sequencing, bacterial RNA-seq and metabolomic analyses. Conclusions Collectively, our data elucidate the shifts in the avocado rhizosphere microbiome after infection by P. cinnamomi and when inoculated with a BCA, and help understand how a pathogen or a beneficial bacterium can alter plant-microbiome crosstalk. Understanding the effect of P. cinnamomi or a BCA on the avocado tree physiology and on the avocado rhizosphere microbiome and metabolome will direct our search for disease biomarkers or potential pathogen antagonists, help identify metabolites related to the recruitment of microorganisms, and assist us in developing integrated disease management strategies.

Avocado (Persea americana Mill.) is a tree crop of great social and economic importance. However, the crop productivity is hindered by fast-spreading diseases, which calls for the search of new biocontrol alternatives to mitigate the impact of avocado phytopathogens. Our objectives were to evaluate the antimicrobial activity of diffusible and volatile organic compounds (VOCs) produced by two avocado rhizobacteria (Bacillus A8a and HA) against phytopathogens Fusarium solani, Fusarium kuroshium, and Phytophthora cinnamomi, and assess their plant growth promoting effect in Arabidopsis thaliana. We found that, in vitro, VOCs emitted by both bacterial strains inhibited mycelial growth of the tested pathogens by at least 20%. Identification of bacterial VOCs by gas chromatography coupled to mass spectrometry (GC–MS) showed a predominance of ketones, alcohols and nitrogenous compounds, previously reported for their antimicrobial activity. Bacterial organic extracts obtained with ethyl acetate significantly reduced mycelial growth of F. solani, F. kuroshium, and P. cinnamomi, the highest inhibition being displayed by those from strain A8a (32, 77, and 100% inhibition, respectively). Tentative identifications carried out by liquid chromatography coupled to accurate mass spectrometry of diffusible metabolites in the bacterial extracts, evidenced the presence of some polyketides such as macrolactins and difficidin, hybrid peptides including bacillaene, and non-ribosomal peptides such as bacilysin, which have also been described in Bacillus spp. for antimicrobial activities. The plant growth regulator indole-3-acetic acid was also identified in the bacterial extracts. In vitro assays showed that VOCs from strain HA and diffusible compounds from strain A8a modified root development and increased fresh weight of A. thaliana. These compounds differentially activated several hormonal signaling pathways involved in development and defense responses in A. thaliana, such as auxin, jasmonic acid (JA) and salicylic acid (SA); genetic analyses suggested that developmental stimulation of the root system architecture by strain A8a was mediated by the auxin signaling pathway. Furthermore, both strains were able to enhance plant growth and decreased the symptoms of Fusarium wilt in A. thaliana when soil-inoculated. Collectively, our results evidence the potential of these two rhizobacterial strains and their metabolites as biocontrol agents of avocado pathogens and as biofertilizers.

Reports on the application of the soft Coulomb hole method to theoretical equilibrium geometries of single diatomic molecules. Optimal parameter calculations; Quantum mechanical calculations; Parameter values; Correlation energy.

Hungatella hathewayi has been observed to be a member of the gut microbiome. Unfortunately, little is known about this organism in spite of being associated with human fatalities; it is important to understand virulence mechanisms and epidemiological prospective to cause disease. In this study, a patient with chronic neurologic symptoms presented to the clinic with subsequent isolation of a strain with phenotypic characteristics suggestive of Clostridium difficile. However, whole-genome sequence found the organism to be H. hathewayi. Analysis including publicly available Hungatella genomes found substantial genomic differences as compared to the type strain, indicating this isolate was not C. difficile. We examined the whole-genome of Hungatella species and related genera, using comparative genomics to fully examine species identification and toxin production. Orthogonal phylogenetic using the 16S rRNA gene and entire genome analyses that included genome distance analyses using Genome-to-Genome Distance (GGDC), Average Nucleotide Identity (ANI), and a pan-genome analysis with inclusion of available public genomes determined the speciation to be Hungatella. Two clearly differentiated groups were identified, one including a reference H. hathewayi genome (strain DSM-13,479) and a second group that was determined to be H. effluvii, which included our clinical isolate. Also, some genomes reported as H. hathewayi were found to belong to other genera, including Clostridium and Faecalicatena. We show that the Hungatella species have an open pan-genome reflecting high genomic diversity. This study highlights the importance of correctly assigning taxonomic identification, particularly in disease-associated strains, to better understand virulence and therapeutic options.

Mycobacteria, like other microorganisms, survive under different environmental variations by expressing an efficient adaptive response, oriented by regulatory elements, such as transcriptional repressors of the TetR family. These repressors in mycobacteria also appear to be related to cholesterol metabolism. In this study, we have evaluated the effect of a fatty acid (oleic–palmitic–stearic)/cholesterol mixture on some phenotypic and genotypic characteristics of a tetR-mutant strain (BCG_2177c mutated gene) of M. bovis BCG, a homologous of Rv2160A of M. tuberculosis. In order to accomplish this, we have analyzed the global gene expression of this strain by RNA-seq and evaluated its neutral-lipid storage capacity and potential to infect macrophages. We have also determined the macrophage response by measuring some pro- and anti-inflammatory cytokine expressions. In comparison with wild-type microorganisms, we showed that the mutation in the BCG_2177c gene did not affect the growth of M. bovis BCG in the presence of lipids but it probably modified the structure/composition of its cell envelope. Compared to with dextrose, an overexpression of the transcriptome of the wild-type and mutant strains was observed when these mycobacteria were cultured in lipids, mainly at the exponential phase. Twelve putative intracellular redox balance maintenance genes and four others coding for putative transcriptional factors (including WhiB6 and three TetR-like) were the main elements repeatedly overexpressed when cultured in the presence of lipids. These genes belonged to the central part of what we called the “genetic lipid signature” for M. bovis BCG. We have also found that all these mycobacteria genotypic changes affected the outcome of BCG-infected macrophages, being the mutant strain most adapted to persist longer inside the host. This high persistence result was also confirmed when mutant-infected macrophages showed overexpression of the anti-inflammatory cytokine TGF-β versus pro-inflammatory cytokines. In summary, the lack of this TetR-like repressor expression, within a lipid environment, may help mycobacteria overcome intracellular redox stress and survive longer inside their host.

Typical enteroaggregative Escherichia coli (tEAEC) is a diarrheagenic E. coli pathotype associated with pediatric and traveler’s diarrhea. Even without diarrhea, EAEC infections in children also lead to increased gut inflammation and growth shortfalls. EAEC strain’s defining phenotype is the aggregative adherence pattern on epithelial cells attributable to the aggregative adherence fimbriae (AAF). EAEC only causes diarrhea in humans; therefore, not much is known of the exact intestinal region of infection and damage or its interactions with intestinal enterocytes in vivo and in situ. This study aimed to develop a new tEAEC mouse model of infection, characterize the microbiota of infected mice, and evaluate in situ the expression of host adherence and surface molecules triggering EAEC infection and the role of the EAEC AAF-II in adherence. Six-week-old C57BL/6 mice, without previous antibiotic treatment, were orally challenged with EAEC 042 strain or EAEC 042 AAF-II mutant (ΔAAF/II) strain, or DAEC-MXR strain (diffusely adherent E. coli clinical isolate), and with saline solution (control group). Paraffin sections of the colon and ileum were stained with H&E and periodic acid-Schiff. ZO-1, β-catenin, MUC1, and bacteria were analyzed by immunofluorescence. EAEC-infected mice, in comparison with DAEC-MXR-infected and control mice, significantly lost weight during the first 3 days. After 7 days post-infection, mucus production was increased in the colon and ileum, ZO-1 localization remained unaltered, and morphological alterations were more pronounced in the ileum since increased expression and apical localization of β-catenin in ileal enterocytes were observed. EAEC-infected mice developed dysbiosis 21 days post-infection. At 4 days post-infection, EAEC strain 042 formed a biofilm on ileal villi and increased the expression and apical localization of β-catenin in ileal enterocytes; these effects were not seen in animals infected with the 042 ΔAAF/II strain. At 3 days post-infection, MUC1 expression on ileal enterocytes was mainly detectable among infected mice and colocalized with 042 strains on the enterocyte surface. We developed a novel mouse model of EAEC infection, which mimics human infection, not an illness, revealing that EAEC 042 exerts its pathogenic effects in the mouse ileum and causes dysbiosis. This model is a unique tool to unveil early molecular mechanisms of EAEC infection in vivo and in situ.

The biogenesis of eukaryotic ribosomes involves the ordered assembly of around 80 ribosomal proteins. Supplying equimolar amounts of assembly-competent ribosomal proteins is complicated by their aggregation propensity and the spatial separation of their location of synthesis and pre-ribosome incorporation. Recent evidence has highlighted that dedicated chaperones protect individual, unassembled ribosomal proteins on their path to the pre-ribosomal assembly site. Here, we show that the co-translational recognition of Rpl3 and Rpl4 by their respective dedicated chaperone, Rrb1 or Acl4, reduces the degradation of the encoding RPL3 and RPL4 mRNAs in the yeast Saccharomyces cerevisiae. In both cases, negative regulation of mRNA levels occurs when the availability of the dedicated chaperone is limited and the nascent ribosomal protein is instead accessible to a regulatory machinery consisting of the nascent-polypeptide-associated complex and the Caf130-associated Ccr4-Not complex. Notably, deregulated expression of Rpl3 and Rpl4 leads to their massive aggregation and a perturbation of overall proteostasis in cells lacking the E3 ubiquitin ligase Tom1. Taken together, we have uncovered an unprecedented regulatory mechanism that adjusts the de novo synthesis of Rpl3 and Rpl4 to their actual consumption during ribosome assembly and, thereby, protects cells from the potentially detrimental effects of their surplus production.

Ubiquitous eukaryotic non-coding circular RNAs regulate transcription and translation. We have reported full-length intronic circular RNAs (flicRNAs) in Entamoeba histolytica with esterified 3′ss and 5′ss. Their 5′ss GU-rich elements are essential for their biogenesis and their suggested role in transcription regulation. Here, we explored whether exonic, exonic-intronic, and intergenic circular RNAs are also part of the E. histolytica and E. invadens ncRNA RNAome and investigated their possible functions. Available RNA-Seq libraries were analyzed with the CIRI-full software in search of circular exonic RNAs (circRNAs). The robustness of the analyses was validated using synthetic decoy sequences with bona fide back splice junctions. Differentially expressed (DE) circRNAs, between the virulent HM1:IMSS and the nonvirulent Rahman E. histolytica strains, were identified, and their miRNA sponging potential was analyzed using the intaRNA software. Respectively, 188 and 605 reverse overlapped circRNAs from E. invadens and E. histolytica were identified. The sequence composition of the circRNAs was mostly exonic although different to human circRNAs in other attributes. 416 circRNAs from E. histolytica were virulent-specific and 267 were nonvirulent-specific. Out of the common circRNAs, 32 were DE between strains. Finally, we predicted that 8 of the DE circRNAs could function as sponges of the bioinformatically reported miRNAs in E. histolytica, whose functions are still unknown. Our results extend the E. histolytica RNAome and allow us to devise a hypothesis to test circRNAs/miRNAs/siRNAs interactions in determining the virulent/nonvirulent phenotypes and to explore other regulatory mechanisms during amoebic encystment.

Infections remain a major cause of morbidity and mortality among hematopoietic stem cell transplant (HSCT) recipients. Unlike Epstein–Barr Virus (EBV) and Human Cytomegalovirus (HCMV), Human Herpesvirus (HHV) 6, HHV7 and HHV8 are not routinely monitored in many centers, especially in the pediatric population of low–medium income countries. We screened EBV, HCMV, HHV6, HHV7 and HHV8 in 412 leukocytes-plasma paired samples from 40 pediatric patients assisted in a tertiary hospital in Mexico. Thirty-two underwent allo-HSCT, whereas eight received auto-HSCT. Overall viral detection frequencies in allo- and auto-HSCT were: EBV = 43.7% and 30.0%, HCMV = 5.0% and 6.7%, HHV6 = 7.9% and 20.0% and HHV7 = 9.7% and 23.3%. HHV8 was not detected in any sample. Interestingly, HHV6 and HHV7 were more frequent in auto-HSCT, and HHV6 was observed in all episodes of multiple detection in auto-HSCT patients. We found EBV DNA in plasma samples, whereas HCMV, HHV6 and HHV7 DNA were predominantly observed in leukocytes, indicative of their expansion in cellular compartments. We also found that IL-1β, IL-2, IL-6 and IL-8 were significantly increased in episodes in which multiple viruses were simultaneously detected, and samples positive for EBV DNA and graft-versus-host disease had a further increase of IL-1β and IL-8. In conclusion, the EBV, HCMV, HHV6 and HHV7 burdens were frequently detected in allo- and auto-HSCT, and their presence associated with systemic inflammation.

Galphimia spp. is popularly used in Mexican traditional medicine. Some populations of Galphimia exert anxiolytic and sedative effects due to the presence of the modified triterpenoids galphimines. However, the galphimine synthesis pathway has not yet been elucidated. Hence, in this study, a comparative transcriptome analysis between two contrasting populations of Galphimia spp., a galphimine-producer, and a non-galphimine-producer, is performed using RNA-Seq in the Illumina Next Seq 550 platform to identify putative candidates genes that encode enzymes of this metabolic pathway. Transcriptome functional annotation was performed using the Blast2GO in levels of gene ontology. For differential expression analysis, edgeR, pheatmap, and Genie3 library were used. To validate transcriptome data, qPCR was conducted. In producer and non-producer plants of both populations of Galphimia spp., most of the transcripts were grouped in the Molecular Function level of gene ontology. A total of 680 differentially expressed transcripts between producer and non-producer plants were detected. In galphimine-producer plants, a larger number of highly expressed transcripts related to acyclic and polycyclic terpene synthesis were identified. As putative candidate genes involved in the galphimine synthesis pathway, P450 family members and enzymes with kinase activity were identified.

The structure and function of rhizosphere microbial communities are affected by the plant health status. In this study, we investigated the effect of root rot on the avocado rhizosphere microbiome, using 16S rDNA and ITS sequencing. Furthermore, we isolated potential fungal pathogens associated with root rot symptoms and assessed their pathogenic activity on avocado. We found that root rot did not affect species richness, diversity or community structure, but induced changes in the relative abundance of several microbial taxa. Root rot increased the proportion of Pseudomonadales and Burkholderiales in the rhizosphere but reduced that of Actinobacteria, Bacillus spp. and Rhizobiales. An increase in putative opportunistic fungal pathogens was also detected in the roots of symptomatic trees; the potential pathogenicity of Mortierella sp., Fusarium spp., Lasiodiplodia sp. and Scytalidium sp., is reported for the first time for the State of Veracruz, Mexico. Root rot also potentially modified the predicted functions carried out by rhizobacteria, reducing the proportion of categories linked with the lipid and amino-acid metabolisms whilst promoting those associated with quorum sensing, virulence, and antibiotic resistance. Altogether, our results could help identifying microbial taxa associated to the disease causal agents and direct the selection of plant growth-promoting bacteria for the development of biocontrol microbial consortia.