Integrating sequence information in microarray data analysis by free energy modeling

Microarray technology is a high-throughput method for genotyping and gene expression profiling. Limited sensitivity and specificity are one of the essential problems for this technology. Most of existing methods of microarray data analysis have an apparent limitation for they merely deal with the numerical part of microarray data and have made little use of gene sequence information. Because it's the gene sequences that precisely define the physical objects being measured by a microarray, it is natural to make the gene sequences an essential part of the data analysis. This dissertation focused on the development of free energy models to integrate sequence information in microarray data analysis. The models were used to characterize the mechanism of hybridization on microarrays and enhance sensitivity and specificity of microarray measurements. ^ Cross-hybridization is a major obstacle factor for the sensitivity and specificity of microarray measurements. In this dissertation, we evaluated the scope of cross-hybridization problem on short-oligo microarrays. The results showed that cross hybridization on arrays is mostly caused by oligo fragments with a run of 10 to 16 nucleotides complementary to the probes. Furthermore, a free-energy based model was proposed to quantify the amount of cross-hybridization signal on each probe. This model treats cross-hybridization as an integral effect of the interactions between a probe and various off-target oligo fragments. Using public spike-in datasets, the model showed high accuracy in predicting the cross-hybridization signals on those probes whose intended targets are absent in the sample. ^ Several prospective models were proposed to improve Positional Dependent Nearest-Neighbor (PDNN) model for better quantification of gene expression and cross-hybridization. ^ The problem addressed in this dissertation is fundamental to the microarray technology. We expect that this study will help us to understand the detailed mechanism that determines sensitivity and specificity on the microarrays. Consequently, this research will have a wide impact on how microarrays are designed and how the data are interpreted. ^

Heterotrimeric G protein -mediated signal transduction in Neurospora crassa

Heterotrimeric G protein-mediated signal transduction is one of numerous means that cells utilize to respond to external stimuli. G proteins consist of α, β andγ subunits. Extracellular ligands bind to seven-transmembrane helix receptors, triggering conformational changes. This is followed by activation of coupled G proteins through the exchange of GDP for GTP on the Gα subunit. Once activated, Gα-GTP dissociates from the βγ dimer. Both of these two moieties can interact with downstream effectors, such as adenylyl cyclase, phospholipase C, phosphodiesterases, or ion channels, leading to a series of changes in cellular metabolism and physiology. ^ Neurospora crassa is a eukaryotic multicellular filamentous fungus, with asexual/vegetative and sexual phases to its life cycle. Three Gα (GNA-1, GNA-2, GNA-3) and one Gβ (GNB-1) proteins have been identified in this organism. This dissertation investigates GNA-1 and GNB-1 mediated signaling pathways in N. crassa. ^ GNA-1 was the first identified microbial Gα that belongs to a mammalian superfamily (Gαi). Deletion of GNA-1 leads to multiple defects in N. crassa. During the asexual cycle, Δgna-1 strains display a slower growth rate and delayed conidiation on solid medium. In the sexual cycle, the Δgna-1 mutant is male-fertile but female-sterile. Biochemical studies have shown that Δ gna-1 strains have lower adenosine 3′–5 ′ cyclic monophosphate (cAMP) levels than wild type under conditions where phenotypic defects are observed. In this thesis work, strains containing one of two GTPase-deficient gna-1 alleles (gna-1 R178C, gna-1Q204L) leading to constitutive activation of GNA-1 have been constructed and characterized. Activation of GNA-1 causes uncontrolled aerial hyphae proliferation, elevated sensitivity to heat and oxidative stresses, and lower carotenoid synthesis. To further study the function of GNA-1, constructs to enable expression of mammalian Gαi superfamily members were transformed into a Δ gna-1 strain, and complementation of Δgna-1 defects investigated. Gαs, which is not a member of Gα i superfamily was used as a control. These mammalian Gα genes were able to rescue the vegetative growth rate defect of the Δ gna-1 strain in the following order: Gαz > Gα o > Gαs > Gαt > Gαi. In contrast, only Gαo was able to complement the sexual defect of a Δgna-1 strain. With regard to the thermotolerance phenotype, none of the mammalian Gα genes restored the sensitivity to a wild type level. These results suggest that GNA-1 regulates two independent pathways during the vegetative and sexual cycles in N. crassa. ^ GNB-1, a G protein β subunit from N. crassa, was identified and its functions investigated in this thesis work. The sequence of the gnb-1 gene predicts a polypeptide of 358 residues with a molecular mass of 39.7 kDa. GNB-1 exhibits 91% identity to Cryphonectria parasitica CPGB-1, and also displays significant homology with human and Dictyostelium Gβ genes (∼66%). A Δ gnb-1 strain was constructed and shown to exhibit defects in asexual spore germination, vacuole number and size, mass accumulation and female fertility. A novel role for GNB-1 in regulation of GNA-1 and GNA-2 protein levels was also demonstrated. ^

