Molecular cytogenetics and gene mapping in sheep (Ovis aries, 2n = 54)

The development of a completely annotated sheep genome sequence is a key need for understanding the phylogenetic relationships and genetic diversity among the many different sheep breeds worldwide and for identifying genes controlling economically and physiologically important traits. The ovine genome sequence assembly will be crucial for developing optimized breeding programs based on highly productive, healthy sheep phenotypes that are adapted to modern breeding and production conditions. Scientists and breeders around the globe have been contributing to this goal by generating genomic and cDNA libraries, performing genome-wide and trait-associated analyses of polymorphism, expression analysis, genome sequencing, and by developing virtual and physical comparative maps. The International Sheep Genomics Consortium (ISGC), an informal network of sheep genomics researchers, is playing a major role in coordinating many of these activities. In addition to serving as an essential tool for monitoring chromosome abnormalities in specific sheep populations, ovine molecular cytogenetics provides physical anchors which link and order genome regions, such as sequence contigs, genes and polymorphic DNA markers to ovine chromosomes. Likewise, molecular cytogenetics can contribute to the process of defining evolutionary breakpoints between related species. The selective expansion of the sheep cytogenetic map, using loci to connect maps and identify chromosome bands, can substantially contribute to improving the quality of the annotated sheep genome sequence and will also accelerate its assembly. Furthermore, identifying major morphological chromosome anomalies and micro-rearrangements, such as gene duplications or deletions, that might occur between different sheep breeds and other Ovis species will also be important to understand the diversity of sheep chromosome structure and its implications for cross-breeding. To date, 566 loci have been assigned to specific chromosome regions in sheep and the new cytogenetic map is presented as part of this review. This review will also summarize the current cytogenomic status of the sheep genome, describe current activities in the sheep cytogenomics research sector, and will discuss the cytogenomics data in context with other major sheep genomics projects.

Ar diffusion in K-feldspar: present and absent

The mechanisms of Ar release from K-feldspar samples in laboratory experiments and during their geological history are assessed here. Modern petrology clearly established that the chemical and isotopic record of minerals is normally dominated by aqueous recrystallization. The laboratory critique is trickier, which explains why so many conflicting approaches have been able to survive long past their expiration date. Current models are evaluated for self-consistency; especially Arrhenian non-linearity leads to paradoxes. The models’ testable geological predictions suggest that temperature-based downslope extrapolations often overestimate observed geological Ar mobility substantially. An updated interpretation is based on the unrelatedness of geological behaviour to laboratory experiments. The isotopic record of K-feldspar in geological samples is not a unique function of temperature, as recrystallisation promoted by aqueous fluids is the predominant mechanism controlling isotope transport. K-feldspar should therefore be viewed as a hygrochronometer. Laboratory degassing proceeds from structural rearrangements and phase transitions such as are observed in situ at high temperature in Na and Pb feldspars. These effects violate the mathematics of an inert Fick’s Law matrix and preclude downslope extrapolation. The similar upward-concave, non-linear shapes of Arrhenius trajectories of many silicates, hydrous and anhydrous, are likely common manifestations of structural rearrangements in silicate structures.

High-resolution sequence stratigraphic correlation in the Upper Jurassic (Kimmeridgian)–Upper Cretaceous (Cenomanian) peritidal carbonate deposits (Western Taurides, Turkey)

