Next Article in Journal
AGC2 (Citrin) Deficiency—From Recognition of the Disease till Construction of Therapeutic Procedures
Previous Article in Journal
Protein–Protein Interactions Mediated by Intrinsically Disordered Protein Regions Are Enriched in Missense Mutations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

From Mouse to Human: Comparative Analysis between Grey and White Matter by Synchrotron-Fourier Transformed Infrared Microspectroscopy

by
Paula Sanchez-Molina
1,2,
Martin Kreuzer
3,
Núria Benseny-Cases
3,
Tony Valente
1,2,
Beatriz Almolda
1,2,
Berta González
1,2,
Bernardo Castellano
1,2 and
Alex Perálvarez-Marín
1,4,*
1
Institute of Neurosciences, Universitat Autònoma de Barcelona, 08193 Bellaterra, Barcelona, Spain
2
Department of Cell Biology, Physiology and Immunology, Universitat Autònoma de Barcelona, 08193 Bellaterra, Barcelona, Spain
3
ALBA Synchrotron Light Source, Carrer de la Llum 2-26, 08290 Cerdanyola del Vallès, Barcelona, Catalonia, Spain
4
Biophysics Unit, Department of Biochemistry and Molecular Biology, Universitat Autònoma de Barcelona, 08193 Bellaterra, Barcelona, Catalonia, Spain
*
Author to whom correspondence should be addressed.
Biomolecules 2020, 10(8), 1099; https://doi.org/10.3390/biom10081099
Submission received: 10 July 2020 / Accepted: 22 July 2020 / Published: 24 July 2020

Abstract

:
Fourier Transform Infrared microspectroscopy (μFTIR) is a very useful method to analyze the biochemical properties of biological samples in situ. Many diseases affecting the central nervous system (CNS) have been studied using this method, to elucidate alterations in lipid oxidation or protein aggregation, among others. In this work, we describe in detail the characteristics between grey matter (GM) and white matter (WM) areas of the human brain by μFTIR, and we compare them with the mouse brain (strain C57BL/6), the most used animal model in neurological disorders. Our results show a clear different infrared profile between brain areas in the lipid region of both species. After applying a second derivative in the data, we established a 1.5 threshold value for the lipid/protein ratio to discriminate between GM and WM areas in non-pathological conditions. Furthermore, we demonstrated intrinsic differences of lipids and proteins by cerebral area. Lipids from GM present higher C=CH, C=O and CH3 functional groups compared to WM in humans and mice. Regarding proteins, GM present lower Amide II amounts and higher intramolecular β-sheet structure amounts with respect to WM in both species. However, the presence of intermolecular β-sheet structures, which is related to β-aggregation, was only observed in the GM of some human individuals. The present study defines the relevant biochemical properties of non-pathological human and mouse brains by μFTIR as a benchmark for future studies involving CNS pathological samples.

1. Introduction

Fourier Transform Infrared microspectroscopy (μFTIR) is an in situ method that offers the possibility of analyzing in detail the biochemical properties of fixed and non-fixed cells and tissues with no needed of staining, homogenization or further manipulations that can alter the nature of the samples. Lipid peroxidation, de/phosphorylation and protein conformations are some of the measurements that the infrared spectrum can provide with a high spatial resolution and sensitivity at the microscopic level [1,2,3].
The intrinsic chemical characteristics of a biological sample is essential information to understand pathophysiological processes. Among others, numerous studies have benefited from this method to elucidate biochemical and molecular features of different pathologies that affect the central nervous system (CNS). The brain is composed of two main areas clearly differentiated: grey matter (GM) and white matter (WM). GM is characterized by the presence of neuronal bodies, whereas WM is characterized by the absence of neuronal bodies and a high presence of myelinated axons. Depending on the disease, the affected brain area may be different and, therefore, biomedical studies are usually focused in one specific anatomical region. μFTIR analyses have been performed in different cerebral areas affected by neurodegenerative disorders, such as Alzheimer’s [4,5,6,7], Parkinson’s [8] and Huntington’s [9] diseases, to detect protein aggregation and lipid peroxidation. Infrared radiation it is also very useful to evaluate demyelinating diseases, such as multiple sclerosis, quantifying the amount and the status of lipids in WM areas. A decrease in the lipid content, together with a high lipid oxidation, is a characteristic feature observed by μFTIR in demyelinated areas [10,11,12]. In addition, some authors have studied by infrared spectroscopy the effect of CNS related pathologies on GM and WM, at the same time finding different alterations by area [9,11,13,14,15]. However, to properly understand the molecular mechanism of pathologies, the characterization of non-pathological tissues needs to be established first. Furthermore, it is important to highlight that most of the studies performed in the CNS are carried out in rodents and, in consequence, it is crucial to determine the differences and the similarities between humans and rodents with the aim of evaluating the reproducibility and applicability of this method, comparing animal models to human.
In the present study, we characterize the protein and lipid composition and their infrared spectral properties in GM and WM areas by μFTIR. Moreover, this work compares the regional-related properties observed by this method in the human brain with those of the mouse, which is the most widely used experimental model in the neuroscience field.

2. Materials and Methods

2.1. Human Samples

Postmortem human brain tissue was obtained from the Neurological Tissue Bank at the Biobanc-Hospital Clínic-IDIBAPS (Barcelona, Spain). The whole procedure was performed in accordance with the Helsinki Declaration, the Convention of the Council of Europe on Human Rights and Biomedicine and approved by the Ethical Committee of the Autonomous University of Barcelona and the Ethical Committee of Clinical Research-Hospital Clinic de Barcelona (A1-OF15016, 05/20/2015). Frontal cortex samples from individuals (n = 7) with no clinical neurological manifestations were used in this study (five women and two men; 78 ± 7.3 years old; 15.2 ± 6.6 h of postmortem delay).

2.2. Animal Samples

C57BL/6 mice (n = 5) were used in the present study (four females and one male). In order to mimic the experimental conditions between species, the age of the mice was 20–22 months old. Animals were maintained in conventional plastic cages with food and water ad libitum, in a 12 h light/dark cycle, at 22 ± 2 °C and 50–60% humidity. All experimental animal work was conducted according to Spanish regulations (Ley 32/2007, Real Decreto 1201/2005, Ley 9/2003 and Real Decreto 178/2004) in agreement with European Union directives (86/609/CEE, 91/628/CEE) and was approved by the Ethical Committee of the Autonomous University of Barcelona.