(Table 2) Annual mass fluxes, 230Th composition and organic geochemistry of sediment trap samples in the South Altantic

The scavenging of 231Pa and 230Th was investigated in the Atlantic Sector of the Southern Ocean by combining results from sediment trap and in situ filtration studies. We present the first high-resolution profile of dissolved 230Th and 231Pa in surface waters across the ACC, showing a dramatic southward increase of both radionuclides around the southern ACC Front at 533S. High dissolved 231Pa/230Th ratios combined with low 230Th/231Pa fractionation factors (F) in these surface waters result in extremely high 231Pa94/230Th94 ratios of material collected in the shallow traps. Particulate 231Pa94/230Th94 ratios in a shallow trap near Bouvet Island increase continuously during the productive period in austral summer, and drop back in the low flux period. This behavior, following the Rayleigh fractionation principle, is interpreted to be due to an increase in the dissolved 231Pa/230Th ratio in the euphotic zone resulting from preferential scavenging of 230Th relative to 231Pa, even in opal-dominated regions. In the post-bloom stage, the depleted radionuclide concentrations are replenished by upwelling of Circumpolar Deep Water. The high particulate 231Pa94/230Th94 signal is weakened during downward transport of the bloom particles in the water column by incorporation of deep suspended particles, which have a lower 231Pa94/230Th94 ratio. It is shown that under the special hydrographic conditions in the Southern Ocean scavenging from the upper water column significantly influences the budgets of 230Th and 231Pa in the sediment. Nevertheless, the budgets are still made up primarily by scavenging from the large standing stock of deep suspended particles.

TEX86 values and temperature estimation for sediment trap samples

The newly introduced temperature proxy, the tetraether index of archaeal lipids with 86 carbon atoms (TEX86), is based on the number of cyclopentane moieties in the glycerol dialkyl glycerol tetraether (GDGT) lipids of marine Crenarchaeota. The composition of sedimentary GDGTs used for TEX86 paleothermometry is thought to reflect sea surface temperature (SST). However, marine Crenarchaeota occur ubiquitously in the world oceans over the entire depth range and not just in surface waters. We analyzed the GDGT distribution in settling particulate organic matter collected in sediment traps from the northeastern Pacific Ocean and the Arabian Sea to investigate the seasonal and spatial distribution of the fluxes of crenarchaeotal GDGTs and the origin of the TEX86 signal transported to the sediment. In both settings the TEX86 measured at all trap deployment depths reflects SST. In the Arabian Sea, analysis of an annual time series showed that the SST estimate based on TEX86 in the shallowest trap at 500 m followed the in situ SST with a 1 to 3 week time delay, likely caused by the relatively low settling speed of sinking particles. This revealed that the GDGT signal that reaches deeper water is derived from the upper water column rather than in situ production of GDGTs. The GDGT temperature signal in deeper traps at 1500 m and 3000 m did not show a seasonal cyclicity observed in the 500 m trap but rather reflected the annual mean SST. This is probably due to a homogenization of the TEX86 SST signal carried by particles as they ultimately reach the interior of the ocean. Our data confirm the use of TEX86 as a temperature proxy of surface ocean waters.

d15N of trap WR2

