- 1Centro de Biología Molecular y Farmacogenética, Núcleo Científico-Tecnológico en Biorecursos BIOREN, Universidad de La Frontera, Temuco, Chile
- 2Tecnología Médica, Facultad de Salud, Universidad Santo Tomás, Temuco, Chile
- 3Centro de Genética y Genómica, Facultad de Medicina, Universidad del Desarrollo, Santiago, Chile
- 4Subdepartamento de Genética Molecular, Instituto de Salud Pública de Chile, Santiago, Chile
- 5Programa Institucional de Fomento a la Investigación, Desarrollo e Innovación, Universidad Tecnológica Metropolitana, Santiago, Chile
- 6Centro de Biotecnología de los Recursos Naturales, Facultad de Ciencias Agrarias y Forestales, Universidad Católica del Maule, Talca, Chile
- 7Advanced Center for Chronic Diseases - ACCDiS, Facultad de Química y Ciencias Farmacéuticas, Universidad de Chile, Santiago, Chile
- 8Departamento de Ciencias Agronómicas y Recursos Naturales, Facultad de Ciencias Agropecuarias y Medioambiente, Temuco, Chile
- 9Departamento de Medicina Interna, Facultad de Medicina, Universidad de La Frontera, Temuco, Chile
Introduction: In recent years, several studies have evidenced the importance of the microbiome to host physiology as metabolism regulator, along with its potential role in triggering various diseases. In this study, we analyzed the gut microbiota in hypercholesterolemic (cases) and normocholesterolemic (controls) individuals to identify characteristic microbial signature for each condition.
Methods: Stool samples were obtained from 57 adult volunteers (27 hypercholesterolemic and 30 controls). The taxonomic profiling of microbial communities was performed using high-throughput sequencing of 16S rRNA V3-V4 amplicons, followed by data analysis using Quantitative Insights Into Microbial Ecology 2 (QIIME2) and linear discriminant analysis (LDA) effect size (LEfSe).
Results: Significant differences were observed in weight, height, body mass index (BMI) and serum levels of triglycerides, total cholesterol and low-density lipoprotein cholesterol (LDL-C) between the groups (p<0.05). LEfSe showed differentially abundant prokaryotic taxa (α=0.05, LDA score > 2.0) in the group of hypercholesterolemic individuals (Methanosphaera, Rothia, Chromatiales, Clostridiales, Bacillaceae and Coriobacteriaceae) and controls (Faecalibacterium, Victivallis and Selenomonas) at various taxonomic levels. In addition, through the application of Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2), the predominance of pathways related to biosynthesis in hypercholesterolemic patients was established, compared to controls in which degradation pathways were predominant. Finally, in the analysis of co-occurrence networks, it was possible to identify associations between the microorganisms present in both studied groups.
Conclusion: Our results point out to unique microbial signatures, which likely play a role on the cholesterol metabolism in the studied population.
Introduction
The human microbiota is composed of cells from all three domains of life, being the prokaryotic assemblages the most well-known. It is estimated that humans harbor approximately 10 bacterial cells per eukaryotic cell, most of which are located in the gastrointestinal tract, where produce vitamins and process dietary components that are not digestible by another way (Bäckhed et al., 2005) . However, in recent years new calculations have been estimated the ratio between bacterial and human cells in 1:1, although the relevance of the microbiota and its impact on the body would not be altered (Sender et al., 2016). The imbalance of bacterial communities, known as dysbiosis, has been associated with the development of various diseases in recent years. Some examples are the association of the microbiome composition with cardiovascular diseases (CVD), diabetes mellitus, obesity, metabolic syndrome, and hypertension (Wu et al., 2015; Mazidi et al., 2016). Furthermore, an important association has been established between diet, gut microbiota composition, and CVDs (Mendelsohn and Larrick, 2013). Moreover, it has been suggested that the characterization of the gut microbiota could be a useful tool in the description of cardiovascular health status in the context of precision medicine (Antman and Loscalzo, 2016). In addition, other studies have associated the presence of gut microorganisms with trimethylamine, which originates from dietary carnitine and lecithin, contributing to development of CVDs (Mendelsohn and Larrick, 2013). Similarly, a possible relation of the intestinal microbiota with the metabolism of circulating lipids in blood should be considered, which could contribute to the development of CVDs, and which may be of interest for the treatment of these diseases (Mazidi et al., 2016) Currently, some studies show the association between microbial groups of the intestinal microbiota with classical lipid parameters and apolipoproteins A1 and B (Yun et al., 2020).
On the other hand, hypercholesterolemia is considered an important modifiable risk factor of CVD, since high levels of low-density lipoprotein cholesterol (LDL-C) are known to increase the risk of coronary artery disease (CAD). In addition, the pharmacological reduction of LDL-C also reduces the incidence and mortality by coronary events (Baigent et al., 2005; Yusuf et al., 2006). Studies in animal models have shown a relationship between the gut microbiota and disturbances in lipid metabolism. Mice lacking their gut microbiota present differences in lipid absorption resulting in variations of serum lipid levels (Wostmann, 1973). Other models have also confirmed the importance of bacteria in cholesterol metabolism, highlighting the relevance of microorganisms in this biological process and their therapeutic potential (Bourgin et al., 2020). Furthermore, specific bacterial groups detected in stool samples have been shown to have beneficial effects on lipid metabolism induced by dietary changes (Martinez et al., 2009). A previous study by our group, using the polymerase chain reaction followed by denaturing gradient gel electrophoresis (PCR-DGGE), showed that subjects with hypercholesterolemia and controls exhibit different gut microbiota profiles (Rebolledo et al., 2017).
Thus, the objective of this study was to describe the microbial communities differentially present in the gut microbiota of hypercholesterolemic and control individuals.
Materials and methods
Subjects
A case-control study was designed, including 27 hypercholesterolemic (cases) and 30 normocholesterolemic individuals (controls). LDL-C levels were considered to define these groups by using the criteria of the National Cholesterol Education Program (NCEP) (Grundy et al., 2004). Participants were recruited by the Center for Cardiological and Internal Medicine Studies of the Universidad de La Frontera (Temuco, Chile). Thus, individuals with concentrations of LDL-C greater than 4.16 mmol/L and less than 2.60 mmol/L were classified as hypercholesterolemic or control, respectively. The presence of inflammatory bowel disease, diabetes mellitus, class 2 and morbid obesity, use of lipid-lowering medications and treatment with antibiotics in the last six months were established as exclusion criteria. For this study, the participants were volunteers and signed an informed consent that was approved by the Scientific Ethics Committee of the Universidad de La Frontera, Temuco, Chile (protocol number 113_2016).
Samples
Blood samples were collected by using standard venipuncture technique and serum was obtained to perform biochemical analysis. Additionally, stool samples were obtained from participants using a DNA/RNA Shield-Fecal Collection Tube (Zymo Research, USA). Then, samples were transported to the laboratory and stored at -80°C until DNA isolation.
Anthropometric, clinical and biochemical parameters
Anthropometric and clinical parameters such as age, height, weight, body mass index (BMI), systolic blood pressure (SBP) and diastolic blood pressure (DBP) were recorded. Serum glucose, triglycerides, total cholesterol, high-density lipoprotein cholesterol (HDL-C) and LDL-C were quantified using enzymatic colorimetric methods. Quality measurements were controlled by using normal and pathological commercial serums (Wiener Lab., Argentina).
DNA extraction and 16S rRNA gene sequencing
DNA was isolated from stool samples using the Fast DNA Spin Kit (MP Biomedicals, USA), following the manufacturer’s instructions. The V3 and V4 region of prokariotic 16S rRNA genes were amplified by polymerase chain reaction (PCR), using 16SF (5’-CCTACGGGNGGCWGCAG-3’) and 16SR (5’-GACTACHVGGGTATCTAATCC-3’) primers (Klindworth et al., 2013). PCR reactions were carried out in 50 µL, with 25 µL 2X KAPA HiFi HotStart ReadyMix (Roche, Switzerland), 50 ng of DNA and primers to final concentration of 200 nM. Thermal cycling consisted of initial denaturation at 95°C for 3 minutes, followed by 25 cycles of denaturation at 95°C for 30 seconds, annealing at 55°C for 30 seconds, and extension 72°C for 30 seconds. A final extension of 72°C for 5 minutes was applied. The products were analyzed by capillary electrophoresis in Fragment Analyzer™ (Agilent Technologies, USA). Then, libraries preparation was obtained following manufacturer’s instructions, quantified in a Quantus Fluorometer (Promega, USA) and analyzed by capillary electrophoresis in Fragment Analyzer™ (Agilent Technologies, USA). Finally, sequencing was performed on a MiSeq sequencer (Illumina, USA), after preparation and loading on a MiSeq Reagent Kit V3 600 cycles following the manufacturer’s instructions.
Community structure analysis
Sequence analysis
The quality control of the sequences was carried out using Divisive Amplicon Denoising Algorithm 2 (DADA2) and reads were demultiplexed and processed using Quantitative Insights Into Microbial Ecology 2 (QIIME2) (Caporaso et al., 2010), using amplicon sequence variants (ASVs) for the analysis. For alpha diversity, Faith’s phylogenetic diversity (PD) index was evaluated. PD and the alpha rarefaction curves were obtained through a QIIME2 plugin. For beta diversity, unweighted UniFrac distance was evaluated, applying PERMANOVA with 999 permutations (alpha = 5%) in addition to the Jaccard distance. Principal Coordinates Analysis (PCoA) 3D graphics were generated using EMPeror (Vázquez-Baeza et al., 2013). The Kolmogorov-Smirnov test was applied to determine the statistical significance (alpha = 5%).
Linear discriminant analysis effect size
LEfSe (Segata et al., 2011) was applied to a table of Features (obtained from QIIME2) to identify potential microorganisms groups that could be associated with hypercholesterolemic individuals and controls. Differentially abundant bacterial groups were defined with a LDA score (log10) over 2. The Kruskal-Wallis test was applied as a factorial test for the classes. The level of significance was alpha = 5%.
Phylogenetic investigation of communities by reconstruction of unobserved states 2
PICRUSt2 allows for the prediction of the functional abundance of bacterial groups based on 16S rRNA sequencing data. For this analysis, a file of aligned Amplicon Sequence Variants (ASVs) was generated using the QIIME data to obtain gene predictions. Finally, gene family profiles and pathway abundances for samples and groups were obtained (Douglas et al., 2020). Statistical Analysis of Metagenomic Profiles (STAMP) tool was used to visualize the results (Parks et al., 2014).
Co-occurrence networks analysis
The co-occurrence analysis was based on the methodology proposed by Bradley (Bradley, 2020), which considers the probabilistic co-occurrence model suggested by Veech (2013), establishing the co-occurrence by pairs and delivering the significant interactions (alpha = 5%) (Veech, 2013). For the latter, the coexistence of pairs of species with lower and higher frequencies than what would be expected to be found by chance is taken into account. Finally, the visualization of the results was carried out using visNetwork, based on R, which takes into account the data frames that describe the nodes and edges of the networks.
Statistical analysis
The data was analyzed using the software GraphPad Prisma v. 8.0 (GraphPad Software Inc., USA). Descriptive statistics were obtained for each variable on the data of all participants. In addition, the biochemical and anthropometric characteristics of the groups were compared using the Mann-Whitney test. The detected phyla of both groups were compared with the Kolmogorov-Smirnov test. Statistical significance was established with a p value less than 0.05.
Results
Anthropometric, clinical and biochemical parameters
Anthropometric, clinical and biochemical parameters are summarized in Table 1. The mean age between the groups was not significantly different (p=0.525). Similarly, there were no significant differences in SBP, DBP, glucose and HDL-C levels. As expected, triglycerides (p=0.003), total cholesterol (p<0.0001) and LDL-C (p<0.0001) serum levels were significantly higher in hypercholesterolemic individuals. However, weight (p<0.0001) and BMI (p=0.004) were significantly higher in control individuals.
Table 1 Anthropometric and biochemical characteristics of hypercholesterolemic individuals and controls.
Community structure analysis
Sequence analysis
The average number of sequences obtained in the MiSeq platform was 318,242 reads per individual. The total number of reads obtained was 18,776,291. There were no significant differences between hypercholesterolemic individuals and controls respect to the evaluation of alpha (Faith’s PD index, p=0.396) and beta (unweighted UniFrac distance, p=0.133) diversity. Jaccard distance, indicated in the Figure 1, shows that for most of the samples analyzed there is no dissimilarity in the microbial communities, with the exception of some cases of hypercholesterolemic participants (red dots), at the bottom of the graph.
Figure 1 Qualitative measurement of the dissimilarity (Jaccard distance) of the microbial communities associated with of hypercholesterolemic individuals and controls. In the Principal Coordinates Analysis (PCoA), produced by EMPeror, hypercholesterolemic individuals are indicated in red, while those control individuals are indicated in blue.
As results of the taxonomic analysis, Figure 2 represents the relative abundance of the different microbial communities of hypercholesterolemic individuals and controls. No significant differences were observed in most of the phylum, excepting Actinobacteria (p=0.025), predominantly present in hypercholesterolemic individuals. The relative abundance of Bacteroidetes was 50.63% in hypercholesterolemic individuals and 47.72% in controls (p=0.845). Firmicutes accounted for 35.57% and 43.66% of the microbial taxa in hypercholesterolemic individuals and controls (p=0.279), respectively. The complete microbial profile is presented in Table 2.
Figure 2 Relative abundance obtained by taxonomic analysis with Quantitative Insights Into Microbial Ecology 2 (QIIME2), which indicate the different phylum present in the deposition samples of hypercholesterolemic individuals and controls.
Table 2 Bacterial phylum detected in fecal samples of hypercholesterolemic individuals and controls.
Differentially abundant bacterial groups present in hypercholesterolemic and controls
The LEfSe comparison between the stool microbiota of hypercholesterolemic individuals and controls suggests the presence of some specific taxa to each group. The genera Faecalibacterium, Victivallis and Selenomonas were likely associated with controls subjects, while the microbiota of participants with hypercholesterolemia contains unique archaeal and bacterial taxa, such as Methanosphaera and Rothia, along with representative taxa from Coriobacteriaceae and Bacillaceae families and Chromatiales and Clostridiales orders (Figure 3).
Figure 3 Microbial taxa with Linear discriminant analysis (LDA) score greater than 2 present in hypercholesterolemic individuals (red) and controls (green). P-value <0.05.
Functional abundance prediction of bacterial groups in hypercholesterolemic and controls
Hypercholesterolemic individuals showed the predominant presence of biosynthesis pathways, including the super-pathway of aromatic amino acid biosynthesis (p=0.00316), chorismate biosynthesis I (p=0.00336), L-lysine biosynthesis I (p=0.00993), 1-4-dihydroxi-6-naphthoate biosynthesis I (p=0.022) and II (p = 0.010), menaquinol-8 biosynthesis II super-pathway (p = 0.014) and cob(II)yrinate a,c-diamide biosynthesis I (early cobalt insertion) (p=0.023). The first three are the pathways with the highest proportion in this study group (close to 0.9%). While in controls, degradation pathways are highlighted, such as 4-deoxy-L-threo-hex-4-enopyranuronate degradation of L-rhamnose degradation I (p=0.022), D-fructuronate degradation (p=0.037), phenylacetate degradation I (aerobic) (p=0.039), and phenylethylamine degradation pathway (p=0.041) (Figure 4).
Figure 4 Pathways obtained through PICRUSt2, which are significantly associated with hypercholesterolemic individuals and controls. The presence of biosynthetic pathways predominates in hypercholesterolemic individuals, while degradation pathways are more relevant in controls. P-value <0.05.
Interactions of the different bacterial groups present in hypercholesterolemic and controls
Co-occurrence networks allow to explore the interactions between different groups of microorganisms. This analysis was carried out at a general level (considering the two studied groups) as well as each study group separately. Figure 5 shows the graphic representation of the co-occurrence considering all the microorganisms of both studied groups.
Figure 5 Co-occurrence networks of microorganisms present in the study groups. Purple nodes represent the presence of microorganisms in both groups, while red nodes are mainly associated with hypercholesterolemic individuals and blue nodes with controls. With the same colors the interactions or edges are evidenced, associating these to the study groups. Continuous lines (edges) refer to co-occurrences that occur between microorganisms with a higher frequency than expected, while dashed lines correspond to interactions with less frequency than expected. P-value <0.05.
Through this analysis, 47 nodes were obtained (Figure 5), presenting co-occurrence interactions mainly associated with microorganisms from the groups Peptostreptococcaceae, Chromatiaceae, Barnesiellaceae, Bifidobacteriaceae, Comamonadaceae, Prevotellaceae, Marinilabiaceae, Fusobacteriaceae, Bacteroidaceae, Victivallaceae, Veillonellaceae, Dethiosulfovibrionaceae, Synergistaceae, Methanobacteriaceae, Erysipelotrichaceae, Clostridiaceae, Thermoanaerobacteraceae, Pasteurellaceae, Micrococcaceae and Leuconostocaceae. For these nodes, co-occurrence relations with a greater tendency in hypercholesterolemic individuals (red edges) and in other cases for controls (blue edges) can be observed. In addition, other microbial groups are presented, but these do not show evident interaction in this analysis.
When performed the analysis separately, Figure 6 showed the microorganisms that were mainly associated with the hypercholesterolemic group, thus obtaining 17 nodes, grouped into a single network, where high and low frequency interactions are present (continuous and dashed lines, respectively). In this study group, we found Prevotellaceae, Bifidobacteriaceae, Barnesiellaceae, Veillonellaceae, Comamonadaceae, Victivallaceae, Methanobacteriaceae, Synergistaceae, Thermoanaerobacteraceae, Micrococcaceae, Leuconostocaceae, Clostridiaceae, Fusobacteriaceae, Pasteurellaceae, Erysipelotrichaceae, Bacteroidaceae and Marinilabiaceae. In some cases, a microorganism can present both types of interactions depending on the group of microorganisms with which it is related, such as the Clostridiaceae group, which has a high frequency co-occurrence when it interacts with Pasteurellaceae or Thermoanaerobacteraceae, while the interaction is low when it is interacting with Fusobacteriaceae.
Figure 6 Co-occurrence networks of microorganisms present in hypercholesterolemic individuals. Continuous lines (edges) refer to co-occurrences that occur between microorganisms with a higher frequency than expected, while dashed lines correspond to interactions with less frequency than expected. P-value <0.05.
Regarding the controls, 17 nodes were also obtained, organizing them in two networks, with the most of microorganisms participating in one of them. In this study group we found Clostridiaceae, Thermoanaerobacteraceae, Victivallaceae, Erysipelotrichaceae, Veillonellaceae, Prevotellaceae, Comamonadaceae, Barnesiellaceae, Peptostreptococcaceae, Marinilabiaceae, Fusobacteriaceae, Bacteroidaceae, Synergistaceae, Dethiosulfovibrionaceae, Methanobacteriaceae, Pasteurellaceae and Micrococcaceae. As in the group of hypercholesterolemic individuals, there are co-occurrences with high and low frequencies of interaction (Figure 7).
Figure 7 Co-occurrence networks of microorganisms present in controls. Continuous lines (edges) refer to co-occurrences that occur between microorganisms with a higher frequency than expected, while dashed lines correspond to interactions with less frequency than expected. P-value <0.05.
Discussion
Considering that hypercholesterolemia is a main modifiable risk factor of CVDs, it is relevant to analyze a potential association of the different microbial taxa present in the intestinal tract with individuals classified as hypercholesterolemic and normocholesterolemic by LDL-C levels. Here, we profiled the stool microbiota of 57 adult subjects in order to identify possible associations between the composition of the microbiota and the LDL-C levels of these individuals.
By analyzing the sequences of the 16S rRNA gene, it was possible to establish that there were no significant differences at alpha and beta diversity between the studied groups (Faith’s PD index, p=0.396 and the unweighted UniFrac distance, p=0.133, together with the Jaccard distance). This contrasts with other works, where the difference in diversity between groups generally stands out, which can also be accompanied by differences in the abundance of microbial groups (Turnbaugh et al., 2009; Hansen et al., 2019).
Through the taxonomic analysis, highlights the presence of phylum Actinobacteria (p=0.025) in hypercholesterolemic individuals. This bacterial group being already indicated in other studies as part of the intestinal microbiota in humans, together with Firmicutes, Bacteroidetes and Proteobacteria phylum (Turnbaugh et al., 2009; Binda et al., 2018). Bacteroidetes and Firmicutes stand out predominantly, and our observations are consistent with this. The Bacteroidetes/Firmicutes relation has been related to the predisposition to certain diseases, with a lower abundance of Bacteroidetes being observed in obese individuals, in which bacteria belonging to the Firmicutes group abound (Ley et al., 2006; Gomes et al., 2018). Although in our study the BMI was increased in controls compared to hypercholesterolemic individuals (p=0.004), no differences were found between Bacteroidetes and Firmicutes (p=0.845 and p=0.279, respectively) by the taxonomic analysis. However, this remains controversial, since The Human Microbiome Project Consortium (2012) found a poor association between the Bacteroidetes/Firmicutes ratio and BMI. Furthermore, a higher abundance of Bacteroidetes in the gut microbiota has been associated with more industrialized populations compared to developing populations (Yatsunenko et al., 2012).
As expected, significant differences in the LDL-C (p<0.0001), triglycerides (p=0.003) and total cholesterol (p<0.0001) levels were observed in hypercholesterolemic compared to controls. Fu et al. (2015), suggest that the gut microbiome plays an important role in blood lipid levels, highlighting significant variations in triglyceride levels (2.74% variation) and HDL-C (3.83% variation). In contrast to our results, these authors showed a little impact of microbial groups on LDL-C (0.01% variation) and total cholesterol (0.01% variation) levels. Even so, the evidence of different studies reveals a close relation of the intestinal microbiota with lipids and lipoproteins, as well as with liver function and some of its alterations (Yu et al., 2019).
The intestinal microbiota obtains its nutrients from the host diet. Thus, Bacteroides, Roseburia, Bifidobacterium, Faecalibacterium, and Enterobacteria have been shown to be capable of fermenting undigested carbohydrates and producing short-chain fatty acids (SCFA), such as acetate, propionate, and butyrate (Al-Lahham et al., 2010). Both butyrate and propionate have low systemic concentrations whereas acetate levels are higher (Schoeler and Caesar, 2019). After absorption, SCFAs can induce lipogenesis and increase triglyceride stores through molecular pathways (Tokarek et al., 2021). Related to this, Granado-Serrano et al. (2019) founded high levels of SCFAs such as isobutyric and isovaleric acids in hypercholesterolemic individuals. In addition, the intestinal microbiota may interact with lipid metabolism through inhibition of hepatic cholesterol synthesis, redistribution of plasma cholesterol to the liver by the action of SCFAs, and/or deconjugation of bile acids by hydrolysis (Sonnenburg et al., 2006). In conjunction with evidence showing that oral supplementation with probiotics can reduce LDL-cholesterol concentration, lowering atherogenic indices and improving glycemic control (Ejtahed et al., 2012). However, we did not have access to data related to the diet of the participants, which we consider a limitation of this study that will be solved in future analyzes, in which the study of SCFAs might be included, considering that several microbial groups present in the studied individuals are related to the profile of these molecules.
In our study, LEfSe shows an enrichment of Faecalibacterium associated with controls individuals. Faecalibacterium prausnitzii is regarded as a butyrate producer in the gut, and it is used to treat dysbiosis in diseases that have an inflammatory component (Ganesan et al., 2018). This bacterial genus is recognized as a marker of gut health (Laursen et al., 2017) and changes in its abundance are associated with ulcerative colitis, colorectal cancer, and other diseases (Benevides et al., 2017). Other bacterial genera enriched in control individuals are Victivallis and Selenomonas. The genus Victivallis is recognized for being naturally present in gut (van Passel et al., 2011), though its abundance has been reported to vary with diet. Zhang et al. (2019), for instance, observed an increasing trend of abundance of this taxon in individuals treated with resistant starch, a prebiotic dietary fiber, also observing a decrease in LDL-C in the studied individuals. This would be concordant with the LDL-C levels and the presence of the genus Victivallis in the studied individuals. Regarding the genus Selenomonas, this has been mainly related to the oral health of children (Mashima et al., 2017), likewise it has been related to long-standing diseases that have an important inflammatory component, such as colorectal cancer and systemic lupus erythematosus (SLE) (Correa et al., 2017; Allali et al., 2018). In addition, some species such as Selenomonas noxia have been related to overweight (Goodson et al., 2009), also observed in the control group.
In hypercholesterolemic subjects, several ASVs were also differently abundant. Chromatiales and Clostridiales orders present genera and species that positively and negatively impact the development of diseases linked to the gut microbiota. Although, Chromatiales can be recognized as belonging to the Proteobacteria phylum, that are related to diverse pathologies including gut inflammation, irritable bowel syndrome and metabolic syndrome (Shin et al., 2015; Litvak et al., 2017). Regarding Clostridiales, it is a diverse group in which the genus Clostridium and some of its species are the most outstanding. Clostridium difficile, is known as a pathogenic agent associated with gut microbiota dysbiosis (Antharam et al., 2016), but other species such as Clostridium coccoides may also occur. Jamar et al. (2018), demonstrated the considerable presence of C. coccoides in individuals with metabolic syndrome, highlighting the effect of microbiota alterations on high levels of triglycerides in obese patients. Considering that hypercholesterolemic subjects in our study also show high levels of triglycerides, we hypothesize that would be related to the enrichment of these taxa in their gut microbiota.
The archaea of the genus Methanosphaera was also significantly associated to hypercholesterolemia in our study. These microorganisms have the capacity to produce methane in the presence of H2 and methanol (Gaci et al., 2014; Lurie-Weinberger and Gophna, 2015). It has been suggested that colonization by archaea is potentially detrimental to host health, due to alterations in intestinal metabolism and synthropic interactions with other microorganisms, which would stimulate the production of certain SCFAs (van de Pol et al., 2017; Canfora et al., 2019). Other studies demonstrate the importance of archaea such as Methanosphaera stadtmanae in the reduction of trimethylamine (TMA) at the gut level, which likely leads to a reduction in trimethylamine-N-oxide (TMAO) concentration (Gaci et al., 2014). In contrast to the groups of microorganisms indicated above, in relation to the genus Rothia, and the families Bacillaceae and Coriobacteriaceae, no data were found to establish a relationship between them and cholesterol metabolism.
The PICRUSt2 analysis showed multiple pathways differentially present in hypercholesterolemic individuals compared to controls. The most significant were the superpathway of aromatic aminoacid biosynthesis (p=0.00316), chorismate biosynthesis I (p=0.00336) and L-lysine biosynthesis I (p=0.00993). These three pathways are associated with the metabolism of amino acids. The production of amino acids and other peptides is part of the molecules generated by bacteria of the gut microbiota on a regular basis, so no evidence has been found that increased production of these molecules is related to hypercholesterolemia or associated pathologies. However, it has been possible to associate the levels of indoxyl sulfate levels with coronary atherosclerosis, as it has also been possible to associate p-cresyl sulfate with the prediction of some cardiovascular event, both metabolic derivatives of aromatic amino acids tryptophan and tyrosine, respectively, both produced by the intestinal microbiota (Wang and Zhao, 2018). Regarding the microorganisms related to the metabolism of amino acids, various groups such as Firmicutes, Bacteroidetes, among others, present a strong association with these biomolecules, such as the Clostridiaceae family, which presents a positive association with the super pathways associated with amino acids, being one of the groups of microorganisms identified in our study (Murga-Garrido et al., 2021).
Regarding pathways significantly associated with controls, we observed differences in 4-deoxy-L-threo-hex-4-enopyranuronate degradation (p=0.016), L-rhamnose degradation I (p=0.022) and D-fructuronate degradation (p=0.037). Unlike the group of hypercholesterolemic individuals, these pathways are mainly related to carbohydrate metabolism. Considering that gut microorganisms have the ability to digest simple and complex carbohydrates (Coker et al., 2021), it is likely that the significant presence of these pathways in this group would be related to a balanced diet in which carbohydrates associated with fiber predominate. However, no data about dietary intake was recorded in the studied individuals.
The application of co-occurrence analysis to the different microorganisms identified in the studied groups allows us to establish possible interactions (edges) among them, although these analyzes correspond to a bioinformatic approximation and require confirmation by other methodologies. In both hypercholesterolemic and control individuals, most co-occurring interactions occur more frequently than expected (continuous lines), compared to those that occur less frequently than expected (dashed lines). An example to highlight in terms of interactions with greater frequency than expected is the case of the Clostridiaceae family (Firmicutes) with other microbial groups (in individuals with hypercholesterolemia). Particularly, in this case, the negative association and absence of direct interaction between this family and the Bacteroidaceae (Bacteroidetes) family should be highlighted, which reflects a similar result described above for these two taxa (Faust et al., 2012). In general, the proper presence and interaction of the microbiome is crucial to maintain the health in hosts, from the first years of life to old age, allowing the development of adequate biological and immunological processes (Thomas et al., 2017).
Given the level of evidence that exists about the microorganisms identified and exposed in this article, and considering their possible association with cholesterol metabolism, would be interesting to investigate them more in depth to understand their impact on the production and regulation of cholesterol in the organism, especially considering that studies in mice have shown that the management of the microbiota together with classical cholesterol therapies can contribute to the prevention of diseases such as non-alcoholic fatty liver disease (NAFLD)-associated hepatocellular carcinoma (X. Zhang et al., 2021), as well as the association between gut microbiota and the effect of the lipid-lowering therapy using statins (Zimmermann et al., 2020); it is of great importance to obtain new evidence that contributes to the diagnosis and treatment of hypercholesterolemia, since its correspond to a one of main modifiable risk factors for cardiovascular diseases.
Data availability statement
The data presented in the study are deposited in the BioProject repository, accesion number PRJNA842179.
Ethics statement
The studies involving human participants were reviewed and approved by Comité Ético Científico, Universidad de La Frontera. The patients/participants provided their written informed consent to participate in this study.
Author contributions
NS and CM conceptualized and wrote the final manuscript. CM performed the experiments and obtained the results of the procedures. MR-H, RA-C, and SC-O analyzed the results using software. NS, GR, and CR. validated the results obtained. VM-C, KS, and PL carried out the formal analysis of the manuscript. NS, FL, and LS managed the resources used. NS, VM-C and PL obtained the financial support for the project. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by FONDECYT-ANID (11160364 to NS; 1211731 to VMC); FONDAP-ANID (15130011 to VMC); and by the ANID PhD fellowship (21181271 to CM).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
CVD, Cardiovascular disease; CAD, Coronary artery disease; LDL-C, Low-density lipoprotein cholesterol; PCR-DGGE, Polymerase chain reaction-gel electrophoresis with denaturing gradient; SBP, Systolic blood pressure; DBP, Diastolic blood pressure; HDL-C, High-density lipoprotein cholesterol; BMI, Body mass index; LDA, Linear Discriminant Analysis; PCoA, Principal Coordinates Analysis; PD, Phylogenetic diversity; SLE, Systemic lupus erythematosus; ASV, Amplicon Sequence Variant; SCFAs, Short-chain fatty acids; TMA, Trimethylamine; TMAO, Trimethylamine-N-oxide; NAFLD, Non-alcoholic fatty liver disease; NCEP, National Cholesterol Education Program; PCR, Polymerase chain reaction; DADA2, Divisive Amplicon Denoising Algorithm 2; LEfSe, Linear Discriminant Analysis Effect Size; QIIME2, Quantitative Insights Into Microbial Ecology 2; PICRUSt2, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2; STAMP, Statistical Analysis of Metagenomic Profiles.
References
Al-Lahham, S. H., Peppelenbosch, M. P., Roelofsen, H., Vonk, R. J., Venema, K. (2010). Biological effects of propionic acid in humans; metabolism, potential applications and underlying mechanisms. Biochim. Biophys. Acta 1801 (11), 1175–1183. doi: 10.1016/j.bbalip.2010.07.007
Allali, I., Boukhatem, N., Bouguenouch, L., Hardi, H., Boudouaya, H. A., Cadenas, M. B., et al. (2018). Gut microbiome of Moroccan colorectal cancer patients. Med. Microbiol. Immunol. 207 (3-4), 211–225. doi: 10.1007/s00430-018-0542-5
Antharam, V. C., McEwen, D. C., Garrett, T. J., Dossey, A. T., Li, E. C., Kozlov, A. N., et al. (2016). An integrated metabolomic and microbiome analysis identified specific gut microbiota associated with fecal cholesterol and coprostanol in clostridium difficile infection. PloS One 11 (2), e0148824. doi: 10.1371/journal.pone.0148824
Antman, E. M., Loscalzo, J. (2016). Precision medicine in cardiology. Nat. Rev. Cardiol. 13 (10), 591–602. doi: 10.1038/nrcardio.2016.101
Bäckhed, F., Ley, R. E., Sonnenburg, J. L., Peterson, D. A., Gordon, J. I. (2005). Host-bacterial mutualism in the human intestine. Science 307 (5717), 1915–1920. doi: 10.1126/science.1104816
Baigent, C., Keech, A., Kearney, P. M., Blackwell, L., Buck, G., Pollicino, C., et al. (2005). Efficacy and safety of cholesterol-lowering treatment: prospective meta-analysis of data from 90,056 participants in 14 randomised trials of statins. Lancet 366 (9493), 1267–1278. doi: 10.1016/s0140-6736(05)67394-1
Benevides, L., Burman, S., Martin, R., Robert, V., Thomas, M., Miquel, S., et al. (2017). New insights into the diversity of the genus faecalibacterium. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.01790
Binda, C., Lopetuso, L. R., Rizzatti, G., Gibiino, G., Cennamo, V., Gasbarrini, A. (2018). Actinobacteria: A relevant minority for the maintenance of gut homeostasis. Dig. Liver. Dis. 50 (5), 421–428. doi: 10.1016/j.dld.2018.02.012
Bourgin, M., Labarthe, S., Kriaa, A., Lhomme, M., Gérard, P., Lesnik, P., et al. (2020). Exploring the bacterial impact on cholesterol cycle: A numerical study. Front. Microbiol. 11. doi: 10.3389/fmicb.2020.01121
Bradley, B. (2020) How to create co-occurrence networks with the r packages ‘cooccur’ and ‘visNetwork’. Available at: https://medium.com/analytics-vidhya/how-to-create-co-occurrence-networks-with-the-r-packages-cooccur-and-visnetwork-f6e1ceb1c523.
Canfora, E. E., Meex, R. C. R., Venema, K., Blaak, E. E. (2019). Gut microbial metabolites in obesity, NAFLD and T2DM. Nat. Rev. Endocrinol. 15 (5), 261–273. doi: 10.1038/s41574-019-0156-z
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7 (5), 335–336. doi: 10.1038/nmeth.f.303
Coker, J. K., Moyne, O., Rodionov, D. A., Zengler, K. (2021). Carbohydrates great and small, from dietary fiber to sialic acids: How glycans influence the gut microbiome and affect human health. Gut. Microbes 13 (1), 1–18. doi: 10.1080/19490976.2020.1869502
Correa, J. D., Calderaro, D. C., Ferreira, G. A., Mendonca, S. M., Fernandes, G. R., Xiao, E., et al. (2017). Subgingival microbiota dysbiosis in systemic lupus erythematosus: Association with periodontal status. Microbiome 5 (34), 1–13. doi: 10.1186/s40168-017-0252-z
Douglas, G. M., Maffei, V. J., Zaneveld, J. R., Yurgel, S. N., Brown, J. R., Taylor, C. M., et al. (2020). PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 38 (6), 685–688. doi: 10.1038/s41587-020-0548-6
Ejtahed, H. S., Mohtadi-Nia, J., Homayouni-Rad, A., Niafar, M., Asghari-Jafarabadi, M., Mofid, V. (2012). Probiotic yogurt improves antioxidant status in type 2 diabetic patients. Nutrition 28 (5), 539–543. doi: 10.1016/j.nut.2011.08.013
Faust, K., Sathirapongsasuti, J. F., Izard, J., Segata, N., Gevers, D., Raes, J., et al. (2012). Microbial co-occurrence relationships in the human microbiome. PloS Comput. Biol. 8 (7), e1002606. doi: 10.1371/journal.pcbi.1002606
Fu, J., Bonder, M. J., Cenit, M. C., Tigchelaar, E. F., Maatman, A., Dekens, J. A., et al. (2015). The gut microbiome contributes to a substantial proportion of the variation in blood lipids. Circ. Res. 117 (9), 817–824. doi: 10.1161/circresaha.115.306807
Gaci, N., Borrel, G., Tottey, W., O’Toole, P. W., Brugere, J. F. (2014). Archaea and the human gut: new beginning of an old story. World J. Gastroenterol. 20 (43), 16062–16078. doi: 10.3748/wjg.v20.i43.16062
Ganesan, K., Chung, S. K., Vanamala, J., Xu, B. (2018). Causal relationship between diet-induced gut microbiota changes and diabetes: A novel strategy to transplant faecalibacterium prausnitzii in preventing diabetes. Int. J. Mol. Sci. 19 (12), 1–28. doi: 10.3390/ijms19123720
Gomes, A. C., Hoffmann, C., Mota, J. F. (2018). The human gut microbiota: Metabolism and perspective in obesity. Gut. Microbes 9 (4), 308–325. doi: 10.1080/19490976.2018.1465157
Goodson, J. M., Groppo, D., Halem, S., Carpino, E. (2009). Is obesity an oral bacterial disease? J. Dent. Res. 88 (6), 519–523. doi: 10.1177/0022034509338353
Granado-Serrano, A. B., Martin-Gari, M., Sanchez, V., Riart Solans, M., Berdun, R., Ludwig, I. A., et al. (2019). Faecal bacterial and short-chain fatty acids signature in hypercholesterolemia. Sci. Rep. 9 (1), 1772. doi: 10.1038/s41598-019-38874-3
Grundy, S. M., Cleeman, J. I., Merz, C. N., Brewer, H. B., Jr., Clark, L. T., Hunninghake, D. B., et al. (2004). Implications of recent clinical trials for the national cholesterol education program adult treatment panel III guidelines. Circulation 110 (2), 227–239. doi: 10.1161/01.CIR.0000133317.49796.0E
Hansen, M. E. B., Rubel, M. A., Bailey, A. G., Ranciaro, A., Thompson, S. R., Campbell, M. C., et al. (2019). Population structure of human gut bacteria in a diverse cohort from rural Tanzania and Botswana. Genome Biol. 20 (16), 1–21. doi: 10.1186/s13059-018-1616-9
Human Microbiome Project Consortium. (2012). Structure, function and diversity of the healthy human microbiome. Nature 486 (7402), 207–214. doi: 10.1038/nature11234
Jamar, G., Santamarina, A. B., Dias, G. C., Masquio, D. C. L., de Rosso, V. V., Pisani, L. P. (2018). Relationship between fatty acids intake and clostridium coccoides in obese individuals with metabolic syndrome. Food Res. Int. 113, 86–92. doi: 10.1016/j.foodres.2018.07.002
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast,, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res 41 (1), e1 (1–11). doi: 10.1093/nar/gks808
Laursen, M. F., Laursen, R. P., Larnkjaer, A., Molgaard, C., Michaelsen, K. F., Frokiaer, H., et al. (2017). Faecalibacterium gut colonization is accelerated by presence of older siblings. mSphere 2 (6), 1–16. doi: 10.1128/mSphere.00448-17
Ley, R. E., Turnbaugh, P. J., Klein, S., Gordon, J. I. (2006). Human gut microbes associated with obesity. Nature 444 (7122), 1022–1023. doi: 10.1038/4441022a
Litvak, Y., Byndloss, M. X., Tsolis, R. M., Baumler, A. J. (2017). Dysbiotic proteobacteria expansion: A microbial signature of epithelial dysfunction. Curr. Opin. Microbiol. 39, 1–6. doi: 10.1016/j.mib.2017.07.003
Lurie-Weinberger, M. N., Gophna, U. (2015). Archaea in and on the human body: Health implications and future directions. PloS Pathog. 11 (6), e1004833. doi: 10.1371/journal.ppat.1004833
Martinez, I., Wallace, G., Zhang, C., Legge, R., Benson, A. K., Carr, T. P., et al. (2009). Diet-induced metabolic improvements in a hamster model of hypercholesterolemia are strongly linked to alterations of the gut microbiota. Appl. Environ. Microbiol. 75 (12), 4175–4184. doi: 10.1128/AEM.00380-09
Mashima, I., Theodorea, C. F., Thaweboon, B., Thaweboon, S., Scannapieco, F. A., Nakazawa, F. (2017). Exploring the salivary microbiome of children stratified by the oral hygiene index. PloS One 12 (9), e0185274. doi: 10.1371/journal.pone.0185274
Mazidi, M., Rezaie, P., Kengne, A. P., Mobarhan, M. G., Ferns, G. A. (2016). Gut microbiome and metabolic syndrome. Diabetes Metab. Syndr. 10 (2 Suppl 1), S150–S157. doi: 10.1016/j.dsx.2016.01.024
Mendelsohn, A. R., Larrick, J. W. (2013). Dietary modification of the microbiome affects risk for cardiovascular disease. Rejuv. Res. 16 (3), 241–244. doi: 10.1089/rej.2013.1447
Murga-Garrido, S. M., Hong, Q., Cross, T. L., Hutchison, E. R., Han, J., Thomas, S. P., et al. (2021). Gut microbiome variation modulates the effects of dietary fiber on host metabolism. Microbiome 9 (1), 117. doi: 10.1186/s40168-021-01061-6
Parks, D. H., Tyson, G. W., Hugenholtz, P., Beiko, R. G. (2014). STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30 (21), 3123–3124. doi: 10.1093/bioinformatics/btu494
Rebolledo, C., Cuevas, A., Zambrano, T., Acuna, J. J., Jorquera, M. A., Saavedra, K., et al. (2017). Bacterial community profile of the gut microbiota differs between hypercholesterolemic subjects and controls. BioMed. Res. Int. 2017, 8127814. doi: 10.1155/2017/8127814
Schoeler, M., Caesar, R. (2019). Dietary lipids, gut microbiota and lipid metabolism. Rev. Endocrine. Metab. Disord. 20 (4), 461–472. doi: 10.1007/s11154-019-09512-0
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12 (6), R60. doi: 10.1186/gb-2011-12-6-r60
Sender, R., Fuchs, S., Milo, R. (2016). Are we really vastly outnumbered? revisiting the ratio of bacterial to host cells in humans. Cell 164 (3), 337–340. doi: 10.1016/j.cell.2016.01.013
Shin, N. R., Whon, T. W., Bae, J. W. (2015). Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 33 (9), 496–503. doi: 10.1016/j.tibtech.2015.06.011
Sonnenburg, J. L., Chen, C. T. L., Gordon, J. I. (2006). Genomic and metabolic studies of the impact of probiotics on a model gut symbiont and host. PloS Biol. 4 (12), e413. doi: 10.1371/journal.pbio.0040413
Thomas, S., Izard, J., Walsh, E., Batich, K., Chongsathidkiet, P., Clarke, G., et al. (2017). The host microbiome regulates and maintains human health: A primer and perspective for non-microbiologists. Cancer Res. 77 (8), 1783–1812. doi: 10.1158/0008-5472.Can-16-2929
Tokarek, J., Gadzinowska, J., Młynarska, E., Franczyk, B., Rysz, J. (2021). What is the role of gut microbiota in obesity prevalence? a few words about gut microbiota and its association with obesity and related diseases. Microorganisms 10 (1), 52. doi: 10.3390/microorganisms10010052
Turnbaugh, P. J., Hamady, M., Yatsunenko, T., Cantarel, B. L., Duncan, A., Ley, R. E., et al. (2009). A core gut microbiome in obese and lean twins. Nature 457 (7228), 480–484. doi: 10.1038/nature07540
van de Pol, J. A., van Best, N., Mbakwa, C. A., Thijs, C., Savelkoul, P. H., Arts, I. C., et al. (2017). Gut colonization by methanogenic archaea is associated with organic dairy consumption in children. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.00355
van Passel, M. W., Kant, R., Palva, A., Lucas, S., Copeland, A., Lapidus, A., et al. (2011). Genome sequence of victivallis vadensis ATCC BAA-548, an anaerobic bacterium from the phylum lentisphaerae, isolated from the human gastrointestinal tract. J. Bacteriol. 193 (9), 2373–2374. doi: 10.1128/JB.00271-11
Vázquez-Baeza, Y., Pirrung, M., Gonzalez, A., Knight, R. (2013). EMPeror: A tool for visualizing high-throughput microbial community data. Gigascience 2 (1), 16. doi: 10.1186/2047-217x-2-16
Veech, J. A. (2013). A probabilistic model for analysing species co-occurrence. Global Ecol. Biogeogr. 22 (2), 252–260. doi: 10.1111/j.1466-8238.2012.00789.x
Wang, Z., Zhao, Y. (2018). Gut microbiota derived metabolites in cardiovascular health and disease. Protein Cell 9 (5), 416–431. doi: 10.1007/s13238-018-0549-0
Wostmann, B. S. (1973). Intestinal bile acids and cholesterol absorption in the germfree rat. J. Nutr. 103 (7), 982–990. doi: 10.1093/jn/103.7.982
Wu, H., Tremaroli, V., Backhed, F. (2015). Linking microbiota to human diseases: A systems biology perspective. Trends Endocrinol. Metab. 26 (12), 758–770. doi: 10.1016/j.tem.2015.09.011
Yatsunenko, T., Rey, F. E., Manary, M. J., Trehan, I., Dominguez-Bello, M. G., Contreras, M., et al. (2012). Human gut microbiome viewed across age and geography. Nature 486 (7402), 222–227. doi: 10.1038/nature11053
Yun, K. E., Kim, J., Kim, M. H., Park, E., Kim, H. L., Chang, Y., et al. (2020). Major lipids, apolipoproteins, and alterations of gut microbiota. J. Clin. Med. 9 (5), 1–16. doi: 10.3390/jcm9051589
Yu, Y., Raka, F., Adeli, K. (2019). The role of the gut microbiota in lipid and lipoprotein metabolism. J. Clin. Med. 8 (12), 1–20. doi: 10.3390/jcm8122227
Yusuf, S., Mehta, S. R., Chrolavicius, S., Afzal, R., Pogue, J., Granger, C. B., et al. (2006). Comparison of fondaparinux and enoxaparin in acute coronary syndromes. N. Engl. J. Med. 354 (14), 1464–1476. doi: 10.1056/NEJMoa055443
Zhang, X., Coker, O. O., Chu, E. S., Fu, K., Lau, H. C. H., Wang, Y. X., et al. (2021). Dietary cholesterol drives fatty liver-associated liver cancer by modulating gut microbiota and metabolites. Gut 70 (4), 761–774. doi: 10.1136/gutjnl-2019-319664
Zhang, L., Ouyang, Y., Li, H., Shen, L., Ni, Y., Fang, Q., et al. (2019). Metabolic phenotypes and the gut microbiota in response to dietary resistant starch type 2 in normal-weight subjects: A randomized crossover trial. Sci. Rep. 9 (1), 4736. doi: 10.1038/s41598-018-38216-9
Keywords: gut microbiota, LDL cholesterol, LEfSe analysis, microbial signature, 16S rRNA sequencing
Citation: Morales C, Rojas G, Rebolledo C, Rojas-Herrera M, Arias-Carrasco R, Cuadros-Orellana S, Maracaja-Coutinho V, Saavedra K, Leal P, Lanas F, Salazar LA and Saavedra N (2022) Characterization of microbial communities from gut microbiota of hypercholesterolemic and control subjects. Front. Cell. Infect. Microbiol. 12:943609. doi: 10.3389/fcimb.2022.943609
Received: 17 May 2022; Accepted: 31 October 2022;
Published: 29 November 2022.
Edited by:
George Seghal Kiran, Pondicherry University, IndiaReviewed by:
Zongfu Hu, Inner Mongolia Minzu University, ChinaRajeev Meora, Inha University, South Korea
Copyright © 2022 Morales, Rojas, Rebolledo, Rojas-Herrera, Arias-Carrasco, Cuadros-Orellana, Maracaja-Coutinho, Saavedra, Leal, Lanas, Salazar and Saavedra. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nicolás Saavedra, nicolas.saavedra@ufrontera.cl