2.3. Sample Preparation

Human brains were frozen at the Neurological Tissue Bank and stored at −80 °C until use. For this study, cerebral frontal cortex samples containing grey matter (GM) and white matter (WM) areas were used. Mice were euthanized under an anesthesia solution of xylazine (30 mg/kg) and ketamine (120 mg/kg) and intracardially perfused for 10 min with 4% paraformaldehyde in 0.1 M phosphate buffer (pH 7.4). The brains were immediately post-fixed in the same solution for 4 h at 4 °C, cryoprotected with 30% sucrose solution in 0.1 M phosphate buffer for 48 h at 4 °C, frozen in ice-cold 2-methylbutane (320404, Sigma-Aldrich, St. Louis, MO, USA) and stored at −80 °C until use. The areas of study in mice were the cerebral cortex and the corpus callosum (bregma between 0.86 mm and −1.22 mm coordinates) as representative areas of GM and WM, respectively. Prior to μFTIR analysis, frozen sections from human and mouse brains containing the areas of interest were cut 8 µm thick on a cryostat (CM3050S Leica) and mounted onto polished calcium fluoride (CaF2) optical windows (CAFP20-1, Crystran, UK). To minimize water contribution, sections on CaF2 slides were air-dried at room temperature and stored in a vacuum protected from light until use.

2.4. μFTIR Data Acquisition

Fourier Transform Infrared microspectroscopy (μFTIR) based on synchrotron radiation was carried out at the MIRAS beamline of ALBA synchrotron light source (Catalonia, Spain) [16]. A Hyperion 3000 microscope equipped with a 36× magnification objective and coupled to a Vertex 70 spectrometer (Bruker, Billerica, MA, USA) was used. The spectra collection was performed in transmission mode at 4 cm−1 spectral resolution, 10 μm × 10 μm aperture dimensions and 128 scans. All spectra were obtained by means of Opus 7.5 software (Bruker, Billerica, MA, USA). The measuring range was 4000–600 cm−1 wavenumbers, and zero filling was performed with fast Fourier transform (FFT), so that we obtained one point every 2 cm−1 in the final spectra. Background spectra were collected from a clean area of each CaF2 window every 10 min. A mercury cadmium telluride (MCT) detector was used, and both the microscope and spectrometer were continuously purged with a flow of dried air. For each studied area, approximately 100 spectra with a step size of 50 μm × 50 μm for human samples and 30 μm × 30 μm for mouse samples were acquired. In order to represent regional differences with high spatial resolution in the tissue, maps of 200 spectra with a step size of 6 μm × 6 μm were performed in one representative sample of each species containing GM and WM areas.

2.5. μFTIR Spectra Analysis

Principal Component Analysis (PCA) was done using the open source software package Orange (Bioinformatics Laboratory of the University of Ljubljana, Version 3.23, Ljubljana, Slovenia) with the spectroscopy add-on (Version 0.4.9) [17]. The analysis was focused on the lipid region (3050–2800 cm−1) and the protein region (1800–1500 cm−1) of the spectra. The PCA were performed after the second derivatives were calculated and unit vector normalizations were performed on each spectrum for each region. The second derivatives were calculated with a polynomial order of 2 and a window size of 5 points for the lipid region and 17 points for the protein region. For a statistical and detailed study of lipid and protein composition, Unscrambler X software (CAMO Software, Oslo, Norway) was used. In this software, second derivation using the Savitsky–Golay algorithm with an eleven-point smoothing filter and a polynomial order of 3 was applied on the spectra to eliminate the baseline contribution and enhance narrow bands. After that, we measured the absorptions at the following wavenumbers related to functional groups of biochemical interest: ~2921 cm−1 (CH2 asymmetric stretching vibrations), ~2962 cm−1 (CH3 asymmetric stretching vibrations), ~3012 cm−1 (C=CH, unsaturated olefinic group), ~1743 cm−1 (C=O, carbonyl group), ~1656 cm−1 (α-helix protein structure), ~1637 cm−1 (intramolecular β-sheet protein structure), ~1625 cm−1 (intermolecular β-sheet protein structure) and ~1548 cm−1 (Amide II). The absorbances at the mentioned wavenumbers, 2 cm−1 before and at 2 cm−1 after, were obtained and averaged to represent the peaks of interest. CH2 and CH3, especially CH2, are functional groups mostly present in lipids, while α-helix and β-sheet structures belong to the Amide I functional group band present in proteins. To account for variations in the tissue thickness, CH3, C=CH and C=O peaks were normalized by the CH2 band (as an approximation of the total lipid content), whereas Amide II and β-sheet peaks were normalized by the total protein content (Amide I or α-helix bands).

2.6. Statistical Analysis

Statistical analysis and graphical representation were performed using the Graph Pad Prism® software. To determine differences between GM and WM, Student’s t-test was carried out on human and mouse data. For statistically significance, we considered p-value < 0.05. All data are expressed as mean values ± standard error of the mean (SEM).

3. Results

In this section, we present the mean spectra in the range 3500–1250 cm−1 wavenumbers and the PCA of the lipid and protein regions of GM and WM samples belonging to humans and mice. Moreover, statistical analyses were applied for the study of CH3 (d2A2962/d2A2921), C=CH (d2A3012/d2A2921) and C=O (d2A1743/d2A2921) compositions with respect to the lipid content, and Amide II (d2A1548/d2A1656+1637), intramolecular β-sheet structure (d2A1637/d2A1656) and intermolecular β-sheet structure (d2A1625/d2A1656+1637) with respect to the protein content, comparing their infrared absorptions between GM and WM in human and murine brain samples.

3.1. Different Chemical Profile of the Brain Tissue

In order to visualize the infrared profile of GM and WM areas, all spectra from each cerebral region were averaged in both species. μFTIR data show a clear different spectrum profile for each brain area characterized by higher CH2 (~2920 cm−1 and ~2850 cm−1) and CH3 (~2960 cm−1) peaks in WM with respect to the GM area; however, Amide I (~1650 cm−1) and Amide II (~1550 cm−1) peaks are very similar in both areas, as shown in Figure 1A,B. CH2 and CH3 groups are mainly related to lipids, which are very abundant in myelin, while Amide I and II groups are related to proteins. Our results demonstrate that lipid/protein (CH2/Amide I) amount is statistically higher in WM compared with GM, as shown in Figure 2A,B. This different chemical composition is also observed by high resolution imaging tissue maps, where GM and WM from the same sample were captured, as shown in Figure 2C,D. All the mentioned observations are present in humans and mice.

