Introduction

The identification of a meta-organismal pathway involving dietary intake, gut microbiota and liver metabolism has raised much interest in the cardiometabolic field of research1. The group of Stanley Hazen provided strong support that metabolism of dietary phosphatidylcholine and L-carnitine by intestinal microbiota, resulting in the formation of the metabolite trimethylamine and its hepatic conversion to trimethylamine N-oxide (TMAO), induces atherosclerosis and that high TMAO levels predict an increased risk of cardiovascular disease (CVD)2,3,4,5. Another group also found elevated TMAO concentrations to independently predict coronary atherosclerosis and mortality in patients with chronic kidney disease6. However, although TMAO levels were found to be increased in hyperglycemic subjects and patients with impaired kidney function in other studies, they were not associated with all-cause mortality, cardiovascular death, or hospitalizations7 or history of myocardial infarction, the angiographically assessed presence of coronary heart disease or incident cardiovascular events during 8 years of follow-up8.

Importantly, Hazen’s group and others found the flavin mono-oxygenase 3 (FMO3) to play a crucial role in the regulation of circulating blood TMAO levels and in the promotion of atherosclerosis and cardiovascular disease in mice9,10,11. Interestingly, very recent studies suggest that FMO3 may also directly and independently of TMAO levels, induce hyperglycemia, hyperlipidemia and atherosclerosis in mice12,13. Given that hepatic FMO3 expression is under the negative control of insulin, hepatic FMO3 gene expression is elevated in obese and insulin resistant mice and, to a lesser extent in obese and patients with diabetes13 and fasting glucose levels were found to correlate positively with TMAO levels in humans2, the question arises to what extent circulating TMAO levels associate with precisely measured insulin sensitivity in humans. Furthermore, to better understand the impact of TMAO on CVD risk it should be clarified whether circulating TMAO levels associate with cardiovascular risk independently of insulin resistance. Finally, it has not been established whether a lifestyle modification, which is considered an important component in the prevention of CVD14,15, can reduce increased circulating fasting TMAO levels in humans.

To address these questions we quantified fasting serum TMAO levels in subjects without manifest CVD or kidney disease, which enables to study physiological and pathophysiological relationships under a minimal level of confounding factors. These subjects were at risk for type 2 diabetes, underwent precise measurements of glucose and lipid metabolism and measurement of the carotid intima-media thickness (cIMT) as an early marker of atherosclerosis and participated in a lifestyle intervention trial.

Materials and Methods

Subjects

Caucasians from the southern part of Germany participated in the Tübingen Lifestyle Intervention Program (TULIP)16,17. Individuals were included into the study when they fulfilled at least one of the following criteria: a family history of type 2 diabetes, a BMI >27 kg/m2, previous diagnosis of impaired glucose tolerance or gestational diabetes. Subjects were considered healthy according to a physical examination and routine laboratory tests. As assessed by means of a standard questionnaire, the participants had no history of liver disease such as hepatitis and did not consume more than >21 drinks (men) and >14 drinks (women) per week, cut-offs that are being considered to represent significant alcohol consumption when evaluating patients with suspected nonalcoholic fatty liver disease (NAFLD)18. A total of 220 subjects had precise measurements of body fat distribution, liver fat content and type 2 diabetes biomarkers at baseline and after 9 months of the lifestyle intervention. To investigate determinants of hepatic FMO3 expression and its relationship with insulin resistance we also studied another group of 55 subjects who underwent liver surgery and donated liver samples for research purposes19. Informed written consent was obtained from all participants and the Ethics Committee of the University of Tübingen had approved the protocol. The methods were carried out in accordance with the approved guidelines.

Lifestyle intervention

After the baseline measurements, individuals underwent dietary counselling and had up to ten sessions with a dietician. Subjects presented a three day food intake protocol to each visit. By individual counselling the aim was to reduce the intake of fat to amounts lower than 30% of total calorie intake, to reduce the intake of saturated fat to amounts lower than 10% of total fat intake, to increase the intake of fibres to 15 g/1000 kcal and to achieve a weight loss >5% during the study. Diet composition was estimated with a validated computer program using two representative days of a 3 day diary (DGE-PC 3.0; Deutsche Gesellschaft für Ernährung, Bonn, Germany). Individuals were asked to perform at least 3 hr of moderate sports per week. Aerobic endurance exercise (e.g. walking, swimming) with only moderate increase in the heart rate was encouraged. The anthropometric, metabolic, physical activity and energy intake measurements were repeated after 9 months.

