Implications of fault constitutive properties for earthquake prediction.

The rate- and state-dependent constitutive formulation for fault slip characterizes an exceptional variety of materials over a wide range of sliding conditions. This formulation provides a unified representation of diverse sliding phenomena including slip weakening over a characteristic sliding distance Dc, apparent fracture energy at a rupture front, time-dependent healing after rapid slip, and various other transient and slip rate effects. Laboratory observations and theoretical models both indicate that earthquake nucleation is accompanied by long intervals of accelerating slip. Strains from the nucleation process on buried faults generally could not be detected if laboratory values of Dc apply to faults in nature. However, scaling of Dc is presently an open question and the possibility exists that measurable premonitory creep may precede some earthquakes. Earthquake activity is modeled as a sequence of earthquake nucleation events. In this model, earthquake clustering arises from sensitivity of nucleation times to the stress changes induced by prior earthquakes. The model gives the characteristic Omori aftershock decay law and assigns physical interpretation to aftershock parameters. The seismicity formulation predicts large changes of earthquake probabilities result from stress changes. Two mechanisms for foreshocks are proposed that describe observed frequency of occurrence of foreshock-mainshock pairs by time and magnitude. With the first mechanism, foreshocks represent a manifestation of earthquake clustering in which the stress change at the time of the foreshock increases the probability of earthquakes at all magnitudes including the eventual mainshock. With the second model, accelerating fault slip on the mainshock nucleation zone triggers foreshocks.

Slip complexity in earthquake fault models.

We summarize studies of earthquake fault models that give rise to slip complexities like those in natural earthquakes. For models of smooth faults between elastically deformable continua, it is critical that the friction laws involve a characteristic distance for slip weakening or evolution of surface state. That results in a finite nucleation size, or coherent slip patch size, h*. Models of smooth faults, using numerical cell size properly small compared to h*, show periodic response or complex and apparently chaotic histories of large events but have not been found to show small event complexity like the self-similar (power law) Gutenberg-Richter frequency-size statistics. This conclusion is supported in the present paper by fully inertial elastodynamic modeling of earthquake sequences. In contrast, some models of locally heterogeneous faults with quasi-independent fault segments, represented approximately by simulations with cell size larger than h* so that the model becomes "inherently discrete," do show small event complexity of the Gutenberg-Richter type. Models based on classical friction laws without a weakening length scale or for which the numerical procedure imposes an abrupt strength drop at the onset of slip have h* = 0 and hence always fall into the inherently discrete class. We suggest that the small-event complexity that some such models show will not survive regularization of the constitutive description, by inclusion of an appropriate length scale leading to a finite h*, and a corresponding reduction of numerical grid size.

The organization of seismicity on fault networks.

Although models of homogeneous faults develop seismicity that has a Gutenberg-Richter distribution, this is only a transient state that is followed by events that are strongly influenced by the nature of the boundaries. Models with geometrical inhomogeneities of fracture thresholds can limit the sizes of earthquakes but now favor the characteristic earthquake model for large earthquakes. The character of the seismicity is extremely sensitive to distributions of inhomogeneities, suggesting that statistical rules for large earthquakes in one region may not be applicable to large earthquakes in another region. Model simulations on simple networks of faults with inhomogeneities of threshold develop episodes of lacunarity on all members of the network. There is no validity to the popular assumption that the average rate of slip on individual faults is a constant. Intermediate term precursory activity such as local quiescence and increases in intermediate-magnitude activity at long range are simulated well by the assumption that strong weakening of faults by injection of fluids and weakening of asperities on inhomogeneous models of fault networks is the dominant process; the heat flow paradox, the orientation of the stress field, and the low average stress drop in some earthquakes are understood in terms of the asperity model of inhomogeneous faulting.

Geometric incompatibility in a fault system.

