Biodiversity and Phylogeny of Planktic Cyanobacteria in Temperate Freshwater Lakes

Currently, the classification used for cyanobacteria is based mainly on morphology. In many cases the classification is known to be incongruent with the phylogeny of cyanobacteria. The evaluation of this classification is complicated by the fact that numerous strains are only described morphologically and have not been isolated. Moreover, the phenotype of many cyanobacterial strains alters during prolonged laboratory cultivation. In this thesis, cyanobacterial strains were isolated from lakes (mainly Lake Tuusulanjärvi) and both morphology and phylogeny of the isolates were investigated. The cyanobacterial community composition in Lake Tuusulanjärvi was followed for two years in order to relate the success of cyanobacterial phenotypes and genotypes to environmental conditions. In addition, molecular biological methods were compared with traditional microscopic enumeration and their ability and usefulness in describing the cyanobacterial diversity was evaluated. The Anabaena, Aphanizomenon, and Trichormus strains were genetically heterogeneous and polyphyletic. The phylogenetic relationships of the heterocytous cyanobacteria were not congruent with their classification. In contrast to heterocytous cyanobacteria, the phylogenetic relationships of the Snowella and Woronichinia strains, which had not been studied before this thesis, reflected the morphology of strains and followed their current classification. The Snowella strains formed a monophyletic cluster, which was most closely related to the Woronichinia strain. In addition, a new cluster of thin, filamentous cyanobacterial strains identified as Limnothrix redekei was revealed. This cluster was not closely related to any other known cyanobacteria. The cyanobacterial community composition in Lake Tuusulanjärvi was studied with molecular methods [denaturant gradient gel electrophoresis (DGGE) and cloning of the 16S rRNA gene], through enumerations of cyanobacteria under microscope, and by strain isolations. Microcystis, Anabaena/Aphanizomenon, and Synechococcus were the major groups in the cyanobacterial community in Lake Tuusulanjärvi during the two-year monitoring period. These groups showed seasonal succession, and their success was related to different environmental conditions. The major groups of the cyanobacterial community were detected by all used methods. However, cloning gave higher estimates than microscopy for the proportions of heterocytous cyanobacteria and Synechococcus. The differences were probably caused by the high 16S rRNA gene copy numbers in heterotrophic cyanobacteria and by problems in the identification and detection of unicellular cyanobacteria.

Stanniocalcin-1 in cell stress and differentiation

