Biomarkers

Introduction – Why do we need ‘biomarkers? Biomarkers of protein oxidation Introduction Major issues/questions Protein carbonyl biomarkers Biochemistry Methods of measurement Storage, stability and limitations in use Protein thiol biomarkers Biochemistry Methods of measurement Storage, stability and limitations on use Aliphatic amino acid biomarkers Biochemistry Methods of measurement Storage, stability and limitations on use Oxidised Tryptophan Biomarkers Biochemistry Method of measurement Storage, stability and limitations on use Oxidised tyrosine biomarkers Biochemistry Methods of measurement Storage, stability and limitations on use Formation of neoepitopes on oxidised proteins Validation of assays for protein oxidation biomarkers Relationship of protein oxidation to disease Modulation of protein oxidation biomarkers by antioxidants Future perspectives Introduction to lipid peroxidation biomarkers Introduction: biochemistry of lipid peroxidation Malondialdehyde Methods of measurement Storage, stability and limitations on use Conjugated dienes Method of measurement Storage, stability and limitations of use LDL lag phase Method of measurement Storage, stability and limitations of use Hydrocarbon gases Biochemistry Method of measurement Storage, stability and limitations on use Lipofuscin Biochemistry Method of measurement Storage, stability and limitation on use Lipid peroxides Biochemistry Method of measurement Storage, stability and limitations on use Isoprostanes Biochemistry Method of measurement Storage, stability and limitations on use Possible new biomarkers of lipid oxidation Relationship of lipid peroxidation to disease Modulation of lipid peroxidation biomarkers by antioxidants Functional consequences of lipid peroxidation Contribution of dietary intake to lipid peroxidation products Biomarkers of DNA oxidation Introduction Confounding factors Units and terminology Nuclear and mitochondrial DNA damage Lymphocytes as surrogate tissues Measurement of DNA damage with the comet assay Practical details Storage, stability, and limitations of the assay Measurement of DNA base oxidation by HPLC Practical details Storage, stability and limitations of the method Measurement of DNA base oxidation by GC–MS Biochemistry of 8-oxoguanine, adenine and fapy derivatives Methods of measurement Storage, stability and limitations of the method Analysis of guanine oxidation products in urine Method of measurement Limitations and criticisms Immunochemical methods Methods of measurement Storage, stability, and limitations of the assay 32P post-labelling Method of measurement Limitations and criticisms Validation of assays for DNA oxidation Oxo-dGuo in lymphocyte DNA Urinary measurements DNA–aldehyde adducts Biochemistry Method of measurement Products of reactive nitrogen species Endpoints arising from oxidative DNA damage Mutations Chromosome aberrations Micronuclei Site-specific DNA damage Relationship of DNA oxidation to disease Modulation of DNA oxidation biomarkers by antioxidants Direct and indirect effects of oxidative stress: measures of total oxidant/antioxidant levels Visualisation of cellular oxidants Biochemistry: histochemical detection of ROS Method of measurement Limitations, storage and stability Measurement of hydrogen peroxide Biochemistry Methods of measurement Storage, stability and limitations on use Measurement of the ratio of antioxidant/oxidised antioxidant Biochemistry Method of measurement Storage, stability and limitations on use Total antioxidant capacity Biochemistry Terminology Methods of measurement Storage, stability and limitations on use Validation of assays for direct oxidant and antioxidant biomarkers Relationship of oxidant/antioxidant measurement to disease Modulation of oxidant/antioxidant biomarkers by dietary antioxidants Induction of genes in response to oxidative stress Background Measurement of antioxidant responsive genes and proteins Effects of antioxidant intake on the activity of antioxidant enzymes

Matrix-free mass spectrometric imaging using laser desorption ionisation Fourier transform ion cyclotron resonance mass spectrometry

Mass spectrometry imaging (MSI) is a powerful tool in metabolomics and proteomics for the spatial localization and identification of pharmaceuticals, metabolites, lipids, peptides and proteins in biological tissues. However, sample preparation remains a crucial variable in obtaining the most accurate distributions. Common washing steps used to remove salts, and solvent-based matrix application, allow analyte spreading to occur. Solvent-free matrix applications can reduce this risk, but increase the possibility of ionisation bias due to matrix adhesion to tissue sections. We report here the use of matrix-free MSI using laser desorption ionisation performed on a 12 T Fourier transform ion cyclotron resonance (FTICR) mass spectrometer. We used unprocessed tissue with no post-processing following thaw-mounting on matrix-assisted laser desorption ionisation (MALDI) indium-tin oxide (ITO) target plates. The identification and distribution of a range of phospholipids in mouse brain and kidney sections are presented and compared with previously published MALDI time-of-flight (TOF) MSI distributions.

