Development of a rapid bioassessment method and index of biotic integrity for estuaries located along the northwest Gulf of Mexico

Recently it has been proposed that the evaluation of effects of pollutants on aquatic organisms can provide an early warning system of potential environmental and human health risks (NRC 1991). Unfortunately there are few methods available to aquatic biologists to conduct assessments of the effects of pollutants on aquatic animal community health. The primary goal of this research was to develop and evaluate the feasibility of such a method. Specifically, the primary objective of this study was to develop a prototype rapid bioassessment technique similar to the Index of Biotic Integrity (IBI) for the upper Texas and Northwestern Gulf of Mexico coastal tributaries. The IBI consists of a series of "metrics" which describes specific attributes of the aquatic community. Each of these metrics are given a score which is then subtotaled to derive a total assessment of the "health" of the aquatic community. This IBI procedure may provide an additional assessment tool for professionals in water quality management.^ The experimental design consisted primarily of compiling previously collected data from monitoring conducted by the Texas Natural Resource Conservation Commission (TNRCC) at five bayous classified according to potential for anthropogenic impact and salinity regime. Standardized hydrological, chemical, and biological monitoring had been conducted in each of these watersheds. The identification and evaluation of candidate metrics for inclusion in the estuarine IBI was conducted through the use of correlation analysis, cluster analysis, stepwise and normal discriminant analysis, and evaluation of cumulative distribution frequencies. Scores of each included metric were determined based on exceedances of specific percentiles. Individual scores were summed and a total IBI score and rank for the community computed.^ Results of these analyses yielded the proposed metrics and rankings listed in this report. Based on the results of this study, incorporation of an estuarine IBI method as a water quality assessment tool is warranted. Adopted metrics were correlated to seasonal trends and less so to salinity gradients observed during the study (0-25 ppt). Further refinement of this method is needed using a larger more inclusive data set which includes additional habitat types, salinity ranges, and temporal variation. ^

NOVEL PHANTOMS AND POST-PROCESSING FOR DIFFUSION SPECTRUM IMAGING

High Angular Resolution Diffusion Imaging (HARDI) techniques, including Diffusion Spectrum Imaging (DSI), have been proposed to resolve crossing and other complex fiber architecture in the human brain white matter. In these methods, directional information of diffusion is inferred from the peaks in the orientation distribution function (ODF). Extensive studies using histology on macaque brain, cat cerebellum, rat hippocampus and optic tracts, and bovine tongue are qualitatively in agreement with the DSI-derived ODFs and tractography. However, there are only two studies in the literature which validated the DSI results using physical phantoms and both these studies were not performed on a clinical MRI scanner. Also, the limited studies which optimized DSI in a clinical setting, did not involve a comparison against physical phantoms. Finally, there is lack of consensus on the necessary pre- and post-processing steps in DSI; and ground truth diffusion fiber phantoms are not yet standardized. Therefore, the aims of this dissertation were to design and construct novel diffusion phantoms, employ post-processing techniques in order to systematically validate and optimize (DSI)-derived fiber ODFs in the crossing regions on a clinical 3T MR scanner, and develop user-friendly software for DSI data reconstruction and analysis. Phantoms with a fixed crossing fiber configuration of two crossing fibers at 90° and 45° respectively along with a phantom with three crossing fibers at 60°, using novel hollow plastic capillaries and novel placeholders, were constructed. T2-weighted MRI results on these phantoms demonstrated high SNR, homogeneous signal, and absence of air bubbles. Also, a technique to deconvolve the response function of an individual peak from the overall ODF was implemented, in addition to other DSI post-processing steps. This technique greatly improved the angular resolution of the otherwise unresolvable peaks in a crossing fiber ODF. The effects of DSI acquisition parameters and SNR on the resultant angular accuracy of DSI on the clinical scanner were studied and quantified using the developed phantoms. With a high angular direction sampling and reasonable levels of SNR, quantification of a crossing region in the 90°, 45° and 60° phantoms resulted in a successful detection of angular information with mean ± SD of 86.93°±2.65°, 44.61°±1.6° and 60.03°±2.21° respectively, while simultaneously enhancing the ODFs in regions containing single fibers. For the applicability of these validated methodologies in DSI, improvement in ODFs and fiber tracking from known crossing fiber regions in normal human subjects were demonstrated; and an in-house software package in MATLAB which streamlines the data reconstruction and post-processing for DSI, with easy to use graphical user interface was developed. In conclusion, the phantoms developed in this dissertation offer a means of providing ground truth for validation of reconstruction and tractography algorithms of various diffusion models (including DSI). Also, the deconvolution methodology (when applied as an additional DSI post-processing step) significantly improved the angular accuracy of the ODFs obtained from DSI, and should be applicable to ODFs obtained from the other high angular resolution diffusion imaging techniques.

