The Soybean Expression Atlas v2: A comprehensive database of over 5000 RNA-seq samples.

Soybean is a crucial crop worldwide, used as a source of food, feed, and industrial products due to its high protein and oil content. Previously, the rapid accumulation of soybean RNA-seq data in public databases and the computational challenges of processing raw RNA-seq data motivated us to develop the Soybean Expression Atlas, a gene expression database of over a thousand RNA-seq samples. Over the past few years, our database has allowed researchers to explore the expression profiles of important gene families, discover genes associated with agronomic traits, and understand the transcriptional dynamics of cellular processes. Here, we present the Soybean Expression Atlas v2, an updated version of our database with a fourfold increase in the number of samples, featuring transcript- and gene-level transcript abundance matrices for 5481 publicly available RNA-seq samples. New features in our database include the availability of transcript-level abundance estimates and equivalence classes to explore differential transcript usage, abundance estimates in bias-corrected counts to increase the accuracy of differential gene expression analyses, a new web interface with improved data visualization and user experience, and a reproducible and scalable pipeline available as an R package. The Soybean Expression Atlas v2 is available at https://soyatlas.venanciogroup.uenf.br/, and it will accelerate soybean research, empowering researchers with high-quality and easily accessible gene expression data.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms.

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae).

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes, and its only congeneric species, P. strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics, and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae).

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

syntenet: an R/Bioconductor package for the inference and analysis of synteny networks.

