Next Article in Journal
The Impact of Dietary Knowledge on Health: Evidence from the China Health and Nutrition Survey
Next Article in Special Issue
Enhancing Nature Connection and Positive Affect in Children through Mindful Engagement with Natural Environments
Previous Article in Journal
Electrophysiological Brain-Cardiac Coupling in Train Drivers during Monotonous Driving
Previous Article in Special Issue
Enhancing Adolescent Girls’ Well-Being in the Arctic—Finding What Motivates Spending Time in Nature
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 18
Issue 7
10.3390/ijerph18073742
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Do Rural Second Homes Shape Commensal Microbiota of Urban Dwellers? A Pilot Study among Urban Elderly in Finland

by
Mika Saarenpää
1,
Marja I. Roslund
1,
Riikka Puhakka
1,
Mira Grönroos
1,
Anirudra Parajuli
1,2,
Nan Hui
1,3,
Noora Nurminen
4,
Olli H. Laitinen
4,
Heikki Hyöty
4,
Ondrej Cinek
5,
Aki Sinkkonen
1,6,* and
the ADELE Research Group
†
1
Ecosystems and Environment Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Niemenkatu 73, 15140 Lahti, Finland
2
Department of Medicine, Karolinska University Hospital, Huddinge, 141 86 Stockholm, Sweden
3
School of Agriculture and Biology, Shanghai Jiao Tong University, 800 Dongchuan RD. Minhang District, Shanghai 200240, China
4
Faculty of Medicine and Health Technology, Tampere University, Arvo Ylpön katu 34, 33520 Tampere, Finland
5
Second Faculty of Medicine, Charles University, V Úvalu 84, 150 06 Prague 5, Czech Republic
6
Natural Resources Institute Finland, Itäinen Pitkäkatu 4 A, 20520 Turku, Finland
*
Author to whom correspondence should be addressed.
†
The ADELE research group: Damiano Cerrone, Mira Grönroos, Nan Hui, Heikki Hyöty, Olli H. Laitinen, Iida Mäkelä, Noora Nurminen, Sami Oikarinen, Anirudra Parajuli, Riikka Puhakka, Juho Rajaniemi, Marja I. Roslund, Mika Saarenpää, Aki Sinkkonen, Laura Soininen, Yan Sun, Raisa Valve, Heli K. Vari.
Int. J. Environ. Res. Public Health 2021, 18(7), 3742; https://doi.org/10.3390/ijerph18073742
Submission received: 4 March 2021 / Revised: 29 March 2021 / Accepted: 31 March 2021 / Published: 2 April 2021
(This article belongs to the Special Issue Well-Being and Health Impacts of Outdoor Recreation in Natural Settings)
Browse Figures
Review Reports Versions Notes

Abstract

:
According to the hygiene and biodiversity hypotheses, increased hygiene levels and reduced contact with biodiversity can partially explain the high prevalence of immune-mediated diseases in developed countries. A disturbed commensal microbiota, especially in the gut, has been linked to multiple immune-mediated diseases. Previous studies imply that gut microbiota composition is associated with the everyday living environment and can be modified by increasing direct physical exposure to biodiverse materials. In this pilot study, the effects of rural-second-home tourism were investigated on the gut microbiota for the first time. Rural-second-home tourism, a popular form of outdoor recreation in Northern Europe, North America, and Russia, has the potential to alter the human microbiota by increasing exposure to nature and environmental microbes. The hypotheses were that the use of rural second homes is associated with differences in the gut microbiota and that the microbiota related to health benefits are more diverse or common among the rural-second-home users. Based on 16S rRNA Illumina MiSeq sequencing of stool samples from 10 urban elderly having access and 15 lacking access to a rural second home, the first hypothesis was supported: the use of rural second homes was found to be associated with lower gut microbiota diversity and RIG-I-like receptor signaling pathway levels. The second hypothesis was not supported: health-related microbiota were not more diverse or common among the second-home users. The current study encourages further research on the possible health outcomes or causes of the observed microbiological differences. Activities and diet during second-home visits, standard of equipment, surrounding environment, and length of the visits are all postulated to play a role in determining the effects of rural-second-home tourism on the gut microbiota.

1. Introduction

The prevalence of immune-mediated diseases, such as allergies, asthma, and autoimmune diseases, have increased significantly in developed countries over the last 30 years [1]. For example, the prevalence of asthma in developed countries can be as high as 21% [2]. The amount of urban developed area is expected to triple in size by 2030 [3], and in 2050, 70% of the world population is predicted to be living in cities [4]. With rapid urbanization, the burden of immune-mediated diseases is expected to increase abruptly in the 21st century.
Urban lifestyle is often characterized by a high level of hygiene, and according to the hygiene hypothesis and its variants (e.g., biodiversity hypothesis), decreased interaction with microbes can explain the high prevalence of immune-mediated diseases in developed countries [5,6,7,8]. Exposure to diverse microbial communities, especially in childhood, is essential for the natural development and functioning of the human immune system [6,8,9,10,11,12]. Kondrashova et al. [13] found the incidence of several immune-mediated diseases in the more urbanized Eastern Finland to be several times higher than in the adjacent Russian Karelia. Populations in both areas have a similar genetic background, and the higher incidence on the Finnish side is therefore believed to be caused by environmental factors, such as high level of hygiene and lack of daily contacts with environmental microbiota.
In addition to the lifestyle, urbanization causes changes in the environment, and urban areas tend to be biologically less diverse than the surrounding rural areas [14]. Pollution, for example, has been found to affect soil microbiota, and soil contamination has been associated with altered commensal microbiota and endocrine disruption [15,16,17]. The biodiversity hypothesis states that contact with natural environments enriches the human microbiota and thus can help prevent the development of immune-mediated diseases [7,18,19,20]. Interaction with agricultural land and forests is considered particularly important [21]. Children living on farms have been found to suffer less from asthma and atopy [10,12,21], and lower biodiversity and amount of green environment around home have been associated with atopic sensitization and altered gut microbiota composition [7,22,23]. Urbanization has decreased contact with nature and diverse environmental microbial communities [24,25]. This can lead to an imbalance in the human microbiota (known as dysbiosis) and eventually to immunological disorders [6,7,10,20].
Gut microbes are associated with multiple systemic and immunomodulatory effects and endocrine signaling pathways, and gut microbiota dysbiosis has been linked with inflammatory bowel disease (IBD), colorectal cancer, obesity, allergy, atopy, and type 1 diabetes, among others [26,27,28,29,30,31]. Different types of gut microbiota dysbiosis are characterized by distinct changes in microbiota composition. For example, decreased total bacterial diversity has been observed in IBD, psoriatic arthritis, and Sjögren’s syndrome [32,33,34], reduced complexity of the phylum Firmicutes in Crohn’s disease [35], and increased ratio of Firmicutes to Bacteroidetes in irritable bowel syndrome [36]. Changes in the butyrate-producing bacteria belonging to, for example, the families Lachnospiraceae and Ruminococcaceae are of special importance in overall gut health and in the development of numerous inflammatory conditions [37,38,39]. The elderly gut microbiota profile differs from the healthy adult one, and studies have reported a decrease in the overall diversity and species diversity for Bacteroides, Prevotella, Bifidobacterium, and Lactobacillus along with an increase in the species number within the Enterobacteriaceae, Clostridium, Proteobacteria, staphylococci, and streptococci [40,41,42,43].
The interaction with nature and environmental materials, for example through outdoor recreation or landscaping materials, offers a natural way of increasing exposure to diverse microbial communities. In recent studies, a diverse environmental microbiota has been incorporated into gardening and health-promoting materials to increase the microbial exposure of urban dwellers. These studies demonstrate for the first time how even short-term direct contacts with soil and plant materials affect the skin and gut microbiota and enhance immune modulation [44,45,46,47].
While deliberate direct exposure to environmental materials and microbes has been found to modify the human microbiota and immune response, there is little research on the possible effects of unimposed nature-based activities, such as outdoor recreation. Outdoor recreation refers to any leisure time activity taking place in a natural environment, such as hiking, fishing, hunting, berry and mushroom picking, camping, swimming, leisure-time gardening, bird watching, and viewing landscapes. Many of these activities increase the physical interaction with natural environments and environmental microbes, which might result in changes in the human microbiota and immune response [23,24,44,45,48].
In Finland, 96% of the adult population reports participating in outdoor activities, on average 2–3 times a week and 170 times a year [49]. Rural-second-home tourism is a popular form of outdoor recreation in Finland and other parts of Northern Europe, as well as in Russia and North America. It provides multiple psychological, social, and physiological health and well-being benefits for all age groups [50]. In Finland, over 40% of the population has been estimated to visit a rural second home regularly, and the amount of rural second homes per capita is one of the highest in the world [51,52]. As a consequence of coronavirus disease 2019 (COVID-19), second homes sales and rentals have increased in Finland and around the world, as people wish to escape the pandemic to the safety of rural second homes [53].
In Finland, visiting rural second homes has been found to increase the rate of participation in many forms of outdoor recreation, such as wild berry and mushroom picking, fishing, swimming, growing own vegetables, chopping firewood, and going for a walk [49,54]. Most Finnish rural second homes are still without an indoor flush toilet, dishwasher, shower, and other amenities, and water for bathing, dishwashing, and drinking often comes from the natural waters (e.g., well, spring, lake) [55]. Rural second home visits characterized by various outdoor activities and a lower level of hygiene have the potential to alter the human microbiota, especially among urban people who have been found to use rural second homes more often than their rural counterparts [51]. The summer vacation season is the main rural-second-home season in Finland, but the standard of equipment is increasing and allows for prolonged second-home seasons. Currently, a third of the Finnish rural second homes are suitable for winter use [55]. Visiting rural second homes during cold, snowy winters might not have the same effect on microbial exposure as in summer, as the transfer of health-related environmental bacteria indoors has been shown to be lower in winter [56].
In this study, the fecal bacterial communities (as a proxy for gut microbiota) were compared between urban elderly visiting rural second homes and elderly without regular access to a rural second home (hereafter referred to as the control group). It was tested whether the use of rural second homes is associated with changes in fecal bacterial composition. It has been shown that the everyday living environment affects commensal microbiota and that the microbiota composition differs between people living in urban and in rural environments [23,24,57]. It is not known whether visiting rural second homes is associated with gut microbiota composition. In the current study, it was hypothesized that the rural-second-home users have differences in their gut microbiota compared to the control group. It was also hypothesized that the gut microbiota associated with health benefits are more diverse or more common among rural-second-home users than in the control group.

