Statistical Hurdle Models for Single Cell Gene Expression: Differential Expression and Graphical Modeling

Thesis (Ph.D.)--University of Washington, 2016-06

Non-linear analysis of the thermo-electro-mechanical behaviour of shear deformable FGM plates with piezoelectric actuators

This paper investigates the non-linear bending behaviour of functionally graded plates that are bonded with piezoelectric actuator layers and subjected to transverse loads and a temperature gradient based on Reddy's higher-order shear deformation plate theory. The von Karman-type geometric non-linearity, piezoelectric and thermal effects are included in mathematical formulations. The temperature change is due to a steady-state heat conduction through the plate thickness. The material properties are assumed to be graded in the thickness direction according to a power-law distribution in terms of the volume fractions of the constituents. The plate is clamped at two opposite edges, while the remaining edges can be free, simply supported or clamped. Differential quadrature approximation in the X-axis is employed to convert the partial differential governing equations and the associated boundary conditions into a set of ordinary differential equations. By choosing the appropriate functions as the displacement and stress functions on each nodal line and then applying the Galerkin procedure, a system of non-linear algebraic equations is obtained, from which the non-linear bending response of the plate is determined through a Picard iteration scheme. Numerical results for zirconia/aluminium rectangular plates are given in dimensionless graphical form. The effects of the applied actuator voltage, the volume fraction exponent, the temperature gradient, as well as the characteristics of the boundary conditions are also studied in detail. Copyright (C) 2004 John Wiley Sons, Ltd.

Whole genome analysis reveals a high incidence of non-optimal codons in secretory signal sequences of Escherichia coli

Translational pausing may occur due to a number of mechanisms, including the presence of non-optimal codons, and it is thought to play a role in the folding of specific polypeptide domains during translation and in the facilitation of signal peptide recognition during see-dependent protein targeting. In this whole genome analysis of Escherichia coli we have found that non-optimal codons in the signal peptide-encoding sequences of secretory genes are overrepresented relative to the mature portions of these genes; this is in addition to their overrepresentation in the 5'-regions of genes encoding non-secretory proteins. We also find increased non-optimal codon usage at the 3' ends of most E. coli genes, in both non-secretory and secretory sequences. Whereas presumptive translational pausing at the 5' and 3' ends of E. coli messenger RNAs may clearly have a general role in translation, we suggest that it also has a specific role in sec-dependent protein export, possibly in facilitating signal peptide recognition. This finding may have important implications for our understanding of how the majority of non-cytoplasmic proteins are targeted, a process that is essential to all biological cells. (C) 2004 Elsevier Inc. All rights reserved.

Expression of the cell surface mucin gene family in adenocarcinomas

Cell surface mucins are complex glycoproteins expressed on the apical membrane surface of mucosal epithelial cells. In malignant epithelial cells they are thought to influence cell adhesion, and are clinical targets for tumor immunotherapy and serum tumor marker assays. We have compared expression of MUC1, MUC3, MUC4, MUC11, MUC12 and MUC13 mRNA in epithelial cancers and/or cell lines with non-malignant tissues. In non-malignant tissues, MUC3, 4, 11, 12 and 13 were expressed at highest levels in gastrointestinal tissues, whereas MUC1 was more widely distributed. Significant down-regulation of the MUC4, MUC12 and MUC13 genes was observed in colonic cancers compared with normal tissue, whereas MUC1 was upregulated. In rectal cancers, levels of all six mucin genes were not significantly different to those in normal rectal tissues. Both MUC1 and MUC4 were down-regulated in gastric cancers, whereas cancer and normal tissue levels were similar for MUC3, 11, 12 and 13. In esophageal cancers there was a general trend toward higher levels than in normal tissue for MUC1, 3, 12 and 13. In ovarian cancers MUC1 levels were very high, whereas only low levels of all other mucins were observed. We also report expression in renal cell carcinomas, bladder carcinomas and breast cancer cell lines. The reported expression profiles of the cell surface mucin gene family will help direct biological and clinical studies of these molecules in mucosal biology, and in malignant and inflammatory diseases of epithelial tissues.