SUMMARY: Interpreting and visualizing synteny relationships across several genomes is a challenging task. We previously proposed a network-based approach for better visualization and interpretation of large-scale microsynteny analyses. Here, we present syntenet, an R package to infer and analyze synteny networks from whole-genome protein sequence data. The package offers a simple and complete framework, including data preprocessing, synteny detection and network inference, network clustering and phylogenomic profiling, and microsynteny-based phylogeny inference. Graphical functions are also available to create publication-ready plots. Synteny networks inferred with syntenet can highlight taxon-specific gene clusters that likely contributed to the evolution of important traits, and microsynteny-based phylogenies can help resolve phylogenetic relationships under debate. AVAILABILITY AND IMPLEMENTATION: syntenet is available on Bioconductor (https://bioconductor.org/packages/syntenet), and the source code is available on a GitHub repository (https://github.com/almeidasilvaf/syntenet). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

hybridexpress: an R/Bioconductor package for comparative transcriptomic analyses of hybrids and their progenitors.

Hybridization, the process of crossing individuals from diverse genetic backgrounds, plays a pivotal role in evolution, biological invasiveness, and crop breeding. At the transcriptional level, hybridization often leads to complex nonadditive effects, presenting challenges for understanding its consequences. Although standard transcriptomic analyses exist to compare hybrids to their progenitors, such analyses have not been implemented in a software package, hindering reproducibility. We introduce hybridexpress, an R/Bioconductor package designed to facilitate the analysis, visualization, and comparison of gene expression patterns in hybrid triplets (hybrids and their progenitors). hybridexpress provides users with a user-friendly and comprehensive workflow that includes all standard comparative analyses steps, including data normalization, calculation of midparent expression values, sample clustering, expression-based gene classification into categories and classes, and overrepresentation analysis for functional terms. We illustrate the utility of hybridexpress through comparative transcriptomic analyses of cotton allopolyploidization and rice root trait heterosis. hybridexpress is designed to streamline comparative transcriptomic studies of hybrid triplets, advancing our understanding of evolutionary dynamics in allopolyploids, and enhancing plant breeding strategies. hybridexpress is freely accessible from Bioconductor (https://bioconductor.org/packages/HybridExpress) and its source code is available on GitHub (https://github.com/almeidasilvaf/HybridExpress).

Maternal exposure to glyphosate-based herbicide causes vascular dysfunction in offspring female rats.

This study analyzed how glyphosate exposure in the gestational period affects vascular function in their female offspring and whether oxidative stress is involved in this effect. To this, pregnant Wistar rats were exposed through drinking water to 0.2% of a glyphosate commercial formulation, and we analyzed the response to acetylcholine and phenylephrine in the aorta from offspring of Glyphosate-based herbicide (O-GBH) and controls (O-CON) rats at six months of age. Relaxation to acetylcholine was reduced in O-GBH than in O-CON. Acute Indomethacin and Apocynin increased relaxation to acetylcholine in O-GBH. The aorta from O-GBH was hyperactive to phenylephrine; the preincubation with N-nitro-L-arginine methyl ester (L-NAME) increased contraction to phenylephrine more in O-CON than O-GBH. TEMPOL similarly reduced phenylephrine response, and L-NAME prevented this effect. The TBARS and GSH levels were increased in O-GBH than in O-CON. Results reinforce the concept that oxidative stress during the perinatal period contributes to the development of vascular changes in adulthood. Results also reveal that oxidative stress parameters altered, and the current levels considered safe for exposure to Glyphosate deserve further investigation, especially in the female gender.

Antitussive, Expectorant and Antipyretic Effect of the Ethanolic Extract of the Leaves of Momordica charantia L.

The traditional use of the M. charantia L. plant to treat coughs, fever and expectoration is widely practiced in different cultures, but its effectiveness and safety still require scientific investigation. This study sought to perform a chemical analysis and evaluate the antitussive, expectorant and antipyretic effects of the ethanolic extract of M. charantia leaves (EEMc) in rats and mice. The EEMc was subjected to chemical analysis by HPLC-DAD, revealing the presence of the flavonoids astragalin and isoquercetin. Acute oral toxicity in mice did not result in deaths, although changes in liver weight and stool consistency were observed. EEMc demonstrated an antitussive effect at doses of 100 and 300 mg/kg in mice subjected to cough induction by citric acid nebulization. Furthermore, it showed expectorant activity at a dose of 300 mg/kg, assessed based on the elimination of the phenol red marker in bronchoalveolar lavage. In the evaluation of antipyretic activity in rats, fever induced by Saccharomyces cerevisiae was reduced at all doses tested during the first hour after treatment. This innovative study identified the presence of astragalin and isoquercetin in EEMc and indicated that the extract has antitussive, expectorant and antipyretic properties. Therefore, EEMc presents itself as a promising option in herbal medicine for the treatment of respiratory symptoms and fever.

Unveiling the Dynamics behind Glioblastoma Multiforme Single-Cell Data Heterogeneity.

Glioblastoma Multiforme is a brain tumor distinguished by its aggressiveness. We suggested that this aggressiveness leads single-cell RNA-sequence data (scRNA-seq) to span a representative portion of the cancer attractors domain. This conjecture allowed us to interpret the scRNA-seq heterogeneity as reflecting a representative trajectory within the attractor's domain. We considered factors such as genomic instability to characterize the cancer dynamics through stochastic fixed points. The fixed points were derived from centroids obtained through various clustering methods to verify our method sensitivity. This methodological foundation is based upon sample and time average equivalence, assigning an interpretative value to the data cluster centroids and supporting parameters estimation. We used stochastic simulations to reproduce the dynamics, and our results showed an alignment between experimental and simulated dataset centroids. We also computed the Waddington landscape, which provided a visual framework for validating the centroids and standard deviations as characterizations of cancer attractors. Additionally, we examined the stability and transitions between attractors and revealed a potential interplay between subtypes. These transitions might be related to cancer recurrence and progression, connecting the molecular mechanisms of cancer heterogeneity with statistical properties of gene expression dynamics. Our work advances the modeling of gene expression dynamics and paves the way for personalized therapeutic interventions.

Oral pregabalin is effective as preemptive analgesia in abdominal hysterectomy-A randomized controlled trial.

Postoperative pain is one of the main negative symptoms resulting from surgery and the use of new methods to control this symptom is of ever-increasing relevance. Opioid-sparing strategies, such as multimodal analgesia, are trends in this scenario. Pregabalin is a well-established treatment for neuropathic pain; however, it is still controversial in the surgical context for postoperative analgesia. This study investigated the effect of pregabalin on postoperative analgesia in patients undergoing abdominal hysterectomy. It is a prospective, randomised, double-blind, placebo-controlled clinical trial. Female patients undergoing abdominal hysterectomy were randomised to use pregabalin (group P1), 300 mg orally 2 h before surgery, or identical placebo pills (group P0). The main outcome includes the postoperative pain index by visual analogue scale (VAS) and McGill's pain questionnaire. Secondary outcomes include opioid consumption and the presence of adverse effects. A value of p < 0.05 was used to reject type I error. Fifty-five patients were randomised amongst the groups. Patients in group P1 had lower pain rates by VAS scale, both at rest and in active motion, than group P0. In McGill's questionnaire, patients from group P1 also had lower pain rates (12 × 28.5). There was approximately twice as much opioid consumption amongst patients in group P0. Regarding side effects, there was a difference between the two groups only for dizziness, being more incident in group P1. This study suggests that pregabalin is an important adjuvant drug in treating postoperative pain in patients with abdominal hysterectomy.

A Strategy Utilizing Protein-Protein Interaction Hubs for the Treatment of Cancer Diseases.

We describe a strategy for the development of a rational approach of neoplastic disease therapy based on the demonstration that scale-free networks are susceptible to specific attacks directed against its connective hubs. This strategy involves the (i) selection of up-regulated hubs of connectivity in the tumors interactome, (ii) drug repurposing of these hubs, (iii) RNA silencing of non-druggable hubs, (iv) in vitro hub validation, (v) tumor-on-a-chip, (vi) in vivo validation, and (vii) clinical trial. Hubs are protein targets that are assessed as targets for rational therapy of cancer in the context of personalized oncology. We confirmed the existence of a negative correlation between malignant cell aggressivity and the target number needed for specific drugs or RNA interference (RNAi) to maximize the benefit to the patient's overall survival. Interestingly, we found that some additional proteins not generally targeted by drug treatments might justify the addition of inhibitors designed against them in order to improve therapeutic outcomes. However, many proteins are not druggable, or the available pharmacopeia for these targets is limited, which justifies a therapy based on encapsulated RNAi.

BioNERO: an all-in-one R/Bioconductor package for comprehensive and easy biological network reconstruction.

Currently, standard network analysis workflows rely on many different packages, often requiring users to have a solid statistics and programming background. Here, we present BioNERO, an R package that aims to integrate all aspects of network analysis workflows, including expression data preprocessing, gene coexpression and regulatory network inference, functional analyses, and intraspecies and interspecies network comparisons. The state-of-the-art methods implemented in BioNERO ensure that users can perform all analyses with a single package in a simple pipeline, without needing to learn a myriad of package-specific syntaxes. BioNERO offers a user-friendly framework that can be easily incorporated in systems biology pipelines.

Chemical Analysis by LC-MS of Cannabis sativa Root Samples from Northeast Brazil and Evaluation of Antitussive and Expectorant Activities.

Cannabis sativa is a millenary medicinal plant. However, contrary to worldwide paradigm-shifting, countries like Brazil still prohibit C. sativa cultivation and its medicinal use, even though many populations use aerial parts and roots of this plant for healthcare. As such, the objective of this work was to identify substances in the samples of the C. sativa roots, tracing a correlation with antitussive and expectorant effects. Therefore, samples of C. sativa roots were donated by the Polícia Federal Brasileira, and its aqueous extract (AECsR) was prepared with subsequent lyophilization, to maintain the material stability. After that, the material was analyzed by LC-MS to observe its chemical profile. Four samples (AECsR-A, B, C, and D) were tested in animal models of citric acid-induced cough (0.4 M) and phenol red expectoration (500 mg/kg). Using LC-MS it was possible to identify 5 molecules in C. sativa roots: p-coumaroyltyramine, tetrahydrocannabinol-C4, feruoiltyramine, anhydrocanabisativine, and cannabisativine. In experimental protocols, male mice (Mus musculus) were treated with samples of AECsR at doses of 12.5, 25, or 50 mg/kg regardless of the pharmacological test. In these tests, all samples showed the potential to treat cough and promote fluid expectoration, differing only in the dose at which these effects were observed. Therefore, the data showed that the C. sativa roots of the Brazilian Northeast showed antitussive and expectorant effects, even with intense secondary metabolites' variation, which alters its potency, but not its effect. This highlights the importance of this medicinal plant for future therapy and corroborates to traditional use.

An updated gene regulatory network reconstruction of multidrug-resistant Pseudomonas aeruginosa CCBH4851.

BACKGROUND: Healthcare-associated infections due to multidrug-resistant (MDR) bacteria such as Pseudomonas aeruginosa are significant public health issues worldwide. A system biology approach can help understand bacterial behaviour and provide novel ways to identify potential therapeutic targets and develop new drugs. Gene regulatory networks (GRN) are examples of in silico representation of interaction between regulatory genes and their targets. OBJECTIVES: In this work, we update the MDR P. aeruginosa CCBH4851 GRN reconstruction and analyse and discuss its structural properties. METHODS: We based this study on the gene orthology inference methodology using the reciprocal best hit method. The P. aeruginosa CCBH4851 genome and GRN, published in 2019, and the P. aeruginosa PAO1 GRN, published in 2020, were used for this update reconstruction process. FINDINGS: Our result is a GRN with a greater number of regulatory genes, target genes, and interactions compared to the previous networks, and its structural properties are consistent with the complexity of biological networks and the biological features of P. aeruginosa. MAIN CONCLUSIONS: Here, we present the largest and most complete version of P. aeruginosa GRN published to this date, to the best of our knowledge.

Biscuits Prepared with Enzymatically-Processed Soybean Meal Are Rich in Isoflavone Aglycones, Sensorially Well-Accepted and Stable during Storage for Six Months.

Soybean meal (SBM) is a co-product of the soybean oil industry that is rich in bioactive compounds, such as isoflavones. We aimed to study the effects of processing SBM by fermentation (Saccharomyces cerevisiae) (FSBM) and enzymatic hydrolysis (CelluMax C, a commercial cellulase) (ESBM) on its chemical composition, with emphasis on isoflavones. Fermentation increased protein content by 9%, ash content by 7%, dietary fiber by 11% and minerals by up to 38%, except for iron, which decreased by 26%. Fermentation completely removed oligosaccharides from SBM, while enzymatic processing decreased oligosaccharides by 45% in SBM. Both processes converted glycosylated isoflavones into the corresponding aglycones, the content of which increased by up to 7.7-fold. Biscuits containing SBM, FSBM and ESBM could be labeled as dietary sources of dietary fibers, potassium, phosphorous, calcium and zinc, as well as high in proteins, copper, iron, manganese and magnesium. While FSBM biscuits had lower sensory scores compared to SBM biscuits, ESBM biscuits had equivalent scores. During storage for 180 days at room temperature, the isoflavone profile of all biscuits remained stable. Moreover, storage did not impair microbiological and sensory qualities of any biscuits. Altogether, ESBM biscuits show great marketing potential.

Systematic analysis of 1298 RNA-Seq samples and construction of a comprehensive soybean (Glycine max) expression atlas.

Soybean (Glycine max [L.] Merr.) is a major crop in animal feed and human nutrition, mainly for its rich protein and oil contents. The remarkable rise in soybean transcriptome studies over the past 5 years generated an enormous amount of RNA-seq data, encompassing various tissues, developmental conditions and genotypes. In this study, we have collected data from 1298 publicly available soybean transcriptome samples, processed the raw sequencing reads and mapped them to the soybean reference genome in a systematic fashion. We found that 94% of the annotated genes (52 737/56 044) had detectable expression in at least one sample. Unsupervised clustering revealed three major groups, comprising samples from aerial, underground and seed/seed-related parts. We found 452 genes with uniform and constant expression levels, supporting their roles as housekeeping genes. On the other hand, 1349 genes showed heavily biased expression patterns towards particular tissues. A transcript-level analysis revealed that 95% (70 963 of 74 490) of the assembled transcripts have intron chains exactly matching those from known transcripts, whereas 3256 assembled transcripts represent potentially novel splicing isoforms. The dataset compiled here constitute a new resource for the community, which can be downloaded or accessed through a user-friendly web interface at http://venanciogroup.uenf.br/resources/. This comprehensive transcriptome atlas will likely accelerate research on soybean genetics and genomics.

Analysis of future climate scenarios for northeastern Brazil and implications for human thermal comfort.

A thermal comfort index for the Northeast of Brazil was analyzed for two scenarios of climatic changes, A1B and A2, for 2021-2080, and compared with the reference period 1961-1990. A technique of regionalization was applied to rainfall, maximum and minimum temperature data from meteorological stations, obtained by statistical downscaling of projections from four global climate models. The results pointed to a significant reduction of rainfall and an increase of temperature for three different climatically homogeneous subregions. Regarding the thermal comfort index, the results point to an increase in days with heat discomfort between 2021 and 2080. In the northern portion, the higher percentage of days with heat discomfort will be significant since the first half of the period under appreciation, i.e., from 2021 to 2050. Conversely, in the eastern of northeastern Brazil, the increase of days with heat discomfort should happen in the period from 2051 to 2080, whereas the central-western part of the region, which, in the reference period, had recorded less than 1% of days with heat discomfort, might see an elevation of that percentage to 7% between 2021 and 2050, potentially reaching 48% of its days made uncomfortable by heat between 2051 and 2080.

Exploring the complexity of soybean (Glycine max) transcriptional regulation using global gene co-expression networks.

MAIN CONCLUSION: We report a soybean gene co-expression network built with data from 1284 RNA-Seq experiments, which was used to identify important regulators, modules and to elucidate the fates of gene duplicates. Soybean (Glycine max (L.) Merr.) is one of the most important crops worldwide, constituting a major source of protein and edible oil. Gene co-expression networks (GCN) have been extensively used to study transcriptional regulation and evolution of genes and genomes. Here, we report a soybean GCN using 1284 publicly available RNA-Seq samples from 15 distinct tissues. We found modules that are differentially regulated in specific tissues, comprising processes such as photosynthesis, gluconeogenesis, lignin metabolism, and response to biotic stress. We identified transcription factors among intramodular hubs, which probably integrate different pathways and shape the transcriptional landscape in different conditions. The top hubs for each module tend to encode proteins with critical roles, such as succinate dehydrogenase and RNA polymerase subunits. Importantly, gene essentiality was strongly correlated with degree centrality and essential hubs were enriched in genes involved in nucleic acids metabolism and regulation of cell replication. Using a guilt-by-association approach, we predicted functions for 93 of 106 hubs without functional description in soybean. Most of the duplicated genes had different transcriptional profiles, supporting their functional divergence, although paralogs originating from whole-genome duplications (WGD) are more often preserved in the same module than those from other mechanisms. Together, our results highlight the importance of GCN analysis in unraveling key functional aspects of the soybean genome, in particular those associated with hub genes and WGD events.

Rainfall, not soil temperature, will limit the seed germination of dry forest species with climate change.

Drylands are predicted to become more arid and saline due to increasing global temperature and drought. Although species from the Caatinga, a Brazilian tropical dry forest, are tolerant to these conditions, the capacity for germination to withstand extreme soil temperature and water deficit associated with climate change remains to be quantified. We aimed to evaluate how germination will be affected under future climate change scenarios of limited water and increased temperature. Seeds of three species were germinated at different temperatures and osmotic potentials. Thermal time and hydrotime model parameters were established and thresholds for germination calculated. Germination performance in 2055 was predicted, by combining temperature and osmotic/salt stress thresholds, considering soil temperature and moisture following rainfall events. The most pessimistic climate scenario predicts an increase of 3.9 °C in soil temperature and 30% decrease in rainfall. Under this scenario, soil temperature is never lower than the minimum and seldomly higher than maximum temperature thresholds for germination. As long as the soil moisture (0.139 cm3 cm3) requirements are met, germination can be achieved in 1 day. According to the base water potential and soil characteristics, the minimum weekly rainfall for germination is estimated to be 17.5 mm. Currently, the required minimum rainfall occurs in 14 weeks of the year but will be reduced to 4 weeks by 2055. This may not be sufficient for seedling recruitment of some species in the natural environment. Thus, in future climate scenarios, rainfall rather than temperature will be extremely limiting for seed germination.

Hospitalisations for mycoses as an indicator of socio-environmental vulnerability in the Brazilian Amazon-Savanna transition region.

BACKGROUND: The infections caused by fungi represent a global concern and an important cause of hospital admissions in endemic areas. The influence of socio-environmental factors in infectious diseases has been documented; however, this phenomenon remains unclear regarding mycoses. OBJECTIVES: This study aimed to analyse the spatio-temporal dynamics of hospitalisations for mycoses (HM) and the association with socio-economic and climate data in the Amazon-Savanna Transition Region in the state of Maranhão, Brazil. METHODS: In this study, Spearman's correlation was applied to determine the correlation between HM, socio-economic and climatic data obtained from national databases in the period from 1998 to 2016. Hospitalisations for mycoses data were spatialised and analysed using the local Moran's index. RESULTS: Our data revealed a negative and significant correlation between HM and socio-economic data regarding population, demographic density, human development index, health facilities and sanitary sewage. Significant correlations were observed between HM and precipitation, maximum temperature and minimum temperature. The main modulating climatic variable was the minimum temperature. The spatial autocorrelation analysis showed the dynamics of HM in municipalities belonging to the different regions of the state influenced by socio-economic conditions. We observed the presence of municipalities with high incidence of HM surrounded by others with low HM cases and vice versa. CONCLUSIONS: Our results indicate that hospitalisations for mycoses represent an important indicator of socio-environmental vulnerability in the Amazon-Savanna transition region in Brazil. We encourage the adoption of measures to mitigate social and environmental impact on these diseases, especially in municipalities with low socio-economic status.