Stanniocalcin-1 (STC-1) is a 56 kD homodimeric protein which was originally identified in bony fish, where it regulates calcium/phosphate homeostasis and protects against toxic hypercalcemia. STC-1 was considered unique to fish until the cloning of cDNA for human STC-1 in 1995 and mouse Stc-1 in 1996. STC-1 is conserved through evolution with human and salmon STC-1 sharing 60% identity and 80% similarity. The surprisingly high homology between mammalian and fish STC-1 and the protective actions of STC-1 in terminally differentiated neurons, originally reported by my colleagues, prompted me to further study the role of STC-1 in cell stress and differentiation. One purpose was to determine whether there is an inter-relationship between terminally differentiated cells and STC-1 expression. The study revealed an accumulation of STC-1 in mature megakaryocytes and adipocytes, i.e. postmitotic cells with limited or lost proliferative capacity. Still proliferating uninduced cells were negative for STC-1 mRNA and protein, whereas differentiating cells accumulated STC-1 in their cytoplasm. Interestingly, in liposarcomas the grade inversely correlated with STC-1 expression. Another aim was to study how STC-1 gene expression is regulated. Given that IL-6 is a cytokine with neuroprotective actions, by unknown mechanisms, we examined whether IL-6 regulates STC-1 gene expression. Treatment of human neural Paju cells with IL-6 induced a dose-dependent upregulation of STC-1 mRNA levels. This induction of STC-1 expression by IL-6 occurred mainly through the MAPK signaling pathway. Furthermore, I studied the role of IL-6-mediated STC-1 expression as a mechanism of cytoprotection conferred by hypoxic preconditioning (HOPC) in brain and heart. My findings show that Stc-1 was upregulated in brain after hypoxia treatment. In the brain of IL-6 deficient mice, however, no upregulation of Stc-1 expression was evident. After induced brain injury the STC-1 response in brains of IL-6 transgenic mice, with IL-6 overexpression in astroglial cells, was stronger than in brains of WT mice. These results indicate that IL-6-mediated expression of STC-1 is one molecular mechanism of HOPC-induced tolerance to brain ischemia. The protection conferred by HOPC in heart occurs during a bimodal time course comprising early and delayed preconditioning. Interestingly, my results showed that the expression of Stc-1 in heart was upregulated in a biphasic manner during HOPC. IL-6 deficient mice did not, however, show a similar biphasic manner of Stc-1 upregulation as did WT mice. Instead, only an early upregulation of Stc-1 expression was evident. The results suggest that the upregulation of Stc-1 during the delayed preconditioning is IL-6-dependent. The upregulated expression of Stc-1 during the early preconditioning, however, is only partly IL-6-dependent and possibly also directly mediated by HIF-1. These findings suggest that STC-1 is a pro-survival protein for terminally differentiated cells and that STC-1 expression may in fact be regulated by stress. In addition, I show that STC-1 gene upregulation, mediated in part by IL-6, is a new mechanism of protection conferred by HOPC in brain and heart. Because of its importance for fundamental biological processes, such as differentiation and cytoprotection, STC-1 may have therapeutic implications for management of stroke, neurodegenerative diseases, cancer, and obesity.

GDP-L-fucose: synthesis and role in inflammation

GDP-L-fucose: synthesis and role in inflammation The migration of leukocytes from intravascular locations to extravascular sites is essential to the immune responses. The initial attachment of leukocytes to the endothelium and the rolling step of the leukocyte extravasation cascade are mediated by selectins, a family of cell adhesion molecules on cell surfaces. Selectins are able to recognize glycoproteins and glycolipids containing the tetrasaccharide sialyl Lewis x (sLex, Neu5Acα2-3Galβ1-4(Fucα1-3)GlcNAc). Several glycosyltransferases are involved in the biosynthesis of sLex, fucosyltransferase VII (Fuc-TVII) being the last enzyme to modify the sLex structure. Fuc-TVII transfers L-fucose from GDP-L-fucose to sialylated N-acetyllactosamine. GDP-L-fucose is synthesized in the cytosol via two different metabolic pathways. The major, constitutively active de novo pathway involves conversion of GDP-α-D-mannose to GDP-β-L-fucose. In the alternative salvage pathway, L-fucokinase synthesizes from free fucose L-fucose-1-phosphate, which is further converted to GDP-L-fucose by GDP-L-fucose pyrophosphorylase. GDP-L-fucose is translocated from the cytosol to Golgi for fucosylation via the GDP-fucose transporter. This thesis involved the study of the synthesis of GDP-L-fucose via the salvage pathway: cloning and expression of murine L-fucokinase and GDP-L-fucose pyrophosphorylase. The gene expression levels of these enzymes were found to be relatively high in various tissues; the mRNA levels were highest in brain, ovary and testis. This study also describes molecular cloning of rat fucosyltransferase VII (FUT7) and its expression as a functional enzyme. Gene expression levels of GDP-L-fucose synthesizing enzymes, GDP-fucose transporter and FUT7 were determined in inflamed tissues as well as cancer cells. Our results revealed a clear upregulation of the enzymes involved in the synthesis of GDP-L-fucose via de novo pathway, GDP-fucose transporter and FUT7 in inflamed tissues and in cancer cells. On the contrary, the GDP-L-fucose salvage pathway was found to be irrelevant in inflammation and in tumorigenesis. Furthermore, our results indicated the transcriptional coregulation of Golgi transporters involved in the synthesis of sulfo sLex, i.e. CMP-sialic acid, GDP-fucose and 3 phosphoadenosine 5 -phosphosulfate transporters, in inflammation.