Habitual physical activity

All individuals completed a standardized self-administered and validated questionnaire to measure physical activity and a habitual physical activity (HPA) score was calculated20.

Total body fat mass, body fat distribution and liver fat content

Measurements of total body- and visceral fat and mass were performed by an axial T1-weighted fast spin echo technique with a 1.5 T whole-body magnetic resonance imager. Liver fat content was measured by localized proton magnetic resonance (1H-MR) spectroscopy as previously described21. NAFLD was defined as liver fat content >5.56%22.

Oral glucose tolerance test

Subjects underwent a frequently sampled 2 hr, 75 g oral glucose tolerance test (OGTT). Venous plasma samples were obtained at 0, 30, 60, 90 and 120 minutes for determination of plasma glucose and insulin levels. Whole body insulin sensitivity was calculated from glucose and insulin values during the OGTT as proposed by Matsuda and DeFronzo23 or by using the homeostatic model assessment of insulin resistance (HOMA-IR)24 in the 55 subjects from whom liver samples were available. Insulin clearance was estimated from the OGTT as C-peptide AUC/insulin AUC.

Euglycemic, hyperinsulinemic clamp

In a subgroup of 167 subjects insulin sensitivity was also determined during a euglycemic, hyperinsulinemic clamp with a primed insulin infusion at a rate of 40 mU · m−2 · min−1 for 2 hours. The insulin sensitivity index measured during the clamp (in μmol·kg−1 · min−1 · pM−1) was calculated as the mean infusion rate of glucose (in μmol · kg−1 · min−1) necessary to maintain euglycemia during the last 40 minutes of the euglycemic hyperinsulinemic clamp divided by the steady state plasma insulin concentration.

Carotid intima-media thickness

The cIMT was measured in the fasting state using a high-resolution ultrasound system (AU5 idea, Esaote Biomedica, Munich, Germany) with an integrated electrocardiography (ECG) package as previously described25.

Liver samples

A total of 55 patients (20 females/35 males, age 62 ± 11 years, bodyweight 76 ± 13 kg, BMI 25.2 ± 3.9 kg ∙ m−2) who underwent liver surgery (e.g. for the resection of solitary hepatic lesions) in the Department of General, Visceral and Transplant Surgery at the University of Tübingen were included in the present study. Patients were fasted overnight prior to collection of the liver samples. Patients tested negative for viral hepatitis and had no liver cirrhosis. Liver samples were taken from normal, non-diseased tissue during surgery, immediately frozen in liquid nitrogen and stored at −80 °C.

RNA isolation, RT-PCR and real-time quantitative PCR analysis of hepatic mRNA expression

Frozen tissue was homogenized in a TissueLyser (Qiagen, Hilden, Germany) and RNA was extracted with the RNeasy Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Total RNA treated with RNase-free DNase I was transcribed into cDNA using first-strand cDNA kit and PCRs were performed in duplicates on a LightCycler480 (Roche Diagnostics, Mannheim, Germany). Data are presented relative to the housekeeping gene Rps13 using the ΔΔCt method. The human primer sequences used are: for RPS13 up 5′-CCCCACTTGGTTGAAGTTGA-3′ and down 5′-ACACCATGTGAATCTCTCAGGA-3′ and for FMO3 5′-cctgctttgagaagagcaatg-3′ and down 5′- tggaaaagactgatttgtaaatgct-3′.

Determination of liver tissue triglyceride content

Tissue samples were homogenized in phosphate-buffered-saline containing 1% Triton X-100 with a TissueLyser (Qiagen, Hilden, Germany). To determine the liver fat content, triglyceride and total protein concentration in the homogenate was quantified using the ADVIA 1800 clinical chemistry analyzer (Siemens Healthcare Diagnostics, Eschborn, Germany) and the results were calculated as mg/100 mg tissue (%) or relative to protein content respectively.

Analytical procedures

