The use of AARS activity as a proxy for zooplankton and ichthyoplankton growth rates

Programa de doctorado en Oceanografía

Spatial and temporal variation in settlement, growth and condition of the rockfish species Sebastes caurinus and Sebastes carnatus in the Channel Islands region, Santa Barbara, CA

Survival during the early life stages of marine species, including nearshore temperate reef fishes, is typically very low, and small changes in mortality rates, due to physiological and environmental conditions, can have marked effects on survival of a cohort and, on a larger scale, on the success of a recruitment season. Moreover, trade offs between larval growth and accumulation of energetic resources prior to settlement are likely to influence growth and survival until this critical period and afterwards. Rockfish recruitment rates are notoriously variable between years and across geographic locations. Monitoring of rates of onshore delivery of pelagic juveniles (defined here as settlement) of two species of nearshore rockfishes, Sebastes caurinus and Sebastes carnatus, was done between 2003-2009 years using artificial collectors placed at San Miguel and Santa Cruz Island, off Southern California coast. I investigated spatiotemporal variation in settlement rate, lipid content, pelagic larval duration and larval growth of the newly settled fishes; I assessed relationships between birth date, larval growth, early life-history characteristics and lipid content at settlement, considering also interspecific differences; finally, I attempt to relate interannual patterns of settlement and of early life history traits to easily accessible, local and regional indices of ocean conditions including in situ ocean temperature and regional upwelling, sea surface temperature (SST) and Chlorophyll-a (Chl-a) concentration. Spatial variations appeared to be of low relevance, while significant interannual differences were detected in settlement rate, pelagic larval duration and larval growth. The amount of lipid content of the newly settled fishes was highly variable in space and time, but did not differ between the two species and did not show any relationships with early life history traits, indicating that no trade off involved these physiological processes or they were masked by high individual variability in different periods of larval life. Significant interspecific differences were found in the timing of parturition and settlement and in larval growth rates, with S. carnatus growing faster and breeding and settling later than S. caurinus. The two species exhibited also different patterns of correlations between larval growth rates and larval duration. S. carnatus larval duration was longer when the growth in the first two weeks post-hatch was faster, while S. caurinus had a shorter larval duration when grew fast in the middle and in the end of larval life, suggesting different larval strategies. Fishes with longer larval durations were longer in size at settlement and exhibited longer planktonic phase in periods of favourable environmental conditions. Ocean conditions had a low explanatory power for interannual variation in early life history traits, but a very high explanatory power for settlement fluctuations, with regional upwelling strength being the principal indicator. Nonetheless, interannual variability in larval duration and growth were related to great phenological changes in upwelling happened during the period of this study and that caused negative consequences at all trophic levels along the California coast. Despite the low explanatory power of the environmental variables used in this study on the variation of larval biological traits, environmental processes were differently related with early life history characteristics analyzed to species, indicating possible species-specific susceptibility to ocean conditions and local environmental adaptation, which should be further investigated. These results have implications for understanding the processes influencing larval and juvenile survival, and consequently recruitment variability, which may be dependent on biological characteristics and environmental conditions.

Impact of ocean acidification on phytoplankton growth and aggregation

