Son and Daughter with Autism Analysis – Part 2

This is a follow up the earlier blog post: Son and Daughter with Autism Analysis from a year ago. There has been a lot of changes of the site and revisions of algorithms.

Comparing Siblings

We know from studies that members in the same family often share about 27% of the same strains. Unfortunately with 16s tests (Biomesight, Ombre), we do not get strain information just species information.

Using the new refactor citizen science symptoms (see New Special Studies on Symptoms ), we are presently surprised! We have many forecast symptoms being the same which supports the observation cited above of share taxa, likely at the strain level.

It does hint that less time with each other and a lot more time with other (ideally normal) children may have benefits to the microbiome. Some of the changes may be connected to gender:

  • About twice as many women as men experience depression [Mayo]
  • Increased inflammation is seen in the periphery in both depression and fatigue. [2019] which agrees with the daughter having a lower Anti inflammatory Bacteria Score
MeasureDaughterSon
Anti inflammatory Bacteria Score25.6%ile89%ile
Buytrate Bacteria Score95.9%ile78%ile
Histamine Producers21.8%ile15.3%
Autism From PubMed53/97 (1%ile)
Prior: 54/97 (1%ile)
73/97 (11%ile)
Prior: 53/97 (1%ile)
Forecast SymptomsOfficial Diagnosis: Depression
28 % match on 7 taxa

DePaul University Fatigue Questionnaire : Blurred Vision
25 % match on 8 taxa
Neurological-Sleep: Inability for deep (delta) sleep
23 % match on 13 taxa

Age: 10-20
17 % match on 23 taxa

DePaul University Fatigue Questionnaire : Forgetting what you are trying to say
16 % match on 31 taxa

Next looking at Percentages of Percentiles, we see significant differences. Unfortunately, we do not have gender and age reference tables, so interpretation is fuzzy.

Potential Medical Conditions Detected had nothing significant for either child. Both are at 95.6%ile on Dr. Jason Hawrelak Recommendations (they were 98.9 and 99.7%iles before) .

Detail Comparison

The thing that stands out is that the Son has a lot more Enzymes out of range (with the resulting substrates(consumers) and products also being out of range).

CriteriaDaughterSon
Lab Read Quality6.710.9
PubMed Bacteria Matches for Autism1%ile (53/97)11%ile (73/97)
Outside Range from JasonH55
Outside Range from Medivere1818
Outside Range from Metagenomics99
Outside Range from MyBioma99
Outside Range from Nirvana/CosmosId2323
Outside Range from XenoGene5151
Outside Lab Range (+/- 1.96SD)1122
Outside Box-Plot-Whiskers9876
Outside Kaltoft-Moldrup132248
Bacteria Reported By Lab842757
Bacteria Over 99%ile1311
Bacteria Over 95%ile2845
Bacteria Over 90%ile6661
Bacteria Under 10%ile66285
Bacteria Under 5%ile22214
Shannon Diversity Index3.0642.807
Simpson Diversity Index0.070.088
Chao1 Index2979124924
Rarely Seen 5%6198
Pathogens3635
Kegg Compounds Low9731001
Kegg Compounds High43162
Kegg Enzymes Low89265
Kegg Enzymes High98381
Kegg Products Low55152
Kegg Products High52209
Kegg Substrates Low46148
Kegg Substrates High58229

Looking at KEGG Derived Probiotic suggestions, the list is full of the soil based bacteria found in Prescript-Assist®/SBO Probiotic or Energybalance / ColoBiotica 28 Colon Support or General Biotics/Equilibrium. There was no probiotic above my usual threshold from the consensus, so the above seems to be worth a try.

KEGG Suggested supplements has nothing significant for the daughter, but for the son we have the following being very significant:

  • Serine
  • Threonine
  • Glutamine
  • Cysteine
  • Arginine

A complex amino-acid supplement may be worth an experiment.

As an experiment (and trying to avoid two different kid diet), I did an uber-consensus from each child’s with tons of prescription medication but only one thing above my usual 50% of highest value.

Son Compared to Prior Sample

We can see the spike in low percentile bacteria. This raises the question, has he had COVID (or a COVID vaccine) prior to the sample being done. These spikes show themselves also via Kaltoft-Moldrup and Box-Plot-Whiskers which are both sensitive to this pattern.

CriteriaCurrent SampleOld Sample
Lab Read Quality10.94.4
Outside Range from JasonH44
Outside Range from Medivere1717
Outside Range from Metagenomics1010
Outside Range from MyBioma1313
Outside Range from Nirvana/CosmosId2727
Outside Range from XenoGene4949
Outside Lab Range (+/- 1.96SD)2222
Outside Box-Plot-Whiskers76100
Outside Kaltoft-Moldrup24889
Bacteria Reported By Lab757708
Bacteria Over 90%ile6182
Bacteria Under 10%ile28526
Shannon Diversity Index2.8072.451
Simpson Diversity Index0.0880.15
Chao1 Index2492419183
Lab: Thryve
Pathogens3530
Condition Est. Over 90%ile20
Kegg Compounds Low10011048
Kegg Compounds High162132
Kegg Enzymes Low265115
Kegg Enzymes High381296
Kegg Products Low15274
Kegg Products High209191
Kegg Substrates Low14869
Kegg Substrates High229212
Anti inflammatory Bacteria Score89.2%ile83.2%ile
Buytrate Bacteria Score77.9%ile90.2%ile
Histamine Producers15.3%ile38.2%ile
Histamine dropping is usually a good sign

From this weekend update of special studies, we can get a count of bacteria shifts strongly associated to symptoms.

  • Old Sample: 32 taxa
  • Latest Sample: 60 taxa

Daughter Compared to Prior Sample

First the numbers which are usually similar to the prior sample.

CriteriaCurrent SampleOld Sample
Lab Read Quality6.73.1
Outside Range from JasonH66
Outside Range from Medivere1919
Outside Range from Metagenomics77
Outside Range from MyBioma1212
Outside Range from Nirvana/CosmosId2626
Outside Range from XenoGene4747
Outside Lab Range (+/- 1.96SD)1161
Outside Box-Plot-Whiskers98203
Outside Kaltoft-Moldrup132134
Bacteria Reported By Lab842852
Bacteria Over 90%ile66202
Bacteria Under 10%ile6610
Shannon Diversity Index3.0643.411
Simpson Diversity Index0.070.028
Chao1 Index2979135210
Lab: Thryve
Pathogens3635
Condition Est. Over 90%ile00
Kegg Compounds Low9731027
Kegg Compounds High4380
Kegg Enzymes Low8944
Kegg Enzymes High98171
Kegg Products Low5529
Kegg Products High5286
Kegg Substrates Low4626
Kegg Substrates High58111
Anti inflammatory Bacteria Score25.5%ile28%ile
Buytrate Bacteria Score95.9%ile74.5%ile
Histamine Producers21.7%ile28.7%ile

From this weekend update of special studies, we can get a count of bacteria shifts strongly associated to symptoms.

  • Old Sample 53
  • Latest Sample: 39

Out of curiosity, I compared the symptom associated outliers. We found 3 are matches (of these 39) and one not matches for the taxa reported for each. That is close to the expected percentage of the same strains for people in the same house.

Bacteria NameDaughterSon
  NegativicutesToo HighToo High
  AcidobacteriiaToo LowToo Low
  Bacteroides eggerthiiToo HighToo Low
  PropionibacterialesToo LowToo Low

Going Forward

Autism has challenges because of its complex nature. This is compounded by a low number of samples to work from for Citizen Science analysis. The shifts reported from PubMed have a high pattern match with people who do not have autism.

I am going to try building a consensus for each by doing two itemsL

  • “Just give Me Suggestions”
  • Doing PubMed Autism on [Changing Microbiome]
  • [All Bacteria identified by special studies]

The rationale is that the last one identify the bacteria that appears to be symptom causing in many people. We have a very poor match from what we do have a match for. This is not surprising because autism is a very wide spectrum.

We then see six sets of suggestion

Son

When I look at the details we have over 150 items with 6 recommended take (i.e. everyone agrees)

The probiotics that have no known adverse risk for any bacteria is below. The high value is 510.

Daughter

When I look at the details we have just 15 items with 6 recommended take (i.e. everyone agrees)

The probiotics that have no known adverse risk for any bacteria are low in computed benefit, so I would ignore them.

  • Their values are low: 16/31 out of a high value of 301

Bottom Line

The failure to find significant matching patterns is a bit of a frustration to me. What we did find had very good agree for the son with 150+ items having each of the size suggestion set agreeing for the take. For the daughter, it was not as strong: 15 for 6 sets being in agreement, and 50 with 5 sets being in agreement.