Kinetic study on thermal decomposition of woods in oxidative environment

The purpose of this work is to gain knowledge on kinetics of biomass decomposition under oxidative atmospheres, mainly examining effect of heating rate on different biomass species. Two sets of experiments are carried out: the first set of experiments is thermal decomposition of four different wood particles, namely aspens, birch, oak and pine under an oxidative atmosphere and analysis with TGA; and the second set is to use large size samples of wood under different heat fluxes in a purpose-built furnace, where the temperature distribution, mass loss and ignition characteristics are recorded and analyzed by a data post-processing system. The experimental data is then used to develop a two-step reactions kinetic scheme with low and high temperature regions while the activation energy for the reactions of the species under different heating rates is calculated. It is found that the activation energy of the second stage reaction for the species with similar constituent fractions tends to converge to a similar value under the high heating rate.

Lithological mapping of Northwest Argentina with remote sensing data using tonal, textural and contextual features

Tonal, textural and contextual properties are used in manual photointerpretation of remotely sensed data. This study has used these three attributes to produce a lithological map of semi arid northwest Argentina by semi automatic computer classification procedures of remotely sensed data. Three different types of satellite data were investigated, these were LANDSAT MSS, TM and SIR-A imagery. Supervised classification procedures using tonal features only produced poor classification results. LANDSAT MSS produced classification accuracies in the range of 40 to 60%, while accuracies of 50 to 70% were achieved using LANDSAT TM data. The addition of SIR-A data produced increases in the classification accuracy. The increased classification accuracy of TM over the MSS is because of the better discrimination of geological materials afforded by the middle infra red bands of the TM sensor. The maximum likelihood classifier consistently produced classification accuracies 10 to 15% higher than either the minimum distance to means or decision tree classifier, this improved accuracy was obtained at the cost of greatly increased processing time. A new type of classifier the spectral shape classifier, which is computationally as fast as a minimum distance to means classifier is described. However, the results for this classifier were disappointing, being lower in most cases than the minimum distance or decision tree procedures. The classification results using only tonal features were felt to be unacceptably poor, therefore textural attributes were investigated. Texture is an important attribute used by photogeologists to discriminate lithology. In the case of TM data, texture measures were found to increase the classification accuracy by up to 15%. However, in the case of the LANDSAT MSS data the use of texture measures did not provide any significant increase in the accuracy of classification. For TM data, it was found that second order texture, especially the SGLDM based measures, produced highest classification accuracy. Contextual post processing was found to increase classification accuracy and improve the visual appearance of classified output by removing isolated misclassified pixels which tend to clutter classified images. Simple contextual features, such as mode filters were found to out perform more complex features such as gravitational filter or minimal area replacement methods. Generally the larger the size of the filter, the greater the increase in the accuracy. Production rules were used to build a knowledge based system which used tonal and textural features to identify sedimentary lithologies in each of the two test sites. The knowledge based system was able to identify six out of ten lithologies correctly.

Sequential pyrolysis and catalytic low temperature reforming of wheat straw

Aim of the work is the implementation of a low temperature reforming (LT reforming) unit downstream the Haloclean pyrolyser in order to enhance the heating value of the pyrolysis gas. Outside the focus of this work was to gain a synthesis gas quality for further use. Temperatures between 400 °C and 500 °C were applied. A commercial pre-reforming catalyst on a nickel basis from Südchemie was chosen for LT reforming. As biogenic feedstock wheat straw has been used. Pyrolysis of wheat straw at 450 °C by means of Haloclean pyrolysis leads to 28% of char, 50% of condensate and 22% of gas. The condensate separates in a water phase and an organic phase. The organic phase is liquid, but contains viscous compounds. These compounds could underlay aging and could lead to solid tars which can cause post processing problems. Therefore, the implementation of a catalytic reformer is not only of interest from an energetic point of view, it is generally interesting for tar conversion purposes after pyrolysis applications. By using a fixed bed reforming unit at 450–490 °C and space velocities about 3000 l/h the pyrolysis gas volume flow could be increased to about 58%. This corresponds to a decrease of the yields of condensates by means of catalysis up to 17%, the yield of char remains unchanged, since pyrolysis conditions are the same. The heating value in the pyrolysis gas could be increased by the factor of 1.64. Hydrogen concentrations up to 14% could be realised.