The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.

Characterisation, cloning and production of industrially interesting enzymes: gluconolactone oxidase of Penicillium cyaneo-fulvum and gluconate 5-dehydrogenase of Gluconobacter suboxydans

The work covered in this thesis is focused on the development of technology for bioconversion of glucose into D-erythorbic acid (D-EA) and 5-ketogluconic acid (5-KGA). The task was to show on proof-of-concept level the functionality of the enzymatic conversion or one-step bioconversion of glucose to these acids. The feasibility of both studies to be further developed for production processes was also evaluated. The glucose - D-EA bioconversion study was based on the use of a cloned gene encoding a D-EA forming soluble flavoprotein, D-gluconolactone oxidase (GLO). GLO was purified from Penicillium cyaneo-fulvum and partially sequenced. The peptide sequences obtained were used to isolate a cDNA clone encoding the enzyme. The cloned gene (GenBank accession no. AY576053) is homologous to the other known eukaryotic lactone oxidases and also to some putative prokaryotic lactone oxidases. Analysis of the deduced protein sequence of GLO indicated the presence of a typical secretion signal sequence at the N-terminus of the enzyme. No other targeting/anchoring signals were found, suggesting that GLO is the first known lactone oxidase that is secreted rather than targeted to the membranes of the endoplasmic reticulum or mitochondria. Experimental evidence supports this analysis, as near complete secretion of GLO was observed in two different yeast expression systems. Highest expression levels of GLO were obtained using Pichia pastoris as an expression host. Recombinant GLO was characterised and the suitability of purified GLO for the production of D-EA was studied. Immobilised GLO was found to be rapidly inactivated during D-EA production. The feasibility of in vivo glucose - D-EA conversion using a P. pastoris strain co-expressing the genes of GLO and glucose oxidase (GOD, E.C. 1.1.3.4) of A. niger was demonstrated. The glucose - 5-KGA bioconversion study followed a similar strategy to that used in the D-EA production research. The rationale was based on the use of a cloned gene encoding a membrane-bound pyrroloquinoline quinone (PQQ)-dependent gluconate 5-dehydrogenase (GA 5-DH). GA 5-DH was purified to homogeneity from the only source of this enzyme known in literature, Gluconobacter suboxydans, and partially sequenced. Using the amino acid sequence information, the GA 5-DH gene was cloned from a genomic library of G. suboxydans. The cloned gene was sequenced (GenBank accession no. AJ577472) and found to be an operon of two adjacent genes encoding two subunits of GA 5-DH. It turned out that GA 5-DH is a rather close homologue of a sorbitol dehydrogenase from another G. suboxydans strain. It was also found that GA 5-DH has significant polyol dehydrogenase activity. The G. suboxydans GA 5-DH gene was poorly expressed in E. coli. Under optimised conditions maximum expression levels of GA 5-DH did not exceed the levels found in wild-type G. suboxydans. Attempts to increase expression levels resulted in repression of growth and extensive cell lysis. However, the expression levels were sufficient to demonstrate the possibility of bioconversion of glucose and gluconate into 5-KGA using recombinant strains of E. coli. An uncharacterised homologue of GA 5-DH was identified in Xanthomonas campestris using in silico screening. This enzyme encoded by chromosomal locus NP_636946 was found by a sequencing project of X. campestris and named as a hypothetical glucose dehydrogenase. The gene encoding this uncharacterised enzyme was cloned, expressed in E. coli and found to encode a gluconate/polyol dehydrogenase without glucose dehydrogenase activity. Moreover, the X. campestris GA 5-DH gene was expressed in E. coli at nearly 30 times higher levels than the G. suboxydans GA 5-DH gene. Good expressability of the X. campestris GA-5DH gene makes it a valuable tool not only for 5-KGA production in the tartaric acid (TA) bioprocess, but possibly also for other bioprocesses (e.g. oxidation of sorbitol into L-sorbose). In addition to glucose - 5-KGA bioconversion, a preliminary study of the feasibility of enzymatic conversion of 5-KGA into TA was carried out. Here, the efficacy of the first step of a prospective two-step conversion route including a transketolase and a dehydrogenase was confirmed. It was found that transketolase convert 5-KGA into TA semialdehyde. A candidate for the second step was suggested to be succinic dehydrogenase, but this was not tested. The analysis of the two subprojects indicated that bioconversion of glucose to TA using X. campestris GA 5-DH should be prioritised first and the process development efforts in future should be focused on development of more efficient GA 5-DH production strains by screening a more suitable production host and by protein engineering.