Since the industrial revolution, the ocean has absorbed around one third of the anthropogenic CO2, which induced a profound alteration of the carbonate system commonly known as ocean acidification. Since the preindustrial times, the average ocean surface water pH has fallen by 0.1 units, from approximately 8.2 to 8.1 and a further decrease of 0.4 pH units is expected for the end of the century. Despite their microscopic size, marine diatoms are bio-geo-chemically a very important group, responsible for the export of massive amount of carbon to deep waters and sediments. The knowledge of the potential effects of ocean acidification on the phytoplankton growth and on biological pump is still at its infancy. This study wants to investigate the effect of ocean acidification on the growth of the diatom Skeletonema marinoi and on its aggregation, using a mechanistic approach. The experiment consisted of two treatments (Present and Future) representing different pCO2 conditions and two sequential experimental phases. During the cell growth phase a culture of S. marinoi was inoculated into transparent bags and the effect of ocean acidification was studied on various growth parameters, including DOC and TEP production. The aggregation phase consisted in the incubation of the cultures into rolling tanks where the sinking of particles through the water column was simulated and aggregation promoted. Since few studies investigated the effect of pH on the growth of S. marinoi and none used pH ranges that are compatible with the OA scenarios, there were no baselines. I have shown here, that OA does not affect the cell growth of S. marinoi, suggesting that the physiology of this species is robust in respect to the changes in the carbonate chemistry expected for the end of the century. Furthermore, according to my results, OA does not affect the aggregation of S. marinoi in a consistent manner, suggesting that this process has a high natural variability but is not influenced by OA in a predictable way. The effect of OA was tested over a variety of factors including the number of aggregates produced, their size and sinking velocity, the algal, bacterial and TEP content. Many of these variables showed significant treatment effects but none of these were consistent between the two experiments.

Harmful algae and their potential impacts on the coastal ecosystem: growth and toxin production dynamics

The main goal of the present thesis was to study some harmful algal species which cause blooms in Italian coastal waters, leading to consequences for human health, coastal ecosystem, fishery and tourism. In particular, in the first part of this thesis the toxicity of Adriatic strains of the raphidophyte Fibrocapsa japonica was investigated. Despite several hypotheses have been proposed for the toxic mechanism of the raphidophytes, especially for the species Chattonella antiqua and C. marina, which have been studied more extensively, just a few studies on the toxic effects of these species for different organisms were reported. Moreover, a careful reading of the literature evidenced as any ichthyotoxic events reported worldwide can be linked to F. japonica blooms. Although recently several studies were performed on F. japonica strains from the USA, Japan, Australia, New Zealand, the Netherlands, Germany, and France in order to characterize their growth and toxicity features, the work reported in this thesis results one of the first investigation on the toxic effects of F. japonica for different organisms, such as bacteria, crustaceans and fish. Mortality effects, together with haemolysis of fish erythrocytes, probably due to the relatively high amount of PUFAs produced by this species, were observed. Mortality for fish, however, was reported only at a high cell density and after a long exposition period (9-10 days); moreover a significant increase of H2O2 obtained in the tanks where sea basses were exposed to F. japonica was also relevant. This result may justify the absence of ichthyotoxic events in the Italian coasts, despite F. japonica blooms detected in these areas were characterized by high cell densities. This work reports also a first complete characterization of the fatty acids produced and extracellularly released by the Adriatic F. japonica, and results were also compared with the fatty acid profile of other strains. The absence of known brevetoxins in F. japonica algal extracts was also highlighted, leading to the hypothesis that the toxicity of F. japonica may be due to a synergic effect of PUFAs and ROS. Another microalgae that was studied in this thesis is the benthic dinoflagellate Ostreopsis cf. ovata. This species was investigated with the aim to investigate the effect of environmental parameters on its growth and toxicity. O. cf. ovata, in fact, shows different blooming periods along the Italian coasts and even the reported toxic effects are variable. The results of this work confirmed the high variability in the growth dynamic and toxin content of several Italian strains which were isolated in recent years along the Adriatic and Tyrrhenian Seas. Moreover, the effects of temperature and salinity on the behaviour of the different isolates are in good agreement with the results obtained from field surveys, which evidence as the environmental parameters are important factors modulating O. cf. ovata proliferation. Another relevant result that was highlighted is the anomaly in the production of palytoxin-like compounds reported by one of the studied isolate, in particular the one isolated in 2008 in Ancona (Adriatic Sea). Only this strain reported the absence of two (ovatoxin-b and –c) of the five ovatoxins so far known in the toxin profile and a different relative abundance of the other toxins. The last aspect that was studied in this thesis regards the toxin biosythesis. In fact, toxins produced (palytoxin-like compounds) or supposed to be produced (brevetoxin-like compounds) by O. cf. ovata and F. japonica, respectively, are polyketides, which are highly oxygenated compounds synthesized by complex enzymes known as polyketide synthase (PKS) enzymes. These enzymes are multi-domain complexes that structurally and functionally resemble the fatty acid synthases (FASs). This work reports the first study of PKS proteins in the dinoflagellates O. cf. ovata, C. monotis and in the raphidophyte F. japonica. For the first time some PKSs were identified in these species, confirming the presence of PKS proteins predicted by the in silico translation of the transcripts found in K. brevis also in other species. The identification of O. cf. ovata PKSs and the localization of the palytoxin-like compounds produced by this dinoflagellate in a similar location (chloroplast) as that observed for other dinoflagellate and cyanobacterial toxins provides some indication that these proteins may be involved in polyketide biosynthesis. However, their potential function as fatty acid synthases cannot be ruled out, as plant fatty acid synthesis also occurs within chloroplasts. This last hypothesis is also supported by the fact that in all the investigated species, and in particular in F. japonica, PKS proteins were present. Therefore, these results provide an important contribution to the study of the polyketides and of the involvement of PKS proteins in the toxin biosynthesis.

