Distribution of C30 carotenoid biosynthesis genes suggests habitat adaptation function in insect-adapted and nomadic Lactobacillaceae

doi:10.21203/rs.3.rs-4637278/v1

Download PDF

Article

Distribution of C30 carotenoid biosynthesis genes suggests habitat adaptation function in insect-adapted and nomadic Lactobacillaceae

https://doi.org/10.21203/rs.3.rs-4637278/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Carotenoids are membrane-bound pigments that are essential for photosynthesizing plants and algae. These specialized metabolites are widely applied in food, feed and cosmetics because of their antioxidant and anti-inflammatory properties. The production of carotenoids, particularly C30 forms, has also been documented in a few non-photosynthetic prokaryotes. However, our understanding of the function, distribution and ecology of these compounds beyond photosynthesizing organisms is limited. In this study, we performed an eco-evolutionary analysis of terpenoid biosynthetic gene clusters in the Lactobacillaceae family, an important family of beneficial bacteria with diverse lifestyles and habitats. Hereto, 4203 dereplicated genomes were screened for terpenoid biosynthesis genes, resulting in detection of crtMN genes in 28/361 (7.7%) species in 14/34 (41.2%) genera. These genes encode the key enzymes that transform two farnesyl pyrophosphate molecules into the C30 carotenoid 4,4’-diaponeurosporene. These crtMN genes appeared to be convergently gained within Fructilactobacillus, and horizontally transferred across species and genera, e.g. from Lactiplantibacillus to Levilactobacillus. Subsequently, in our in-house Lactobacillaceae culture collection from various habitats (n = 575), the phenotype was confirmed in 87% of the isolates predicted to have the crtMN genes (27/31). Nomadic and insect-adapted species, particularly those isolated from vegetable fermentations, such as Lactiplantibacillus, and floral habitats, such as Fructilactobacillus, were capable of C30 carotenoid biosynthesis, while vertebrate-associated species, such as those from the vagina, lacked this trait. This habitat association aligned with the observations that C30 carotenoid-producing strains were more resistant to oxidative and UV-stress. Taken together, this study revealed that C30 carotenoid biosynthesis plays a role in habitat adaptation and is scattered across Lactobacillaceae in line with this habitat adaptation.

Biological sciences/Microbiology/Environmental microbiology

Biological sciences/Ecology/Microbial ecology

Biological sciences/Microbiology/Microbial genetics/Bacterial genes

Carotenoids are membrane-bound specialized metabolites found in all photosynthesizing organisms, such as plants, algae and cyanobacteria. They play vital roles in protecting producer cells against excess light energy and scavenging reactive oxygen species^1–3. However, these molecules also perform functions that are unrelated to photosynthesis, such as being precursors to phytohormones^4,5 or attracting pollinators and promoting seed dispersal by accumulating in flowers and fruits⁶. Moreover, carotenoids are found in non-photosynthetic organisms, for example in the pathogenic bacterium Staphylococcus aureus, where the carotenoid staphyloxanthin increases virulence by protecting the bacteria against reactive oxygen species and inactivating host neutrophils⁷. In leaf-dwelling non-photosynthetic bacteria, such as Clavibacter, Pseudomonas and Methylobacteria, carotenoids protect bacterial cells and, in some cases, the plant host from UV radiation^8,9. Animals require but are generally unable to synthesize carotenoids and consequently supplementation depends on food intake, with the notable exception of aphids, which have acquired the ability to produce carotenoids as fungal carotenoid biosynthetic genes have been taken up into their genome¹⁰. In humans, dietary intake of carotenoids is essential for vitamin A biosynthesis¹¹ and is associated with a reduction in inflammatory markers¹² and a reduced risk of prostate cancer and breast cancer¹³, although also opposite associations have been found such as for lung cancer¹⁴. Because mostly positive health benefits have been found, carotenoids are increasingly used as functional ingredients in food and feed, where they can also have additional technological benefits, for instance as natural pigments and antioxidants¹.

The best-studied carotenoids are those found in photosynthetic organisms. These specialized metabolites are mostly C40-carotenoids, which have 40 carbon atoms within their backbone¹⁵. More recently, the full scope of carotenoid diversity and evolution has gained increasing attention^15,16. Phylogenetic analyses revealed that carotenoid biosynthesis is an ancient process that evolved prior to photosynthesis, likely under increased UV radiation conditions^16–18. According to Santana-Molina et al.¹⁹, the earliest evolved group of carotenoids might be composed of 30 carbon atoms, or triterpenoids. These specialized compounds are found in plants under various forms, such as triterpenoid saponins²⁰ or oleanolic acids²¹, and they are often involved in the plant’s defense system. However, C30 carotenoid biosynthesis originated in prokaryotes^16,18 and has been detected in the photosynthesizing bacteria Heliobacter²², and several nonphotosynthetic bacteria such as Lactiplantibacillus plantarum²³, Bacillus subtilis²⁴, Enterococcus faecium (formerly Streptococcus faecium)²⁵, Staphylococcus aureus²⁶ and Methylobacterium rhodinum (formerly Pseudomonas rhodos)²⁷. In addition to these phenotypic observations, Santana and colleagues¹⁸ showed that both 4,4’-diapophytoene (synthesis encoded by the crtM gene) as well as squalene (encoded by Sqs or the HpnCDE cluster) can be precursors to C30 carotenoids, and that these pathways are now found scattered across prokaryotes, ranging from Firmicutes to Planctomycetes to Archaea and found to be mobile via horizontal gene transfer.

While carotenoid biosynthesis originated in, and is now scattered across nonphotosynthesizing prokaryotes, the function and ecological role of these specialized metabolites in these organisms have not been systematically studied. Within the Firmicutes, Lactobacillaceae are an interesting family of nonphotosynthesizing bacteria, with great potential for use in food, feed and pharma^28,29. Two phenotypical screenings of Lactobacillaceae, focusing on isolates from fermented foods, identified Lp. plantarum as a C30 carotenoid producer via the 4,4'-diapophytoene pathway^23,30. However, the genomic and ecological diversity of the biosynthetic and phenotypic potential of this family has not yet been investigated.