Dissipative nonlinear structures in optical fiber systems II : self-similar parabolic pulses in active fibers

In this second talk on dissipative structures in fiber applications, we overview theoretical aspects of the generation, evolution and characterization of self-similar parabolic-shaped pulses in fiber amplifier media. In particular, we present a perturbation analysis that describes the structural changes induced by third-order fiber dispersion on the parabolic pulse solution of the nonlinear Schrödinger equation with gain. Promising applications of parabolic pulses in optical signal post-processing and regeneration in communication systems are also discussed.

Fast pyrolysis and nitrogenolysis of biomass and biogenic residues:production of a sustainable slow release fertiliser

The production of agricultural and horticultural products requires the use of nitrogenous fertiliser that can cause pollution of surface and ground water and has a large carbon footprint as it is mainly produced from fossil fuels. The overall objective of this research project was to investigate fast pyrolysis and in-situ nitrogenolysis of biomass and biogenic residues as an alternative route to produce a sustainable solid slow release fertiliser mitigating the above stated problems. A variety of biomasses and biogenic residues were characterized by proximate analysis, ultimate analysis, thermogravimetric analysis (TGA) and Pyrolysis – Gas chromatography – Mass Spectroscopy (Py–GC–MS) for their potential use as feedstocks using beech wood as a reference material. Beech wood was virtually nitrogen free and therefore suitable as a reference material as added nitrogen can be identified as such while Dried Distillers Grains with Solubles (DDGS) and rape meal had a nitrogen content between 5.5wt.% and 6.1wt.% qualifying them as high nitrogen feedstocks. Fast pyrolysis and in-situ nitrogenolysis experiments were carried out in a continuously fed 1kg/h bubbling fluidized bed reactor at around 500°C quenching the pyrolysis vapours with isoparaffin. In-situ nitrogenolysis experiments were performed by adding ammonia gas to the fast pyrolysis reactor at nominal nitrogen addition rates between 5wt.%C and 20wt.%C based on the dry feedstock’s carbon content basis. Mass balances were established for the processing experiments. The fast pyrolysis and in-situ nitrogenolysis products were characterized by proximate analysis, ultimate analysis and GC– MS. High liquid yields and good mass balance closures of over 92% were obtained. The most suitable nitrogen addition rate for the in-situ nitrogenolysis experiments was determined to be 12wt.%C on dry feedstock carbon content basis. However, only a few nitrogen compounds that were formed during in-situ nitrogenolysis could be identified by GC–MS. A batch reactor process was developed to thermally solidify the fast pyrolysis and in-situ nitrogenolysis liquids of beech wood and Barley DDGS producing a brittle solid product. This was obtained at 150°C with an addition of 2.5wt% char (as catalyst) after a processing time of 1h. The batch reactor was also used for modifying and solidifying fast pyrolysis liquids derived from beech wood by adding urea or ammonium phosphate as post processing nitrogenolysis. The results showed that this type of combined approach was not suitable to produce a slow release fertiliser, because the solid product contained up to 65wt.% of highly water soluble nitrogen compounds that would be released instantly by rain. To complement the processing experiments a comparative study via Py–GC–MS with inert and reactive gas was performed with cellulose, hemicellulose, lignin and beech wood. This revealed that the presence of ammonia gas during analytical pyrolysis did not appear to have any direct impact on the decomposition products of the tested materials. The chromatograms obtained showed almost no differences between inert and ammonia gas experiments indicating that the reaction between ammonia and pyrolysis vapours does not occur instantly. A comparative study via Fourier Transformed Infrared Spectroscopy of solidified fast pyrolysis and in-situ nitrogenolysis products showed that there were some alterations in the spectra obtained. A shift in frequencies indicating C=O stretches typically related to the presence of carboxylic acids to C=O stretches related to amides was observed and no double or triple bonded nitrogen was detected. This indicates that organic acids reacted with ammonia and that no potentially harmful or non-biodegradable triple bonded nitrogen compounds were formed. The impact of solid slow release fertiliser (SRF) derived from pyrolysis and in-situ nitrogenolysis products from beech wood and Barley DDGS on microbial life in soils and plant growth was tested in cooperation with Rothamsted Research. The microbial incubation tests indicated that microbes can thrive on the SRFs produced, although some microbial species seem to have a reduced activity at very high concentrations of beech wood and Barley DDGS derived SRF. The plant tests (pot trials) showed that the application of SRF derived from beech wood and barley DDGS had no negative impact on germination or plant growth of rye grass. The fertilizing effect was proven by the dry matter yields in three harvests after 47 days, 89 days and 131 days. The findings of this research indicate that in general a slow release fertiliser can be produced from biomass and biogenic residues by in-situ nitrogenolysis. Nevertheless the findings also show that additional research is necessary to identify which compounds are formed during this process.