Interdependence between geometry of a fault system, its kinematics, and seismicity is investigated. Quantitative measure is introduced for inconsistency between a fixed configuration of faults and the slip rates on each fault. This measure, named geometric incompatibility (G), depicts summarily the instability near the fault junctions: their divergence or convergence ("unlocking" or "locking up") and accumulation of stress and deformations. Accordingly, the changes in G are connected with dynamics of seismicity. Apart from geometric incompatibility, we consider deviation K from well-known Saint Venant condition of kinematic compatibility. This deviation depicts summarily unaccounted stress and strain accumulation in the region and/or internal inconsistencies in a reconstruction of block- and fault system (its geometry and movements). The estimates of G and K provide a useful tool for bringing together the data on different types of movement in a fault system. An analog of Stokes formula is found that allows determination of the total values of G and K in a region from the data on its boundary. The phenomenon of geometric incompatibility implies that nucleation of strong earthquakes is to large extent controlled by processes near fault junctions. The junctions that have been locked up may act as transient asperities, and unlocked junctions may act as transient weakest links. Tentative estimates of K and G are made for each end of the Big Bend of the San Andreas fault system in Southern California. Recent strong earthquakes Landers (1992, M = 7.3) and Northridge (1994, M = 6.7) both reduced K but had opposite impact on G: Landers unlocked the area, whereas Northridge locked it up again.

Long-term in vitro correction of alpha-L-iduronidase deficiency (Hurler syndrome) in human bone marrow.

Allogeneic bone marrow transplantation is the most effective treatment for Hurler syndrome but, since this therapy is not available to all patients, we have considered an alternative approach based on transfer and expression of the normal gene in autologous bone marrow. A retroviral vector carrying the full-length cDNA for alpha-L-iduronidase has been constructed and used to transduce bone marrow from patients with this disorder. Various gene-transfer protocols have been assessed including the effect of intensive schedules of exposure of bone marrow to viral supernatant and the influence of growth factors. With these protocols, we have demonstrated successful gene transfer into primitive CD34+ cells and subsequent enzyme expression in their maturing progeny. Also, by using long-term bone marrow cultures, we have demonstrated high levels of enzyme expression sustained for several months. The efficiency of gene transfer has been assessed by PCR analysis of hemopoietic colonies as 25-56%. No advantage has been demonstrated for the addition of growth factors or intensive viral exposure schedules. The enzyme is secreted into the medium and functional localization has been demonstrated by reversal of the phenotypic effects of lysosomal storage in macrophages. This work suggests that retroviral gene transfer into human bone marrow may offer the prospect for gene therapy of Hurler syndrome in young patients without a matched sibling donor.

Targeted gene correction of episomal DNA in mammalian cells mediated by a chimeric RNA.DNA oligonucleotide.

An experimental strategy to facilitate correction of single-base mutations of episomal targets in mammalian cells has been developed. The method utilizes a chimeric oligonucleotide composed of a contiguous stretch of RNA and DNA residues in a duplex conformation with double hairpin caps on the ends. The RNA/DNA sequence is designed to align with the sequence of the mutant locus and to contain the desired nucleotide change. Activity of the chimeric molecule in targeted correction was tested in a model system in which the aim was to correct a point mutation in the gene encoding the human liver/bone/kidney alkaline phosphatase. When the chimeric molecule was introduced into cells containing the mutant gene on an extrachromosomal plasmid, correction of the point mutation was accomplished with a frequency approaching 30%. These results extend the usefulness of the oligonucleotide-based gene targeting approaches by increasing specific targeting frequency. This strategy should enable the design of antiviral agents.

Selection of transduced CD34+ progenitors and enzymatic correction of cells from Gaucher patients, with bicistronic vectors.