Comparison of Tumor Shrinkage and Cumulative Dose Distribution for Lung Cancers

Background: The physical characteristic of protons is that they deliver most of their radiation dose to the target volume and deliver no dose to the normal tissue distal to the tumor. Previously, numerous studies have shown unique advantages of proton therapy over intensity-modulated radiation therapy (IMRT) in conforming dose to the tumor and sparing dose to the surrounding normal tissues and the critical structures in many clinical sites. However, proton therapy is known to be more sensitive to treatment uncertainties such as inter- and intra-fractional variations in patient anatomy. To date, no study has clearly demonstrated the effectiveness of proton therapy compared with the conventional IMRT under the consideration of both respiratory motion and tumor shrinkage in non-small cell lung cancer (NSCLC) patients. Purpose: This thesis investigated two questions for establishing a clinically relevant comparison of the two different modalities (IMRT and proton therapy). The first question was whether or not there are any differences in tumor shrinkage between patients randomized to IMRT versus passively scattered proton therapy (PSPT). Tumor shrinkage is considered a standard measure of radiation therapy response that has been widely used to gauge a short-term progression of radiation therapy. The second question was whether or not there are any differences between the planned dose and 5D dose under the influence of inter- and intra-fractional variations in the patient anatomy for both modalities. Methods: A total of 45 patients (25 IMRT patients and 20 PSPT patients) were used to quantify the tumor shrinkage in terms of the change of the primary gross tumor volume (GTVp). All patients were randomized to receive either IMRT or PSPT for NSCLC. Treatment planning goals were identical for both groups. All patients received 5 to 8 weekly repeated 4-dimensional computed tomography (4DCT) scans during the course of radiation treatments. The original GTVp contours were propagated to T50 of weekly 4DCT images using deformable image registration and their absolute volumes were measured. Statistical analysis was performed to compare the distribution of tumor shrinkage between the two population groups. In order to investigate the difference between the planned dose and the 5D dose with consideration of both breathing motion and anatomical change, we re-calculated new dose distributions at every phase of the breathing cycle for all available weekly 4DCT data sets which resulted 50 to 80 individual dose calculations for each of the 7 patients presented in this thesis. The newly calculated dose distributions were then deformed and accumulated to T50 of the planning 4DCT for comparison with the planned dose distribution. Results: At the end of the treatment, both IMRT and PSPT groups showed mean tumor volume reductions of 23.6% ( 19.2%) and 20.9% ( 17.0 %) respectively. Moreover, the mean difference in tumor shrinkage between two groups is 3% along with the corresponding 95% confidence interval, [-8%, 14%]. The rate of tumor shrinkage was highly correlated with the initial tumor volume size. For the planning dose and 5D dose comparison study, all 7 patients showed a mean difference of 1 % in terms of target coverage for both IMRT and PSPT treatment plans. Conclusions: The results of the tumor shrinkage investigation showed no statistically significant difference in tumor shrinkage between the IMRT and PSPT patients, and the tumor shrinkage between the two modalities is similar based on the 95% confidence interval. From the pilot study of comparing the planned dose with the 5D dose, we found the difference to be only 1%. Overall impression of the two modalities in terms of treatment response as measured by the tumor shrinkage and 5D dose under the influence of anatomical change that were designed under the same protocol (i.e. randomized trial) showed similar result.