The intermediate pyrolysis and catalytic steam reforming of brewers spent grain

Brewers spent grain (BSG) is a widely available feedstock representing approximately 85% of the total by-products generated in the brewing industry. This is currently either disposed of to landfill or used as cattle feed due to its high protein content. BSG has received little or no attention as a potential energy resource, but increasing disposal costs and environmental constraints are now prompting the consideration of this. One possibility for the utilisation of BSG for energy is via intermediate pyrolysis to produce gases, vapours and chars. Intermediate pyrolysis is characterised by indirect heating in the absence of oxygen for short solids residence times of a few minutes, at temperatures of 350-450 °C. In the present work BSG has been characterised by chemical, proximate, ultimate and thermo-gravimetric analysis. Intermediate pyrolysis of BSG at 450 °C was carried out using a twin coaxial screw reactor known as Pyroformer to give yields of char 29%, 51% of bio-oil and 19% of permanent gases. The bio-oil liquid was found to separate in to an aqueous phase and organic phase. The organic phase contained viscous compounds that could age over time leading to solid tars that can present problems in CHP application. The quality of the pyrolysis vapour products before quenching can be upgraded to achieve much improved suitability as a fuel by downstream catalytic reforming. A Bench Scale batch pyrolysis reactor has then been used to pyrolyse small samples of BSG under a range of conditions of heating rate and temperature simulating the Pyroformer. A small catalytic reformer has been added downstream of the reactor in which the pyrolysis vapours can be further cracked and reformed. A commercial reforming nickel catalyst was used at 500, 750 and 850 °C at a space velocity about 10,000 L/h with and without the addition of steam. Results are presented for the properties of BSG, and the products of the pyrolysis process both with and without catalytic post-processing. Results indicate that catalytic reforming produced a significant increase in permanent gases mainly (H2 and CO) with H2 content exceeding 50 vol% at higher reforming temperatures. Bio-oil yield decreased significantly as reforming temperature increased with char remaining the same as pyrolysis condition remained unchanged. The process shows an increase in heating value for the product gas ranging between 10.8-25.2 MJ/m as reforming temperature increased. © 2012 Elsevier B.V. All rights reserved.

Pilot investigation of the changes in cortical activation during facial affect recognition with lamotrigine monotherapy in bipolar disorder

Background: Bipolar disorder is associated with dysfunction in prefrontal and limbic areas implicated in emotional processing. Aims: To explore whether lamotrigine monotherapy may exert its action by improving the function of the neural network involved in emotional processing. Method: We used functional magnetic resonance imaging to examine changes in brain activation during a sad facial affect recognition task in 12 stable patients with bipolar disorder when medication-free compared with healthy controls and after 12 weeks of lamotrigine monotherapy. Results: At baseline, compared with controls, patients with bipolar disorder showed overactivity in temporal regions and underactivity in the dorsal medial and right ventrolateral prefrontal cortex, and the dorsal cingulate gyrus. Following lamotrigine monotherapy, patients demonstrated reduced temporal and increased prefrontal activation. Conclusions: This preliminary evidence suggests that lamotrigine may enhance the function of the neural circuitry involved in affect recognition.

Dissipative nonlinear structures in optical fiber systems II:self-similar parabolic pulses in active fibers

