Plasma and Fecal Metabolite
Profiles
in
Autism Spectrum
Disorder
Brittany D.
Needham
1
,
Mark D.
Adame
1
,
Gloria
Serena
2
,
Destanie
R.
Rose
3
,4
,
Gregory M.
Preston
5
,
Mary C.
Conrad
5
,
A.
Stewart
Campbell
5
,
David H.
Donabedian
5
,
Alessio
F
asano
2
,
Paul
Ashwood
3
,4
,
and Sarkis K.
Mazmanian
1
1
Biology and Biological Engineering, California Institute of Technology, Pasadena
,
CA
,
91125, USA
2
Division of Pediatric Gastroenterology and Nutrition,
Mucosal Immunology and Biology Research
Center
, Massachusetts General Hospital for Children, Boston, MA, 02114, USA
3
Department of Medical Microbiology and Immunology,
University of California Davis
, Da
vis, CA,
95616, USA
4
The M.I.N.D. Institute, University of California, Davis, Sacramento, CA, 95817, USA
5
Axial Biotherapeutics
,
Waltham
, MA, 02
451
, USA
Correspondence:
bneedham@caltech.edu
and
sarkis@caltech.edu
Running title:
Plasma
and Fecal Metabolite Profiles in ASD
Keywords:
Autism Spectrum Disorder, ASD, metabolomics,
plasma metabolites, fecal metabolites,
steroid hormones, mitochondrial dysfunction,
phenolic metabolites
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
2
ABSTRACT
Autism Spectrum Disorder (ASD) is a neurodevelopmental condition with hallmark behavior
al
manifestations
including
impaired
social communication
and
restricted
repetitive behavior.
In
addition,
many
affected
individuals display
metabolic imbalances, immune dysregulation
,
gastrointestinal (GI) dysfunction
, and altered gut microbiome compositions
. We sought to
better
understand non
-
behavioral
features
of ASD by
determining
molecular
signatures
in
peripheral
tiss
ues
. Herein, we present
the
untargeted
metabolom
e
of 231 plasma and 97
fecal
samples from
a
large
cohort of
children with ASD and
typically developing
(
TD
)
controls.
D
ifferences in lipid,
amino acid, and xenobiotic metabolism
discriminate
ASD and
TD
sample
s.
We
reveal
correlations between
specific
metabolite profiles and clinical behavior score
s
,
and
identify
metabolites
particularly
associated with
GI
dysfunction in ASD.
These findings
support a
connection between GI
physiology
, metabolism
,
and
complex
behavior
al traits
, and may advance
discovery and development of molecular biomarkers for ASD.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
3
INTRODUCTION
Many
diseases
are associated with
informative
metabolic signatures
,
or
biomarkers
,
that
enable
diagnoses, predict disease course, and
guide
treatment
strategies
.
In contrast, a
utism spectrum
disorder (ASD) is diagnosed based on
observational
evaluation
of
behavioral symptoms,
including reduced social interaction and repetitive/stereotyped behaviors
(
1
)
.
The average age of
ASD
diagnos
is is
between
3
-
4 years
old
(
2
)
, at which time children can receive behavioral
therapy
,
the gold standard
treatment.
Because
earlier diagnosis
imp
roves
efficacy of behavioral
therapies
(
3
,
4
)
, molecular biomarkers represent an attractive approach
for
identif
ying
‘at
-
risk’
populations
and may
aid
deve
lopment of
personalized
therap
ies
.
This prospect
is
increasingly
important given the rising rate of
ASD
diagnoses
, which currently stands at
up to
1 in
59
children
in the United States
(
2
)
, with no FDA approved drugs for core
behavioral
symptoms
.
Metabolic
abnormalities have been
reported
in ASD
(
5
)
, though most have measured onl
y
a
small subset of metabolites and many outcomes have not reproduced between cohorts.
Mitochondrial d
isease, which heavily
influences systemic metabolism,
is
estimated to be higher
in ASD compared to controls (5% vs ~0.01%)
(
6
)
,
and
genes crucial for
mitochondrial
function
are risk factors for ASD in humans and
rodent models
(
7
)
.
The metabolic abnormalities
associated with mitochondrial dysfunction
in ASD
affect cellular energy, oxidative stress, and
neurotransmission in the gut and the brain
(
7
–
19
)
.
Othe
r
metaboli
c
profiles
in
ASD
implicate
a
romatic and phenolic metabolites
,
including
derivatives of
nicotinic, amino
acid, and hippurate
metabolism
(
20
–
30
)
.
Various amino acids are detected at differential levels across studies and
across sample types, but any consistent
patterns are difficult to
discern
(
15
,
24
,
26
,
29
,
31
–
34
)
.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
4
Some of the discrepancies between these studies are likely due to differences in
sample number
,
tissue analyzed
,
and methodology.
Other sources of variability include differences in
environmental
factors, such as
diet and gut
bacteria
,
which
differ
between ASD and TD
populations
(
35
–
37
)
. Diet is a major source of
circulating
metabolites
, impacting
the metabolome
directly or indirectly
through
chemical
transform
ation
by t
he
trillions of
gut
microbes
,
the
microbio
me
,
which
has been proposed
to
modulate complex behaviors in animal models and
humans
(
38
–
40
)
.
Such p
roposed environmental modulators of ASD
may integrate with genetic
risks to impact behavioral endpoints through the actions of small molecules produc
ed in
peripheral tissues outside the brain.
Herein
,
we
present
a
comprehensive
comparison of
an extensive
panel of
identified
metabolites
in
human
plasma an
d
feces
from a large cohort of
matched
ASD and
TD
children
.
We identified
differential
levels
of
metabolites ranging from
hormone
s
, amino acid
s
, xenobiotic
s
, and l
ipid
s
,
many of which
correlate
with clinical
behavior
and GI
scores
.
To our knowledge, this is the first
study to concurrently evaluate paired intestinal and systemic metabolomes in a high
-
powered
analysis with a large number of identified metabolites, allowing direct associations between
metabolites previously highlighted
in ASD samples
and
discovering new metabolites of interest.
These findings
support the emerging concept of
evaluating non
-
behavioral features
in the
diagnosis of ASD
and its GI comorbidities
.
METHODS AND MATERIALS
Participants:
Samples for this study
, age
d 3
-
12 years old,
were collected through the UC Davis
MIND institute
(
41
,
42
)
. ASD diagnosis
was confirmed at the MIND Institute by trained staff
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
5
using the Autism Diagnostic Observation Schedule (ADOS), and the Autism Diagnostic
Interview
-
Revised (ADI
-
R).
S
ubjects were diagnosed prior to 2013 based on DSM IV. Typically
developing (TD) participant
s were screened using the Social Communication Questionnaire
(SCQ). Participants in the TD group
score
d
within the typical range, i.e. below 15, on the SCQ
and above 70 on the
Mullen Scales of Early Learning (MSEL) and Vineland Adaptive Behavior
Score (VAB
S). Ninety
-
seven of the participants, who also provided the stool samples, completed
an additional evaluation to determine gastrointestinal (GI) symptoms. GI status was determined
using the GI symptom survey (GISS), based upon Rome III Diagnostic Questionn
aire for the
Pediatric Functional GI Disorders
(
43
)
and described in detail elsewhere
(
42
)
.
S
ee Supplemental materials
for additional methods
.
RESULTS
Plasma and Fecal Metabolomes Differ between ASD and
TD
Plasma samples from 130 ASD and 101
TD
children
were analyzed along with fecal samples
collected from a subset of these same ASD (n=57) and
TD
individuals
(n=40) (Fig
ures
S
1A
-
S
1G
)
.
S
amples
and metrics of behavioral and GI scores
were
obtained from the
UC Davis MIND
institute
(
41
,
42
)
.
T
he ASD group
was stratified
into subsets of children with
GI symptoms
(ASD+GI) or without
GI symptoms
(ASD
-
GI)
,
to
explore
potential
effects
of
comorbid
intestinal
dysfunction in ASD
(
40 out of 130 ASD
samples were +GI)
. This stratification was based on
symptoms associated with ASD including diarrhea, constipation, and irritable bowel syndrome
-
like symptoms
.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
6
S
amples were analyzed by the globa
l metabolite panel (
plasma
and fecal samples) and the
complex lipid panel
(
plasma
samples
) by Metabolon
, Inc.
(Durham, NC)
,
which
identified
a total
panel of 1,611
plasma
and 814 fecal metabolites
(Table
s
S
1
-
S
2)
.
Overall, we
discovered that
310
metabolites
were differentially abundant between the ASD and
TD
groups
in plasma
and 112
metabolites
were
differentially abundant in feces
(Fig
ure
s
1
A
-
1
B
,
S1C, and
Tables
S
1
-
S
2
)
.
Using a quantitative assay for a targeted panel of metabolites
,
we observed a
high correlation
between relative abundance and
precise
concentration
(Fig
ure S
2A
)
.
Overall, t
hese data
expand
on previous evidence
that the metabolomic profile
of
ASD
and
TD
population
s
display
differences
not only in the gut compartment, but also in
circulation
, which
may
affect the levels
of metabolites
throughout the body,
in
cluding
the brain
(
44
,
45
)
.
To
appreciate the biological relevance of
differen
t
metabolomes
between ASD and
TD
,
individual metabolite
s
were integrated
into
biochemical pathways
for
pathway enrichment
analysis
, revealing
the
degree of change
within each
.
Here, we identified large scale changes
,
mostly
in
lipid, xenobiotic
,
and nucleoti
de pathways
associated with
diverse
physiological
process
es
(Fig
ures
1
C
-
1
D
,
S2
B
)
.
We then used
Random Forest
machine learning analysis
to
determine if
metabolite profiles can unbiasedly predict whether the sample came from an
ASD
and TD
donor
. Overall, the modest predictive accuracy of this machine learning approach was
68% for plasma and 67% for feces.
T
o test whether
focusing on
the
most
discriminating
metabolites
would
improve th
e
prediction
, we repeated the Random Forest analysis using the
top
30 metabolites
and found that the predictive accuracy
for plasma
improved slightly
to 69%
when
using all ASD samples
,
and
to 73%
when using only ASD
-
GI samples
.
For fecal samples, the
predictive accuracy improved
to 75% using all ASD and to 73%
using only ASD
-
GI.
The top 30
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
7
metabolites
,
calculated b
y measuring the mean decrease in accuracy of
the
machine learning
algorithm
,
are useful to describe the
strongest
drivers of overall metabolic differences between
ASD and TD populations. These most dis
criminatory metabolites were
primarily
from the
lipid,
amino acid,
xenobiotic, and
cofacto
r
/
vitamin
super pathways
(Fig
ures
1
E
-
1
F
)
.
Several
of these
metabolites
have been previou
s
ly linked
to ASD
,
such as
steroids, bile acids,
acyl
-
carnitines
, and
nicotinamide
metabolites
(
46
–
49
)
.
Further, m
ultiple molecules known to be produced or
manipulated by the gut microbiota also feature
d
prominently, including 4
-
ethylphenyl sulfate
(4EPS
),
which is elevated in a
n ASD
mouse model,
and
indolelactate
, a
microbe
-
derived
trypt
ophan metabolite
(Fig
ure
1
E
)
(
50
,
51
)
.
The
two most discriminatory molecules in plasma
(Figures
1
G
-
1
H
)
and
feces
(Figures
1
I
-
1
J
) are
depicted
.
M
etabolites
correlating
the strongest
with these discriminatory metabolites are closely related on a structural and metabolic level
(Figure
s
1K
-
1L).
Global
Metabolite Levels Correlate with
Clinical Behavioral Scores
Using
clinical
metadata
for
ASD
individual
s
, we
correlate
d
the levels of
individual
metabolite
s
to the verbal, social
,
and nonverbal scores of
standard diagnostic measures
:
the
Autism
Diagnostic Interview, Revised (
ADIR
)
,
a parent questionnaire,
and
the cumulative A
utism
D
iagnostic
O
bservation
S
chedule
severity score
(
ADOS
-
SS)
, conducted by trained health
professionals
(
1
)
(
Table
s
S1
-
S
2
)
.
Next
,
group
ed
correlations between
the
behavior metrics and
entire
metabolite pathways
were calculated
(Figure
1
M
).
We
found that v
erbal and social scores
primarily correlate with lipid
metabolism
pathways
and
that
n
onverbal
scores have
the fewest
correlations
.
T
he
ADOS
-
SS correlated
with
diverse metabolite pathways, including
amino acid
s
and
food
/
plant component
pathways
, which may be
partly
due to
the
diverse array of symptoms
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
8
integrated into
the ADOS
-
SS
score
.
Additionally, when stratified by severity of behavior
al
score
s
,
we observe modest
clustering according to
plasma
metabolites correlated to the dis
order
between the
top and bottom quartiles of behavior severity groups for verbal, social, and ADOSS
-
SS metrics, but not for nonverbal
, and less so with fecal samples
(Figure
s
1
N
,
S1
H
-
S1
I
).
These
correlations
support
the involvement of lipid, amino acid, a
nd xenobiotic metabolism in the
etiology of ASD,
as previously described
(
8
,
48
,
52
)
, and
reveal new candidates for
ASD
biomarkers
that correlate with symptom severity
.
Transfer of ASD Fecal Microbiota into Mice is
Accompanied by Metabolic Signatures
Since microbial metabolites ranked highly in the
Random Forest
machine learning analysis
, w
e
tested if any of the observed metabolite differences in humans could be transferred to mice via
fecal microbial transplant. We
selected 4 male donor samples from
each of
the ASD and
TD
group
s
and colonized
2
-
3 male
germ
-
free
mice per donor for three weeks
before collecting
plasma and fecal samples for metabolite profiling and bacterial DNA sequencing
(Figure S2C
-
S2D)
, respectively
.
Global metabolomic analysis revealed
that
c
olonized mice
modestly
cluster
by donor and
group
when statistically significant metabolites are considered in PCA analysis
(Figure
S2E
).
We
selected the human
donors based on 4EPS
levels
(Figure S2
F
)
,
due to
its
involvement in an ASD mouse model
(
50
)
and
dysregulation of similar phenolic compounds in
human ASD
(
22
,
24
,
28
,
29
,
53
)
.
4EPS is not produced by the host and is strictly a bacterial
metabolite
(
50
,
54
)
. Surprisingly, we observed
4EPS
levels
in
a bimodal distribution in mouse
samples (Figure
S2
G
-
S2H
).
In spite of the surprising results with 4EPS,
many of t
he metabolites
with the highest fold change and lowest
p
-
value are
indeed
other
phenolic molecules such as
metabolites of hippurate, tyrosine, and diet
-
derived phenols
(Figure
s
S2
I
,
Table S5).
W
hile
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
9
preliminary,
these studies
reveal
microbiome
-
mediated
effects on xenobiotic pathways and
phenolic
metabolites
dysregulated in ASD.
Lipid and Xenobiotic Metabolites
Correlate
with GI Symptoms
Considering the prevalence of intestinal issues in ASD,
GI
-
dependent
m
etabolite
analysis
between
individuals
may
provide
useful
insights
.
Accordingly,
we found a total of 87 and 2
4
metabolites in plasma and feces, respectively, that were differentially abundant within
the
ASD+GI compared to
ASD
-
GI
individuals
(Figure
s
2
A
-
2
B
).
At the pathway level
, we
found
differences
in
free fatty acid
and xenobiotic metabolites
in the food component and plant
pathways
between ASD+GI vs ASD
-
GI
plasma
(Figure
2
C
)
,
with
no broad
GI
-
dependent
pathway alterations in fecal samples
(
Table S2).
In plasma,
free
fatty
acids
of
multiple chain
length
s
we
re lower in ASD+GI
compared to ASD
-
GI samples
, including monounsaturated,
saturated, and polyunsaturated fatty acids (PUFAs)
(Figure
2
D
).
PUFA
s
are anti
-
inflammatory
and lower levels of these fatty acids may contribute to
GI dysfunction directly
(
12
)
.
Several metabolites from the food component and plant pathway also
discriminate ASD+GI from
ASD
-
GI
, such as
lower circulating
piperine metabolites
in
ASD+GI plasma samples
(
Figure
s
2
E
-
2
G
)
.
A
lower level of piperine
metabolites
was
also
associ
ated with
higher
ADOS
-
SS
(Fig
ure
2
H
) and worse nonverbal scores
(
Table
S
1
).
Observed alterations to piperine metabolite levels
are
supported by the fact that
oral administration of
piperine has been
successfully
used as a
treatment in preclini
cal ASD
model
s
, presumably
due to
its antioxidant quality
(
55
)
.
While these
correlations are interesting,
the implications of
correlating an
altered metabolome and GI
symptoms in ASD remain to be determined.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
10
Steroid H
ormone Levels
are Elevated in ASD
Multiple
human
ASD
studies
have examined the levels of specific steroid metabolites within
androgenic, pregnenolone and progesterone metabolism, with some finding
aberrant
levels,
positive correlations with ASD
severity
, and behavioral improvement
following
treatment
s
th
at
lower levels
of certain
hormon
es
(
47
,
56
–
64
)
.
On the other hand, in a recent clinical t
rial
of
ASD
childre
n
given an
antioxidant treatment, levels of
pregnenolones and androgens increased and
correlated with improved behavior
(
65
)
.
In our dataset
, we found
robust
increase
s
of
al
most
every
detected
metabolite within the
pregnenolone,
androgen, progesterone
,
and corticosteroid
pathways
in
the
plasma of ASD children (Fig
ures
1
C
,
3
A
,
Table
S
1
).
Similarly,
some
androgenic
steroid pathway
metabolites
we
re elevated in ASD fecal samples (Fig
ures
3
A
,
Table S2
). This is
a strong indicator that the physiological pathways associated with
the downstream metabolism of
chole
s
terol
are significantly altered between ASD and
TD
populations
(Figure
3
B
)
.
There does
not appear to be a global change in steroid meta
bolism, as most primary bile acid
and sterol
me
tabolites
we
re unaffected (
Table
s
S
1
-
S
2
).
W
e observe
d
some
elevat
ion of
these hormone
levels
independent of sex
, which is notable considering the male
bias in ASD
, reflected in our
primarily male sample
set
(7
-
15% female)
(
66
)
(
Fig
ure S3
A
-
S3B
).
Because a cluster of
our
samples
are from older individuals
in the ASD group, and
to account for
age
-
dependent
increases
in androgen
s
,
we stratified by age
and
still
observed
heightened
androgenic and pregnenolone
metabolite
levels in ASD
sub
populations
(
Figure
S3
C
).
Taken together, these data indicate that
steroidal hormone
metabolism
may be
altered in the ASD population relative to
TD
samples and
that these differences are n
ot driven
sole
l
y
by sex or age differences in our cohort.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
11
Lipid Metabolite Levels
Differ
in ASD
Lipids are crucial for energy storage, cellular membrane integrity and cell signaling. They play a
variety of roles in the central nervous system, and their
dysregulation
has been linked to
ASD
(
9
,
11
)
.
P
hospholipids have been measured at lower levels, while long chain fatty acids
are
reportedly
elevated but PUFAs have been measured at both higher and lower levels
,
depending
on the cohort
(
9
,
12
,
13
)
.
Here,
we performed an
untargeted
, quantitative
metabolite
analysis on
complex lipid
s
.
T
he concentration
s
of
999 lipids were quantified
, and 13.7% of these were
significantly
different
in the ASD samples, with 4.6% increased and 9.1% decreased compared to
TD
samples.
Many of these differentially abundant lipids included phospholipids, cholesterol
esters and glycerolipids.
In general, shorter (14
-
18 chain length) saturated fatty acids wer
e
less
abundant
in ASD throughout the lipid classes
(Figures
3
C, S
4
A
-
S
4
B
)
.
PUFA
lipid
levels were
also
different
, with elevations
in diacylglyce
rols and free fatty acids, and
an enrichment
of the
18:2
(linolenic)
chain length in most lipid classes
(Figure
3
C)
.
F
ecal
samples generally
trended in the opposite direction from plasma
in PUFA lipids
(Figure S4
C
,
Table S2
).
Intriguingly, m
ultiple lipids with
linolenic
and
linoleic
(
18:3
)
chains we
re correla
ted
with social behavior (Tables
S1
and
S2
).
PUFA lipids containing linolenic and linoleic acids are
precursors to the important PUFAs (arachidonic 20:4 and docosahexaenoic 22:6 acids) for brain
development,
function,
and
structural i
n
tegrity
(
12
)
.
Future studies are needed to determine if
these specific changes in lipid
levels
contribute to ASD symptoms.
ASD Correlates with
Cellular Energy and Oxidative Stress Metabolite
s
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
12
Many lipids are
markers of
mitochondrial and oxidative stress
,
offering a snapshot into cellular
metabolic state
s
. These markers include
acyl
-
carnitines
,
which
have been
highlighted
in
various
ASD studies
and are established indicators of
mitochondrial dysfunction
(
11
,
12
,
14
–
18
,
24
,
25
,
33
,
46
,
58
,
67
–
71
)
.
A
cyl
-
carnitines
are formed
to allow transport
of lipids
across the
mitochondrial membranes for beta
-
oxidation
, and abnormal levels of these conjugated lipids
accumulate to higher levels with a
decrease
in beta
-
oxidation
.
Interestingly, high levels of short
acyl
-
carnitines are found in
rodent
mode
ls
where
ASD
-
l
ike
behavior
s
are
induced with the short
chain fatty acid
s
,
valproic and propionic acids
(
72
)
.
We fou
nd
differential
lev
els of various acyl
-
carnitines
in ASD, creating a pattern of more
abundant short chain acyl
-
carnitines and less abundant long chain acyl
-
carnitines in
the ASD
-
GI
samples
compared to
TD
samples
(Fig
ure
4
A
)
.
Acyl
-
carnitines a
re positively
correlated with
more severe social
defects
,
an
effect
driven by
stru
c
tures
with
s
horter moieties (C2
-
C14) (Fig
ure
1
M
,
Table
S
1
)
. In fecal samples, a
cetylcarnitine (C2) and free car
nitine we
re elevated in ASD
(Figures
4
B
-
4
C)
and
were highly
discriminatory
(Figure
1
F
).
Other mitochondrial markers
in
both plasma and feces
we
re
also
differentially abundant
in
ASD
and are summarized in
Figure
4
D along with
markers of
phospholipid
metabolism
, which occurs largely in the mitochondria
and
was
significant
ly altered
in
fecal
samples
(Fig
ures
4
D
,
1
D
).
Additionally
, the 5
plasma
metabolites most
positively correlated with ADOS
-
SS are all involved in
these
cellular energy
pathways
(Figure
S4D
,
Table
S
1
). These observed differences in energy markers and lipids could
have a
neurodevelopmental effect during
periods
when
the
high lipid
and energy requirement
in
the brain is crucial
(
73
–
75
)
, and alterations to levels of tricarboxylic acid (TCA) cycle
intermediates have been observed in human ASD prefrontal cortex samp
les
(
49
)
.
Such d
ef
ects in
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
13
cellular metabolism
su
pport the theory
that mitochondrial dysfunction
may not only be comorbid
with ASD
but also a potential
contributing factor, as suggested by numerous previous reports
(
6
,
13
,
67
,
69
,
72
)
.
A
mino acid
signatures
of
oxidative stress and mitochondria
l function
are
also
present
in our
dataset. Dysregulated amino acid degradation, homeostasis, and
import into the
brain ha
ve
been
implicated as a cause of neuronal stress in ASD, and supporting metab
olomic data has shown
perturbations of various amino acid pathways
, such as glutamate, methionine, glutathione, and
gamma
-
glutamyl metabolites
(
10
,
17
,
24
,
26
,
29
,
33
,
49
,
52
,
76
)
, which exhibit
differences as
well (Figure 4E)
.
We
also demonstrate correlations between pathways
of
oxidative stress
(
cysteine, methionine, SAM and glutathione
pathways)
with
the ADOS
-
SS
(Figure
4
F)
.
S
ome of
these molecules we
re found in higher levels in the ASD feces and at lower levels in the ASD
plasm
a
, such as hypotaurine
(Figures
4
E,
4
G
-
4
H
), which might indicate altered
fecal
production
,
excretion
or differential uptake into the plasma potentially through vari
ed
intestinal permeability.
Similar to
hypotaurine,
levels of its precursor,
taurine
,
were significantly increased in ASD fecal
samples, although not altered in plasma (Figure S4E and S4F).
Taurine
plays
many roles
throughout the host, and
has previously been measured at altered levels in ASD, although with
little consensus
(
17
,
18
,
24
–
26
,
33
,
65
,
71
,
77
)
.
Hypota
urine and taurine
deficiency has
been
shown to lead to defects in
cell
differentiation
in the brain
(
77
)
and their dysregulation could
alter neuronal signali
ng
(
78
)
.
O
xidative stress
-
related
glutathione pathway
precursors, gamma
-
glutamyl amino acids
,
are
relevant to
ASD
through their
influenc
e on
levels of neurotransmitters such as gamma
-
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
14
aminobutyric acid (GABA)
(
78
,
79
)
, which is
widely thought to play a role in ASD
. In a recent
ASD study,
an experimental treatment that led to increased
gamma
-
glutamyl AAs and other
redox pathway metabolites correlated with improved behavior metrics
in children
(
65
)
.
A
lmost
every gamma
-
glutamyl amino acid is negatively correlated with ADOS
-
SS
(Fig
ure
4
I
).
We
also
observe
perturbations
in the urea cycle, which process
es
the amino
group of amino acids for
excretion in urine
(Figure S4G
-
S4H
)
.
Abnormalities in this pathway can
be indicative of altered
amino acid degradation
observed in ASD
and can
lead to
neuro
tox
ic accumulation of nitrogen
-
containing compounds in the bloo
d
(
80
)
.
Together, these
results corroborate and extend a
growing body of research into altered mitochondrial metabolism and oxidative stress in ASD.
Differential Phenolic Xenobiotic Metabolite Levels
in ASD
P
henolic m
etabolites
, a diverse structural class
comprised of
thousand
s of molecules
containing
a
phenol
moiety
,
c
ome from dietary ingredients or from biotransformation of aromatic amino acids
by the gut microbiota. In their free phenolic forms they can be readily absorbed through
intestinal tissues; however, microbial modif
ication of these molecules can significant
ly
alter
their
absorption, bioavailability, and bioactivity
(
81
)
, leading to various benefits or harm to the host
(
81
,
82
)
.
In fact, altered levels of phenolic molecules
have
been
hig
hlighted
in
many
ASD
metabolomic studie
s
(
18
,
20
–
26
,
28
,
31
,
32
,
53
,
65
,
83
,
84
)
.
However, a consensus of enrichment
or depletion o
f these metabolites across ASD vs
.
TD
groups has yet to be reached, and most
studies
have
only
measured
a
small
subset of this structural class
of metabolites
.
W
e
observe
d
altered levels of
phenolic metabolite
s
belonging to several
interrelated pathways,
including tyrosine, benzoate, and food component and plant metabolites
(Fig
ure
5
A
).
S
ome
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
15
molecules
,
such as
homovanillate and tyramine,
are involved in neurotransmitter metabolism
(
18
)
.
Oth
er
s
include
hippurate derivatives
,
which we
measured
at
lower
levels
in the ASD
plasma
samples
(Figure
5
A)
.
Phytoestrogens such as daidzein, genistein, and equol derivatives
are elevated
in our ASD plasma samples
(Figure
5
A)
.
These metabolites
have purported health
benefits but are also known disruptors of endocrine signaling
(
85
)
.
Additionally,
we detected
altered levels of
molecules
with
structural similarity to
para
-
cresol sulfate,
a toxic molecule
elevated in the urine of young ASD children
(
20
,
22
)
. These
include
, among others,
4
-
ethylphenyl sulfate
(
4EPS
)
, 2
-
ethylphenylsulfate,
cresol derivatives, 4
-
allylpheny
l sulfate,
and 4
-
methylbenze
ne
sulfonate
, the latter of which was
elevated
a remarkable
60
-
fold
in a small subset
of ASD samples
(
Figure 5B
-
5D,
Table
S
1)
.
Some
changes
in the levels of
these
phenolic
molecules
are
also
observed
in animal models of ASD
(
50
,
86
)
. P
reviously
,
4EPS was
observed
at a
46
-
fold
elevated level
in the maternal immune activation mouse model, and daily
administration of synthetic 4EPS to wild type mice was sufficient to induce an anxiety
-
like
phenotype
(
50
)
.
Here, i
n ASD plasma samples,
4EPS
levels
were
increased
6.9
-
fold
(
Figure
5
D)
.
Other phenolic metabolite levels correlate with 4EPS in ASD samples, including
4
-
acetylpheny
l
sulfate,
a
derivative of 4EPS
, and others
.
(Fig
ure
5
E
, S5A
-
S5B
)
.
These observed
alterations in
phenolic molecules
, combined with
mounting evidence from
previous
studies
,
suggest that
phenolic
structural metabolites
may
play a role in ASD
.
DISCUSSION
Changes in the metabolome have been linked to a number of
neurodevelopmental and
neurodegenerative disorders
(
87
,
88
)
. The
current
study
includes
a
comprehensive
profiling of
the
metabolites from
23
1
plasma and
97 matched
fecal samples from
ASD and
TD
individuals
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
16
who have been extensively behavior tested
.
We observed that
specific
individual
metabolites and
metabolic
pathways in
lipid, amino acid, and xenobiotic metabol
ism
were altered between ASD
and
TD
groups
,
and that many of these same
molecules and
pathways correlated to severity of
ASD
behavioral
scores
.
The core findings are summarized in Figure 6.
Various
s
teroid hormone
metabolite levels
were elevated
in ASD sample
s
.
PUFA levels, and short chain acyl
-
carnitines
were generally elevated in ASD plasma, while saturated fatty acid levels and long
-
chain acyl
-
carnitines were decreased,
contributing to
a
general
picture of dysregulated cellular energy and
oxidative state
, with potential connections to
mitochondrial dysfunction in ASD.
Some of the
lipid and xenobiotic metabolites also showed interesting changes w
ithin the ASD group when
stratified by the presence of GI symptoms.
Finally, m
any phenolic metabolites, derived largely
from host and bacterial metabolism of amino acids, plant polyphenols, and other food
components were
detected
at differential levels in
plasma and feces
between the comparison
groups
.
ASD
is
diagnosed by behavioral tests,
with
extensive
heterogeneity of symptoms, severity, and
etiology
between individuals
and little consensus on molecular mechanisms
.
This enigmatic
spectrum
has a strong
but complex
genetic basis, with
hundreds
of
reported
risk genes
(
89
)
.
The
contributions of risk alleles to behavior is a vast and active area of research
. In addition to
gene
tics
, understanding
of
altered metabolite levels in blood, feces, and brain
s
of ASD
individuals
may prov
ide a glimpse into physiologic aspects of the disorder and hold the potential
to advance diagnosis and/or stratification of sub
-
populations of ASD.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
17
Interestingly,
m
ultiple mouse models support the notion that gut metabolites
are
associated with
brain development and
function
(
48
,
90
,
91
)
. Some metabolites are closely linked with
neurological disorders, either
to positive or negative outcomes
(
92
,
93
)
.
Several
exa
mples have
been
reported
of gut metabolites entering circulation and directly affecting the brain
,
as well as
cases where metabolites stimulate pathways in the gut, immune, or autonomic nervous system
and exert changes to
the
brain and
to
behavior
(
94
–
97
)
.
Comprehensive metabolic profiling in
human
s
and animal model
s
provide
s
insight into the molecular status of disease and how genetic
factors and environmental
risks
interact.
Deeper
analysis of
our
dataset
along with
additional
studies
,
with future
empirical studies to validate the relevance of
our
observations, could
illumin
ate
aspects of
ASD
pathophysiology
. The
intriguing
correlations between ASD behavior
s
,
altered levels of
fecal and plasma
metabolites,
and
GI symptoms contribute to the concept that
ASD
may be viewed as a
whole
-
body
condition, and argue for increased inves
tigation into
peripheral aspects of disease that may lead to advances in diagnosis and improved stratification
of ASD populations.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
18
ACKNOWLEDGMENTS
This work was supported by
grants from
Autism Speaks
(Grant #7567 to P.A., A.F. and S.K.M.)
;
the Johnson Foundation (to P.A
.)
;
the Brain Foundation (to P.A
.)
;
Genome, Environment,
Microbiome, and Metabolome in Autism (GEMMA)
(#
825033 to AF
)
; Axial
Biotherapeutics
(to
S.K.M.)
; and
the National Institutes of Health (
HD090214 to P.A.
) and (
MH100556
to S.K.M.)
.
We would like to thank the MIND Institute study staff, and the dedication and commitment of
the families who took part in these studies is gratefully acknowledged.
Declaration of Interests
A.S.C., D.H.D
.
and S.K.M.
ha
ve
financial interest
in Axial Biotherapeutics.
A.F. has financial
interest in Alba Therapeutics.
G.M.P. and M.C.C. are employed by Axial Biotherapeutics.
B.D.N., M.D.A., G.S., D.R.R., and P.A. report no financial conflicts of interest.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
19
REFERENCES
1.
E. Doernberg, E. Hollander, Neurodevelopmental Disorders (ASD and ADHD): DSM
-
5,
ICD
-
10, and ICD
-
11.
CNS Spectr.
21
, 295
–
299 (2016).
2.
J. Baio, L. Wiggins, D. L. Christensen,
M. J. Maenner, J. Daniels, Z. Warren, M. Kurzius
-
Spencer, W. Zahorodny, C. Robinson Rosenberg, T. White, M. S. Durkin, P. Imm, L.
Nikolaou, M. Yeargin
-
Allsopp, L.
-
C. Lee, R. Harrington, M. Lopez, R. T. Fitzgerald, A.
Hewitt, S. Pettygrove, J. N. Constantin
o, A. Vehorn, J. Shenouda, J. Hall
-
Lande, K. Van
Naarden Braun, N. F. Dowling, Prevalence of Autism Spectrum Disorder Among Children
Aged 8 Years
—
Autism and Developmental Disabilities Monitoring Network, 11 Sites,
United States, 2014.
MMWR Surveill. Summ
.
67
, 1
–
23 (2018).
3.
G. Dawson, S. Rogers, J. Munson, M. Smith, J. Winter, J. Greenson, A. Donaldson, J.
Varley, Randomized, controlled trial of an intervention for toddlers with autism: the Early
Start Denver Model.
Pediatrics
.
125
, e17
-
23 (2010).
4.
J
. Virués
-
Ortega, Applied behavior analytic intervention for autism in early childhood: meta
-
analysis, meta
-
regression and dose
-
response meta
-
analysis of multiple outcomes.
Clin.
Psychol. Rev.
30
, 387
–
399 (2010).
5.
A. Azhari, F. Azizan, G. Esposito, A sys
tematic review of gut
-
immune
-
brain mechanisms in
Autism Spectrum Disorder.
Dev. Psychobiol.
(2018), doi:10.1002/dev.21803.
6.
D. A. Rossignol, R. E. Frye, Mitochondrial dysfunction in autism spectrum disorders: a
systematic review and meta
-
analysis.
Mol.
Psychiatry
.
17
, 290
–
314 (2012).
7.
N. Cheng, J. M. Rho, S. A. Masino, Metabolic Dysfunction Underlying Autism Spectrum
Disorder and Potential Treatment Approaches.
Front. Mol. Neurosci.
10
(2017),
doi:10.3389/fnmol.2017.00034.
8.
A. El
-
Ansary, L. Al
-
Ayad
hi, Lipid mediators in plasma of autism spectrum disorders.
Lipids
Health Dis.
11
, 160 (2012).
9.
A. K. El
-
Ansary, A. G. B. Bacha, L. Y. Al
-
Ayahdi, Impaired plasma phospholipids and
relative amounts of essential polyunsaturated fatty acids in autistic pat
ients from Saudi
Arabia.
Lipids Health Dis.
10
, 63 (2011).
10.
S. J. James, P. Cutler, S. Melnyk, S. Jernigan, L. Janak, D. W. Gaylor, J. A. Neubrander,
Metabolic biomarkers of increased oxidative stress and impaired methylation capacity in
children with
autism.
Am. J. Clin. Nutr.
80
, 1611
–
1617 (2004).
11.
Q.
-
Q. Lv, C. You, X.
-
B. Zou, H.
-
Z. Deng, Acyl
-
carnitine, C5DC, and C26 as potential
biomarkers for diagnosis of autism spectrum disorder in children.
Psychiatry Res.
267
,
277
–
280 (2018).
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
20
12.
G. A. Most
afa, L. Y. AL
-
Ayadhi, Reduced levels of plasma polyunsaturated fatty acids and
serum carnitine in autistic children: relation to gastrointestinal manifestations.
Behav.
Brain Funct. BBF
.
11
(2015), doi:10.1186/s12993
-
014
-
0048
-
2.
13.
E. Pastural, S. Ritchi
e, Y. Lu, W. Jin, A. Kavianpour, K. Khine Su
-
Myat, D. Heath, P. L.
Wood, M. Fisk, D. B. Goodenowe, Novel plasma phospholipid biomarkers of autism:
mitochondrial dysfunction as a putative causative mechanism.
Prostaglandins Leukot.
Essent. Fatty Acids
.
81
,
253
–
264 (2009).
14.
H. Wang, S. Liang, M. Wang, J. Gao, C. Sun, J. Wang, W. Xia, S. Wu, S. J. Sumner, F.
Zhang, C. Sun, L. Wu, Potential serum biomarkers from a metabolomics study of autism.
J. Psychiatry Neurosci.
41
, 27
–
37 (2016).
15.
P. R. West, D. G.
Amaral, P. Bais, A. M. Smith, L. A. Egnash, M. E. Ross, J. A. Palmer, B.
R. Fontaine, K. R. Conard, B. A. Corbett, G. G. Cezar, E. L. R. Donley, R. E. Burrier,
Metabolomics as a Tool for Discovery of Biomarkers of Autism Spectrum Disorder in the
Blood Pla
sma of Children.
PLoS ONE
.
9
, e112445 (2014).
16.
R. Cozzolino, L. De Magistris, P. Saggese, M. Stocchero, A. Martignetti, M. Di Stasio, A.
Malorni, R. Marotta, F. Boscaino, L. Malorni, Use of solid
-
phase microextraction coupled
to gas chromatography
–
mass
spectrometry for determination of urinary volatile organic
compounds in autistic children compared with healthy controls.
Anal. Bioanal. Chem.
406
,
4649
–
4662 (2014).
17.
X. Ming, T. P. Stein, V. Barnes, N. Rhodes, L. Guo, Metabolic Perturbance in Autism
Spectrum Disorders: A Metabolomics Study.
J. Proteome Res.
11
, 5856
–
5862 (2012).
18.
A. Noto, V. Fanos, L. Barberini, D. Grapov, C. Fattuoni, M. Zaffanello, A. Casanova, G.
Fenu, A. De Giacomo, M. De Angelis, C. Moretti, P. Papoff, R. Ditonno, R. Francavi
lla,
The urinary metabolomics profile of an Italian autistic children population and their
unaffected siblings.
J. Matern.
-
Fetal Neonatal Med. Off. J. Eur. Assoc. Perinat. Med. Fed.
Asia Ocean. Perinat. Soc. Int. Soc. Perinat. Obstet.
27 Suppl 2
,
46
–
52 (2014).
19.
W. Shaw, E. Kassen, E. Chaves, Increased urinary excretion of analogs of Krebs cycle
metabolites and arabinose in two brothers with autistic features.
Clin. Chem.
41
, 1094
–
1104 (1995).
20.
L. Altieri, C. Neri, R. Sacco, P. Curatolo, A.
Benvenuto, F. Muratori, E. Santocchi, C.
Bravaccio, C. Lenti, M. Saccani, R. Rigardetto, M. Gandione, A. Urbani, A. M. Persico,
Urinary p
-
cresol is elevated in small children with severe autism spectrum disorder.
Biomark. Biochem. Indic. Expo. Response Sus
ceptibility Chem.
16
, 252
–
260 (2011).
21.
P. Emond, S. Mavel, N. Aïdoud, L. Nadal
-
Desbarats, F. Montigny, F. Bonnet
-
Brilhault, C.
Barthélémy, M. Merten, P. Sarda, F. Laumonnier, P. Vourc’h, H. Blasco, C. R. Andres,
GC
-
MS
-
based urine metabolic profiling of
autism spectrum disorders.
Anal. Bioanal.
Chem.
405
, 5291
–
5300 (2013).
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
21
22.
S. Gabriele, R. Sacco, S. Cerullo, C. Neri, A. Urbani, G. Tripi, J. Malvy, C. Barthelemy, F.
Bonnet
-
Brihault, A. M. Persico, Urinary p
-
cresol is elevated in young French children
with autism spectrum disorder: a replication study.
Biomark. Biochem. Indic. Expo.
Response Susceptibility Chem.
19
, 463
–
470 (2014).
23.
A. W. Lis, I. Mclaughlin, R. K. Mpclaughlin, E. W. Lis, E. G. Stubbs, Profiles of ultraviolet
-
absorbing components of
urine from autistic children, as obtained by high
-
resolution ion
-
exchange chromatography.
Clin. Chem.
22
, 1528
–
1532 (1976).
24.
M. Lussu, A. Noto, A. Masili, A. C. Rinaldi, A. Dessì, M. D. Angelis, A. D. Giacomo, V.
Fanos, L. Atzori, R. Francavilla, The u
rinary 1H
-
NMR metabolomics profile of an italian
autistic children population and their unaffected siblings.
Autism Res.
10
, 1058
–
1066
(2017).
25.
L. Nadal
-
Desbarats, N. Aïdoud, P. Emond, H. Blasco, I. Filipiak, P. Sarda, F. Bonnet
-
Brilhault, S. Mavel, C.
R. Andres, Combined
1
H
-
NMR and
1
H
–
13
C HSQC
-
NMR to
improve urinary screening in autism spectrum disorders.
The Analyst
.
139
, 3460
–
3468
(2014).
26.
I. K. S. Yap, M. Angley, K. A. Veselkov
, E. Holmes, J. C. Lindon, J. K. Nicholson, Urinary
metabolic phenotyping differentiates children with autism from their unaffected siblings
and age
-
matched controls.
J. Proteome Res.
9
, 2996
–
3004 (2010).
27.
J. B. Adams, L. J. Johansen, L. D. Powell, D.
Quig, R. A. Rubin, Gastrointestinal flora and
gastrointestinal status in children with autism
--
comparisons to typical children and
correlation with autism severity.
BMC Gastroenterol.
11
, 22 (2011).
28.
M. D. Angelis, M. Piccolo, L. Vannini, S. Siragusa,
A. D. Giacomo, D. I. Serrazzanetti, F.
Cristofori, M. E. Guerzoni, M. Gobbetti, R. Francavilla, Fecal Microbiota and
Metabolome of Children with Autism and Pervasive Developmental Disorder Not
Otherwise Specified.
PLOS ONE
.
8
, e76993 (2013).
29.
D.
-
W. Kan
g, Z. E. Ilhan, N. G. Isern, D. W. Hoyt, D. P. Howsmon, M. Shaffer, C. A.
Lozupone, J. Hahn, J. B. Adams, R. Krajmalnik
-
Brown, Differences in fecal microbial
metabolites and microbiota of children with autism spectrum disorders.
Anaerobe
.
49
,
121
–
131 (2018
).
30.
L. Wang, C. T. Christophersen, M. J. Sorich, J. P. Gerber, M. T. Angley, M. A. Conlon,
Elevated Fecal Short Chain Fatty Acid and Ammonia Concentrations in Children with
Autism Spectrum Disorder.
Dig. Dis. Sci.
57
, 2096
–
2102 (2012).
31.
T. Bitar, S
. Mavel, P. Emond, L. Nadal
-
Desbarats, A. Lefèvre, H. Mattar, M. Soufia, H.
Blasco, P. Vourc’h, W. Hleihel, C. R. Andres, Identification of metabolic pathway
disturbances using multimodal metabolomics in autistic disorders in a Middle Eastern
population.
J
. Pharm. Biomed. Anal.
152
, 57
–
65 (2018).
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
22
32.
F. Gevi, L. Zolla, S. Gabriele, A. M. Persico, Urinary metabolomics of young Italian autistic
children supports abnormal tryptophan and purine metabolism.
Mol. Autism
.
7
(2016),
doi:10.1186/s13229
-
016
-
0109
-
5.
33.
H. Kuwabara, H. Yamasue, S. Koike, H. Inoue, Y. Kawakubo, M. Kuroda, Y. Takano, N.
Iwashiro, T. Natsubori, Y. Aoki, Y. Kano, K. Kasai, Altered metabolites in the plasma of
autism spectrum disorder: a capillary electrophoresis time
-
of
-
flight mass spect
roscopy
study.
PloS One
.
8
, e73814 (2013).
34.
J. S. Orozco, I. Hertz
-
Picciotto, L. Abbeduto, C. M. Slupsky, Metabolomics analysis of
children with autism, idiopathic
-
developmental delays, and Down syndrome.
Transl.
Psychiatry
.
9
, 243 (2019).
35.
B. J. F
erguson, K. Dovgan, D. Severns, S. Martin, S. Marler, K. Gross Margolis, M. L.
Bauman, J. Veenstra
-
VanderWeele, K. Sohl, D. Q. Beversdorf, Lack of Associations
Between Dietary Intake and Gastrointestinal Symptoms in Autism Spectrum Disorder.
Front. Psychia
try
.
10
, 528 (2019).
36.
M. R. Sanctuary, J. N. Kain, K. Angkustsiri, J. B. German, Dietary Considerations in Autism
Spectrum Disorders: The Potential Role of Protein Digestion and Microbial Putrefaction
in the Gut
-
Brain Axis.
Front. Nutr.
5
(2018), doi:1
0.3389/fnut.2018.00040.
37.
S. A. Cermak, C. Curtin, L. G. Bandini, Food selectivity and sensory sensitivity in children
with autism spectrum disorders.
J. Am. Diet. Assoc.
110
, 238
–
246 (2010).
38.
B. D. Needham, W. Tang, W.
-
L. Wu, Searching for the gut
microbial contributing factors to
social behavior in rodent models of autism spectrum disorder.
Dev. Neurobiol.
78
, 474
–
499 (2018).
39.
J. Pulikkan, A. Mazumder
, T. Grace, Role of the Gut Microbiome in Autism Spectrum
Disorders.
Adv. Exp. Med. Biol.
1118
, 253
–
269 (2019).
40.
H. E. Vuong, J. M. Yano, T. C. Fung, E. Y. Hsiao, The Microbiome and Host Behavior.
Annu. Rev. Neurosci.
40
, 21
–
49 (2017).
41.
I. Hertz
-
Pi
cciotto, L. A. Croen, R. Hansen, C. R. Jones, J. van de Water, I. N. Pessah, The
CHARGE Study: An Epidemiologic Investigation of Genetic and Environmental Factors
Contributing to Autism.
Environ. Health Perspect.
114
, 1119
–
1125 (2006).
42.
D. R. Rose, H.
Yang, G. Serena, C. Sturgeon, B. Ma, M. Careaga, H. K. Hughes, K.
Angkustsiri, M. Rose, I. Hertz
-
Picciotto, J. Van de Water, R. L. Hansen, J. Ravel, A.
Fasano, P. Ashwood, Differential immune responses and microbiota profiles in children
with autism spectr
um disorders and co
-
morbid gastrointestinal symptoms.
Brain. Behav.
Immun.
70
, 354
–
368 (2018).
43.
D. A. Drossman, The Functional Gastrointestinal Disorders and the Rome III Process.
Gastroenterology
.
130
, 1377
–
1390 (2006).
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
23
44.
R. S. Kelly, A. Boulin, N.
Laranjo, K. Lee
-
Sarwar, S. H. Chu, A. P. Yadama, V. Carey, A.
A. Litonjua, J. Lasky
-
Su, S. T. Weiss, Metabolomics and Communication Skills
Development in Children; Evidence from the Ages and Stages Questionnaire.
Metabolites
.
9
, 42 (2019).
45.
A. Weiler,
A. Volkenhoff, H. Hertenstein, S. Schirmeier, Metabolite transport across the
mammalian and insect brain diffusion barriers.
Neurobiol. Dis.
107
, 15
–
31 (2017).
46.
R. E. Frye, S. Melnyk, D. F. MacFabe, Unique acyl
-
carnitine profiles are potential
biomark
ers for acquired mitochondrial disease in autism spectrum disorder.
Transl.
Psychiatry
.
3
, e220 (2013).
47.
C. Gillberg, E. Fernell, E. Kočovská, H. Minnis, T. Bourgeron, L. Thompson, C. S. Allely,
The role of cholesterol metabolism and various steroid ab
normalities in autism spectrum
disorders: A hypothesis paper.
Autism Res. Off. J. Int. Soc. Autism Res.
10
, 1022
–
1044
(2017).
48.
A. V. Golubeva, S. A. Joyce, G. Moloney, A. Burokas, E. Sherwin, S. Arboleya, I. Flynn, D.
Khochanskiy, A. Moya
-
Pérez, V. Pet
erson, K. Rea, K. Murphy, O. Makarova, S.
Buravkov, N. P. Hyland, C. Stanton, G. Clarke, C. G. M. Gahan, T. G. Dinan, J. F. Cryan,
Microbiota
-
related Changes in Bile Acid & Tryptophan Metabolism are Associated with
Gastrointestinal Dysfunction in a Mouse M
odel of Autism.
EBioMedicine
.
24
, 166
–
178
(2017).
49.
I. Kurochkin, E. Khrameeva, A. Tkachev, V. Stepanova, A. Vanyushkina, E.
Stekolshchikova, Q. Li, D. Zubkov, P. Shichkova, T. Halene, L. Willmitzer, P. Giavalisco,
S. Akbarian, P. Khaitovich, Metabolome
signature of autism in the human prefrontal
cortex.
Commun. Biol.
2
, 234 (2019).
50.
E. Y. Hsiao, S. W. McBride, S. Hsien, G. Sharon, E. R. Hyde, T. McCue, J. A. Codelli, J.
Chow, S. E. Reisman, J. F. Petrosino, P. H. Patterson, S. K. Mazmanian, Microbio
ta
Modulate Behavioral and Physiological Abnormalities Associated with
Neurodevelopmental Disorders.
Cell
.
155
, 1451
–
1463 (2013).
51.
H. M. Roager, T. R. Licht, Microbial tryptophan catabolites in health and disease.
Nat.
Commun.
9
, 1
–
10 (2018).
52.
H.
-
F
. Zheng, W.
-
Q. Wang, X.
-
M. Li, G. Rauw, G. B. Baker, Body fluid levels of neuroactive
amino acids in autism spectrum disorders: a review of the literature.
Amino Acids
.
49
, 57
–
65 (2017).
53.
B. Diémé, S. Mavel, H. Blasco, G. Tripi, F. Bonnet
-
Brilhault, J.
Malvy, C. Bocca, C. R.
Andres, L. Nadal
-
Desbarats, P. Emond, Metabolomics Study of Urine in Autism Spectrum
Disorders Using a Multiplatform Analytical Methodology.
J. Proteome Res.
14
, 5273
–
5282 (2015).
54.
L. Santamaría, I. Reverón, F. L. de Felipe, B.
de las Rivas, R. Muñoz, Ethylphenol
Formation by Lactobacillus plantarum: Identification of the Enzyme Involved in the
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
24
Reduction of Vinylphenols.
Appl. Environ. Microbiol.
84
(2018),
doi:10.1128/AEM.01064
-
18.
55.
B. Pragnya, J. S. L. Kameshwari, B.
Veeresh, Ameliorating effect of piperine on behavioral
abnormalities and oxidative markers in sodium valproate induced autism in BALB/C
mice.
Behav. Brain Res.
270
, 86
–
94 (2014).
56.
B. Auyeung, J. Ahluwalia, L. Thomson, K. Taylor, G. Hackett, K. J. O’Don
nell, S. Baron
-
Cohen, Prenatal versus postnatal sex steroid hormone effects on autistic traits in children
at 18 to 24 months of age.
Mol. Autism
.
3
, 17 (2012).
57.
S. Baron
-
Cohen, B. Auyeung, B. Nørgaard
-
Pedersen, D. M. Hougaard, M. W. Abdallah, L.
Melga
ard, A. S. Cohen, B. Chakrabarti, L. Ruta, M. V. Lombardo, Elevated fetal
steroidogenic activity in autism.
Mol. Psychiatry
.
20
, 369
–
376 (2015).
58.
M. H. Hassan, T. Desoky, H. M. Sakhr, R. H. Gabra, A. H. Bakri, Possible Metabolic
Alterations among Autis
tic Male Children: Clinical and Biochemical Approaches.
J. Mol.
Neurosci. MN
(2019), doi:10.1007/s12031
-
018
-
1225
-
9.
59.
M. P. Mamidala, A. Polinedi, P. T. V. P. Kumar, N. Rajesh, O. R. Vallamkonda, V. Udani,
N. Singhal, V. Rajesh, Maternal hormonal interv
entions as a risk factor for Autism
Spectrum Disorder: an epidemiological assessment from India.
J. Biosci.
38
, 887
–
892
(2013).
60.
P. M. Whitaker
-
Azmitia, M. Lobel, A. Moyer, Low maternal progesterone may contribute to
both obstetrical complications and
autism.
Med. Hypotheses
.
82
, 313
–
318 (2014).
61.
M. D. Majewska, M. Hill, E. Urbanowicz, P. Rok
-
Bujko, P. Bieńkowski, I. Namysłowska, P.
Mierzejewski, Marked elevation of adrenal steroids, especially androgens, in saliva of
prepubertal autistic children.
Eur. Child Adolesc. Psychiatry
.
23
, 485
–
498 (2014).
62.
J. Willing, C. K. Wagner, Exposure to the Synthetic Progestin, 17α
-
Hydroxyprogesterone
Caproate During Development Impairs Cognitive Flexibility in Adulthood.
Endocrinology
.
157
, 77
–
82 (2016).
63.
A. Quartier, L. Chatrousse, C. Redin, C. Keime, N. Haumess
er, A. Maglott
-
Roth, L. Brino,
S. Le Gras, A. Benchoua, J.
-
L. Mandel, A. Piton, Genes and Pathways Regulated by
Androgens in Human Neural Cells, Potential Candidates for the Male Excess in Autism
Spectrum Disorder.
Biol. Psychiatry
.
84
, 239
–
252 (2018).
64.
D. A. Geier, J. K. Kern, P. G. King, L. K. Sykes, M. R. Geier, An evaluation of the role and
treatment of elevated male hormones in autism spectrum disorders.
Acta Neurobiol. Exp.
(Warsz.)
.
72
, 1
–
17 (2012).
65.
S. Bent, B. Lawton, T. Warren, F. Widjaja,
K. Dang, J. W. Fahey, B. Cornblatt, J. M.
Kinchen, K. Delucchi, R. L. Hendren, Identification of urinary metabolites that correlate
with clinical improvements in children with autism treated with sulforaphane from
broccoli.
Mol. Autism
.
9
(2018), doi:10.1
186/s13229
-
018
-
0218
-
4.
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
25
66.
L. Kanner, Autistic disturbances of affective contact.
Acta Paedopsychiatr.
35
, 100
–
136
(1968).
67.
K. K. Griffiths, R. J. Levy, Evidence of Mitochondrial Dysfunction in Autism: Biochemical
Links, Genetic
-
Based Associations, an
d Non
-
Energy
-
Related Mechanisms.
Oxid. Med.
Cell. Longev.
2017
(2017), doi:10.1155/2017/4314025.
68.
A. Legido, R. Jethva, M. J. Goldenthal, Mitochondrial dysfunction in autism.
Semin. Pediatr.
Neurol.
20
, 163
–
175 (2013).
69.
J. Lombard, Autism: a mitoch
ondrial disorder?
Med. Hypotheses
.
50
, 497
–
500 (1998).
70.
M. Coleman, J. P. Blass, Autism and lactic acidosis.
J. Autism Dev. Disord.
15
, 1
–
8 (1985).
71.
S. Mavel, L. Nadal
-
Desbarats, H. Blasco, F. Bonnet
-
Brilhault, C. Barthélémy, F. Montigny,
P. Sarda,
F. Laumonnier, P. Vourc′h, C. R. Andres, P. Emond, 1H
–
13C NMR
-
based urine
metabolic profiling in autism spectrum disorders.
Talanta
.
114
, 95
–
102 (2013).
72.
D. F. MacFabe, Enteric short
-
chain fatty acids: microbial messengers of metabolism,
mitochondria,
and mind: implications in autism spectrum disorders.
Microb. Ecol. Health
Dis.
26
(2015), doi:10.3402/mehd.v26.28177.
73.
C.
-
Y. Chang, D.
-
S. Ke, J.
-
Y. Chen, Essential fatty acids and human brain.
Acta Neurol.
Taiwanica
.
18
, 231
–
241 (2009).
74.
C. W. Kuz
awa, H. T. Chugani, L. I. Grossman, L. Lipovich, O. Muzik, P. R. Hof, D. E.
Wildman, C. C. Sherwood, W. R. Leonard, N. Lange, Metabolic costs and evolutionary
implications of human brain development.
Proc. Natl. Acad. Sci. U. S. A.
111
, 13010
–
13015 (2014).
75.
J. S. O’Brien, E. L. Sampson, Lipid composition of the normal human brain: gray matter,
white matter, and myelin.
J. Lipid Res.
6
, 537
–
544 (1965).
76.
A. M. Smith, J. J. King, P. R. West, M. A. Ludwig, E. L. R. Donley, R. E. Burrier, D. G.
Amaral, A
mino Acid Dysregulation Metabotypes: Potential Biomarkers for Diagnosis and
Individualized Treatment for Subtypes of Autism Spectrum Disorder.
Biol. Psychiatry
.
85
,
345
–
354 (2019).
77.
H. Ripps, W. Shen, Review: Taurine: A “very essential” amino acid.
Mol
. Vis.
18
, 2673
–
2686 (2012).
78.
A. Dickinson, M. Jones, E. Milne, Measuring neural excitation and inhibition in autism:
Different approaches, different findings and different interpretations.
Brain Res.
1648
,
277
–
289 (2016).
79.
C. A. Olson, H. E. Vuong
, J. M. Yano, Q. Y. Liang, D. J. Nusbaum, E. Y. Hsiao, The Gut
Microbiota Mediates the Anti
-
Seizure Effects of the Ketogenic Diet.
Cell
.
173
, 1728
-
1741.e13 (2018).
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
26
80.
A. Liu, W. Zhou, L. Qu, F. He, H. Wang, Y. Wang, C. Cai, X. Li, W. Zhou, M. Wang,
Alter
ed Urinary Amino Acids in Children With Autism Spectrum Disorders.
Front. Cell.
Neurosci.
13
(2019), doi:10.3389/fncel.2019.00007.
81.
A. Duda
-
Chodak, T. Tarko, P. Satora, P. Sroka, Interaction of dietary compounds, especially
polyphenols, with the intest
inal microbiota: a review.
Eur. J. Nutr.
54
, 325
–
341 (2015).
82.
L. Valdés, A. Cuervo, N. Salazar, P. Ruas
-
Madiedo, M. Gueimonde, S. González, The
relationship between phenolic compounds from diet and microbiota: impact on human
health.
Food Funct.
6
, 242
4
–
2439 (2015).
83.
J. Kałużna
-
Czaplińska, E. Żurawicz, W. Struck, M. Markuszewski, Identification of organic
acids as potential biomarkers in the urine of autistic children using gas
chromatography/mass spectrometry.
J. Chromatogr. B Analyt. Technol. Biom
ed. Life. Sci.
966
, 70
–
76 (2014).
84.
W. J. Johannsen, S. H. Friedman, E. I. Feldman, A. Negrete, A Re
-
examination of the
Hippuric Acid
--
Anxiety Relationship.
Psychosom. Med.
24
, 569 (1962).
85.
H. B. Patisaul
, W. Jefferson, The pros and cons of phytoestrogens.
Front. Neuroendocrinol.
31
, 400
–
419 (2010).
86.
J. C. Naviaux, M. A. Schuchbauer, K. Li, L. Wang, V. B. Risbrough, S. B. Powell, R. K.
Naviaux, Reversal of autism
-
like behaviors and metabolism in adult
mice with single
-
dose
antipurinergic therapy.
Transl. Psychiatry
.
4
, e400
–
e400 (2014).
87.
M. Grochowska, T. Laskus, M. Radkowski, Gut Microbiota in Neurological Disorders.
Arch.
Immunol. Ther. Exp. (Warsz.)
.
67
, 375
–
383 (2019).
88.
G. Sharon, T. R. Samp
son, D. H. Geschwind, S. K. Mazmanian, The Central Nervous
System and the Gut Microbiome.
Cell
.
167
, 915
–
932 (2016).
89.
L. M. Iakoucheva, A. R. Muotri, J. Sebat, Getting to the Cores of Autism.
Cell
.
178
, 1287
–
1298 (2019).
90.
L. Desbonnet, G. Clarke, A
. Traplin, O. O’Sullivan, F. Crispie, R. D. Moloney, P. D. Cotter,
T. G. Dinan, J. F. Cryan, Gut microbiota depletion from early adolescence in mice:
Implications for brain and behaviour.
Brain. Behav. Immun.
48
, 165
–
173 (2015).
91.
K. M. Neufeld, N. Kang
, J. Bienenstock, J. A. Foster, Reduced anxiety
-
like behavior and
central neurochemical change in germ
-
free mice.
Neurogastroenterol. Motil. Off. J. Eur.
Gastrointest. Motil. Soc.
23
, 255
–
264, e119 (2011).
92.
M. Kujawska, J. Jodynis
-
Liebert, Polyphenols
in Parkinson’s Disease: A Systematic Review
of In Vivo Studies.
Nutrients
.
10
, 642 (2018).
93.
T. R. Sampson, J. W. Debelius, T. Thron, S. Janssen, G. G. Shastri, Z. E. Ilhan, C. Challis, C.
E. Schretter, S. Rocha, V. Gradinaru, M.
-
F. Chesselet, A. Keshav
arzian, K. M. Shannon,
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
27
R. Krajmalnik
-
Brown, P. Wittung
-
Stafshede, R. Knight, S. K. Mazmanian, Gut
Microbiota Regulate Motor Deficits and Neuroinflammation in a Model of Parkinson’s
Disease.
Cell
.
167
, 1469
-
1480.e12 (2016).
94.
M. Jaglin, M. Rhimi, C. Phil
ippe, N. Pons, A. Bruneau, B. Goustard, V. Daugé, E. Maguin,
L. Naudon, S. Rabot, Indole, a Signaling Molecule Produced by the Gut Microbiota,
Negatively Impacts Emotional Behaviors in Rats.
Front. Neurosci.
12
(2018),
doi:10.3389/fnins.2018.00216.
95.
J.
C. O’Connor, M. A. Lawson, C. André, M. Moreau, J. Lestage, N. Castanon, K. W.
Kelley, R. Dantzer, Lipopolysaccharide
-
induced depressive
-
like behavior is mediated by
indoleamine 2,3
-
dioxygenase activation in mice.
Mol. Psychiatry
.
14
, 511
–
522 (2009).
96.
P. Tian, G. Wang, J. Zhao, H. Zhang, W. Chen, Bifidobacterium with the role of 5
-
hydroxytryptophan synthesis regulation alleviates the symptom of depression and related
microbiota dysbiosis.
J. Nutr. Biochem.
66
, 43
–
51 (2019).
97.
D. Wang, L. Ho, J. Fait
h, K. Ono, E. M. Janle, P. J. Lachcik, B. R. Cooper, A. H. Jannasch,
B. R. D’Arcy, B. A. Williams, M. G. Ferruzzi, S. Levine, W. Zhao, L. Dubner, G. M.
Pasinetti, Role of intestinal microbiota in the generation of polyphenol
-
derived phenolic
acid mediated
attenuation of Alzheimer’s disease β
-
amyloid oligomerization.
Mol. Nutr.
Food Res.
59
, 1025
–
1040 (2015).
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint
28
FIGURE TITLES AND LEGENDS
Figure 1. Plasma and Fecal Metabolomes Differ between ASD and TD and Global
Metabolite Levels Correlate with Clinical Behavioral Scores
(A
-
B) The number of significantly elevated and decreased metabolites (pval<0.05) in ASD
samples compared to the TD contr
ol group by ANOVA contrasts in plasma and feces,
respectively. Samples are stratified by all samples or samples without or with GI symptoms (
-
GI,
+GI).
See also Table S1 and S2.
(
C
-
D
) Pathway analysis results of human plasma and fecal
comparisons (all samp
les), indicating which metabolomic pathways are the most altered between
groups, with enrichment value plotted and p
-
value to the right of each bar. Metabolites within
each pathway could be observed at either higher or lower levels, as this plot only indic
ates
significant changes. (
E
-
F
) Top 30 most distinguishing metabolites between each group in plasma
and feces by random forest analysis, with mean decrease accuracy along the x
-
axis. Metabolites
known to be produced by (asterisks) or influenced by (triangl
es) the gut microbiota are denoted.
The super pathway to which each metabolite belongs to is indicated by color of sphere and
defined in the legend.
See also Figure S1.
(
G
-
H
) Scaled intensity values indicating relative levels
of the most distinguishing mol
ecules between ASD and TD (all samples) in plasma. Asterisks
indicate significance in ANOVA contrasts performed on total metabolomics dataset. p values:
2.87E
-
09(
G
) and 1.57E
-
06(
H
). Data are represented as mean ± SEM. (
I
-
J
) Scaled intensity
values
indicating relative levels of the most distinguishing molecules between ASD and TD (all
samples) in feces. Asterisks indicate significance in ANOVA contrasts performed on total
metabolomics dataset. p values: 1.37E
-
05(
I
) and 0.002(
J
). Data are represented
as mean ± SEM.
(K) Top correlated plasma metabolites that covary with margaroylcarnitine and indolelactate. All
p
-
values <0.0001. (L) Top correlated fecal metabolites that covary with nicotinamide and 9
-
HOTrE. All p
-
values <0.0001.
(M) Spearman correlation
s between behavior scores of ASD
children in the ADIR diagnostic test (Verbal, Social, and Nonverbal metrics) and the ADOS
-
SS
and metabolite pathways in ASD samples. Directionality of correlation is indicated in the legend
at bottom. Colors of pathways are
defined at the top left of the chart. A split box means that both
positive and negative correlations occur with metabolites within that pathway.
See also Tables S1
and S2.
(N) PCA plot comparing the metabolic profile least and most severe ~quartiles withi
n
verbal, social and ADOS
-
SS scores. PCA input included all metabolites significantly associated
with the behavior. Clustering is denoted by ellipses of the 95% confidence interval.
LPC,
lysophosphatidylcholine; CE, cholesterol ester; FFA, free fatty acid;
androst., androstane;
hydroxypreg, hydroxypregnenalone; PFOS, perfluorooctanesulfonic acid; hydroxy
-
CMPF,
hydroxy
-
3
-
carboxy
-
4
-
methyl
-
5
-
propyl
-
2
-
furanpropionate; DHEA
-
s, dehydroepiandrosterone
sulfate; 9
-
HOTrE, 9S
-
hydroxy
-
10E,12Z,15Z
-
octadecatrienoic acid
; AMP, adenosine
monophosphate; HExCer, Hexosylceramide; PC, phosphatidylcholine; DAG, diacylglycerol;
TAG, triacylglycerol; MFA, monounsaturated fatty acid; PE plas, phosphatidylethanolamine
plasmalogens; Endocann, endocannabinoid; Met, methionine; Cys, c
ysteine; SAM, s
-
adenosyl
methionine; Ala, alanine; Asp, aspartate;
ɣ
-
glutAA, gamma
-
glutamyl amino acids; Second. Bile,
secondary bile; MAG, monoacylglycerol; Gly, glycine; Ser, serine; thr, threonine.
See also
Figures S1 and S2.
Figure 2.
Lipid and Xenobi
otic Metabolites Correlate with GI Symptoms
was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission.
The copyright holder for this preprint (which
this version posted May 19, 2020.
.
https://doi.org/10.1101/2020.05.17.098806
doi:
bioRxiv preprint