Expression and function of $\beta$-tubulin isotypes in Chinese hamster ovary (CHO) cells

The availability of isotype specific antisera for $\beta$-tubulin, coupled with genetic and biochemical analysis, has allowed the determination of $\beta$-tubulin isotype expression and distribution in Chinese hamster ovary (CHO) cells. Using genetic manipulations involving selection for colcemid resistance followed by reversion and reselection for drug resistance, we have succeeded in isolating cell lines that exhibit three major and one minor $\beta$-tubulin spots by two-dimensional gel electrophoresis. In concert with isotype specific antibodies, analysis of these mutants demonstrates that CHO cells express two copies of isotype I, at least one copy of isotype IV, and very small amounts of isotype V. Their stoichiometry is approximately 1:1:0.7:0.2. All three isotypes assemble into both cytoplasmic and spindle microtubules, and are similar in their responses to cold, colcemid, and calcium induced depolymerization. They have comparable turnover rates and are equally sensitive to depression of synthesis upon colchicine treatment. These results suggest that $\beta$-tubulin isotypes are used interchangeably to assemble microtubule structures in CHO cells. However, of 18 colcemid resistant mutants with a demonstrable alteration in $\beta$-tubulin, all were found to have the alteration in isotype I, thus leaving open the possibility that subtle differences in isotype properties may exist. Under various conditions of the cell growth, the relative proportion of each expressed isotype does not significantly seem to change except in the early G1 phase of the cell cycle. At this time the synthesis of isotype V increases more than two fold relative to isotype I and IV, while at the same time, total $\beta$-tubulin synthesis is decreased about 60-70%. ^

Cell surface galactosyltransferase structure and function during mesenchymal cell migration

In order to more fully understand the function of surface GalTase on mesenchymal cells, anti-GalTase IgG was used to (a) examine the role of surface GalTase during mouse mesenchymal cell migration on laminin and fibronectin; (b) define the plasma membrane distribution of GalTase by indirect immunofluorescence on migrating cells; (c) quantitate the level of surface GalTase on migrating cells; and (d) determine whether GalTase is associated with the cytoskeleton.^ Results show that anti-GalTase IgG was able to inhibit migration (48-80% as compared to basal rate) when cells were migrating on laminin-containing matrices. Monovalent Fab fragments inhibited migration on laminin by 90% after 4 hours. On the other hand, anti-GalTase IgG had no effect on cells migrating on fibronectin. This illustrates the substrate specificity of GalTase mediated-migration. When anti-GalTase IgG was used to localize surface GalTase on cells migratory on laminin, the enzyme was restricted to the leading and trailing edges of the cell. Assays indicate that GalTase is elevated approximately 3-fold when cells are migrating on laminin-containing matrices as compared to migratory cells on plastic or fibronectin, or as compared to stationary cells on any substrate. Laminin appears to recruit GalTase from preexisting intracellular pools to the growing lamellipodia.^ Double-label indirect immunofluorescence studies indicate that there is an apparent co-localization between some of the surface GalTase and some actin filaments. This relationship was explored by extracting cells prelabeled with anti-GalTase IgG and quantitated by radiolabeled second antibodies. Results show that 79% of the surface GalTase is associated with the cytoskeleton (as judged by detergent insolubility) when monovalent antibodies (Fab) are used. However virtually all (80-100%) of the surface GalTase can be induced to associate with the cytoskeleton when cross-linked with bivalent antibodies. Furthermore, when cells in suspension are incubated with divalent antibodies, an additional 66% of the surface GalTase can be induced to associate with the cytoskeleton. The elevated levels of surface GalTase detectable on cells migrating on laminin also appear to be associated with the cytoskeleton.^ Several lines of evidence suggest that GalTase is associated with F-actin. Data suggest that laminin induces the expression of surface GalTase to the growing lamellipodia where it becomes associated with the cytoskeleton leading to cell spreading and migration. (Abstract shortened with permission of author.) ^

