Comparison of Routine Brain Imaging at 3 T and 7 T.

OBJECTIVE The aim of this study was to compare quantitative and semiquantitative parameters (signal-to-noise ratio [SNR], contrast-to-noise ratio [CNR], image quality, diagnostic confidence) from a standard brain magnetic resonance imaging examination encompassing common neurological disorders such as demyelinating disease, gliomas, cerebrovascular disease, and epilepsy, with comparable sequence protocols and acquisition times at 3 T and at 7 T. MATERIALS AND METHODS Ten healthy volunteers and 4 subgroups of 40 patients in total underwent comparable magnetic resonance protocols with standard diffusion-weighted imaging, 2D and 3D turbo spin echo, 2D and 3D gradient echo and susceptibility-weighted imaging of the brain (10 sequences) at 3 T and 7 T. The subgroups comprised patients with either lesional (n = 5) or nonlesional (n = 4) epilepsy, intracerebral tumors (n = 11), demyelinating disease (n = 11) (relapsing-remitting multiple sclerosis [MS, n = 9], secondary progressive MS [n = 1], demyelinating disease not further specified [n = 1]), or chronic cerebrovascular disorders [n = 9]). For quantitative analysis, SNR and CNR were determined. For a semiquantitative assessment of the diagnostic confidence, a 10-point scale diagnostic confidence score (DCS) was applied. Two experienced radiologists with additional qualification in neuroradiology independently assessed, blinded to the field strength, 3 pathology-specific imaging criteria in each of the 4 disease groups and rated their diagnostic confidence. The overall image quality was semiquantitatively assessed using a 4-point scale taking into account whether diagnostic decision making was hampered by artifacts or not. RESULTS Without correction for spatial resolution, SNR was higher at 3 T except in the T2 SPACE 3D, DWI single shot, and DIR SPACE 3D sequences. The SNR corrected by the ratio of 3 T/7 T voxel sizes was higher at 7 T than at 3 T in 10 of 11 sequences (all except for T1 MP2RAGE 3D).In CNR, there was a wide variation between sequences and patient cohorts, but average CNR values were broadly similar at 3 T and 7 T.DCS values for all 4 pathologic entities were higher at 7 T than at 3 T. The DCS was significantly higher at 7 T for diagnosis and exclusion of cortical lesions in vascular disease. A tendency to higher DCS at 7 T for cortical lesions in MS was observed, and for the depiction of a central vein and iron deposits within MS lesions. Despite motion artifacts, DCS values were higher at 7 T for the diagnosis and exclusion of hippocampal sclerosis in mesial temporal lobe epilepsy (improved detection of the hippocampal subunits). Interrater agreement was 69.7% at 3 T and 93.3% at 7 T. There was no significant difference in the overall image quality score between 3 T and 7 T taking into account whether diagnostic decision making was hampered by artifacts or not. CONCLUSIONS Ultra-high-field magnetic resonance imaging at 7 T compared with 3 T yielded an improved diagnostic confidence in the most frequently encountered neurologic disorders. Higher spatial resolution and contrast were identified as the main contributory factors.

Structural brain correlates of defective gesture performance in schizophrenia

INTRODUCTION The neural correlates of impaired performance of gestures are currently unclear. Lesion studies showed variable involvement of the ventro-dorsal stream particularly left inferior frontal gyrus (IFG) in gesture performance on command. However, findings cannot be easily generalized as lesions may be biased by the architecture of vascular supply and involve brain areas beyond the critical region. The neuropsychiatric syndrome of schizophrenia shares apraxic-like errors and altered brain structure without macroanatomic lesions. Schizophrenia may therefore qualify as a model disorder to test neural correlates of gesture impairments. METHODS We included 45 schizophrenia patients and 44 healthy controls in the study to investigate the structural brain correlates of defective gesturing in schizophrenia using voxel based morphometry. Gestures were tested in two domains: meaningful gestures (transitive and intransitive) on verbal command and imitation of meaningless gestures. Cut-off scores were used to separate patients with deficits, patients without deficits and controls. Group differences in gray matter (GM) volume were explored in an ANCOVA. RESULTS Patients performed poorer than controls in each gesture category (p < .001). Patients with deficits in producing meaningful gestures on command had reduced GM predominantly in left IFG, with additional involvement of right insula and anterior cingulate cortex. Patients with deficits differed from patients without deficits in right insula, inferior parietal lobe (IPL) and superior temporal gyrus. CONCLUSIONS Impaired performance of meaningful gestures on command was linked to volume loss predominantly in the praxis network in schizophrenia. Thus, the behavioral similarities between apraxia and schizophrenia are paralleled by structural alterations. However, few associations between behavioral impairment and structural brain alterations appear specific to schizophrenia.

Geomorphology of the Imhotep region on comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Context. Since August 2014, the OSIRIS Narrow Angle Camera (NAC) onboard the Rosetta spacecraft has acquired high spatial resolution images of the nucleus of comet 67P/Churyumov-Gerasimenko, down to the decimeter scale. This paper focuses on the Imhotep region, located on the largest lobe of the nucleus, near the equator. Aims. We map, inventory, and describe the geomorphology of the Imhotep region. We propose and discuss some processes to explain the formation and ongoing evolution of this region. Methods. We used OSIRIS NAC images, gravitational heights and slopes, and digital terrain models to map and measure the morphologies of Imhotep. Results. The Imhotep region presents a wide variety of terrains and morphologies: smooth and rocky terrains, bright areas, linear features, roundish features, and boulders. Gravity processes such as mass wasting and collapse play a significant role in the geomorphological evolution of this region. Cometary processes initiate erosion and are responsible for the formation of degassing conduits that are revealed by elevated roundish features on the surface. We also propose a scenario for the formation and evolution of the Imhotep region; this implies the presence of large primordial voids inside the nucleus, resulting from its formation process.

Two independent and primitive envelopes of the bilobate nucleus of comet 67P

The factors shaping cometary nuclei are still largely unknown, but could be the result of concurrent effects of evolutionary(1,2) and primordial processes(3,4). The peculiar bilobed shape of comet 67P/Churyumov-Gerasimenko may be the result of the fusion of two objects that were once separate or the result of a localized excavation by outgassing at the interface between the two lobes(5). Here we report that the comet's major lobe is enveloped by a nearly continuous set of strata, up to 650 metres thick, which are independent of an analogous stratified envelope on the minor lobe. Gravity vectors computed for the two lobes separately are closer to perpendicular to the strata than those calculated for the entire nucleus and adjacent to the neck separating the two lobes. Therefore comet 67P/Churyumov-Gerasimenko is an accreted body of two distinct objects with 'onion-like' stratification, which formed before they merged. We conclude that gentle, low-velocity collisions occurred between two fully formed kilometre-sized cometesimals in the early stages of the Solar System. The notable structural similarities between the two lobes of comet 67P/Churyumov-Gerasimenko indicate that the early-forming cometesimals experienced similar primordial stratified accretion, even though they formed independently.

Size-frequency distribution of boulders ≥7 m on comet 67P/Churyumov-Gerasimenko

Aims. We derive for the first time the size-frequency distribution of boulders on a comet, 67P/Churyumov-Gerasimenko (67P), computed from the images taken by the Rosetta/OSIRIS imaging system. We highlight the possible physical processes that lead to these boulder size distributions. Methods. We used images acquired by the OSIRIS Narrow Angle Camera, NAC, on 5 and 6 August 2014. The scale of these images (2.44−2.03 m/px) is such that boulders ≥7 m can be identified and manually extracted from the datasets with the software ArcGIS. We derived both global and localized size-frequency distributions. The three-pixel sampling detection, coupled with the favorable shadowing of the surface (observation phase angle ranging from 48° to 53°), enables unequivocally detecting boulders scattered all over the illuminated side of 67P. Results. We identify 3546 boulders larger than 7 m on the imaged surface (36.4 km2), with a global number density of nearly 100/km2 and a cumulative size-frequency distribution represented by a power-law with index of −3.6 +0.2/−0.3. The two lobes of 67P appear to have slightly different distributions, with an index of −3.5 +0.2/−0.3 for the main lobe (body) and −4.0 +0.3/−0.2 for the small lobe (head). The steeper distribution of the small lobe might be due to a more pervasive fracturing. The difference of the distribution for the connecting region (neck) is much more significant, with an index value of −2.2 +0.2/−0.2. We propose that the boulder field located in the neck area is the result of blocks falling from the contiguous Hathor cliff. The lower slope of the size-frequency distribution we see today in the neck area might be due to the concurrent processes acting on the smallest boulders, such as i) disintegration or fragmentation and vanishing through sublimation; ii) uplifting by gas drag and consequent redistribution; and iii) burial beneath a debris blanket. We also derived the cumulative size-frequency distribution per km2 of localized areas on 67P. By comparing the cumulative size-frequency distributions of similar geomorphological settings, we derived similar power-law index values. This suggests that despite the selected locations on different and often opposite sides of the comet, similar sublimation or activity processes, pit formation or collapses, as well as thermal stresses or fracturing events occurred on multiple areas of the comet, shaping its surface into the appearance we see today.

Shape model, reference system definition, and cartographic mapping standards for comet 67P/Churyumov-Gerasimenko Stereo-photogrammetric analysis of Rosetta/OSIRIS image data

We analyzed more than 200 OSIRIS NAC images with a pixel scale of 0.9-2.4 m/pixel of comet 67P/Churyumov-Gerasimenko (67P) that have been acquired from onboard the Rosetta spacecraft in August and September 2014 using stereo-photogrammetric methods (SPG). We derived improved spacecraft position and pointing data for the OSIRIS images and a high-resolution shape model that consists of about 16 million facets (2 m horizontal sampling) and a typical vertical accuracy at the decimeter scale. From this model, we derive a volume for the northern hemisphere of 9.35 km(3) +/- 0.1 km(3). With the assumption of a homogeneous density distribution and taking into account the current uncertainty of the position of the comet's center-of-mass, we extrapolated this value to an overall volume of 18.7 km(3) +/- 1.2 km(3), and, with a current best estimate of 1.0 X 10(13) kg for the mass, we derive a bulk density of 535 kg/m(3) +/- 35 kg/m(3). Furthermore, we used SPG methods to analyze the rotational elements of 67P. The rotational period for August and September 2014 was determined to be 12.4041 +/- 0.0004 h. For the orientation of the rotational axis (z-axis of the body-fixed reference frame) we derived a precession model with a half-cone angle of 0.14 degrees, a cone center position at 69.54 degrees/64.11 degrees (RA/Dec J2000 equatorial coordinates), and a precession period of 10.7 days. For the definition of zero longitude (x-axis orientation), we finally selected the boulder-like Cheops feature on the big lobe of 67P and fixed its spherical coordinates to 142.35 degrees right-hand-rule eastern longitude and -0.28 degrees latitude. This completes the definition of the new Cheops reference frame for 67P. Finally, we defined cartographic mapping standards for common use and combined analyses of scientific results that have been obtained not only within the OSIRIS team, but also within other groups of the Rosetta mission.

Learning in the cerebellar nucleus is necessary for acquisition of Pavlovian eyelid conditioning

The cumulative work presented here supports the hypothesis that plasticity in the cerebellar cortex and cerebellar nuclei mediates a simple associative form of motor teaming-Pavlovian eyelid conditioning. It was previously demonstrated that focal ablative lesions of cerebellar anterior lobe or pharmacological block of the cerebellar cortex output disrupted the timing of the conditioned eyeblink response, unmasking a response with a relatively fixed and very short latency to onset. The results of this thesis demonstrate that the short-latency responses are due to associative learning. Unpaired training does not support the acquisition of short-latency responses while the rate of acquisition of short-latency responses during paired training is approximately the same as that of timed conditioned responses. The acquisition of short-latency responses is dependent on an intact cerebellar cortex. Both ablative lesions of the cerebellar cortex and inactivation of cerebellar cortex output with picrotoxin block the acquisition of short-latency responses. However, once the short-latency responses are acquired neither disconnection of cerebellar cortex nor inactivation of the cerebellar nucleus block reacquisition. The results are consistent with the proposal that plasticity in the cerebellar cortex is necessary for learning the timing of conditioned responses, plasticity in the interpositus nucleus mediates the short latency responses, and cerebellar cortical output and mossy fiber input are necessary for the acquisition of short latency responses. ^

The contributions of the amygdala, orbital frontal cortex and hippocampal formation to primate social cognition

Deficits in social cognition are prominent symptoms of many human psychiatric disorders, but the origin of such deficits remains largely unknown. To further current knowledge regarding the neural network mediating social cognition, the present research program investigated the individual contributions of two temporal lobe structures, the amygdala and hippocampal formation, and one frontal lobe region, the orbital frontal cortex (Areas 11 and 13), to primate social cognition. Based on previous research, we hypothesized that the amygdala, hippocampal formation and orbital frontal cortex contribute significantly to the formation of new social relationships, but less to the maintenance of familiar ones. ^ Thirty-six male rhesus macaques (Macaca mulatta) served as subjects, and were divided into four experimental groups: Neurotoxic amygdala lesion (A-ibo, n = 9), neurotoxic or aspiration orbital frontal cortex lesion (O, n = 9), neurotoxic hippocampal formation lesion (H-ibo, n = 9) or sham-operated control (C, n = 9). Six social groups (tetrads) were created, each containing one member from each experimental group. The effect of lesion on established social relationships was assessed during pre- and post-surgical unrestrained social interactions, whereas the effect of lesion on the formation of new relationships was assessed during an additional phase of post-surgical testing with shuffled tetrad membership. Results indicated that these three neural structures each contribute significantly to both the formation and maintenance of social relationships. Furthermore, the amygdala appears to primarily mediate normal responses to threatening social signals, whereas the orbital frontal cortex plays a more global role in social cognition by mediating responses to both threatening and affiliative social signals. By contrast, the hippocampal formation seems to contribute to social cognition indirectly by providing access to previous experience during social judgments. ^ These conclusions were further investigated with three experiments that measured behavioral and physiological (stress hormone) reactivity to threatening stimuli, and three additional experiments that measured subjects' ability to flexibly alter behavioral responses depending on the incentive value of a food reinforcer. Data from these six experiments further confirmed and strengthened the three conclusions originating from the social behavior experiments and, when combined with the current literature, helped to formulate a simple, but testable, theoretical model of primate social cognition. ^

Evaluation of Mycobacterium avium complex lung disease in women

Mycobacterium avium complex (MAC) is a ubiquitous organism responsible for most pulmonary and disseminated disease caused by non-tuberculosis (NTM) mycobacteria. Though MAC lung disease without predisposing factors is uncommon, in recent years it has been increasingly described in middle-aged and elderly women. Recognition and correct diagnosis, is often delayed due to the indolent nature of the disease. It is unclear if these women have significant clinical disease as or if their airways are simply colonized by the bacterium. This study describes the clinical presentation, identifies risk factors, and describes the clinical significance of MAC lung disease in HIV-negative women aged 50 or greater. ^ A hybrid study design utilizing both cross-sectional and case-control methodologies was used. A comparison population was selected from previously identified tuberculosis suspects found throughout Harris County. The study population had at least one acid fast bacillus pulmonary culture performed between 1/1/1998 and 12/31/2000 from a pulmonary source. Clinical presentation and symptoms were analyzed using a cross-sectional design. Past medical history and other risk factors were evaluated using a traditional case-control study design. Differences in categorical variables were estimated with the Chi Square or Fisher's Exact test as appropriate. Odds ratios and 95% confidence intervals were utilized to evaluate associations. Multivariate logistic regression was used to identify predictive factors for MAC. All statistical tests were two-sided and P-values <0.05 were considered statistically significant. ^ Culture confirmed MAC pulmonary cases were more likely to be white, have bronchiectasis, scoliosis, evidence of cavitation and pleural changes on chest radiography and granulomas on histopathologic examination than women whose pulmonary cultures were AFB negative. After controlling for selected risk factors, white race continued to be significantly associated with MAC lung disease (OR = 4.6, 95% CI = 2.3, 9.2). In addition, asthma history, smoking history and alcohol use were less likely to be evident among MAC cases in a multivariate analysis. Right upper and right middle lobe disease was further noted among clinically significant cases. Based on population data, MAC lung disease appears to represent a significant clinical syndrome in HIV-negative women thus supporting the theory of the Lady Windermere Syndrome. ^

Contribution of transmembrane and cytoplasmic domains to membrane protein association

Membrane proteins are critical to every aspect of cell physiology, with their association mediating important biological functions. The transmembrane and cytoplasmic domains are known to be important for their association. In order to characterize their role in detail, we have applied different biophysical techniques in detergent micelles to two model systems. The first study involves FcRγ, a single transmembrane domain protein existing as a disulfide linked homodimer. We investigated the role of a conserved transmembrane polar residue and the cytoplasmic tail in FcRγ homo-interactions. Our results by various biophysical techniques including SDS-PAGE, circular dichroism and sedimentation equilibrium in detergent micelles indicate importance of both the transmembrane polar residue and cytoplasmic tail in maintaining proper conformation for FcRγ homo-interactions. A contrasting second study was on L-selectin, another single transmembrane domain protein with a large extracellular domain and a short cytoplasmic tail. Previous cross-linking experiments indicate its possible dimerization. However, the purified fragment of L-selectin and corresponding mutants did not dimerize when analyzed by TOXCAT assay, sedimentation equilibrium and fluorescence resonance energy transfer. It was likely that the presence of juxtamembrane positively charged residues led to decreased migrational rates in SDS PAGE. In conclusion, complementary biophysical techniques should be used with care when studying membrane protein association in detergent micelles. As an extension to our study on L-selectin, we also investigated its interaction with Calmodulin (CaM) in detergent micelles. CaM was found to interact with different detergents. We applied fluorescence and NMR spectroscopy to characterize the interaction of both the apo and Ca 2+ bound form of CaM, with commonly used detergents, below and above their respective critical micelle concentrations. The interaction of apo-CaM with detergents was found to vary with the nature of the detergent head group, whereas Ca2+-CaM interacted with individual detergent molecules irrespective of the nature of their head group. NMR titration experiments of CaM with detergents indicated involvement of specific residues from the N-lobe, linker and C-lobe of CaM. ^

Identification of a conserved cluster in the RH domain of GRK critical for activation by GPCRs

One of the most critical aspects of G Protein Coupled Receptors (GPCRs) regulation is their rapid and acute desensitization following agonist stimulation. Phosphorylation of these receptors by GPCR kinases (GRK) is a major mechanism of desensitization. Considerable evidence from studies of rhodopsin kinase and GRK2 suggests there is an allosteric docking site for the receptor distinct from the GRK catalytic site. While the agonist-activated GPCR appears crucial for GRK activation, the molecular details of this interaction remain unclear. Recent studies suggested an important role for the N- and C-termini and domains in the small lobe of the kinase domain in allosteric activation; however, neither the mechanism of action of that site nor the RH domain contributions have been elucidated. To search for the allosteric site, we first indentified evolutionarily conserved sites within the RH and kinase domains presumably deterministic of protein function employing evolutionary trace (ET) methodology and crystal structures of GRK6. Focusing on a conserved cluster centered on helices 3, 9, and 10 in the RH domain, key residues of GRK5 and 6 were targeted for mutagenesis and functional assays. We found that a number of double mutations within helices 3, 9, and 10 and the N-terminus markedly reduced (50–90%) the constitutive phosphorylation of the β-2 Adrenergic Receptor (β2AR) in intact cells and phosphorylation of light-activated rhodopsin (Rho*) in vitro as compared to wild type (WT) GRK5 or 6. Based on these results, we designed peptide mimetics of GRK5 helix 9 both computationally and through chemical modifications with the goal of both confirming the importance of helix 9 and developing a useful inhibitor to disrupt the GPCR-GRK interaction. Several peptides were found to block Rho* phosphorylation by GRK5 including the native helix 9 sequence, Peptide Builder designed-peptide preserving only the key ET residues, and chemically locked helices. Most peptidomimetics showed inhibition of GRK5 activity greater than 80 % with an IC50 of ∼ 30 µM. Alanine scanning of helix 9 has further revealed both essential and non-essential residues for inhibition. Importantly, substitution of Arg 169 by an alanine in the native helix 9-based peptide gave an almost complete inhibition at 30 µM with an IC50 of ∼ 10 µM. In summary we report a previously unrecognized crucial role for the RH domain of GRK5 and 6, and the subsequent identification of a lead peptide inhibitor of protein-protein interaction with potential for specific blockade of GPCR desensitization. ^

CONFORMATIONAL DYNAMICS OF K-RAS AND H-RAS PROTEINS: IS THERE FUNCTIONAL SPECIFICITY AT THE CATALYTIC DOMAIN?

Ras proteins serve as crucial signaling modulators in cell proliferation through their ability to hydrolyze GTP and exist in a GTP “on” state and GTP “off” state. There are three different human Ras isoforms: H-ras, N-ras and K-ras (4A and 4B). Although their sequence identity is very high at the catalytic domain, these isoforms differ in their ability to activate different effectors and hence different signaling pathways. Much of the previous work on this topic has attributed this difference to the hyper variable region of Ras proteins, which contains most of the sequence variance among the isoforms and encodes specificity for differential distribution in the membrane. However, we hypothesize that sequence variation on lobe II of Ras catalytic domain alters dynamics and leads to differential preference for different effectors or modulators. In this work, we used all atom molecular dynamics to analyze the dynamics in the catalytic domain of H-ras and K-ras. We have also analyzed the dynamics of a transforming mutant of H-ras and K-ras and further studied the dynamics of an effectorselective mutant of H-ras. Collectively we have determined that wild type K-ras is more dynamic than H-ras and that the structure of the effector binding loop more closely resembles that of the T35S Raf-selective mutant, possibly giving us a new view and insight into the v mode of effector specificity. Furthermore we have determined that specific mutations at the same location perturb the conformational equilibrium differently in H-ras and K-ras and that an enhanced oncogenic potential may arise from different structural perturbations for each point mutation of a specific isoform.

Clay mineraology and grain size composition of surface sediments from South Pacific Ocean

The reconstruction of low-latitude ocean-atmosphere interactions is one of the major issues of (paleo-)environmental studies. The trade winds, extending over 20° to 30° of latitude in both hemispheres, between the subtropical highs and the intertropical convergence zone, are major components of the atmospheric circulation and little is known about their long-term variability on geological time-scales, in particular in the Pacific sector. We present the modern spatial pattern of eolian-derived marine sediments in the eastern equatorial and subtropical Pacific (10°N to 25°S) as a reference data set for the interpretation of SE Pacific paleo-dust records. The terrigenous silt and clay fractions of 75 surface sediment samples have been investigated for their grain-size distribution and clay-mineral compositions, respectively, to identify their provenances and transport agents. Dust delivered to the southeast Pacific from the semi- to hyper-arid areas of Peru and Chile is rather fine-grained (4-8 µm) due to low-level transport within the southeast trade winds. Nevertheless, wind is the dominant transport agent and eolian material is the dominant terrigenous component west of the Peru-Chile Trench south of ~ 5°S. Grain-size distributions alone are insufficient to identify the eolian signal in marine sediments due to authigenic particle formation on the sub-oceanic ridges and abundant volcanic glass around the Galapagos Islands. Together with the clay-mineral compositions of the clay fraction, we have identified the dust lobe extending from the coasts of Peru and Chile onto Galapagos Rise as well as across the equator into the doldrums. Illite is a very useful parameter to identify source areas of dust in this smectite-dominated study area.

Lead isotopes of feldspar in North Atlantic Ocean sediments

The provenance of ice-rafted debris (IRD) deposited in the North Atlantic before, during, and after Heinrich event 2 has been determined through measuring the lead isotopic composition of single feldspar grains and multiple-grain composites from the larger than 150-µm size fraction, from cores from the eastern and western North Atlantic and from the Labrador Sea. Single-grain analyses are used to identify the specific continental sources of the IRD, whereas composite samples are used to assess the relative IRD contributions from different sources. All single grains from Heinrich layer 2 (H 2) as well as H 2 composites plot along a correlation line on a 207Pb/204Pb versus 206Pb/204Pb diagram characteristic of the Churchill province of the Canadian shield. This is yet another strong piece of evidence that this Heinrich event was dominated by a massive iceberg discharge of the Laurentide ice sheet lobe located over Hudson Bay. In contrast, single grains from the ambient glacial sediment (above and below H 2) have multiple sources: many of them also lie along the correlation line with H 2 grains, but many others have Pb signatures consistent with derivation from the Grenville province and the Appalachian range in North America and possibly from Scandinavia and Greenland. Composites from the ambient sediment generally lie well to the right of the H 2 reference line in agreement with the results of the single-grain analyses. The evidence provided by lead isotopes regarding the dominant role played by the Hudson Bay lobe of the Laurentide ice sheet in the development of the Heinrich events lends support to the binge/purge model advanced by MacAyeal [1993a, b] that invokes trapping of geothermal heat by the base of the icecap and subsequent basal melting as the mechanism that triggered the Heinrich events.