Global changes in gene expression observed at the transition from growth to stationary phase in Listeria monocytogenes ScottA batch culture

Listeria monocytogenes is a food-borne Gram-positive bacterium that is responsible for a variety of infections (worldwide) annually. The organism is able to survive a variety of environmental conditions and stresses, however, the mechanisms by which L. monocytogenes adapts to environmental change are yet to be fully elucidated. An understanding of the mechanism(s) by which L. monocytogenes survives unfavourable environmental conditions will aid in developing new food processing methods to control the organism in foodstuffs. We have utilized a proteomic approach to investigate the response of L. monocytogenes batch cultures to the transition from exponential to stationary growth phase. Proteomic analysis showed that batch cultures of L. monocytogenes perceived stress and began preparations for stationary phase much earlier (approximately A(600) = 0.75, mid-exponential) than predicted by growth characteristics alone. Global analysis of the proteome revealed that the expression levels of more than 50% of all proteins observed changed significantly over a 7-9 h period during this transition phase. We have highlighted ten proteins in particular whose expression levels appear to be important in the early onset of the stationary phase. The significance of these findings in terms of functionality and the mechanistic picture are discussed.

Insertion of an E-coli lacZ gene in Acetobacter xylinus for the production of cellulose in whey

A mini-Tn10:lacZ: kan was inserted into a wild-type strain of Acetobacter xylinus by random transposon mutagenesis, generating a lactose-utilising and cellulose-producing mutant strain designated ITz3. Antibiotic selection plate assays and Southern hybridisation revealed that the lacZ gene was inserted once into the chromosome of strain ITz3 and was stably maintained in non-selective medium after more than 60 generations. The modified strain had, on the average, a 28-fold increase in cellulose production and a 160-fold increase in beta-galactosidase activity when grown in lactose medium. beta-Galactosidase activity is present in either lactose or sucrose medium indicating that the gene is constitutively expressed. Cellulose and beta-galactosidase production by the modified strain was also evaluated in pure and enriched whey substrates. Utilisation of lactose in whey substrate by ITz3 reached 17 g l(-1) after 4 days incubation. (C) 2004 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.

Correlating gene expression with larval competence and the effect of age and parentage on metamorphosis in the tropical abalone Haliotis asinina

The non-geniculate crustose coralline alga (CCA) Mastophora pacifica can induce the metamorphosis of competent Haliotis asinina (Vetigastropoda) larvae. The ability to respond to this natural cue varies considerably with larval age, with a higher proportion of older larvae (e.g. 90 h) able to metamorphose in response to M. pacifica than younger larvae (e.g. 66 h). Here we document the variation in time to acquisition of competence within a larval age class. For example, after 18 h of exposure to M. pacifica, approximately 15 and 36% of 84 and 90-h-old H. asinina larvae had initiated metamorphosis, respectively. This age-dependent response to M. pacifica is also observed when different aged larvae are exposed to CCA for varying periods. A higher proportion of older larvae require shorter periods of exposure to CCA than younger larvae in order to initiate metamorphosis. In this experiment, as in the previous, a small proportion of young larvae were able to respond to brief periods of CCA exposure, suggesting that they had developed the same state of competency as the majority of their older counterparts. Comparisons of the proportions of larvae undergoing metamorphosis between families reveals that parentage also has a significant (P < 0.05) affect on whether an individual will initiate metamorphosis at a given age. These familial differences are more pronounced when younger, largely pre-competent larvae (i.e. 66 h old) are exposed to M. pacifica, with proportions of larvae undergoing metamorphosis differing by as much as 10 fold between families. As these data suggest that variation in the rate of development of the competent state has a genetic basis, and as a first step towards identifying the molecular basis to this variation, we have identified numerous genes that are differentially expressed later in larval development using a differential display approach. Spatial expression analysis of these genes suggests that they may be directly involved in the acquisition of competence, or may play a functional role in the postlarva following metamorphosis.

Cloning and expression of the large zebrafish protocadherin gene, Fat