Blood glucose was determined using a bedside glucose analyzer (glucose-oxidase method; YSI, Yellow Springs Instruments, Yellow Springs, OH). Plasma insulin was determined on an ADVIA Centaur XP and all other blood parameters on an ADVIA 1800 clinical chemistry system (Siemens Healthcare systems, Erlangen, Germany). The estimated glomerular filtration rate (eGFR) was quantified by using the MDRD-4 equation26.

TMAO was measured in serum samples that were stored at −80 °C. Stability studies revealed that TMAO is stable under these storage conditions for several years27. For the quantification of serum TMAO levels from samples first sample preparation was performed. The extraction solvent consisted of methanol (MeOH)/acetonitirile (CAN) (v/v, 1/1) containing d9-TMAO (1 μM) and d3-carnitine (20 μM). 200 μL extraction solvent was added to 50 μL serum, vortexed for 1 min and centrifuged at 13,000 g for 10 min. Then 100 μL supernatant was freeze-dried and before analysis redissolved in 1 mL 70% ACN, sonicated for 10 min, then transferred to a sample vial for ultra-high performance liquid chromatography (UHPLC)-MS/MS analysis. The chromatographic separation was performed with a Nexera X2 UHPLC system (Shimadzu, Japan) equipped with a silica hydride column (150 × 2.1 mm, 4 μm, MicroSolv Technology Corporation, USA). The following settings were used: oven temperature 40 °C, flow rate 0.4 mL/min, autosampler temperature 10 °C, injection volume 1 μL. Mobile phase A was 50% ACN containing 0.1% formic acid and 10 mM ammonium formate, mobile phase B was 90% ACN containing 0.1% formic acid. The gradient program started with 50% A, increased to 85% A within 3 min, held 85% A for 2 min, increased to 100% A in 3 min, held this for 2 min, finally decreased back to 50% A in 0.5 min for an equilibration of 2.5 min. Mass spectrometric analysis was performed using a triple quadrupole 8050 electrospray ionization (ESI)-MS (Shimadzu, Japan) in the positive ion mode. The ESI source parameters were: nebulizing gas flow 3 L/min, heating gas flow 10 L/min, interface temperature 300 °C, desolvation line temperature 250 °C, heat block temperature 400 °C, drying gas flow 10 L/min. Multiple Reaction Monitoring (MRM) was used. TMAO and carnitine was quantified by deuterated internal standards.

Statistical analyses

Data that were not normally distributed (Shapiro-Wilk W test; e.g. cIMT, insulin sensitivity, body fat mass, liver fat, content, blood biomarkers) were logarithmically transformed. Differences between baseline and follow-up were tested using the matched pairs t test. Univariate associations between parameters were tested using Pearson correlation analyses. To adjust the effects of covariates and identify independent relationships, stepwise and multivariate linear regression analyses were used. The statistical software package JMP 11.0 (SAS Institute Inc, Cary, NC, USA) was used.

Results

Subject characteristics

Anthropometrics and metabolic characteristics of the 220 subjects, who were phenotyped at baseline and at follow-up, are shown in the Table 1. A total of 90 males and 130 females with a mean age of 46 years and a mean BMI of 29.5 kg ∙ m−2 were studied. Applying ultra-high performance liquid chromatography mass spectrometry, as was done in previous studies2,4 mean ± SD baseline fasting levels of serum TMAO of 2.83 ± 1.6 μmol/l in our study were very similar in magnitude as was reported in those studies. There was a large variability in serum TMAO levels among the subjects ranging from 0.77 to 11.51 μmol/l.

Table 1 Anthropometric and metabolic parameters at baseline and at follow-up.

Cross-sectional relationships at baseline

At baseline circulating TMAO levels were not different between males (2.82 ± 1.60 μmol/l) and females (2.85 ± 1.64 μmol/l, p = 0.89). They correlated positively with age, BMI, total- and LDL- cholesterol, apolipoprotein B levels and tumor necrosis factor (TNF) alpha levels (N = 144), but not with total body fat mass, HDL-cholesterol or high sensitivity C-reactive protein (hs-CRP) levels (Table 2).

Table 2 Relationships of cIMT and fasting TMAO levels with selected parameters at baseline.

