Abstract
Background
Endogenous Alu RNAs form double-stranded RNAs recognized by double-stranded RNA sensors and activate IRF and NF-kB transcriptional paths and innate immunity. Deamination of adenosines to inosines by the ADAR family of enzymes, a process termed A-to-I editing, disrupts double-stranded RNA structure and prevents innate immune activation. Innate immune activation is observed in Alzheimer's disease, the most common form of dementia. We have previously reported loss of A-to-I editing in hippocampus vasculature, but no change in cortex or cortex vasculature, associated with Alzheimer's disease.
Objective
Here, we investigated the status of Alu RNA A-to-I editing in cortex extracellular vesicles in Alzheimer's disease.
Methods
We used existing RNA-seq data sets and the SPRINT software package to determine levels of Alu RNA A-to-I editing in cortex extracellular vesicles in Alzheimer's disease and control groups and compared these editing profiles to those found in both total cortex and hippocampus vasculature.
Results
We find substantial loss of Alu A-to-I editing in cortex extracellular vesicles in Alzheimer's disease. By measuring editing patterns on a gene-by-gene basis, we determined that editing patterns in cortex extracellular vesicles resemble editing patterns in hippocampus vasculature rather than total cortex.
Conclusions
We conclude that hippocampus vasculature unedited Alu RNAs are packaged in extracellular vesicles, travel to the cortex, deliver their cargo and stimulate innate immunity and alter other basic biological processes contributing to Alzheimer's disease progression.
Keywords
Introduction
Deamination of adenosines in double-stranded RNA to yield inosines catalyzed by the adenosine deaminase acting on RNA (ADAR) family of enzymes is a common post-transcriptional modification of RNA in eukaryotic cells.1–6 This process is commonly referred to as ‘A-to-I editing’. In humans, the most frequent editing sites are found in RNAs transcribed from Alu retrotransposons.7–9 Alu genomic elements are primate specific and arose from insertion into the genome of a head-to-tail fusion of 7SL RNA. 10 About 1 million copies of Alu genomic sequences, which are about 300 nucleotides in length, exist in the human genome, thus Alu elements make up about 10% of the human genome. 11 RNAs containing Alu elements are transcribed by RNA polymerase 2 as part of pre-mRNAs or by RNA polymerase 3 as part of their normal life cycle. Alu RNAs transcribed from non-coding genome space, gene introns, intergenic space, 3’ and 5’ untranslated regions (UTR) are the most common sites of A-to-I editing while editing of Alu RNAs in exons is relatively less frequent.12,13
In general terms RNA A-to-I editing affects several biological processes.7,8,14 Editing of RNAs transcribed from introns can produce alternative splicing of mRNAs. Further, looping of Alu containing enhancers to promoters can alter promoter-enhancer selectivity and impact expression of cognate genes. 15 Editing of Alu containing elements in 3’ UTRs may alter microRNA binding sights that may also affect mRNA stability. 16 Editing may also affect additional functions encoded by nucleotide sequences within the 5’ UTR such as ribosome recruitment to the mRNA and start codon choices. 17
Since Alu elements arose from a head-to-tail fusion of 7SL RNA,11–13 Alu RNAs can form double-stranded RNA (dsRNA) structures capable of activating host dsRNA sensors including RIG-I, MDA5, and TLR3.18–25 In addition, Alu genomic elements inserted close in the genome in opposite orientations are common and when transcribed, exist as two Alu RNA elements joined by a linker sequence, also referred to as Alu inverted repeats (IR/Alu), and IR/Alu RNAs are potent activators of the host dsRNA sensor, MDA5.18–20 Numerous lines of evidence indicate that A-to-I editing of Alu RNAs reduces their dsRNA character and prevents activation of dsRNA sensors. Activation of dsRNA sensors causes innate immune activation resulting from induction of interferon-stimulated genes (ISGs), cytokines, and other inflammatory mediators because of activation of transcription factors, interferon regulatory factors (IRF) and nuclear factor-kappa B (NF-kB). In humans, inactivating mutations in ADAR result in one form of Aicardi-Goutières syndrome, characterized by chronic innate immune activation resulting in severe encephalopathy that can lead to death or existence in a vegetative state. 26
Extracellular vesicles (EVs), also referred to as exosomes are small vesicles (30–150 nanometers in diameter) produced in the endosomal compartment and secreted by most eukaryotic cells.27–31 EVs contain proteins, nucleic acids, lipids, and metabolites and are released into a variety of biological fluids, including blood, urine, saliva, breast milk, synovial fluid, amniotic liquid, and cerebrospinal fluid. As such, they mediate both near and long-distance intercellular communication in both health and disease and affect a variety of aspects of cell biology.
Loss of Alu RNA A-to-I editing in blood occurs as a result of relapsing remitting multiple sclerosis (MS) and is accompanied by accumulation of Alu dsRNAs and induction of IRF and NF-kB regulated genes.25,32 Loss of Alu RNA A-to-I editing also results from COVID-19 disease (blood and lung) and severe influenza disease (blood) and in cell culture following infection with SARS-CoV-2.33,34 Thus, loss of endogenous Alu RNA A-to-I editing may contribute to both acute and chronic inflammation observed in these conditions.
Alzheimer's disease (AD) is a progressive disease and the most common form of dementia. 35 Its causes are incompletely understood but AD is associated with accumulation of aggregated misfolded proteins including amyloid-β plaques and neurofibrillary tangles. We find loss of Alu RNA A-to-I editing in hippocampus vasculature (HPCvasc) associated with AD. 36 In contrast, AD-dependent changes in editing are not observed in prefrontal cortex (CTX) or CTX vasculature. Here, we examined levels of A-to-I edited Alu RNAs in EVs isolated from CTX and find reduced levels of Alu RNA editing in AD samples compared to normal cognition (NCI) samples. Gene-specific editing patterns in CTX EVs are more like patterns observed in HPCvasc than CTX leading us to conclude that CTX EVs may have arisen from HPCvasc and travelled to CTX. CTX EVs carrying reduced levels of edited Alu RNAs may fuse with cells of the CTX to activate dsRNA sensors to create an inflammatory state or affect other cellular processes in CTX.
Methods
RNA-Seq FASTQ files
RNA-seq FASTQ files measuring RNA levels in EVs isolated from CTX used here have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE197505: AD (N = 8), and age and gender-matched controls (normal cognition, NCI) (N = 10). 37 RNA-seq FASTQ files measuring RNA levels in HPCvasc are accessible through accession number GSE163577 and have been previously described. 38 RNA-seq FASTQ files measuring RNA levels in total pre-frontal CTX are accessible via accession number GSE53699 and have been previously described: AD group; N = 8, and NCI group; N = 8. 39
Briefly, EVs were purified from pre-frontal CTX from postmortem brain tissues using standard procedures in accordance with the Minimal Information for Studies of Extracellular Vesicles 2018 guidelines, total RNA was purified and subjected to RNA sequencing using standard procedures. The human brain tissues in this study were provided by the National Human Brain Bank for Development and Function, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China, from neuropathologically classified individuals with AD or NCI at the time of death and brain donation. 37
Analysis of RNA A-to-I editing
Our methods for analysis of RNA A-to-I editing have been described.33–35 Briefly, bulk RNA-seq FASTQ files were trimmed, and we used the following workflow to identify endogenous RNA A-to-I–editing sites from paired FASTQ sequencing files. The main identification tool was a python-based package called the SPRINT toolkit accepting sequence files and producing text files with the following information for each edit site: (1) genomic location; (2) type of edit (e.g., A-to-G; T-to-C), strand (+ or −); (3) number of edits per site and total number of reads per site.40,41 Mathematica programs were developed to synthesize data: numbers of samples in groups with shared editing sites, mean numbers of total reads and edits for each editing site, editing sites common and unique to groups, and editing sites per gene. This information was tied to an Alu database to annotate each site: gene locations (intronic, noncoding RNA, intergenic, 3′ and 5′ untranslated regions (UTRs), exons) and if editing sites were within Alu or non-Alu elements. 42 To avoid sequencing errors, editing sites were only included if total reads were greater than 5 and edit/read ratios were greater than 0.05. To create genome-wide A-to-I editing indices (AEI), we identified all A-to-I editing sites present in one study group and summed edit/read ratios for all editing sites across the genome among study groups. Gene expression levels from RNA-seq FASTQ files were determined using the DESeq2 R package as described. 43 Both raw p-values and p-values after adjusting for multiple testing, false discovery rate (FDR), were determined.
Results
AD-associated loss of a-to-i editing of Alu RNAs in pre-frontal cortex EVs
We obtained RNA-seq FASTQ files derived from a study that comprised isolation of EVs from frontal CTX of brain from individuals with normal cognition (NCI, N = 10) or AD (N = 8) at the time of death and brain donation (GSE197505). 36 Groups were age and gender matched and did not differ in the post-mortem time between brain harvest and RNA processing. We used the SPRINT software package to determine levels of RNA A-to-I editing among samples using RNA-seq data. Briefly, this package scans the genome for A-to-G mismatches compared to the reference genome after removing known single nucleotide polymorphisms and reports proportion of edits to total reads at each edited site. We required the proportion of edits/reads to exceed 0.05 and the number of total reads to exceed 5 to guard against sequencing errors. We also identified the number of times a given nucleotide was edited in each sample from separate groups to examine both rare and more common editing sites. To express these data, we calculated both total number of edited sites per group and created an A-to-I editing index (AEI) by multiplying average proportion of edits/reads at each site times the total number of samples in which a given site was edited in the NCI group or the AD group.
To initiate these studies, we used the SPRINT software package to determine the total number of expressed Alu RNA elements (reads ≥ 1) and the average number reads of all expressed Alu RNA elements in NCI and AD groups. We did not find these two measurements to be substantially different between the two groups (Table 1). We also employed the DESeq2 software package to determine total RNA counts in the two groups. Again, we did not find this measurement to be substantially different between NCI and AD groups. Thus, total number of Alu RNA elements detected and average expression levels of Alu RNA elements were similar in NCI and AD groups. Further, levels of total RNA detected by DESeq2 analysis were similar between NCI and AD groups.
Number and expression levels of RNA Alu elements in CTX EVs from NCI and AD groups determined from RNA-seq data.
Total # of Alu RNA elements and average ± standard deviation Alu RNA reads determined using the SPRINT software package in NCI and AD groups. Total RNA counts determined by DESeq2.
We found reduced numbers of total A-to-I edited sites in EVs isolated from the AD group compared to the NCI group (Figure 1(a)). We also found reduced overall AEI, which also takes into count the frequency of individual editing sites, in the AD group compared to the NCI group (Figure 1(b)). Overall, the majority of editing sites were found in introns followed by intergenic space, 3’ UTRs, noncoding RNAs (ncRNA), 5’ UTRs, and synonymous (syn) and nonsynonymous edits found in exons (Figure 1(c), Supplemental File 1). Loss of editing in the AD group compared to the NCI group was similar in RNAs transcribed from these different genomic regions. We conclude from these analyses that overall levels of A-to-I edited Alu RNAs in EVs were reduced in the AD group compared to the NCI group.