2. Materials and Methods

2.1. Study Area and Participants

The participants—25 elderly people (65–79 years) living in urban apartment houses in the city of Lahti, Finland—were chosen from a large prospective study called GOAL (Good Aging in Lahti region, 2002–2012). The original GOAL study aimed to find connections between the living environment and chronic diseases and functional disabilities in the elderly and retired population [58]. In 2015, 60 GOAL participants were invited to participate in a follow-up study (2015–2016). In the follow-up study, stool samples for microbial analysis and survey data were collected [23,25]. Twenty-five participants from the follow-up study living in urban areas were included in this pilot study, and 10 of them were rural-second-home users. The second-home users were asked to report the number of years they had been using a second home, the number of visiting days per year, and how much time was usually spent outdoors vs. indoors (Table 1).
All participants with at least one of the following noncommunicable diseases (NCDs, also known as chronic diseases) affecting the immune response were excluded in the 2015 follow-up study: celiac disease, rheumatoid arthritis, chronic obstructive pulmonary disease, diabetes, psoriasis, multiple sclerosis, corticosteroid-treated asthma, chronic obstructive pulmonary disease, and actively treated or metastatic cancer. Daily smokers and participants with dementia or treated with corticosteroid and other immunosuppressive medication were also excluded. If the participants were treated with antibiotics within the last six months before giving the stool samples, they were excluded. None of the participants owned pets with access to the outdoors.

2.2. Sample Collection, DNA Extraction, Amplification, and Sequencing

In this study, stool samples collected in August 2015 were analyzed (see [23], Supplementary Figure S2 for a flowchart describing the sample collection and subsequent processes). Samples collected in August were used because the summer vacation season is the peak rural-second-home season in Finland and because exposure to environmental microbiota among aging Finns is high during the summer months [56]. Sample collection, DNA extraction, and amplicon sequencing were described in our previous study [44]. In short, participants took the samples independently using a collection kit and stored them at −20 °C until the samples were collected by the study personnel. The samples were transferred in dry ice and stored at −80 °C until further analyses. DNA was extracted from 30 to 60 mg of frozen, unprocessed stool, and bacterial community analyses were based on the amplification of the V4 region of 16S rRNA. Raw sequence reads are available in the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra accessed on 31 March 2021) with accession numbers SAMN08991885–SAMN08992045.

2.3. Bioinformatics