The gene transfer efficiency of human hematopoietic stem cells is still inadequate for efficient gene therapy of most disorders. To overcome this problem, a selectable retroviral vector system for gene therapy has been developed for gene therapy of Gaucher disease. We constructed a bicistronic retroviral vector containing the human glucocerebrosidase (GC) cDNA and the human small cell surface antigen CD24 (243 bp). Expression of both cDNAs was controlled by the long terminal repeat enhancer/promoter of the Molony murine leukemia virus. The CD24 selectable marker was placed downstream of the GC cDNA and its translation was enhanced by inclusion of the long 5' untranslated region of encephalomyocarditis virus internal ribosomal entry site. Virus-producing GP+envAM12 cells were created by multiple supernatant transductions to create vector producer cells. The vector LGEC has a high titer and can drive expression of GC and the cell surface antigen CD24 simultaneously in transduced NIH 3T3 cells and Gaucher skin fibroblasts. These transduced cells have been successfully separated from untransduced cells by fluorescence-activated cell sorting, based on cell surface expression of CD24. Transduced and sorted NIH 3T3 cells showed higher GC enzyme activity than the unsorted population, demonstrating coordinated expression of both genes. Fibroblasts from Gaucher patients were transduced and sorted for CD24 expression, and GC enzyme activity was measured. The transduced sorted Gaucher fibroblasts had a marked increase in enzyme activity (149%) compared with virgin Gaucher fibroblasts (17% of normal GC enzyme activity). Efficient transduction of CD34+ hematopoietic progenitors (20-40%) was accomplished and fluorescence-activated cell sorted CD24(+)-expressing progenitors generated colonies, all of which (100%) were vector positive. The sorted, CD24-expressing progenitors generated erythroid burst-forming units, colony-forming units (CFU)-granulocyte, CFU-macrophage, CFU-granulocyte/macrophage, and CFU-mix hematopoietic colonies, demonstrating their ability to differentiate into these myeloid lineages in vitro. The transduced, sorted progenitors raised the GC enzyme levels in their progeny cells manyfold compared with untransduced CD34+ progenitors. Collectively, this demonstrates the development of high titer, selectable bicistronic vectors that allow isolation of transduced hematopoietic progenitors and cells that have been metabolically corrected.

Fixing human factor IX (fIX): correction of a cryptic RNA splice enables the production of biologically active fIX in the mammary gland of transgenic mice.

Transgenic mice and sheep secrete only low levels of human factor IX in their milk because of an aberrant splicing of the transgene RNA in the mammary gland. Removal of the cryptic 3' splice site prevents this splicing and leads to the production of relatively high levels of factor IX. The purified protein is fully active showing that the mammary gland is capable of the efficient post-translational modification of this protein and that transgenic animals are a suitable means of its production.

Phenotype correction in retinal pigment epithelium in murine mucopolysaccharidosis VII by adenovirus-mediated gene transfer.

We have studied the use of adenovirus-mediated gene transfer to reverse the pathologic changes of lysosomal storage disease caused by beta-glucuronidase deficiency in the eyes of mice with mucopolysaccharidosis VII. A recombinant adenovirus carrying the human beta-glucuronidase cDNA coding region under the control of a non-tissue-specific promoter was injected intravitreally or subretinally into the eyes of mice with mucopolysaccharidosis VII. At 1-3 weeks after injection, the treated and control eyes were examined histochemically for beta-glucuronidase expression and histologically for phenotypic correction of the lysosomal storage defect. Enzymatic expression was detected 1-3 weeks after injection. Storage vacuoles in the retinal pigment epithelium (RPE) were still present 1 week after gene transfer but were reduced to undetectable levels by 3 weeks in both intravitreally and subretinally injected eyes. There was minimal evidence of ocular pathology associated with the viral injection. These data indicate that adenovirus-mediated gene transfer to the eye may provide for adjunctive therapy for lysosomal storage diseases affecting the RPE in conjunction with enzyme replacement and/or gene therapies for correction of systemic disease manifestations. The data also support the view that recombinant adenovirus may be useful as a gene therapy vector for retinal degenerations that result from a primary genetic defect in the RPE cells.

Correction of chromosomal instability and sensitivity to diverse mutagens by a cloned cDNA of the XRCC3 DNA repair gene.

The mutagen-sensitive CHO line irs1SF was previously isolated on the basis of hypersensitivity to ionizing radiation and was found to be chromosomally unstable as well as cross-sensitive to diverse kinds of DNA-damaging agents. The analysis of somatic cell hybrids formed between irs1SF and human lymphocytes implicated a human gene (defined as XRCC3; x-ray repair cross-complementing), which partially restored mitomycin C resistance to the mutant. A functional cDNA that confers mitomycin C resistance was transferred to irs1SF cells by transforming them with an expression cDNA library and obtaining primary and secondary transformants. Functional cDNA clones were recovered from a cosmid library prepared from a secondary transformant. Transformants also showed partial correction of sensitivity to cisplatin and gamma-rays, efficient correction of chromosomal instability, and substantially improved plating efficiency and growth rate. The XRCC3 cDNA insert is approximately 2.5 kb and detects an approximately 3.0-kb mRNA on Northern blots. The cDNA was mapped by fluorescence in situ hybridization to human chromosome 14q32.3, which was consistent with the chromosome concordance data of two independent hybrid clone panels.