Parameter estimation in a growth model for a biological population

The motivating problem concerns the estimation of the growth curve of solitary corals that follow the nonlinear Von Bertalanffy Growth Function (VBGF). The most common parameterization of the VBGF for corals is based on two parameters: the ultimate length L∞ and the growth rate k. One aim was to find a more reliable method for estimating these parameters, which can capture the influence of environmental covariates. The main issue with current methods is that they force the linearization of VBGF and neglect intra-individual variability. The idea was to use the hierarchical nonlinear model which has the appealing features of taking into account the influence of collection sites, possible intra-site measurement correlation and variance heterogeneity, and that can handle the influence of environmental factors and all the reliable information that might influence coral growth. This method was used on two databases of different solitary corals i.e. Balanophyllia europaea and Leptopsammia pruvoti, collected in six different sites in different environmental conditions, which introduced a decisive improvement in the results. Nevertheless, the theory of the energy balance in growth ascertains the linear correlation of the two parameters and the independence of the ultimate length L∞ from the influence of environmental covariates, so a further aim of the thesis was to propose a new parameterization based on the ultimate length and parameter c which explicitly describes the part of growth ascribable to site-specific conditions such as environmental factors. We explored the possibility of estimating these parameters characterizing the VBGF new parameterization via the nonlinear hierarchical model. Again there was a general improvement with respect to traditional methods. The results of the two parameterizations were similar, although a very slight improvement was observed in the new one. This is, nevertheless, more suitable from a theoretical point of view when considering environmental covariates.

High-resolution reconstruction of Holocene climate variability and environmental changes in the North Pacific using bivalve shells