In this second talk on dissipative structures in fiber applications, we overview theoretical aspects of the generation, evolution and characterization of self-similar parabolic-shaped pulses in fiber amplifier media. In particular, we present a perturbation analysis that describes the structural changes induced by third-order fiber dispersion on the parabolic pulse solution of the nonlinear Schrödinger equation with gain. Promising applications of parabolic pulses in optical signal post-processing and regeneration in communication systems are also discussed.

Permanent matching of coupled optical bottle resonators with better than 0.16 GHz precision

The fabrication precision is one of the most critical challenges to the creation of practical photonic circuits composed of coupled high Q-factor microresonators. While very accurate transient tuning of microresonators based on local heating has been reported, the record precision of permanent resonance positioning achieved by post-processing is still within 1 and 5 GHz. Here we demonstrate two coupled bottle microresonators fabricated at the fiber surface with resonances that are matched with a better than 0.16 GHz precision. This corresponds to a better than 0.17 Å precision in the effective fiber radius variation. The achieved fabrication precision is only limited by the resolution of our optical spectrum analyzer and can be potentially improved by an order of magnitude.

Computational fluid dynamics (CFD) based approach to consequence assessment of accidental release of hydrocarbon during storage and transportation

This thesis investigated the risk of accidental release of hydrocarbons during transportation and storage. Transportation of hydrocarbons from an offshore platform to processing units through subsea pipelines involves risk of release due to pipeline leakage resulting from corrosion, plastic deformation caused by seabed shakedown or damaged by contact with drifting iceberg. The environmental impacts of hydrocarbon dispersion can be severe. Overall safety and economic concerns of pipeline leakage at subsea environment are immense. A large leak can be detected by employing conventional technology such as, radar, intelligent pigging or chemical tracer but in a remote location like subsea or arctic, a small chronic leak may be undetected for a period of time. In case of storage, an accidental release of hydrocarbon from the storage tank could lead pool fire; further it could escalate to domino effects. This chain of accidents may lead to extremely severe consequences. Analyzing past accident scenarios it is observed that more than half of the industrial domino accidents involved fire as a primary event, and some other factors for instance, wind speed and direction, fuel type and engulfment of the compound. In this thesis, a computational fluid dynamics (CFD) approach is taken to model the subsea pipeline leak and the pool fire from a storage tank. A commercial software package ANSYS FLUENT Workbench 15 is used to model the subsea pipeline leakage. The CFD simulation results of four different types of fluids showed that the static pressure and pressure gradient along the axial length of the pipeline have a sharp signature variation near the leak orifice at steady state condition. Transient simulation is performed to obtain the acoustic signature of the pipe near leak orifice. The power spectral density (PSD) of acoustic signal is strong near the leak orifice and it dissipates as the distance and orientation from the leak orifice increase. The high-pressure fluid flow generates more noise than the low-pressure fluid flow. In order to model the pool fire from the storage tank, ANSYS CFX Workbench 14 is used. The CFD results show that the wind speed has significant contribution on the behavior of pool fire and its domino effects. The radiation contours are also obtained from CFD post processing, which can be applied for risk analysis. The outcome of this study will be helpful for better understanding of the domino effects of pool fire in complex geometrical settings of process industries. The attempt to reduce and prevent risks is discussed based on the results obtained from the numerical simulations of the numerical models.

AWI Bathymetric Chart of the Gakkel Ridge (AWI BCGR) (Scale 1:150,000)

The 1 : 1,500,000 AWI Bathymetric Chart of the Gakkel Ridge (AWI BCGR) has been developed from multibeam data measured during the Arctic Mid-Ocean Ridge Expedition in 2001 (AMORE 2001, ARK-XVII/2). This expedition was conducted to investigate the Gakkel Ridge in the Arctic Ocean and was carried out by the icebreaking research vessels RV Polarstern and USCGC Healy. Polarstern is equipped with the multibeam sonar system Hydrosweep DS-2, whereas Healy carries Seabeam 2112. During the expedition an area of 8890 km length and 18 - 46 km width, situated between 82°N/8°W and 87°N/75°E, was surveyed simultaneously by both vessels. Water depths ranged from 566 to 5673 meters. Dense sea ice cover derogated the sonar measurements and decreased data quality. Data errors were corrected in an extensive post-processing. The data of two different sonar systems had to be consolidated in order to derive a high resolution bathymetry of the Gakkel Ridge. Final result was a digital terrain model (DTM) with a grid spacing of 100 meters, which was utilized for generating the map series AWI Bathymetric Chart of the Gakkel Ridge, consisting of ten map sheets.