3.2. Principal Component Analysis of Brain Areas and Species

Principal Component Analysis (PCA) was applied to calculate Euclidean distances between each spectral dataset and to cluster the spectral properties of the four sample groups (WM and GM for humans and mice). The PCA scores plots directly unravel differences and similarities among the spectra, while the PCA loadings detail the spectral differences between the clusters. The analysis is shown in Figure 3—left panels show the comparison of the C-H bond region (3050–2800 cm−1). This region is mainly dominated by the infrared absorption of lipids in the system. Figure 3A shows the score plot and Figure 3B shows the corresponding PC1 (59%) and PC2 (17%) loadings. The score plot distribution shows that PC1 clearly differentiates human and mouse data, while PC2 differentiates between GM and WM. PC1 loading indicates that slight shifts characterize the differences between human and mouse samples, while PC2 indicates that differences between 3012 cm−1 (C=CH bond), 2921 cm−1 (CH2 group) and 2962 cm−1 (CH3 group) intensities differentiate between WM and GM. For a better illustration of differences, the average spectra of the second derivative of each group are plotted in Figure 3C. Right panels in Figure 3 show the PCA corresponding to the 1800–1500 cm−1 region where carbonyl groups and the Amide I and II absorb. The score plot in Figure 3D shows that, as in the C–H bond region, PC1 (54%) discriminates between humans and mice and PC2 (17%) between GM and WM. The loadings in Figure 3E indicate the wavelengths that allow for distinguishing between groups. The strongest contribution arises at 1625 cm−1 for PC1 and at 1637 cm−1 for PC2, indicating more intermolecular β-sheet structures for human samples, especially in GM, and more intramolecular β-sheet structures in human and mouse GM samples. Moreover, the contribution at 1685 cm−1 of PC1 suggests differences in the antiparallel β-sheet structures for humans and mice. As also observed in the representation of the average spectra of the second derivative of each group, as shown in Figure 3F, small changes in protein secondary structure discriminate between human and mouse spectra. Differently, changes in the intensity of the carbonyl band with a maximum at 1741 cm−1 discriminate between WM and GM.

3.3. Protein and Lipid Properties Differ by Cerebral Region

To thoroughly analyze the typical properties of lipids and proteins in brain areas, ratios for different infrared wavenumber absorptions after the second derivative were performed. For lipids, we quantified carbonyl (C=O), unsaturated olefinic (C=CH) and methyl (CH3) group absorbances with respect to the total lipidic amount in the sample, normalizing by the CH2 absorbance. Our results show that lipids in the GM area of humans and mice contain more carbonyl (d2A1743/d2A2921), olefinic (d2A3012/d2A2921) and methyl (d2A2962/d2A2921) groups with respect to the WM area, as shown in Figure 4A–C. Regarding proteins, we quantified Amide II, intramolecular β-sheet and intramolecular β-sheet absorbances with respect to the total protein amount normalizing by the Amide I absorbance. Compared to GM, our data show that proteins of WM present more absorption in the Amide II groups (d2A1548/d2A1656+1637), that account for the N–H and C–H absorption, and less intramolecular β-sheet structures (d2A1637/d2A1656) in humans and mice, as shown in Figure 4D,E. Interestingly, intermolecular β-sheet structures (d2A1625/d2A1656+1637) seem not to be present in mice, whereas in humans we can observe their presence in the GM area belonging to four of the seven subjects studied, as shown in Figure 4F. Attending to the mean values of the samples, this work also provides evidence of higher heterogeneity, especially in the C–H region, in GM than in WM, as shown in Figure 4.

4. Discussion