Bivalve mollusk shells are useful tools for multi-species and multi-proxy paleoenvironmental reconstructions with a high temporal and spatial resolution. Past environmental conditions can be reconstructed from shell growth and stable oxygen and carbon isotope ratios, which present an archive for temperature, freshwater fluxes and primary productivity. The purpose of this thesis is the reconstruction of Holocene climate and environmental variations in the North Pacific with a high spatial and temporal resolution using marine bivalve shells. This thesis focuses on several different Holocene time periods and multiple regions in the North Pacific, including: Japan, Alaska (AK), British Columbia (BC) and Washington State, which are affected by the monsoon, Pacific Decadal Oscillation (PDO) and El Niño/Southern Oscillation (ENSO). Such high-resolution proxy data from the marine realm of mid- and high-latitudes are still rare. Therefore, this study contributes to the optimization and verification of climate models. However, before using bivalves for environmental reconstructions and seasonality studies, life history traits must be well studied to temporally align and interpret the geochemical record. These calibration studies are essential to ascertain the usefulness of selected bivalve species as paleoclimate proxy archives. This work focuses on two bivalve species, the short-lived Saxidomus gigantea and the long-lived Panopea abrupta. Sclerochronology and oxygen isotope ratios of different shell layers of P. abrupta were studied in order to test the reliability of this species as a climate archive. The annual increments are clearly discernable in umbonal shell portions and the increments widths should be measured in these shell portions. A reliable reconstruction of paleotemperatures may only be achieved by exclusively sampling the outer shell layer of multiple contemporaneous specimens. Life history traits (e.g., timing of growth line formation, duration of the growing season and growth rates) and stable isotope ratios of recent S. gigantea from AK and BC were analyzed in detail. Furthermore, a growth-temperature model based on S. gigantea shells from Alaska was established, which provides a better understanding of the hydrological changes related to the Alaska Coastal Current (ACC). This approach allows the independent measurement of water temperature and salinity from variations in the width of lunar daily growth increments of S. gigantea. Temperature explains 70% of the variability in shell growth. The model was calibrated and tested with modern shells and then applied to archaeological specimens. The time period between 988 and 1447 cal yrs BP was characterized by colder (~1-2°C) and much drier (2-5 PSU) summers, and a likely much slower flowing ACC than at present. In contrast, the summers during the time interval of 599-1014 cal yrs BP were colder (up to 3°C) and fresher (1-2 PSU) than today. The Aleutian Low may have been stronger and the ACC was probably flowing faster during this time.

The reconstruction of paleoclimate variability in Northern Germany, during marine isotope stage 5

In order to obtain a better understanding about the influence of post-depositional diagenesis on speleothem 230Th/U-ages and paleoclimate variability during Marine Isotope Stage (MIS) 5 in northern Germany, four stalagmites from the Riesenberghöhle (RBH) were investigated by thin section analysis, 230Th/U-dating as well as stable oxygen and carbon isotope and laser ablation inductively coupled mass spectrometry (LA-ICPMS) trace element analysis. The RBH is located in the Weser Hills and is one of the northernmost limestone caves in Germany.rnMulti collector (MC) ICPMS 230Th/U-ages and thin section analysis of the RBH stalagmites shows that some growth phases of the stalagmites were diagenetically altered after their deposition. The impact of post-depositional diagenesis (PDD) on the 230Th/U-ages is modeled, and potential processes leading to PDD are discussed. In this context, it is suggested that PDD may be induced by rapid climate change at the inception of the GIS.rnDespite of the dating uncertainties resulting from PDD, 230Th/U-dating shows that the RBH stalagmites grew during the Eemian and most of the Greenland Interstadials (GIS) during MIS 5. Thus, the growth phases of the RBH stalagmites might be related to a reorganization of the Atlantic Meridional Overturning Circulation (AMOC). The stable isotope (δ13C and δ18O) and the trace element variability of the stalagmites reflects rapid changes of past temperature and precipitation on millennial and sub-millennial timescales. These past climate changes can be amplified by orbitally forced variations of the July solar insolation at 65°N.

Growth hormone receptor polymorphism and growth hormone therapy response in children: a bayesian meta-analysis

Recombinant human growth hormone (rhGH) therapy is used in the long-term treatment of children with growth disorders, but there is considerable treatment response variability. The exon 3-deleted growth hormone receptor polymorphism (GHR(d3)) may account for some of this variability. The authors performed a systematic review (to April 2011), including investigator-only data, to quantify the effects of the GHR(fl-d3) and GHR(d3-d3) genotypes on rhGH therapy response and used a recently established Bayesian inheritance model-free approach to meta-analyze the data. The primary outcome was the 1-year change-in-height standard-deviation score for the 2 genotypes. Eighteen data sets from 12 studies (1,527 children) were included. After several prior assumptions were tested, the most appropriate inheritance model was codominant (posterior probability = 0.93). Compared with noncarriers, carriers had median differences in 1-year change-in-height standard-deviation score of 0.09 (95% credible interval (CrI): 0.01, 0.17) for GHR(fl-d3) and of 0.14 (95% CrI: 0.02, 0.26) for GHR(d3-d3). However, the between-study standard deviation of 0.18 (95% CrI: 0.10, 0.33) was considerable. The authors tested by meta-regression for potential modifiers and found no substantial influence. They conclude that 1) the GHR(d3) polymorphism inheritance is codominant, contrasting with previous reports; 2) GHR(d3) genotypes account for modest increases in rhGH effects in children; and 3) considerable unexplained variability in responsiveness remains.