Temporal changes in d15N values of sinking particles collected with sediment traps in the Benguela upwelling regime off southwest Africa mirrored variations in the input of inorganic nitrogen to the surface water. Reductions in d15N (to as low as 2.5 per mil) corresponded to low sea surface temperatures during austral spring and late austral autumn/early winter, indicating increased nitrate availability due to the presence of recently upwelled water. High particulate fluxes accompanied the low d15N values and sea surface temperatures, reflecting increased productivity, fueled by the upwelled nutrients. High d15N values (up to 13.1 per mil) coincided with high sea surface temperatures and low particle fluxes. In this area, the seaward extension of upwelling filaments, which usually occurs twice yearly, brings nutrient-rich water to the euphotic zone and leads to elevated productivity and relatively lower d15N values of the particulate nitrogen. Satellite images of ocean chlorophyll show that productivity variations coincide with d15N changes. The observed isotopic pattern does not appear to have been caused by variations in the species composition of the phytoplankton assemblage. Calculations based on d15N of the sinking particulate nitrogen show that the surface nitrate pool was more depleted during late austral summer/early fall and mid-winter and that supply exceeded demand during the intense spring bloom and in late austral fall. The main uncertainty associated with these estimates is the effect of diagenesis on d15N and possible variability in preservation of the isotope signal between periods of high and low particle flux.

Stable isotope, seismic sequence boundaries and bioevents in ODP Hole 166-1006A

During ODP Leg 166, the recovery of cores from a transect of drill sites across the Bahamas margin from marginal to deep basin environments was an essential requirement for the study of the response of the sedimentary systems to sea-level changes. A detailed biostratigraphy based on planktonic foraminifera was performed on ODP Hole 1006A for an accurate stratigraphic control. The investigated late middle Miocene-early Pliocene sequence spans the interval from about 12.5 Ma (Biozone N12) to approximately 4.5 Ma (Biozone N19). Several bioevents calibrated with the time scale of Berggren et al. (1995a,b) were identified. The ODP Site 1006 benthic oxygen isotope stratigraphy can be correlated to the corresponding deep-water benthic oxygen isotope curve from ODP Site 846 in the Eastern Equatorial Pacific (Shackleton et al., 1995. Proc. ODP Sci. Res. 138, 337-356), which was orbitally tuned for the entire Pliocene into the latest Miocene at 6.0 Ma. The approximate stratigraphic match of the isotopic signals from both records between 4.5 and 6.0 Ma implies that the paleoceanographic signal from the Bahamas is not simply a record of regional variations but, indeed, represents glacio-eustatic fluctuations. The ODP Site 1006 oxygen and carbon isotope record, based on benthic and planktonic foraminifera, was used to define paleoceanographic changes on the margin, which could be tied to lithostratigraphic events on the Bahamas carbonate platform using seismic sequence stratigraphy. The oxygen isotope values show a general cooling trend from the middle to late Miocene, which was interrupted by a significant trend towards warmer sea-surface temperatures (SST) and associated sea-level rise with decreased ice volume during the latest Miocene. This trend reached a maximum coincident with the Miocene/Pliocene boundary. An abrupt cooling in the early Pliocene then followed the warming which continued into the earliest Pliocene. The late Miocene paleoceanographic evolution along the Bahamas margin can be observed in the ODP Site 1006 delta13C values, which support other evidence for the beginning of the closure of the Panama gateway at 8 Ma followed by a reduced intermediate water supply of water from the Pacific into the Caribbean at about 5 Ma. A general correlation of lower sedimentation rates with the major seismic sequence boundaries (SSBs) was observed. Additionally, the SSBs are associated with transitions towards more positive oxygen isotope excursions. This observed correspondence implies that the presence of a SSB, representing a density impedance contrast in the sedimentary sequence, may reflect changes in the character of the deposited sediment during highstands versus those during lowstands. However, not all of the recorded oxygen isotope excursions correspond to SSBs. The absence of a SSB in association with an oxygen isotope excursion indicates that not all oxygen isotope sea-level events impact the carbonate margin to the same extent, or maybe even represent equivalent sea-level fluctuations. Thus, it can be tentatively concluded that SSBs produced on carbonate margins do record sea-level fluctuations but not every sea-level fluctuation is represented by a SSB in the sequence stratigraphic record.

Gene discovery in the wood-forming tissues of poplar: Analysis of 5,692 expressed sequence tags