The present work shows that synchrotron-based μFTIR is a very sensitive and reproducible method for biochemical analyses of GM and WM areas of both human and mouse brain tissues. Furthermore, the molecular features reported by μFTIR based on synchrotron light can be simultaneously compared with histological features observed in the same examined sample. Depending on the tissue sample processing (e.g., paraffin embedding, freezing, chemical fixation, etc.) and the data analysis (e.g., spectral correction, normalization, derivation, etc.), the infrared values obtained may vary. Although Mazur and collaborators showed no significant differences in the infrared output of unfixed versus formalin-fixed cells [18,19], Hackett and collaborators demonstrated some slight differences in the infrared spectrum of the cryo-fixed tissue before and after formalin fixation [20]. However, our data in both cryo-fixed and paraformaldehyde-fixed tissue show a specific composition between brain areas already reported in rats [9,15,21] and humans [11], regardless of histological sample processing. A special note of caution is due for the analysis of specific wavenumbers, such as phosphate and sugars, depending on the chemistry of sample processing (phosphate buffer, sugar-based cryoprotectants, etc.). Because of differential sample processing between human and murine samples in this study, we did not analyze the phosphate region to avoid the substantial contribution of the buffer in murine samples compared to human samples.
Our infrared spectra confirm the different biochemical compositions between brain areas characterized by a higher amount of lipid content in WM compared to GM, showing a significantly higher lipid/protein (CH2/Amide I) ratio in WM with respect to GM. We can considerer the lipid/protein ratio to distinguish between brain areas by the μFTIR method, considering that second derivative absorption values < 1.5 correspond with measures in GM, whereas values > 1.5 correspond with measures in WM. Despite the different sample processing used between species, this discrimination can be perfectly applied to human and mouse samples. Thus, this lipid/protein 1.5 threshold value is an interesting GM versus WM discrimination factor, that may allow unsupervised μFTIR large data acquisition in human or murine whole brain samples. As a potential diagnostic method, values greater than 1.5 in GM or values lower than 1.5 in WM could be indicative of lipid and/or protein alterations in the brains of mice and humans. A similar lipid/protein threshold value (1.05) between GM and WM was established by Bonda and collaborators in rat brain cryo-fixed samples [9], supporting the use of the μFTIR method to discriminate between brain regions in humans, mice and rats. However, our data in cryo-fixed human samples show values that are slightly different to those reported in the study of Bonda et al. [9], but very similar to those of paraformaldehyde-fixed murine samples, which may, albeit unlikely, mean that mouse is more similar to human than rat as an animal model. A more likely explanation is that differential sample processing (fixation procedure) and data processing (water subtraction, calculation of area, derivation, etc.) can account for variability between studies.
In addition, we are able to clearly visualize WM and GM in the tissue by high resolution infrared maps based on lipid/protein ratio with no need of other histological or biochemical procedures. Among the many biological differences between GM and WM, one of the most significant is that GM is predominantly composed of rich-protein neuronal somas and WM by myelinated neuronal axons [22]. Since WM is mainly constituted by multilamellar myelin sheaths and the lipid composition of myelin is ~70% dry weight [23], the higher content of lipid in WM can be attributed to this multilamellar structure surrounding the axons. In addition to neurons and myelin, we have to considerer that both GM and WM are composed by glial cells [24,25]. Cellular density and intrinsic features of glial cells differ in each brain area and, therefore, they probably also participate in the specific tissue chemical composition described in this study.
The PCA allows us to identify spectral differences and similarities between the four sample groups. On the one hand, we observed that humans and mice share common spectral features (3012 cm−1 and 2962 cm−1), while the distinction can be made between GM and WM. On the other hand, strong differences for the PC1 of the protein region (1741 cm−1 and 1625 cm−1) were found between humans and mice.
Our results show that lipid characteristics vary between GM and WM in humans and mice. In both species, CH3/CH2, C=CH/CH2 and C=O/CH2 ratios were higher in GM compared to WM. Since lipids are composed mostly by long CH2 chains and to a lesser extent by methyl groups, CH2 groups are used to measure the lipid content of a sample, while CH3 groups could be useful to determine the lipid chain length. Cholesterol and cerebrosides are the most abundant lipids in myelin, constituting approximately 27% and 19%, respectively, of the lipids in human WM [21,26]. Whereas cholesterol is also highly present in GM, high concentration of cerebroside is specific for WM but not for GM. As an example, galactocerebrosides, the main cerebrosides present in the WM, are composed of two chains of numerous CH2 groups and only two CH3 groups at the end of the chains, arguing for a differential observation depending on lipid tail length in the ratio CH3/CH2. However, as it is shown in Benseny-Cases et al. [27], proteins have a similar number of CH2 and CH3 groups in their lateral chains, whereas lipids have significantly more CH2 than CH3 in their hydrophobic tail. Thus, changes in the lipid/protein ratio also implies changes in the C–H region (3000–2000 cm−1). Absorptions of olefinic C=CH functional group at 3012 cm−1 and C=O functional group at 1743 cm−1 were used to evaluate lipid unsaturation and lipid carbonyl content, respectively. C=CH and C=O bands examined by μFTIR are absent in cholesterol, galactocerebroside, sphingomyelin and sulfatide lipids, which are characteristic of WM, and present in some phospholipids [28]. Thus, our results show higher C=CH/CH2 and C=O/CH2 ratios in GM with respect to WM and are in agreement with these data. However, Krafft et al. showed a C=CH peak in galactocerebroside, sphingomyelin and sulfatide lipids by Raman spectroscopy [29]. On the other hand, it is known that free radicals mainly react with carbon double bond sites present in unsaturated lipids, resulting in carbonyl compound formation and lipid chain fragmentation [30,31,32,33]. Therefore, the high expression of C=O groups is related to the stable end-products of lipid peroxidation and, in consequence, the analyses of lipid oxidation status, measuring C=O/CH2 content by μFTIR, have been performed in different studies [15,27,34]. Taking these data into account, our observations suggest that the lipids in GM could be more oxidized than in WM. Since carbon double bonds are the main target of free radicals, lipid oxidation should result in less olefinic content, however, we found a higher C=CH/CH2 ratio in GM than WM. Some authors have attributed this high C=CH absorption to the accumulation of lipid peroxidation end-products, which contain numerous double bonds [4,35,36]. Altogether, the lower lipid content associated with the higher olefinic and carbonyl content in GM, could be interpreted as a higher oxidative status of GM with respect to WM. Brain tissue, especially WM, is very vulnerable to oxidative stress due to the presence of a rich lipid content, high oxidative metabolic activity and low antioxidant mechanisms [37]. However, in non-pathological samples, we observed more markers related to oxidation in GM than WM. The possible accumulation of oxidized lipids in GM could be due to the greater vulnerability to oxidative stress that neurons present [38].
Regarding protein composition, we also found differences between GM and WM areas. Analyzing the ratio Amide II/Amide I, we observed higher Amide II absorption (1548 cm−1) in human and mouse WM with respect to GM. Mathematical methods, such as deconvolution or the second derivation of the Amide I band, allow us to detect α-helix at ~1656 cm−1 and intramolecular β-sheet at ~1637 cm−1 protein secondary conformations by μFTIR [1,39,40]. Our Amide I analysis reveals that WM is composed by fewer intramolecular β-sheet structures (d2A1637/d2A1656) compared with GM. These results were observed in both human and mouse brains. In addition to the principal β-sheet structure at ~1637 cm−1, more β-components have been detected in the infrared spectra after resolution enhancement. The presence of a band at ~1625 cm−1 represents an unusual intermolecular β-sheet structure [41] that has been observed in β-amyloid plaques [4,42,43,44], α-synuclein [45] and prion proteins [46]. In consequence, this band is used as a marker of protein β-aggregation in several pathologies [5,6,9,45,46,47]. Infrared absorption intensity of the ratio d2A1625/d2A1656+1637 after second derivative shows negative values in human WM and in both mouse brain areas. This low vibrational intensity at 1625 cm−1 could be due to no expression of aggregated proteins in these areas. Interestingly, some humans show positive values at this wavelength in GM. Considering studied samples belong to aged humans, the apparition of this peak could be due to the age-dependent onset of protein aggregation [48,49,50]. Murine samples also belong to aged mice and there is no expression of this band. This finding may indicate that in mice there is no production of β-aggregation in non-pathological conditions, while humans are susceptible to produce it in GM. Furthermore, a weak band at ~1685 cm−1, which is assigned to antiparallel β-sheet structures [41,51], was also observed in the GM of humans. The antiparallel β-sheet structure has been described as a signature of β-amyloid oligomers being absent in fibrils [42,52]. Thereby, infrared absorptions at both 1685 cm−1 and 1625 cm−1 in the GM of some individuals could suggest the presence of oligomeric β-amyloid aggregates. This assumption is consistent with the fact that only genetically modified mice are capable to accumulate β-amyloid [53], and on the contrary, has been reported that healthy elder humans can present β-amyloid plaques in GM areas [48,49,50]. Furthermore, the β-aggregation, observed specifically in the GM of some individuals, is in agreement with the weaker mechanisms that present neurons to avoid protein aggregation with respect to other cells types of the CNS [54,55]. Nevertheless, β-aggregation in healthy aged individuals and in pathophysiology, both in human and murine samples, raises an issue that requires further research.
In summary, this work demonstrates by μFTIR: (i) a specific lipid/protein composition of GM and WM areas; (ii) a different presence of lipid and protein functional groups by cerebral area; (iii) that the regional differences observed in the human brain are fairly reproduced in mice, with the exception of β-aggregation content. These similarities between species in non-pathological conditions make the mouse a good experimental animal model. µFTIR comparative biology studies, but also comparative sample and data processing studies [18,19,20] are essential to assess and understand differences and similarities between species, especially for disease-related animal models. Specifically, the study of lipid-related disorders in mice could be very reliable in humans due to the high reproducibility observed in the lipid characteristics of GM and WM; however, the study of protein-related pathologies in mouse models requires a proper assessment sample and data processing variables.