Isolated autosomal dominant growth hormone deficiency: stimulating mutant GH-1 gene expression drives GH-1 splice-site selection, cell proliferation, and apoptosis

The majority of mutations that cause isolated GH deficiency type II (IGHD II) affect splicing of GH-1 transcripts and produce a dominant-negative GH isoform lacking exon 3 resulting in a 17.5-kDa isoform, which further leads to disruption of the GH secretory pathway. A clinical variability in the severity of the IGHD II phenotype depending on the GH-1 gene alteration has been reported, and in vitro and transgenic animal data suggest that the onset and severity of the phenotype relates to the proportion of 17.5-kDa produced. The removal of GH in IGHD creates a positive feedback loop driving more GH expression, which may itself increase 17.5-kDa isoform productions from alternate splice sites in the mutated GH-1 allele. In this study, we aimed to test this idea by comparing the impact of stimulated expression by glucocorticoids on the production of different GH isoforms from wild-type (wt) and mutant GH-1 genes, relying on the glucocorticoid regulatory element within intron 1 in the GH-1 gene. AtT-20 cells were transfected with wt-GH or mutated GH-1 variants (5'IVS-3 + 2-bp T->C; 5'IVS-3 + 6 bp T->C; ISEm1: IVS-3 + 28 G->A) known to cause clinical IGHD II of varying severity. Cells were stimulated with 1 and 10 mum dexamethasone (DEX) for 24 h, after which the relative amounts of GH-1 splice variants were determined by semiquantitative and quantitative (TaqMan) RT-PCR. In the absence of DEX, only around 1% wt-GH-1 transcripts were the 17.5-kDa isoform, whereas the three mutant GH-1 variants produced 29, 39, and 78% of the 17.5-kDa isoform. DEX stimulated total GH-1 gene transcription from all constructs. Notably, however, DEX increased the amount of 17.5-kDa GH isoform relative to the 22- and 20-kDa isoforms produced from the mutated GH-1 variants, but not from wt-GH-1. This DEX-induced enhancement of 17.5-kDa GH isoform production, up to 100% in the most severe case, was completely blocked by the addition of RU486. In other studies, we measured cell proliferation rates, annexin V staining, and DNA fragmentation in cells transfected with the same GH-1 constructs. The results showed that that the 5'IVS-3 + 2-bp GH-1 gene mutation had a more severe impact on those measures than the splice site mutations within 5'IVS-3 + 6 bp or ISE +28, in line with the clinical severity observed with these mutations. Our findings that the proportion of 17.5-kDa produced from mutant GH-1 alleles increases with increased drive for gene expression may help to explain the variable onset progression, and severity observed in IGHD II.

Whole-body MRI of multiple myeloma: comparison of different MRI sequences in assessment of different growth patterns

PURPOSE: To determine sensitivity, specificity and inter-observer variability of different whole-body MRI (WB-MRI) sequences in patients with multiple myeloma (MM). METHODS AND MATERIALS: WB-MRI using a 1.5T MRI scanner was performed in 23 consecutive patients (13 males, 10 females; mean age 63+/-12 years) with histologically proven MM. All patients were clinically classified according to infiltration (low-grade, n=7; intermediate-grade, n=7; high-grade, n=9) and to the staging system of Durie and Salmon PLUS (stage I, n=12; stage II, n=4; stage III, n=7). The control group consisted of 36 individuals without malignancy (25 males, 11 females; mean age 57+/-13 years). Two observers independently evaluated the following WB-MRI sequences: T1w-TSE (T1), T2w-TIRM (T2), and the combination of both sequences, including a contrast-enhanced T1w-TSE with fat-saturation (T1+/-CE/T2). They had to determine growth patterns (focal and/or diffuse) and the MRI sequence that provided the highest confidence level in depicting the MM lesions. Results were calculated on a per-patient basis. RESULTS: Visual detection of MM was as follows: T1, 65% (sensitivity)/85% (specificity); T2, 76%/81%; T1+/-CE/T2, 67%/88%. Inter-observer variability was as follows: T1, 0.3; T2, 0.55; T1+/-CE/T2, 0.55. Sensitivity improved depending on infiltration grade (T1: 1=60%; 2=36%; 3=83%; T2: 1=70%; 2=71%; 3=89%; T1+/-CE/T2: 1=50%; 2=50%; 3=89%) and clinical stage (T1: 1=58%; 2=63%; 3=79%; T2: 1=58%; 2=88%; 3=100%; T1+/-CE/T2: 1=50%; 2=63%; 3=100%). T2w-TIRM sequences achieved the best reliability in depicting the MM lesions (65% in the mean of both readers). CONCLUSIONS: T2w-TIRM sequences achieved the highest level of sensitivity and best reliability, and thus might be valuable for initial assessment of MM. For an exact staging and grading the examination protocol should encompass unenhanced and enhanced T1w-MRI sequences, in addition to T2w-TIRM.

Variability of SCCmec in the Zurich area

A periodic survey of methicillin-resistant Staphylococcus aureus (MRSA) in Zurich in 2004 and 2006 revealed a consistently low prevalence of MRSA. SCCmec and ccr typing showed fluctuations in the proportions of SCCmec types and in the carriage of mobile virulence determinants. Together with the presence of variant SCCmecs these findings suggest a high clonal diversity and level of SCCmec recombination. The prevalence of a local "drug clone", associated with low-level methicillin resistance and rapid growth, significantly decreased. This clone had spread among intraveneous drug users, steadily increasing from 1994 to 2001 and was dominant in 2001. Apparently, changes in the management of the Zurich drug scene have restricted the spread of this clone.

Climate-induced interannual variability of marine primary and export production in three global coupled climate carbon cycle models

Fully coupled climate carbon cycle models are sophisticated tools that are used to predict future climate change and its impact on the land and ocean carbon cycles. These models should be able to adequately represent natural variability, requiring model validation by observations. The present study focuses on the ocean carbon cycle component, in particular the spatial and temporal variability in net primary productivity (PP) and export production (EP) of particulate organic carbon (POC). Results from three coupled climate carbon cycle models (IPSL, MPIM, NCAR) are compared with observation-based estimates derived from satellite measurements of ocean colour and results from inverse modelling (data assimilation). Satellite observations of ocean colour have shown that temporal variability of PP on the global scale is largely dominated by the permanently stratified, low-latitude ocean (Behrenfeld et al., 2006) with stronger stratification (higher sea surface temperature; SST) being associated with negative PP anomalies. Results from all three coupled models confirm the role of the low-latitude, permanently stratified ocean for anomalies in globally integrated PP, but only one model (IPSL) also reproduces the inverse relationship between stratification (SST) and PP. An adequate representation of iron and macronutrient co-limitation of phytoplankton growth in the tropical ocean has shown to be the crucial mechanism determining the capability of the models to reproduce observed interactions between climate and PP.

Phenotypic variability in familial combined pituitary hormone deficiency caused by a PROP1 gene mutation resulting in the substitution of Arg-->Cys at codon 120 (R120C).

As pituitary function depends on the integrity of the hypothalamic-pituitary axis, any defect in the development and organogenesis of this gland may account for a form of combined pituitary hormone deficiency (CPHD). A mutation in a novel, tissue-specific, paired-like homeodomain transcription factor, termed Prophet of Pit-1 (PROP1), has been identified as causing the Ames dwarf (df) mouse phenotype, and thereafter, different PROP1 gene alterations have been found in humans with CPHD. We report on the follow-up of two consanguineous families (n = 12), with five subjects affected with CPHD (three males and two females) caused by the same nucleotide C to T transition, resulting in the substitution of Arg-->Cys in PROP1 at codon 120. Importantly, there is a variability of phenotype, even among patients with the same mutation. The age at diagnosis was dependent on the severity of symptoms, ranging from 9 months to 8 yr. Although in one patient TSH deficiency was the first symptom of the disorder, all patients became symptomatic by exhibiting severe growth retardation and failure to thrive, which was mainly caused by GH deficiency (n = 4). The secretion of the pituitary-derived hormones (GH, PRL, TSH, LH, and FSH) declined gradually with age, following a different pattern in each individual; therefore, the deficiencies developed over a variable period of time. All of the subjects entered puberty spontaneously, and the two females also experienced menarche and periods before a replacement therapy was necessary.

Heart rate variability in response to pain stimulus in VLBW infants followed longitudinally during NICU stay.

The objective of this longitudinal study, conducted in a neonatal intensive care unit, was to characterize the response to pain of high-risk very low birth weight infants (<1,500 g) from 23 to 38 weeks post-menstrual age (PMA) by measuring heart rate variability (HRV). Heart period data were recorded before, during, and after a heel lanced or wrist venipunctured blood draw for routine clinical evaluation. Pain response to the blood draw procedure and age-related changes of HRV in low-frequency and high-frequency bands were modeled with linear mixed-effects models. HRV in both bands decreased during pain, followed by a recovery to near-baseline levels. Venipuncture and mechanical ventilation were factors that attenuated the HRV response to pain. HRV at the baseline increased with post-menstrual age but the growth rate of high-frequency power was reduced in mechanically ventilated infants. There was some evidence that low-frequency HRV response to pain improved with advancing PMA.

Inherited and noninherited sources of variation in metastasis and the growth rate of tumors: Implications for tumor evolution and therapy

I have undertaken measurements of the genetic (or inherited) and nongenetic (or noninherited) components of the variability of metastasis formation and tumor diameter doubling time in more than 100 metastatic lines from each of three murine tumors (sarcoma SANH, sarcoma SA4020, and hepatocarcinoma HCA-I) syngeneic to C3Hf/Kam mice. These lines were isolated twice from lung metastases and analysed immediately thereafter to obtain the variance to spontaneous lung metastasis and tumor diameter doubling time. Additional studies utilized cells obtained from within 4 passages of isolation. Under the assumption that no genetic differences in metastasis formation or diameter doubling time existed among the cells of a given line, the variance within a line would estimate nongenetic variation. The variability derived from differences between lines would represent genetic origin. The estimates of the genetic contribution to the variation of metastasis and tumor diameter doubling time were significantly greater than zero, but only in the metastatic lines of tumor SANH was genetic variation the major source of metastatic variability (contributing 53% of the variability). In the tumor cell lines of SA4020 and HCA-I, however, the contribution of nongenetic factors predominated over genetic factors in the variability of the number of metastasis and tumor diameter doubling time. A number of other parameters examined, such as DNA content, karyotype, and selection and variance analysis with passage in vivo, indicated that genetic differences existed within the cell lines and that these differences were probably created by genetic instability. The mean metastatic propensity of the lines may have increased somewhat during their isolation and isotransplantation, but the variance was only slightly affected, if at all. Analysis of the DNA profiles of the metastatic lines of SA4020 and HCA-I revealed differences between these lines and their primary parent tumors, but not among the SANH lines and their parent tumor. Furthermore, there was a direct correlation between the extent of genetic influence on metastasis formation and the ability of the tumor cells to develop resistance to cisplatinum. Thus although nongenetic factors might predominate in contributing to metastasis formation, it is probably genetic variation and genetic instability that cause the progression of tumor cells to a more metastatic phenotype and leads to the emergence of drug resistance. ^