Reduced Alu RNA A-to-I editing in AD compared to NCI CTX EVs. (a) Number of editing sites (X-axis) found in CTX EV RNAs obtained from individuals with NCI (N = 10) or AD (N = 8) at the time of death and brain donation determined by SPRINT analysis of RNA-Seq data. Y-axis is the proportion of edits/reads at each site, p < 0.0001, Welch's t-test. (b) As in (a), except AEI was calculated for each editing site, p < 0.0001, Welch's t-test. (c) Genomic locations of Alu RNA A-to-I editing sites.
As a second approach to our analysis, we determined frequencies that individual editing sites were detected in two or more EV samples in NCI and AD groups. Overall, number of individual editing sites we detected declined as frequency of presence of editing sites increased from two to ten. Overall, we found loss of editing in the AD group compared to the NCI group at all frequencies (Figure 2(a)). This was also the case when we calculated AEI (Figure 2(b)). Thus, loss of editing was found when we examined editing sites that are rare in frequency and editing sites that are common in frequency in AD compared to NCI groups.

Number of times a given Alu RNA A-to-I editing site was detected in NCI or AD CTX EVs. (a) Frequency (≥2, X-axis) of editing sites in CTX EVs from NCI or AD samples determined from RNA-seq data, p < 0.0001, 2-way ANOVA. (b) As in (a), except AEI was calculated, p < 0.0001, 2-way ANOVA.
A-to-I editing patterns in CTX EVs resemble HPCvasc editing patterns
We determined AEI per gene in CTX EV samples from NCI and AD groups. For this analysis we included all editing sites traversing the 5’ UTR to the 3’ UTR including both exons and introns. We identified those genes that exhibited high levels of editing in EVs in the NCI group but were completely unedited in the AD group (Figure 3(a)). We have previously reported that A-to-I editing in total hippocampus vasculature (HPCvasc) is reduced in the AD group compared to the NCI group. 36 In contrast, we did not find differences in A-to-I editing in total frontal cortex between NCI and AD groups. Therefore, we compared editing patterns of these same genes in HPCvasc samples between NCI and AD groups. We found that editing patterns of these genes were uniformly lower in the HPCvasc AD group compared to the NCI group (Figure 3(b)). We performed a similar comparison of editing in frontal cortex. Here, we did not detect markedly different editing patterns of this group of genes when NCI and AD groups were compared (Figure 3(c)). In fact, most of these genes were un-edited in CTX samples in both NCI and AD groups. We conclude that the pattern of A-to-I editing of this group of genes in CTX EVs was similar to the editing pattern found in total HPCvasc and dissimilar from the editing pattern of these genes in total CTX.