5. Conclusions

The understanding of the brain biochemical properties in non-pathological conditions shown in this work will help us to study the molecular dynamics in different pathologies affecting GM, WM or both. Interestingly, this study provides evidence of mice (strain C57BL/6) as good animal models to reproduce human properties of the brain. Thus, μFTIR is a powerful tool to progress in the knowledge of CNS physiology and pathology performed in mice, by which data can be extrapolated to humans.

Author Contributions

Conceptualization, all authors; investigation, P.S.-M., M.K. and N.B.-C.; resources, A.P.-M., B.C. and T.V.; writing—original draft preparation, P.S.-M., M.K., N.B.-C. and A.P.-M.; writing—review and editing, all authors; funding acquisition, A.P.-M. and B.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Spanish Ministry of Economy and Business through Plan Estatal BFU2017-87843-R to A.P.-M. and B.C.

Acknowledgments

These experiments were performed at MIRAS beamline at ALBA Synchrotron with the collaboration of ALBA staff. We would like to thank Ellen Gelpí from the Neurological Tissue Bank of the Biobanc (Barcelona) for providing the human samples.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

References

  1. Kong, J.; Yu, S. Fourier transform infrared spectroscopic analysis of protein secondary structures. Acta Biochim. Biophys. Sin. (Shanghai) 2007, 39, 549–559. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Lamba, O.P.; Lal, S.; Yappert, M.C.; Lou, M.F.; Borchman, D. Spectroscopic detection of lipid peroxidation products and structural changes in a sphingomyelin model system. Biochim. Biophys. Acta 1991, 1081, 181–187. [Google Scholar] [CrossRef]
  3. Sarver, R.W., Jr.; Krueger, W.C. Protein secondary structure from Fourier transform infrared spectroscopy: A data base analysis. Anal. Biochem. 1991, 194, 89–100. [Google Scholar] [CrossRef]
  4. Benseny-Cases, N.; Klementieva, O.; Cotte, M.; Ferrer, I.; Cladera, J. Microspectroscopy (muFTIR) reveals co-localization of lipid oxidation and amyloid plaques in human Alzheimer disease brains. Anal. Chem. 2014, 86, 12047–12054. [Google Scholar] [CrossRef] [PubMed]
  5. Benseny-Cases, N.; Alvarez-Marimon, E.; Aso, E.; Carmona, M.; Klementieva, O.; Appelhans, D.; Ferrer, I.; Cladera, J. In situ structural characterization of early amyloid aggregates in Alzheimer’s disease transgenic mice and Octodon degus. Sci. Rep. 2020, 10, 5888. [Google Scholar] [CrossRef] [Green Version]
  6. Klementieva, O.; Willen, K.; Martinsson, I.; Israelsson, B.; Engdahl, A.; Cladera, J.; Uvdal, P.; Gouras, G.K. Pre-plaque conformational changes in Alzheimer’s disease-linked Abeta and APP. Nat. Commun. 2017, 8, 14726. [Google Scholar] [CrossRef] [PubMed]
  7. Liao, C.R.; Rak, M.; Lund, J.; Unger, M.; Platt, E.; Albensi, B.C.; Hirschmugl, C.J.; Gough, K.M. Synchrotron FTIR reveals lipid around and within amyloid plaques in transgenic mice and Alzheimer’s disease brain. Analyst 2013, 138, 3991–3997. [Google Scholar] [CrossRef] [Green Version]
  8. Szczerbowska-Boruchowska, M.; Dumas, P.; Kastyak, M.Z.; Chwiej, J.; Lankosz, M.; Adamek, D.; Krygowska-Wajs, A. Biomolecular investigation of human substantia nigra in Parkinson’s disease by synchrotron radiation Fourier transform infrared microspectroscopy. Arch. Biochem. Biophys. 2007, 459, 241–248. [Google Scholar] [CrossRef]
  9. Bonda, M.; Perrin, V.; Vileno, B.; Runne, H.; Kretlow, A.; Forro, L.; Luthi-Carter, R.; Miller, L.M.; Jeney, S. Synchrotron infrared microspectroscopy detecting the evolution of Huntington’s disease neuropathology and suggesting unique correlates of dysfunction in white versus gray brain matter. Anal. Chem. 2011, 83, 7712–7720. [Google Scholar] [CrossRef] [Green Version]
  10. Heraud, P.; Caine, S.; Campanale, N.; Karnezis, T.; McNaughton, D.; Wood, B.R.; Tobin, M.J.; Bernard, C.C. Early detection of the chemical changes occurring during the induction and prevention of autoimmune-mediated demyelination detected by FT-IR imaging. Neuroimage 2010, 49, 1180–1189. [Google Scholar] [CrossRef]
  11. Choo, L.P.; Jackson, M.; Halliday, W.C.; Mantsch, H.H. Infrared spectroscopic characterisation of multiple sclerosis plaques in the human central nervous system. Biochim. Biophys. Acta 1993, 1182, 333–337. [Google Scholar] [CrossRef]
  12. LeVine, S.M.; Wetzel, D.L. Chemical analysis of multiple sclerosis lesions by FT-IR microspectroscopy. Free Radic. Biol. Med. 1998, 25, 33–41. [Google Scholar] [CrossRef]
  13. Ali, M.H.M.; Rakib, F.; Abdelalim, E.M.; Limbeck, A.; Mall, R.; Ullah, E.; Mesaeli, N.; McNaughton, D.; Ahmed, T.; Al-Saad, K. Fourier-Transform Infrared Imaging Spectroscopy and Laser Ablation -ICPMS New Vistas for Biochemical Analyses of Ischemic Stroke in Rat Brain. Front. Neurosci. 2018, 12, 647. [Google Scholar] [CrossRef] [PubMed]
  14. Leskovjan, A.C.; Kretlow, A.; Miller, L.M. Fourier transform infrared imaging showing reduced unsaturated lipid content in the hippocampus of a mouse model of Alzheimer’s disease. Anal. Chem. 2010, 82, 2711–2716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Cakmak, G.; Miller, L.M.; Zorlu, F.; Severcan, F. Amifostine, a radioprotectant agent, protects rat brain tissue lipids against ionizing radiation induced damage: An FTIR microspectroscopic imaging study. Arch. Biochem. Biophys. 2012, 520, 67–73. [Google Scholar] [CrossRef]
  16. Yousef, I.; Ribó, L.; Crisol, A.; Šics, I.; Ellis, G.; Ducic, T.; Kreuzer, M.; Benseny-Cases, N.; Quispe, M.; Dumas, P.; et al. MIRAS: The Infrared Synchrotron Radiation Beamline at ALBA. Synchrotron Radiat. News 2017, 30, 4–6. [Google Scholar] [CrossRef]
  17. Demšar, J.; Curk, T.; Erjavec, A.; Gorup, Č.; Hočevar, T.; Milutinovič, M.; Možina, M.; Polajnar, M.; Toplak, M.; Starič, A.; et al. Orange: Data mining toolbox in python. J. Mach. Learn. Res. 2013, 14, 2349–2353. [Google Scholar]
  18. Mazur, A.I.; Marcsisin, E.J.; Bird, B.; Miljković, M.; Diem, M. Evaluating different fixation protocols for spectral cytopathology, part 1. Anal. Chem. 2012, 84, 1259–1266. [Google Scholar] [CrossRef] [Green Version]
  19. Mazur, A.I.; Marcsisin, E.J.; Bird, B.; Miljković, M.; Diem, M. Evaluating different fixation protocols for spectral cytopathology, part 2: Cultured cells. Anal. Chem. 2012, 84, 8265–8271. [Google Scholar] [CrossRef] [Green Version]
  20. Hackett, M.J.; McQuillan, J.A.; El-Assaad, F.; Aitken, J.B.; Levina, A.; Cohen, D.D.; Siegele, R.; Carter, E.A.; Grau, G.E.; Hunt, N.H.; et al. Chemical alterations to murine brain tissue induced by formalin fixation: Implications for biospestroscopic imaging and mapping studies of disease pathogenesis. Analyst 2011, 136, 2941–2952. [Google Scholar] [CrossRef]
  21. Levine, S.M.; Wetzel, D.L.B. Analysis of Brain Tissue by FT-IR Microspectroscopy. Appl. Spectrosc. Rev. 1993, 28, 385–412. [Google Scholar] [CrossRef]
  22. Purves, D.A.; Augustine, G.J.; Fitzpatrick, D.; Hall, W.C.; LaMantia, A.S.; McNamara, J.O.; White, L.E. Neuroscience, 3th ed.; Sinauer Associates: Sunderland, MA, USA, 2004. [Google Scholar]
  23. Norton, W.T.; Poduslo, S.E. Myelination in rat brain: Changes in myelin composition during brain maturation. J. Neurochem. 1973, 21, 759–773. [Google Scholar] [CrossRef] [PubMed]
  24. von Bartheld, C.S.; Bahney, J.; Herculano-Houzel, S. The search for true numbers of neurons and glial cells in the human brain: A review of 150 years of cell counting. J. Comp. Neurol. 2016, 524, 3865–3895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Herculano-Houzel, S. The glia/neuron ratio: How it varies uniformly across brain structures and species and what that means for brain physiology and evolution. Glia 2014, 62, 1377–1391. [Google Scholar] [CrossRef] [PubMed]
  26. Olsson, N.U.; Harding, A.J.; Harper, C.; Salem, N., Jr. High-performance liquid chromatography method with light-scattering detection for measurements of lipid class composition: Analysis of brains from alcoholics. J. Chromatogr. B Biomed. Appl. 1996, 681, 213–218. [Google Scholar] [CrossRef]
  27. Benseny-Cases, N.; Alvarez-Marimon, E.; Castillo-Michel, H.; Cotte, M.; Falcon, C.; Cladera, J. Synchrotron-Based Fourier Transform Infrared Microspectroscopy (muFTIR) Study on the Effect of Alzheimer’s Abeta Amorphous and Fibrillar Aggregates on PC12 Cells. Anal. Chem. 2018, 90, 2772–2779. [Google Scholar] [CrossRef]
  28. Dreissig, I.; Machill, S.; Salzer, R.; Krafft, C. Quantification of brain lipids by FTIR spectroscopy and partial least squares regression. Spectrochim. Acta A Mol. Biomo.l Spectrosc. 2009, 71, 2069–2075. [Google Scholar] [CrossRef]
  29. Krafft, C.; Neudert, L.; Simat, T.; Salzer, R. Near infrared Raman spectra of human brain lipids. Spectrochim. Acta A Mol. Biomo.l Spectrosc. 2005, 61, 1529–1535. [Google Scholar] [CrossRef]
  30. Rice-Evans, C.A. Formation of free radicals and mechanisms of action in normal biochemical processes and pathological states. In Free Radical Damage and its Control; Rice-Evans, C.A., Burdon, R.H., Eds.; Elsevier Science: Amsterdam, The Netherlands, 1994; Volume 28, pp. 131–153. [Google Scholar]
  31. de Zwart, L.L.; Meerman, J.H.; Commandeur, J.N.; Vermeulen, N.P. Biomarkers of free radical damage applications in experimental animals and in humans. Free Radic. Biol. Med. 1999, 26, 202–226. [Google Scholar] [CrossRef]
  32. Niki, E.; Yoshida, Y.; Saito, Y.; Noguchi, N. Lipid peroxidation: Mechanisms, inhibition, and biological effects. Biochem. Biophys. Res. Commun. 2005, 338, 668–676. [Google Scholar] [CrossRef]
  33. Yin, H.; Xu, L.; Porter, N.A. Free radical lipid peroxidation: Mechanisms and analysis. Chem. Rev. 2011, 111, 5944–5972. [Google Scholar] [CrossRef] [PubMed]
  34. Rakib, F.; Ali, C.M.; Yousuf, M.; Afifi, M.; Bhatt, P.R.; Ullah, E.; Al-Saad, K.; Ali, M.H.M. Investigation of Biochemical Alterations in Ischemic Stroke Using Fourier Transform Infrared Imaging Spectroscopy-A Preliminary Study. Brain Sci. 2019, 9, 293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Severcan, F.; Gorgulu, G.; Gorgulu, S.T.; Guray, T. Rapid monitoring of diabetes-induced lipid peroxidation by Fourier transform infrared spectroscopy: Evidence from rat liver microsomal membranes. Anal. Biochem. 2005, 339, 36–40. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, K.B.; Bose, R.; Mantsch, H.H. Infrared spectroscopic study of diabetic platelets. Vib. Spectrosc. 2002, 28, 131–136. [Google Scholar] [CrossRef]
  37. Halliwell, B. Reactive oxygen species and the central nervous system. J. Neurochem. 1992, 59, 1609–1623. [Google Scholar] [CrossRef] [PubMed]
  38. Lassmann, H.; van Horssen, J. Oxidative stress and its impact on neurons and glia in multiple sclerosis lesions. Biochim. Biophys. Acta 2016, 1862, 506–510. [Google Scholar] [CrossRef]
  39. Byler, D.M.; Susi, H. Examination of the secondary structure of proteins by deconvolved FTIR spectra. Biopolymers 1986, 25, 469–487. [Google Scholar] [CrossRef]
  40. Dong, A.; Huang, P.; Caughey, W.S. Protein secondary structures in water from second-derivative amide I infrared spectra. Biochemistry 1990, 29, 3303–3308. [Google Scholar] [CrossRef]
  41. Arrondo, J.L.; Young, N.M.; Mantsch, H.H. The solution structure of concanavalin A probed by FT-IR spectroscopy. Biochim. Biophys. Acta 1988, 952, 261–268. [Google Scholar] [CrossRef]
  42. Peralvarez-Marin, A.; Barth, A.; Graslund, A. Time-resolved infrared spectroscopy of pH-induced aggregation of the Alzheimer Abeta(1-28) peptide. J. Mol. Biol. 2008, 379, 589–596. [Google Scholar] [CrossRef]
  43. Benseny-Cases, N.; Cocera, M.; Cladera, J. Conversion of non-fibrillar beta-sheet oligomers into amyloid fibrils in Alzheimer’s disease amyloid peptide aggregation. Biochem. Biophys. Res. Commun. 2007, 361, 916–921. [Google Scholar] [CrossRef] [PubMed]
  44. Rak, M.; Del Bigio, M.R.; Mai, S.; Westaway, D.; Gough, K. Dense-core and diffuse Abeta plaques in TgCRND8 mice studied with synchrotron FTIR microspectroscopy. Biopolymers 2007, 87, 207–217. [Google Scholar] [CrossRef] [PubMed]
  45. Conway, K.A.; Harper, J.D.; Lansbury, P.T., Jr. Fibrils formed in vitro from alpha-synuclein and two mutant forms linked to Parkinson’s disease are typical amyloid. Biochemistry 2000, 39, 2552–2563. [Google Scholar] [CrossRef] [PubMed]
  46. Kretlow, A.; Wang, Q.; Beekes, M.; Naumann, D.; Miller, L.M. Changes in protein structure and distribution observed at pre-clinical stages of scrapie pathogenesis. Biochim. Biophys. Acta 2008, 1782, 559–565. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Miller, L.M.; Bourassa, M.W.; Smith, R.J. FTIR spectroscopic imaging of protein aggregation in living cells. Biochim. Biophys. Acta 2013, 1828, 2339–2346. [Google Scholar] [CrossRef] [Green Version]
  48. Price, J.L.; Morris, J.C. Tangles and plaques in nondemented aging and “preclinical” Alzheimer’s disease. Ann. Neurol. 1999, 45, 358–368. [Google Scholar] [CrossRef]
  49. Rodrigue, K.M.; Kennedy, K.M.; Devous, M.D., Sr.; Rieck, J.R.; Hebrank, A.C.; Diaz-Arrastia, R.; Mathews, D.; Park, D.C. beta-Amyloid burden in healthy aging: Regional distribution and cognitive consequences. Neurology 2012, 78, 387–395. [Google Scholar] [CrossRef] [Green Version]
  50. Spies, P.E.; Verbeek, M.M.; van Groen, T.; Claassen, J.A. Reviewing reasons for the decreased CSF Abeta42 concentration in Alzheimer disease. Front. Biosci. 2012, 17, 2024–2034. [Google Scholar] [CrossRef] [Green Version]
  51. Kubelka, J.; Keiderling, T.A. Differentation of Beta-Sheet Forming Structures: Ab Initio-Based Simulations of IR Absorption nd Vibrational CD for Model Peptide and Protein Beta-Sheets. J. Am. Chem. Soc. 2001, 123, 12048–12058. [Google Scholar] [CrossRef]
  52. Cerf, E.; Sarroukh, R.; Tamamizu-Kato, S.; Breydo, L.; Derclaye, S.; Dufrene, Y.F.; Narayanaswami, V.; Goormaghtigh, E.; Ruysschaert, J.M.; Raussens, V. Antiparallel beta-sheet: A signature structure of the oligomeric amyloid beta-peptide. Biochem. J. 2009, 421, 415–423. [Google Scholar] [CrossRef]
  53. Elder, G.A.; Gama Sosa, M.A.; De Gasperi, R. Transgenic mouse models of Alzheimer’s disease. Mt. Sinai. J. Med. 2010, 77, 69–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kundra, R.; Dobson, C.M.; Vendruscolo, M. A Cell- and Tissue-Specific Weakness of the Protein Homeostasis System Underlies Brain Vulnerability to Protein Aggregation. iScience 2020, 23, 100934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Fu, H.; Possenti, A.; Freer, R.; Nakano, Y.; Hernandez Villegas, N.C.; Tang, M.; Cauhy, P.V.M.; Lassus, B.A.; Chen, S.; Fowler, S.L.; et al. A tau homeostasis signature is linked with the cellular and regional vulnerability of excitatory neurons to tau pathology. Nat. Neurosci. 2019, 22, 47–56. [Google Scholar] [CrossRef] [PubMed]