Questions

  1. I assume higher anti-inflammatory score is better – Daughter was 25% and Son was 89%
  2. Deep Sleep with Son – 23% match that he has deep sleep issues is pretty strong?
    • Does not jump out, but indicates that microbiome is playing a role.
  3. Son – lot more missed enzymes – what is that do you believe and probiotics help with that?
    • I avoid the word “believe”. A rational assumption is that disruption of enzymes compare to others impacts how the cells (including brain cells) react.
  4. Spike in low percentage bacteria – likely long Covid for Son means he has less good bacteria now?
    • I avoid the words “good” or “bad” bacteria. Any bacteria far enough from typical values become bad; disruption to the microbiome and the body. Theses spikes are typically seen (pattern matching) with two conditions: Long COVID and ME/CFS. A common symptom of these two issues are cognitive issues – for example: memory, ability to learn, etc.
  5. Histamine – Higher percentage is worse correct?  Daughter was 21% and son was 15% 
  6. Higher Butryate percentage is better?  95% daughter / 78% son

The Journey Begins with your microbiome

Thanks for joining me!

This is a companion site to the analysis site at: 

https://microbiomeprescription.com/

The intent of this site to assist people with health issues that are, or could be, microbiome connected. There are MANY conditions known to have the severity being a function of the microbiome dysfunction, including Autism, Alzheimer’s, Anxiety and Depression. See this list of studies from the US National Library of Medicine. Individual symptoms like brain fog, anxiety and depression have strong statistical association to the microbiome. A few of them are listed here.

The base rule of the site is to avoid speculation, keep to facts from published studies and to facts from statistical analysis(with the source data available for those wish to replicate the results). Internet hearsay is avoid like the plague it is.

Open data and Open source are our mottos!

Continue reading “The Journey Begins with your microbiome” →Posted on  by lassesenEdit

Biomesight #4 Sample: IBS and COVID

We have a varied history with some storms blowing us off courses. Here’s a list of the tests and prior blog posts:

His comments are short:

  • I would say some small subjective improvements since last time, but no major changes. Reminder: I have a friendly MD in terms of antibiotics.
  • Metronidazole was on top in the last samples, I did it back then.
    • Comment: Metronidazole is no longer at the top but dropped down to 16% of the highest value. It appears to have done its magic in reducing the bacteria pointing to it as a tool..

Base Analysis

When people have multiple samples, I like to do side-by-side comparisons, especially when someone has been doing some of the suggestions suggested. The suggestions are computed and may not always work. Expert Systems and AI are not perfect; they typically do better than a person with only a few years of experience that has training in the discipline (better consistency, remember more facts, etc). How are we doing objectively?

Scores

We see two positive shifts in the latest sample: Increase of Anti inflammatory Bacteria Score and decrease of Histamine Producers.

Percentages of Percentiles

We see a lot of bouncing around between samples. The middle two images matches the typical pattern seen with ME/CFS and Long COVID. Those shifts have faded over the last 3 months with a different pattern appearing indicating a different dialect of gut dysfunction.

Multi-Vector Comparison

The main numbers are below. The take away, less bacteria that are in the high percentile range (at 95%ile, 10 -> 28 -> 23 -> 8). The numbers bounce around with the middle two being similar and the other two also similar. There are no really clear shift in these measures.

Criteria11/18/20215/20/20226/22/20239/4/2023
Lab Read Quality8.15.54.77.2
Outside Range from JasonH6699
Outside Range from Medivere16161515
Outside Range from Metagenomics8877
Outside Range from MyBioma5566
Outside Range from Nirvana/CosmosId20202323
Outside Range from XenoGene29293535
Outside Lab Range (+/- 1.96SD)76173
Outside Box-Plot-Whiskers36695438
Outside Kaltoft-Moldrup93484788
Bacteria Reported By Lab652508542558
Bacteria Over 99%ile7462
Bacteria Over 95%ile1028238
Bacteria Over 90%ile29423622
Bacteria Under 10%ile2084150175
Bacteria Under 5%ile180198157
Shannon Diversity Index1.8531.8261.2721.556
Simpson Diversity Index0.0560.0380.0870.09
Rarely Seen 1%2271
Rarely Seen 5%145218
Pathogens41242936

From Special Studies

The top match was the same on all of the samples, with an increase when there was actually COVID.

Criteria11/18/20215/20/20226/22/20239/4/2023
COVID19 (Long Hauler)28%ile33%ile41%ile28%ile
Next one:15%ile26%ile20%ile13%ile

The “next one” dropping implies some reduction of dysbiosis

Health Analysis

Using Dr. Jason Hawrelak Recommendations, there are many items on the edge of being in range with some items of interest (I strike out those that are unlikely to be of great concern):

  • Faecalibacterium prausnitzii at 27% of the microbiome or 96%ile
  • Akkermansia — 0.009 % of the microbiome or 35%ile
  • Bifidobacterium 0.016 % of the microbiome or 16%ile
  • Bacteroides – 27% of microbiome, or 64%ile

Additionally, two indicate increased risk of Candida (new feature just added)

  • Phocaeicola dorei at 10% of the microbiome or 91%ile
  • Faecalibacterium prausnitzii at 27% of the microbiome or 96%ile

I would suggest a test for candida to be safe. The data suggests a risk. If confirmed, candida would contribute significantly to gut dysbiosis [The interplay between gut bacteria and the yeast Candida albicans[2021]). I did a “back-flip” check of the top prescription items, and all of them reduces Candida (studies cited below).

Addendum – Predicted Symptoms

This was just added to the site today as a further refactor based on New Special Studies on Symptoms data. These are from [My Profile Tab]

Criteria11/18/20215/20/20226/22/20239/4/2023
Forecast Major SymptomsNeurological: Cognitive/Sensory Overload
40 % match on 25 taxa

DePaul University Fatigue Questionnaire : Racing heart
38 % match on 13 taxa

DePaul University Fatigue Questionnaire : Difficulty falling asleep
37 % match on 27 taxa

DePaul University Fatigue Questionnaire : Difficulty finding the right word
35 % match on 20 taxa
Autonomic Manifestations: urinary frequency dysfunction
66 % match on 6 taxa
Immune Manifestations: Bloating
37 % match on 45 taxa

Neurological-Audio: hypersensitivity to noise
35 % match on 28 taxa
NoneNeurological-Sleep: Chaotic diurnal sleep rhythms (Erratic Sleep)
50 % match on 18 taxa

Neurological: Spatial instability and disorientation
37 % match on 16 taxa

This can be helpful for judging possible severity (and potential improvement of some symptoms), for example: Neurological: Cognitive/Sensory Overload. See [Special Studies] tab.

  • 2021 – 40% matches
  • 2022- 24% matches
  • 6/22/23 – 16% matches
  • 9/4/2023 – 4% matches

Going Forward

COVID has had quite an impact on this microbiome. I am going to just go with the “Just Give Me Suggestions” option with the addition of what matched his diagnosis:

  • Irritable Bowel Syndrome  (68 %ile) 7 of 68

To explain a bit more. First I click the button below

And then click I could click the consensus report to see what the top items are:

Which are shown below.

In this case, I want to add Irritable Bowel Syndrome suggestions (on the Changing Microbiome Tab)

Instead of the usual 4 packages of suggestions, we have 5

When we look at the consensus report we see the same items there, but the values have increased.

The intent is put a little bias on the numbers towards specific conditions of greatest concern.

PDF Suggestions

I tend to favor the PDF suggestions because it simplifies things for many readers. Also the PDF gives a good list of citations (never complete) used to make the citations to persuade MDs to see that the suggestions are based on studies — a lot of studies.

The PDF suggestions are below (using the consensus view is another option for those more technically orientated). I clip from the PDF to keep the blog simpler for the typical reader.

This is a little longer list than usual, so I went to the consensus report to get priority data. Top value was 618, so 309 is the 50% threshold.

These appear to be of low influence with the exception of l.bulgaricus:

Minor note: quercetin with resveratrol is an avoid, quercetin is a take. resveratrol by itself is a negative (-113). At times, you need to look at the technical details/consensus to clarify things; the data we are using is incomplete and sparse…. If clearly contradictory suggestions appear, then don’t do them (thing an abundance of caution).

Because he has an antibiotic friendly MD, the following are the TOP antibiotics with notes:

CFS Antibiotics are also above the threshold. Since the prior sample had a strong Long COVID or ME/CFS Profile, I would be inclined to include one of those below in the antibiotic rotation. The microbiome cannot make a diagnosis of most things, with most ME/CFS microbiomes there is a particular pattern which you had in your last sample but which has disappeared from your current sample which looks more like your first sample. I read this as recovering from ME/CFS….  in likely a fragile state since relapse is very common with ME/CFS.

My own experience is that it is better to overcure ME/CFS and when there are signs of recovery…. no backflips of joy or running marathons; keep doing slow walks that becomes a bit further each week for 6-12 months. Your microbiome is fragile and can quickly slip back.

I prefer to use the strategy of going for prescription items that are both suggested from the microbiome and been shown to help with one or more of the diagnosis conditions. This usually encounter low resistance from physicians — they are clueless for the microbiome, but very accepting of published studies. An antibiotic that is used as a prophylaxis usually encounter little resistance.

KEGG Suggestions

This is done by using information from the bacteria found with some fudge factors. I am in discussion with some Ph.D. candidates to build this concept directly from the FASTQ files and will hopefully have this as an added feature next year.

The KEGG probiotics is the usual pattern for ME/CFS and Long COVID with the top one being the usual, with the top reasonably available ones for other families shown below. I usually like to compare the values with those from consensus to minimum risk (i.e. two thumbs up, we do; mixed, we skipped)

KEGG Supplements

From the list, we will look only at those with a z-score (statistical significance) over 2. After each we put the consensus value (if it is listed)

Only two items are with high confidence.

How to Proceed Suggestions

The suggestions should be thought as influencers. The human population is often a good analogy or parable for the microbiome population. Each influencer shifts the population in the desired direction. Based on Cecile Jadin’s work and several studies, I am a firm believer in short duration (1-2 weeks) of each influencers. Just as with human influencers, people stop listening if the same person just keeps droning on and on. If a different person starts speaking, you get persuaded more. If a mob start to shout, yet a different human behavior will occur. In terms of the microbiome, “stop listening” means mutations that are resistant to the item will start to increase. Items line vitamins and minerals can be taken continuously; items that are likely to have bacteria resistance developed should be taken for a week and then another item replace it.

The items to rotate:

  • Antibiotics listed above
  • Probiotics: lactobacillus salivarius and lactobacillus bulgaricus
  • Herbs and spices: cinnamon, ginger, black cumin, thyme, rosemary, quercetin (suggests just before each antibiotic with a few days of overlap because it has potential synergistic activity with antibiotics [2020], [2016],[2018] )

Remember our goal is to destabilize a stable microbiome dysfunction.

Questions and Answers

While there has not been significant changes in many of the vectors between this sample and the prior sample from a few months earlier, there has been two significant objective changes:

  • Significant improvement of Anti inflammatory Bacteria Score (higher) and Histamine Producers (Lower).
  • The lost of the ME/CFS – Long COVID spike in the 0-9%ile

Q: Do you/should I use the colored list now instead of the consensus list?

  • Either are fine, the color list (from PDF) is what I tend to use in post because it is easier for new readers to understand (and automatically sent on new uploads). The consensus page is more complex but allows people to apply their own logic and priorities.

Q: “Quercetin (suggests just before each antibiotic with a few days of overlap because it has potential synergistic activity with antibiotics”

Q: I just did Mutaflor for 8 days and felt really tired all the time (but in the end I also got a flu/cold, so maybe that was the reason and not mutaflor). Nevertheless, if it was a herx reaction, I wonder if I should have taken it for longer until the reaction disappeared? (I stopped it 4 days ago.) Not sure if this question even makes sense.

  • My personal choice would be to keep taking it for at least a week (perhaps 2). Remember that the traditional pattern for a herx is feeling bad for X hours and then things get better. The duration of the feeling bad usually decrease from day to day. Catching a cold makes interpretation challenging.

Postscript – and Reminder

I am not a licensed medical professional and there are strict laws where I live about “appearing to practice medicine”.  I am safe when it is “academic models” and I keep to the language of science, especially statistics. I am not safe when the explanations have possible overtones of advising a patient instead of presenting data to be evaluated by a medical professional before implementing.

I can compute items to take, those computations do not provide information on rotations etc.

I cannot tell people what they should take or not take. I can inform people items that have better odds of improving their microbiome as a results on numeric calculations. I am a trained experienced statistician with appropriate degrees and professional memberships. All suggestions should be reviewed by your medical professional before starting.

The answers above describe my logic and thinking and is not intended to give advice to this person or any one. Always review with your knowledgeable medical professional.

Technical Note: Prevalence, Average and Not Reported

In reviewing many papers with the microbiome I noticed that often the researchers restrict their examinations to the taxa that is reported in all samples. I suspect this is due to a lack of sufficient statistical training and/or not understanding the natures of the microbiome.

Recently I came across these papers that uses an approach that I often have used, working off relative frequency of detection a.k.a. prevalence.

This post is going to use samples available at Microbiome Prescription Citizen Science site. We are going to restrict to one lab source and divide the data into two groups based on their self-declare symptoms and diagnosis.

  • Patients with Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) declared [Obs: 271]
  • Patients without Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) declared and other status declared (for example: “Asymptomatic” [Obs:569]

Naive First Pass

We are going to take the average count for each group ignoring no values reported. We are going to restrict it to taxa where we have at least 30 non-zero values [1,564 taxa]. We found some 77 taxa with a t-score over 2.81 (p < 0.005)

taxa nametaxa rankShiftT_score
Prevotella coprispecieslow in ME/CFS-5.27
Prevotellagenuslow in ME/CFS-4.52
Sporolactobacillaceaefamilylow in ME/CFS-4.2
Sporolactobacillus putidusspecieslow in ME/CFS-4.19
Sporolactobacillusgenuslow in ME/CFS-4.19
Prevotellaceaefamilylow in ME/CFS-4.1
Firmicutesphylumhigh in ME/CFS3.94
Blautiagenushigh in ME/CFS3.91
Cetobacterium cetispecieshigh in ME/CFS3.89
Cetobacteriumgenushigh in ME/CFS3.84

Deeming Not Reported to be Zero

In this case we have 78 taxa with a t-score over 2.81 with slight changes of t-scores.

taxa nametaxa rankShiftT_score
Prevotella coprispecieslow in ME/CFS-5.31
Sporolactobacillaceaefamilylow in ME/CFS-4.63
Sporolactobacillus putidusspecieslow in ME/CFS-4.62
Sporolactobacillusgenuslow in ME/CFS-4.62
Prevotellagenuslow in ME/CFS-4.5
Prevotella oulorumspecieslow in ME/CFS-4.35
Prevotellaceaefamilylow in ME/CFS-4.08
Bifidobacterium gallicumspecieslow in ME/CFS-3.97
Firmicutesphylumhigh in ME/CFS3.94
Blautiagenushigh in ME/CFS3.91

Prevalence

We followed the same process as above and limited things to a Chi-2 probability of < 0.005 (as used above) We ended up with 65 taxa.

tax_NameTax_RankPrevalence
in MECFS %
Prevalence
Control %
DifferenceChi2FoldChange
Deferribacteresphylum33.62013.5141.7
Erysipelothrix inopinataspecies2110.710.3142
Deferribacteralesorder33.62013.5141.7
Deferribacteraceaefamily33.62013.5141.7
Deferribacteresclass33.62013.5141.7
Mogibacterium vescumspecies27.715.811.9131.8
Haploplasma cavigenitaliumspecies8.52.85.7133
Haploplasmagenus8.52.85.7133
Gluconobactergenus15.16.98.3132.2
Prosthecobacter fluviatilisspecies7.72.55.3123.1

Comparing these two lists, we found only 6 taxa in common

  • Bifidobacterium angulatum
  • Propionigenium modestum
  • Pseudomonas viridiflava
  • Cetobacterium ceti
  • Cetobacterium
  • Propionigenium

The next result is that we have 78+65 – 6 = 137 statistically significant bacteria with p < 0.005.

Bottom Line

There are at least two different statistical ways of determining significance. IMHO, the prevalence approach is likely to be a superior tool for diagnosis purposes because it is possible to compute the probability of a match to the above patterns despite some bacteria not being reported.

The full list of bacteria is listed here.

Technical Notes: Statistics and Diversity Indices and PofP

A reader raised a valid question which actually triggers other related questions.

You seem to like the “percentage of percentiles” measurement, but I’m not convinced it’s being analyzed appropriately. As I understand it, you first convert to percentiles, getting numbers in [0, 100]. I think this is fine. Then you histogram these percentiles. Because each lab will perform the same measurements every time, I think this is also fine. However, the result is compositional data in the sense of Aitchison, and it should be analyzed in a manner consistent with that. For compositional data, a chi^2 test is inappropriate because it relies on the number of species (or genera) measured.

My suggestion is to apply a centered logratio transform to each person’s percentages and fit a normal distribution to the transformed data. To determine whether someone’s microbiome deviates significantly, calculate a multivariate normal tail probability. Beware that the covariance matrix will be rank deficient (you’re in a ten-dimensional space, but there are only nine parameters because percentages sum to 100). You may want a robust fit because it’s reasonable to expect that the microbiome of someone ill might be an outlier.

For more information about compositional data, see Aitchison, J., “The Statistical Analysis of Compositional Data,” Journal of the Royal Statistical Society. Series B (Methodological) Vol. 44, No. 2 (1982), pp. 139-177; Aitchison, J., “The Statistical Analysis of Compositional Data,” Chapman & Hall, London, 1986; and Aitchison, J. “A Concise Guide to Compositional Data Analysis,” unpublished manuscript, 2005, available online (just Google). For other approaches to compositional data analysis, see Greenacre, Michael; Grunsky, Eric; Bacon-Shone, John; Erb, Ionas; Quinn, Thomas, “Aitchison’s Compositional Data Analysis 40 Years On: A Reappraisal,” arXiv:2201.05197, 13 Jan 2022, to appear in Statistical Science.

What is the statistical basis for other Diversity Indices?

How to calculate these numbers is well determined — they seem to be brilliant ideas tossed out there that seems to fit the data for some study. For some background, see this page. The problem is a lack of rigor, especially statistical rigor.

Diversity indices, particularly the Shannon-Wiener index, have extensively been used in analyzing patterns of diversity at different geographic and ecological scales. These indices have serious conceptual and statistical problems which make comparisons of species richness or species abundances across communities nearly impossible. 

Conceptual and statistical problems associated with the use of diversity indices in ecology [2009]

The problem is an absence of a native statistical model. For example, it does not fit the usual ones.

The key question is simple, what is the distribution underlying diversity Indices? We read ” In the literature of biodiversity, according to Ricotta (2005), there are a “jungle” of biological measures of diversity.”[2017]. Zheng’s A new diversity estimator[2017] in Journal of Statistical Distributions and Applications where he states “There are many other open problems built on this connection between birthday problem and diversity measures. ” The problem is this, the birthday problems deals with 366 discrete well defined boxes that are well defined. Dealing with the microbiome, we lack these boxes. Consider a measure of a microbiome sample in 2000, there are a large number of different bacteria species in Lactobacillus. Today, we have these species no longer placed in 1 genus, but 25 genus [2020] including:

  • Acetilactobacillus,
  • Agrilactobacillus,
  • Amylolactobacillus,
  • Apilactobacillus,
  • Bombilactobacillus,
  • Companilactobacillus,
  • Dellaglioa,
  • Fructilactobacillus,
  • Furfurilactobacillus,
  • Holzapfelia,
  • Lacticaseibacillus,
  • Lactiplantibacillus,
  • Lapidilactobacillus,
  • Latilactobacillus,
  • Lentilactobacillus,
  • Levilactobacillus,
  • Ligilactobacillus,
  • Limosilactobacillus,
  • Liquorilactobacillus,
  • Loigolactobacilus,
  • Paucilactobacillus,
  • Schleiferilactobacillus,  
  • Secundilactobacillus.

With the same strains/species, our diversity indices will be very different because our boxes are arbitrary and “soft” unlike the days of the year or the roll of a dice.

Back to percentage of percentiles

While I show genus and species in the table for ease of understanding of the typical reader, I originally did it solely with the lowest identifiable levels (the “atoms” or the microbiome) – species. At the species level, it is not compositional. There is no composition! Looking at the data that was actually received, I noticed many genus had no species listed. In some cases, the genus had species, but none of the known ones were detected. In other cases, the test did not report any species in over 3000 test results.

On this basis I decided to use try using both species and genus. I soon discovered that they almost always exhibit a similar pattern and chi^2. At this point, I opted for benefiting my readers and not as much rigor as some would like. We could do the lowest taxonomical level reporting across the hierarchy as one solution.

This approach ends up with us side-stepping the classification issues cited above. We are dealing with distinctive, non-overlapping events (a bacteria being identified) and then convert them to percentile giving use a continuous uniform distribution for each of these independent events. IMHO, at this point we have a good model to chi^2 test. We are not dealing with measuring a population, just a sample.

In answer to “a chi^2 test is inappropriate because it relies on the number of species (or genera) measured.” is missing the point. If I get two bags of coins from the bank and then flip them to determine if they are biased — whether the bags contains 1000 or 100,000 coins is significant only on the ability to determine the margin of error. The number of species/genus is only significant in that sense. If there is a strong bias with a small number, then having more will not change the bias.

Technical Note: Lab Quality Versus Bacteria Reported

For samples coming from uBiome, Ombre/Thryve and Biomesight there are two important numbers reported. [Count] and [Count_Norm].

  • [Count] is the number coming from the lab equipment, the bacteria detected.
  • [Count_Norm] is the above number scaled to be out of one million (1,000,000)

The question arises, if you have low lab quality and the number of bad bacteria also dropped: Is this an actual improvement or a false improvement due to low lab quality?

To partially answer that question, I pulled biomesight samples (biggest collection) and plotted the data. Shown below:

  • Vertical axis is number of bacteria reported
  • Horizontal axis is lab quality measure

As is shown, there is a relationship.

Using this data and sample values of 4.3 and 8.4, we compute 546 and 643 for expected bacteria (just use the formula on the chart). This means that with 4.3 we expect only 85% (546/643 * 100) of the count seen with 8.4.

So we see the changes below are likely true improvements: (Left is 4.3, right is 8.4)

For bacteria reported by lab, the numbers suggests the left column has less odd bacteria and the gut microbiome may becoming more uniform.

For anyone interested in doing their own charting and analysis, the raw data is https://citizenscience.microbiomeprescription.com/

The foolishness of Numeric Health Measures for the Microbiome

Today I was asked if the numbers shown on various sights like the one below, ” if this result is correlated with the severity of the patient?”. The answer is No. The usual reason that such numbers appear on sites is to satisfy customers asking for such numbers. Simple, easy to understand number.

Businesses want to make customers happy — so they literally cook-up a number to show on their reports. There is no research supporting any of the magical numbers that I have seen. Some one put together some numeric formula to generate the numbers.

For Biomesight, the logic is shown on the page with the right display

If you sum up the values for each dial, you get the total. When you view the percentage, it is not so obvious: 100%, 69%, 88%, 85% looks like a complex formula is being used.

Analogy: Give me a Health Measure for a Person

The microbiome is very much like a person. How would you create a single value for a person?

  • Probiotics –> Income
  • Commensals –> Savings
  • Pathobionts –> Debts
  • Diversity –> Health?
    • Married and marriage status
    • Chronic Conditions
    • Health Status
    • Height
    • Weight
    • Gender
    • etc

Would a person with no debt, good savings, a low income and married with 5 kids be a higher or lower measure than someone with moderate debt, high income, small savings and no relationships?

Are there any studies?

What about the ratios used in the literature?

  • Bacteroides/Bifidobacterium Ratio
  • Bacteroides/Clostridium Ratio
  • Blautia/Bacteroides Ratio
  • Firmicutes/Bacteroidetes Ratio 
  • Prevotella/Bacteroides Ratio

Well, the sample that I am looking at have 43%ile, 98%ile, 7.2%ile, 58%ile and 19%ile. The numbers are all over the place!

We also have two studies:

Their latest study states: ” accuracy of 80% in distinguishing healthy (no disease) from non-healthy (diseased)” So 1 in five times, they will misclassify a healthy person as having a disease.

An example of the algorithm being used is below.


This is just trying to determine if the person is healthy or not — nothing about severity.

My approach is simple, I simply flagged the bacteria deemed to be unhealthy. Just list them.

How to use it:

Same Sample – 2 Labs: 16s vs Shotgun

Back Story

Latest microbiome results are in. Seems like my microbiome is stubborn and stuck these last few tests. Wondering if I should just use this test suggestions exclusively or combine with my prior Thorne test 

We have two sample – one via BiomeSight and one via Thorne. This post is going to do two things:

  • Look at Suggestions – by combining both sets of suggestions using the Uber Consensus
  • Look at the differences between the reports.

We also review “which is better”. My focus is clinical application to individuals — not research papers; answer at bottom.

Uber Consensus

The process has become very simple — “Just give me Suggestions!” on both samples and then going to uber consensus as illustrated by the video below.

The result was excellent agreement on suggestion between each set of results. The CSV files are attached below.

Differences between Reports

I compared two things between the reports:

  • Percentage of the bacteria in the microbiome
  • Percentile of the bacteria in the microbiome

At the Phylum Level

Items less than 100 should be ignored (accuracy of measurement limits). There are a few dramatic differences.

Bacteria NameThorne CountBiomeSight Count
Firmicutes396799529540
Actinobacteria606102100
Bacteroidetes461289448230
Proteobacteria609518150
Chlorobi36429
Acidobacteria35100
Cyanobacteria8320
Spirochaetes8530
Verrucomicrobia5910
Chloroflexi7750
Tenericutes5430
Deinococcus-Thermus4830
Fibrobacteres410
Synergistetes1720
By Count

Looking at Percentiles next

Bacteria NameThorne %ileBiomeSight %ile
Chlorobi2584
Actinobacteria8533
Acidobacteria3481
Spirochaetes8136
Cyanobacteria311
Deinococcus-Thermus5529
Firmicutes1437
Chloroflexi6750
Verrucomicrobia141
Tenericutes132
Proteobacteria1018
Synergistetes64
Bacteroidetes5556
Fibrobacteres10
By Percentile ranking

We have Bacteroidetes in agreement with both — but for the rest…

At the genus level

Bacteria NameThorne CountBiomeSight Count
Bacteroides180054397640
Blautia16470107220
Roseburia1679373640
Faecalibacterium109196152890
Corynebacterium43413820
Ruminococcus917744170
Phocaeicola223209199669
Parabacteroides1185531940
Phascolarctobacterium610123980
Dorea3613000
Sutterella1611339
Oscillospira08250
Coprococcus612012589
Eggerthella6491760
Pseudobutyrivibrio1495790
Lachnospira115936230
Prevotella9544260
Anaerostipes93036310
Clostridium20394960
Pedobacter462410
Odoribacter40772060
Bifidobacterium27831019
Escherichia751610
Porphyromonas1372150
Mediterraneibacter1483113629
Bilophila61110
Veillonella751160
Desulfovibrio19001250
Streptococcus1477840
Acetivibrio33470
Chlorobaculum6429
Finegoldia1339920
Gemella17400
Enterococcus585220
Paenibacillus37620
Mogibacterium39370
Acetobacterium15340
Serratia47350
Eubacterium517240
Megasphaera35290
Selenomonas52290
Bacillus24810
Caldicellulosiruptor11240
Campylobacter23510
Slackia16240
Sphingobacterium48270
Caloramator10190
Staphylococcus18110
Hathewaya8170
Peptoniphilus656800
Peptostreptococcus6150
Microbacterium12510
Adlercreutzia525620
Rhodothermus690
Erysipelothrix1290
Acidaminococcus1290
Hymenobacter8010
Negativicoccus11550
Collinsella7410
Rhodococcus6710
Dialister2580
Anaerococcus336390
Pseudoclostridium860
Moorella960
Vibrio6010
Caldilinea150
Brochothrix250
Mycobacterium6720
Neisseria5710
Pectinatus750
Thermoclostridium1650
Alkaliphilus940
Shewanella3160
Lactobacillus5730
Leptospira430
Deinococcus3510
Tetragenococcus530
Ethanoligenens3410
Weissella1030
Gulosibacter120
Pseudoclavibacter220
Kocuria2810
Meiothermus220
Stenotrophomonas2810
Symbiobacterium320
Devosia420
Dysgonomonas3420
Azoarcus2110
Leuconostoc920
Glaciecola110
Turicibacter2130
Pelotomaculum110
Parascardovia210
Lentibacillus210
Actinopolyspora210
Kitasatospora210
MLOs310
Ochrobactrum310
Rickettsia310
Luteibacter310
Fibrobacter410
Pediococcus1420
Halanaerobium610
Dyadobacter1410
Mycoplasma1720
Thauera910
Lysobacter1110
By Counts

Looking at the percentile rankings — the absolute numbers may vary greatly, but what about relative percentiles?

Bacteria NameThorne %ileBiomesight %ile
Ochrobactrum22
Actinopolyspora11
Halanaerobium11
MLOs11
Glaciecola11
Lentibacillus11
Pelotomaculum11
Parascardovia11
Luteibacter11
Phocaeicola8989
Rickettsia10
Pediococcus109
Fibrobacter20
Mycoplasma53
Alkaliphilus13
Finegoldia8588
Kitasatospora30
Thauera51
Streptococcus5550
Turicibacter1217
Peptoniphilus6458
Hathewaya18
Clostridium1811
Desulfovibrio6169
Eubacterium3846
Symbiobacterium19
Enterococcus8879
Sphingobacterium1323
Pseudoclavibacter111
Anaerococcus7283
Eggerthella9886
Gulosibacter012
Lactobacillus2311
Bifidobacterium5543
Leuconostoc214
Shewanella3547
Prevotella5063
Corynebacterium9986
Collinsella130
Oscillospira016
Faecalibacterium4965
Meiothermus117
Caloramator119
Coprococcus3957
Lysobacter180
Odoribacter8163
Adlercreutzia6381
Pedobacter1331
Dyadobacter201
Dysgonomonas244
Mediterraneibacter6990
Devosia122
Acetivibrio527
Thermoclostridium932
Ethanoligenens251
Dialister1135
Veillonella1641
Pectinatus127
Porphyromonas8862
Moorella128
Negativicoccus6639
Lachnospira5121
Rhodothermus132
Tetragenococcus132
Acetobacterium334
Anaerostipes6596
Bilophila133
Ruminococcus1447
Weissella235
Parabacteroides4275
Acidaminococcus439
Pseudoclostridium137
Leptospira142
Serratia3475
Slackia445
Phascolarctobacterium5697
Erysipelothrix446
Sutterella146
Bacteroides3987
Roseburia4391
Escherichia2877
Selenomonas2173
Deinococcus541
Megasphaera1872
Brochothrix156
Kocuria582
Mogibacterium1774
Stenotrophomonas633
Azoarcus610
Caldilinea061
Caldicellulosiruptor264
Mycobacterium8724
Hymenobacter681
Blautia573
Paenibacillus8719
Neisseria690
Pseudobutyrivibrio2595
Campylobacter751
Gemella482
Peptostreptococcus181
Chlorobaculum184
Staphylococcus850
Vibrio912
Bacillus921
Rhodococcus910
Dorea193
Microbacterium941
By Percentile

Bottom Line

I have never had much belief in the absolute accuracy of the bacteria named or the count. Why? Simple, I understand the statistical process being used and its weakness. I will leave arguments over “which is better” and “which is accurate” to others.

See The taxonomy nightmare before Christmas… for more information.

My focus and concern is to improve the microbiome. With sparse data and the great complexity involved, I am actually very pleased that the suggestions are in agreement. The suggestions are computed using fuzzy logic expert systems. The noise in the data and the statistical processes involved seem to be smoothed out by this Artificial Intelligence engine approach.

Score: Labs: -2, Microbiome Prescription 2

Which is better?. My focus is clinical application to individuals, both give similar suggestions using the Fuzzy Logic Expert System. There is no difference in that sense.

ME/CFS recovery short circuited by COVID

I have been doing periodic review of this person’s sample. He just got his latest results and it was a shocker (of the wrong type!). This is worth a review.

Comparing Samples Overtime

At the typical analysis level there has been no change in these broad criteria since the last sample:

  • Outside Range from JasonH
  • Outside Range from Medivere
  • Outside Range from Metagenomics
  • Outside Range from MyBioma
  • Outside Range from Nirvana/CosmosId
  • Outside Range from XenoGene

Why compare over 90%ile to under 10%ile? The reason is probability – we are converting the data to a uniform distribution for all of the bacteria. This allows for reliable statistical significance to be determined for all of the bacteria. If things are “normal” then the ratio should be 1.0 The further from 1, the more abnormal. This is independent of any assumptions on bacteria distributions.

Looking at over 90% and under 10%, our expected count are 64 for both (10% of 639)

  • We have 20 over 90%ile, so we have under representation of dominant
  • We have 273 under 10%ile, the typical over representation of low levels of many bacteria seenwith most ME/CFS people
  • The ratios is higher at 13.7 compare to prior ratios (8.2, 11,3,3,6.5)

Where we see differences

  • Outside Kaltoft-Moldrup count returned to the size of the very first sample.
  • Compounds over 90%ile to under 10%ile (which should be 1.0 theoretically) had been close to 1.0 on the prior 3 samples, jumped up to 8.8:1. This was not as bad as the first sample with a 12.9:1 ration.
  • Enzymes over 90%ile to under 10%ile, continue to be bias towards low with a 2.7:1 ratio (prior 3.5, 1.7, 2.9, 3.2)
  • Conditions: jumped from none over 90%ile to 13!

My general impression is that ground has been lost. This is the first time that subsequent results appear to be worse. WHAT HAPPENED!???!!???!!!

Criteria2/22/20228/11/20223/25/202212/3/20218/31/2021
Lab Read Quality9.75.56.23.67.8
Bacteria Reported By Lab639461593445551
Bacteria Over 99%ile433515
Bacteria Over 95%ile1113112423
Bacteria Over 90%ile2023214035
Bacteria Under 10%ile273189237123227
Bacteria Under 5%ile21910714366192
Bacteria Under 1%ile17523449167
Lab: BiomeSight
Rarely Seen 1%671423
Rarely Seen 5%22143379
Pathogens3732463138
Outside Range from JasonH77446
Outside Range from Medivere1515151519
Outside Range from Metagenomics88667
Outside Range from MyBioma77778
Outside Range from Nirvana/CosmosId2323181821
Outside Range from XenoGene3232363639
Outside Lab Range (+/- 1.96SD)786914
Outside Box-Plot-Whiskers3833385841
Outside Kaltoft-Moldrup210111123100211
Condition Est. Over 99%ile50007
Condition Est. Over 95%ile900015
Condition Est. Over 90%ile1300029
Enzymes Over 99%ile3510301972
Enzymes Over 95%ile1006821982162
Enzymes Over 90%ile191183296126192
Enzymes Under 10%ile520645514369616
Enzymes Under 5%ile375423264186450
Enzymes Under 1%ile219864937272
Compounds Over 99%ile2347622844
Compounds Over 95%ile7225423112786
Compounds Over 90%ile12633829830798
Compounds Under 10%ile11043082972271265
Compounds Under 5%ile10681732241111241
Compounds Under 1%ile10456567471206

What Happened?

This person sent the following notes

  • I’ve got COVID in October, I feel as I have fully recovered.
  • I have a little bit more energy than before.
  • My body feels extremely stressed, I have started to get a pressure over the neck / thyroid when I get totally stressed out.
  • My sleep is much better. I have been following Andrew Hubermans protocol for sleep, which had a great impact on me. 
  • Would be great if I could get some recommendations for food, supplements, antibiotics etc.

Ah, the person feels like he has recovered but his microbiome is still recovering. We have a clean explanation for the regression! Our goal is now to try to stop Long COVID from setting in.

The Conditions matches include: hypercholesterolemia (High Cholesterol), Hyperlipidemia (High Blood Fats), Hypertension (High Blood Pressure, Nonalcoholic Fatty Liver Disease (nafld) Nonalcoholic and Atherosclerosis. None of those are concerning –they were not matches last time and thus should be viewed as transient red herrings. Looking at PUBMED Long COVID explicitly, we see good news: Long COVID   (29 %ile) 37 of 212. However when we go over to Special Studies, the very top one is

  •  41 % match COVID19 (Long Hauler), the next match was 21% –this really sticks out!
    • All Prior samples also had Long Hauler at the top too.

Going Forward

To build our consensus, we will do the usual and toss in our top Special Studies one.

Remember that we have a massive over representation of low %ile bacteria — hence some selection methods produce must larger numbers.

The suggestions downloads are below.

The top items are two similar probiotic mixtures from studies:

Looking at probiotics and the components above — bifidobacterium breve (probiotic) and propionibacterium freudenreichii appear to be the best of the components. The last one is used to make Emmental and Jarlsberg cheese or is available as a probiotic: Nutricology/Securil.

Given the general hostility between lactobacillus and E.Coli, plus the risk or lactic acid issue with lactobacillus, I would suggest avoiding lactobacillus casei initially, perhaps try it in a later cycle. A similar contradiction happened with different forms of cranberry as shown below. We want absolutely clean positive choices.

If different results from similar versions of something — then AVOID.

I reran the suggestions with everything — in case any prescription items may be of special interest.

The top antibiotics included:

As always, I prefer the Cecile Jadin approach of taking a single course, take a break and then take a different antibiotics.

Looking at Foods

The top food suggested was Protein powder, whey based, protein >70%, unfortified (ignore the unfortified – go for the standard used by ME/CFS Physicians: Undenatured Whey). This is a recommended with good results by many ME/CFS physicians like Teitelbaum and Dr. Paul Cheney. The second choice was duck liver which is not easy to obtain in many places (with other liver below). #3 items was Rye, whole grain flour, and near the top of the list: Confectionery, peanut, chocolate-coated which I chuckled over: See my 2012 post, Honestly! Chocolate!!! (Yes 11 years ago!) and 2013 post, Peanut Butter – a complex food?

I personally have always love Liver pâté! There was something that always felt so good after eating. Liverwurst is essentially the same food.

So for a Scandinavian, this is almost going to a church social!! Liver pâté on Dark 100% Rye bread with Jarlsberg cheese also on Dark 100% Rye bread!

Postscript – and Reminder

I am not a licensed medical professional and there are strict laws where I live about “appearing to practice medicine”.  I am safe when it is “academic models” and I keep to the language of science, especially statistics. I am not safe when the explanations have possible overtones of advising a patient instead of presenting data to be evaluated by a medical professional before implementing.

I can compute items to take, those computations do not provide information on rotations etc.

I cannot tell people what they should take or not take. I can inform people items that have better odds of improving their microbiome as a results on numeric calculations. I am a trained experienced statistician with appropriate degrees and professional memberships. All suggestions should be reviewed by your medical professional before starting.

The answers above describe my logic and thinking and is not intended to give advice to this person or any one. Always review with your knowledgeable medical professional.

Beyond simple naïve microbiome analysis?

Foreword

This post started out with a title of “Post-Acute COVID-19 Syndrome vs Myalgic Encephalomyelitis – Similarities and Differences“. It scope was pretty obvious — compare microbiome shifts from these two sibling conditions. Expectations was a bland informational review.

The result was calling into question the typical assumption that we could isolate symptoms and conditions to specific bacteria. I follow the statistics and discovered that you can get a magnitude better statistical significance by moving beyond bacteria. In coming weeks, I hope to code up suggestions AI based on this enlightenment.

My starting point

In my last post, Long COVID – an update, I did a comparison between the citizen science data and the literature published on the US Library of Medicine. In terms of symptoms, these two appear the same — but at the microbiome (and enzyme levels) how similar are they?

  • Post-Acute COVID-19 Syndrome (PCAS), also known as Long COVID
  • Myalgic Encephalomyelitis (ME), also known as Chronic Fatigue Syndrome (CFS)

One frustrating aspect of many studies on the US Library of Medicine for many conditions is simple: results are not replicated in subsequent studies for the same condition. Historically I have viewed this as a result of different equipment and different reference libraries. In many cases the bacteria deemed significant are often different and when they did report the same bacteria, they report opposite shifts!

This post explores some of these issues, and came to an interesting conclusion.

Study Caveats

The studies on the US Library of Medicine compare people with the condition to healthy controls. With the citizen science data that is almost impossible to do. If a person has gotten a microbiome test, they likely have some condition(s) and thus are not healthy controls!!

This is not all bad. It means that when we find things that are statistically significant they are differentiators against other people with microbiome issues. That is, how are people with ME different than people with FM and IBS. Conceptually, we are more likely to identify the key features for these conditions and not key features for auto-immune conditions or a gut disturbance in general. It is a nuisance difference, but may be a very important nuisance.

Comparison that we will review are from:

  • US Pubmed — bacteria reported by both with direction
  • KEGG Enzymes shifts from Citizen Science (using only Biomesight data)
  • Bacteria shifts from Citizen Science (using only Biomesight data)
Citizen Science Samples
Studies on US National Library of Medicine

For citizen science we may have many uploaded samples annotated both with PCAS and ME. To resolve this conflict, ME will contain only samples with ME and without PCAS. Both ME and PCAS have many, many comorbid symptoms which may also come into play. Many of the pure ME samples are before COVID swept the world, hence relatively clean. PCAS are more recent samples.

For PACS citizen science data, we have only significance difference identified from Biomesight data, hence we will compare those only.

ScopeMEPACSSame
US National Library of Medicine6823325
Enzymes – Citizen Science with p < 0.00122819931
Bacteria – Citizen Science109360
Entities reported as significant or found significant

I must admit that finding no bacteria in common with the same lab and the same reference library was a little bit of a surprise. One explanation is that microbiome dysfunctions evolve over time. People with PACS have had it less then 3 years, likely an average of just 1 year. People with ME has had it often for 30+ years. Comparing the two may be similar to comparing a one bottle of grape juice to a bottle of vintage wine.

Details for Common Bacteria from US National Library of Medicine

In the table below: H indicates High, L indicates Low.

Note that Bacteroides are reported high and low in different studies, suggesting there are subsets of each condition

tax_ranktax_NameDirection
classBacteroidiaH
familyBacteroidaceaeH
familyClostridiaceaeH
familyLachnospiraceaeL
genusAnaerostipesL
genusBacteroidesH
genusBacteroidesL
genusBifidobacteriumL
genusCoprobacillusH
genusCoprococcusL
genusDoreaL
genusEggerthellaH
genusEnterococcusH
genusFaecalibacteriumL
genusLactobacillusL
genusStreptococcusH
genusTuricibacterH
orderEubacterialesL
phylumBacteroidetesH
phylumFirmicutesL
phylumFusobacteriaH
speciesAnaerobutyricum halliiL
speciesEnterocloster bolteaeH
speciesFaecalibacterium prausnitziiL
speciesRuminococcus gnavusH
From https://microbiomeprescription.com/Library/PubMed

Details for Shared Enzymes with p < 0.001

In recent posts for conditions comorbid with ME, PACS, I found that enzyme analysis had greater statistical significance than bacteria. All of these posts reported higher enzyme levels were significant with these conditions.

The result for items shared that had p < 0.001 was almost overwhelming!

ECKeyEnzymeName
1.1.1.2921,5-anhydro-D-mannitol:NADP+ oxidoreductase
1.12.98.4H2:polysulfide oxidoreductase
1.7.2.2ammonia:ferricytochrome-c oxidoreductase
1.8.7.3CoB,CoM:ferredoxin oxidoreductase
1.8.98.4CoB,CoM,ferredoxin:coenzyme F420 oxidoreductase
1.8.98.5CoB,CoM,ferredoxin:H2 oxidoreductase
1.8.98.6coenzyme B,coenzyme M,ferredoxin:formate oxidoreductase
2.3.1.201acetyl-CoA:UDP-2-acetamido-3-amino-2,3-dideoxy-alpha-D-glucuronate N-acetyltransferase
2.7.1.2271-phosphatidyl-1D-myo-inositol:a very-long-chain (2’R)-2′-hydroxy-phytoceramide phosphoinositoltransferase
2.7.8.12CDP-glycerol:4-O-[(2R)-glycerophospho]-N-acetyl-beta-D-mannosaminyl-(1->4)-N-acetyl-alpha-D-glucosaminyl-diphospho-ditrans,octacis-undecaprenol glycerophosphotransferase
2.7.8.36UDP-N,N’-diacetylbacillosamine:tritrans,heptacis-undecaprenyl-phosphate N,N’-diacetylbacillosamine transferase
3.1.1.114methyl acetate acetohydrolase
3.1.3.27phosphatidylglycerophosphate phosphohydrolase
3.1.6.6choline-sulfate sulfohydrolase
3.1.6.8cerebroside-3-sulfate 3-sulfohydrolase
3.10.1.1N-sulfo-D-glucosamine sulfohydrolase
3.2.1.116-alpha-D-glucan 6-glucanohydrolase
3.2.1.152mannosylglycoprotein endo-beta-mannosidase
3.2.1.197beta-1,2-D-mannoside mannohydrolase
3.2.1.24alpha-D-mannoside mannohydrolase
3.4.21.26prolyl oligopeptidase
4.1.99.1L-tryptophan indole-lyase (deaminating; pyruvate-forming)
4.2.2.20chondroitin-sulfate-ABC endolyase
4.2.2.21chondroitin-sulfate-ABC exolyase
4.2.2.3alginate beta-D-mannuronate—uronate lyase
4.2.2.8heparin-sulfate lyase
4.3.1.7ethanolamine ammonia-lyase (acetaldehyde-forming)
5.1.1.20L-alanyl-D-glutamate epimerase
5.1.3.11cellobiose 2-epimerase
5.3.1.22hydroxypyruvate aldose-ketose-isomerase
6.1.1.13D-alanine:poly(phosphoribitol) ligase (AMP-forming)

One of them caught my eye, heparin-sulfate lyase, because micro-clots and “thick blood” are associated with these conditions with good results reported from the use of heparin for some patients.

 HSGAGs are widely distributed on the cell surface and extracellular cell matrix of virtually every mammalian cell type and play critical role in regulating numerous functions of blood vessel wall, blood coagulation, inflammation response and cell differentiation.

Microbial heparin/heparan sulphate lyases: potential and applications [2012]

Bacteria – Citizen Science

This blew me away — we have over 150 people with PCAS and over 250 with ME giving us superior sample sizes. We have 145 bacteria deemed significant for one or the other. We had NONE that was in common. This gut punch gives two main options: denial and look for an excuse to disregard. or roll with the punch and roll with enzymes.

Bottom Line

The enzyme aspect is the strongest association between PACS and ME. The count is higher, but more important, we are dealing with p < 0.001 data making false positives unlikely. This leads to a model that there is no ‘magical list of bacteria involved’ but a diverse array of bacteria that may be different for each person — but united in the over production of enzymes. This shifts the microscope of research into a different light spectrum. This is very interesting and may require some new brain cells to be used.

Using this information to improve..

If Enzymes estimate trumps bacteria levels (in a statistical sense), then we need to look at the enzyme levels and deduce for each one of concern, which collection of bacteria needs reduction — even when they are in the normal range. It is the aggregation of bacteria where the issue arises, not individual bacteria or specific subsets of bacteria.

A few examples may illustrate things a little

Example: (S)-3-hydroxy-3-methylglutaryl-CoA acetoacetate-lyase (acetyl-CoA-forming) a.k.a. EC 4.1.3.4, which was the most significant for PACS in the post: Long COVID – an update. There are some 2000+ taxon associated with it. We look at the averages for these below.

RankTax_NameWith PACSWithout PACSTScoreDF
speciesFaecalibacterium prausnitzii1381511096043.796775667
speciesPseudomonas viridiflava53252.62810832
speciesComamonas kerstersii125402.60038854
speciesPseudomonas aeruginosa62311.82464443
speciesEmticicia oligotrophica23039671.727619455
speciesDenitratisoma oestradiolicum42241.65065722
speciesGranulicella tundricola29211.6122548
speciesBacillus subtilis40191.37443117
speciesNiabella soli31240.96394316
speciesRalstonia insidiosa53380.91487436
speciesOligella ureolytica51320.8976819
speciesGlaciecola nitratireducens27240.6789966
speciesBacillus halotolerans32280.37121958
speciesAcidaminococcus intestini7496240.357674146
speciesAcinetobacter guillouiae67630.09895618
Key Contributors to EC 4.1.3.4,

For another one, we see the pattern stronger. Below we see the difference of Faecalibacterium prausnitzii is around 30,000 units. Looking at the other contributors, we see an additional 40,000 units. These extra units doubles the shift (and thus significance) of the enzyme above that of a single bacteria. Some of these are deemed healthy usually, for example: Akkermansia muciniphila which was at the 78%ile for Without PACS and 83%ile with PCAS. Neither would be deem to be outliers.

RankTaxon NameWith
PACS
Without
PACS
t-scoreDF
speciesFaecalibacterium prausnitzii1381511096043.796775667
speciesSutterella wadsworthensis962667722.380718452
speciesAliarcobacter skirrowii3756212.22360217
speciesAkkermansia muciniphila19096122901.896922547
speciesDesulfovibrio desulfuricans14234691.76977232
speciesEmticicia oligotrophica23039671.727619455
speciesEnterococcus casseliflavus1965811.63566638
speciesPorphyromonas asaccharolytica13502541.59988186
speciesBacteroides fragilis808055951.523991489
speciesBifidobacterium dentium14544611.433823239
speciesPhocaeicola dorei35482290751.396731649
speciesCorynebacterium aurimucosum11054071.27523496
speciesBacteroides eggerthii14379103451.108857263
speciesCorynebacterium jeikeium18977230.85828270
speciesPhocaeicola coprophilus649636420.856783152
speciesDesulfovibrio piger203215340.848976141
speciesMegamonas funiformis167711300.62028390
speciesHathewaya histolytica289027290.467066660
speciesHaemophilus parainfluenzae134312500.282656500
speciesMesoplasma entomophilum118210690.230055294
speciesPhocaeicola vulgatus51403512130.034398665
Key Contributors for EC6.1.1.6

Going Forward

The logical approach is simple to describe. For a person with the symptom, determine the enzymes which are abnormal. Determine the bacteria that are too high (even if only a little). Then use the suggestions AI Engine to determine the substances that will affect the greatest number of these bacteria to shift in the desired direction without encouraging other bacteria that could contribute to these enzymes to increase.

Now, the mathematics and complexities of this computation is a different matter but well within the power of today’s computer.

German CFS Patient got COVID….

This request came from the person discussed in Follow up to: A German CFS Patient Experience and Analysis.

I don’t know, if you remember me, we did two reports together, and your suggestions really helped to get my microbiome back on track (which shows in the samples).

And then I got COVID in November 2021-December 2021. But I felt better with it, but unfortunately I couldn’t give up my sample while having COVID.

Anyway I wanted to ask you whether you may be interested in my case, because I had a huge, irreversible it seems, crash from 20-30 on Bell CFIDS disability scale to now under 10 and my microbiome crashed along with me . (The crash also resulted in a high number of Lorazepam intake from which I’m slowly withdrawing now. But I didn’t get a clear idea of the effect of Lorazepam on the microbiome, other then they make the slowed gut motility worse of course.)

I have a very severe and have a progressive form of ME/CFS in the way that whenever I really crash I always go down to a lower baseline and do not recover. And with most crashes I loose about 50 % of my functionality, so it just took me one year to go from very mild to very severe.

From my lowest point onwards I’ve always had to take about 4 Lorazepam to guarantee a minimum of a bearable quality of life in bed. I succeeded for 4 years not to have a major crash and did did not build up a tolerance towards Lorazepam in that time.

A lot of things seem to have reversed, what I should take before are often things that I now should avoid. What Biomesight says seems to contradict slightly from what I can gather from your site. (Yes, I know you explained why there can be contradictory results).

And I have difficulties getting the suggestions for the handpicked criteria to show. Of course I would be super glad, if you could help, but I understand if you have more interesting projects to work on. (I would of course donate for your effort, as this is the only or easiest way to say thank you),

I believe one of the differences between Microbiome Prescription(MP) and Biomesight is simply the number of studies used to make suggestions. At present, we have over 11,000 studies coded into MP, I do not know the number that Biomesight uses, but I expect less than 1%. Also, MP suggestions was written by a person that has worked professionally in Artificial Intelligence. I suspect Biomesight lacks that skill set for development. Regardless, put items not in disagreement as first priority.

Analysis — The Numbers

There is no magic number that answers questions about the microbiome. Usually, I look for abnormalities. Since the earlier post, she had 6 more microbiome samples done periodically and shown below. She is wise to regularly monitor and ideally take moderate steps (diet and supplements) to counter any concerning trends.

CriteriaSep21Mar22May22Sep22Jan23May23
Shannon Diversity Index78.294.367.153.998.984.70
Simpson Diversity Index30.740.744.417.542.948.90
Chao1 Index53.666.881.536.765.161.90
Lab Read Quality4.87.37.75.256.5
Bacteria Reported By Lab612653717536636642
Bacteria Over 99%ile2159684
Bacteria Over 95%ile42050333212
Bacteria Over 90%ile294469506838
Bacteria Under 10%ile44181181405343
Bacteria Under 5%ile1216416592014
Bacteria Under 1%ile1140148130
Rarely Seen 1%537022
Rarely Seen 5%161121101915
Pathogens282938313234
Outside Range from JasonH556688
Outside Range from Medivere121219191919
Outside Range from Metagenomics99101066
Outside Range from MyBioma666699
Outside Range from Nirvana/CosmosId202014142121
Outside Range from XenoGene363636363939
Outside Lab Range (+/- 1.96SD)2122416189
Outside Box-Plot-Whiskers67831069410658
Outside Kaltoft-Moldrup641832188710675
Condition Est. Over 99%ile000000
Condition Est. Over 95%ile020000
Condition Est. Over 90%ile035000
Enzymes Over 99%ile0021000
Enzymes Over 95%ile19066151736
Enzymes Over 90%ile68131193427118
Enzymes Under 10%ile302852039420080
Enzymes Under 5%ile13225130418127
Enzymes Under 1%ile1164802211
Compounds Over 99%ile1017000
Compounds Over 95%ile1803531018
Compounds Over 90%ile49573131764
Compounds Under 10%ile78987696511241135998
Compounds Under 5%ile77984892710921057959
Compounds Under 1%ile77383290410691018930
Sep21Sep21Mar22Mar22May22May22Sep22Sep22Jan23Jan23May23May23
PercentileGenus%Genus%Genus%Genus%Genus%Genus%
0 – 974%4626%4323%75%116%96%
10-191911%137%95%1712%2715%2113%
20 – 292615%148%158%1611%148%2314%
30 – 39138%137%168%128%159%1811%
40 – 49148%137%147%1712%1810%159%
50 – 59148%169%147%107%159%1610%
60 – 692213%2011%189%107%1710%1811%
70 – 792313%158%2212%1913%169%138%
80 – 892313%1911%189%1812%2112%159%
90 – 99116%116%2212%1913%2112%138%
Total172180191145175161
Sep21Sep21Mar22Mar22May22May22Sep22Sep22Jan23Jan23May23May23
Percentile%Species%Species%Species%Species%Species%Species
0 – 95%1028%6027%707%1410%228%16
10-1913%265%115%1312%2315%3316%34
20 – 2913%278%186%169%1810%2312%25
30 – 397%154%86%1512%2310%2211%23
40 – 499%198%169%248%158%1810%21
50 – 5912%259%2012%309%1810%2212%25
60 – 698%1712%258%217%148%176%12
70 – 7910%209%197%1911%229%2010%21
80 – 8913%279%197%1916%309%219%19
90 – 997%158%1712%308%1512%287%15
201213257192226211

We lack any data on Lorazepam and other Benzodiazepines impacts on the microbiome which complicates interpretations. I did a search on the US National Library of Medicine and found nothing useful.

More History of Patient

I discovered Pregabalin in March 22 which brought me from Bell CFIDS disability scale below 10 up to nearly bell 20-30. The ditch in the curve around May 22nd is probably me taking too much Lorazepam and Pregabalin, as I for the first time I could take care of things that needed to betaken care of.

At that time I took about 500mg metformin (which did give me energy) most of the day and stayed with my Thorne Fibre mend, Inulin (in the beginning Inulin from the Argave helped dramatically with nausea and headaches)and Acacia Fibre, sometime an Amino Acid complex, but they make me jittery. Usually completely constipated I suddenly developed a strong diarrhea along with an unbearable itching of my whole skin in August / September for which Famotidine(Pepcid) and Cromoglicic acid (Cromolyn – prescription in US) worked best. That was a time where I ate lots of cake and carbohydrates and would take Metformin (I am always hovering around the entrance point to prediabetic) afterwards. That seemed to be too much sugar, my body couldn’t deal with. After I stopped the cake , eat more vegetables again, it went away. 

I have got restless legs, which are kind of turned on or off with every mayor crash. Now unfortunately they are turned on, and the only thing apart from medication that helps is when I eat complex carbohydrates lie brown rice, whole food, pasta, oat flakes etc, when I don’t I use Pramipexole.

Pregabalin been used with Fibromyalgia, a sibling condition for some, and suggested by the American Family Physician journal in 2023. Pregabalin with Lorazepam has known interactions: ” increase side effects such as dizziness, drowsiness, confusion, and difficulty concentrating.”[Src] so she is right about her loss of effectiveness.

Of the many items cited, we know what a few of them likely shifts. Others we lack data.

This missing data illustrates the challenge of trying to manipulate the microbiome — an absence of data. For antibiotics we have a reasonable amount of information, thus we can negotiate with MDs between their desired goal for the antibiotic and our goal of improving the microbiome to find a mutually acceptable compromise.

Going Forward

As part of my learning process, I evaluated each against the “Just Give Me Suggestions” consensus to see it that provide any insight. I also looked at the top items in three other classes.

CriteriaSep21Mar22May22Sep22Jan23May23
cromolyn disodium salt275.5243.9368.9504.3391.1393.3
famotidine275.5248.8378.5504.3379.7393.3
metformin146.4163.8234.1249.9293-17
inulin-231.6-79.1-207.5-203.3-333-157
Total399587725806631491
Best Probioticlactobacillus caseilactobacillus caseilactobacillus caseilactobacillus caseilactobacillus caseilactobacillus casei
Best Amino Acidpolymannuronic acidpolymannuronic acidmelatonin supplementmelatonin supplementmelatonin supplementpolymannuronic acid
Best Vitamin/MineralVitamin B7Vitamin CVitamin CVitamin B1Vitamin B1Vitamin B-12
Great Consistency across the samples!

This helps us evaluate possible (we do not know for certain) impact on various microbiome.

I am not a medical professional and have no clinical experience, so picking items tend to be arbitrary in most cases. I am familiar with the literature for ME/CFS and if the person has ME/CFS, I will tend to pick items that studies reporting helping.

My preference is simple.

My suggestion (given all of the fuzziness and items being taken) is to persist with the prescription items — they help both her symptoms and her microbiome! I would suggest adding the following items (see Dosages for Supplements for literature on dosage):

  • lactobacillus casei – at least 48 BCFU/day — this is the suggested serving size from Custom Probiotics product. Or a Yakult bottle with each meal (each bottle is 20 BCFU). Depending on availability and cost.
  • melatonin – 10 mg/d – in three dosages, i.e. one with each meal.
  • Vitamin B1, B12, C7 and C. (see above for dosages)

One additional item that I would suggest, being prediabetic is to take the Pendulum Akkermansia muciniphila probiotic. This may be a challenge to obtain in Germany (if someone is visiting the US, that may be a backdoor to get it).

Postscript – and Reminder

I am not a licensed medical professional and there are strict laws where I live about “appearing to practice medicine”.  I am safe when it is “academic models” and I keep to the language of science, especially statistics. I am not safe when the explanations have possible overtones of advising a patient instead of presenting data to be evaluated by a medical professional before implementing.

I can compute items to take, those computations do not provide information on rotations etc.

I cannot tell people what they should take or not take. I can inform people items that have better odds of improving their microbiome as a results on numeric calculations. I am a trained experienced statistician with appropriate degrees and professional memberships. All suggestions should be reviewed by your medical professional before starting.

The answers above describe my logic and thinking and is not intended to give advice to this person or any one. Always review with your knowledgeable medical professional.

The Transcribed Tests – a New Option: Condition Matching

As a result of doing an analysis for a 19 month old toddler, I added a new option that can also be used with Transcribed tests. This post applied to the following tests:

You must save your input.

Process

Here’s an example

When you logged in, you will see your saved tests, CLICK ON Review.

And then we have the details you entered below with an important column, taxon number.

Below this are conditions where your pattern matches at least 5 shifts reported in Published Studies.

There may be many items listed. This is by pattern matching and is not predictive.

If you have any of these conditions, or suspect you may have. Just click the appropriate button.

An example is below. These are tuned safest-suggestions for the matches. What do I mean by safest? It means the items are not reported in any study in the database to adversely impact any of bacteria listed. Many substances have contradictory reports on shifts — this substances are excluded.

Not Listed Condition?

This person believes they may have Autoimmune, so going to https://microbiomeprescription.com/Library/PubMed we find that it is listed.

If it is not listed, search for bacteria shifts reported and use those (please send me the studies so I may add them).

The bacteria are shown in a tree. You have to manually match between the two.

In this we have:

  • Escherichia     ⬇️  but our sample is high,
  • Roseburia intestinalis   ⬆️  - we are high on Roseburia, we will include it

We have only one match — this tests with limited reporting is not a good fit for this condition. Doing a test like Biomesight, Xenogene, Thorne or Ombre is likely the best choice.

We just copy the taxon number into the form at the bottom of the page, and then click suggestions.

In this case, we get a short list. Remember, doing a single bacteria means you are ignoring a lot of interactions and factors. The suggestions could feed other bacteria that are too high.