The white-rot fungi Phlebia radiata and Dichomitus squalens in wood-based cultures: expression of laccases, lignin peroxidases, and oxalate decarboxylase

Basidiomycetous white-rot fungi are the only organisms that can efficiently decompose all the components of wood. Moreover, white-rot fungi possess the ability to mineralize recalcitrant lignin polymer with their extracellular, oxidative lignin-modifying enzymes (LMEs), i.e. laccase, lignin peroxidase (LiP), manganese peroxidase (MnP), and versatile peroxidase (VP). Within one white-rot fungal species LMEs are typically present as several isozymes encoded by multiple genes. This study focused on two effi cient lignin-degrading white-rot fungal species, Phlebia radiata and Dichomitus squalens. Molecular level knowledge of the LMEs of the Finnish isolate P. radiata FBCC43 (79, ATCC 64658) was complemented with cloning and characterization of a new laccase (Pr-lac2), two new LiP-encoding genes (Pr-lip1, Pr-lip4), and Pr-lip3 gene that has been previously described only at cDNAlevel. Also, two laccase-encoding genes (Ds-lac3, Ds-lac4) of D. squalens were cloned and characterized for the first time. Phylogenetic analysis revealed close evolutionary relationships between the P. radiata LiP isozymes. Distinct protein phylogeny for both P. radiata and D. squalens laccases suggested different physiological functions for the corresponding enzymes. Supplementation of P. radiata liquid culture medium with excess Cu2+ notably increased laccase activity and good fungal growth was achieved in complex medium rich with organic nitrogen. Wood is the natural substrate of lignin-degrading white-rot fungi, supporting production of enzymes and metabolites needed for fungal growth and the breakdown of lignocellulose. In this work, emphasis was on solid-state wood or wood-containing cultures that mimic the natural growth conditions of white-rot fungi. Transcript analyses showed that wood promoted expression of all the presently known LME-encoding genes of P. radiata and laccase-encoding genes of D. squalens. Expression of the studied individual LME-encoding genes of P. radiata and D. squalens was unequal in transcript quantities and apparently time-dependent, thus suggesting the importance of several distinct LMEs within one fungal species. In addition to LMEs, white-rot fungi secrete other compounds that are important in decomposition of wood and lignin. One of these compounds is oxalic acid, which is a common metabolite of wood-rotting fungi. Fungi produce also oxalic-acid degrading enzymes of which the most widespread is oxalate decarboxylase (ODC). However, the role of ODC in fungi is still ambiguous with propositions from regulation of intra and extracellular oxalic acid levels to a function in primary growth and concomitant production of ATP. In this study, intracellular ODC activity was detected in four white-rot fungal species, and D. squalens showed the highest ODC activity upon exposure to oxalic acid. Oxalic acid was the most common organic acid secreted by the ODC-positive white-rot fungi and the only organic acid detected in wood cultures. The ODC-encoding gene Ds-odc was cloned from two strains of D. squalens showing the first characterization of an odc-gene from a white-rot polypore species. Biochemical properties of the D. squalens ODC resembled those described for other basidiomycete ODCs. However, the translated amino acid sequence of Ds-odc has a novel N-terminal primary structure with a repetitive Ala-Ser-rich region of ca 60 amino acid residues in length. Expression of the Ds-odc transcripts suggested a constitutive metabolic role for the corresponding ODC enzyme. According to the results, it is proposed that ODC may have an essential implication for the growth and basic metabolism of wood-decaying fungi.