Editing patterns of hyper-edited genes in CTX EV compared to total HPCvasc and total CTX. (a) AEI/gene was determined for indicated genes in NCI CTX EVs; editing of these genes was undetectable in AD CTX EVs. (b) As in (a) except AEI/gene in HPCvasc was determined. (c) As in (a) except AEI/gene in total CTX was determined.
Studies to analyze HPCvasc gene expression in NCI and AD groups employed single cell RNA sequencing (scRNA-seq). 38 Investigators used this approach to identify genes selectively expressed by different cell lineages within HPCvasc including astrocytes, neurons, oligodendrocytes, endothelial cells, different myeloid cells, T cells and others. We previously compared AEI of these lineage-defining genes between NCI and AD groups to ask if differences in editing in HPCvasc may be attributed to specific lineages. 36 In general terms, we found that AEI of lineage defining genes was reduced in AD groups compared to NCI groups leading us to conclude that loss of editing observed in HPCvasc was not cell-type specific but rather a general feature of cell lineages present within HPCvasc. Frontal CTX contains many of these same cell lineages, astrocytes, neurons, oligodendrocytes, etc. In general terms, EVs are packaged intracellularly, released from cells and migrate across extracellular space to reach target cells or tissues. Therefore, we asked if editing patterns of lineage-defining genes found in CTX EVs resembled editing patterns of these genes in HPCvasc or total CTX in NCI and AD groups. To illustrate our approach, we determined editing patterns of astrocyte lineage defining genes in EV, HPCvasc and CTX tissues by determining AEI and expressed data as the AD/NCI ratio, log2 (Figure 4(a)). Gene AEIs were segregated into three categories; 1) log2 ratio < -0.4 = reduced editing in AD compared to NCI, 2) ratio −0.4 to +0.4 = no difference in editing between AD and NCI groups, and 3) ratio > 0.4 = increased editing in AD compared to NCI groups. We determined if this ratio fell in the same category in cortex EVs and HPCvasc or in cortex EVs and CTX, 1 = ‘agree’ or 0 = ‘disagree’.