In this study, an in-depth pangenome and evolutionary analysis of the presence of terpenoid biosynthetic gene clusters was performed on a large dataset of 4203 unique genomes of the Lactobacillaceae family. These genomic analyses were complemented with phenotypic confirmation using a biobank of Lactobacillaceae from various habitats by high-throughput screening based on absorbance maxima and high-performance liquid chromatography (HPLC) and high-resolution quadrupole time-of-flight mass spectrometry (Q-TOF MS). Additionally, the role of C30 carotenoids in UV and oxidative stress resilience was assessed for selected strains. Finally, associations between the presence of carotenoid biosynthesis genes in a species and its lifestyle were established, using assigned lifestyles from literature^31,32, and new isolation and metabolic analyses from our in-house library of strains.

Pangenome and phylogenetic analysis of the terpenoid biosynthetic gene clusters.

All publicly available Lactobacillaceae genomes were downloaded from the Genome Taxonomy Database (GTDB, gtdb.ecogenomic.org, version r214). With checkM, any incomplete (< 90%) or contaminated (> 5%) genomes were excluded³³. This resulted in a dataset of 6259 public genomes. To avoid pseudoreplication, the sample function of SCARAP (github.com/SWittouck/SCARAP) was used with an average nucleotide identity (ANI) cutoff of 99.99%, retaining only genomes with ANI values below this cutoff. Second, the pan genomes of this family were inferred using the SCARAP tool³⁴. As a start, orthogroups of detected genes, i.e.4,4’-diapophytoene desaturase (encoded by the crtN gene) and of 4,4’-diapophytoene synthase (encoded by the crtM gene) were determined (Table 1).³³ For reference, the genes from two experimentally confirmed species, Lp. plantarum WCFS1²³ and Latilactobacillus fragifolii AMBP162 (experimentally confirmed in this study), were used as queries. Based on the results of this pangenome analysis, 37 in-house isolates were selected for whole genome sequencing, and were included in the subsequent analysis. Finally, all genomes were checked for other terpene biosynthesis genes using Antismash 6.0³⁵ focusing on terpenes and terpene pathways.

Table 1

**References to experimentally confirmed** **crtMN** **genes** used to identify similar sequences in the pangenome.
Strain	crtM 4,4’-diapophytoene synthase (882 nt)	crtN 4,4’-diapophytoene desaturase (1497 nt)	Reference
Lactiplantibacillus plantarum WCFS1	NC_004567.2_2720	NC_004567.2_2719	²³
Latilactobacillus fragifolii AMBP162	NZ_CAJOCF010000002.1_110	NZ_CAJOCF010000002.1_111	Carotenoid biosynthesis was first observed during the description of the species³⁶ and further validated in this study

The core tree of all Lactobacillaceae was constructed following the method described in Eilers et al.³⁷. The core genome was inferred using SCARAP, which consists of 296 core genes. Protein sequences were extracted, aligned with MAFFT and trimmed with trimAL with a gap threshold of 10%. Aligned and trimmed core proteins were subsequently inferred with IQ-TREE using LG + F + G4.

First, the prevalence of crtMN genes within all species of the family Lactobacillaceae was calculated and integrated with lifestyle data from Zheng et al., 2020³², and metadata from the GTDB using tidygenomes (github.com/SWittouck/tidygenomes) based on ggtree packages in R. Genera were collapsed when no species of this genus contained the crtMN genes.

Second, to examine the phylogeny of the crtM and crtN within the Lactobacillaceae family, a gene tree was inferred at the amino acid level using a similar procedure to the species tree. Since the crtM and crtN tree were highly similar, the crtN tree was used as a model due to its larger size and thus resolution. To reduce the number of branches, sequences were first clustered with cd-hit³⁸ with a 95% similarity threshold. In instances where sequences from different species were present in the same CD-hit cluster, they are shown in the adjacent table. The biosynthetic gene clusters obtained from antiSMASH were visualized using Bigscape³⁹ based on dereplicated genomes, and the clusters were mapped onto the crtN tree.

Strains used in this study

Phenotypic characterization started from the crtMN carriers within 575 in-house Lactobacillaceae isolates and three publicly available strains (Supplementary Table S1). These strains were previously isolated from vegetable fermentations³⁹; the human vagina⁴⁰; the human respiratory tract⁴¹, the phyllosphere³⁶; anthosphere⁴²; and liquid compost fermentations⁴³. Forty-eight isolates, taxonomically related to carotenoid producers based on the pangenome analysis, were phenotypically screened. To complement the publicly available data, the genomes of 25 crtMN and 17 non-crtMN containing in-house isolates were sequenced and included in the pangenome analysis (Supplementary Table S1, study number PRJEB57255). UV and oxidative stress assays were performed using nine Lp. plantarum strains, as indicated in Supplementary Table S1.

C30 carotenoid extraction and identification

Lipophilic compounds were extracted, and their absorption spectra were measured based on the methods of Garrido-Fernández and colleagues²³. In brief, cells were harvested from a 50 ml overnight culture in Weissella Medium Broth (WMB) by centrifugation for 15 min at 2000 × g. The cells were subsequently washed with 50 ml of sterile distilled water. Afterwards, 10 ml of N,N-dimethylformamide was added to the washed cells, which were incubated for 15 min at 65°C. The cell debris was separated by centrifugation at 3000 × g for 10 min, after which the supernatant was transferred to a separator funnel. The extraction of the remaining cell debris with N,N-dimethylformamide was repeated four times. All the extracts were pooled and mixed with 100 ml of diethyl ether, and 10% NaCl was added to aid in the separation of the liquid phases. The organic phase was dried with anhydrous Na₂SO₄, followed by solvent evaporation in a rotary evaporator. The resulting residue was dissolved in 2 ml of methanol/tert-butyl methyl ether (1:1 v:v) containing 1% BHT. The absorption spectrum was measured between 550 nm and 300 nm using a spectrophotometer to determine the characteristic absorption maxima. Carotenoids were further purified and identified using high-performance liquid chromatography (HPLC) coupled inline to a diode array detector (DAD) and high-resolution quadrupole time-of-flight mass spectrometry (Q-TOF MS), a procedure conducted by the RIC group. For high-throughput identification, a similar procedure was used with different extraction solvents and without further purification. The washed cells were extracted with 2 ml of molecular biology grade ethanol, vortexed thoroughly and incubated for 20 min at 65°C. The cell debris was separated by centrifugation at 8603 × g for 5 min, after which the supernatant was transferred to a 15 ml tube. The ethanol extraction step was repeated once on the remaining cell debris. Both extracts were pooled, and 2 ml of heptane was added and mixed vigorously for 1 min. Next, 2 ml of distilled water was added and mixed. The hydrophilic and lipophilic phases were then separated by centrifuging the tubes for 5 min at 4000 × g. The top layer (lipophilic phase) was cautiously transferred to a clean 2 ml tube and dried with anhydrous Na₂SO₄. Finally, the absorption peaks between 300 and 550 nm were measured using a spectrophotometer to detect characteristic absorption peaks.