Of particular interest were possible relationships of circulating TMAO levels with renal function, body fat distribution, liver fat content, glycemia, insulin sensitivity and cIMT. TMAO levels correlated negatively with the eGFR, tended to associate positively with visceral fat mass and did not associate with liver fat content. While a positive correlation of TMAO levels was found with fasting glycemia, no relationship was observed with 2 hr glucose levels and insulin sensitivity, measured during the OGTT or the euglycemic, hyperinsulinemic clamp (Table 2). There was also no relationship of TMAO levels with insulin clearance. Furthermore, after adjustment for age, sex and BMI no relationship of TMAO levels with insulin sensitivity (OGTT, p = 0.25; clamp, p = 0.69) was found. However, TMAO levels correlated positively with cIMT (Fig. 1).

Figure 1
figure 1

Cross-sectional relationship of fasting serum TMAO levels with cIMT.

Univariate cross-sectional relationship (Pearson correlation coefficient and p-value) of fasting serum trimethylamine N-Oxide (TMAO) levels with carotid intima-media thickness (cIMT) in 220 subjects.

When we investigated possible relationships of cIMT with anthropometrics and metabolic characteristics we found that a somewhat other picture emerged compared to the relationships of TMAO levels with these parameters. While similar to TMAO levels, cIMT correlated with age, BMI, fasting glucose levels, total- and LDL-cholestreol levels and apolipoprotein B levels, cIMT also correlated with sex (women: 0.54 ± 0.11; men: 0.60 ± 0.13 mm, p = 0.0009), visceral fat mass, liver fat content and insulin sensitivity (Table 2).

We next investigated whether higher TMAO levels associated with increased cIMT independently of established determinants of CVD risk. For this we first ran a forward stepwise linear regression analysis including all parameters that correlated with cIMT. In this analysis age was found to be the strongest determinant of cIMT (p < 0.0001) and was followed by visceral fat mass (p = 0.0001) and TMAO levels (p = 0.03), while all the other parameters were not independent determinants of cIMT. In a multivariate linear regression model including age, visceral fat mass and additionally sex, which strongly associated with cIMT, higher TMAO associated with elevated cIMT, independently of these parameters (Table 3).

Table 3 Determinants of cIMT in multivariate linear regression models.

Longitudinal relationships

During the lifestyle intervention most anthropometric and metabolic parameters, which are associated with CVD risk, as well as cIMT improved. Interestingly, mean TMAO levels did not change (Table 1). Mean total energy (from 2014 ± 527 kcal to 1855 ± 393 kcal, p < 0.0001), total fat (from 33 ± 6% to 31 ± 5%, p < 0.0001) and saturated fat (percentage of total fat; from 44 ± 6% to 41 ± 4%, p < 0.0001) intake decreased, yet the changes of these parameters did not correlate with the change of TMAO levels (all p values ≥0.31). Furthermore, no significant relationship was observed between the change of TMAO levels with the change of cIMT in the whole population. However, when we stratified subjects into tertiles of change in TMAO levels (Table 4) we found that cIMT decreased significantly only in subjects in the tertile with the largest decrease of TMAO levels (>20%) (Table 4 and Fig. 2). Interestingly, while in the other two tertiles most parameters considerably improved (e.g. liver fat content, insulin sensitivity) such changes were weaker or absent in the tertile with the largest decrease in TMAO levels (Table 4).

Table 4 Anthropometric and metabolic parameters in tertiles of changes of TMAO levels during the lifestyle intervention.
Figure 2
figure 2

Change of TMAO levels and change of cIMT during 9 months of lifestyle intervention.

Individual carotid intima-media thickness (cIMT) before and after 9 months of lifestyle intervention according to tertiles of changes in fasting serum trimethylamine N-Oxide (TMAO) levels (p-value from the paired t test) in 220 subjects.

FMO3 mRNA expression in liver samples

Considering that TMAO levels are not only influenced by the diet but also by FMO3 activity and that hepatic FMO3 expression is under the control of insulin signalling, we also investigated whether hepatic FMO3 mRNA expression associated with anthropometric and metabolic parameters. In 55 liver samples FMO3 mRNA expression levels were higher in females than in males (Fig. 3, panel A). FMO3 mRNA expression was not associated with BMI or eGFR (both p > 0.45). However, there was a negative correlation between hepatic FMO3 mRNA expression and age (Fig. 3, panel B). Liver fat content as well as insulin resistance, measured by HOMA-IR, were not significantly associated with hepatic FMO3 mRNA expression (both p > 0.26). After adjustment for age and gender, HOMA-IR tended to positively associate with hepatic FMO3 mRNA expression (Fig. 3, panel C).