Editing patterns of astrocyte specific genes in CTX EVs, total CTX and HPCvasc. (a) Editing patterns (Y-axis) of indicated genes (X-axis) selectively expressed by astrocytes were determined in CTX EVs, HPCvasc and total CTX. (b) Gene AEIs were segregated into three categories; 1) log2 ratio < −0.4 = reduced editing in AD compared to NCI, 2) ratio −0.4 to +0.4 = no difference in editing between AD and NCI groups, and 3) ratio > 0.4 = increased editing in AD compared to NCI groups. We determined if this ratio fell in the same category in cortex EVs and HPCvasc or in cortex EVs and CTX, 1 = ‘agree’ or 0 = ‘disagree’. p = 0.003 determined by chi-square analysis assuming 33% of gene AEIs fell into each category by chance.
We determined statistical significance by chi-square analysis assuming 33% of gene AEIs should fall into each category by chance (Figure 4(b)). This comparison indicated that the editing pattern of astrocyte lineage defining genes in CTX EVs was similar to the editing patterns observed in HPCvasc but not CTX editing patterns. We performed a similar analysis of genes that define other cell lineages and obtained similar results (Table 2, Supplemental File 2). These studies support the notion that editing patterns of RNA in CTX EVs are more similar to RNA editing patterns in HPCvasc than CTX and raise the possibility that CTX EVs may not have been packaged in CTX but another brain region such as HPCvasc.
Comparison of lineage specific gene AEIs in CTX EVs to HPCvasc or total CTX.
Lineage = different cell types in brain. AEI per lineage-defining gene was determined from CTX EVs, HPCvasc, and CTX RNA-seq data representing 3 independent studies. AD AEI/NCI, log2, was determined for each gene per group. Gene AEIs were segregated into 3 categories: 1) ratio < -0.4 = reduced editing in AD compared to NCI groups, 2) ratio −0.4 to +0.4 = no difference in editing between AD and NCI groups, and 3) ratio > 0.4 = increased editing in AD compared to NCI groups. We determined if this ratio fell in the same category in CTX EVs and HPCvasc groups or in CTX EV and CTX groups = ‘agree’ or fell into one of the other categories = ‘disagree’. Numbers in the table are # of genes from each lineage that ‘agree’ or ‘disagree’ between groups. p values were determined by chi-square analysis assuming 33% of gene AEIs should fall into each category by chance. * = ‘# agree’ lower than expected by chance
Gene expression and pathways analysis of hyper-edited genes
Enhancer RNAs (eRNAs) transcribed from distal genomic enhancers regulate transcription of target genes by processes that are incompletely understood.44–50 One view is that eRNAs participate in chromatin looping via RNA:DNA;DNA hybrids, termed R-loops, to recruit enhancers to promoters resulting in gene activation. New studies demonstrate that Alu sequences in eRNAs interact with Alu sequences in complimentary promoters of target genes and alter rates of gene transcription. 15 We reasoned that A-to-I editing of Alu RNAs may affect ability of eRNAs to interact with their target promoters, either positively or negatively, and this may change transcript levels of target genes. Further, eRNAs packaged in EVs may be delivered to distal tissues and alter mRNA levels of target genes in tissues in which EVs travel to. As an initial test of this notion, we identified genes in EVs that exhibited the greatest level of editing in CTX EVs in the NCI group, hyper-edited genes, AEI > 50 (Figure 5(a)). We found that editing of these hyper-edited genes was uniformly reduced in the AD group. The vast majority of these edited Alu RNAs were derived from introns (Figure 1(c)).

