Literature DB >> 33313377

Evaluation of a visual interpretation method for tau-PET with 18F-flortaucipir.

Ida Sonni1,2, Orit H Lesman Segev3,4, Suzanne L Baker1, Leonardo Iaccarino3, Deniz Korman5, Gil D Rabinovici1,3,6, William J Jagust1,5, Susan M Landau5, Renaud La Joie3.   

Abstract

INTRODUCTION: Positron emission tomography targeting tau (tau-PET) is a promising diagnostic tool for the identification of Alzheimer's disease (AD). Currently available data rely on quantitative measures, and a visual interpretation method, critical for clinical translation, is needed.
METHODS: We developed a visual interpretation method for 18F-flortaucipir tau-PET and tested it on 274 individuals (cognitively normal controls, patients with mild cognitive impairment [MCI], AD dementia, and non-AD diagnoses). Two readers interpreted 18F-flortaucipir PET using two complementary indices: a global visual score and a visual distribution pattern.
RESULTS: Global visual scores were reliable, correlated with global cortical 18F-flortaucipir standardized uptake value ratio (SUVR) and were associated with clinical diagnosis and amyloid status. The AD-like 18F-flortaucipir pattern had good sensitivity and specificity to identify amyloid-positive patients with AD dementia or MCI. DISCUSSION: This 18F-flortaucipir visual rating scheme is associated with SUVR quantification, clinical diagnosis, and amyloid status, and constitutes a promising approach to tau measurement in clinical settings.
© 2020 The Authors. Alzheimer's & Dementia: Diagnosis, Assessment & Disease Monitoring published by Wiley Periodicals, LLC on behalf of Alzheimer's Association.

Entities:  

Keywords:  Alzheimer's disease; flortaucipir; qualitative assessment; tau positron emission tomography; visual assessment

Year:  2020        PMID: 33313377      PMCID: PMC7699207          DOI: 10.1002/dad2.12133

Source DB:  PubMed          Journal:  Alzheimers Dement (Amst)        ISSN: 2352-8729


BACKGROUND

Positron emission tomography (PET) has transformed the way we study cognitive decline and dementia, offering the possibility to visualize age‐associated pathophysiological changes, with major consequences for diagnosis, disease monitoring, and screening strategies for clinical trials. The importance of biomarkers of amyloid beta (Aβ) deposition has been recognized by the National Institute on Aging and the Alzheimer Association (NIA‐AA) in 2011. The advancement of Aβ‐PET with C‐Pittsburgh compound B (PiB), and subsequent fluorinated counterparts , , was rapid. Aβ‐PET radiotracers are now clinically available and interpreted visually using a binary (positive/negative) read. , The successful experience of Aβ‐PET has prompted the development of PET radiopharmaceuticals targeting the other aggregated protein implicated in Alzheimer's disease (AD) and other neurodegenerative diseases: tau. Multiple tau‐targeted radiopharmaceuticals are being evaluated in humans. The most widely used tracer, 18F‐flortaucipir (FTP; previously 18F‐T807 and 18F‐AV1451), , , binds with high affinity to paired helical filament tau. FTP‐SUVR (standardized uptake value ratio) values correlate with cognitive deficits, phosphorylated tau in CSF , and plasma, and post mortem tau burden in neuropathologically proven AD cases. , FTP signal strongly correlates with neurodegeneration , , and discriminates AD patients from healthy controls and patients with other neurodegenerative diseases with high accuracy. , Altogether, FTP‐PET could be a promising diagnostic tool, although most published work relies on SUVR quantification and a visual interpretation method, critical for clinical translation, is needed. Several features of FTP‐PET complicate the definition of a simple binary read similar to Aβ‐PET: strong between‐patient variations in signal intensity and regional distribution, mild‐to‐moderate signal in non‐AD associated conditions, and the presence of “off‐target” binding. Given these challenges, our goals were: (1) to develop a reliable approach to guide visual interpretation and reporting of FTP‐PET scans, (2) to evaluate the ability of FTP‐PET visual reads to distinguish Aβ‐positive patients with AD dementia (ADdem) and mild cognitive impairment (MCI) from other participants, (3) to compare visual reads to a previously developed definition of tau‐positivity using temporal SUVRs.

MATERIALS AND METHODS

Study design

Development phase

A pilot phase was first led by a nuclear medicine physician (IS) and a radiologist (OLS) to develop a visual rating scheme, using a set of 60 FTP‐PET scans. This phase allowed the readers to (1) gain familiarity with the different FTP patterns across clinical diagnoses, (2) choose a color scale to analyze images, and (3) develop and test a visual read procedure before proceeding to a larger scale analysis. Subjects were selected from the Berkeley Aging Cohort Study/University of California San Francisco (BACS/UCSF) cohort and included 20 cognitively normal controls, 20 patients with a clinical diagnosis of MCI or probable AD dementia, and 20 patients with non‐AD neurodegenerative syndromes. This phase was conducted using FTP SUVR images, intensity normalized by mean inferior cerebellar cortex signal.

RESEARCH IN CONTEXT

Systematic review: Brain imaging has become increasingly important in characterizing patients with cognitive decline in the context of research studies, clinical diagnosis, and clinical trials. Positron emission tomography (PET) radiotracers targeting tau pathology were recently developed and are currently limited to investigational settings. After reviewing the growing literature on flortaucipir‐PET, we identified tau‐PET as a promising diagnostic tool although a visual interpretation method, critical for clinical translation, is needed. Interpretation: Our proposed flortaucipir‐PET visual read approach provides ratings that are strongly associated with standardized uptake value ratio quantification, clinical diagnosis, and amyloid beta status across two cohorts with distinct characteristics. Inter‐ and intra‐rater reliability for our main indices are strong, suggesting that our visual approach is reproducible and can represent a promising alternative to quantitative tau‐PET measurements in clinical settings. Future directions: Future investigations including an autopsy‐based gold standard are needed to evaluate the exact potential and limitations of tau‐PET and the proposed visual approach.

Testing phase

The study was conducted on a separate set of 274 non‐intensity‐normalized, unitless FTP‐PET images to make the visual interpretation independent from image processing, which is not often used in clinical settings. Two co‐authors (RLJ and SML) selected 274 subjects from two independent cohorts (137 from BACS/UCSF, 137 from Alzheimer's Disease Neuroimaging Initiative [ADNI]), including various clinical syndromes, Aβ‐positive and negative cases. Two PET readers (IS and OLS) were unaware of the sample composition. Scans were anonymized, randomized, and independently rated by both readers blinded to clinical diagnosis and any other clinical, radiological, or demographic information using the visual rating method described below. Ten percent of the scans (15 from BACS/UCSF, 15 from ADNI) were randomly selected, assigned a second anonymized number, and read twice to evaluate intra‐reader reliability.

Study participants

The institutional review boards of all participating institutions approved the study and informed consent was obtained from all participants or authorized representatives. Cohort characteristics are summarized in Table 1.
TABLE 1

Cohort characteristics

BACS/UCSFADNI
YCOCMCIADdemnon‐AD a OCMCIADdem
N115183631604730
Age39.2 ± 16.076.8 ± 6.171.0 ± 6.062.6 ± 8.066.6 ± 7.675.5 ± 7.076.8 ± 7.678.5 ± 9.7
Age (range)21–5960–9363–8048–8246–7959–9460–9256‐94
Females (%)10 (91%)24 (46%)3 (38%)23 (64%)15 (48%)29 (48%)31 (66%)14 (47%)
Years of education16.8 ± 2.016.7 ± 1.918.3 ± 3.216.9 ± 3.016.4 ± 3.116.5 ± 2.316.6 ± 2.915.1 ± 2.6
MMSE29.0 ± 1.528.8 ± 1.227.9 ± 2.221.3 ± 5.024.1 ± 5.729.1 ± 1.128.0 ± 2.021.7 ± 4.9
Aβ‐positive (%) b 0 (0%) b 23 (45%)8 (100%)36 (100%)7 (23%) b 16 (27%)20 (43%)27 (90%)
APOE ε4 carriers (%) c 4 (44%) b 15 (29%)5 (63%)21 (58%)8 (26%)19 (32%)14 (30%)13 (43%)

Abbreviations: Aβ, amyloid beta; ADdem, patients with a clinical diagnosis of Alzheimer's disease dementia; ADNI, Alzheimer's Disease Neuroimaging Initiative; APOE, apolipoprotein E; BACS, Berkeley Aging Cohort Study; MCI, patients with a clinical diagnosis of mild cognitive impairment; MMSE, Mini‐Mental State Examination; non‐AD, patients with a clinical diagnosis of non‐AD neurodegenerative syndrome; OC, older controls; PET, positron emission tomography; PiB, Pittsburgh compound B; UCSF, University of California San Francisco; YC, young controls.

Notes: Clinical diagnoses were independent from Aβ status. For continuous variables, mean ± SD is indicated unless specified otherwise.

Non AD cases: eight with behavioral variant frontotemporal dementia, seven with non‐fluent primary progressive aphasia, six with corticobasal syndrome, four with Parkinson's disease, three with progressive supranuclear palsy, two with semantic variant primary progressive aphasia.

Missing data in the BACS/UCSF data: Aβ status for 7 YC and 1 non‐AD (patient with behavioral variant frontotemporal dementia who carried a V337M MAPT mutation known to cause Alzheimer's like 3R/4R paired helical filaments tau), and APOE genotype for 2 YC. The Aβ tracer used in BACS/UCSF was PiB, while in ADNI florbetapir was used for 128 subjects and florbetaben for 9. PiB‐PET was missing for a UCSF patient with clinical AD dementia who died 36 months after FTP‐PET and showed high Aβ burden at autopsy (Thal 5 and CERAD frequent); this patient was considered Aβ‐positive in all tables, figures, and analyses.

APOE ε4 indicates the percentage of patients with APOE ε4/ε2, APOE ε4/ε3 or APOE ε4/ε4.

Cohort characteristics Abbreviations: Aβ, amyloid beta; ADdem, patients with a clinical diagnosis of Alzheimer's disease dementia; ADNI, Alzheimer's Disease Neuroimaging Initiative; APOE, apolipoprotein E; BACS, Berkeley Aging Cohort Study; MCI, patients with a clinical diagnosis of mild cognitive impairment; MMSE, Mini‐Mental State Examination; non‐AD, patients with a clinical diagnosis of non‐AD neurodegenerative syndrome; OC, older controls; PET, positron emission tomography; PiB, Pittsburgh compound B; UCSF, University of California San Francisco; YC, young controls. Notes: Clinical diagnoses were independent from Aβ status. For continuous variables, mean ± SD is indicated unless specified otherwise. Non AD cases: eight with behavioral variant frontotemporal dementia, seven with non‐fluent primary progressive aphasia, six with corticobasal syndrome, four with Parkinson's disease, three with progressive supranuclear palsy, two with semantic variant primary progressive aphasia. Missing data in the BACS/UCSF data: Aβ status for 7 YC and 1 non‐AD (patient with behavioral variant frontotemporal dementia who carried a V337M MAPT mutation known to cause Alzheimer's like 3R/4R paired helical filaments tau), and APOE genotype for 2 YC. The Aβ tracer used in BACS/UCSF was PiB, while in ADNI florbetapir was used for 128 subjects and florbetaben for 9. PiB‐PET was missing for a UCSF patient with clinical AD dementia who died 36 months after FTP‐PET and showed high Aβ burden at autopsy (Thal 5 and CERAD frequent); this patient was considered Aβ‐positive in all tables, figures, and analyses. APOE ε4 indicates the percentage of patients with APOE ε4/ε2, APOE ε4/ε3 or APOE ε4/ε4.

Sample 1—BACS/UCSF

Participants were selected from BACS, and the UCSF Memory and Aging Center research cohorts. The sample included cognitively normal individuals, patients with clinical diagnosis of MCI, AD dementia, or non‐Alzheimer's neurodegenerative syndromes. Cognitively normal participants included 11 young controls (YC, aged 20 to 60 years) and 55 older controls (OC, age>60).

Sample 2—ADNI

ADNI is a multi‐site longitudinal study. Subjects selected from this cohort included OC and patients with amnestic MCI and probable AD dementia.

Magnetic resonance imaging

T1‐weighted magnetic resonance imaging (MRI) was used to define the regions of interest (ROI) used for FTP quantification and were not available to the readers. All scans were processed with FreeSurfer 5.3 (http://surfer.nmr.mgh.harvard.edu/) to derive ROIs in each subject's native space (see below).

PET acquisition and processing

Details of the PET acquisition are provided in the supporting information. Briefly, the FTP‐PET data were acquired as four 5‐minute frames between 80 and 100 minutes after the injection of ≈370 MBq of radiotracer. Frames were realigned and averaged; resulting images were used for the visual read to make the visual interpretation entirely qualitative. For quantification purposes, SUVR images were created using the inferior cerebellar cortex as the reference region. FreeSurfer segmentation was used to extract the average SUVR value in two ROIs: (1) the entire cortex and (2) a temporal meta‐ROI (including bilateral amygdala, entorhinal, parahippocampal, fusiform, and inferior and middle temporal cortices). We used the latter as a quantification‐based criteria to define tau‐positivity, using the 1.27 threshold derived from independent samples in Ossenkoppele et al. Aβ‐PET was obtained using PiB‐PET (135 BACS/UCSF participants), 18F‐florbetapir (n = 128, ADNI), or 18F‐florbetaben (n = 9, ADNI). Aβ‐PET were processed as previously described and Aβ status was defined from tracer/pipeline specific thresholds: 1.065 DVR for PiB, 1.11 for 18F‐florbetapir, and 1.08 for 18F‐florbetaben). , , , , ,

FTP‐PET visual interpretation method: visual score and distribution pattern

Readers assessed FTP‐PET scans with no access to anatomical imaging, using the image viewer software Mango (version 4.0.1–Research Imaging Institute, UTHSCSA; http://ric.uthscsa.edu/mango) and the “Spectrum” color scale. The color scale was manually adjusted so that the predominant color in the inferior cerebellar cortex would be set as the midpoint of the color scale (Figure 1A). Areas of frequent off‐target FTP binding (Figure 1A) were taken into consideration. The visual interpretation method provides a report for each FTP scan describing two parameters detailed below.
FIGURE 1

Overview of the visual read approach and global score. A, The image intensity has to be manually adjusted by each reader, fixing the value of the inferior cerebellar cortex value to the mid‐range of the color scale (green area). For both visual indices (global visual score and flortaucipir pattern), readers considered tracer binding beyond the common areas of off‐target (non‐specific) binding, as illustrated here in amyloid‐negative controls. B, The global visual score is based on seven regions of interest (ROIs), only shown on the left to simplify the display, each scored on a 0 to 2 scale (white arrows illustrate the scale applied to the inferior temporal lobe), resulting in a 0 to 14 global scale

Overview of the visual read approach and global score. A, The image intensity has to be manually adjusted by each reader, fixing the value of the inferior cerebellar cortex value to the mid‐range of the color scale (green area). For both visual indices (global visual score and flortaucipir pattern), readers considered tracer binding beyond the common areas of off‐target (non‐specific) binding, as illustrated here in amyloid‐negative controls. B, The global visual score is based on seven regions of interest (ROIs), only shown on the left to simplify the display, each scored on a 0 to 2 scale (white arrows illustrate the scale applied to the inferior temporal lobe), resulting in a 0 to 14 global scale

Global visual score

FTP‐PET images were assessed following a predefined regional cortical binding system (0 = no binding; 1 = mild binding; 2 = intense binding) for seven brain regions (Figure 1B). In the case of asymmetrical uptake within the ROI, that is, higher on one side than on the contralateral, the highest signal intensity was chosen. The global visual score was obtained by summing each score given to the seven cortical ROIs to provide a measure of whole cortex FTP intensity (ranging from 0, no uptake in any of the ROIs, to 14, intense uptake in all ROIs).

Distribution pattern

Four FTP distribution patterns were described (Figure 2):
FIGURE 2

Four distribution patterns. Three cases are provided to illustrate each pattern; for all 12 examples shown here, both readers independently assigned the same flortaucipir visual pattern. For each case, four slices are provided, including an axial slice at the level of the inferior cerebellum, and axial, coronal, and sagittal sections throughout the brain

Pattern I (negative scan): absence of FTP signal in any brain area beyond background or off‐target binding. Pattern II (mild temporal binding only): mild to moderate elevation of FTP signal limited to the medial temporal cortex and fusiform gyrus. This pattern could be consistent with early (ie, I to III) Braak stage tau pathology that can be seen in older individuals with or without cognitive deficits. , , , Pattern III (AD‐like binding): FTP distribution consistent with the neuropathological distribution of tau in Braak stage ≥ IV, , , , extending beyond the medial temporal/fusiform area (eg, with lateral temporal, parietal, or frontal cortex binding). Pattern IV (non‐AD‐like): atypical FTP signal distribution, not following the known distribution of tau tangles in AD by Braak et al. , , , This pattern includes, but is not limited to, predominant white matter or subcortical binding as previously described in non‐AD syndromes. Four distribution patterns. Three cases are provided to illustrate each pattern; for all 12 examples shown here, both readers independently assigned the same flortaucipir visual pattern. For each case, four slices are provided, including an axial slice at the level of the inferior cerebellum, and axial, coronal, and sagittal sections throughout the brain

Statistical analysis

Correlations were quantified using Spearman's rho coefficient. The effects of Aβ status and clinical diagnostic on the global visual score was assessed in a general linear model including a group x Aβ interaction term. Relationships between FTP distribution patterns and temporal meta‐ROI SUVR or global visual scores were tested using Kruskall Wallis analysis of variance, followed by post hoc pairwise Dwass‐Steel‐Critchlow‐Fligner tests. Linear weighted kappa (κ) coefficients were used to evaluate inter‐rater and intra‐rater agreement on the global visual scores; non‐weighted κ to evaluate inter‐ and intra‐rater reliability for the identification of the four distribution patterns. When analyzing mismatches between distribution patterns, we distinguished major disagreements (between “negative” and “AD‐like” or “AD‐like” and “non‐AD‐like” patterns) from minor disagreements (other combinations). Additional reliability analyses were conducted to specifically focus on reliability of the “AD‐like” pattern, versus any of the other patterns.

RESULTS

Visual scores: reliability and association with cortical SUVR

Global visual scores from the two readers were strongly correlated (rho = .805, P <.001) and inter‐rater agreement was good (linear weighted κ = .73, 95% confidence interval [CI; .68 to .78]; Figure S1A in supporting information). Intra‐rater agreement was high: linear weighted κreader1 = 0.79 [.68 to .91] and κreader2 = 0.84 [.7 to .93] (Figure S1B). The average visual score was significantly associated with cortical FTP–SUVR (rho = .71, P <.001; Figure 3A). Inspection of the scatter plot showed a non‐linear relationship and differences in dynamic range between the two variables: SUVR values appear little related to the visual score in the low range, while the visual score reaches a ceiling for cortical SUVR >1.6. Complementary analyses were run in three non‐overlapping subgroups of similar visual score span, showing that the correlation between SUVR and visual score was significant in low values (rho = .29, P <.001 when restricted to scans with a score ≤4), but gradually strengthened with increasing visual scores (rho = .60, for scores between 5 and 9, inclusive; rho = .78 for scores ≥10).
FIGURE 3

Global visual score: association with global cortical standardized uptake value ratio (SUVR), clinical diagnosis, and amyloid status. A, Associations between flortaucipir cortical SUVR and global visual score (average of both readers). See Figure S2 in supporting information for similar plots using each reader's score or each cohort separately. The scatter plot shows a non‐linear relationship and differences in dynamic range between the two variables: SUVR values appear little related to the visual score in the low range, while the visual score reaches a ceiling for cortical SUVR >1.6. B, Association of clinical diagnosis and amyloid status on the global visual score (average of both readers). See Figure S3 in supporting information for similar plots using each reader's score. The YC group was excluded from the general linear model results described on top (overall R2 = .68). YC, young controls; OC, older controls; MCI, patients with a clinical diagnosis of mild cognitive impairment; ADdem, patients with a clinical diagnosis of Alzheimer's disease dementia; non‐AD, patients with a clinical diagnosis of non‐AD neurodegenerative syndrome. Clinical diagnoses were independent from amyloid beta (Aβ) status. C, Receiver operating curve (ROC) testing on the ability of the average global visual score and two SUVR measures to distinguish Aβ‐positive MCI/ADdem (n = 90) versus all other participants (n = 184). AUROC, area under the receiver operating curve

Global visual score: association with global cortical standardized uptake value ratio (SUVR), clinical diagnosis, and amyloid status. A, Associations between flortaucipir cortical SUVR and global visual score (average of both readers). See Figure S2 in supporting information for similar plots using each reader's score or each cohort separately. The scatter plot shows a non‐linear relationship and differences in dynamic range between the two variables: SUVR values appear little related to the visual score in the low range, while the visual score reaches a ceiling for cortical SUVR >1.6. B, Association of clinical diagnosis and amyloid status on the global visual score (average of both readers). See Figure S3 in supporting information for similar plots using each reader's score. The YC group was excluded from the general linear model results described on top (overall R2 = .68). YC, young controls; OC, older controls; MCI, patients with a clinical diagnosis of mild cognitive impairment; ADdem, patients with a clinical diagnosis of Alzheimer's disease dementia; non‐AD, patients with a clinical diagnosis of non‐AD neurodegenerative syndrome. Clinical diagnoses were independent from amyloid beta (Aβ) status. C, Receiver operating curve (ROC) testing on the ability of the average global visual score and two SUVR measures to distinguish Aβ‐positive MCI/ADdem (n = 90) versus all other participants (n = 184). AUROC, area under the receiver operating curve Associations were comparable in strength and shape when looking at each reader's or cohort separately (Figure S2 in supporting information).

Visual scores across diagnoses and Aβ status

All but one YC had an average visual score of 0; this group was not included in the following analyses (Figure 3B). In the rest of the cohort, clinical group, Aβ status, and their interaction accounted for 68% of the variance in visual scores. All clinical groups differed from each other (all Ps <.02 using Holm Bonferroni correction), except for the OC and the non‐AD group (P = .65). Within Aβ‐negative participants, no clinical group difference was significant (all Ps >.5), whereas differences existed between clinical subgroups in Aβ‐positive cases (ADdem> MCI > OC and non‐AD). Within the OC, MCI, and ADdem clinical groups, visual scores were higher in Aβ‐positive patients compared to their Aβ‐negative counterparts (P = .023, P <.001 and P <.001, respectively). Receiver operating curve (ROC) analyses were run to assess the ability of the visual score to distinguish the Aβ‐negative from Aβ‐positive cases; area under the ROC curve (AUROC) were .62 in OC, .77 in MCI, and .97 in ADdem. The visual score was not related to Aβ status in the non‐AD clinical group (AUROC = .47, p = .79). Group comparisons were comparable when looking at each reader separately or considering the global cortical SUVR (Figure S3 in supporting information). ROC analyses were used to determine the ability of the global visual score to distinguish Aβ‐positive patients with MCI and AD dementia from all other cases (Figure 3C). The average visual score (AUROC = .92 [.88 to .94]; optimal threshold value ≥7, resulting in 78% sensitivity and 96% specificity) performed significantly better (P = .003) than global cortical SUVR (AUROC = .83 [.78 to .88]), but not different (P = .39) from the temporal meta‐ROI (AUROC = .90 [.85, .93]).

FTP patterns: reliability

When attributing one of the four visual patterns to the 274 scans, the two readers agreed in 66% of cases (compared to a chance level of 25% based on a four‐pattern classification), and the inter‐rater reliability was moderate (κ = .53 [.45 to .61]; Figure 4A and Figure S4 in supporting information). Most disagreement cases (26% of all scans) were “minor.” “Major” disagreements (ie, between AD‐like patterns and either negative or non‐AD patterns) happened in 22 (8%) of all scans (Figure S5 in supporting information). Out of these 18 discrepant cases, 7 were OC (6 Aβ‐negative), 10 were MCI (9 Aβ‐negative) and 1 was an Aβ‐negative case with a non‐AD diagnosis; 16/18 were ADNI cases.
FIGURE 4

Visual flortaucipir patterns: inter‐reader agreement and association with temporal meta‐region of interest (ROI) standardized uptake value ratio (SUVR). A, Inter‐reader agreement based on the 274 cases read by both readers (see Results section and Figure S4 in supporting information for statistics or intra‐rater agreement). Numbers in the table indicate raw numbers of cases for each combination of patterns, not percentages. B, Association between temporal meta‐ROI FTP‐SUVR and visual patterns defined by each reader in the 274 cases. The dotted line represents the independently defined 1.27 threshold; small bar graph inserts indicate the proportion of case above the 1.27 threshold within each visual pattern group. Black bars show median and interquartile range. P values correspond to post hoc tests following significant Kruskal Wallis test (see Result section). C, Similar to panel (B), but distinguishing cases for which the two readers assigned similar or different visual patterns. The right panel shows the discriminative characteristics of the visual AD‐like pattern to separate the two groups.

Visual flortaucipir patterns: inter‐reader agreement and association with temporal meta‐region of interest (ROI) standardized uptake value ratio (SUVR). A, Inter‐reader agreement based on the 274 cases read by both readers (see Results section and Figure S4 in supporting information for statistics or intra‐rater agreement). Numbers in the table indicate raw numbers of cases for each combination of patterns, not percentages. B, Association between temporal meta‐ROI FTP‐SUVR and visual patterns defined by each reader in the 274 cases. The dotted line represents the independently defined 1.27 threshold; small bar graph inserts indicate the proportion of case above the 1.27 threshold within each visual pattern group. Black bars show median and interquartile range. P values correspond to post hoc tests following significant Kruskal Wallis test (see Result section). C, Similar to panel (B), but distinguishing cases for which the two readers assigned similar or different visual patterns. The right panel shows the discriminative characteristics of the visual AD‐like pattern to separate the two groups. When separating the AD‐like pattern versus any of the other three, the two readers agreed on 237/274 cases (86% agreement, κ = .71 [.62 to .79]). Other patterns had low inter‐reader agreement (κ values between .39 and .49, Figure S4). Intra‐rater reliability was good for the definition of the four patterns (κreader1 = .76 and κreader2 = .79) and for the distinction of the AD‐like pattern versus any other (κreader1 = .73, κreader2 = .92; Figure S4).

FTP patterns: relationship to visual scores and SUVR quantification

Negative patterns had very minimal global visual scores, showing a floor effect (>95% of cases with a score of 0). In contrast, AD‐like patterns had median visual scores of 12 (reader 1) or 13 (reader 2), and reached a ceiling effect >40% cases being assigned the maximal score (Figure S6 in supporting information). The other patterns were associated with mildly elevated visual scores with median values around 2 to 4. Figure 4B shows that FTP temporal meta‐ROI SUVR strongly varied across visually determined patterns (Kruskal‐Wallis test: P <.001, ε2 = .49 for reader 1; P <.001, ε2 = .62 for reader 2). For both readers, the AD‐like pattern was associated with higher SUVR values than all other patterns (all Ps <.001). Most patients classified as AD‐like had temporal SUVR values > 1.27 (77% for reader 1, 96% for reader 2). The mild temporal and the non‐AD‐like patterns showed significantly increased temporal meta‐ROI SUVR values compared to the negative pattern at the group level, but most cases remained <1.27 (Figure 4B). When readers disagreed, reader 2 was usually more conservative than reader 1 (Figure 4A), and 20% to 23% of these scans had temporal meta‐ROI SUVR >1.27 threshold (Figure 4C).

FTP patterns across diagnoses and Aβ status

Both readers assigned the AD‐like pattern to the vast majority of Aβ‐positive ADdem, around half Aβ‐positive MCI, and some Aβ‐positive OC (reader 1: 12/31, reader 2: 8/31; Figure 5A). Reader 1 also classified 24 Aβ‐negative participants (10/72 controls, 11/28 MCI, 1/3 ADdem, and 2/23 non‐AD) as AD‐like. The mild temporal binding visual pattern was found across clinical and Aβ subgroups. The non‐AD‐like FTP pattern was assigned to a large proportion of patients with a non‐AD diagnosis (11/30 for reader 1, 16/30 for reader 2) and some cases from other subgroups. Similarly, tau‐positive cases defined by temporal meta‐ROI SUVR > 1.27, were found in the majority of Aβ‐positive patients with MCI and ADdem, but also in other subgroups (Figure 5B).
FIGURE 5

Visual flortaucipir patterns and temporal meta‐ region of interest (ROI) standardized uptake value ratio (SUVR)‐based positivity: association with clinical diagnosis and amyloid beta (Aβ) status. A, FTP visual patterns assigned by each reader (top: reader 1, bottom: reader 2; see Figure S7 in supporting information for the combination of both readers). Left panel shows the distribution of visual patterns (y‐axis) in each clinical/Aβ subgroup (x‐axis). Middle panel illustrates the distribution of visual patterns (specific n indicated on top), and in Aβ‐positive MCI/ADdem (circled numbers, n = 90 total), versus all other participants (n = 184). Right panel shows the diagnostic properties of the AD‐like pattern to distinguish Aβ‐positive MCI/ADdem versus all other participants. (B) Distribution of temporal SUVR‐based tau‐positive cases in each clinical/Aβ subgroup (left) or in Aβ‐positive MCI/ADdem versus all other participants (middle). Right panel shows the diagnostic properties of the meta‐ROI‐based FTP positivity to distinguish Aβ‐positive MCI/ADdem versus all other participants. YC, young controls; OC, older controls; MCI, patients with a clinical diagnosis of mild cognitive impairment; ADdem, patients with a clinical diagnosis of Alzheimer's disease dementia; non‐AD, patients with a clinical diagnosis of non‐AD neurodegenerative syndrome. Clinical diagnoses were independent from Aβ status. Se, sensitivity; Sp, specificity; PLR, positive likelihood ratio; NLR, negative likelihood ratio

Visual flortaucipir patterns and temporal meta‐ region of interest (ROI) standardized uptake value ratio (SUVR)‐based positivity: association with clinical diagnosis and amyloid beta (Aβ) status. A, FTP visual patterns assigned by each reader (top: reader 1, bottom: reader 2; see Figure S7 in supporting information for the combination of both readers). Left panel shows the distribution of visual patterns (y‐axis) in each clinical/Aβ subgroup (x‐axis). Middle panel illustrates the distribution of visual patterns (specific n indicated on top), and in Aβ‐positive MCI/ADdem (circled numbers, n = 90 total), versus all other participants (n = 184). Right panel shows the diagnostic properties of the AD‐like pattern to distinguish Aβ‐positive MCI/ADdem versus all other participants. (B) Distribution of temporal SUVR‐based tau‐positive cases in each clinical/Aβ subgroup (left) or in Aβ‐positive MCI/ADdem versus all other participants (middle). Right panel shows the diagnostic properties of the meta‐ROI‐based FTP positivity to distinguish Aβ‐positive MCI/ADdem versus all other participants. YC, young controls; OC, older controls; MCI, patients with a clinical diagnosis of mild cognitive impairment; ADdem, patients with a clinical diagnosis of Alzheimer's disease dementia; non‐AD, patients with a clinical diagnosis of non‐AD neurodegenerative syndrome. Clinical diagnoses were independent from Aβ status. Se, sensitivity; Sp, specificity; PLR, positive likelihood ratio; NLR, negative likelihood ratio The AD‐like pattern was far more frequent in Aβ‐positive patients with MCI or ADdem (reader 1: 73/90, reader 2: 70/90) than in other participants (reader 1: 36/184, reader 2: 12/184). Very few Aβ‐positive MCI/ADdem patients were attributed a negative FTP pattern (reader 1: 3/90, reader 2: 7/90). The AD‐like pattern had good sensitivity (78% to 81%) and good to excellent specificity (80% to 93%). The AD‐like pattern had discrimination performances comparable to tau‐positivity based on the temporal meta‐ROI (82% sensitive, 84% specific). When combining patterns from both readers, the sensitivity of the AD‐like pattern (assigned by both readers) was 76%, but the specificity reached 95% (Figure S7 in supporting information).

DISCUSSION

We developed and evaluated a visual rating scheme for the qualitative interpretation of FTP tau PET scans. We tested this visual approach in a large (n = 274) multi‐center sample derived from two independent cohorts with the intent to encompass the heterogeneity of a clinical population, including various clinical diagnoses, and using PET data acquired from multiple sites. Results indicate that the interpretation of FTP‐PET using our visual approach strongly correlates with SUVR quantification and is associated with both clinical diagnosis and Aβ status. The two visual indices we developed (a global score and an AD‐like visual pattern) allowed the differentiation of Aβ‐positive patients with MCI or AD dementia from other participants with ≈78% to 81% sensitivity and ≈80% to 96% specificity, similar to the SUVR quantification of temporal signal in the same patients (82% sensitivity, 84% specificity, Figure 4). Regardless of the visual or SUVR‐based definition of tau‐positivity, “false negative” cases were more likely to be patients with MCI rather than dementia, suggesting that all the proposed approaches might be missing lower tau levels. Alternatively, it is possible that these FTP‐negative Aβ‐positive patients did not harbor any tau pathology, and that their cognitive deficits might be due to another, non‐AD etiology. Within the recently proposed NIA‐AA research framework these cases would correspond to “Alzheimer's pathologic change [rather than full AD, ie, Aβ + tau] with concomitant non‐AD changes.” Inversely, a visual score >7 or an AD‐like visual pattern was frequent in Aβ‐positive OC, suggesting that the visual indices might capture tau pathology at the preclinical stage. Additional “false positives” visual FTP‐PET were found in Aβ‐negative subjects (Figure 4A and Figure S4). The absence of autopsy data does not make it possible to clarify whether these are true false positive (ie, positive FTP‐PET without underlying tau pathology, for instance due to off‐target binding related to non‐tau related pathology), or represent neurofibrillary pathology without Aβ such as primary age‐related tauopathy (PART). , In parallel to our effort, Avid Radiopharmaceuticals developed an independent visual read method consisting of separating FTP‐PET patterns consistent with AD (with an additional distinction between moderate and advanced patterns) from negative tau patterns. Although no formal head‐to‐head comparison was conducted between the two protocols, our approach and Avid's largely overlap, and equivalences can be drawn. Avid's “negative” patterns encompasses our own negative pattern (pattern I, no activity beyond off‐target binding), scans with isolated mesial temporal signal (our pattern II), or increased neocortical activity isolated to the anterolateral temporal, and/or frontal regions (which would likely be classified as our pattern IV). The identification of Avid's “AD tau pattern” relied on the presence of signal in the postero‐lateral temporal or occipital regions, and additional parietal/precuneus or frontal binding to be considered “advanced” AD tau; these scans would be classified as “AD‐like” (pattern III) in our visual scheme. The Avid protocol was used in an end‐of‐life study including 64 terminally ill participants (50 cognitively impaired, 14 cognitively unimpaired) to test the association with post mortem findings. An AD‐tau pattern was able to detect patients with Braak stages V‐VI (sensitivity: 92.3% to 100% based on five independent readers; specificity: 52% to 92%) or with high AD neuropathologic changes (ie, Braak V/VI with high amyloid burden; sensitivity: 94.7% to 100%; specificity: 50% to 92.3%). In patients with lower stages of tau pathology, most scans were read as negative (75% for Braak IV, 80% for Braak III, and 83% for Braak II), and were all negative based on image quantification. In spite of differences in the visual approach and study design, our data support the Avid results as both studies show that FTP‐PET visual reads are feasible and can identify advanced cases of clinically diagnosed AD or AD pathology with high accuracy, supporting its translation to clinical settings. Both the global visual score and the identification of an AD‐like pattern showed satisfactory intra‐ and inter‐reader agreement (κs≥.71), underlying the reliability of the main indices of interest derived from our visual rating scheme. However, the inter‐reader agreement on the four‐pattern distinction was only moderate (κ = .53), and the specific identification of the other three patterns were suboptimal (κs <.5, Figure S4). The difficulty in identifying these patterns can be attributed to the lower amount of tau‐PET signal associated with them and to some choices we made in the development of our visual scheme. Indeed, our approach focused on identifying a clear AD pattern based on a limited amount of data and data pre‐processing so that it could be widely used in clinical settings. Thus, visual reads were based on non‐intensity normalized images to make this method accessible without a dedicated neuro‐software. Image contrast had to be adjusted by each reader, resulting in intrinsic variability in the reading process. Alternative approaches, including reading an SUVR image or access to more advanced quantification methods might help improve pattern identification, similar to recent propositions related to amyloid‐PET. , , Morever, FTP scans were read without any anatomical (computed tomography or MRI) imaging while it could help PET interpretation by highlighting macroscopic abnormalities (eg, atrophy, presence of vascular anomalies or meningiomas ) and better distinguishing cortical from white matter signal. Anatomical information could be particularly helpful to identify the “mild temporal binding” pattern due to the multiple sources of off‐target binding in areas surrounding the areas of interest (eg, choroid plexus, extra‐brain binding medial to the entorhinal/parahippocampal areas). The main challenge to an FTP‐PET visual approach relates to the presence of off‐target binding. Contrary to amyloid imaging, in which white matter predominant off‐target binding is a constant feature across Aβ‐radiotracers and is key to the interpretation of PET scan, , , FTP‐PET off‐target signal is highly variable across individuals both in terms of severity and location, therefore hampering FTP interpretability. Various factors are suspected to cause FTP off‐target signal: neuromelanin, brain hemorrhagic lesions, as well as the enzymes monoamine oxidase (MAO)‐A and MAO‐B, widely distributed throughout the brain in both neurons and glia. , , Sufficient knowledge of FTP off‐target binding is fundamental to the accurate interpretation of the PET scan. In conclusion, our study shows that the proposed visual ratings of FTP tau PET are strongly associated with FTP global quantification, clinical diagnosis, and Aβ status across two cohorts with distinct characteristics. Inter‐ and intra‐rater reliability for our main indices are strong, suggesting that the proposed visual approach is reproducible and can therefore represent a promising alternative to quantitative tau‐PET measurements in clinical settings. Future investigations including an autopsy‐based gold standard are needed to evaluate the exact potential and limitations of tau‐PET and the proposed visual read approach.

FUNDING INFORMATION

Data collection and sharing for this project was funded by the Alzheimer's Disease Neuroimaging Initiative (ADNI; National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH‐12‐2‐0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie; Alzheimer's Association; Alzheimer's Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol‐Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann‐La Roche Ltd and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC; Johnson & Johnson Pharmaceutical Research & Development LLC; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer's Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California. The BACS/UCSF cohorts receive funding from National Institute on Aging grants (R01‐AG045611.; R01‐AG034570; R01‐AG062542; P50‐AG023501; P01‐AG19724 to) and the Rainwater Charitable Foundation. Renaud La Joie n funded by the Alzheimer's Association (AARF:16‐443577) and the National Institute of Aging (K99AG065501). Avid Radiopharmaceuticals enabled use of the [18F]Flortaucipir tracer by providing precursor, but did not provide direct funding and was not involved in data analysis or interpretation.

ETHICAL APPROVAL

The Institutional Review Boards of all participating institutions approved the study and informed consent was obtained from all participants or authorized representatives.  The study procedures were in accordance with the institutional ethical standards and with the 1964 Helsinki declaration and its later amendments.

CONFLICTS OF INTEREST

Suzanne L. Baker serves as a consultant for Genentech. No other potential conflicts of interest relevant to this article exist. Gil D. Rabinovici has research support from Avid Radiopharmaceuticals, GE Healthcare, Eli Lilly, and Life Molecular Imaging; has served as a consultant for Axon Neurosciences, GE Healthcare, Eisai, and Merck; and is an Associate Editor for JAMA Neurology. No other potential conflicts of interest relevant to this article exist. William J. Jagust serves as a consultant to Bioclinica, Genentech, CuraSen, Biogen, and Grifols and owns equity in OptoCeutics. No other potential conflicts of interest relevant to this article exist. Susan M. Landau has served as a consultant for NeuroVision and Cortexyme. No other potential conflicts of interest relevant to this article exist. No potential conflicts of interest relevant to this article exist for all other authors. Supporting Information Click here for additional data file.
  50 in total

1.  The diagnosis of dementia due to Alzheimer's disease: recommendations from the National Institute on Aging-Alzheimer's Association workgroups on diagnostic guidelines for Alzheimer's disease.

Authors:  Guy M McKhann; David S Knopman; Howard Chertkow; Bradley T Hyman; Clifford R Jack; Claudia H Kawas; William E Klunk; Walter J Koroshetz; Jennifer J Manly; Richard Mayeux; Richard C Mohs; John C Morris; Martin N Rossor; Philip Scheltens; Maria C Carrillo; Bill Thies; Sandra Weintraub; Creighton H Phelps
Journal:  Alzheimers Dement       Date:  2011-04-21       Impact factor: 21.566

2.  Association of lifetime cognitive engagement and low β-amyloid deposition.

Authors:  Susan M Landau; Shawn M Marks; Elizabeth C Mormino; Gil D Rabinovici; Hwamee Oh; James P O'Neil; Robert S Wilson; William J Jagust
Journal:  Arch Neurol       Date:  2012-05

3.  Correlation of In Vivo [18F]Flortaucipir With Postmortem Alzheimer Disease Tau Pathology.

Authors:  Ruben Smith; Moa Wibom; Daria Pawlik; Elisabet Englund; Oskar Hansson
Journal:  JAMA Neurol       Date:  2019-03-01       Impact factor: 18.302

4.  Regional profiles of the candidate tau PET ligand 18F-AV-1451 recapitulate key features of Braak histopathological stages.

Authors:  Adam J Schwarz; Peng Yu; Bradley B Miller; Sergey Shcherbinin; James Dickson; Michael Navitsky; Abhinay D Joshi; Michael D Devous; Mark S Mintun
Journal:  Brain       Date:  2016-03-02       Impact factor: 13.501

Review 5.  Strategic roadmap for an early diagnosis of Alzheimer's disease based on biomarkers.

Authors:  Giovanni B Frisoni; Marina Boccardi; Frederik Barkhof; Kaj Blennow; Stefano Cappa; Konstantinos Chiotis; Jean-Francois Démonet; Valentina Garibotto; Panteleimon Giannakopoulos; Anton Gietl; Oskar Hansson; Karl Herholz; Clifford R Jack; Flavio Nobili; Agneta Nordberg; Heather M Snyder; Mara Ten Kate; Andrea Varrone; Emiliano Albanese; Stefanie Becker; Patrick Bossuyt; Maria C Carrillo; Chiara Cerami; Bruno Dubois; Valentina Gallo; Ezio Giacobini; Gabriel Gold; Samia Hurst; Anders Lönneborg; Karl-Olof Lovblad; Niklas Mattsson; José-Luis Molinuevo; Andreas U Monsch; Urs Mosimann; Alessandro Padovani; Agnese Picco; Corinna Porteri; Osman Ratib; Laure Saint-Aubert; Charles Scerri; Philip Scheltens; Jonathan M Schott; Ida Sonni; Stefan Teipel; Paolo Vineis; Pieter Jelle Visser; Yutaka Yasui; Bengt Winblad
Journal:  Lancet Neurol       Date:  2017-07-11       Impact factor: 44.182

6.  Alzheimer's pathology targets distinct memory networks in the ageing brain.

Authors:  Anne Maass; David Berron; Theresa M Harrison; Jenna N Adams; Renaud La Joie; Suzanne Baker; Taylor Mellinger; Rachel K Bell; Kaitlin Swinnerton; Ben Inglis; Gil D Rabinovici; Emrah Düzel; William J Jagust
Journal:  Brain       Date:  2019-08-01       Impact factor: 13.501

7.  Tau pathology and neurodegeneration contribute to cognitive impairment in Alzheimer's disease.

Authors:  Alexandre Bejanin; Daniel R Schonhaut; Renaud La Joie; Joel H Kramer; Suzanne L Baker; Natasha Sosa; Nagehan Ayakta; Averill Cantwell; Mustafa Janabi; Mariella Lauriola; James P O'Neil; Maria L Gorno-Tempini; Zachary A Miller; Howard J Rosen; Bruce L Miller; William J Jagust; Gil D Rabinovici
Journal:  Brain       Date:  2017-12-01       Impact factor: 13.501

8.  Alzheimer's Disease Neuroimaging Initiative (ADNI): clinical characterization.

Authors:  R C Petersen; P S Aisen; L A Beckett; M C Donohue; A C Gamst; D J Harvey; C R Jack; W J Jagust; L M Shaw; A W Toga; J Q Trojanowski; M W Weiner
Journal:  Neurology       Date:  2009-12-30       Impact factor: 9.910

9.  Augmenting Amyloid PET Interpretations With Quantitative Information Improves Consistency of Early Amyloid Detection.

Authors:  Nicholas R Harn; Suzanne L Hunt; Jacqueline Hill; Eric Vidoni; Mark Perry; Jeffrey M Burns
Journal:  Clin Nucl Med       Date:  2017-08       Impact factor: 7.794

10.  Existing Pittsburgh Compound-B positron emission tomography thresholds are too high: statistical and pathological evaluation.

Authors:  Sylvia Villeneuve; Gil D Rabinovici; Brendan I Cohn-Sheehy; Cindee Madison; Nagehan Ayakta; Pia M Ghosh; Renaud La Joie; Samia Kate Arthur-Bentil; Jacob W Vogel; Shawn M Marks; Manja Lehmann; Howard J Rosen; Bruce Reed; John Olichney; Adam L Boxer; Bruce L Miller; Ewa Borys; Lee-Way Jin; Eric J Huang; Lea T Grinberg; Charles DeCarli; William W Seeley; William Jagust
Journal:  Brain       Date:  2015-05-06       Impact factor: 13.501

View more
  6 in total

1.  Diagnostic value of amyloid-PET and tau-PET: a head-to-head comparison.

Authors:  Daniele Altomare; Camilla Caprioglio; Frédéric Assal; Gilles Allali; Aline Mendes; Federica Ribaldi; Kelly Ceyzeriat; Marta Martins; Szymon Tomczyk; Sara Stampacchia; Alessandra Dodich; Marina Boccardi; Christian Chicherio; Giovanni B Frisoni; Valentina Garibotto
Journal:  Eur J Nucl Med Mol Imaging       Date:  2021-02-27       Impact factor: 9.236

2.  Evaluation of a visual interpretation method for tau-PET with 18F-flortaucipir.

Authors:  Ida Sonni; Orit H Lesman Segev; Suzanne L Baker; Leonardo Iaccarino; Deniz Korman; Gil D Rabinovici; William J Jagust; Susan M Landau; Renaud La Joie
Journal:  Alzheimers Dement (Amst)       Date:  2020-11-28

3.  Deep learning improves utility of tau PET in the study of Alzheimer's disease.

Authors:  James Zou; David Park; Aubrey Johnson; Xinyang Feng; Michelle Pardo; Jeanelle France; Zeljko Tomljanovic; Adam M Brickman; Devangere P Devanand; José A Luchsinger; William C Kreisl; Frank A Provenzano
Journal:  Alzheimers Dement (Amst)       Date:  2021-12-31

Review 4.  Tau Biomarkers in Dementia: Positron Emission Tomography Radiopharmaceuticals in Tauopathy Assessment and Future Perspective.

Authors:  Maria Ricci; Andrea Cimini; Riccardo Camedda; Agostino Chiaravalloti; Orazio Schillaci
Journal:  Int J Mol Sci       Date:  2021-11-30       Impact factor: 5.923

5.  [18F]THK-5351 PET Patterns in Patients With Alzheimer's Disease and Negative Amyloid PET Findings.

Authors:  Minyoung Oh; Jungsu S Oh; Seung Jun Oh; Sang Ju Lee; Jee Hoon Roh; Woo Ram Kim; Ha-Eun Seo; Jae Myeong Kang; Sang Won Seo; Jae-Hong Lee; Duk L Na; Young Noh; Jae Seung Kim
Journal:  J Clin Neurol       Date:  2022-07       Impact factor: 2.566

6.  Visually Identified Tau 18F-MK6240 PET Patterns in Symptomatic Alzheimer's Disease.

Authors:  Natasha Krishnadas; Kun Huang; Stephanie A Schultz; Vincent Doré; Pierrick Bourgeat; Anita M Y Goh; Fiona Lamb; Svetlana Bozinovski; Samantha C Burnham; Joanne S Robertson; Simon M Laws; Paul Maruff; Colin L Masters; Victor L Villemagne; Christopher C Rowe
Journal:  J Alzheimers Dis       Date:  2022       Impact factor: 4.160

  6 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.