Figure 1. μFTIR profile of grey matter (GM) and white matter (WM). Spectra of GM (blue line) and WM (red line) areas were collected and average for human (A) and mouse (B) brain samples. Note the increased absorptions of CH2 and CH3 peaks (*) in WM with respect to GM. a.u.—arbitrary units.
Figure 1. μFTIR profile of grey matter (GM) and white matter (WM). Spectra of GM (blue line) and WM (red line) areas were collected and average for human (A) and mouse (B) brain samples. Note the increased absorptions of CH2 and CH3 peaks (*) in WM with respect to GM. a.u.—arbitrary units.
Biomolecules 10 01099 g001
Figure 2. Lipid/Protein composition by μFTIR in grey matter (GM) and white matter (WM). Quantification of CH2/Amide I absorbance after second derivative (d2A2921/d2A1656) in humans (A) and mice (B) confirms higher lipid amount in WM with respect to GM. Each dot in the graphs correspond to the mean of approximately 100 infrared measurements from a subject. Representative in situ heat maps from human (C) and murine (D) cerebral tissue showing CH2/Amide I ratio. Dashed lines in the tissue images represent the limit between GM and WM. **** p-value < 0.0001. Scale bar = 50 μm.
Figure 2. Lipid/Protein composition by μFTIR in grey matter (GM) and white matter (WM). Quantification of CH2/Amide I absorbance after second derivative (d2A2921/d2A1656) in humans (A) and mice (B) confirms higher lipid amount in WM with respect to GM. Each dot in the graphs correspond to the mean of approximately 100 infrared measurements from a subject. Representative in situ heat maps from human (C) and murine (D) cerebral tissue showing CH2/Amide I ratio. Dashed lines in the tissue images represent the limit between GM and WM. **** p-value < 0.0001. Scale bar = 50 μm.
Biomolecules 10 01099 g002
Figure 3. Principal Component Analysis (PCA) of grey matter (GM) and white matter (WM). PCA of the normalized second derivatives of the lipid region (3050–2800 cm−1) and protein region (1800–1500 cm−1), showing the scores plots (A,D), the resulting loadings for PC-1 and PC-2 (B,E) and the average spectra (C,F) of the two regions, respectively. Spectra shadowing in C and F indicates standard deviation. a.u.—arbitrary units.
Figure 3. Principal Component Analysis (PCA) of grey matter (GM) and white matter (WM). PCA of the normalized second derivatives of the lipid region (3050–2800 cm−1) and protein region (1800–1500 cm−1), showing the scores plots (A,D), the resulting loadings for PC-1 and PC-2 (B,E) and the average spectra (C,F) of the two regions, respectively. Spectra shadowing in C and F indicates standard deviation. a.u.—arbitrary units.
Biomolecules 10 01099 g003
Figure 4. Chemical properties of lipids and proteins in grey matter (GM) and white matter (WM). Graphs represents absorbance second derivative (d2A) ratios of functional groups related to lipids (AC) and proteins (DF) in human and murine brain samples. Each dot in the graphs correspond to the mean of approximately 100 infrared measurements from a subject. * p-value < 0.05; ** p-value < 0.01; *** p-value < 0.001; **** p-value < 0.0001.
Figure 4. Chemical properties of lipids and proteins in grey matter (GM) and white matter (WM). Graphs represents absorbance second derivative (d2A) ratios of functional groups related to lipids (AC) and proteins (DF) in human and murine brain samples. Each dot in the graphs correspond to the mean of approximately 100 infrared measurements from a subject. * p-value < 0.05; ** p-value < 0.01; *** p-value < 0.001; **** p-value < 0.0001.
Biomolecules 10 01099 g004