Upper Jurassic (Kimmeridgian)±Upper Cretaceous (Cenomanian) inner platform carbonates in the Western Taurides are composed of metre-scale upward-shallowing cyclic deposits (parasequences) and important karstic surfaces capping some of the cycles. Peritidal cycles (shallow subtidal facies capped by tidal-¯at laminites or fenestrate limestones) are regressive- and transgressive-prone (upward-deepening followed by upward-shallowing facies trends). Subtidal cycles are of two types and indicate incomplete shallowing. Submerged subtidal cycles are composed of deeper subtidal facies overlain by shallow subtidal facies. Exposed subtidal cycles consist of deeper subtidal facies overlain by shallow subtidal facies that are capped by features indicative of prolonged subaerial exposure. Subtidal facies occur characteristically in the Jurassic, while peritidal cycles are typical for the Lower Cretaceous of the region. Within the foraminiferal and dasyclad algal biostratigraphic framework, four karst breccia levels are recognized as the boundaries of major second-order cycles, introduced for the ®rst time in this study. These levels correspond to the Kimmeridgian±Portlandian boundary, mid-Early Valanginian, mid-Early Aptian and mid-Cenomanian and represent important sea level falls which affected the distribution of foraminiferal fauna and dasyclad ¯ora of the Taurus carbonate platform. Within the Kimmeridgian±Cenomanian interval 26 third-order sequences (types 1 and 2) are recognized. These sequences are the records of eustatic sea level ¯uctuations rather than the records of local tectonic events because the boundaries of the sequences representing 1±4 Ma intervals are correlative with global sea level falls. Third-order sequences and metre-scale cyclic deposits are the major units used for long-distance, high-resolution sequence stratigraphic correlation in the Western Taurides. Metre-scale cyclic deposits (parasequences) in the Cretaceous show genetical stacking patterns within third-order sequences and correspond to fourth-order sequences representing 100±200 ka. These cycles are possibly the E2 signal (126 ka) of the orbital eccentricity cycles of the Milankovitch band. The slight deviation of values, calculated for parasequences, from the mean value of eccentricity cycles can be explained by the currently imprecise geochronology established in the Cretaceous and missed sea level oscillations when the platform lay above fluctuating sea level.

The effect of idursulfase on growth in patients with Hunter syndrome: data from the Hunter Outcome Survey (HOS)

Hunter syndrome (mucopolysaccharidosis type II) is a rare and life-limiting multisystemic disorder with an X-linked recessive pattern of inheritance. Short stature is a prominent feature of this condition. This analysis aimed to investigate the effects of enzyme replacement therapy with idursulfase on growth in patients enrolled in HOS - the Hunter Outcome Survey which is a multinational observational database. As of Jan 2012, height data before treatment were available for 567 of 740 males followed prospectively after HOS entry. Cross-sectional analysis showed that short stature became apparent after approximately 8 years of age; before this, height remained within the normal range. Age-corrected standardized height scores (z-scores) before and after treatment were assessed using piecewise regression model analysis in 133 patients (8-15 years of age at treatment start; data available on ≥ 1 occasion within +/-24 months of treatment start; growth hormone-treated patients excluded). Results showed that the slope after treatment (slope=-0.005) was significantly improved compared with before treatment (slope=-0.043) (difference=0.038, p=0.004). Analysis of covariates (age at treatment start, cognitive involvement, presence of puberty at the start of ERT, mutation type, functional classification), showed a significant influence on growth of mutation type (height deficit in terms of z-scores most pronounced in patients with deletions/large rearrangements/nonsense mutations, p<0.0001) and age (most pronounced in the 12-15-year group, p<0.0001). Cognitive involvement, pubertal status at the start of ERT and functional classification were not related to the growth deficit or response to treatment. In conclusion, the data showed an improvement in growth rate in patients with Hunter syndrome following idursulfase treatment.

Does Sex Speed Up Evolutionary Rate and Increase Biodiversity?