The cadherin superfamily members play an important role in mediating cell-cell contact and adhesion (Takeichi, M., 1991. Cadherin cell adhesion receptors as a morphogenetic regulator. Science 251, 1451-1455). A distinct subfamily, neither belonging to the classical or protocadherins includes Fat, the largest member of the cadherin super-family. Fat was originally identified in Drosophila. Subsequently, orthologues of Fat have been described in man (Dunne, J., Hanby, A. M., Poulsom, R., Jones, T. A., Sheer, D., Chin, W. G., Da, S. M., Zhao, Q., Beverley, P. C., Owen, M. J., 1995. Molecular cloning and tissue expression of FAT, the human homologue of the Drosophila fat gene that is located on chromosome 4q34-q35 and encodes a putative adhesion molecule. Genomics 30, 207-223), rat (Ponassi, M., Jacques, T. S., Ciani, L., ffrench, C. C., 1999. Expression of the rat homologue of the Drosophila fat tumour suppressor gene. Mech. Dev. 80, 207-212) and mouse (Cox, B., Hadjantonakis, A. K., Collins, J. E., Magee, A. I., 2000. Cloning and expression throughout mouse development of mfat 1, a homologue of the Drosophila tumour suppressor gene fat [In Process Citation]. Dev. Dyn. 217, 233-240). In Drosophila, Fat has been shown to play an important role in both planar cell polarity and cell boundary formation during development. In this study we describe the characterization of zebrafish Fat, the first non-mammalian, vertebrate Fat homologue to be identified. The Fat protein has 64% amino acid identity and 80% similarity to human FAT and an identical domain structure to other vertebrate Fat proteins. During embryogenesis fat mRNA is expressed in the developing brain, specialised epithelial surfaces the notochord, ears, eyes and digestive tract, a pattern similar but distinct to that found in mammals. (c) 2005 Elsevier B.V. All rights reserved.

Risk of non-small cell lung cancer and the cytochrome P4501A1 Ile462Val polymorphism

Objective: The Ile462Val substitution in the cytochrome P450 1A1 gene (CYP1A1) results in increased enzymatic activity. Preliminary data suggesting a link between this polymorphism and lung cancer risk in Caucasians are inconsistent, reflecting small sample sizes and the relatively low frequency of the variant. Methods: The data set consisted of 1050 primary non-small cell lung cancer cases and 581 controls, a large homogenous population designed specifically to address previous inconsistencies. Patients were genotyped using a PCR-RFLP technique. Results: Carriers of the valine allele, CYP1A1*2C, (Ile/Val or Val/Val genotypes) were significantly over-represented in non-small cell lung cancer compared to controls (OR=1.9; 95% CI=1.2-2.9; p=0.005) when adjusted for confounders, particularly in women (OR=4.6; 95% CI=1.7-12.4; p=0.003). The valine variant was statistically significantly over-represented in cases of lung cancer younger than the median age (64 years) (OR=2.5; 95% CI=1.3-4.8; p=0.005) and cases with less than the median cumulative tobacco-smoke exposure (46 pack-years) (OR=2.4; 95% CI=1.3-4.7; p=0.007). Conclusions: These new data establish an association between the CYP1A1 Ile462Val polymorphism and the risk of developing non-small cell lung cancer, especially among women.

Inherent size constraints on prokaryote gene networks due to "accelerating" growth

Networks exhibiting accelerating growth have total link numbers growing faster than linearly with network size and either reach a limit or exhibit graduated transitions from nonstationary-to-stationary statistics and from random to scale-free to regular statistics as the network size grows. However, if for any reason the network cannot tolerate such gross structural changes then accelerating networks are constrained to have sizes below some critical value. This is of interest as the regulatory gene networks of single-celled prokaryotes are characterized by an accelerating quadratic growth and are size constrained to be less than about 10,000 genes encoded in DNA sequence of less than about 10 megabases. This paper presents a probabilistic accelerating network model for prokaryotic gene regulation which closely matches observed statistics by employing two classes of network nodes (regulatory and non-regulatory) and directed links whose inbound heads are exponentially distributed over all nodes and whose outbound tails are preferentially attached to regulatory nodes and described by a scale-free distribution. This model explains the observed quadratic growth in regulator number with gene number and predicts an upper prokaryote size limit closely approximating the observed value. (c) 2005 Elsevier GmbH. All rights reserved.

The human stress-activated protein kinase-interacting 1 gene encodes JNK-binding proteins

The orthologous proteins of the stress-activated protein kinase-interacting 1 (Sin1) family have been implicated in several different signal transduction pathways. In this study, we have investigated the function of the full-length human Sin1 protein and a C-terminally truncated isoform, Sin 1 alpha, which is produced by alternative splicing. Immunoblot analysis using an anti-Sin 1 polyclonal antibody showed that full-length Sin I and several smaller isoforms are widely expressed. Sin 1 was demonstrated to bind to c-Jun N-terminal kinase (JNK) in vitro and in vivo, while no interaction with p38- or ERK1/2-family MAPKs was observed. The Sin1 alpha isoform could also form a complex with JNK in vivo. Despite localizing in distinct compartments within the cell, both Sin1 and Sin1 alpha co-localized with JNK, suggesting that the Sin1 proteins could recruit JNK. Over-expression of full-length Sin1 inhibited the activation of JNK by UV-C in DG75 cells, as well as basal JNK-activity in HEK293 cells. These data suggest that the human Sin1 proteins may act as scaffold molecules in the regulation of signaling by JNK. (c) 2004 Elsevier Inc. All rights reserved.

Association between polymorphisms in the progesterone receptor gene and endometriosis

The progesterone receptor (PR) is a candidate gene for the development of endometriosis, a complex disease with strong hormonal features, common in women of reproductive age. We typed the 306 base pair Alu insertion (AluIns) polymorphism in intron G of PR in 101 individuals, estimated linkage disequilibrium (LD) between five single-nucleotide polymorphisms (SNPs) across the PR locus in 980 Australian triads (endometriosis case and two parents) and used transmission disequilibrium testing (TDT) for association with endometriosis. The five SNPs showed strong pairwise LD, and the AluIns was highly correlated with proximal SNPs rs1042839 ({Delta}2 = 0.877, D9 = 1.00, P < 0.0001) and rs500760 ({Delta}2 = 0.438, D9 = 0.942, P < 0.0001). TDT showed weak evidence of allelic association between endometriosis and rs500760 (P = 0.027) but not in the expected direction. We identified a common susceptibility haplotype GGGCA across the five SNPs (P = 0.0167) in the whole sample, but likelihood ratio testing of haplotype transmission and non-transmission of the AluIns and flanking SNPs showed no significant pattern. Further, analysis of our results pooled with those from two previous studies suggested that neither the T2 allele of the AluIns nor the T1/T2 genotype was associated with endometriosis.

Experimental validation of the regulated expression of large numbers of non-coding RNAs from the mouse genome

Recent large-scale analyses of mainly full-length cDNA libraries generated from a variety of mouse tissues indicated that almost half of all representative cloned sequences did flat contain ail apparent protein-coding sequence, and were putatively derived from non-protein-coding RNA (ncRNA) genes. However, many of these clones were singletons and the majority were unspliced, raising the possibility that they may be derived from genomic DNA or unprocessed pre-rnRNA contamination during library construction, or alternatively represent nonspecific transcriptional noise. Here we Show, using reverse transcriptase-dependent PCR, microarray, and Northern blot analyses, that many of these clones were derived from genuine transcripts Of unknown function whose expression appears to be regulated. The ncRNA transcripts have larger exons and fewer introns than protein-coding transcripts. Analysis of the genomic landscape around these sequences indicates that some cDNA clones were produced not from terminal poly(A) tracts but internal priming sites within longer transcripts, only a minority of which is encompassed by known genes. A significant proportion of these transcripts exhibit tissue-specific expression patterns, as well as dynamic changes in their expression in macrophages following lipopolysaccharide Stimulation. Taken together, the data provide strong support for the conclusion that ncRNAs are an important, regulated component of the mammalian transcriptome.