Oxidative and UV stress resistance assays

Cells were harvested by centrifugation (10 min at 1500 × g) from a 10 mL two-day culture, grown in MRS medium at 28°C with shaking, washed with sterile phosphate-buffered saline (PBS) and resuspended to an optical density (OD) of 0.16 at 600 nm. For oxidative stress, 100 µL of diluted culture was added to a suspension containing a final concentration of 0.03% H₂O₂ and incubated while shaking for 10, 20, or 30 min. For UV stress, 500 µL of each suspension and 8 mL of sterile PBS were dispensed into small Petri plates and placed on an orbital shaker with a gentle swirling motion inside a laminar air flow cabinet, after which the lids were removed. The swirling Petri plates were exposed to UV treatment, and samples were collected at four time points: before UV exposure, after 10 sec, after 20 sec, and after 30 sec of UV exposure. The number of colony-forming units (CFUs) was determined for all the samples through serial dilution and plating out in triplicate on MRS agar. The tested strains were tested in triplicates and the entire assay was repeated three times. The rate of decrease in CFU counts was determined as a measure of susceptibility to oxidative and UV stress and involved transforming the CFU counts to a logarithmic scale and fitting a linear regression to infer the slope of the decrease over time.

Association of C30 carotenoid biosynthesis capacity and habitat

To further associate the carotenoids and habitats, multiple statistical analyses were performed on the dereplicated dataset in R with the packages dunn.test, caper, nlme, and geiger to assess the correlation between crtMN prevalence and lifestyle and genome size, taking phylogenetic effects into account. A gene was considered a core trait of a species when more than 95% of the genomes within a species had this trait and an accessory trait when the prevalence was between 5% and 95%.

Data availability

The publicly available Lactobacillaceae genomes were downloaded from the Genome Taxonomy Database (GTDB, version r214). All newly sequenced genomes from in-house isolates are available in the ENA database in project PRJEB57255 using the accession numbers given in Supplementary Table S1.

C30 carotenoid biosynthesis genes are scattered across the Lactobacillaceae family and show frequent horizontal transfer

To explore the carotenoid biosynthetic gene cluster potential within the Lactobacillaceae family, the presence of all known enzymes involved in carotenoid biosynthesis was evaluated. For this purpose, 6297 publicly available genomes were screened. Due to the high similarity among genomes in this dataset, genomes were dereplicated based on pairwise ANI < 99.99%, resulting in 4179 unique genomes, from which a pangenome was constructed. Enzymes involved in the 4,4'-diapophytoene pathway for the production of C30 carotenoids were present; 4,4'-diapophytoene synthase (encoded by the crtM gene) and 4,4’-diapophytoene desaturase (crtN) (Fig. 1). The orthogroups to which each gene belonged were the same in the experimentally confirmed carotenoid-producers, Lp. plantarum WCFS1²³ and Lt. fragifolii AMBP162^T (phenotype first observed during the description of the species³⁶ and confirmed in this paper), validating the use of these orthogroups to detect crtMN genes across the Lactobacillaceae pangenome. Additionally, analysis with antiSMASH confirmed that no other known clusters for carotenoid biosynthesis, such as the squalene pathway, were present in the Lactobacillaceae. The crtMN C30 carotenoid biosynthesis gene families showed to be scattered across 28 species from 14 genera out of a total of 361 species and 34 genera analyzed within the Lactobacillaceae family (i.e., 7.7% species and 41.2% genus prevalence) (Fig. 2). We selected 37 in-house strains (Supplementary Table S1) for genome sequencing to complement the public dataset based on the likelihood for presence of crtMN genes (29 strains) and underrepresented species (8 strains). After genomic dereplication, 24 remained (Supplementary Table S1) and were included in the subsequent analysis (4203 unique genomes). The detected crtMN genes clustered into two clades (Fig. 3). One clade included Lactiplantibacillus, Fructilactobacillus, Latilactobacillus and Companilactobacillus, while the other exhibited broader diversity, including genera such as Leuconostoc, Oenococcus, Holzapfelia and other taxa within Fructilactobacillus. Notably, the gene trees for crtM and crtN were highly similar, indicating that these genes were typically transferred or inherited together. The crtN tree was preferred for visualization because it was more accurate due to the greater length of the gene (i.e. 499 amino acids). Notably, the Fructilactobacillus genus was present in both clades, indicating two independent acquisition events of the crtMN genes in this genus and showing that this trait had been gained convergently in this genus. The crtMN phylogenetic trees did not align with the species tree based on the core genome (Fig. 2), suggesting that these genes were not only inherited through vertical transmission but also horizontally transferred across species and even genera. Clear examples of recent horizontal gene transfer (HGT) events were observed within clusters of multiple species that contain nearly identical crtMN genes (up to 100% similarity at the protein level), for example amongst Lp. plantarum, Levilactobacillus buchneri, Levilactobacillus brevis, Lactiplantibacillus pentosus and Pediococcus pentosaceus, as well as amongst several Leuconostoc species. Of note, one cluster present in Apilactobacillus ozensis and Apilactobacillus xinyiensis seemed to contain two genes that were both annotated as crtN.

Identification of C30 carotenoid-producing species in Lactobacillaceae biobank