Share and Cite

MDPI and ACS Style

Sanchez-Molina, P.; Kreuzer, M.; Benseny-Cases, N.; Valente, T.; Almolda, B.; González, B.; Castellano, B.; Perálvarez-Marín, A. From Mouse to Human: Comparative Analysis between Grey and White Matter by Synchrotron-Fourier Transformed Infrared Microspectroscopy. Biomolecules 2020, 10, 1099. https://doi.org/10.3390/biom10081099

AMA Style

Sanchez-Molina P, Kreuzer M, Benseny-Cases N, Valente T, Almolda B, González B, Castellano B, Perálvarez-Marín A. From Mouse to Human: Comparative Analysis between Grey and White Matter by Synchrotron-Fourier Transformed Infrared Microspectroscopy. Biomolecules. 2020; 10(8):1099. https://doi.org/10.3390/biom10081099

Chicago/Turabian Style

Sanchez-Molina, Paula, Martin Kreuzer, Núria Benseny-Cases, Tony Valente, Beatriz Almolda, Berta González, Bernardo Castellano, and Alex Perálvarez-Marín. 2020. "From Mouse to Human: Comparative Analysis between Grey and White Matter by Synchrotron-Fourier Transformed Infrared Microspectroscopy" Biomolecules 10, no. 8: 1099. https://doi.org/10.3390/biom10081099

APA Style

Sanchez-Molina, P., Kreuzer, M., Benseny-Cases, N., Valente, T., Almolda, B., González, B., Castellano, B., & Perálvarez-Marín, A. (2020). From Mouse to Human: Comparative Analysis between Grey and White Matter by Synchrotron-Fourier Transformed Infrared Microspectroscopy. Biomolecules, 10(8), 1099. https://doi.org/10.3390/biom10081099

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop