2. Materials and Methods
2.1. General Context
Yakuts represent the largest ethnic group of the Sakha Republic of Eastern Siberia (Yakutia), with almost half a million individuals. Although their history remains largely unknown, they descend from a population that migrated northwards from the Lake Baikal area, following the expansion of the Great Mongolian Empire in the 13th century, and the present-day Yakuts are the descendants of these ancient Yakuts [
8]. Except for the Arctic and Antarctic, Yakutia is the coldest country in the world, with winter temperatures dropping below −70 °C and annual temperature amplitudes over 100 °C. People and cattle adapted to extreme climatic conditions, and they have developed a traditional lifestyle based on livestock herding, where horses, cattle and reindeer provided means of transport and sources of meat, milk, and clothing [
9]. The first contacts of the autochthonous population with the Russians date back to 1620 AD, but there was almost no trace of these contacts in material culture before 1689/1700 AD [
10]. Trappers and itinerant merchants made the first contact, but then trade developed. The inhabitants of Yakutia traded with Russia and arrivals from the East as well as China via the trading posts of Nertchinsk and then of Kiakhta, south on the Chinese and Mongolian borders. From the end of the 18th century, the ports of Okhotsk and Ayan on the Pacific coast became more and more critical for imports [
5]. Samples recovered from burials are not always random because, before the 19th century, most individuals were not buried. In the 17th century, only male clan chiefs were buried, whereas after 1700 AD, as many men as women were buried. They were mainly from the elite, and the assumed identity of several subjects and/or their status could be inferred from the associated artefacts, modes of burial, and sometimes local legends that specified who they are (chiefs, shamans) [
5]. Assimilation in the Russian Orthodox culture was fully accomplished after 1800 AD. Russian colonization had a far-reaching impact on the social sphere, including the rise of an economic elite in the early 18th century [
8], the reduction of seasonal mobility and the settlement in villages, and the later progressive decline of earlier shamanic religious practices in favor of Christianity [
5]. The exposure of the immunologically-naive Yakut population to Russian germs also resulted in massive epidemiological outbreaks of smallpox [
10], tuberculosis [
11] and some changes in cultural practices and diet that affected dental health [
12]. According to our historical research, it does not appear that cholera affected Yakutia.
2.2. Data Collection
We have excavated more than 150 frozen bodies in four geographical regions (31,000 km²) of the Sakha Republic (Yakutia) (
Figure 1) [
13]. We classified them into five chronological periods: before 1700, 1700–1750, 1750–1800, 1800–1850, and after 1850 AD [
14]. We categorized and analyzed the characteristics of these ancient subjects and studied their familial relationships, social status, gender, and age [
15]. There are some multiple graves with relatives who died within a short time of one another which will allow us to examine whether there are intra-family variations. Among these, 47 had well-preserved head hair not in contact with the sediment. The sample was genetically homogeneous [
8], and all subjects had the same hair color (black).
Hair samples were collected and stored during the autopsy of the subjects following procedures ruling out the possibility of outside contamination. In the laboratory, hair samples (n = 47) were photographed and decontaminated in 2 × 2 min baths in water, followed by 2 × 1 min baths in dichloromethane. After drying proximal hair segments (2 to 3 cm length), we cut them into small pieces (<1 mm length). We analyzed proximal hair segments (2 to 3 cm in length) close to the skull corresponding to the last few weeks of life. It is challenging to know if these few cm represent individual consumption accurately because, at certain times, there could be supply difficulties during part of the year. One element of the discussion will therefore concern the relationship between the analyses and other archaeological elements to specify to what extent our results can be considered representative.
2.3. Hair Analysis
Concentrations of the main compounds related to tea and tobacco consumption were analysed. In one case, we found remains of tea leaves in a teapot associated with a woman, and we analysed them (
Figure 2). For extraction, we incubated 20 mg of hair samples in methanol, and we used the solution obtained after centrifugation for analysis. Previous studies on nicotine and caffeine analyses in hair and the use of liquid chromatography with high resolution mass spectrometry detection (LC-HRMS) or with high tandem mass spectrometry detection (LC-MS/MS) tools for toxicological investigations in hair [
16,
17] supported our decision to use these two analytical methods to detect tobacco and tea products in our ancient population sample. We combined these two techniques according to previously published methods [
18,
19,
20] screening for drugs and toxic compounds, psychoactive drugs including natural alkaloids, and tea methylxanthines (caffeine, theobromine, theophylline) and tobacco alkaloids (nicotine) and metabolites (cotinine).
1—A large screening for drugs and toxics compounds was performed in hair samples using LC-HRMS. We used a Waters ultra-performance liquid chromatography system consisting of two binary solvent manager LC pumps, a sample manager autosampler, and an Acquity column manager oven. Mass spectrometry data were acquired using a XEVO G2-XS QTOF instrument (Waters, Manchester, UK) controlled with MassLynx 4.1 software. After the addition of 50 µL of internal standards (methyl-clonazepam and β-OH-ethyltheophyllin) to 100 µL of the methanolic hair extract, the obtained mixture was subsequently evaporated to dryness at +30 °C under a gentle stream of nitrogen. The dry residue was reconstituted using 50 µL of a mixture of ammonium formate buffer 5 mM, pH 3/acetonitrile in 1% formic acid (87/13; v/v): 10 µL were injected in the chromatographic system. Chromatographic separation was performed using an ACQUITY HSS C18 column (150 × 2.1 mm, 1.8 µm, Waters) in an oven at 50 °C, and mobile phases including ammonium formate buffer 5 mM at pH 3 and acetonitrile in 0.1% formic acid. For detection, mass spectrometric conditions were as follows: positive electrospray ionization interface (ESI+), ion spray voltage set at 20 V, source temperature set at 140 °C and desolvation temperature at 500 °C with a desolvation gas flow rate of 900 L/h, nitrogen as desolvation gas, and argon as collision gas. The data were processed using ChromaLynx, TargetLynx, MassFragment, and MetaboLynx associated software (Waters) using a homemade database of more than 1400 substances.
2—Using LC-MS/MS (XEVO TQ-S system, Waters, Manchester, UK), another screening for psychoactive drugs including natural alkaloid compounds and related metabolites was performed in hair samples. After the addition of 50 µL of internal standards (methyl-clonazepam and β-OH-ethyltheophyllin) to 50 µL of the methanolic hair extract, the obtained mixture was subsequently evaporated to dryness at +30 °C under a gentle stream of nitrogen. The dry residue was reconstituted using 100 µL of ammonium formate buffer; 5 mM at pH 3: 10 µL was injected for separation using an Acquity UPLC HSS C18 column (150 × 2.1 mm, 1.8 µm, Waters) and a gradient including ammonium formate buffer 5 mM at pH 3 and acetonitrile in 0.1% formic acid. A Xevo TQ-S tandem mass spectrometer was used for detection after positive electrospray ionization mode in the MRM mode (using two transitions for each analyte) for more than 300 compounds and metabolites.
3—Lastly, the detection in hair of caffeine, theobromine, theophylline, nicotine, and cotinine was performed using the same LC-MS/MS device (XEVO TQ-S system, Waters, Manchester, UK). Calibration curves were extemporaneously achieved in the 0.1 to 20 ng/mg range for caffeine, theobromine, and theophylline, and in the 0.2 to 20 ng/mg range for nicotine and cotinine using a blank head hair sample. After the addition of 50 µL of internal standards (nicotine-D4 and cotinine-D3 at 5 µg/L) to 100 µL of the methanolic hair extract, the obtained mixture was subsequently evaporated to dryness at +30 °C under a gentle stream of nitrogen. The dry residue was reconstituted using 100 µL of ammonium formate buffer 5 mM at pH 3. Chromatographic separation of this extract was performed using an Acquity UPLC HSS C18 column (150 × 2.1 mm, 1.8 µm, Waters) and a gradient of (A) ammonium formate buffer 5 mM at pH 3 and (B) ACN/0.1% formic acid as mobile phase at a flow-rate of 0.4 mL/min during 5 min. The oven temperature was set at 50 °C and the injection volume was 10 µL. The gradient elution started with 100% of A and decreased to 50% at 3 min, and to 5% at 3.3 min. The washing step with 95% of solution B was held from 3.3 min to 4 min and the initial condition was applied from 4.1 min to 5 min. A Xevo TQ-S tandem mass spectrometer was used for detection after positive electrospray ionization mode in the MRM mode using the following transitions:
m/
z 195.2 to 138.0 (for quantitation) and
m/
z 195.2 to 110.0 for caffeine;
m/
z 181.1 to 68.9 (for quantitation) and
m/
z 181.1 to 95.9 for theobromine;
m/
z 181.1 to 124.0 (for quantitation) and
m/
z 181.1 to 68.9 for theophylline;
m/
z 163.2 to 132.0 (for quantitation) and
m/
z 163.2 to 117.0 for nicotine;
m/
z 177.1 to 79.9 (for quantitation) and
m/
z 177.1 to 98.0 for cotinine;
m/
z 167.2 to 136.0 (for quantitation) and
m/
z 167.2 to 134.1 for nicotine-D4; and
m/
z 180.1 to 80.0 (for quantitation) and
m/
z 180.1 to 100.9 for cotinine-D3. The validation procedure applied for this analytical method complied with both the French Analytical Toxicology Society (SFTA) and international recommendations for the validation of new analytical methods [
21,
22] including linearity, limit of detection (LOD), lower limit of quantification (LLOQ), accuracy, precision, and matrix effects. Calibration curves, estimated using 1/x weighted quadratic regression, were considered acceptable if the coefficient of determination (r
2) was at least 0.99. The LLOQ was the lowest concentration with the two transitions’ presence and an intra-assay precision CV% and a relative bias lower than 25%. Intra-day and inter-day precision were calculated by analysing 3 concentration levels in five replicates on five different days. Relative standard deviation (RSD) and percentage deviation of the average concentration from the corresponding nominal value were used to estimate precision and accuracy, respectively. These parameters were considered acceptable when they were lower than 25% at the LLOQ, and lower than 20% for other levels. Ion suppression phenomenon was studied following the experimental system previously proposed [
23]. Briefly, a standard solution containing the compounds of interest (at 100 µg/L) was continuously and directly infused into the mass spectrometer interface. A simultaneous LC flow containing either a pure mobile phase or a blank biological extract (blank hair) was introduced through a tee. Evolution of the signal of the transitions at the retention times of the corresponding compounds of interest was studied to determine the presence and intensity of ion suppression.
Chemical and reagents β-OH-ethyltheophyllin, methyl-clonazepam, 5-sulfosalicylic acid, ammonium formate, and formic acid, were purchased from Sigma-Aldrich (Saint-Quentin-Fallavier, France). LC-MS grade water and acetonitrile were purchased from Biosolve (Dieuze, France), while acetonitrile and methanol HPLC grade and 30% hydrochloric acid were purchased from VWR Prolabo (Fontenay-sous-Bois, France).
2.4. Lungs and Left Breast Sample Preparation
From our sample of 150 bodies, we were able to study lung histology from 14 subjects. In 7 out of these 14 cases, we checked hair and lungs concomitantly to establish whether high nicotine use was associated with anthracosis. For one woman (shamanic tree, #2), we endeavoured to ascertain if she smoked while breastfeeding because she presented high cotinine and nicotine levels, and one of the children—deceased between one and three years—that we excavated in the same grave was her son [
15] (whose hair was not preserved).
Lungs were rehydrated in Ruffer I solution for 3 h and then immersed in 10% neutral buffered formalin for 7 days. After fixation, tissues were dehydrated through increasing grades of ethylic alcohol, cleared in xylene and embedded in paraffin wax. Tissue sections of 5 μm thickness were stained with hematoxylin and eosin (H&E) to assess the general morphology.
Left breast: Tissues were naturally mummified by freezing and remained in this state until they were removed from the funeral chamber and autopsied. Upon reception in the laboratory tissues were rehydrated, given their very dry aspect, before the histological processing. Tissues were rehydrated in Ruffer I solution. Rehydration was completed after 3 h.
Once rehydrated, tissues were fixed in 10% buffered formalin for 4 days. One third of the tissues was deep frozen and the remaining pieces dehydrated through increasing grades of ethylic alcohol, cleared in xylene, and embedded in paraffin wax. Serial sections of 4 μm in thickness were cut and stained with: (i) Hematoxylin and Eosin (H&E) to assess the general morphology on paraffin sections, (ii) Giemsa stain for the demonstration of parasites on paraffin sections, and (iii) Sudan III and IV for the demonstration of lipids on frozen sections.
2.5. Relationship between Methylxanthines and Substances
To compare our results with historical data, we conducted an ethnopharmacological and historical study to determine what substances could have been consumed by the Yakuts at different periods. We then searched the literature in English and Russian for the composition of the various substances that could be the natural source of the alkaloid compounds of interest in this study. Our search strategies used a combination of standardized terms related to vegetal material (e.g., plant, tree, fruit, flowers) or the exact name of the vegetal, and keywords that we implemented in (i) NCBI PubMed (1900–present) and Google Scholar (1900–present); (ii) the data library of the N.N. Vorozhtsov Novosibirsk Institute of Organic Chemistry; and (iii) old monographs and articles of Belarus where alkaloids were very much studied but not available on the Internet. We also consider publications that were not found in the literature search but cited in retrieved publications. Unfortunately, the absence of tests on living organisms or significant bibliographic data concerning the relationships between the substances ingested and hair patterns of methylxanthines will oblige us to present these data in the form of hypotheses that we will compare with historical data.
2.6. Statistical Analyses
Principal Component Analysis was computed over the levels of the five pooled hair analytes, and outliers were defined qualitatively on the PCA projection.
2.6.1. Logistic Regression for Hair Concentration throughout Time
Multivariate logistic-regression analysis was used to assess a putative shift in the use of tobacco and tea with time. The dependent variable was binary coded: 0 for individuals dated before 1800 and 1 for individuals dated after 1800. The independent variables were the hair concentrations in the five studied xenobiotics and potential confounding factors that were adjusted for in the multivariable analyses included sex, age (0 to 14, 15 to 29, 30 to 49, and >50), social status (suspected shaman or not), and presence of an associated pipe (absence, presence: simple one, special one, and imported one). Computations were made on data including and excluding outliers.
2.6.2. Association between the Number of Metabolites and Time
The number of tea components detected in hair varies from 0 to 3 methylxanthines (theophylline, theobromine, caffeine) depending on the subject. We tested the variations over time (before and after 1800) for three possible combinations: theobromine alone versus all others, theobromine + caffeine versus theobromine + caffeine + theophylline, and all three components versus nothing.
2.6.3. Correlations Nicotine/Cotinine among Age Classes throughout the Time Range
Correlation by age class: Spearman correlation coefficient was computed between nicotine and cotinine concentrations within each age class (0–15, 15–30, 30–50, and over 50 years old).
Co-correlation between age classes: Co-correlation between age classes was assessed between the children age class (0–15 y/o) and each one of the others. The computation was performed with data including and not including the 4 outliers (#14, 29, 31, 41). A one-sided Fischer Z test was conducted with the null hypothesis being that correlation is lesser in the child group than in adult group and the alternative being that data for the children display a lesser correlation coefficient. Zou’s confidence intervals [
24] were also computed. Computations were made using the Cocor R package [
25].
2.6.4. Geographic Variation of Xenobiotic Hair Concentrations and Pipes
The association between xenobiotic hair concentrations and geography was assessed using several permutations-based methods. First, a permutation ANOVA was performed using the four archeological fields (North, East, West, and South) as groups [
26]. We then assessed correlation (Spearman’s rank test) between analyte concentrations and distances from the three main trading posts and/or fairs in the east through which the product was likely to transit and then the known entry areas (North-East and East) through which it was likely to arrive if it had been shipped to one of the ports on the Pacific coast (Okhotsk or Ayan) or the west, south-east and north-west, routes of arrival in Yakutia of products from the southern trading posts (on the Mongolian or Chinese border, Kiakhta or Nerchinsk) and/or from Russia. Strength of association was tested using a permutation over individual concentration.
As the associations between the levels of cotinine and/or nicotine with the pipes are statistically significant (9/10), we implemented the same methods of association on the totality of the pipes (for a majority of subjects the hair was not preserved), before and after 1800 AD. We present only the one before 1800 AD (20 pipes), the only positive one.
2.6.5. Interpolation & Mapping of PCA by Product
Tea and tobacco use were then approximated by combining the correspondent chemical xenobiotics by PCA (methylxanthines for tea, nicotine and cotinine for tobacco). Resulting PCA first axis values were interpolated over a rhomboid area comprised between the four studied archaeological fields using spline interpolation and triangulation based on Renkas tripack [
27] for linear interpolation as implemented in the “akima” R package. The interpolated values were overlaid according to their geographic coordinates over the map as black transparency (values normalized between 0 and 0.7). Maps were produced using the “marmap” R package and the National Oceanic and Atmospheric Administration bathymetrics and global relief data (resolution of 5 min) and then drawn following an orthographic projection for 100° longitude, 60° latitude, and 0 m elevation.
4. Discussion
Our study is the first to analyse many hair samples collected in the minutes or hours following the excavation of frozen bodies. Considering that the growth rate of human scalp hair ranges from 0.7 to 1.4 cm/month with an average of about 1 cm/month [
28] we analysed samples of hair corresponding to approximately two to three months’ growth before the subject’s death, and much less for the children for whom the hair was not very long.
Classically, researchers find that xenobiotics are well preserved in hair [
29], and several associations suggest the robustness of our results: (i) The ten subjects buried with a pipe had traces and/or significant hair concentrations of cotinine and/or nicotine, and one had no trace of nicotine and traces of cotinine. However, in this last case, it was an imported European earthenware pipe—the only one known in Yakutia—which indicates that the pipe was more of a prestige object than an object for smoking (
Figure 7). (ii) The nicotine/cotinine correlation increased with age and the correlation was significantly different between the young (0–15 years old) and the oldest quartile (over 50 years old). One of the two subjects with the lowest cotinine level was a child aged 6 to 12 months, and for four subjects where cotinine was detected alone, these levels were low. Present-day studies provide subjects with low levels of cotinine not associated with nicotine [
30]. (iii) An older woman from the end of the 18th century was the only subject with a level of nicotine associated with the highest cotinine level; she died with anthracosis pigments associated with emphysema. (iv) In graves dated before 1800 AD, we found a sub-significant association between the presence of a pipe and the distance to Yakutsk or the distance to the southwest of the studied area.
The measurement of nicotine hair concentration can be an informative tool to assess tobacco smoke exposure (whether active or passive) and related illnesses. The main pitfall of nicotine hair concentration interpretation is the impossibility to discriminate between active and passive exposure [
31,
32,
33]. Indeed, a positive result for nicotine in hair (usually greater than 0.2 ng/mg) could not confirm with any certainty an active use of tobacco by the deceased during their lifetime. Nicotine hair contamination from external sources (such as smoking by other members of the communities in which the deceased individuals lived) cannot be ruled out. Several authors reported the inability of this method to differentiate between active and passive exposure [
31,
32,
33]. In fact, observed hair concentrations for smokers range usually between 0.4 and 11 ng/mg and between 0.9 and 11 ng/mg for nicotine and cotinine, respectively. However, comparable hair values are reported in non-smoker populations, especially in the case of passive exposure [
34]. All in all, the determination of nicotine (and/or cotinine) in mummy hair samples cannot definitively confirm active nicotine use during the lifetime of the deceased, as external contamination and exposure to second-hand smoke during their lifetime cannot be excluded [
35].
Evolution over time of tobacco smoking: The only individual of the 17th century (Djoussulen) presented only a few traces of cotinine and nicotine in hair, perhaps reflecting passive exposure due to contact with Europeans as some of his grave goods demonstrate. Less than a century after contact, proximity to the major tobacco import trading post of Kiakhta was a determining factor in consumption, as were direct relationships of individuals with this spot. The development of economic circuits affected consumption, which largely flouted laws—as early as 1634 AD, smoking was prohibited in Russia on punishment of death (1634), and it remained so for autochthonous people until 1822 AD [
36]. The development of smoking was linked to the economic circuits in the 18th century, while in the 19th century tobacco consumption became generalized and economic circuits were detectable only for tea which was imported via the Pacific harbors. Between 1700 AD and 1800 AD, two men of the elite (#29 and #30) with a Mongolian haircut and a southern fashion for their clothes were heavy consumers of tobacco. One of them was found with a beautiful mammoth ivory pipe (
Figure 7). The manufacturing of highly attractive smoking accessories was instrumental in the dissemination of tobacco. The custom then was to pass the pipe from one person to the other to smoke thus encouraging family transmission. These men probably travelled to a trading post on the Mongolian border in contact with China where tobacco, tea, and furs were traded [
5]. Two other contemporary heavy smokers, females in this case (#12 and #26), were related to one of them and one was breastfeeding when she died. The mother of #29 and #30, who was not supposed to be living under the same roof as them when she died, was not a heavy tobacco user and drank undoubtedly a different tea than her daughter (see below). Similar family habits were also observed in the 19th century.
Currently, in most individuals, tea and caffeine are present at very high concentrations compared to the two other components [
37]. Theophylline and theobromine hair concentrations range between 0.1 and 0.6 ng/mg and between 0.3 and 10 ng/mg, respectively [
37]. In case of a positive result for caffeine, theophylline, or theobromine in hair (usually greater than 0.1 ng/mg), there is no data in the literature allowing us to assess the level of tea drinking by the subject (even if our sample the teapot was associated with the highest tea consumer). Tea contains some 3-methylxanthines (caffeine, theophylline, and theobromine), but its chemical composition is significantly affected by tea processing [
38], leaf maturation, botanical variety, geographical origin, agricultural practice [
39], and preparation [
40]. The influence of physiological variations from one subject to another on the incorporation of these methylxanthines in hair has not, to the best of our knowledge, been described. While all three 3-methylxanthines are present in most black teas, theophylline is often absent in green teas [
41] and some black teas [
42] (
Table 9). Following these observations, we could hypothesise that black tea consumption, which appeared punctually before 1750 AD, became widespread in the 19th century, which is consistent with historical data [
43].
The type of drink has been assigned to each subject according to the levels of theobromine, theophylline, and caffeine measured in the hair. Herbal tea is associated with the presence of theobromine only. Green tea was assigned when both theobromine and caffeine were present. Coffee is assigned to subjects with caffeine only. Black tea is assigned to subjects with all three xenobiotics associated with tea consumption (theobromine, caffeine, theophylline).
The level of tobacco consumption has been defined by taking into account the average of the two compounds associated with tobacco (nicotine and cotinine), the orders of magnitude of the two measures being too different to analyse them separately. A distribution curve of these averages has made it possible to define four groups of consumption: low consumer (nicotine and cotinine under 200 ng/mg), medium consumer (between 200 and 501 ng/mg), high consumer (between 501 and 2500 ng/mg), and very high consumer (over 2500 ng/mg).
In addition to these imported substances, we identified 19 other local plants that may have been drunk or smoked during historical times in Yakutia. Only a few contain the methylxanthines involved. An herbal tea based on (or including) the roots of
Rhodiola rosea, which contains caffeine in small quantities [
44], is known in the Siberian people’s pharmacopoeia for its anti-fatigue, anti-ageing effect. Three senior women had low levels of caffeine, and we excavated two of them in Villuy, one of the only places in Yakutia where it grows [
45]. Ivan tea, a fermented herbal tea from
Chamerion latifolium and
Chamerion angustifolium, which was drunk at that time, contains only theobromine in large quantities [
45]. It could have been consumed by five subjects of various ages from the 17th and 18th centuries; in the 19th century, we only suspected Ivan tea in three subjects, two of whom were less than one year old (
Table 9). Two subjects had only traces of theophylline; three senior women had only traces of caffeine, and three had traces of both caffeine and theophylline. This association could be related to herbal teas that we could not identify or passive exposure to tea. In the summer, Yakuts boil water before drinking it, and the container used to boil the water may have been used as a teapot beforehand.
It is highly challenging to give hypotheses about the substances’ origins because several substances could be drunk or mixed. Thus, two individuals from the 18th-century with a high level of theobromine compared to theophylline and caffeine levels in their hair were able to consume green tea and/or black tea and at the same time as Ivan tea. For example, the remains of tea leaf from the teapot tested positive for caffeine (+++), theophylline (++), and theobromine (+), whereas the hair of the woman buried with the teapot tested positive for these three substances but with a concentration of theobromine higher than that of theophylline. The tea preparation for the funeral meal may have been different from that consumed during the last months of the deceased’s life.
All periods combined, we found a few heavy tobacco consumers, men and women, amongst light or passive consumers, but tobacco-related co-morbidity was only recorded after 1750 AD (in a woman) at a time when a large part of the population, even older children, became addicted according to travelers’ accounts [
46,
47]. Heavy tea users were only detected in the 19th century and the heaviest users of the two substances dated from this century. In the 19th century, while the Yakut population assimilated into the Russian orthodox way of life, several contemporary trends were detected: all subjects seemed to become at least light consumers of tea, and moderate tea users were also tobacco users. The generalization of tea and tobacco is linked to economic circuits (smoking in the 18th century and tea in the 19th century). In our study area, tobacco use spread faster than tea use, probably because tobacco is intrinsically more addictive than tea, and because of the absence of a previous tradition in the use of these substances. A variety of herbal teas—especially Ivan tea—have been used by the Yakuts for centuries and it is not impossible that we have detected it. The use of green tea by women in the 18th century also supports this hypothesis, and the generalization of the use of tea follows the development of tuberculosis [
11] for which green tea was used as a treatment, perhaps with some results [
48]. In the 19th century, tea had made the transition from a beverage used in medicine and spirituality-based socialization (shamans?) to being used by elites with the tea ceremony led by women of fashion and means. The ceremony dictating that tea must be brewed and drunk using specialized equipment soon became established as a more generalized practice [
44]. For tobacco, there were no previous traditions; socio-economic and family contacts seem to have played a decisive role very early with this substance, which is more addictive than tea. In contemporary populations, numerous epidemiologic studies have consistently demonstrated positive associations between the use of tobacco, alcohol, and caffeine [
49], and it would be interesting to determine whether the same association existed among the Yakuts. To address this question, we have only historical records. The popularity of alcohol seems to have been delayed until the early 19th century [
47,
50]. It was banned from sale to the local population in 1822 [
51]. However, similar to tobacco before it, it continued to be illicitly traded, and the authors note the increasing consumption of tea and alcohol in the 19th century [
52,
53].