Having shown the taxonomic spread of C30 carotenoid biosynthesis genes within Lactobacillaceae, we subsequently aimed to phenotypically substantiate the carotenoid biosynthesis capacity in strains containing crtMN genes. First, C30 carotenoid biosynthesis was confirmed in Latilactobacillus fragifolii AMBP162^T, originally isolated from the phyllosphere of a strawberry plant³⁶ and compared to that of the known C30 carotenoid producer Lp. plantarum WCFS1²³. The pellets of both strains appeared yellow, providing the first indication of carotenoid biosynthesis (Fig. 4A). The extract was purified using HPLC coupled inline to a diode-array detector (DAD) and high-resolution quadrupole time-of-flight mass spectrometry (Q-TOF MS). The highest peak in the UV chromatogram of the HPLC eluent (Fig. 4B) had absorption maxima at 467, 438 and 414 nm (Fig. 4C), similar to the peaks prior to purification (465, 435 and 409 nm) and corresponding to the absorption maxima of 4,4’-diaponeurosporene²³. The peak also had a mass-to-charge ratio (m/z) of 403.3338 (Fig. 4D), corresponding to the formula C₃₀H₄₂, which is identical to that of 4,4’-diaponeurosporene.

Subsequently, a high-throughput method was developed to phenotypically test C30 carotenoids, using less solvent and excluding unnecessary purification of the extract (see Methods). Starting from 575 in-house isolates and five publicly sourced strains (Supplementary Table S1) based on taxonomy and isolation source, 45 in-house and 5 public Lactobacillaceae were tested phenotypically. This selection included 17 related noncarriers to assess the absence of phenotype. In total, the biosynthesis of 4,4’-diaponeurosporene-like carotenoids was substantiated in 27 strains belonging to the crtMN-harboring species Lt. fragifolii, Lp. plantarum, Leuconostoc citreum, Lc. pseudomesenteroides and Holzapfelia floricola, and absent in all tested noncarrier strains (17) (Fig. 3, Supplementary Table S1). For only six of the isolates tested, the observed phenotype did not match the genotype. Specifically, the two tested strains of the species Lactiplantibacillus mudanjiangensis (AMBF-0197 and AMBF-0209) and one strain of Lc. mesenteroides (LMG 6893) did not appear to express the phenotype under the tested conditions (Fig. 1C), whereas two Lp. plantarum strains (AMBP-0214 and AMBP-0424) appeared to harbor inactive 4,4’-diapophytoene desaturase caused by the deletion of 522 bp in the crtN gene. These latter strains were considered natural nonproducing mutants of carotenoid production in Lp. plantarum, a characteristic employed in subsequent physiological tests. Finally, the Apilactobacillus xinyiensis(AMBP-0461) containing the unusual cluster with the duplicated crtN did not appear to express the phenotype under the tested conditions.

C30 carotenoids are associated with oxidative and UV stress resistance in Lp. plantarum strains

To assess whether carotenoid biosynthesis in Lactobacillaceae plays a role in resistance to oxidative and UV stress, we compared seven carotenoid-producing Lp. plantarum strains to the natural nonproducing strains AMBP-0214 and AMBP-0424. A phylogenetic tree of the core genomes of all the strains tested showed that all the strains were closely related and that the nonproducers were not an outgroup (Supplementary Figure S1). Carotenoid-negative strains were more susceptible to oxidative and UV-induced stress than carotenoid-positive strains, in line with the protective role of carotenoid biosynthesis in Lp. plantarum strains. The degree of oxidative and UV stress resistance varied among the strains tested (Fig. 5A). A linear model including time, different experiments and the production of carotenoids showed that the viability of the strains decreased with longer oxidative stress treatment (Est − 0.32, p < 0.001), with increased survival for the carotenoid-producing strains (indicated by a less negative slope) (Est. 0.06, p < 0.05). The same protection was observed for UV-induced stress, as phenotypic carotenoid biosynthesis in these strains was associated with increased survival in all trials (Est. 0.024, p < 0.05).

C30 carotenoid biosynthesis is associated with nomadic and insect-adapted lifestyles

We subsequently evaluated whether the C30 carotenoid biosynthesis genes were associated with a particular Lactobacillaceae lifestyle and habitat adaptation using the categories defined by Zheng et al.Click or tap here to enter text.as free-living, vertebrate-adapted, invertebrate-adapted or nomadic ³². Such lifestyle assignments were available for 170 out of 361 species in our public genome datasets. The presence of C30 carotenoid biosynthesis genes (crtMN) was most common in nomadic and insect-adapted lactobacilli, where they appeared as core (e.g. in nomadic Lp. plantarum, insect-adapted Fructilactobacillus lindneri and uncategorized Lc. citreum) and accessory traits. In contrast, crtMN genes were completely absent in vertebrate-associated lactobacilli such as Lactobacillus crispatus and Limosilactobacillus reuteri. They were also found in a few genomes representing free-living species such as Lv. buchneri (Fig. 5B). Consequently, a significant association between the proportion of genomes with C30 carotenoid biosynthesis and the lifestyle of a species (p < 0.001) was found when the Kruskal-Wallis test was used. Pairwise comparison with a Dunn test with Bonferroni correction for multiple testing (Supplementary Figure S2A) revealed the strongest association with nomadic and insect-adapted lifestyles. When considering the phylogenetic background with phylogenetic generalized least squares (PGLS), the same trends were found, but only insect-adapted versus free living was near significant (p = 0.059) (Supplementary Figure S2B).

Additionally, nomadic and insect-associated crtMN-harboring strains differed in genome size, with the nomadic strains having large genomes, similar to other nomadic Lactobacillaceae, and the insect-adapted strains having the smallest genomes within the family. These findings point to positive selection of these genes after horizontal gain in species with nomadic and insect-adapted lifestyles (Fig. 5C). To further investigate the lifestyle association, we also screened the isolation sources of 575 in-house isolates and three publicly sourced isolates. We correlated this with the C30 carotenoid phenotype and genotype. This analysis pointed towards leaves and flowers carrying a high number of carotenoid producers belonging to the species Lt. fragifolii (n = 13), Lc. pseudomesenteroides (n = 1), Lc. citreum (n = 3) and H. floricola (n = 1) (Supplementary Figure S3). A second major habitat of carotenoid producers showed to be plant-based fermentations, which harbored carotenoid-producing Lc. citreum (n = 2), Lp. plantarum (n = 6) and Lactiplantibacillus paraplantarum (n = 2). In contrast, in our laboratory, carotenoid-producing Lactobacillaceae strains were only sporadically isolated from vertebrate habitats (e.g., the human vagina⁴⁰ or respiratory tract⁴¹) and were all nomadic Lp. plantarum strains.