From actin monomers to bundles: The roles of twinfilin and α-actinin in Drosophila melanogaster development

The actin cytoskeleton is essential for a large variety of cell biological processes. Actin exists in either a monomeric or a filamentous form, and it is very important for many cellular functions that the local balance between these two actin populations is properly regulated. A large number of proteins participate in the regulation of actin dynamics in the cell, and twinfilin, one of the proteins examined in this thesis, belongs to this category. The second level of regulation involves proteins that crosslink or bundle actin filaments, thereby providing the cell with a certain shape. α-Actinin, the second protein studied, mainly acts as an actin crosslinking protein. Both proteins are conserved in organisms ranging from yeast to mammals. In this thesis, the roles of twinfilin and α-actinin in development were examined using Drosophila melanogaster as a model organism. Twinfilin is an actin monomer binding protein that is structurally related to cofilin. In vitro, twinfilin reduces actin polymerisation by sequestering actin monomers. The Drosophila twinfilin (twf) gene was identified and found to encode a protein functionally similar to yeast and mammalian twinfilins. A strong hypomorphic twf mutation was identified, and flies homozygous for this allele were viable and fertile. The adult twf mutant flies displayed reduced viability, a rough eye phenotype and severely malformed bristles. The shape of the adult bristle is determined by the actin bundles that are regularly spaced around the perimeter of the developing pupal bristles. Examination of the twf pupal bristles revealed an increased level of filamentous actin, which in turn resulted in splitting and displacement of the actin bundles. The bristle defect was rescued by twf overexpression in developing bristles. The Twinfilin protein was localised at sites of actin filament assembly, where it was required to limit actin polymerisation. A genetic interaction between twinfilin and twinstar (the gene encoding Cofilin) was detected, consistent with the model predicting that both proteins act to limit the amount of filamentous actin. α-Actinin has been implicated in several diverse cell biological processes. In Drosophila, the only function for α-actinin yet known is in the organisation of the muscle sarcomere. Muscle and non-muscle cells utilise different α-actinin isoforms, which in Drosophila are produced by alternative splicing of a single gene. In this work, novel α-actinin deletion alleles, including ActnΔ233, were generated, which specifically disrupted the transcript encoding the non-muscle α-actinin isoform. Nevertheless, ActnΔ233 homozygous mutant flies were viable and fertile with no obvious defects. By comparing α-actinin protein distribution in wild type and ActnΔ233 mutant animals, it could be concluded that non-muscle α-actinin is the only isoform expressed in young embryos, in the embryonic central nervous system and in various actin-rich structures of the ovarian germline cells. In the ActnΔ233 mutant, α-actinin was detected not only in muscle tissue, but also in embryonic epidermal cells and in certain follicle cell populations in the ovaries. The population of α-actinin protein present in non-muscle cells of the ActnΔ233 mutant is referred to as FC-α-actinin (Follicle Cell). The follicular epithelium in the Drosophila ovary is a well characterised model system for studies on patterning and morphogenesis. Therefore, α-actinin expression, regulation and function in this tissue were further analysed. Examination of the α-actinin localisation pattern revealed that the basal actin fibres of the main body follicle cells underwent an organised remodelling during the final stages of oogenesis. This involved the assembly of a transient adhesion site in the posterior of the cell, in which α-actinin and Enabled (Ena) accumulated. Follicle cells genetically manipulated to lack all α-actinin isoforms failed to remodel their cytoskeleton and translocate Ena to the posterior of the cell, while the actin fibres as such were not affected. Neither was epithelial morphogenesis disrupted. The reorganisation of the basal actin cytoskeleton was also disturbed following ectopic expression of Decapentaplegic (Dpp) or as a result of a heat shock. At late oogenesis, the main body follicle cells express both non-muscle α-actinin and FC-α-actinin, while the dorsal anterior follicle cells express only non-muscle α-actinin. The dorsal anterior cells are patterned by the Dpp and Epidermal growth factor receptor (EGFR) signalling pathways, and they will ultimately secrete the dorsal appendages of the egg. Experiments involving ectopic activation of EGFR and Dpp signalling showed that FC-α-actinin is negatively regulated by combined EGFR and Dpp signalling. Ubiquitous overexpression of the adult muscle-specific α-actinin isoform induced the formation of aberrant actin bundles in migrating follicle cells that did not normally express FC-α-actinin, provided that the EGFR signalling pathway was activated in the cells. Taken together, this work contributes new data to our knowledge of α-actinin function and regulation in Drosophila. The cytoskeletal remodelling shown to depend on α-actinin function provides the first evidence that α-actinin has a role in the organisation of the cytoskeleton in a non-muscle tissue. Furthermore, the cytoskeletal remodelling constitutes a previously undescribed morphogenetic event, which may provide us with a model system for in vivo studies on adhesion dynamics in Drosophila.