Figure 3
figure 3

Relationships of hepatic FMO3 mRNA expression with sex, age and adjusted HOMA-IR.

Relationship of hepatic FMO3 mRNA expression (normalized for the mRNA expression of the housekeeping gene RSP13) with gender (A) and age (B) and of hepatic FMO3 mRNA expression, adjusted for sex and gender, with the homeostatic model assessment of insulin resistance (HOMA-IR) in 55 patients who donated liver samples (C).

Discussion

Because CVD is the leading cause of deaths in industrialized countries28 better prediction and prevention of the disease is of major importance for the individual and the society. In the search for the major pathomechanisms promoting CVD the identification of the meta-organismal pathway involving dietary intake, gut microbiota and liver metabolism has raised much interest in the field of CVD research. Particularly increased generation of TMAO is thought to represent a major pathomechanism contributing to increased CVD risk1. While in vitro and animal data support that TMAO may play an important role in the pathogenesis of atherosclerosis29,30, in humans such a role has only been estimated from correlation studies2,3,4,5. Recently other studies could not confirm that elevated fasting blood TMAO levels predict increased CVD risk7,8. It was hypothesized that a positive relationship of TMAO levels with CVD risk may in part be confounded by impaired kidney function and poor metabolic control8. Furthermore, because hepatic FMO3, the enzyme that is involved in the generation of TMAO, is under strong control of insulin signalling13 and because TMAO is thought to regulate insulin signalling in mice31, insulin sensitivity may represent another important confounding factor in such a relationship.

To address these issues we measured fasting serum TMAO levels in our study population that was phenotyped for glucose and lipid metabolism, body fat distribution, ectopic lipid deposition in the liver and cIMT, an early marker of atherosclerosis. In agreement with other studies2,8 we found TMAO levels to correlate positively with age, BMI, fasting glycemia and blood lipids. Furthermore, although TMAO levels did not associate with hs-CRP levels, there was a relative strong and positive correlation of TMAO levels with TNF-alpha levels, a finding that may warrant further investigation. TMAO levels also correlated positively with cIMT. Interestingly, while insulin sensitivity also associated with cIMT, no relationship was found between TMAO levels and insulin sensitivity, measured by the OGTT or the clamp.

To our knowledge this is the first study investigating possible relationships of TMAO levels with precisely measured insulin sensitivity and early atherosclerosis in a large group of subjects. Our findings support the notion that TMAO levels predict CVD risk independently of insulin sensitivity, a major pathomechanism of cardiometabolic diseases32. Furthermore, we provide information about possible effects of insulin sensitivity on endogenous production of TMAO. Studies in insulin receptor knock-out mice or in patients with defects in insulin signalling can help to understand the impact of insulin signalling on gene expression33. In this respect recently Miao et al. showed that in Liver Insulin Receptor Knockout (LIRKO) mice FMO3 was the second most highly upregulated hepatic transcript and that TMAO was the most strongly upregulated metabolite in these animals, compared to controls13. Furthermore, in that study insulin was found to suppress FMO3 mRNA expression in primary hepatocytes and it was shown that FMO3 knockdown improved glucose tolerance in mice. Finally, hepatic FMO3 mRNA expression was found to be increased in females and in heavily obese and hyperglycemic patients, compared to leaner and normoglycemic subjects13. In agreement with Miao et al.13 in our normal weight patients who donated liver samples we also found hepatic FMO3 mRNA expression to be elevated in females. Furthermore, we found FMO3 expression to negatively associate with age and that FMO3 expression tended to associate positively with adjusted insulin resistance. These data and the fact that in our larger group of 220 subjects TMAO levels correlated positively with fasting glycemia, an estimate of hepatic insulin resistance, indicates that also in humans fasting TMAO levels appear to be regulated to a considerable extent by hepatic FMO3. However, because TMAO levels did not correlate with whole-body insulin resistance TMAO does not appear to be a major regulator of glucose metabolism in humans.