In this study, we investigated the biosynthesis of C30 carotenoids in the Lactobacillaceae family, the largest family of beneficial bacteria known to date, and linked the crtMN-positive genotype to their lifestyles. For this purpose, we used an integrated comparative genomic approach combined with a phenotypic screening of a diverse in-house Lactobacillaceae strain collection and an assessment of the functional role of carotenoids in stress resistance.

Pangenome analysis of Lactobacillaceae revealed that crtMN-mediated carotenoid biosynthesis was a rare and scattered trait within the family: it occurred in 41% of all the genera but only in a few species or strains within each genus. It was found to be a core property in some nomadic species, such as Lp. plantarum; insect-adapted species, such as Fl. lindneri; and accessory in many others. Phylogenetic analyses indicated that the biosynthetic crtMN genes are frequently transferred horizontally across species and genera, for example, from Lp. plantarum to Lv. brevis and within the Leuconostoc genus. This finding is consistent with the high mobility of carotenoid pathways observed at higher taxonomic levels^16–18. Moreover, the trait appeared to have been gained convergently in Fructilactobacillus, with this genus having acquired crtMN genes from distinct donors during distinct events. Using a high-throughput extraction and analysis method, the crtMN genotype was matched with the synthesis of 4,4’-diaponeurosporene in five species. Functionally, we showed that Lp. plantarum strains that do not produce carotenoids were less resistant to UV and oxidative stress, in line with the general knowledge on carotenoids¹⁶ and the previously observed protection against oxidative stress conferred by 4,4’-diaponeurosporene biosynthesis in Lp. plantarum⁴⁴. Finally, the scattered distribution and mobility of this trait across the family, coupled with its advantages in UV and oxidative stress, prompted us to systematically investigate the link between carotenoid biosynthesis and the lifestyle and ecology of Lactobacillaceae.

When testing for associations between two features associated with a set of species, as we have done here for carotenoids, phylogeny needs to be considered. This was done in our study by applying a phylogenetic generalized least squares (PGLS) approach⁴⁵. Such an approach effectively takes into account that closely related species will likely share similarity in any two traits (such as crtMN prevalence and lifestyle studied here) because of ”phylogenetic inertia”, not necessarily because the traits are correlated⁴⁶. Since the lifestyles of Lactobacillaceae species are mostly conserved at genus level, this implies that the independent units of information are in general the genera rather than the species, resulting in a lower effective sample size and limitations in our dataset concerning statistical power. In addition, another limitation of our work is that a lifestyle assessment has not yet be attributed by ^31,32 to 191 of the 361 Lactobacillaceae species in our dataset. This is due to lack of sufficient data on the isolation sources, metabolic potential, and related properties for these species^31,32. For many of these species, only a single strain has been isolated from a single source. Repetitive isolation of species from the same environment, as well as substantiation with specific metabolic and experimental validation is required to attribute lifestyles, as the environment of isolation does not reflect the environment of adaptation (niche) for various reasons, such as random dispersal events and increasing anthropogenic effects on the biosphere. Such detailed information will have to be collected for these 191 Lactobacillaceae species to also be able to attribute a lifestyle in the future and further substantiate our analyses.

Despite these shortcomings in public data and taking into account the phylogeny, a near significant association was found between carotenoid biosynthesis genes and an insect-adapted lifestyle (p = 0.056). However, excluding phylogenetic relatedness, we found that crtMN genes were mostly absent in free-living species and completely absent in vertebrate-associated species, such as the L. crispatus, which is dominant in the human vagina⁴⁰ and Lm. reuteri which typically colonizes the vertebrate gut³¹. The complete absence of carotenoids in well-studied vertebrate associated Lactobacillaceae suggests that carotenoid production is not selected for in mucosal and low-oxygen habitats, such as the gut and vagina. In contrast, carotenoid biosynthesis genes were strongly associated with nomadic and insect-adapted Lactobacillaceae indicating that oxygen- and UV-rich environments encountered by nomadic and insect-adapted strains can select for this trait. The habitat-adaptation strategy of nomadic carotenoid producers appeared to differ from that of insect-adapted species. Nomadic Lactobacillaceae species typically have large genomes (between 2.4 and 3.6 Mbp), making them metabolically versatile and adaptable to various environments. In other studies, nomadic species have been found to typically occur in low numbers in oligotrophic environmental niches, such as plant surfaces^29,52, and in high abundances once carbohydrates become more available, such as in vegetable fermentation products⁵³. Such fermentations are characterized by intense microbial competition and high-salt concentration⁵³, possibly leading to osmotic and oxidative stress. In such fermentations and outdoor oligotrophic environments, we speculate based on the data obtained here that C30 carotenoids could provide a fitness advantage to nomadic lactobacilli by reducing susceptibility to oxidative stress.

Insect-adapted lactobacilli have smaller genomes (1.2–2.2 Mbp) with a concomitant decrease in carbohydrate metabolic capacity ^31,32. This difference was also observed in our present study among crtMN carriers, as the insect-adapted crtMN carriers had remarkably small genome sizes (between 1.6 and 2.1 Mbp), but they still contained crtMN genes, indicating an evolutionary advantage. Interestingly, the insect-adapted crtMN carriers were all part of one clade within the Fructilactobacillus genus, a genus known to be transferred between pollinators via the environment, with flowers serving as key hubs^31,47,48. Our data presented here indicate that the biosynthesis of C30 carotenoids, and its associated protection against UV radiation and oxidative stress, could be a significant advantage for these environmentally dispersed, insect-adapted Lactobacillaceae species. Our hypothesis is in line with the adaptation strategy previously described for leaf-dwelling bacteria, such as Clavibacter and Pseudomonas, which produce C40 carotenoids and other pigments to increase their survival in this UV-stressed environment⁹. Notably, carotenoid biosynthesis was absent in Lactobacillaceae associated with social pollinators, such as Lactobacillus apis and Bombilactobacillus. These bacteria are vertically passed down to offspring within the hive⁴⁹, where they are protected from UV- and oxidative stress. An environmental survival strategy does not seem required for these bacteria⁴⁹. An exception to this is the Apilactobacillus genus, which is also known to be dispersed among solitary bees via flower⁵⁰. This genus was shown in our study here to have an unusual putative terpenoid cluster, characterized by a duplicated crtN and an absent crtM gene. However, it remains to be substantiated whether this duplication is associated with a particular phenotype or pigment.