Most empirical and theoretical studies have shown that sex increases the rate of evolution, although evidence of sex constraining genomic and epigenetic variation and slowing down evolution also exists. Faster rates with sex have been attributed to new gene combinations, removal of deleterious mutations, and adaptation to heterogeneous environments. Slower rates with sex have been attributed to removal of major genetic rearrangements, the cost of finding a mate, vulnerability to predation, and exposure to sexually transmitted diseases. Whether sex speeds or slows evolution, the connection between reproductive mode, the evolutionary rate, and species diversity remains largely unexplored. Here we present a spatially explicit model of ecological and evolutionary dynamics based on DNA sequence change to study the connection between mutation, speciation, and the resulting biodiversity in sexual and asexual populations. We show that faster speciation can decrease the abundance of newly formed species and thus decrease long-term biodiversity. In this way, sex can reduce diversity relative to asexual populations, because it leads to a higher rate of production of new species, but with lower abundances. Our results show that reproductive mode and the mechanisms underlying it can alter the link between mutation, evolutionary rate, speciation and biodiversity and we suggest that a high rate of evolution may not be required to yield high biodiversity.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor "silent" germline micronuclear genome by a process of "unscrambling" and fragmentation. The tiny macronuclear "nanochromosomes" typically encode single, protein-coding genes (a small portion, 10%, encode 2-8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

Mutagenicity of herpes simplex virus type 1 in a shuttle vector

The shuttle vector plasmid pZ189 was used to find the kinds of mutations that are induced by herpes simplex virus type-1 (HSV-1). In cells infected by HSV-1 the frequency of mutation in supF gene, the mutagenesis marker, was increased over background by from two- to seven-fold, reaching 0.14-0.45%. No increase was induced by infection by vaccinia virus under the same conditions. Mutagenesis was an early event, showing a four-fold increase in mutation frequency at only two hours after infection, and peaking at a seven-fold increase at four hours after infection. DNA sequencing and gel electrophoresis analysis were performed on 105 HSV-1 induced mutants and 65 spontaneous mutants and provided the following information: (1) A change in plasmid size was seen in 54% of HSV-1 related mutants, compared with only 37% of spontaneous mutants. (2) Among point mutations, the predominant type was G:C to A:T transition, which accounted for 51% of point mutations in mutants isolated from cells infected with HSV-1, and 32% of point mutations in spontaneous mutants. (3) Deletions of DNA were seen in HSV-1 related mutants at a frequency of 40%, compared with 29% in spontaneous mutants. The HSV-1 related deletions were about half the length of spontaneous mutants and three contained short filler sequences. (4) Fifteen (15%) of HSV-1 induced mutants revealed the altered restriction patterns on agarose gel electrophoresis analysis and were due either to rearrangements of plasmid DNA, and/or to insertion of sequences derived from chromosomal DNA (seven plasmids). No insertions of DNA from HSV-1 were detected. Among spontaneous mutants, only 5 (7.7%) were rearrangements and none had inserted chromosomal DNA. (5) DNA sequence analysis of seven plasmids with inserted chromosomal DNA revealed that four cases had repetitive DNA sequences integrated and the other three were unidentified sequences from the GenBank database. Three repetitive DNA included $\alpha$ satellite, Alu and KpnI family sequences. The other sequence was identified as tRNA-like component. The observed mutations have implications for the mechanism of malignant transformation of cells by HSV-1. ^

Regulation of Spec genes in aboral ectoderm cells of sea urchin {\it Strongylocentrotus purpuratus}

The Spec genes of the sea urchin Stronylocentrotus purpuratus serves as an excellent model for studying cell type-specific gene expression during early embryogenesis. The Spec1/Spec2 genes encode cytosolic calcium-binding proteins related to the calmodulin/troponin C/myosin light chain superfamily. Members of the Spec gene family are activated shortly after the sixth cleavage as the lineage-specific founder cells giving rise to aboral ectoderm are established, and the accumulation of the Spec mRNAs is limited exclusively to aboral ectoderm cell lineages. In this dissertation, the transcriptional regulation of the Spec genes was studied. Sequence comparisons of the Spec gene 5$\sp\prime$ flanking regions showed that a DNA block of approximately 800 bp from the 3$\sp\prime$ end of the first exon to the 5$\sp\prime$ end of a repetitive DNA element, termed RSR, was highly conserved. In Spec2a, the conserved region was a continuous stretch of DNA, but in Spec1 and Spec2c, DNA insertions interrupt the conserved sequence block and alter the relative placement of the RSR element and other 5$\sp\prime$ flanking DNA. Thus, drastic rearrangements have occurred within the putative control regions of the Spec genes. In vivo expression experiments using the sea urchin embryo gene-transfer system showed that while the 5$\sp\prime$ flanking regions of all three Spec genes conferred proper temporal activation to the reporter CAT gene, only the Spec2a 5$\sp\prime$ flanking region could restrict lacZ gene expression to aboral ectoderm cells. However, the Spec2a conserved region alone was not sufficient to confer proper spatial expression, suggesting that negative spatial elements are also associated with the proper activation of Spec2a. A major positive regulatory region, defined as the RSR enhancer, was identified between base pairs $-$631 and $-$443 on Spec2a. The RSR enhancer was essential for maximal activity and conferred preferential aboral ectoderm expression to a lacZ reporter gene. DNaseI footprinting and band-shift analysis of the RSR enhancer revealed multiple DNA-elements. One of the elements, an A/T-rich sequence called the A/T palindrome was studied in detail. This element binds a single 45-kDa nuclear protein, the A/T palindrome binding protein (A/TBP), whose DNA-binding specificity suggests a possible relationship with the bicoid-class homeodomain proteins. Mutated A/T palindromes are incapable of binding the 45-kDa protein and lower promoter activity by 8-fold. DNA-binding activity for A/TBP is low in unfertilized eggs, increases by the 16-cell stage and continues rising in blastulae. These data suggest that A/TBP plays a major role in the activation of the Spec2a gene in aboral ectoderm cells. ^

The physical map of {\it Escherichia coli\/} K-12 strain MG1655: Construction, comparison with other strains and implications for mechanisms controlling oligonucleotide frequency

Complete NotI, SfiI, XbaI and BlnI cleavage maps of Escherichia coli K-12 strain MG1655 were constructed. Techniques used included: CHEF pulsed field gel electrophoresis; transposon mutagenesis; fragment hybridization to the ordered $\lambda$ library of Kohara et al.; fragment and cosmid hybridization to Southern blots; correlation of fragments and cleavage sites with EcoMap, a sequence-modified version of the genomic restriction map of Kohara et al.; and correlation of cleavage sites with DNA sequence databases. In all, 105 restriction sites were mapped and correlated with the EcoMap coordinate system.^ NotI, SfiI, XbaI and BlnI restriction patterns of five commonly used E. coli K-12 strains were compared to those of MG1655. The variability between strains, some of which are separated by numerous steps of mutagenic treatment, is readily detectable by pulsed-field gel electrophoresis. A model is presented to account for the difference between the strains on the basis of simple insertions, deletions, and in one case an inversion. Insertions and deletions ranged in size from 1 kb to 86 kb. Several of the larger features have previously been characterized and some of the smaller rearrangements can potentially account for previously reported genetic features of these strains.^ Some aspects of the frequency and distribution of NotI, SfiI, XbaI and BlnI cleavage sites were analyzed using a method based on Markov chain theory. Overlaps of Dam and Dcm methylase sites with XbaI and SfiI cleavage sites were examined. The one XbaI-Dam overlap in the database is in accord with the expected frequency of this overlap. The occurrence of certain types of SfiI-Dcm overlaps are overrepresented. Of the four subtypes of SfiI-Dcm overlap, only one has a partial inhibitory effect on the activity of SfiI. Recognition sites for all four enzymes are rarer than expected based on oligonucleotide frequency data, with this effect being much stronger for XbaI and BlnI than for NotI and SfiI. The latter two enzyme sites are rare mainly due to apparent negative selection against GGCC (both) and CGGCCG (NotI). The former two enzyme sites are rare mainly due to effects of the VSP repair system on certain di-tri- and tetranucleotides, most notably CTAG. Models are proposed to explain several of the anomalies of oligonucleotide distribution in E. coli, and the biological significance of the systems that produce these anomalies is discussed. ^

Defining the regions of homology between human chromosome 16p12--13 and mouse chromosomes

In this study, the evolutionary relationship between human chromosome 16p12-p13 and mouse chromosomes was investigated by determining the order of marker loci in the region and then identifying the chromosomal locations of the homologous loci in mice. Eighteen genes from human 16 were mapped to fifteen subchromosomal regions by a variety of mapping approaches.^ Thirteen of the genes were mapped in the mouse. Linkage analysis with backcross mice and segregation analysis in a mouse - Chinese Hamster Ovary (CHO) somatic cell hybrid panel informative for different regions of mouse genome were used. The results assigned the thirteen genes to three different mouse chromosomes.^ A group of six genes on mouse 16 was found to be closely linked to Scid. The order of Myh11 and Mrp remains ambiguous since no recombination was detected in backcross analysis. Their relative position in human is also uncertain since they were shown to be very close to each other. For the other mouse loci, an unambiguous gene order could be determined and was found to be identical to that in human. Therefore, they comprise a new conserved linkage group between the two species. The orientation of the group was inverted relative to the centromeres, i.e. the proximal loci in one species become distal in another. The size of the group was estimated to be from 4.4 to 8 Mb and 10 to 32 cM in human. In mouse, it was about 21 cM in the backcross analysis. The two boundaries of the conserved linkage were defined within a 1 Mb range. It is now possible to predict the locations of mouse homologs for some human disease genes based on their locations on human 16p.^ The six human 16p genes that map to MMU7 showed a different gene order in mouse than in human. No recombination was found between Crym and Umod while Crym was distal to D16S79A and proximal to D16S92. The location of Stp and Cdr2 with respect to the above four loci was not determined since they were not mapped in the same set of backcross mice. These genes greatly expanded an existing conserved synteny group between the human 16p12-p13 region and the MMU7. It now consists of eleven loci that span a region of probably more than 10 Mb in human. The gene order derived from this study provided further evidence for chromosomal rearrangements within the conserved synteny. (Abstract shortened by UMI.) ^

Molecular analysis of the human $\gamma$-chain T cell receptor genes

In our studies we have focused on the issue of variability and diversity of the $\gamma$ (or $\delta)$ chain T cell receptor (TCR) genes by studying cDNA transcripts in peripheral blood mononuclear cells or $\gamma\delta$ TCR+ T cell clones. The significance of these studies lies in the better understanding of the molecular biology of the $\gamma\delta$ T cell receptor as well as in answering the question whether certain molecular forms predominate in $\gamma\delta$ T cells exhibiting specific immunologic functions. We establish that certain $\gamma$-chain TCR genes exhibit particular patterns of rearrangements in cDNA transcripts in normal individuals. V$\gamma$I subgroup were shown to preferentially rearrange to J$\gamma$2C$\gamma$2 gene segments. These preferential VJC rearrangements, may have implications regarding the potential for diversity and polymorphism of the $\gamma$-chain TCR gene. In addition, the preferential association of V$\gamma$I genes with J$\gamma$2C$\gamma$2, which encode a non-disulfide-linked $\gamma\delta$ TCR, suggests that $\gamma$ chains utilizing V$\gamma$I are predominantly expressed as non-disulfide-linked $\gamma\delta$ TCR heterodimers. The implications of this type of expression remain to be determined. We identified two alternative splicing events of the $\gamma$-chain TCR genes occurring in high frequency in all the normal individuals examined. These events may suggest additional mechanisms of regulation and control as well as diversification of $\gamma\delta$ TCR gene expression. The question whether particular forms of $\gamma$ or $\delta$-chain TCR genes are involved in HLA Class I recognition by specific $\gamma\delta$ cytotoxic T cell clones was addressed. Our results indicated that the T cell clones expressed identical $\gamma$ but distinct $\delta$-chains suggesting that the specificity for recognition of HLA-A2 or HLA-A3 may be conferred by the $\delta$-chain TCR. The issue of the degree of diversity and polymorphism of the $\delta$-chain TCR genes in a patient with a primary immunodeficiency (Omenn's syndrome) was addressed. A limited pattern of rearrangements in peripheral blood transcripts was found, suggesting that a limited $\gamma\delta$ TCR repertoire may be expressed in this particular primary immunodeficiency syndrome. Overall, our findings suggest that $\delta$-chain TCR genes exhibit the potential for significant diversity and that there are certain preferential patterns of expression that may be associated with particular immunologic functions. ^

The expression levels of three raft-associated molecules in cultivated vascular cells are dependent on culture conditions

Relaying a signal across the plasma membrane requires functional connections between the partner molecules. Membrane microdomains or lipid rafts provide an environment in which such specific interactions can take place. The integrity of these sites is often taken for granted when signalling pathways are investigated in cell culture. However, it is well known that smooth muscle and endothelial cells undergo cytoskeletal rearrangements during monolayer culturing. Likewise affected--and with potentially important consequences for signalling events--is the organization of the plasma membrane. The expression levels of three raft markers were massively upregulated, and raft-associated 5'-nucleotidase activity increased in conventional monolayer cultures as compared with a spheroidal coculture model, shown to promote the differentiation of endothelial cells. Our data point to a shift of raft components in monolayer cultures and demonstrate potential advantages of the spheroid coculture system for investigation of raft-mediated signalling events in endothelial cells.

GENETIC EVENTS PREDISPOSING TO WILMS' TUMOR

Wilms' tumor (WT) is a childhood embryonic tumor of the kidney. In some cases, WT has been associated with a chromosome deletion in the region 11p13. The majority of WT cases, however, have normal karyotypes with no discernable deletions or rearrangements of chromosome 11.^ To study the genetic events predisposing to the development of WT, I have used a number of gene markers specific for chromosome 11. Gene probes for human catalase and apolipoprotein A1 were localized to chromosome 11 by in situ hybridization. A number of other probes previously mapped to chromosome 11 were also used. Nine WT patients who were heterozygous for at least one 11p marker were shown to lose heterozygosity in their tumor DNA. Gene dosage experiments demonstrated that two chromosomes 11 were present although loss of heterozygosity had occurred in all but two cases. By using gene probes from the short and long arms of chromosome 11, I discerned that loss of heterozygosity was due to somatic recombination in four cases, chromosome deletion in two cases, and chromosome loss and reduplication or somatic recombination in these cases. Examination of DNAs from the parents of six of these patients indicated that the alleles that were lost in tumor tissues were alleles inherited from the mother. In sporadic WT cases one would expect the loss of alleles to be random. These data suggest that the loss of alleles resulting in the development of WT is not a random event, however, the significance of this is not known. ^

The Draft Assembly of the Radically Organized Stylonychia lemnae Macronuclear Genome.

Stylonychia lemnae is a classical model single-celled eukaryote, and a quintessential ciliate typified by dimorphic nuclei: A small, germline micronucleus and a massive, vegetative macronucleus. The genome within Stylonychia's macronucleus has a very unusual architecture, comprised variably and highly amplified "nanochromosomes," each usually encoding a single gene with a minimal amount of surrounding noncoding DNA. As only a tiny fraction of the Stylonychia genes has been sequenced, and to promote research using this organism, we sequenced its macronuclear genome. We report the analysis of the 50.2-Mb draft S. lemnae macronuclear genome assembly, containing in excess of 16,000 complete nanochromosomes, assembled as less than 20,000 contigs. We found considerable conservation of fundamental genomic properties between S. lemnae and its close relative, Oxytricha trifallax, including nanochromosomal gene synteny, alternative fragmentation, and copy number. Protein domain searches in Stylonychia revealed two new telomere-binding protein homologs and the presence of linker histones. Among the diverse histone variants of S. lemnae and O. trifallax, we found divergent, coexpressed variants corresponding to four of the five core nucleosomal proteins (H1.2, H2A.6, H2B.4, and H3.7) suggesting that these ciliates may possess specialized nucleosomes involved in genome processing during nuclear differentiation. The assembly of the S. lemnae macronuclear genome demonstrates that largely complete, well-assembled highly fragmented genomes of similar size and complexity may be produced from one library and lane of Illumina HiSeq 2000 shotgun sequencing. The provision of the S. lemnae macronuclear genome sets the stage for future detailed experimental studies of chromatin-mediated, RNA-guided developmental genome rearrangements.

Influence of the Central American Seaway and Drake Passage on ocean circulation and neodymium isotopes: A model study

The sensitivity of the neodymium isotopic composition (ϵNd) to tectonic rearrangements of seaways is investigated using an Earth System Model of Intermediate Complexity. The shoaling and closure of the Central American Seaway (CAS) is simulated, as well as the opening and deepening of Drake Passage (DP). Multiple series of equilibrium simulations with various intermediate depths are performed for both seaways, providing insight into ϵNd and circulation responses to progressive throughflow evolutions. Furthermore, the sensitivity of these responses to the Atlantic Meridional Overturning Circulation (AMOC) and the neodymium boundary source is examined. Modeled ϵNd changes are compared to sediment core and ferromanganese (Fe-Mn) crust data. The model results indicate that the North Atlantic ϵNd response to the CAS shoaling is highly dependent on the AMOC state, i.e., on the AMOC strength before the shoaling to shallow depths (preclosure). Three scenarios based on different AMOC forcings are discussed, of which the model-data agreement favors a shallow preclosure (Miocene) AMOC (∼6 Sv). The DP opening causes a rather complex circulation response, resulting in an initial South Atlantic ϵNd decrease preceding a larger increase. This feature may be specific to our model setup, which induces a vigorous CAS throughflow that is strongly anticorrelated to the DP throughflow. In freshwater experiments following the DP deepening, ODP Site 1090 is mainly influenced by AMOC and DP throughflow changes, while ODP Site 689 is more strongly influenced by Southern Ocean Meridional Overturning Circulation and CAS throughflow changes. The boundary source uncertainty is largest for shallow seaways and at shallow sites.