Cytochrome P450 in mammalian brain: Purification, anatomical distribution and regulation of the cytochrome P450 monooxygenase system in rat brain

The cytochrome P450 (P450) monooxygenase system plays a major role in metabolizing a wide variety of xenobiotic as well as endogenous compounds. In performing this function, it serves to protect the body from foreign substances. However, in a number of cases, P450 activates procarcinogens to cause harm. In most animals, the highest level of activity is found in the liver. Virtually all tissues demonstrate P450 activity, though, and the role of the P450 monooxygenase system in these other organs is not well understood. In this project I have studied the P450 system in rat brain; purifying NADPH-cytochrome P450 reductase (reductase) from that tissue. In addition, I have examined the distribution and regulation of expression of reductase and P450 in various anatomical regions of the rat brain.^ NADPH-cytochrome P450 reductase was purified to apparent homogeneity and cytochrome P450 partially purified from whole rat brain. Purified reductase from brain was identical to liver P450 reductase by SDS-PAGE and Western blot techniques. Kinetic studies utilizing cerebral P450 reductase reveal Km values in close agreement with those determined with enzyme purified from rat liver. Moreover, the brain P450 reductase was able to function successfully in a reconstituted microsomal system with partially purified brain cytochrome P450 and with purified hepatic P4501A1 as measured by 7-ethoxycoumarin and 7-ethoxyresorufin O-deethylation. These results indicate that the reductase and P450 components may interact to form a competent drug metabolism system in brain tissue.^ Since the brain is not a homogeneous organ, dependent upon the well orchestrated interaction of numerous parts, pathology in one nucleus may have a large impact upon its overall function. Hence, the anatomical distribution of the P450 monooxygenase system in brain is important in elucidating its function in that organ. Related to this is the regulation of P450 expression in brain. In order to study these issues female rats--both ovariectomized and not--were treated with a number of xenobiotic compounds and sex steroids. The brains from these animals were dissected into 8 discrete regions and the presence and relative level of message for P4502D and reductase determined using polymerase chain reaction. Results of this study indicate the presence of mRNA for reductase and P4502D isoforms throughout the rat brain. In addition, quantitative PCR has allowed the determination of factors affecting the expression of message for these enzymes. ^

Structure and function study of prostaglandin H synthase

Prostaglandin H synthase (PGHS) is a key enzyme in biosynthesis of prostaglandins, thromboxane, and prostacyclin. It has two activities, cyclooxygenase and peroxidase. "PGHS" means PGHS-1. A current hypothesis considers the cyclooxygenase reaction to be a free radical chain reaction, initiated by interaction of the synthase peroxidase with hydroperoxides leading to the production of a tyrosyl free radical. According to this hypothesis, tyrosyl residue(s) may play a key role in the cyclooxygenase reaction. Tetranitromethane (TNM) can relatively selectively nitrate tyrosines at pH 8.0. The effect of TNM on both cyclooxygenase activity and peroxidase activity has been examined: reaction of the synthase holoenzyme with TNM at pH 8.0 led to inactivation of both activities, with the cyclooxygenase activity being lost rapidly and completely, while the peroxidase activity was lost more slowly. Indomethacin, a non-steroidal anti-inflammatory agent, can protect the synthase from the inactivation of TNM. Amino acid analyses indicated that a loss of tyrosine and formation of nitrotyrosine residues occurred during reaction with TNM, and that TNM-reacted holoenzyme with $<$10% residual cyclooxygenase activity had about 2.0 nitrotyrosine/subunit.^ PGH synthase is known to be an endoplasmic reticulum membrane-associated protein. Antibodies directed at particular PGHS peptide segments and indirect immunofluorescence have been used to characterize the membrane topology of crucial portions of PGHS. PGHS was expressed in COS-1 cells transfected with the appropriate cDNA. Stably-transfected human endothelial cells were also used for the topology study. The cells were treated with streptolysin-O, which selectively permeabilizes the plasma membrane, or with saponin to achieve general membrane disruption, before incubation with the antipeptide antibodies. Bound antipeptide antibody was stained by FITC-labelled secondary antibody and visualized by fluorescence microscopy. With the antipeptide antibodies against residues 51-66, 156-170 or 377-390, there was a significant reticular and perinuclear pattern of staining in cells permeabilized with saponin but not in cells permeabilized with SLO alone. Antibodies directed against the endogenous C-terminal peptide or against residues 271-284 produced staining in cells permeabilized with saponin, and also in a lower, but significant fraction of cells permeabilized with SLO. Similar results were obtained when COS-1 cells expressing recombinant PGHS with a viral reporter peptide inserted at the C-terminus were stained with antibody against the reporter epitope.^ The PGHS C-terminal sequence is similar to that of the consensus KDEL ER retention signal. The potential function of the PGHS C-terminus segment in ER retention was examined by mutating this segment and analyzing the subcellular distribution of the mutants expressed in COS-1 cells. None of the mutants had an altered subcellular distribution, although some had greatly diminished the enzyme activities. (Abstract shortened by UMI.) ^

Survival analysis of longevity in siblings

Do siblings of centenarians tend to have longer life spans? To answer this question, life spans of 184 siblings for 42 centenarians have been evaluated. Two important questions have been addressed in analyzing the sibling data. First, a standard needs to be established, to which the life spans of 184 siblings are compared. In this report, an external reference population is constructed from the U.S. life tables. Its estimated mortality rates are treated as baseline hazards from which the relative mortality of the siblings are estimated. Second, the standard survival models which assume independent observations are invalid when correlation within family exists, underestimating the true variance. Methods that allow correlations are illustrated by three different methods. First, the cumulative relative excess mortality between siblings and their comparison group is calculated and used as an effective graphic tool, along with the Product Limit estimator of the survival function. The variance estimator of the cumulative relative excess mortality is adjusted for the potential within family correlation using Taylor linearization approach. Second, approaches that adjust for the inflated variance are examined. They are adjusted one-sample log-rank test using design effect originally proposed by Rao and Scott in the correlated binomial or Poisson distribution setting and the robust variance estimator derived from the log-likelihood function of a multiplicative model. Nether of these two approaches provide correlation estimate within families, but the comparison with the comparison with the standard remains valid under dependence. Last, using the frailty model concept, the multiplicative model, where the baseline hazards are known, is extended by adding a random frailty term that is based on the positive stable or the gamma distribution. Comparisons between the two frailty distributions are performed by simulation. Based on the results from various approaches, it is concluded that the siblings of centenarians had significant lower mortality rates as compared to their cohorts. The frailty models also indicate significant correlations between the life spans of the siblings. ^

Effects of informative censoring on the hazard function: Application to the proportional hazards regression model

The standard analyses of survival data involve the assumption that survival and censoring are independent. When censoring and survival are related, the phenomenon is known as informative censoring. This paper examines the effects of an informative censoring assumption on the hazard function and the estimated hazard ratio provided by the Cox model.^ The limiting factor in all analyses of informative censoring is the problem of non-identifiability. Non-identifiability implies that it is impossible to distinguish a situation in which censoring and death are independent from one in which there is dependence. However, it is possible that informative censoring occurs. Examination of the literature indicates how others have approached the problem and covers the relevant theoretical background.^ Three models are examined in detail. The first model uses conditionally independent marginal hazards to obtain the unconditional survival function and hazards. The second model is based on the Gumbel Type A method for combining independent marginal distributions into bivariate distributions using a dependency parameter. Finally, a formulation based on a compartmental model is presented and its results described. For the latter two approaches, the resulting hazard is used in the Cox model in a simulation study.^ The unconditional survival distribution formed from the first model involves dependency, but the crude hazard resulting from this unconditional distribution is identical to the marginal hazard, and inferences based on the hazard are valid. The hazard ratios formed from two distributions following the Gumbel Type A model are biased by a factor dependent on the amount of censoring in the two populations and the strength of the dependency of death and censoring in the two populations. The Cox model estimates this biased hazard ratio. In general, the hazard resulting from the compartmental model is not constant, even if the individual marginal hazards are constant, unless censoring is non-informative. The hazard ratio tends to a specific limit.^ Methods of evaluating situations in which informative censoring is present are described, and the relative utility of the three models examined is discussed. ^

A Bayesian approach to estimating the regression coefficients of a multinomial logit model

A Bayesian approach to estimation of the regression coefficients of a multinominal logit model with ordinal scale response categories is presented. A Monte Carlo method is used to construct the posterior distribution of the link function. The link function is treated as an arbitrary scalar function. Then the Gauss-Markov theorem is used to determine a function of the link which produces a random vector of coefficients. The posterior distribution of the random vector of coefficients is used to estimate the regression coefficients. The method described is referred to as a Bayesian generalized least square (BGLS) analysis. Two cases involving multinominal logit models are described. Case I involves a cumulative logit model and Case II involves a proportional-odds model. All inferences about the coefficients for both cases are described in terms of the posterior distribution of the regression coefficients. The results from the BGLS method are compared to maximum likelihood estimates of the regression coefficients. The BGLS method avoids the nonlinear problems encountered when estimating the regression coefficients of a generalized linear model. The method is not complex or computationally intensive. The BGLS method offers several advantages over Bayesian approaches. ^

Expression and function of $\beta$-1,4-galactosyltransferase during wild-type and T/{\it t\/} -mutant spermatogenesis

In the mouse, gamete recognition is mediated in part by the binding of sperm surface $\beta$1,4 galactosyltransferase (GalTase) to specific oligosaccharide residues on the zona pellucida ZP3. The expression of GalTase on the sperm surface is regulated by alleles within the distal segment of the T/t complex and results in a haploid-specific increase in GalTase expression on spermatids and sperm from t-bearing males, suggesting that differences in sperm GalTase activity may contribute to t-sperm transmission ratio distortion. In this study, the expression of GalTase RNA during wild-type and T/t-mutant spermatogenesis was characterized and the role of GalTase was analyzed in transmission ratio distortion. It was found that spermatogenic cells predominantly express the long form of the GalTase RNA, which encodes the GalTase protein that is preferentially targeted to the cell surface in somatic cells. In wild-type testes, GalTase RNA accumulates during the maturation of primary spermatocytes, reaches peak levels prior to meiosis, and decreases and meiosis. GalTase RNA accumulates to similar levels during the maturation of +/t and t/t primary spermatocytes, but unlike wild-type, the level of GalTase RNA in t-spermatocytes remains elevated during meiotic division. Consequently, spermatids in t-mutant testes inherit higher levels of GalTase RNA than do wild-type spermatids, which likely accounts for the haploid-specific increase in surface GalTase activity characteristic of spermatids from t-bearing mice.^ The functional significance of the increased GalTase activity during t-sperm transmission ratio distortion was determined by examining the distribution of GalTase RNA and surface GalTase protein in haploid spermatids from +/t males. Results show that +- and t-spermatids have similar levels of both GalTase RNA and protein, indicating that transmission ratio distortion in +/t mice is not likely due to haploid-specific differences in sperm surface GalTase activity.^ The presence of GalTase on the surface of an early spermatogenic cells before it is required on the mature sperm to perform its function during gamete binding suggests a separate function for GalTase in Sertoli-germ cell adhesion. Studies indicate that cell surface GalTase partly mediates the initial adhesion of pachytene spermatocytes, but not haploid spermatids, to Sertoli cells. ^