Our hypothesis that Lactobacillaceae bacteria that are dispersed to plants via insects have a competitive advantage expressing C30 carotenoids was supported by our extensive culture approach and collection. A diverse array of carotenoid-producing Lactobacillaceae were isolated from flowers and leaves, with a relative high prevalence found for Lc. citreum in flowers. This species has – to the best of our knowledge - not yet been assigned to a certain lifestyle butour data presented here add support to an insect or flower-adapted lifestyle. In contrast, Lactobacillaceae isolates from vertebrate habitats studied here (mainly from the human vagina and respiratory tract), showed to be predominantly non-producers. Among the producing strains isolated, the nomadic Lp. plantarum was the most predominant species.

In addition to the ecological role, the presence of crtMN genes in Lactobacillaceae is of interest from an applied perspective, especially considering the beneficial properties and lack of virulence factors in this family of bacteria. For example, incorporating carotenoid-producing bacteria into food fermentations could add additional functional properties to these foods. In fact, 4,4’-diaponeurosporene is already present in many vegetable fermentations, as Lp. plantarum, a core producer, dominates the later stages of most typical vegetable fermentations, and Leuconostoc species generally dominate in the early stages⁵³. Although the added benefits of 4,4’-diaponeurosporene in these fermented food ecosystems have not yet been studied, this metabolite has been connected to health-promoting effects via immune modulation. For instance, the introduction of crtMN genes originating from Staphylococcus aureus into Bacillus subtilis has been shown to reduce colitis in mice⁵⁴ and increase resistance to Salmonella typhimurium infection⁵⁵. Furthermore, in piglets, heterologously crtMN expressing-B. subtilis bacteria have been shown to improve the mucosal immune system of the gut⁵⁶ and the respiratory tract⁵⁷. These studies were carried out with genetically modified bacteria and are thus unlikely to reach large market applications, especially in Europe. In contrast, our results presented here indicate that natural carotenoid-producing Lactobacillaceae constitute an interesting alternative. Moreover, microbial biosynthesis can offer advantages over traditional production methods at it can be safer and less reliant on fossil resources than chemical synthesis and less influenced by seasonality or climate than plant-based biosynthesis⁵⁸ and applied to C30 carotenoids⁵⁹.

In summary, this study on the ecology and evolution of carotenoid biosynthesis in the Lactobacillaceae family revealed a scattered distribution of crtMN-mediated C30 carotenoid 4,4’-diaponeurosporene biosynthesis across 28 species and 14 genera and highlighted the mobility of this trait. C30 carotenoid biosynthesis appears to have emerged as a core property in several species, notably Lp. plantarum, Lc. citreum and Fl. lindneri. Furthermore, carotenoid biosynthesis was strongly associated with nomadic and insect-adapted lifestyles, where it offers an advantage via protection from UV and oxidative stress.

Conflict of interest

S.L. received funding from different probiotic companies that were not involved in this research. T.E. is partially funded through an industrial research VLAIO grant not related to this work. M.L. is employed part-time by Biobest Group NV, this company was not involved in this research. I.V. and K.S. are employed by the RIC and their work was compensated via a service contract. P.A.B. is an independent consultant for several companies in the food and pharmaceutical industry bound by confidentiality agreements.

Funding

This research has received funding from the following funding bodies: the European Research Council (ERC; starting grant Lacto-Be 852600 of S.L.), the Belgian Science Policy Office (BELSPO; BRAIN-be 2.0; B@SEBALL (B2/191/P3/B@SEBALL)), the Scientific Research Foundation – Flanders (FWO; doctoral grant of J.T. 1SC3623N and postdoctoral grant of S.W. 12AZ624N), the Industrial Research Fund of the University of Antwerp (IOF-SEP 48530), VLAIO (HBC.2022.1000), and from the revenue of our microbiome platform services.

Acknowledgments

The authors would like to thank Tim Van Rillaer for his help with the bioinformatic pipelines and Ines Tuyaerts, Nele Van de Vliet and Sam Bakelants for their help with the UV and oxidative stress tests.