A rapidly growing area of genome research is the generation of expressed sequence tags (ESTs) in which large numbers of randomly selected cDNA clones are partially sequenced. The collection of ESTs reflects the level and complexity of gene expression in the sampled tissue. To date, the majority of plant ESTs are from nonwoody plants such as Arabidopsis, Brassica, maize, and rice. Here, we present a large-scale production of ESTs from the wood-forming tissues of two poplars, Populus tremula L. × tremuloides Michx. and Populus trichocarpa ‘Trichobel.’ The 5,692 ESTs analyzed represented a total of 3,719 unique transcripts for the two cDNA libraries. Putative functions could be assigned to 2,245 of these transcripts that corresponded to 820 protein functions. Of specific interest to forest biotechnology are the 4% of ESTs involved in various processes of cell wall formation, such as lignin and cellulose synthesis, 5% similar to developmental regulators and members of known signal transduction pathways, and 2% involved in hormone biosynthesis. An additional 12% of the ESTs showed no significant similarity to any other DNA or protein sequences in existing databases. The absence of these sequences from public databases may indicate a specific role for these proteins in wood formation. The cDNA libraries and the accompanying database are valuable resources for forest research directed toward understanding the genetic control of wood formation and future endeavors to modify wood and fiber properties for industrial use.

Identification and analysis of the plant peroxisomal targeting signal 1 receptor NtPEX5

Protein translocation into peroxisomes takes place via recognition of a peroxisomal targeting signal present at either the extreme C termini (PTS1) or N termini (PTS2) of matrix proteins. In mammals and yeast, the peroxisomal targeting signal receptor, Pex5p, recognizes the PTS1 consisting of -SKL or variants thereof. Although many plant peroxisomal matrix proteins are transported through the PTS1 pathway, little is known about the PTS1 receptor or any other peroxisome assembly protein from plants. We cloned tobacco (Nicotiana tabacum) cDNAs encoding Pex5p (NtPEX5) based on the protein’s interaction with a PTS1-containing protein in the yeast two-hybrid system. Nucleotide sequence analysis revealed that the tobacco Pex5p contains seven tetratricopeptide repeats and that NtPEX5 shares greater sequence similarity with its homolog from humans than from yeast. Expression of NtPEX5 fusion proteins, consisting of the N-terminal part of yeast Pex5p and the C-terminal region of NtPEX5, in a Saccharomyces cerevisiae pex5 mutant restored protein translocation into peroxisomes. These experiments confirmed the identity of the tobacco protein as a PTS1 receptor and indicated that components of the peroxisomal translocation apparatus are conserved functionally. Two-hybrid assays showed that NtPEX5 interacts with a wide range of PTS1 variants that also interact with the human Pex5p. Interestingly, the C-terminal residues of some of these peptides deviated from the established plant PTS1 consensus sequence. We conclude that there are significant sequence and functional similarities between the plant and human Pex5ps.

The protein kinases of Caenorhabditis elegans: A model for signal transduction in multicellular organisms

Caenorhabditis elegans should soon be the first multicellular organism whose complete genomic sequence has been determined. This achievement provides a unique opportunity for a comprehensive assessment of the signal transduction molecules required for the existence of a multicellular animal. Although the worm C. elegans may not much resemble humans, the molecules that regulate signal transduction in these two organisms prove to be quite similar. We focus here on the content and diversity of protein kinases present in worms, together with an assessment of other classes of proteins that regulate protein phosphorylation. By systematic analysis of the 19,099 predicted C. elegans proteins, and thorough analysis of the finished and unfinished genomic sequences, we have identified 411 full length protein kinases and 21 partial kinase fragments. We also describe 82 additional proteins that are predicted to be structurally similar to conventional protein kinases even though they share minimal primary sequence identity. Finally, the richness of phosphorylation-dependent signaling pathways in worms is further supported with the identification of 185 protein phosphatases and 128 phosphoprotein-binding domains (SH2, PTB, STYX, SBF, 14-3-3, FHA, and WW) in the worm genome.

Molecular uncoupling of C-C chemokine receptor 5-induced chemotaxis and signal transduction from HIV-1 coreceptor activity

The C-C chemokine receptor 5 (CCR5) plays a crucial role in facilitating the entry of macrophage-tropic strains of the HIV-1 into cells, but the mechanism of this phenomenon is completely unknown. To explore the role of CCR5-derived signal transduction in viral entry, we introduced mutations into two cytoplasmic domains of CCR5 involved in receptor-mediated function. Truncation of the terminal carboxyl-tail to eight amino acids or mutation of the highly conserved aspartate-arginine-tyrosine, or DRY, sequence in the second cytoplasmic loop of CCR5 effectively blocked chemokine-dependent activation of classic second messengers, intracellular calcium fluxes, and the cellular response of chemotaxis. In contrast, none of the mutations altered the ability of CCR5 to act as an HIV-1 coreceptor. We conclude that the initiation of signal transduction, the prototypic function of G protein coupled receptors, is not required for CCR5 to act as a coreceptor for HIV-1 entry into cells.

Inhibition of mRNA export in vertebrate cells by nuclear export signal conjugates

Leucine-rich nuclear export signals (NESs) are recognized by the NES receptor exportin 1 and are central to the export of multiple shuttling proteins and RNAs. The export of messenger RNA in vertebrates was, however, thought to occur by a different pathway, because inhibition by injection of a synthetic Rev NES conjugate could not be demonstrated. Here we find that peptide conjugates composed of the NES of either protein kinase A inhibitor protein (PKI) or the HIV-1 Rev protein, when coupled to human serum albumin, are potent inhibitors of mRNA and small nuclear RNA export. These results provide direct evidence that mRNA export in vertebrates depends on interactions between an NES and its cognate NES receptors. PKI NES conjugates are significantly more efficient at inhibiting RNA export than are REV NES conjugates, indicating that different NESs may have different abilities to promote protein and RNA export. Surprisingly, an expected control conjugate containing the mutant Rev NES sequence M10 strongly inhibited the export of intronless dihydrofolate reductase mRNA. Nuclear injection of NES peptide conjugates led to mislocalization to the nucleus of 10–20% of the cytoplasmic Ran GTPase-binding protein (RanBP1) indicating that RanBP1 shuttles between the nucleus and the cytoplasm via an NES pathway. These results demonstrate that in vertebrates the export of mRNA, like that of small nuclear RNA, 5S rRNA, and transport factors such as RanBP1, employs NES-mediated molecular machinery.

Genomic and cDNA sequence analysis of the cell matrix adhesion regulator gene

The cell matrix adhesion regulator (CMAR) gene has been suggested to be a signal transduction molecule influencing cell adhesion to collagen and, through this, possibly involved in tumor suppression. The originally reported CMAR cDNA was 464 bp long with a tyrosine phosphorylation site at the extreme 3′ end, which mutagenesis studies had shown to be central to the function of this gene. Since the discovery of a 4-bp insertion polymorphism within the originally reported coding region, further sequence information has been obtained. The cDNA has been extended 5′ by ≈2 kb revealing a 559-bp region showing strong homology to the proposed 5′ untranslated sequence of a murine protein kinase receptor family member, variant in kinase (vik). CMAR genomic sequencing has shown the presence of an intron, the intron/exon boundary lying within this region of homology. An RNA transcript for CMAR of ≈2.5 kb has also been identified. The data suggest complex mechanisms for control of expression of two closely associated genes, CMAR and the vik- associated sequence.

Efficient gene tagging in Arabidopsis thaliana using a gene trap approach

Large quantities of DNA sequence information about plant genes are rapidly accumulating in public databases, but to progress from DNA sequence to biological function a mutant allele for each of the genes ideally should be available. Here we describe a gene trap construct that allowed us to disrupt transcribed genes with a high efficiency in Arabidopsis thaliana. In the T-DNA vector used, the expression of a bacterial reporter gene coding for neomycin phosphotransferase II (nptII) depends on the in vivo generation of a translation fusion upon the T-DNA integration into the Arabidopsis genome. Analysis of 20 selected transgenic lines showed that 12 lines are T-DNA insertion mutants. The disrupted genes analyzed encoded ribosomal proteins (three lines), aspartate tRNA synthase, DNA ligase, basic-domain leucine zipper DNA binding protein, ATP-binding cassette transporter, and five proteins of unknown function. Four tagged genes were new for Arabidopsis. The results presented here suggest that gene trapping, using nptII as a reporter gene, can be as high as 80% and opens novel perspectives for systematic gene tagging in A. thaliana.