RF MIMO Systems for Wide-Area Indoor Human Motion Monitoring

Human motion monitoring is an important function in numerous applications. In this dissertation, two systems for monitoring motions of multiple human targets in wide-area indoor environments are discussed, both of which use radio frequency (RF) signals to detect, localize, and classify different types of human motion. In the first system, a coherent monostatic multiple-input multiple-output (MIMO) array is used, and a joint spatial-temporal adaptive processing method is developed to resolve micro-Doppler signatures at each location in a wide-area for motion mapping. The downranges are obtained by estimating time-delays from the targets, and the crossranges are obtained by coherently filtering array spatial signals. Motion classification is then applied to each target based on micro-Doppler analysis. In the second system, multiple noncoherent multistatic transmitters (Tx's) and receivers (Rx's) are distributed in a wide-area, and motion mapping is achieved by noncoherently combining bistatic range profiles from multiple Tx-Rx pairs. Also, motion classification is applied to each target by noncoherently combining bistatic micro-Doppler signatures from multiple Tx-Rx pairs. For both systems, simulation and real data results are shown to demonstrate the ability of the proposed methods for monitoring patient repositioning activities for pressure ulcer prevention.

Assessment of the Spatial and Temporal Distribution of Functional Connectivity in Resting-State BOLD fMRI

Recent research into resting-state functional magnetic resonance imaging (fMRI) has shown that the brain is very active during rest. This thesis work utilizes blood oxygenation level dependent (BOLD) signals to investigate the spatial and temporal functional network information found within resting-state data, and aims to investigate the feasibility of extracting functional connectivity networks using different methods as well as the dynamic variability within some of the methods. Furthermore, this work looks into producing valid networks using a sparsely-sampled sub-set of the original data.

In this work we utilize four main methods: independent component analysis (ICA), principal component analysis (PCA), correlation, and a point-processing technique. Each method comes with unique assumptions, as well as strengths and limitations into exploring how the resting state components interact in space and time.

Correlation is perhaps the simplest technique. Using this technique, resting-state patterns can be identified based on how similar the time profile is to a seed region’s time profile. However, this method requires a seed region and can only identify one resting state network at a time. This simple correlation technique is able to reproduce the resting state network using subject data from one subject’s scan session as well as with 16 subjects.

Independent component analysis, the second technique, has established software programs that can be used to implement this technique. ICA can extract multiple components from a data set in a single analysis. The disadvantage is that the resting state networks it produces are all independent of each other, making the assumption that the spatial pattern of functional connectivity is the same across all the time points. ICA is successfully able to reproduce resting state connectivity patterns for both one subject and a 16 subject concatenated data set.

Using principal component analysis, the dimensionality of the data is compressed to find the directions in which the variance of the data is most significant. This method utilizes the same basic matrix math as ICA with a few important differences that will be outlined later in this text. Using this method, sometimes different functional connectivity patterns are identifiable but with a large amount of noise and variability.

To begin to investigate the dynamics of the functional connectivity, the correlation technique is used to compare the first and second halves of a scan session. Minor differences are discernable between the correlation results of the scan session halves. Further, a sliding window technique is implemented to study the correlation coefficients through different sizes of correlation windows throughout time. From this technique it is apparent that the correlation level with the seed region is not static throughout the scan length.

The last method introduced, a point processing method, is one of the more novel techniques because it does not require analysis of the continuous time points. Here, network information is extracted based on brief occurrences of high or low amplitude signals within a seed region. Because point processing utilizes less time points from the data, the statistical power of the results is lower. There are also larger variations in DMN patterns between subjects. In addition to boosted computational efficiency, the benefit of using a point-process method is that the patterns produced for different seed regions do not have to be independent of one another.

This work compares four unique methods of identifying functional connectivity patterns. ICA is a technique that is currently used by many scientists studying functional connectivity patterns. The PCA technique is not optimal for the level of noise and the distribution of the data sets. The correlation technique is simple and obtains good results, however a seed region is needed and the method assumes that the DMN regions is correlated throughout the entire scan. Looking at the more dynamic aspects of correlation changing patterns of correlation were evident. The last point-processing method produces a promising results of identifying functional connectivity networks using only low and high amplitude BOLD signals.