Paired-end sequence data (.fastq) from the stool samples were processed using Mothur (version 1.39.5, http://www.mothur.org accessed on 31 March 2021) [59], as described in our earlier studies [23,44], following the protocols by Schloss and Westcott [60] and Kozich et al. [61]. The sequences were aligned using the Mothur version of SILVA bacterial reference (version 132) [62]. Less abundant (≤10 sequences across all experimental units) operational taxonomic units (OTUs) were removed to avoid PCR or sequencing artifacts [63]. Each sample was subsampled to 4024 sequences to control for the varying number of sequences.
The 16S rRNA OTU data were picked against the Greengenes Database [64] according to 97% similarity, and the OTU table was used to predict metabolic functions by referencing the Kyoto Encyclopedia of Genes and Genomes (KEGG) ortholog database (release 89.0) in PICRUSt (http://picrust.github.io/picrust/ accessed on 31 March 2021) [65].

2.4. Statistical Analyses

Statistical analyses were conducted for the whole bacterial community at different taxonomic levels (i.e., OTU, family, order, class, and phylum) and for the functional orthologs. At the OTU level, analyses were conducted within the most abundant dominant gut phyla and predominant colon families (Table 2); p-values were corrected with Benjamini–Hochberg correction to account for the multiple testing [66].
The Shannon and Simpson diversity indices, observed richness, and predicted functional orthologs of the bacterial communities were compared between the rural-second-home users and the control group using the Student’s t-test or Mann–Whitney U test. The Student’s t-test was used when the data were normally distributed based on the Shapiro–Wilk test, and the Mann–Whitney U test was used when the data were not normally distributed.
Differences in the bacterial community composition between the groups were compared using permutational multivariate analysis of variance (PERMANOVA; 999 permutations). The bacterial community compositions were visualized with non-metric multidimensional scaling (NMDS). Both PERMANOVA and NMDS were based on the Bray–Curtis dissimilarity. All statistical analyses and data visualization were performed using the R statistical software environment (version 4.0.3, R Foundation, Vienna, Austria) [68] and the packages phyloseq (version 1.32.0) [69], ggplot2 (version 3.3.2) [70], and vegan (version 2.5-6) [71].

3. Results

3.1. Characterization of the Bacterial Communities

Based on Illumina MiSeq sequencing of the bacterial 16S rRNA gene, the most common bacterial phyla in the stool samples of the rural-second-home users and the control group were Firmicutes, Bacteroidetes, Proteobacteria, and Actinobacteria (Figure 1).

3.2. Diversity of the Bacterial Communities

The Shannon diversity of the fecal microbiota was higher in the control group than among the rural-second-home users at the OTU, family, order, class, and phylum levels, for the phylum Firmicutes and for the families Ruminococcaceae and Lachnospiraceae (Figure 2, Table 3). For the phyla Bacteroidetes, Actinobacteria, and Proteobacteria and for the family Bacteroidaceae, no differences in Shannon diversity were observed between the rural-second-home users and the control group. The observed richness followed a similar pattern and was higher in the control group than among the rural-second-home users at the OTU, family, order, class, and phylum levels, for the phylum Firmicutes and for the family Ruminococcaceae. For the phyla Bacteroidetes, Actinobacteria, and Proteobacteria and for the families Lachnospiraceae and Bacteroidaceae, no differences in the observed richness were detected between the groups. The Simpson diversity differed from the Shannon diversity and observed richness: only the diversity of the phylum Firmicutes was observed to be higher in the control group than among the rural-second-home users, and at all other taxonomic levels, no differences were detected between the groups.

3.3. Bacterial Community Composition and Functional Profile

The permutational analysis of variance (PERMANOVA) revealed that the fecal bacterial community compositions did not differ between the rural-second-home users and the control group at the OTU (Figure 3) or any other taxonomic level (Table 4).
According to the predicted gut metagenomic functions, the fecal microbiota of the rural-second-home users showed lower levels of RIG-I-like receptor signaling pathway than the control group (p = 0.019). No other differences in the functional orthologs were observed (p > 0.1, except for endocytosis (p > 0.006); data not shown).

4. Discussion

In this pilot study, the fecal bacterial communities of urban elderly with and without regular access to a rural second home were compared to see whether the use of second homes is associated with changes in the gut microbiota. Previous studies have shown how the everyday living environment affects the human microbiota [7,22,24,25], but this is the first time when the effects of rural-second-home tourism, a form of outdoor recreation, on the gut microbiota have been investigated. The hypotheses were that the rural-second-home users have differences in their gut microbiota compared to the control group and that the gut microbiota associated with health benefits are more diverse or more common among the second-home users. Based on Shannon diversity and observed richness, this pilot study indicates that the fecal bacterial diversity of urban elderly rural-second-home users is lower than in the control group. Differences in the fecal bacterial diversity were observed at the OTU, family, order, class, and phylum levels, and for the OTUs, within the phylum Firmicutes and the families Lachnospiraceae and Ruminococcaceae. Interestingly, when the Simpson diversity indices were compared between the groups, differences in fecal bacterial diversity were observed only for the OTUs within the phylum Firmicutes. The Simpson diversity has been seen to emphasize the evenness of the community, whereas the Shannon diversity is more affected by the amount of rare species. The results were somewhat surprising, as previous studies have shown that direct contact with nature diversifies the gut microbiota [44,47]. The researchers originally expected that rural-second-home users regularly engage in activities that increase exposure to nature. A diverse gut microbiota is often associated with health benefits and eubiosis. While the observed differences in α-diversity confirm that the gut microbiota of the rural-second-home users is different compared to that of the control group, the microbiota associated with health benefits (Lachnospiraceae and Ruminococcaceae) were not more diverse or common among the second-home users. Furthermore, the bacterial community composition did not differ between the rural-second-home users and the control group.
Previous studies have found that the gut microbiota of elderly people with health issues is characterized by, for example, a reduction in the overall diversity, Bacteroides-Prevotella group, bifidobacteria, lactobacilli, and some clostridia, and an increase in enterobacteria, proteobacteria, staphylococci, and streptococci [72,73,74]. A reduction in the phylum Firmicutes, a major gut phylum, has also been observed in the elderly [75]. The reduction of Firmicutes has also been linked to colorectal cancer [76], type 2 diabetes [77], and Crohn’s disease [35]. On the other hand, elevated Firmicutes levels, especially of bacteria belonging to the family Lachnospiraceae, have been associated with irritable bowel syndrome (IBS) [78]. The families Ruminococcaceae and Lachnospiraceae, both belonging to the phylum Firmicutes, contain many species that produce butyrate, an important bacterial metabolite with positive anti-inflammatory, epithelial barrier-protective, and cell-regulatory effects [38]. The decreased abundance of both Ruminococcaceae and Lachnospiraceae has been linked with microbial dysbiosis in IBD [27,76,79] and even with a higher risk of stroke [80]. A reduction in the overall diversity of the gut microbiota has also been observed in IBD [32], as well as in psoriatic arthritis and Sjögren syndrome [33,34].
Based on the predicted gut metagenomic functions, the RIG-I-like receptor signaling pathway was lower among the rural-second-home users than in the control group. RIG-I-like receptors are part of the eukaryotic intracellular pattern recognition, and bacterial and viral infections are known to activate the RIG signaling pathway [81,82]. The results do not suggest that the detected bacteria express RIG-I-like receptor signaling pathway genes but that they might express molecules that could act as part of the pathway. The effects of RIG-I-like receptors are mediated by the induction of interferons [83]. As RIG-I-like receptors are able to induce a strong proinflammatory response through interferon-stimulated elements, this pathway is tightly regulated [84,85] to avoid inflammation-associated overreactions. Intriguingly, RIG-I-mediated interferon responses have been detected without viral or bacterial triggers in association with autoimmune diseases [86], such as relapsing–remitting multiple sclerosis, Rheumatoid arthritis, Sjögren’s syndrome, and systemic lupus erythematosus [87,88,89,90,91]. The current results suggest that the RIG-I-like receptors in elderly spending time at rural second homes might be better regulated, leading to a lower risk of immune-mediated diseases, such as allergies [22].
Since specific changes in the gut microbiota can have both health- and illness-related consequences, it is difficult to determine whether the observed differences are associated with health outcomes. This is even more true in the case of the elderly, whose gut microbiota display greater inter-individual variability than that of younger adults [92]. This variability is caused by many factors, such as changes in digestion, bowel function, medication, everyday living environment, mobility, and, perhaps most importantly, diet [41,42,43,92]. Patterns in the microbiota composition also differ between geographical regions, which highlights the danger of generalizing health-related observations in the elderly gut microbiota [42,43,93]. Because defining a healthy elderly gut microbiota is challenging, it cannot be concluded that the lower gut microbiota diversity of rural-second-home users is a sign of gut dysbiosis, especially since the community composition did not differ between the studied groups. Rural-second-home tourism comes in many forms, and it can be postulated that there is a myriad of factors that determine the effects of rural-second-home tourism on the gut microbiota.
In a previous study with the same study population and stool samples, the effects of ethnicity, socioeconomic status, high parasite load, access to modern health care, and use of antibiotics were ruled out, and the importance of the everyday interaction with natural environments in the homeostasis of elderly gut microbiota was underlined (see [23], Supplementary Table S1). The type and quality of the environment surrounding homes are important determinants of the commensal microbiota [22,24,48], and for example, broad-leaved forests and diverse yard vegetation have health-related effects on the gut microbiota [17,23]. Unfortunately, information on the surrounding environment of the second homes was not collected in this study. Therefore, the details of the non-microbial diversity surrounding the second homes are unknown. However, the coverage of built area surrounding the rural second homes can be assumed to be low, i.e., resembling the land cover structure in sparsely populated rural areas [23,25]. Information on the activities during the second home visits was also missing. However, most second-home users did report that they tend to spend more time outdoors while visiting the second homes (Table 1). Having direct physical contact with nature (e.g., chopping firewood, picking wild berries or mushrooms, swimming, or gardening) probably has a different effect on the commensal microbiota than sitting on a patio or going for a walk.
While no information on the standard of equipment of the second homes is available, the equipment plausibly varies based on the reported number of visiting days. On average, rural-second-home users reported spending 67 days a year at their second homes (Table 1), but the timing and exact number of visiting days in 2015 is unknown. Those rural second homes where study participants reported spending 100–150 days per year are probably fit for winter use and resemble ordinary detached homes. However, most rural second homes in Finland, and presumably also in the current study, are typical lakeside summer cottages with latrines and no running water [51,55]. The standard of equipment can be assumed to affect the exposure to microbes through differing activities and levels of hygiene. For example, if the second home lacks electricity and electrical appliances, many activities and chores might take place outdoors (e.g., open-fire grilling). The lack of running tap water often means that water for bathing, washing dishes, and drinking comes from natural waters. Natural water sources often contain more organic material than chlorinated tap water, which can affect human exposure to microbes [94]. As washing dishes by hand is associated with a low incidence of allergic diseases [95], contact with natural water sources can be expected to have a similar health-related effect.
The effect of diet on the elderly gut microbiota composition was ruled out in a previous study with the same study population [23]. Therefore, information on the dietary habits during the second-home visits was not collected in this study. In the previous study [23], diet-related questions targeted only the everyday dietary habits, and these habits might be different during rural-second-home visits due to possibly more rudimentary conditions. While some studies have observed that rural-second-home users favor homegrown and local products [96], heavily grilled food (e.g., sausages) and alcohol are also an important component of the rural-second-home life for many [50]. The elderly gut microbiota has been strongly linked with diet [92], and diets rich in fat and poor in dietary fiber along with overconsumption of alcohol have been associated with disadvantageous changes in the gut microbiota [97,98,99,100,101,102]. Interestingly, butyrate-producing Firmicutes, such as Lachnospiraceae, are sensitive to changes in dietary fiber [103]. Dietary choices at rural second homes might be imbalanced and non-diversified due to, for example, long distances to grocery stores.
While lower gut microbiota diversity and RIG-I-like receptor signaling pathway levels were observed in urban elderly rural-second-home users, the small sample size limits the possibilities to draw conclusions on the nature or causes of these differences. Arguably, rural-second-home tourism is not automatically associated with only health- or illness-related changes in the gut microbiota. Instead, individual choices concerning diet and activities are potentially crucial explanatory factors. There are as many ways of using rural second homes as there are users. Therefore, larger studies with different age groups and geographical regions are needed. It would also be advisable to compare the fecal bacterial communities before and after the rural-second-home visits. To overcome the uncertainties of the current findings, future research should take into account the standard of equipment of the rural second homes, land cover and biodiversity around the homes, and dietary habits during the second-home visits.

5. Conclusions

It was hypothesized that there are differences in the gut microbiota between rural-second-home users and control group and that microbiota associated with health benefits are more diverse or common among rural-second-home users. The findings of this pilot study suggest that there are differences in the gut microbiota: the use of rural second homes is associated with lower gut microbiota diversity and RIG-I-like receptor signaling pathway levels when compared with non-use. However, the microbiota associated with health benefits were not more diverse or common among the rural-second-home users. Although the observed differences have been associated with health outcomes in previous studies, making any conclusions related to these is not meaningful in the context of this study. Larger studies taking into account the standard of equipment of the rural second homes, the type and quality of the surrounding environment, the choice of activities, the length and season of the visits, and dietary habits are encouraged to understand the role of rural-second-home tourism in shaping the gut microbiota structure and function and their potential association with health.

Author Contributions

Conceptualization, R.P., M.S., M.I.R., M.G. and A.S.; methodology, R.P., M.S., M.I.R., M.G., A.S., A.P. and N.H.; validation, R.P., M.I.R. and A.S.; formal analysis, M.S., M.I.R., M.G. and A.P.; investigation, R.P., M.S., M.I.R., M.G., A.P., N.H. and N.N.; resources, A.S., O.C. and H.H.; data curation, R.P., M.S., M.I.R., M.G., A.P., N.H., N.N. and O.C.; writing—original draft preparation, M.S.; writing—review and editing, R.P., M.S., M.I.R., A.S., N.H., O.H.L. and O.C.; visualization, M.S.; supervision, R.P., M.I.R., M.G. and A.S.; project administration, A.S.; funding acquisition, M.S. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Finnish Funding Agency for Innovation (TEKES, currently known as Business Finland), grant number 40333/14; by Child and Nature Foundation (Lapsi ja Luonto Säätiö), personal grant to Mika Saarenpää; and by Onni and Hilja Tuovinen Foundation (Onni ja Hilja Tuovisen Säätiö), personal grant to Mika Saarenpää.

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Tampere University Hospital (approval code R15044).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Raw sequence reads are available in the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra, accessed on 31 March 2021) with accession numbers SAMN08991885-SAMN08992045.

Acknowledgments

We would like to thank all the ADELE team members for their support, volunteers for participating in the study, CSC—IT Center for Science, Finland, for computational resources, and personnel at the Environmental Laboratory at the University of Helsinki for laboratory services. BEI Resources (American Type Culture Collection, ATCC) are also gratefully acknowledged for the Genomic DNA from Microbial Mock Community (HM-276D).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lerner, A.; Jeremias, P.; Matthias, T. The World Incidence and Prevalence of Autoimmune Diseases Is Increasing. IJCD 2016, 3, 151–155. [Google Scholar] [CrossRef] [Green Version]
  2. To, T.; Stanojevic, S.; Moores, G.; Gershon, A.S.; Bateman, E.D.; Cruz, A.A.; Boulet, L.-P. Global Asthma Prevalence in Adults: Findings from the Cross-Sectional World Health Survey. BMC Public Health 2012, 12, 204. [Google Scholar] [CrossRef] [Green Version]
  3. Seto, K.C.; Güneralp, B.; Hutyra, L.R. Global Forecasts of Urban Expansion to 2030 and Direct Impacts on Biodiversity and Carbon Pools. Proc. Natl. Acad. Sci. USA 2012, 109, 16083–16088. [Google Scholar] [CrossRef] [Green Version]
  4. United Nations, Department of Economic and Social Affairs, Population Division. World Urbanization Prospects: The 2018 Revision (ST/ESA/SER.A/420); United Nations: New York, NY, USA, 2019. [Google Scholar]
  5. Strachan, D.P. Hay Fever, Hygiene, and Household Size. BMJ 1989, 299, 1259–1260. [Google Scholar] [CrossRef] [Green Version]
  6. Stiemsma, L.T.; Reynolds, L.A.; Turvey, S.E.; Finlay, B.B. The Hygiene Hypothesis: Current Perspectives and Future Therapies. Immunotargets Ther. 2015, 4, 143–157. [Google Scholar] [CrossRef] [Green Version]
  7. Hanski, I.; von Hertzen, L.; Fyhrquist, N.; Koskinen, K.; Torppa, K.; Laatikainen, T.; Karisola, P.; Auvinen, P.; Paulin, L.; Mäkelä, M.J.; et al. Environmental Biodiversity, Human Microbiota, and Allergy Are Interrelated. Proc. Natl. Acad. Sci. USA 2012, 109, 8334–8339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Rook, G.A.W.; Adams, V.; Hunt, J.; Palmer, R.; Martinelli, R.; Brunet, L.R. Mycobacteria and Other Environmental Organisms as Immunomodulators for Immunoregulatory Disorders. Springer Semin. Immun. 2004, 25, 237–255. [Google Scholar] [CrossRef]
  9. Noverr, M.C.; Huffnagle, G.B. The “microflora Hypothesis” of Allergic Diseases. Clin. Exp. Allergy 2005, 35, 1511–1520. [Google Scholar] [CrossRef] [PubMed]
  10. Ege, M.J.; Mayer, M.; Schwaiger, K.; Mattes, J.; Pershagen, G.; van Hage, M.; Scheynius, A.; Bauer, J.; von Mutius, E. Environmental Bacteria and Childhood Asthma. Allergy 2012, 67, 1565–1571. [Google Scholar] [CrossRef] [PubMed]
  11. Valkonen, M.; Wouters, I.M.; Täubel, M.; Rintala, H.; Lenters, V.; Vasara, R.; Genuneit, J.; Braun-Fahrländer, C.; Piarroux, R.; von Mutius, E.; et al. Bacterial Exposures and Associations with Atopy and Asthma in Children. PLoS ONE 2015, 10, e0131594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Stein, M.M.; Hrusch, C.L.; Gozdz, J.; Igartua, C.; Pivniouk, V.; Murray, S.E.; Ledford, J.G.; Marques dos Santos, M.; Anderson, R.L.; Metwali, N.; et al. Innate Immunity and Asthma Risk in Amish and Hutterite Farm Children. N. Engl. J. Med. 2016, 375, 411–421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Kondrashova, A.; Seiskari, T.; Ilonen, J.; Knip, M.; Hyöty, H. The ‘Hygiene Hypothesis’ and the Sharp Gradient in the Incidence of Autoimmune and Allergic Diseases between Russian Karelia and Finland. APMIS 2013, 121, 478–493. [Google Scholar] [CrossRef]
  14. Sala, O.E.; Chapin, F.S., III; Armesto, J.J.; Berlow, E.; Bloomfield, J.; Dirzo, R.; Huber-Sanwald, E.; Huenneke, L.F.; Jackson, R.B. Global Biodiversity Scenarios for the Year 2100. Science 2000, 287, 1770–1774. [Google Scholar] [CrossRef] [PubMed]
  15. Parajuli, A.; Grönroos, M.; Kauppi, S.; Płociniczak, T.; Roslund, M.I.; Galitskaya, P.; Laitinen, O.H.; Hyöty, H.; Jumpponen, A.; Strömmer, R.; et al. The Abundance of Health-Associated Bacteria Is Altered in PAH Polluted Soils—Implications for Health in Urban Areas? PLoS ONE 2017, 12, e0187852. [Google Scholar] [CrossRef] [Green Version]
  16. Roslund, M.I.; Rantala, S.; Oikarinen, S.; Puhakka, R.; Hui, N.; Parajuli, A.; Laitinen, O.H.; Hyöty, H.; Rantalainen, A.-L.; Sinkkonen, A.; et al. Endocrine Disruption and Commensal Bacteria Alteration Associated with Gaseous and Soil PAH Contamination among Daycare Children. Environ. Int. 2019, 130, 104894. [Google Scholar] [CrossRef] [PubMed]
  17. Vari, H.K.; Roslund, M.I.; Oikarinen, S.; Nurminen, N.; Puhakka, R.; Parajuli, A.; Grönroos, M.; Siter, N.; Laitinen, O.H.; Hyöty, H.; et al. Associations between Land Cover Categories, Gaseous PAH Levels in Ambient Air and Endocrine Signaling Predicted from Gut Bacterial Metagenome of the Elderly. Chemosphere 2020, 128965. [Google Scholar] [CrossRef]
  18. von Hertzen, L.; Hanski, I.; Haahtela, T. Natural Immunity. EMBO Rep. 2011, 12, 1089–1093. [Google Scholar] [CrossRef] [PubMed]
  19. Haahtela, T.; Laatikainen, T.; Alenius, H.; Auvinen, P.; Fyhrquist, N.; Hanski, I.; von Hertzen, L.; Jousilahti, P.; Kosunen, T.U.; Markelova, O.; et al. Hunt for the Origin of Allergy—Comparing the Finnish and Russian Karelia. Clin. Exp. Allergy 2015, 45, 891–901. [Google Scholar] [CrossRef] [PubMed]
  20. Haahtela, T. A Biodiversity Hypothesis. Allergy 2019, 74, 1445–1456. [Google Scholar] [CrossRef] [Green Version]
  21. von Hertzen, L.; Haahtela, T. Disconnection of Man and the Soil: Reason for the Asthma and Atopy Epidemic? J. Allergy Clin. Immunol. 2006, 117, 334–344. [Google Scholar] [CrossRef] [PubMed]
  22. Ruokolainen, L.; von Hertzen, L.; Fyhrquist, N.; Laatikainen, T.; Lehtomäki, J.; Auvinen, P.; Karvonen, A.M.; Hyvärinen, A.; Tillmann, V.; Niemelä, O.; et al. Green Areas around Homes Reduce Atopic Sensitization in Children. Allergy 2015, 70, 195–202. [Google Scholar] [CrossRef]
  23. Parajuli, A.; Hui, N.; Puhakka, R.; Oikarinen, S.; Grönroos, M.; Selonen, V.A.O.; Siter, N.; Kramna, L.; Roslund, M.I.; Vari, H.K.; et al. Yard Vegetation Is Associated with Gut Microbiota Composition. Sci. Total Environ. 2020, 713, 136707. [Google Scholar] [CrossRef] [PubMed]
  24. Lehtimäki, J.; Karkman, A.; Laatikainen, T.; Paalanen, L.; von Hertzen, L.; Haahtela, T.; Hanski, I.; Ruokolainen, L. Patterns in the Skin Microbiota Differ in Children and Teenagers between Rural and Urban Environments. Sci. Rep. 2017, 7, 45651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Parajuli, A.; Grönroos, M.; Siter, N.; Puhakka, R.; Vari, H.K.; Roslund, M.I.; Jumpponen, A.; Nurminen, N.; Laitinen, O.H.; Hyöty, H.; et al. Urbanization Reduces Transfer of Diverse Environmental Microbiota Indoors. Front. Microbiol. 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  26. Macpherson, A.J.; Harris, N.L. Interactions between Commensal Intestinal Bacteria and the Immune System. Nat. Rev. Immunol. 2004, 4, 478–485. [Google Scholar] [CrossRef]
  27. Frank, D.N.; Amand, A.L.S.; Feldman, R.A.; Boedeker, E.C.; Harpaz, N.; Pace, N.R. Molecular-Phylogenetic Characterization of Microbial Community Imbalances in Human Inflammatory Bowel Diseases. Proc. Natl. Acad. Sci. USA 2007, 104, 13780–13785. [Google Scholar] [CrossRef] [Green Version]
  28. Round, J.L.; Mazmanian, S.K. The Gut Microbiota Shapes Intestinal Immune Responses during Health and Disease. Nat. Rev. Immunol. 2009, 9, 313–323. [Google Scholar] [CrossRef]
  29. Sekirov, I.; Russell, S.L.; Antunes, L.C.M.; Finlay, B.B. Gut Microbiota in Health and Disease. Physiol. Rev. 2010, 90, 859–904. [Google Scholar] [CrossRef] [Green Version]
  30. Clarke, G.; Stilling, R.M.; Kennedy, P.J.; Stanton, C.; Cryan, J.F.; Dinan, T.G. Minireview: Gut Microbiota: The Neglected Endocrine Organ. Mol. Endocrinol. 2014, 28, 1221–1238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Meng, X.; Zhang, G.; Cao, H.; Yu, D.; Fang, X.; de Vos, W.M.; Wu, H. Gut Dysbacteriosis and Intestinal Disease: Mechanism and Treatment. J. Appl. Microbiol. 2020, 129, 787–805. [Google Scholar] [CrossRef] [Green Version]
  32. Ott, S.J.; Musfeldt, M.; Wenderoth, D.F.; Hampe, J.; Brant, O.; Fölsch, U.R.; Timmis, K.N.; Schreiber, S. Reduction in Diversity of the Colonic Mucosa Associated Bacterial Microflora in Patients with Active Inflammatory Bowel Disease. Gut 2004, 53, 685–693. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Scher, J.U.; Ubeda, C.; Artacho, A.; Attur, M.; Isaac, S.; Reddy, S.M.; Marmon, S.; Neimann, A.; Brusca, S.; Patel, T.; et al. Decreased Bacterial Diversity Characterizes the Altered Gut Microbiota in Patients with Psoriatic Arthritis, Resembling Dysbiosis in Inflammatory Bowel Disease. Arthritis Rheumatol. 2015, 67, 128–139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. de Paiva, C.S.; Jones, D.B.; Stern, M.E.; Bian, F.; Moore, Q.L.; Corbiere, S.; Streckfus, C.F.; Hutchinson, D.S.; Ajami, N.J.; Petrosino, J.F.; et al. Altered Mucosal Microbiome Diversity and Disease Severity in Sjögren Syndrome. Sci. Rep. 2016, 6, 23561. [Google Scholar] [CrossRef]
  35. Manichanh, C.; Rigottier-Gois, L.; Bonnaud, E.; Gloux, K.; Pelletier, E.; Frangeul, L.; Nalin, R.; Jarrin, C.; Chardon, P.; Marteau, P.; et al. Reduced Diversity of Faecal Microbiota in Crohn’s Disease Revealed by a Metagenomic Approach. Gut 2006, 55, 205–211. [Google Scholar] [CrossRef] [Green Version]
  36. Rajilić–Stojanović, M.; Biagi, E.; Heilig, H.G.H.J.; Kajander, K.; Kekkonen, R.A.; Tims, S.; de Vos, W.M. Global and Deep Molecular Analysis of Microbiota Signatures in Fecal Samples From Patients With Irritable Bowel Syndrome. Gastroenterology 2011, 141, 1792–1801. [Google Scholar] [CrossRef] [PubMed]
  37. Louis, P.; Flint, H.J. Diversity, Metabolism and Microbial Ecology of Butyrate-Producing Bacteria from the Human Large Intestine. Fems. Microbiol. Lett. 2009, 294, 1–8. [Google Scholar] [CrossRef] [Green Version]
  38. Plöger, S.; Stumpff, F.; Penner, G.; Schulzke, J.; Gäbel, G.; Martens, H.; Shen, Z.; Günzel, D.; Aschenbach, J. Microbial Butyrate and Its Role for Barrier Function in the Gastrointestinal Tract. Ann. N. Y. Acad. Sci. 2012, 1258, 52–59. [Google Scholar] [CrossRef] [Green Version]
  39. Geirnaert, A.; Calatayud, M.; Grootaert, C.; Laukens, D.; Devriese, S.; Smagghe, G.; De Vos, M.; Boon, N.; Van de Wiele, T. Butyrate-Producing Bacteria Supplemented in Vitro to Crohn’s Disease Patient Microbiota Increased Butyrate Production and Enhanced Intestinal Epithelial Barrier Integrity. Sci. Rep. 2017, 7, 11450. [Google Scholar] [CrossRef] [Green Version]
  40. Greenhalgh, K.; Meyer, K.M.; Aagaard, K.M.; Wilmes, P. The Human Gut Microbiome in Health: Establishment and Resilience of Microbiota over a Lifetime. Environ. Microbiol. 2016, 18, 2103–2116. [Google Scholar] [CrossRef]
  41. Rinninella, E.; Raoul, P.; Cintoni, M.; Franceschi, F.; Miggiano, G.A.D.; Gasbarrini, A.; Mele, M.C. What Is the Healthy Gut Microbiota Composition? A Changing Ecosystem across Age, Environment, Diet, and Diseases. Microorganisms 2019, 7, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. O’Toole, P.W.; Claesson, M.J. Gut Microbiota: Changes throughout the Lifespan from Infancy to Elderly. Int. Dairy J. 2010, 20, 281–291. [Google Scholar] [CrossRef]
  43. Nagpal, R.; Mainali, R.; Ahmadi, S.; Wang, S.; Singh, R.; Kavanagh, K.; Kitzman, D.W.; Kushugulova, A.; Marotta, F.; Yadav, H. Gut Microbiome and Aging: Physiological and Mechanistic Insights. Nutr. Healthy Aging 2018, 4, 267–285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Nurminen, N.; Lin, J.; Grönroos, M.; Puhakka, R.; Kramna, L.; Vari, H.K.; Viskari, H.; Oikarinen, S.; Roslund, M.; Parajuli, A.; et al. Nature-Derived Microbiota Exposure as a Novel Immunomodulatory Approach. Future Microbiol. 2018, 13, 737–744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Grönroos, M.; Parajuli, A.; Laitinen, O.H.; Roslund, M.I.; Vari, H.K.; Hyöty, H.; Puhakka, R.; Sinkkonen, A. Short-Term Direct Contact with Soil and Plant Materials Leads to an Immediate Increase in Diversity of Skin Microbiota. MicrobiologyOpen 2019, 8, e00645. [Google Scholar] [CrossRef]
  46. Hui, N.; Grönroos, M.; Roslund, M.; Parajuli, A.; Vari, H.K.; Soininen, L.; Laitinen, O.H.; Sinkkonen, A. Diverse Environmental Microbiota as a Tool to Augment Biodiversity in Urban Landscaping Materials. Front. Microbiol. 2019, 10, 536. [Google Scholar] [CrossRef] [PubMed]
  47. Roslund, M.I.; Puhakka, R.; Grönroos, M.; Nurminen, N.; Oikarinen, S.; Gazali, A.M.; Cinek, O.; Kramná, L.; Siter, N.; Vari, H.K.; et al. Biodiversity Intervention Enhances Immune Regulation and Health-Associated Commensal Microbiota among Daycare Children. Sci. Adv. 2020, 6, eaba2578. [Google Scholar] [CrossRef]
  48. Lehtimäki, J.; Sinkko, H.; Hielm-Björkman, A.; Salmela, E.; Tiira, K.; Laatikainen, T.; Mäkeläinen, S.; Kaukonen, M.; Uusitalo, L.; Hanski, I.; et al. Skin Microbiota and Allergic Symptoms Associate with Exposure to Environmental Microbes. Proc. Natl. Acad. Sci. USA 2018, 115, 4897–4902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Sievänen, T.; Neuvonen, M. Luonnon virkistyskäytön kysyntä 2010 ja kysynnän muutos. In Luonnon Virkistyskäyttö 2010; Sievänen, T., Neuvonen, M., Eds.; Working Papers of the Finnish Forest Research Institute: Vantaa, Finland, 2011; Volume 212, pp. 37–79. [Google Scholar]
  50. Pitkänen, K.; Lehtimäki, J.; Puhakka, R. How Do Rural Second Homes Affect Human Health and Well-Being? Review of Potential Impacts. Int. J. Environ. Res. Public Health 2020, 17, 6748. [Google Scholar] [CrossRef]
  51. Adamiak, C.; Vepsäläinen, M.; Strandell, A.; Hiltunen, M.; Pitkänen, K.; Hall, M.; Rinne, J.; Hannonen, O.; Paloniemi, R.; Åkerlund, U. Second Home Tourism in Finland—Perceptions of Citizens and Municipalities on the State and Development of Second Home Tourism; Suomen Ympäristökeskus: Helsinki, Finland, 2015. [Google Scholar]
  52. Statistics Finland. Free-Time Residences 2019. Helsinki. 2020. Available online: http://www.stat.fi/til/rakke/2019/rakke_2019_2020-05-27_kat_001_en.html (accessed on 2 December 2020).
  53. Pitkänen, K.; Hannonen, O.; Toso, S.; Gallent, N.; Hamiduddin, I.; Halseth, G.; Hall, M.C.; Müller, D.K.; Treivish, A.; Nefedova, T. Second Homes during Corona—Safe or Unsafe Haven and for Whom? Reflections from Researchers around the World. Matkailututkimus 2020, 16, 20–39. [Google Scholar] [CrossRef]
  54. Sievänen, T.; Pouta, E.; Neuvonen, M. Recreational Home Users—Potential Clients for Countryside Tourism? Scand. J. Hosp. Tour. 2007, 7, 223–242. [Google Scholar] [CrossRef]
  55. FCG Finnish Consulting Group Oy. Mökkibarometri 2016; Saaristoasiain Neuvottelukunta, Maa-ja Metsätalousministeriö: Helsinki, Finland, 2016; Available online: https://valtioneuvosto.fi/documents/1410837/1880296/Mokkibarometri+2016/7b69ab48-5859-4b55-8dc2-5514cdfa6000 (accessed on 22 December 2020).
  56. Hui, N.; Parajuli, A.; Puhakka, R.; Grönroos, M.; Roslund, M.I.; Vari, H.K.; Selonen, V.A.O.; Yan, G.; Siter, N.; Nurminen, N.; et al. Temporal Variation in Indoor Transfer of Dirt-Associated Environmental Bacteria in Agricultural and Urban Areas. Environ. Int. 2019, 132, 105069. [Google Scholar] [CrossRef]
  57. Kondrashova, A.; Reunanen, A.; Romanov, A.; Karvonen, A.; Viskari, H.; Vesikari, T.; Ilonen, J.; Knip, M.; Hyöty, H. A Six-fold Gradient in the Incidence of Type 1 Diabetes at the Eastern Border of Finland. Ann. Med. 2005, 37, 67–72. [Google Scholar] [CrossRef] [PubMed]
  58. Fogelholm, M.; Valve, R.; Absetz, P.; Heinonen, H.; Uutela, A.; Patja, K.; Karisto, A.; Konttinen, R.; Mäkelä, T.; Nissinen, A.; et al. Rural—Urban Differences in Health and Health Behaviour: A Baseline Description of a Community Health-Promotion Programme for the Elderly. Scand. J. Public Health 2006, 34, 632–640. [Google Scholar] [CrossRef] [PubMed]
  59. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; et al. Introducing Mothur: Open-Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial Communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Schloss, P.D.; Westcott, S.L. Assessing and Improving Methods Used in Operational Taxonomic Unit-Based Approaches for 16S RRNA Gene Sequence Analysis. Appl. Environ. Microbiol. 2011, 77, 3219–3226. [Google Scholar] [CrossRef] [Green Version]
  61. Kozich, J.J.; Westcott, S.L.; Baxter, N.T.; Highlander, S.K.; Schloss, P.D. Development of a Dual-Index Sequencing Strategy and Curation Pipeline for Analyzing Amplicon Sequence Data on the MiSeq Illumina Sequencing Platform. Appl. Environ. Microbiol. 2013, 79, 5112–5120. [Google Scholar] [CrossRef] [Green Version]
  62. Pruesse, E.; Quast, C.; Knittel, K.; Fuchs, B.M.; Ludwig, W.; Peplies, J.; Glöckner, F.O. SILVA: A Comprehensive Online Resource for Quality Checked and Aligned Ribosomal RNA Sequence Data Compatible with ARB. Nucleic Acids Res. 2007, 35, 7188–7196. [Google Scholar] [CrossRef] [Green Version]
  63. Oliver, A.K.; Brown, S.P.; Callaham, M.A.; Jumpponen, A. Polymerase Matters: Non-Proofreading Enzymes Inflate Fungal Community Richness Estimates by up to 15%. Fungal Ecol. 2015, 15, 86–89. [Google Scholar] [CrossRef] [Green Version]
  64. DeSantis, T.Z.; Hugenholtz, P.; Larsen, N.; Rojas, M.; Brodie, E.L.; Keller, K.; Huber, T.; Dalevi, D.; Hu, P.; Andersen, G.L. Greengenes, a Chimera-Checked 16S RRNA Gene Database and Workbench Compatible with ARB. Appl. Environ. Microbiol. 2006, 72, 5069–5072. [Google Scholar] [CrossRef] [Green Version]
  65. Langille, M.G.I.; Zaneveld, J.; Caporaso, J.G.; McDonald, D.; Knights, D.; Reyes, J.A.; Clemente, J.C.; Burkepile, D.E.; Vega Thurber, R.L.; Knight, R.; et al. Predictive Functional Profiling of Microbial Communities Using 16S RRNA Marker Gene Sequences. Nat. Biotechnol. 2013, 31, 814–821. [Google Scholar] [CrossRef] [PubMed]
  66. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  67. Donaldson, G.P.; Lee, S.M.; Mazmanian, S.K. Gut Biogeography of the Bacterial Microbiota. Nat. Rev. Microbiol. 2016, 14, 20–32. [Google Scholar] [CrossRef] [Green Version]
  68. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 22 December 2020).
  69. McMurdie, P.J.; Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [Green Version]
  70. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016; ISBN 978-0-387-98141-3. [Google Scholar]
  71. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R Package Version 2.5-6. 2019. Available online: https://CRAN.R-project.org/package=vegan (accessed on 22 December 2020).
  72. Bartosch, S.; Fite, A.; Macfarlane, G.T.; McMurdo, M.E.T. Characterization of Bacterial Communities in Feces from Healthy Elderly Volunteers and Hospitalized Elderly Patients by Using Real-Time PCR and Effects of Antibiotic Treatment on the Fecal Microbiota. Appl. Environ. Microbiol. 2004, 70, 3575–3581. [Google Scholar] [CrossRef] [Green Version]
  73. Woodmansey, E.J.; McMurdo, M.E.T.; Macfarlane, G.T.; Macfarlane, S. Comparison of Compositions and Metabolic Activities of Fecal Microbiotas in Young Adults and in Antibiotic-Treated and Non-Antibiotic-Treated Elderly Subjects. Appl. Environ. Microbiol. 2004, 70, 6113–6122. [Google Scholar] [CrossRef] [Green Version]
  74. Zwielehner, J.; Liszt, K.; Handschur, M.; Lassl, C.; Lapin, A.; Haslberger, A.G. Combined PCR-DGGE Fingerprinting and Quantitative-PCR Indicates Shifts in Fecal Population Sizes and Diversity of Bacteroides, Bifidobacteria and Clostridium Cluster IV in Institutionalized Elderly. Exp. Gerontol. 2009, 44, 440–446. [Google Scholar] [CrossRef] [PubMed]
  75. Mäkivuokko, H.; Tiihonen, K.; Tynkkynen, S.; Paulin, L.; Rautonen, N. The Effect of Age and Non-Steroidal Anti-Inflammatory Drugs on Human Intestinal Microbiota Composition. Br. J. Nutr. 2010, 103, 227–234. [Google Scholar] [CrossRef] [Green Version]
  76. Kostic, A.D.; Gevers, D.; Pedamallu, C.S.; Michaud, M.; Duke, F.; Earl, A.M.; Ojesina, A.I.; Jung, J.; Bass, A.J.; Tabernero, J.; et al. Genomic Analysis Identifies Association of Fusobacterium with Colorectal Carcinoma. Genome Res. 2012, 22, 292–298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Larsen, N.; Vogensen, F.K.; van den Berg, F.W.J.; Nielsen, D.S.; Andreasen, A.S.; Pedersen, B.K.; Al-Soud, W.A.; Sørensen, S.J.; Hansen, L.H.; Jakobsen, M. Gut Microbiota in Human Adults with Type 2 Diabetes Differs from Non-Diabetic Adults. PLoS ONE 2010, 5, e9085. [Google Scholar] [CrossRef]
  78. Krogius-Kurikka, L.; Lyra, A.; Malinen, E.; Aarnikunnas, J.; Tuimala, J.; Paulin, L.; Mäkivuokko, H.; Kajander, K.; Palva, A. Microbial Community Analysis Reveals High Level Phylogenetic Alterations in the Overall Gastrointestinal Microbiota of Diarrhoea-Predominant Irritable Bowel Syndrome Sufferers. BMC Gastroenterol. 2009, 9, 95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Matsuoka, K.; Kanai, T. The Gut Microbiota and Inflammatory Bowel Disease. Semin. Immunopathol. 2015, 37, 47–55. [Google Scholar] [CrossRef] [Green Version]
  80. Zeng, X.; Gao, X.; Peng, Y.; Wu, Q.; Zhu, J.; Tan, C.; Xia, G.; You, C.; Xu, R.; Pan, S.; et al. Higher Risk of Stroke Is Correlated With Increased Opportunistic Pathogen Load and Reduced Levels of Butyrate-Producing Bacteria in the Gut. Front. Cell. Infect. Microbiol. 2019, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Chiu, Y.-H.; Macmillan, J.B.; Chen, Z.J. RNA Polymerase III Detects Cytosolic DNA and Induces Type I Interferons through the RIG-I Pathway. Cell 2009, 138, 576–591. [Google Scholar] [CrossRef] [Green Version]
  82. Martínez, I.; Oliveros, J.C.; Cuesta, I.; de la Barrera, J.; Ausina, V.; Casals, C.; de Lorenzo, A.; García, E.; García-Fojeda, B.; Garmendia, J.; et al. Apoptosis, Toll-like, RIG-I-like and NOD-like Receptors Are Pathways Jointly Induced by Diverse Respiratory Bacterial and Viral Pathogens. Front. Microbiol. 2017, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Laitinen, O.H.; Svedin, E.; Kapell, S.; Nurminen, A.; Hytönen, V.P.; Flodström-Tullberg, M. Enteroviral Proteases: Structure, Host Interactions and Pathogenicity. Rev. Med. Virol. 2016, 26, 251–267. [Google Scholar] [CrossRef] [PubMed]
  84. Oshiumi, H. Recent Advances and Contradictions in the Study of the Individual Roles of Ubiquitin Ligases That Regulate RIG-I-Like Receptor-Mediated Antiviral Innate Immune Responses. Front. Immunol. 2020, 11, 1296. [Google Scholar] [CrossRef] [PubMed]
  85. Ren, Z.; Ding, T.; Zuo, Z.; Xu, Z.; Deng, J.; Wei, Z. Regulation of MAVS Expression and Signaling Function in the Antiviral Innate Immune Response. Front. Immunol. 2020, 11, 1030. [Google Scholar] [CrossRef]
  86. Sharma, S.; Fitzgerald, K.A.; Cancro, M.P.; Marshak-Rothstein, A. Nucleic Acid-Sensing Receptors: Rheostats of Autoimmunity and Autoinflammation. J. Immunol. 2015, 195, 3507–3512. [Google Scholar] [CrossRef] [Green Version]
  87. Baechler, E.C.; Batliwalla, F.M.; Karypis, G.; Gaffney, P.M.; Ortmann, W.A.; Espe, K.J.; Shark, K.B.; Grande, W.J.; Hughes, K.M.; Kapur, V.; et al. Interferon-Inducible Gene Expression Signature in Peripheral Blood Cells of Patients with Severe Lupus. Proc. Natl. Acad. Sci. USA 2003, 100, 2610–2615. [Google Scholar] [CrossRef] [Green Version]
  88. Bennett, L.; Palucka, A.K.; Arce, E.; Cantrell, V.; Borvak, J.; Banchereau, J.; Pascual, V. Interferon and Granulopoiesis Signatures in Systemic Lupus Erythematosus Blood. J. Exp. Med. 2003, 197, 711–723. [Google Scholar] [CrossRef] [Green Version]
  89. Lübbers, J.; Brink, M.; van de Stadt, L.A.; Vosslamber, S.; Wesseling, J.G.; van Schaardenburg, D.; Rantapää-Dahlqvist, S.; Verweij, C.L. The Type I IFN Signature as a Biomarker of Preclinical Rheumatoid Arthritis. Ann. Rheum. Dis. 2013, 72, 776–780. [Google Scholar] [CrossRef] [PubMed]
  90. Postal, M.; Vivaldo, J.F.; Fernandez-Ruiz, R.; Paredes, J.L.; Appenzeller, S.; Niewold, T.B. Type I Interferon in the Pathogenesis of Systemic Lupus Erythematosus. Curr. Opin. Immunol. 2020, 67, 87–94. [Google Scholar] [CrossRef] [PubMed]
  91. Tossberg, J.T.; Heinrich, R.M.; Farley, V.M.; Crooke, P.S.; Aune, T.M. Adenosine-to-Inosine RNA Editing of Alu Double-Stranded (Ds)RNAs Is Markedly Decreased in Multiple Sclerosis and Unedited Alu DsRNAs Are Potent Activators of Proinflammatory Transcriptional Responses. J. Immunol. 2020, 205, 2606–2617. [Google Scholar] [CrossRef]
  92. Claesson, M.J.; Jeffery, I.B.; Conde, S.; Power, S.E.; O’Connor, E.M.; Cusack, S.; Harris, H.M.B.; Coakley, M.; Lakshminarayanan, B.; O’Sullivan, O.; et al. Gut Microbiota Composition Correlates with Diet and Health in the Elderly. Nature 2012, 488, 178–184. [Google Scholar] [CrossRef] [PubMed]
  93. Mueller, S.; Saunier, K.; Hanisch, C.; Norin, E.; Alm, L.; Midtvedt, T.; Cresci, A.; Silvi, S.; Orpianesi, C.; Verdenelli, M.C.; et al. Differences in Fecal Microbiota in Different European Study Populations in Relation to Age, Gender, and Country: A Cross-Sectional Study. Appl. Environ. Microbiol. 2006, 72, 1027–1033. [Google Scholar] [CrossRef] [Green Version]
  94. Hertzen, L.V.; Laatikainen, T.; Pitkänen, T.; Vlasoff, T.; Mäkelä, M.J.; Vartiainen, E.; Haahtela, T. Microbial Content of Drinking Water in Finnish and Russian Karelia—Implications for Atopy Prevalence. Allergy 2007, 62, 288–292. [Google Scholar] [CrossRef]
  95. Hesselmar, B.; Hicke-Roberts, A.; Wennergren, G. Allergy in Children in Hand Versus Machine Dishwashing. Pediatrics 2015, 135, e590–e597. [Google Scholar] [CrossRef] [Green Version]
  96. Turkki, P. Lähiruokaa Mökille: Etelä-Savon Vapaa-Ajan Asukkaiden Ruokaostoskäyttäytyminen ja Lähiruoan Saatavuus; Tutkimuksia ja Raportteja—Research Reports 106; Mikkelin Ammattikorkeakoulu: Mikkeli, Finland, 2016; ISBN 978-951-588-553-1. [Google Scholar]
  97. Martínez, I.; Kim, J.; Duffy, P.R.; Schlegel, V.L.; Walter, J. Resistant Starches Types 2 and 4 Have Differential Effects on the Composition of the Fecal Microbiota in Human Subjects. PLoS ONE 2010, 5, e15046. [Google Scholar] [CrossRef] [Green Version]
  98. David, L.A.; Maurice, C.F.; Carmody, R.N.; Gootenberg, D.B.; Button, J.E.; Wolfe, B.E.; Ling, A.V.; Devlin, A.S.; Varma, Y.; Fischbach, M.A.; et al. Diet Rapidly and Reproducibly Alters the Human Gut Microbiome. Nature 2014, 505, 559–563. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Venkataraman, A.; Sieber, J.R.; Schmidt, A.W.; Waldron, C.; Theis, K.R.; Schmidt, T.M. Variable Responses of Human Microbiomes to Dietary Supplementation with Resistant Starch. Microbiome 2016, 4, 33. [Google Scholar] [CrossRef] [Green Version]
  100. Lin, D.; Peters, B.A.; Friedlander, C.; Freiman, H.J.; Goedert, J.J.; Sinha, R.; Miller, G.; Bernstein, M.A.; Hayes, R.B.; Ahn, J. Association of Dietary Fibre Intake and Gut Microbiota in Adults. Br. J. Nutr. 2018, 120, 1014–1022. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Bjørkhaug, S.T.; Aanes, H.; Neupane, S.P.; Bramness, J.G.; Malvik, S.; Henriksen, C.; Skar, V.; Medhus, A.W.; Valeur, J. Characterization of Gut Microbiota Composition and Functions in Patients with Chronic Alcohol Overconsumption. Gut Microbes 2019, 10, 663–675. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Wan, Y.; Wang, F.; Yuan, J.; Li, J.; Jiang, D.; Zhang, J.; Li, H.; Wang, R.; Tang, J.; Huang, T.; et al. Effects of Dietary Fat on Gut Microbiota and Faecal Metabolites, and Their Relationship with Cardiometabolic Risk Factors: A 6-Month Randomised Controlled-Feeding Trial. Gut 2019, 68, 1417–1429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Flint, H.J.; Duncan, S.H.; Scott, K.P.; Louis, P. Links between Diet, Gut Microbiota Composition and Gut Metabolism. Proc. Nutr. Soc. 2015, 74, 13–22. [Google Scholar] [CrossRef] [Green Version]
Ijerph 18 03742 g001 550
Figure 1. Relative abundances of bacterial phyla in stool samples of 10 rural-second-home users (cottage; blue color) and 15 control group participants (no cottage; orange color) (values are expressed as mean ± SE).
Figure 1. Relative abundances of bacterial phyla in stool samples of 10 rural-second-home users (cottage; blue color) and 15 control group participants (no cottage; orange color) (values are expressed as mean ± SE).
Ijerph 18 03742 g001
Ijerph 18 03742 g002a 550Ijerph 18 03742 g002b 550Ijerph 18 03742 g002c 550
Figure 2. The Shannon diversity of fecal microbiota was higher in the control group (no cottage; red color) than among rural-second-home users (cottage; blue color) at all tested taxonomic levels (operational taxonomic units (OTU) (a), Family (b), Order (c), Class (d), Phylum (e) and for OTUs within the phylum Firmicutes (f) and the families Lachnospiraceae (g) and Ruminococcaceae (h). The observed richness was higher in the control group than among rural-second-home users at all tested taxonomic levels and within the phylum Firmicutes and the family Ruminococcaceae. Simpson diversity was higher in the control group only within the phylum Firmicutes. Boxplots show medians (thick line), upper and lower hinges (box), values 1.5 times IQR (whiskers), and values outside hinges (data points, outliers).
Figure 2. The Shannon diversity of fecal microbiota was higher in the control group (no cottage; red color) than among rural-second-home users (cottage; blue color) at all tested taxonomic levels (operational taxonomic units (OTU) (a), Family (b), Order (c), Class (d), Phylum (e) and for OTUs within the phylum Firmicutes (f) and the families Lachnospiraceae (g) and Ruminococcaceae (h). The observed richness was higher in the control group than among rural-second-home users at all tested taxonomic levels and within the phylum Firmicutes and the family Ruminococcaceae. Simpson diversity was higher in the control group only within the phylum Firmicutes. Boxplots show medians (thick line), upper and lower hinges (box), values 1.5 times IQR (whiskers), and values outside hinges (data points, outliers).
Ijerph 18 03742 g002aIjerph 18 03742 g002bIjerph 18 03742 g002c
Ijerph 18 03742 g003 550
Figure 3. Fecal bacterial community composition did not differ between rural-second-home users and control group (p = 0.139, stress = 0.178, F = 1.327, R2 = 0.055). Non-metric dimensional scaling ordination based on Bray–Curtis dissimilarity of bacteria in stool samples of urban rural-second-home users (cottage; blue triangles) and control group (no cottage; orange dots) at the operational taxonomic unit (OTU) level is shown.
Figure 3. Fecal bacterial community composition did not differ between rural-second-home users and control group (p = 0.139, stress = 0.178, F = 1.327, R2 = 0.055). Non-metric dimensional scaling ordination based on Bray–Curtis dissimilarity of bacteria in stool samples of urban rural-second-home users (cottage; blue triangles) and control group (no cottage; orange dots) at the operational taxonomic unit (OTU) level is shown.
Ijerph 18 03742 g003
Table
Table 1. Rural second homes had been used on average for 36 years, and the number of visiting days per year was on average 67. Most participants reported spending more time outdoors than indoors while visiting rural second homes.
Table 1. Rural second homes had been used on average for 36 years, and the number of visiting days per year was on average 67. Most participants reported spending more time outdoors than indoors while visiting rural second homes.
Years UsedVisiting Days per YearTime Spent Outdoors vs. Indoors
4640Outdoors
2310Outdoors
3035Outdoors
2790Equally
-100Outdoors
4430Outdoors
29150Equally
5330Outdoors
61-Outdoors
9120Equally
Table
Table 2. Dominant bacterial phyla in the gut and predominant bacterial families in the colon [67].
Table 2. Dominant bacterial phyla in the gut and predominant bacterial families in the colon [67].
TaxaRole and GeographyNote
BacteroidetesDominant phylum in the gut
FirmicutesDominant phylum in the gut
ActinobacteriaDominant phylum in the gut
ProteobacteriaDominant phylum in the gut
VerrucomicrobiaDominant phylum in the gutNot analyzed due to frequent zero-valued observations
RuminococcaceaePredominant family in the inter-fold regions of the lumen
LachnospiraceaePredominant family in the inter-fold regions of the lumen
BacteroidaceaePredominant family in the digesta
PrevotellaceaePredominant family in the digestaNot analyzed due to frequent zero-valued observations
RikenellaceaePredominant family in the digestaNot analyzed due to frequent zero-valued observations
Table
Table 3. Summary of bacterial diversity results. p-values, Benjamini–Hochberg critical values, significance after the correction, and type of the statistical test used (U or t-test) are given for each group.
Table 3. Summary of bacterial diversity results. p-values, Benjamini–Hochberg critical values, significance after the correction, and type of the statistical test used (U or t-test) are given for each group.
Diversity MeasureLevelp-ValueBenjamini–Hochberg Critical Value [66]Significant with FDR 0.1Test
ShannonPhyl. Firmicutes0.0010.008YesU test
Fam. Ruminococcaceae0.0140.017Yest-test
Fam. Lachnospiraceae0.0170.025YesU test
Phylum0.0310.033YesU test
Class0.0310.033YesU test
Order0.0310.033YesU test
OTU0.0350.042YesU test
Family0.0350.042YesU test
Phyl. Proteobacteria0.1540.050Not-test
Phyl. Bacteroidetes0.2290.058Not-test
Fam. Bacteroidaceae0.5420.067NoU test
Phyl. Actinobacteria0.8680.075Not-test
SimpsonPhyl. Firmicutes0.0040.008YesU test
Fam. Ruminococcaceae0.0270.017NoU test
Fam. Lachnospiraceae0.0310.025NoU test
OTU0.1200.033NoU test
Phylum0.1200.033NoU test
Class0.1200.033NoU test
Order0.1200.033NoU test
Family0.1200.033NoU test
Phyl. Bacteroidetes0.2670.042NoU test
Fam. Bacteroidaceae0.5790.050NoU test
Phyl. Proteobacteria0.6800.058Not-test
Phyl. Actinobacteria0.8680.067NoU test
Observed richnessFam. Ruminococcaceae0.0030.008Yest-test
OTU0.0080.017Yest-test
Phyl. Firmicutes0.0090.025Yest-test
Family0.0180.033Yest-test
Class0.0190.042Yest-test
Order0.0190.05Yest-test
Phylum0.0240.058Yest-test
Phyl. Proteobacteria0.0530.067Not-test
Fam. Lachnospiraceae0.1260.075NoU test
Phyl. Actinobacteria0.5480.083Not-test
Phyl. Bacteroidetes0.7600.092NoU test
Fam. Bacteroidaceae0.8890.1Not-test
Table
Table 4. Summary of results from PERMANOVA analysis of fecal bacterial community composition.
Table 4. Summary of results from PERMANOVA analysis of fecal bacterial community composition.
Levelp-ValueFR2
Phyl. Proteobacteria0.0681.4410.059
Phyl. Actinobacteria0.0801.6610.067
OTU0.1391.3270.055
Fam. Bacteroidaceae0.1721.3760.056
Phyl. Bacteroidetes0.2541.2050.050
Phylum0.3001.1170.046
Order0.3181.0880.045
Class0.3241.0890.045
Fam. Lachnospiraceae0.3281.1040.046
Family0.3431.0710.044
Fam. Ruminococcaceae0.4080.9900.041
Phyl. Firmicutes0.4271.0240.043
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saarenpää, M.; Roslund, M.I.; Puhakka, R.; Grönroos, M.; Parajuli, A.; Hui, N.; Nurminen, N.; Laitinen, O.H.; Hyöty, H.; Cinek, O.; et al. Do Rural Second Homes Shape Commensal Microbiota of Urban Dwellers? A Pilot Study among Urban Elderly in Finland. Int. J. Environ. Res. Public Health 2021, 18, 3742. https://doi.org/10.3390/ijerph18073742

AMA Style

Saarenpää M, Roslund MI, Puhakka R, Grönroos M, Parajuli A, Hui N, Nurminen N, Laitinen OH, Hyöty H, Cinek O, et al. Do Rural Second Homes Shape Commensal Microbiota of Urban Dwellers? A Pilot Study among Urban Elderly in Finland. International Journal of Environmental Research and Public Health. 2021; 18(7):3742. https://doi.org/10.3390/ijerph18073742

Chicago/Turabian Style

Saarenpää, Mika, Marja I. Roslund, Riikka Puhakka, Mira Grönroos, Anirudra Parajuli, Nan Hui, Noora Nurminen, Olli H. Laitinen, Heikki Hyöty, Ondrej Cinek, and et al. 2021. "Do Rural Second Homes Shape Commensal Microbiota of Urban Dwellers? A Pilot Study among Urban Elderly in Finland" International Journal of Environmental Research and Public Health 18, no. 7: 3742. https://doi.org/10.3390/ijerph18073742

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop