CHROMOSOMAL MAPPING, LINKAGE RELATIONSHIPS, AND EVOLUTIONARY STUDIES OF THE HUMAN ANONYMOUS DNA CLONE D1S1 (IN SITU HYBRIDIZATION, PRIMATES, GENE MAPPING, AUTOSOMAL DOMINANT RETINITIS PIGMENTOSA, GENE DUPLICATION)

D1S1, an anonymous human DNA clone originally called (lamda)Ch4-H3 or (lamda)H3, was the first single copy mapped to a human chromosome (1p36) by in situ hybridization. The chromosomal assignment has been confirmed in other laboratories by repeating the in situ hybridization but not by another method. In the present study, hybridization to a panel of hamster-human somatic cell hybrids revealed copies of D1S1 on both chromosomes 1 and 3. Subcloning D1S1 showed that the D1S1 clone itself is from chromosome 3, and the sequence detected by in situ hybridization is at least two copies of part of the chromosome 3 copy. This finding demonstrates the importance of verifying gene mapping with two methods and questions the accuracy of in situ hybridization mapping.^ Non-human mammals have only one copy of D1S1, and the non-human primate D1S1 map closely resembles the human chromosome 3 copy. Thus, the human chromosome 1 copies appear to be part of a very recent duplication that occurred after the divergence between humans and the other great apes.^ A moderately informative HindIII D1S1 RFLP was mapped to chromosome 3. This marker and 12 protein markers were applied to a linkage study of autosomal dominant retinitis pigmentosa (ADRP). None of the markers proved linkage, but adding the three families examined to previously published data raises the ADRP:Rh lod score to 1.92 at (THETA) = 0.30. ^

Determination of Δ9-tetrahydrocannabinolic acid A (Δ9-THCA-A) in whole blood and plasma by LC–MS/MS and application in authentic samples from drivers suspected of driving under the influence of cannabis

Delta-9-tetrahydrocannabinolic acid A (THCA-A) is the biosynthetic precursor of delta-9-tetrahydrocannabinol (THC) in cannabis plants, and has no psychotropic effects. THCA-A can be detected in blood and urine, and several metabolites have been identified. THCA-A was also shown to be incorporated in hair by side stream smoke to a minor extent, but incorporation via blood stream or sweat seems unlikely. The detection of THCA-A in biological fluids may serve as a marker for differentiating between the intake of prescribed THC medication – containing only pure THC – and cannabis products containing THC besides THC-acid A and other cannabinoids. However, the knowledge about its usefulness in forensic cases is very limited. The aim of the present work was the development of a reliable method for THCA-A determination in human blood or plasma using LC–MS/MS and application to cases of driving under the influence of drugs. Fifty eight (58) authentic whole blood and the respective plasma samples were collected from drivers suspected of driving under the influence of cannabis from the region of Bern (Switzerland). Samples were first tested for THC, 11-OH-THC and THC-COOH, and then additionally for THCA-A. For this purpose, the existing LC–MS/MS method was modified and validated, and found to be selective and linear over a range of 1.0 to 200 ng/mL (the correlation coefficients were above 0.9980 in all validation runs). Limit of detection (LOD) and limit of quantification (LOQ) were 0.3 ng/mL and 1.0 ng/mL respectively. Intra- and inter-assay accuracy were equal or better than 90% and intra- and inter-assay precision were equal or better than 11.1%. The mean extraction efficiencies were satisfactory being equal or higher than 85.4%. THCA-A was stable in whole blood samples after 3 freeze/thaw cycles and storage at 4 °C for 7 days. Re-injection (autosampler) stability was also satisfactory. THC was present in all blood samples with levels ranging from 0.7 to 51 ng/mL. THCA-A concentrations ranged from 1.0 to 496 ng/mL in blood samples and from 1.4 to 824 ng/mL in plasma samples. The plasma:blood partition coefficient had a mean value of 1.7 (±0.21, SD). No correlation was found between the degree of intoxication or impairment stated in the police protocols or reports of medical examinations and the detected THCA-A-concentration in blood.

TECPR2 Associated Neuroaxonal Dystrophy in Spanish Water Dogs.

Clinical, pathological and genetic examination revealed an as yet uncharacterized juvenile-onset neuroaxonal dystrophy (NAD) in Spanish water dogs. Affected dogs presented with various neurological deficits including gait abnormalities and behavioral deficits. Histopathology demonstrated spheroid formation accentuated in the grey matter of the cerebral hemispheres, the cerebellum, the brain stem and in the sensory pathways of the spinal cord. Iron accumulation was absent. Ultrastructurally spheroids contained predominantly closely packed vesicles with a double-layered membrane, which were characterized as autophagosomes using immunohistochemistry. The family history of the four affected dogs suggested an autosomal recessive inheritance. SNP genotyping showed a single genomic region of extended homozygosity of 4.5 Mb in the four cases on CFA 8. Linkage analysis revealed a maximal parametric LOD score of 2.5 at this region. By whole genome re-sequencing of one affected dog, a perfectly associated, single, non-synonymous coding variant in the canine tectonin beta-propeller repeat-containing protein 2 (TECPR2) gene affecting a highly conserved region was detected (c.4009C>T or p.R1337W). This canine NAD form displays etiologic parallels to an inherited TECPR2 associated type of human hereditary spastic paraparesis (HSP). In contrast to the canine NAD, the spinal cord lesions in most types of human HSP involve the sensory and the motor pathways. Furthermore, the canine NAD form reveals similarities to cases of human NAD defined by widespread spheroid formation without iron accumulation in the basal ganglia. Thus TECPR2 should also be considered as candidate gene for human NAD. Immunohistochemistry and the ultrastructural findings further support the assumption, that TECPR2 regulates autophagosome accumulation in the autophagic pathways. Consequently, this report provides the first genetic characterization of juvenile canine NAD, describes the histopathological features associated with the TECPR2 mutation and provides evidence to emphasize the association between failure of autophagy and neurodegeneration.

Traditional linkage analysis in admixed families

Currently there is no general method to study the impact of population admixture within families on the assumptions of random mating and consequently, Hardy-Weinberg equilibrium (HWE) and linkage equilibrium (LE) and on the inference obtained from traditional linkage analysis. ^ First, through simulation, the effect of admixture of two populations on the log of the odds (LOD) score was assessed, using Prostate Cancer as the typical disease model. Comparisons between simulated mixed and homogeneous families were performed. LOD scores under both models of admixture (within families and within a data set of homogeneous families) were closest to the homogeneous family scores of the population having the highest mixing proportion. Random sampling of families or ascertainment of families with disease affection status did not affect this observation, nor did the mode of inheritance (dominant/recessive) or sample size. ^ Second, after establishing the effect of admixture on the LOD score and inference for linkage, the presence of induced disequilibria by population admixture within families was studied and an adjustment procedure was developed. The adjustment did not force all disequilibria to disappear but because the families were adjusted for the population admixture, those replicates where the disequilibria exist are no longer affected by the disequilibria in terms of maximization for linkage. Furthermore, the adjustment was able to exclude uninformative families or families that had such a high departure from HWE and/or LE that their LOD scores were not reliable. ^ Together these observations imply that the presence of families of mixed population ancestry impacts linkage analysis in terms of the LOD score and the estimate of the recombination fraction. ^

The genetic basis of thoracic aortic aneurysms and dissections: Genetic heterogeneity and mapping of TAAD1 and TAAD2 loci

Thoracic aortic aneurysms leading to aortic dissections (TAAD) are a major cause of morbidity and mortality in the United States. TAAD is a complication of some known genetic disorders, such as Marfan syndrome and Turner syndrome, but the majority of familial cases are not due to a known genetic syndrome. Previous studies by our group have established that nonsyndromic, familial TAAD is inherited in an autosomal dominant manner with decreased penetrance and variable expression. Using one large family with multiple members with TAAD for the genome wide scan, a major locus for familial TAAD was mapped to 5q13–14 (TAAD1). Nine out of 15 families studied were linked to this locus, establishing that TAAD1 was a major locus, and that there was genetic heterogeneity for the condition. Mapping of TAAD2 locus was accomplished using a single large family with multiple members with TAAD not linked to known loci of aneurysm formation. This established a second novel locus for familial TAAD on 3p24–25 (LOD score of 4.3), termed the TAAD2 locus. Two putative loci with suggestive LOD scores were mapped on 4q and 12q through a genome scan carried out using three families. TAAD phenotype in 12 families did not segregate with known loci, indicating further genetic heterogeneity. An STS-tagged BAC based contig was constructed for 7.8Mb and 25Mb critical interval of TAAD1 and TAAD2 respectively and characterized to identify the defective gene. The hypothesis that the defective genes responsible for the TAAD1 and TAAD2 encoded extracellular matrix (ECM) proteins, the major components of the elastic fiber system in the aortic media was tested. Four genes encoding ECM proteins, versican, thrombospondin-3, CRTL1, on TAAD1 and FBLN2 at TAAD2 were sequenced, but no disease-causing mutations were identified. Studies to identify the defective gene are initiated through the positional candidate gene approach using combination of bioinformatics and expression studies. The identification of the TAAD susceptibility genes will allow for presymptomatic diagnosis of individuals at risk for this life threatening disease. The identification of the molecular defects that contribute to TAAD will also further our understanding of the proteins that provide structural integrity to the aortic wall. ^

The genetic contribution of the major histocompatibility complex to bleomycin-induced lung injury

Pulmonary fibrosis (PF) is the result of a variety of environmental and cancer treatment related insults and is characterized by excessive deposition of collagen. Gas exchange in the alveoli is impaired as the normal lung becomes dense and collapsed leading to a loss of lung volume. It is now accepted that lung injury and fibrosis are in part genetically regulated. ^ Bleomycin is a chemotherapeutic agent used for testicular cancer and lymphomas that induces significant pulmonary toxicity. We delivered bleomycin to mice subcutaneously via a miniosmotic pump in order to elicit lung injury (LI) and quantified the %LI morphometrically using video imaging software. We previously identified a quantitative trait loci, Blmpf-1(LOD=17.4), in the Major Histocompatibility Complex (MHC), but the exact genetic components involved have remained unknown. ^ In the current studies, Blmpf-1 was narrowed to an interval spanning 31.9-32.9Mb on Chromosome 17 using MHC Congenic mice. This region includes the MHC Class II and III genes, and is flanked by the TNF-alpha super locus and MHC Class I genes. Knockout mice of MHC Class I genes (B2mko), MHC Class II genes (Cl2ko), and TNF-alpha (TNF-/-) and its receptors (p55-/-, p75-/-, and p55/p75-/-) were treated with bleomycin in order to ascertain the role of these genes in the pathogenesis of lung injury. ^ Cl2ko mice had significantly better survival and %LI when compared to treated background BL/6 (B6, P<.05). In contrast, B2mko showed no differences in survival or %LI compared to B6. This suggests that the MHC Class II locus contains susceptibility genes for bleomycin-induced lung injury. ^ TNF-alpha, a Class III gene, was examined and it was found that TNF-/- and p55-/- mice had higher %LI and lower survival when compared to B6 (P<.05). In contrast, p75-/- mice had significantly reduced %LI when compared to TNF-/-, p55-/-, and B6 mice as well as higher survival (P<.01). These data contradict the current paradigm that TNF-alpha is a profibrotic mediator of lung injury and suggest a novel and distinct role for the p55 and p75 receptors in mediating lung injury. ^

Identification of genetic variation influencing apolipoprotein E levels: Follow-up of genome-wide linkage scans in the GENOA study

Apolipoprotein E (ApoE) plays a major role in the metabolism of high density and low density lipoproteins (HDL and LDL). Its common protein isoforms (E2, E3, E4) are risk factors for coronary artery disease (CAD) and explain between 16 to 23% of the inter-individual variation in plasma apoE levels. Linkage analysis has been completed for plasma apoE levels in the GENOA study (Genetic Epidemiology Network of Atherosclerosis). After stratification of the population by lipoprotein levels and body mass index (BMI) to create more homogeneity with regard to biological context for apoE levels, Hispanic families showed significant linkage on chromosome 17q for two strata (LOD=2.93 at 104 cM for a low cholesterol group, LOD=3.04 at 111 cM for a low cholesterol, high HDLC group). Replication of 17q linkage was observed for apoB and apoE levels in the unstratified Hispanic and African-American populations, and for apoE levels in African-American families. Replication of this 17q linkage in different populations and strata provides strong support for the presence of gene(s) in this region with significant roles in the determination of inter-individual variation in plasma apoE levels. Through a positional and functional candidate gene approach, ten genes were identified in the 17q linked region, and 62 polymorphisms in these genes were genotyped in the GENOA families. Association analysis was performed with FBAT, GEE, and variance-component based tests followed by conditional linkage analysis. Association studies with partial coverage of TagSNPs in the gene coding for apolipoprotein H (APOH) were performed, and significant results were found for 2 SNPs (APOH_20951 and APOH_05407) in the Hispanic low cholesterol strata accounting for 3.49% of the inter-individual variation in plasma apoE levels. Among the other candidate genes, we identified a haplotype block in the ACE1 gene that contains two major haplotypes associated with apoE levels as well as total cholesterol, apoB and LDLC levels in the unstratified Hispanic population. Identifying genes responsible for the remaining 60% of inter-individual variation in plasma apoE level, will yield new insights into the understanding of genetic interactions involved in the lipid metabolism, and a more precise understanding of the risk factors leading to CAD. ^

Gene by smoking interaction in blood pressure

Hypertension (HT) is mediated by the interaction of many genetic and environmental factors. Previous genome-wide linkage analysis studies have found many loci that show linkage to HT or blood pressure (BP) regulation, but the results were generally inconsistent. Gene by environment interaction is among the reasons that potentially explain these inconsistencies between studies. Here we investigate influences of gene by smoking (GxS) interaction on HT and BP in European American (EA), African American (AA) and Mexican American (MA) families from the GENOA study. A variance component-based method was utilized to perform genome-wide linkage analysis of systolic blood pressure (SBP), diastolic blood pressure (DBP), and HT status, as well as bivariate analysis for SBP and DBP for smokers, non-smokers, and combined groups. The most significant results were found for SBP in MA. The strongest signal was for chromosome 17q24 (LOD = 4.2), increased to (LOD = 4.7) in bivariate analysis but there was no evidence of GxS interaction at this locus (p = 0.48). Two signals were identified only in one group: on chromosome 15q26.2 (LOD = 3.37) in non-smokers and chromosome 7q21.11 (LOD = 1.4) in smokers, both of which had strong evidence for GxS interaction (p = 0.00039 and 0.009 respectively). There were also two other signals, one on chromosome 20q12 (LOD = 2.45) in smokers, which became much higher in the combined sample (LOD = 3.53), and one on chromosome 6p22.2 (LOD = 2.06) in non-smokers. Neither peak had very strong evidence for GxS interaction (p = 0.08 and 0.06 respectively). A fine mapping association study was performed using 200 SNPs in 30 genes located under the linkage signals on chromosomes 15 and 17. Under the chromosome 15 peak, the association analysis identified 6 SNPs accounting for a 7 mmHg increase in SBP in MA non-smokers. For the chromosome 17 linkage peak, the association analysis identified 3 SNPs accounting for a 6 mmHg increase in SBP in MA. However, none of these SNPs was significant after correcting for multiple testing, and accounting for them in the linkage analysis produced very small reductions in the linkage signal. ^ The linkage analysis of BP traits considering the smoking status produced very interesting signals for SBP in the MA population. The fine mapping association analysis gave some insight into the contribution of some SNPs to two of the identified signals, but since these SNPs did not remain significant after multiple testing correction and did not explain the linkage peaks, more work is needed to confirm these exploratory results and identify the culprit variations under these linkage peaks. ^

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

FZD6, MATN2 AND SLC25A32, POSSIBLE CANDIDATE GENES IN NONSYNDROMIC CLEFT LIP AND PALATE

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect with a multifactorial etiology. Despite decades of research, the genetic underpinnings of NSCLP still remain largely unexplained. A genome wide association study (GWAS) of a large NSCLP African American family with seven affected individuals across three generations found evidence for linkage at 8q21.3-24.12 (LOD = 2.98). This region contained three biologically relevant candidate genes: Frizzled-6 (FZD6) (LOD = 2.8), Matrilin-2 (MATN2) (LOD = 2.3), and Solute Carrier Family 25, Member 32 (SLC26A32) (LOD = 1.6). Sequencing of the coding regions and the 5’ and 3’ UTRs of these genes in two affected family members identified a rare intronic variant, rs138557689 (c.-153+432A>C), in FZD6. The rs138557689/C allele segregated with the NSCLP phenotype; in silico analysis predicted and EMSA analysis showed that the 138557689/C allele creates new DNA binding sites. FZD6 is part of the WNT pathway, which is involved in craniofacial development, including midface development and upper lip fusion. Our novel findings suggest that an alteration in FZD6 gene regulation may perturb this tightly controlled biological pathway and in turn contribute to the development of NSCLP in this family. Studies are underway to further define how the rs138557689/C variant affects expression of FZD6.

Plio-Pleistocene diatom biostratigraphy from ODP Leg 177, Atlantic sector of the Southern Ocean

Seven sites were drilled during Ocean Drilling Program Leg 177 in the Atlantic sector of the Southern Ocean (SO) on a transect over the Antarctic Circumpolar Current from the Subantarctic to the Antarctic Zone. At four sites sediments were recovered with a Pliocene/Pleistocene sediment package of up to 580 m allowing the refinement of previous diatom zonation concepts. Samples were analyzed on stratigraphic distribution and abundance of diatom species. A refined diatom biozonation tied to the geomagnetic polarity record is proposed. For the middle and late Pleistocene two zonations applicable to the northern and southern area of the SO were constructed, considering different latitudinal distributions of biostratigraphic diatom marker species. The southern zonation for the Pleistocene relies on the occurrence of species of the genus Rouxia, R. leventerae and R. constricta n. sp. as well as on a revised last occurrence datum (LOD) of Actinocyclus ingens (0.38 Ma, late marine isotope stage (MIS) 11). The use of these new stratigraphic marker species refines the temporal resolution for biostratigraphic age assignment to up to 0.1 Myr. In particular the LOD of R. leventerae as an indicator for the MIS 6/5 boundary (Termination II) will improve future dating of carbonate-free Antarctic sediments. These new data were obtained from sediments of Sites 1093 and 1094 (Antarctic Zone). The northern zonation for the middle and late Pleistocene time interval is based on the Pleistocene abundance pattern of Hemidiscus karstenii which was already proposed by previous investigations (e.g. Gersonde and Barcena, 1998). One new species (R. constricta) and two new combinations (Fragilariopsis clementia, Fragilariopsis reinholdii) are proposed in this study.

A refined Pliocene to early Pleistocene chronostratigraphic frame at ODP Hole 653A (West Mediterranean)

Two planktonic foraminiferal oxygen isotope records of ODP Hole 653A (Tyrrhenian Sea) are presented for the time period extending from approximately 0.8 to 3.0 Ma. Six, generally accepted, synchronous bioevents were used to precise the oxygen isotope chronology and to identify the oxygen isotope stages 22 down to 114. Subsequently, this oxygen isotope chronology was used to determine the synchronism or diachronism of various other biostratigraphic events with those recorded in the Singa and Ficarazzi land sections (Italy) and those in other DSDP/ODP sites. New results concern the diachronity of the FOD of the planktonic foraminiferal species N. atlantica, G.truncatulinoides truncatulinoides and G. inflata between ODP Hole 653A and the Italian landsections. Because many species entered the Mediterranean in short term fluxes, strongly related to the southward migration of cool North Atlantic surface waters, their time distribution through the Pliocene-Pleistocene generally corresponds to alternated intervals of presence and absence. This should explain most of the apparently diachronous appearances and disappearances. Alternating presence-absence patterns are of less importance for the various nannofossil events. The LOD of D. surculus occurs during the transition of stage 100 to 101 in both ODP Hole 653A and the Singa section, which is in perfect agreement with the disappearance of this species from the open ocean. The LOD of D. pentaradiatus in the Mediterranean occurs in stages 100-99, which seems to be consistent with the extinction of this species in the southern Hemisphere. G. oceanica, which corresponds to the 4 µm < Gephyrocapsa spp <5.5 µm is recorded in stages 65 to 64 at ODP Hole 653A. The Gephyrocapsa spp. >5.5 µm first occurred in stage 51 at Hole 653A, which fits within the uncertainty interval for this event stretching from stage 51 to 47 in the open ocean and seems therefore a useful tool for conventional biostratigraphy in the Mediterranean.

Calcareous nannofossils in sediments of the Japan Sea

Preliminary results of the biostratigraphic analysis of calcareous nannofossils recovered from Ocean Drilling Program Leg 128, Sites 798 and 799, provide clues to the Quaternary oceanography of the Japan Sea. The distribution of calcareous nannofossils from the Quaternary sediments at Site 798 (903 m water depth) may record the position of an Oceanographic frontal boundary between warm water derived from a branch of the Kuroshio Current as it entered the Japan Sea through the Tsushima Straits to the south, and colder water introduced into the western portion of the Japan Sea derived from the winter chilling of northern Japan Sea surface waters. This Oceanographic front probably oscillated north-south over Site 798 in response to glacial/interglacial cycles, or perhaps to some other climatic event or combination of events unique to the Japan Sea. During the last 1.5 m.y., six major intervals are recognized when the Oceanographic front may have been north of Site 798 separated by five major intervals when the frontal boundary may have been south of the site. These migrations were centered around approximately 0.125, 0.29, 0.56, 0.62, 0.85, 0.91, 0.98, 1.0, 1.11, and 1.5 Ma, which correspond to the boundaries separating nannofossil-rich sediments from barren or nearly barren, low-carbonate intervals. Nannofossil-rich intervals may represent times when the frontal boundary was north of Site 798, and the site was above the CCD. Barren or nearly barren intervals represent times when the frontal boundary may have been south of Site 798 and the CCD was probably higher. The distribution of calcareous nannofossils at Site 799 (2073 m water depth) appears to be controlled more by the depth of the CCD than by any climatic effects. The FOD (first occurrence datum) of Emiliania huxleyi, the LOD (last occurrence datum) of Psuedoemiliania lacunosa, Helicosphaera sellii, Calcidiscus macintyrei (10 ?m), and the FOD and LOD of Reticulofenestra asanoi are recognized from Site 798 cores. The LOD of P. lacunosa is observed in sediments from Site 799. Only in the sediments younger than 1.5 Ma are the nannofossils from Sites 798 and 799 preserved well enough and sufficiently numerous for age dating and paleoceanographic conjecture. In-situ dissolution in older sediments at both sites precludes any dating or paleoenvironmental interpretations.

Age model of ODP Hole 107-654A

It has long been speculated that glacio-eustatic sea level drop may have caused or contributed to the isolation and consequent desiccation of the Mediterranean in the late Miocene (the 'Messinian salinity crisis'). Ocean Drilling Program site 654 on the Sardinia margin of the Tyrrhenian Sea is the first deep-sea drill site to penetrate through upper Messinian evaporites into lower Messinian/upper Tortonian open marine sediments, and thus offers a unique opportunity to date the onset of the salinity crisis. A reexamination of the magnetostratigraphic, biostratigraphic, and stable isotope-stratigraphic constraints on the preevaporite sediments of site 654 has yielded two possible ages for the contact between salinity crisis sediments and the underlying normal marine sediments. One magnetostratigraphic interpretation plus the biostratigraphically determined position of the Tortonian/Messinian boundary imply a date of about 6.2 Ma for the youngest presalinity crisis sediments. An alternative magnetostratigraphic interpretation plus the carbon isotope stratigraphy imply a date of about 5.2 Ma. The younger of these dates coincides with a delta18O spike in open ocean sediments [Keigwin, 1987 doi:10.1029/PA002i006p00639], which is attributed to increased ice volume.

Oxygen and carbon isotope stratigraphy of ODP Holes 107-653A and 107-654A

Stable isotope analysis of two species (or groups of species) of planktonic foraminifers: Globigerinoides ruber (or G. obliquus and G. obliquus extremus) and Globigerina bulloides (or G. falconensis and G. obesa) from ODP Hole 653A and Site 654 in the Tyrrhenian basin, records the Pliocene-Pleistocene glacial history of the Northern Hemisphere. The overall increase in mean d18O values through the interval 4.6-0.08 Ma is 1.7 per mil for G. bulloides and 1.5 per mil for G. ruber. The time interval 3.1-2.5 Ma corresponds to an important phase of 18O enrichment for planktonic foraminifers. In this interval, glacial d18O values of both species G. bulloides and G. ruber increase by about l per mil, this increase being more progressive for G. ruber than for G. bulloides. The increase of interglacial d18O values is higher for G. bulloides (1.5 per mil) than for the Gruber group (1 per mil). These data suggest a more pronounced seasonal stratification of the water masses during interglacial phases. Large positive d18O fluctuations of increasing magnitude are also recorded at 2.25 and 2.15 Ma by G bulloides and appear to be diachronous with those of Site 606 in the Atlantic Ocean. Other events of increasing d18O values are recorded between 1.55 and 1.3 Ma, at 0.9 Ma, 0.8 Ma, and near 0.34 Ma. In the early Pliocene the d18O variability recorded by the planktonic species G. bulloides was higher in the Mediterranean than in the Atlantic at the same latitude. This suggests that important cyclic variations in the water budget of the Mediterranean occurred since that time. Step increases in the d18O variability are synchronous with those of the open ocean at 0.9 and 0.34 Ma. The higher variability as well as the higher amplitude of the peaks of 18O enrichment may be partly accounted for by increase of dryness over the Mediterranean area. In particular the high amplitude d18O fluctuations recorded between 3.1 and 2.1 Ma are correlated with the onset of a marked seasonal contrast and a summer dryness, revealed by pollen analyses. Strong fluctuations towards d13C values higher than modern ones are recorded by the G. ruber group species before 1.7 Ma and suggest a high production of phytoplankton. When such episodes of high primary production are correlated with episodes of decreasing 13C content of G. bulloides, they are interpreted as the consequence of a higher stratification of the upper water masses resulting itself from a marked seasonality. Such episodes occur between 4.6 and 4.05 Ma, 3.9 and 3.6 Ma, and 3.25 and 2.66 Ma. The interval 2.66-1.65 Ma corresponds to a weakening of the stratification of the upper water layers. This may be related to episodes of cooling and increasing dryness induced by the Northern Hemisphere Glaciations. The Pleistocene may have been a less productive period. The transition from highly productive to less productive surface waters also coincides with a new step increase in dryness and cooling, between 1.5 and 1.3 Ma. The comparison of the 13C records of G ruber and G. bulloides in fact suggests that a high vertical convection became a dominant feature after 2.6 Ma. Increases in the nutrient input and the stratification of the upper water masses may be suspected, however, during short episodes near 0.86 Ma (isotopic stage 25), 0.57-0.59 Ma (isotopic stage 16), 0.49 Ma (isotopic stage 13), 0.4-0.43 Ma (isotopic stage 11), and 0.22 and 0.26 Ma (part of isotopic stage 7 and transition 7/8). In fact, changes in the C02 balance within the different water masses of the Tyrrhenian basin as well as in the local primary production did not follow the general patterns of the open ocean.