AEI and expression profiles of hyper-edited genes. (a) Differences in AEI of hyper-edited genes between NCI and AD EVs, Y-axis is AEI/gene and X-axis is rank, p < 0.0001 determined by 2-way ANOVA. (b) NCI and AD expression differences of hyper-edited genes in total CTX, Y-axis shows expression ratios of hyper-edited genes (a), NCI/AD, log2, X-axis is rank, p < 0.0001 determined by 2-way ANOVA. (c) Pathways analysis of hyper-edited genes determined by NIAID DAVID Functional Annotation Bioinformatics Microarray Analysis program. Input was hyper-edited genes (a) and indicated pathways were over-represented. X-axis is Benjamini p value, -log10. Dashed line = p 0.05.
We next asked if mRNA levels of these same genes were different in total CTX among NCI and AD groups reasoning that CTX EVs may deliver edited or unedited Alu-containing eRNAs to CTX and influence levels of target gene mRNAs in CTX. Overall, we found most of the genes within this group exhibited greater levels of expression in total CTX in the NCI group compared to the AD group (Figure 5(b), Supplemental File 3). We submitted this gene list to the NIAID DAVID Functional Annotation Bioinformatics Microarray Analysis program for pathways analysis. 51 Pathways found to be significantly over-represented were 1) zinc finger proteins, 2) microtubules, 3) actin binding, 4) protein K48-linked ubiquitination, 5) ubiquitin conjugation, 6) isopeptide bonds, and 7) cross-inking glycyl-lysine isopeptide bonds (Figure 5(c)). Although we found that differences in expression levels of individual genes were not huge, it seems that small differences in expression of the large number of genes involved in these pathways may impact functions of these pathways.
Discussion
Our studies demonstrate elevated levels of unedited Alu RNAs in EVs isolated from CTX in AD samples compared to NCI samples. Reduced editing in AD is seen at both rare and common editing sites. Our previous studies show that reduced editing of Alu RNAs associated with AD is seen in HPCvasc but not total CTX or CTXvasc. 36 Our analysis of A-to-I editing of individual genes shows that gene editing patterns are more similar to editing patterns observed in HPCvasc than in CTX. Similarly, editing patterns of CTX EV cell lineage-specific genes are more similar to the editing patterns found in HPCvasc than in total CTX in both AD and NCI groups. A model consistent with these results is that EVs found in CTX are packaged with Alu RNAs from HPCvasc and travel to CTX to deliver their cargo.
Two lines of evidence support our notion that EVs isolated from CTX develop in HPCvasc and travel to CTX. First, genes that have high levels of editing in NCI CTX EVs but no detectable levels of editing in AD CTX EVs display a similarly high editing level in NCI HPCvasc and a reduced or undetectable editing level in AD HPCvasc. In contrast, these same genes are largely unedited in both NCI CTX and AD CTX, suggesting RNAs in CTX EVs do not come from CTX but rather come from HPCvasc. Second, editing patterns of cell lineage specific genes observed in CTX EVs, both from NCI and AD samples, look more like the gene editing patterns observed in HPCvasc than gene editing patterns observed in CTX. Further studies will be required to determine if RNAs from EVs located in other brain regions resemble HPCvasc editing patterns or exhibit unique or region-specific RNA editing patterns.
Neuroinflammation is implicated in the pathogenesis of AD.52–55 Loss of editing of Alu RNAs results in activation of IRF and NF-kB transcriptional paths and innate immune responses. Activation of these same transcriptional paths and induction of IRF and NF-kB responsive genes is also observed in AD.53–55 Thus, packaging of unedited Alu RNAs in CTX EVs associated with AD may contribute to the neuroinflammation observed in AD.
As outlined above, editing of Alu elements within eRNAs may affect eRNA function that, in turn, will affect gene transcript levels. 15 Therefore, it seems that a second possibility is that CTX EVs deliver unedited, in the case of AD, and edited, in the case of NCI, eRNAs to CTX that may alter gene expression programs. In support of this idea, we find that hyper-edited genes in CTX EVs display uniformly reduced editing in AD samples compared to NCI samples. Expression of this group of genes in CTX is largely reduced in AD samples compared to NCI samples. Pathway analysis shows that these genes fall into groups whose functions have been implicated in AD pathogenesis, such as microtubule and actin function, 56 ubiquitination, 57 and others. 58 Interestingly, our studies in multiple sclerosis result in a similar finding 25 ; reduced editing on a per gene basis is associated with reduced gene expression raising the possibility that this may be a more general finding.
In general terms, EVs are released from most eukaryotic cells and carry nucleic acids, proteins, lipids, and various other cellular products as cargo.27–31 EVs can deliver this cargo to either adjacent cells or cells at a distance and this is thought to represent a form of intercellular communication. Our results support the idea that CTX EVs may come from the HPCvasc given the similar A-to-I gene editing patterns that are observed. Loss of volume or atrophy of the HPC is an early sign of AD.59,60 Further, breakdown of the blood-brain barrier in the HPC is not only an early sign of AD but is also associated with mild cognitive impairment.60–64 Although somewhat speculative, it seems that accumulation of unedited Alu RNAs in HPCvasc, their packaging in EVs, and spreading to other regions of the brain may contribute to initiation and progression of AD-associated cognitive decline.
Footnotes
Acknowledgments
The authors have no acknowledgments to report.
Author contributions
Philip Crooke, 3rd (Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Validation; Writing – original draft; Writing – review & editing); John Tossberg (Data curation; Investigation; Validation; Writing – review & editing); Thomas Martin Aune (Conceptualization; Data curation; Formal analysis; Funding acquisition; Investigation; Project administration; Writing – original draft; Writing – review & editing).
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by grants from the NIH (R21AI144193 and R21AI164107 to T.M. Aune.
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Data availability
Supplemental data will be deposited in PubMed Central.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