Sun, T., Rao, S., Zhou, X., Li, L.: Plant carotenoids: recent advances and future perspectives. Molecular Horticulture vol. 2 Preprint at (2022). https://doi.org/10.1186/s43897-022-00023-2
Havaux, M., And, Niyogi, K.K.: The Violaxanthin Cycle Protects Plants from Photooxidative Damage by More than One Mechanism. Plant Biology vol. 96 (1999).
Fraser, N.J., Hashimoto, H., Cogdell, R.J.: Carotenoids and Bacterial Photosynthesis: The Story so Far. Photosynth. Res. vol. 70 (2001)
Schwartz, S.H., Qin, X., Loewen, M.C.: The biochemical characterization of two carotenoid cleavage enzymes from Arabidopsis indicates that a carotenoid-derived compound inhibits lateral branching. J. Biol. Chem. 279, 46940–46945 (2004)
Booker, J., et al.: MAX3/CCD7 Is a Carotenoid Cleavage Dioxygenase Required for the Synthesis of a Novel Plant Signaling Molecule. Curr. Biol. 14, 1232–1238 (2004)
Schemske, D.W., Bradshaw, H.D.: Pollinator Preference and the Evolution of Floral Traits in Monkeyflowers (Mimulus).
Liu, G.Y., et al.: Staphylococcus aureus golden pigment impairs neutrophil killing and promotes virulence through its antioxidant activity. J. Exp. Med. 202, 209–215 (2005)
Mohanty, S.R., et al.: Methylotroph bacteria and cellular metabolite carotenoid alleviate ultraviolet radiation-driven abiotic stress in plants. Front. Microbiol. 13, 1–18 (2023)
Jacobs, J.L., Carroll, T.L., Sundin, G.W.: The role of pigmentation, ultraviolet radiation tolerance, and leaf colonization strategies in the epiphytic survival of phyllosphere bacteria. Microb. Ecol. 49, 104–113 (2005)
Moran, N.A., Jarvik, T.: Lateral transfer of genes from fungi underlies carotenoid production in aphids. Sci. (1979). 328, 624–627 (2010)
Karrer, P.: Carotenoids, Flavins and Vitamin A and B2. (1937). https://www.nobelprize.org/prizes/chemistry/1937/karrer/lecture/
Hajizadeh-Sharafabad, F., Zahabi, E.S., Malekahmadi, M., Zarrin, R., Alizadeh, M.: Carotenoids supplementation and inflammation: a systematic review and meta-analysis of randomized clinical trials. Critical Reviews in Food Science and Nutrition vol. 62 8161–8177 Preprint at (2022). https://doi.org/10.1080/10408398.2021.1925870
Dehnavi, M.K., Ebrahimpour-Koujan, S., Lotfi, K., Azadbakht, L.: The Association between Circulating Carotenoids and Risk of Breast Cancer: A Systematic Review and Dose–Response Meta-Analysis of Prospective Studies. Advances in Nutrition vol. 15 Preprint at (2024). https://doi.org/10.1016/j.advnut.2023.10.007
Zhang, Y., Yang, J., Na, X., Zhao, A.: Association between β-carotene supplementation and risk of cancer: a meta-analysis of randomized controlled trials. Nutr. Rev. 81, 1118–1130 (2023)
Yabuzaki, J.: Carotenoids Database: structures, chemical fingerprints and distribution among organisms. Database 1–11 (2017). (2017)
Klassen, J.L.: Phylogenetic and evolutionary patterns in microbial carotenoid biosynthesis are revealed by comparative genomics. PLoS One 5, (2010)
Armstrong,’, G.A., Hearstt, J.E.: Genetics and molecular biology of carotenoid pigment biosynthesis. FASEB J. 10, 228–237 (1996)
Santana-Molina, C., Henriques, V., Hornero-Méndez, D., Devos, D.P., Rivas-Marin, E.: The squalene route to C30 carotenoid biosynthesis and the origins of carotenoid biosynthetic pathways. Proc. Natl. Acad. Sci. U S A 119, (2022)
Santana-Molina, C., Henriques, V., Hornero-Méndez, D., Devos, D.P., Rivas-Marin, E.: Supplementary Materials for The Squalene Route to C30 Carotenoid Biosynthesis and the Origins of Carotenoid Biosynthetic Pathways
Augustin, J.M., Kuzina, V., Andersen, S.B., Bak, S.: Molecular activities, biosynthesis and evolution of triterpenoid saponins. Phytochemistry vol. 72 435–457 Preprint at (2011). https://doi.org/10.1016/j.phytochem.2011.01.015
Pollier, J., Goossens, A.: Oleanolic acid. Phytochemistry vol. 77 10–15 Preprint at (2012). https://doi.org/10.1016/j.phytochem.2011.12.022
Takaichi, S., et al.: The major carotenoid in all known species of heliobacteria is the C30 carotenoid 4,4’-diaponeurosporene, not neurosporene. Arch. Microbiol. 168, 277–281 (1997)
Garrido-Fernández, J., Maldonado-Barragán, A., Caballero-Guerrero, B., Hornero-Méndez, D., Ruiz-Barba, J.L.: Carotenoid production in Lactobacillus plantarum. Int. J. Food Microbiol. 140, 34–39 (2010)
Pramastya, H., Song, Y., Elfahmi, E.Y., Sukrasno, S., Quax, W.J.: Positioning Bacillus subtilis as terpenoid cell factory. Journal of Applied Microbiology vol. 130 1839–1856 Preprint at (2021). https://doi.org/10.1111/jam.14904
Taylor, R.F., Davies, B.H.: Triterpenoid carotenoids and related lipids. Triterpenoid carotenoid aldehydes from Streptococcus faecium UNH 564P. Biochem. J. 153, 233–239 (1976)
Pelz, A., et al.: Structure and biosynthesis of staphyloxanthin from Staphylococcus aureus. J. Biol. Chem. 280, 32493–32498 (2005)
Kleinig, H., Schmitt, R., Meister, W., Englert, G., Thommen, H.: New C30-Carotenoic Acid Glucosyl Esters from Pseudomonas rhodos. Z. fur Naturforschung - Sect. C J. Biosci. 34, 181–185 (1979)
Lebeer, S., Vanderleyden, J., De Keersmaecker, S.C.J.: Genes and Molecules of Lactobacilli Supporting Probiotic Action. Microbiol. Mol. Biol. Rev. 72, 728–764 (2008)
Yu, A.O., Leveau, J.H.J., Marco, M.L.: Abundance, diversity and plant-specific adaptations of plant-associated lactic acid bacteria. Environmental Microbiology Reports vol. 12 16–29 Preprint at (2020). https://doi.org/10.1111/1758-2229.12794
Turpin, W., et al.: PCR of crtNM combined with analytical biochemistry: An efficient way to identify carotenoid producing lactic acid bacteria. Syst. Appl. Microbiol. 39, 115–121 (2016)
Duar, R.M., et al.: Lifestyles in transition: evolution and natural history of the genus Lactobacillus. FEMS Microbiol. Rev. 41 (2017). https://doi.org/10.1093/FEMSRE/FUX030 S27–S48 Preprint at
Zheng, J., et al.: A taxonomic note on the genus Lactobacillus: Description of 23 novel genera, emended description of the genus Lactobacillus Beijerinck 1901, and union of Lactobacillaceae and Leuconostocaceae. Int. J. Syst. Evol. Microbiol. 70, 2782–2858 (2020)
Parks, D.H., Imelfort, M., Skennerton, C.T., Hugenholtz, P., Tyson, G.W.: CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055 (2015)
Wittouck, S., Wuyts, S., Meehan, C.J., van Noort, V., Lebeer, S.: A Genome-Based Species Taxonomy of the Lactobacillus Genus Complex. mSystems 4, (2019)
Blin, K., et al.: AntiSMASH 6.0: Improving cluster detection and comparison capabilities. Nucleic Acids Res. 49, W29–W35 (2021)
Legein, M., Wittouck, S., Lebeer, S.: Latilactobacillus fragifolii sp. nov., isolated from leaves of a strawberry plant (Fragaria x ananassa). Int. J. Syst. Evol. Microbiol. 72, 005193 (2022)
Eilers, T., Dillen, J., Van de Vliet, N., Wittouck, S., Lebeer, S.: Lactiplantibacillus carotarum AMBF275T sp. nov. isolated from carrot juice fermentation. Int. J. Syst. Evol. Microbiol. 73, (2023)
Fu, L., Niu, B., Zhu, Z., Wu, S., Li, W.: CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 28, 3150–3152 (2012)
Navarro-Muñoz, J.C., et al.: A computational framework to explore large-scale biosynthetic diversity. Nature Chemical Biology 2019 16:1 16, 60–68 (2019)
Lebeer, S., et al.: Citizen-science map of the vaginal microbiome. 10.21203/rs.3.rs-1350465/v1
De Boeck, I., et al.: Lactobacilli Have a Niche in the Human Nose. Cell. Rep. 31, (2020)
Temmermans, J., et al.: The biocontrol agent Lactiplantibacillus plantarum AMBP214 is dispersible to plants via bumblebees. Appl. Environ. Microbiol. (2023). 10.1128/aem.00950-23
Legein, M., et al.: Compost Teas: market survey and microbial analysis 1. Biorxiv (2022). 10.1101/2022.08.05.503013
Kim, M., Jung, D., Seo, D., Park, Y., Seo, M.: 4,4′-Diaponeurosporene from Lactobacillus plantarum subsp. plantarum KCCP11226: Low Temperature Stress- Induced Production Enhancement and In Vitro Antioxidant Activity. J. Microbiol. Biotechnol. 31, 63–69 (2021)
Grafen, A.: The phylogenetic regression. Philosophical Trans. Royal Soc. Lond. B Biol. Sci. 326, 119–157 (1989)
Felsenstein, J.: PHYLOGENIES AND THE COMPARATIVE METHOD. Am. Nat. 125, 1–15 (1985)
McFrederick, Q.S., et al.: Flowers and Wild Megachilid Bees Share Microbes. Microb. Ecol. 73, 188–200 (2017)
McFrederick, Q.S., et al.: Environment or kin: Whence do bees obtain acidophilic bacteria? Mol. Ecol. 21, 1754–1768 (2012)
Kwong, W.K., Moran, N.A.: Gut microbial communities of social bees. Nature Reviews Microbiology vol. 14 374–384 Preprint at (2016). https://doi.org/10.1038/nrmicro.2016.43
Vuong, H.Q., Mcfrederick, Q.S., Angert, E.: Comparative Genomics of Wild Bee and Flower Isolated Lactobacillus Reveals Potential Adaptation to the Bee Host. Genome Biol. Evol. 11, 2151–2161 (2019)
Martino, M.E., et al.: Nomadic lifestyle of Lactobacillus plantarum revealed by comparative genomics of 54 strains isolated from different habitats. Environ. Microbiol. 18, 4974–4989 (2016)
Miller, E.R., et al.: Establishment limitation constrains the abundance of lactic acid bacteria in the Napa cabbage phyllosphere. Appl. Environ. Microbiol. 85, (2019)
Wuyts, S., et al.: Carrot juice fermentations as man-made microbial ecosystems dominated by lactic acid bacteria. Appl. Environ. Microbiol. 84, 2021 (2018)
Jing, Y., Liu, H., Xu, W., Yang, Q.: Amelioration of the DSS-induced colitis in mice by pretreatment with 4,4’-diaponeurosporene-producing Bacillus subtilis. Exp. Ther. Med. 14, 6069–6073 (2017)
Liu, H., Xu, W., Yu, Q., Yang, Q.: 4,4’-diaponeurosporene-producing Bacillus subtilis increased mouse resistance against Salmonella typhimurium infection in a CD36-dependent manner. Front. Immunol. 8, 483 (2017)
Jing, Y., Liu, H., Xu, W., Yang, Q.: 4,4′-Diaponeurosporene-Producing Bacillus subtilis Promotes the Development of the Mucosal Immune System of the Piglet Gut. Anat. Rec. 302, 1800–1807 (2019)
Zhang, P., Huang, L., Zhang, E., Yuan, C., Yang, Q.: Oral administration of Bacillus subtilis promotes homing of CD3 + T cells and IgA-secreting cells to the respiratory tract in piglets. Res. Vet. Sci. 136, 310–317 (2021)
Cravens, A., Payne, J., Smolke, C.D.: Synthetic biology strategies for microbial biosynthesis of plant natural products. Nature Communications vol. 10 Preprint at (2019). https://doi.org/10.1038/s41467-019-09848-w
Siziya, I.N., Hwang, C.Y., Seo, M.J.: Antioxidant Potential and Capacity of Microorganism-Sourced C30 Carotenoids—A Review. Antioxidants 11, (2022)

Yes there is potential Competing Interest. S.L. received funding from different probiotic companies that were not involved in this research. T.E. is partially funded through an industrial research VLAIO grant not related to this work. M.L. is employed part-time by Biobest Group NV, this company was not involved in this research. I.V. and K.S. are employed by the RIC and their work was compensated via a service contract. P.A.B. is an independent consultant for several companies in the food and pharmaceutical industry bound by confidentiality agreements.

Download PDF

Version 1

posted

You are reading this latest preprint version

Distribution of C30 carotenoid biosynthesis genes suggests habitat adaptation function in insect-adapted and nomadic Lactobacillaceae

Status:

Version 1

Abstract

Figures

Introduction

Materials and methods

Strains used in this study

C30 carotenoid extraction and identification

Oxidative and UV stress resistance assays

Association of C30 carotenoid biosynthesis capacity and habitat

Data availability

Results

C30 carotenoid biosynthesis is associated with nomadic and insect-adapted lifestyles

Discussion

Declarations

Conflict of interest

Funding

Acknowledgments

References

Additional Declarations

Supplementary Files

Status:

Version 1