Positional cloning and pathway analysis of the asthma susceptibility gene, NPSR1

In the present study, we identified a novel asthma susceptibility gene, NPSR1 (neuropeptide S receptor 1) on chromosome 7p14.3 by the positional cloning strategy. An earlier significant linkage mapping result among Finnish Kainuu asthma families was confirmed in two independent cohorts: in asthma families from Quebec, Canada and in allergy families from North Karelia, Finland. The linkage region was narrowed down to a 133-kb segment by a hierarchial genotyping method. The observed 77-kb haplotype block showed 7 haplotypes and a similar risk and nonrisk pattern in all three populations studied. All seven haplotypes occur in all three populations at frequences > 2%. Significant elevated relative risks were detected for elevated total IgE (immunoglobulin E) or asthma. Risk effects of the gene variants varied from 1.4 to 2.5. NPSR1 belongs to the G protein-coupled receptor (GPCR) family with a topology of seven transmembrane domains. NPSR1 has 9 exons, with the two main transcripts, A and B, encoding proteins of 371 and 377 amino acids, respectively. We detected a low but ubiquitous expression level of NPSR1-B in various tissues and endogenous cell lines while NPSR1-A has a more restricted expression pattern. Both isoforms were expressed in the lung epithelium. We observed aberrant expression levels of NPSR1-B in smooth muscle in asthmatic bronchi as compared to healthy. In an experimental mouse model, the induced lung inflammation resulted in elevated Npsr1 levels. Furthermore, we demonstrated that the activation of NPSR1 with its endogenous agonist, neuropeptide S (NPS), resulted in a significant inhibition of the growth of NPSR1-A overexpressing stable cell lines (NPSR1-A cells). To determine which target genes were regulated by the NPS-NPSR1 pathway, NPSR1-A cells were stimulated with NPS, and differentially expressed genes were identified using the Affymetrix HGU133Plus2 GeneChip. A total of 104 genes were found significantly up-regulated and 42 down-regulated 6 h after NPS administration. The up-regulated genes included many neuronal genes and some putative susceptibility genes for respiratory disorders. By Gene Ontology enrichment analysis, the biological process terms, cell proliferation, morphogenesis and immune response were among the most altered. The expression of four up-regulated genes, matrix metallopeptidase 10 (MMP10), INHBA (activin A), interleukin 8 (IL8) and EPH receptor A2 (EPHA2), were verified and confirmed by quantitative reverse-transcriptase-PCR. In conclusion, we identified a novel asthma susceptibility gene, NPSR1, on chromosome 7p14.3. NPS-NPSR1 represents a novel pathway that regulates cell proliferation and immune responses, and thus may have functional relevance in the pathogenesis of asthma.