In agreement with the findings from our cross-sectional analysis also in our longitudinal analyses the change of insulin sensitivity did not associate with the change of TMAO levels. Importantly, while insulin sensitivity and other cardiometabolic risk parameters improved, mean TMAO levels did not change. When we investigated whether the change of macronutrient intake and the intake of saturated fat, which is mainly derived from animal sources, correlates with the change of TMAO levels, no such relationship was found. Previous studies investigating the impact of diet modification on fasting TMAO levels revealed that a 2-week high-fat diet did not change fasting TMAO levels34, while a small increase of TMAO levels was found during a 4-week high-fat diet intervention35. Together with our data these findings indicate that fasting TMAO levels are not strongly regulated by a moderate modification of the diet that is currently being recommended for the prevention of cardiometabolic diseases36. However, we have no solid information about the dietary intake of choline or carnitine in our study and, therefore, cannot comment whether a specific modification of foods rich in choline or carnitine occurred in some of our subjects. Furthermore, we have no information about the gut microbial composition and function of our subjects. It is important to study these aspects in future studies, because Hazen’s group very recently elegantly showed in mice that the microbial TMA lyase activity is very important for the conversion of choline and carnitine into TMA37. Individual differences in the gut microbiome and the nutrient intake may have resulted in that we found no significant mean change of fasting TMAO levels during our intervention, yet a large variability in the change of TMAO levels, ranging from a decrease of 82% to an increase of 730%.

When we stratified our subjects into tertiles of the change of TMAO levels we found that cIMT decreased significantly only in subjects in the tertile with the largest decrease of TMAO levels. Interestingly, specifically in this tertile, except for TMAO levels, the improvement of many metabolic parameters (e.g. insulin sensitivity and glycemia), liver fat content and viscveral fat mass was smallest, indicating that the improvement of cIMT in this tertile was most probably not strongly affected by the improvement of these parameters. Otherwise, one may expect cIMT not to decrease, or even increase, in the tertile where TMAO levels increased. However, particularly in this tertile the improvement of insulin sensitivity, glycemia, liver fat content visceral fat mass and serum lipids was largest, which may have counteracted the effect of increasing TMAO levels on cIMT. Of note, most studies investigating the change of cIMT during lifestyle or pharmacological interventions, or during a natural follow-up, found that, based on low reversibility of arterial structure, cIMT did not decrease, but that the intervention or a healthy lifestyle resulted in a reduced progression of cIMT38,39. The fact that mean cIMT decreased during the intervention in our study indicates that lifestyle modification is effective to revert atherosclerosis, particularly in early stages of its pathophysiological process.

In summary, with this report we provide novel information, particularly from a metabolic angle, that may help to better understand the role of TMAO in the natural history of cardiometabolic diseases. Based on the relationships of fasting TMAO levels with cIMT in the cross-sectional and longitudinal analyses, our data support the conclusion that derived from published animal and human data, that TMAO levels have an important role in the pathophysiology of CVD. Furthermore, we found that fasting TMAO levels associated positively with fasting glycemia and that hepatic FMO3 mRNA expression tended to associate positively with the HOMA-IR. These findings indicate that, similar as shown in animals and in obese, hyperglycemic humans13 FMO3 expression and TMAO levels appear to be increased in hepatic insulin resistance. However and in contrast to other liver-derived proteins, such as the hepatokine fetuin-A40, TMAO does not appear to be a major regulator of whole-body insulin sensitivity. Finally, the fact that mean TMAO levels did not change during a standard lifestyle intervention, which is effective to improve most of the studied cardiometabolic risk parameteres, indicates that it is important to better understand mechanisms regulating TMAO blood levels.

In conclusion, we found that increased fasting serum TMAO levels associate with increased cIMT, independently of established cardiovascular risk markers, including insulin resistance, visceral obesity and fatty liver. Furthermore, the decrease of cIMT during a lifestyle intervention appears to be related to the decrease of TMAO levels. Our data support the idea that further research into this meta-organismal pathway may help to better understand the natural history of CVD.

Additional Information

How to cite this article: Randrianarisoa, E. et al. Relationship of Serum Trimethylamine N-Oxide (TMAO) Levels with early Atherosclerosis in Humans. Sci. Rep. 6, 26745; doi: 10.1038/srep26745 (2016).