The role of pollen in the changing environmental conditions of Scots pine

Rapid change in climate is challenge for the adaptation of forest trees in the future. In wind pollinated tree species pollen mediated long distance gene flow may provide alleles that are (pre)adapted to a future climate. In order to examine the long distance pollen flow in Scots pine (Pinus sylvestris L.), we measured the amount and viability of airborne pollen and flowering phenology in central, northern, and northernmost Finland during four years. Viable airborne pollen grains were detected during female flowering and before local pollen shedding in all study sites. The situation when there was nonlocal pollen in the air lasted from one to four days depending on the year and study site. The amount of nonlocal airborne pollen varied also between years and study sites, the total amount of nonlocal viable pollen in the air was 2.3% from all detected viable pollen grains. The effect of pollen origin on seeds siring ability was studied with artificial pollination experiments. Pollen genotypes originating from southern Finland sired 76% and 48 % of the analysed seeds in competition studies where both pollen origin were introduced simultaneously into the female strobili. We examined the importance of arrival order of pollen grains in to the strobili in a study where pollen genotypes of different origin were introduced in two hours interval. Northern genotypes sired 76% of the analysed seeds when it was injected first, but in the "southern first" experiment both pollen types sired equal amount of seeds. The first pollen grain in the pollen chamber do not always fertilizes the ovum, instead there likely is more complex way of competition between pollen grains. To examine chemically mediated pollen-pollen interactions we conducted in vitro germination experiment where different pollen genotypes had chemical but not physical contact. Both positive and negative effects of interactions were found. We found highly negative effects in germinability of northern pollen grains when they were germinating with southern pollen, and increase in the germinability of southern pollen. There were no variation in the size of the dry pollen grains between pollen origins, and minor variation between different genotypes. After hydration and germination northern pollen grains were larger than southern pollen. Pollen genotypes having high hydration rates had low germinability and tube growth rate, however, germinated pollen grains were larger in size than nongerminated. This supports the suggestion that the early germination and growth of pollen tube is dependent on pollen storage materialsand less dependent on water intake and hydration. Long distance pollen movements and good competition ability of southern pollen makes gene flow possible, although rising temperature and timing of pollen movements may affect pollen competition and the amount of gene flow.

Mutation analysis of TGF-β signaling pathway genes among Finnish patients with primary pulmonary hypertension and hereditary hemorrhagic telangiectasia

Primary pulmonary hypertension (PPH), or according to the recent classification idiopathic pulmonary hypertension (IPAH), is a rare, progressive disease of pulmonary vasculature leading to pulmonary hypertension and right heart failure. Most of the patients are sporadic but in about 6% of cases the disease is familial (FPPH). In 2000 two different groups identified the gene predisposing to PPH. This gene, Bone morphogenetic protein receptor type 2 (BMPR2), encodes a subunit of transforming growth factor β (TGF-β) receptor complex. There is a genetic connection between PPH and hereditary hemorrhagic telangiectasia (HHT), a bleeding disorder characterized by local telangiectasias and sometimes with pulmonary hypertension. In HHT, mutations in ALK1 (activin like kinase type 1) and Endoglin, another members of the TGF-β signaling pathway are found. In this study we identified all of the Finnish PPH patients for the years 1986-1999 using the hospital discharge registries of Finnish university hospitals. During this period we found a total of 59 confirmed PPH patients: 55 sporadic and 4 familial representing 3 different families. In 1999 the prevalence of PPH was 5.8 per million and the annual incidence varied between 0.2-1.3 per million. Among 28 PPH patients studied, heterozygous BMPR2 mutations were found in 12% (3/26) of sporadic patients and in 33% of the PPH families (1/3). All the mutations found were different. Large deletions of BMPR2 were excluded by single-stranded chain polymomorphism analysis. As a candidate gene approach we also studied ALK1, Endoglin, Bone Morphogenetic Receptor Type IA (BMPR1A or ALK3), Mothers Against Decapentaplegic Homolog 4 (SMAD4) and Serotonine Transporter Gene (SLC6A4) using single-strand conformational polymorphism (SSCP) analysis and direct sequencing. Among patients and family members studied, we found two mutations in ALK1 in two unrelated samples. We also identified all the HHT patients treated at the Department of Otorhinolaryngology at Helsinki University Central Hospital between the years of 1990-2005 and 8 of the patients were studied for Endoglin and ALK1 mutations using direct sequencing. A total of seven mutations were found and all the mutations were different. The absence of a founder mutation in the Finnish population in both PPH and HHT was somewhat surprising. This suggests that the mutations of BMPR2, ALK1 and Endoglin are quite young and the older mutations have been lost due to repetitive genetic bottlenecks and/or negative selection. Also, other genes than BMPR2 may be involved in the pathogenesis of PPH. No founder mutations were found in PPH or HHT and thus no simple genetic test is available for diagnostics.

The DISC1 Pathway in the Genetic Etiology of Schizophrenia

Schizophrenia is a severe psychotic disorder affecting 0.5-1 % of the population. The disorder is characterized by hallucinations; delusions; disorganized behavior and speech; avolition; anhedonia; flattened affect and cognitive deficits. The etiology of the disorder is complex with evidence for multiple genes contributing to the onset of the disorder along with environmental factors. DISC1 is one of the most promising candidate genes for schizophrenia. It codes for a protein which takes part in numerous molecular interactions along several pathways. This network, termed as the DISC1 pathway, is evidently important for the development and maturation of the central nervous system from the embryo until young adulthood. Disruption at these pathways is thought to predispose schizophrenia. In the present study, we have studied the DISC1 pathway in the etiology of schizophrenia in the Finnish population. We have utilized large Finnish samples; the schizophrenia family sample where DISC1 was originally shown to associate with schizophrenia and the Northern Finland birth cohort 1966 (NFBC66). Several DISC1 binding partners displayed evidence for association in the family sample along with DISC1. Through a genome-wide linkage study, we found a significant linkage signal to a locus where a DISC1 binding partner NDE1 is located at the carriers of a certain DISC1 risk variant. In a follow-up study, genetic markers in NDE1 displayed significant evidence for association with schizophrenia. Further exploration of association between 11 genes of the DISC1 pathway and schizophrenia led to recognition of novel variants in NDEL1, PDE4B and PDE4D that significantly either increased or decreased the risk for schizophrenia. Further, we found evidence that DISC1 itself has a significant role in the human mental functioning even in the healthy population. Variants in DISC1 had a significant effect on anhedonia which is a trait present at everybody but is in its severe form one of the main symptoms of schizophrenia and correlates with the risk of developing the disorder. Further, utilizing genome-wide marker data, we recognized three genes; MIR620; CCDC141 and LCT; that are closely related to the DISC1 pathway but which effects on anhedonia were observable only at the individuals who carried these specific DISC1 variants. Our findings significantly add up to the previous evidence for the involvement of DISC1 and the DISC1 pathway in the etiology of schizophrenia and psychosis. Our results support the concept of a number of DISC1 pathway related genes contributing in the etiology of schizophrenia along with DISC1 and provide new candidates for the studies of schizophrenia. Our findings also significantly increase the importance of DISC1 itself as having a role in psychological functioning in the general population.

Autoxidation of conjugated linoleic acid methyl ester in the presence of α-tocopherol: the hydroperoxide pathway

The autoxidation of conjugated linoleic acid (CLA) is poorly understood in spite of increasing interest in the beneficial biological properties of CLA and growing consumption of CLA-rich foods. In this thesis, the autoxidation reactions of the two major CLA isomers, 9-cis,11-trans-octadecadienoic acid and 10-trans,12-cis-octadecadienoic acid, are investigated. The results contribute to an understanding of the early stages of the autoxidation of CLA methyl ester, and provide for the first time a means of producing and separating intact CLA methyl ester hydroperoxides as well as basic knowledge on lipid hydroperoxides and their hydroxy derivatives. Conjugated diene allylic monohydroperoxides were discovered as primary autoxidation products formed during autoxidation of CLA methyl esters in the presence and absence of α-tocopherol. This established that one of the autoxidation pathways of CLA methyl ester is the hydroperoxide pathway. Hydroperoxides were produced from the two major CLA methyl esters by taking advantage of the effect of α-tocopherol to promote hydroperoxide formation. The hydroperoxides were analysed and separated first as methyl hydroxyoctadecadienoates and then as intact hydroperoxides by HPLC. The isolated products were characterized by UV, GC-MS, and NMR techniques. In the presence of a high amount of α-tocopherol, the autoxidation of CLA methyl ester yields six kinetically-controlled conjugated diene monohydroperoxides and is diastereoselective in favour of one particular geometric isomer as a pair of enantiomers. The primary autoxidation products produced from the two major CLA isomers include new positional isomers of conjugated diene monohydroperoxides, the 8-, 10-, 12-, and 14-hydroperoxyoctadecadienoates. Furthermore, two of these new positional isomers have an unusual structure for a cis,trans lipid hydroperoxide where the allylic methine carbon is adjacent to the cis instead of the usual trans double bond. The 1H and 13C NMR spectra of nine isomeric methyl hydroxyoctadecadienoates and of ten isomeric methyl hydroperoxyoctadecadienoates including the unusual cis,trans hydroperoxides, i.e. Me 8-OOH-9c,11t and Me 14-OOH-10t,12c, were fully assigned with the aid of 2D NMR spectroscopy. The assigned NMR data enabled determination of the effects of the hydroxyl and hydroperoxyl groups on the carbon chemical shifts of CLA isomers, identification of diagnostic signals, and determination of chemical shift differences of the olefinic resonances that may help with the assignment of structure to as yet unknown lipid hydroperoxides either as hydroxy derivatives or as intact hydroperoxides. A mechanism for the hydroperoxide pathway of CLA autoxidation in the presence of a high amount of α-tocopherol was proposed based on the characterized primary products, their relative distribution, and theoretical calculations. This is an important step forward in CLA research, where exact mechanisms for the autoxidation of CLA have not been presented before. Knowledge of these hydroperoxide formation steps is of crucial importance for understanding the subsequent steps and the different pathways of the autoxidation of CLA. Moreover, a deeper understanding of the autoxidation mechanisms is required for ensuring the safety of CLA-rich foods. Knowledge of CLA oxidation and how it differs from the oxidation of nonconjugated polyunsaturated fatty acids may also be the key to understanding the biological mechanisms of CLA activity.

Lymphatic vessels in health and disease: Role of the VEGF-C/VEGFR-3 pathway and the transcription factor FOXC2

The circulatory system consists of two vessel types, which act in concert but significantly differ from each other in several structural and functional aspects as well as in mechanisms governing their development. The blood vasculature transports oxygen, nutrients and cells to tissues whereas the lymphatic vessels collect extravasated fluid, macromolecules and cells of the immune system and return them back to the blood circulation. Understanding the molecular mechanisms behind the developmental and functional regulation of the lymphatic system long lagged behind that of the blood vasculature. Identification of several markers specific for the lymphatic endothelium, and the discovery of key factors controlling the development and function of the lymphatic vessels have greatly facilitated research in lymphatic biology over the past few years. Recognition of the crucial importance of lymphatic vessels in certain pathological conditions, most importantly in tumor metastasis, lymphedema and inflammation, has increased interest in this vessel type, for so long overshadowed by its blood vascular cousin. VEGF-C (Vascular Endothelial Growth Factor C) and its receptor VEGFR-3 are essential for the development and maintenance of embryonic lymphatic vasculature. Furthermore, VEGF-C has been shown to be upregulated in many tumors and its expression found to positively correlate with lymphatic metastasis. Mutations in the transcription factor FOXC2 result in lymphedema-distichiasis (LD), which suggests a role for FOXC2 in the regulation of lymphatic development or function. This study was undertaken to obtain more information about the role of the VEGF-C/VEGFR-3 pathway and FOXC2 in regulating lymphatic development, growth, function and survival in physiological as well as in pathological conditions. We found that the silk-like carboxyterminal propeptide is not necessary for the lymphangiogenic activity of VEGF-C, but enhances it, and that the aminoterminal propeptide mediates binding of VEGF-C to the neuropilin-2 coreceptor, which we suggest to be involved in VEGF-C signalling via VEGFR-3. Furthermore, we found that overexpression of VEGF-C increases tumor lymphangiogenesis and intralymphatic tumor growth, both of which could be inhibited by a soluble form of VEGFR-3. These results suggest that blocking VEGFR-3 signalling could be used for prevention of lymphatic tumor metastasis. This might prove to be a safe treatment method for human cancer patients, since inhibition of VEGFR-3 activity had no effect on the normal lymphatic vasculature in adult mice, though it did lead to regression of lymphatic vessels in the postnatal period. Interestingly, in contrast to VEGF-C, which induces lymphangiogenesis already during embryonic development, we found that the related VEGF-D promotes lymphatic vessel growth only after birth. These results suggest, that the lymphatic vasculature undergoes postnatal maturation, which renders it independent of ligand induced VEGFR-3 signalling for survival but responsive to VEGF-D for growth. Finally, we show that FOXC2 is necessary for the later stages of lymphatic development by regulating the morphogenesis of lymphatic valves, as well as interactions of the lymphatic endothelium with vascular mural cells, in which it cooperates with VEGFR-3. Furthermore, our study indicates that the absence of lymphatic valves, abnormal association of lymphatic capillaries with mural cells and an increased amount of basement membrane underlie the pathogenesis of LD. These findings have given new insight into the mechanisms of normal lymphatic development, as well as into the pathogenesis of diseases involving the lymphatic vasculature. They also reveal new therapeutic targets for the prevention and treatment of tumor metastasis and lymphatic vascular failure in certain forms of lymphedema. Several interesting questions were posed that still need to be addressed. Most importantly, the mechanism of VEGF-C promoted tumor metastasis and the molecular nature of the postnatal lymphatic vessel maturation remain to be elucidated.

Engineering the pentose phosphate pathway of Saccharomyces cerevisiae for production of ethanol and xylitol

The baker s yeast Saccharomyces cerevisiae has a long tradition in alcohol production from D-glucose of e.g. starch. However, without genetic modifications it is unable to utilise the 5-carbon sugars D-xylose and L arabinose present in plant biomass. In this study, one key metabolic step of the catabolic D-xylose pathway in recombinant D-xylose-utilising S. cerevisiae strains was studied. This step, carried out by xylulokinase (XK), was shown to be rate-limiting, because overexpression of the xylulokinase-encoding gene XKS1 increased both the specific ethanol production rate and the yield from D xylose. In addition, less of the unwanted side product xylitol was produced. Recombinant D-xylose-utilizing S. cerevisiae strains have been constructed by expressing the genes coding for the first two enzymes of the pathway, D-xylose reductase (XR) and xylitol dehydrogenase (XDH) from the D-xylose-utilising yeast Pichia stipitis. In this study, the ability of endogenous genes of S. cerevisiae to enable D-xylose utilisation was evaluated. Overexpression of the GRE3 gene coding for an unspecific aldose reductase and the ScXYL2 gene coding for a xylitol dehydrogenase homologue enabled growth on D-xylose in aerobic conditions. However, the strain with GRE3 and ScXYL2 had a lower growth rate and accumulated more xylitol compared to the strain with the corresponding enzymes from P. stipitis. Use of the strictly NADPH-dependent Gre3p instead of the P. stipitis XR able to utilise both NADH and NADPH leads to a more severe redox imbalance. In a S. cerevisiae strain not engineered for D-xylose utilisation the presence of D-xylose increased xylitol dehydrogenase activity and the expression of the genes SOR1 or SOR2 coding for sorbitol dehydrogenase. Thus, D-xylose utilisation by S. cerevisiae with activities encoded by ScXYL2 or possibly SOR1 or SOR2, and GRE3 is feasible, but requires efficient redox balance engineering. Compared to D-xylose, D-glucose is a cheap and readily available substrate and thus an attractive alternative for xylitol manufacture. In this study, the pentose phosphate pathway (PPP) of S. cerevisiae was engineered for production of xylitol from D-glucose. Xylitol was formed from D-xylulose 5-phosphate in strains lacking transketolase activity and expressing the gene coding for XDH from P. stipitis. In addition to xylitol, ribitol, D-ribose and D-ribulose were also formed. Deletion of the xylulokinase-encoding gene increased xylitol production, whereas the expression of DOG1 coding for sugar phosphate phosphatase increased ribitol, D-ribose and D-ribulose production. Strains lacking phosphoglucose isomerase (Pgi1p) activity were shown to produce 5 carbon compounds through PPP when DOG1 was overexpressed. Expression of genes encoding glyceraldehyde 3-phosphate dehydrogenase of Bacillus subtilis, GapB, or NAD-dependent glutamate dehydrogenase Gdh2p of S. cerevisiae, altered the cellular redox balance and enhanced growth of pgi1 strains on D glucose, but co-expression with DOG1 reduced growth on higher D-glucose concentrations. Strains lacking both transketolase and phosphoglucose isomerase activities tolerated only low D-glucose concentrations, but the yield of 5-carbon sugars and sugar alcohols on D-glucose was about 50% (w/w).

Identification of the fungal catabolic D-galacturonate pathway

Pectin is a natural polymer consisting mainly of D-galacturonic acid monomers. Microorganisms living on decaying plant material can use D-galacturonic acid for growth. Although bacterial pathways for D-galacturonate catabolism had been described previously, no eukaryotic pathway for D-galacturonate catabolism was known at the beginning of this work. The aim of this work was to identify such a pathway. In this thesis the pathway for D-galacturonate catabolism was identified in the filamentous fungus Trichoderma reesei. The pathway consisted of four enzymes: NADPH-dependent D-galacturonate reductase (GAR1), L-galactonate dehydratase (LGD1), L-threo-3-deoxy-hexulosonate aldolase (LGA1) and NADPH-dependent glyceraldehyde reductase (GLD1). In this pathway D-galacturonate was converted to pyruvate and glycerol via L-galactonate, L-threo-3-deoxy-hexulosonate and L-glyceraldehyde. The enzyme activities of GAR1, LGD1 and LGA1 were present in crude mycelial extract only when T. reesei was grown on D-galacturonate. The activity of GLD1 was equally present on all the tested carbon sources. The corresponding genes were identified either by purifying and sequencing the enzyme or by expressing genes with homology to other similar enzymes in a heterologous host and testing the activities. The new genes that were identified were expressed in Saccharomyces cerevisiae and resulted in active enzymes. The GAR1, LGA1 and GLD1 were also produced in S. cerevisiae as active enzymes with a polyhistidine-tag, and purified and characterised. GAR1 and LGA1 catalysed reversible reactions, whereas only the forward reactions were observed for LGD1 and GLD1. When gar1, lgd1 or lga1 was deleted in T. reesei the deletion strain was unable to grow with D-galacturonate as the only carbon source, demonstrating that all the corresponding enzymes were essential for D-galacturonate catabolism and that no alternative D-galacturonate pathway exists in T. reesei. A challenge for biotechnology is to convert cheap raw materials to useful and more valuable products. Filamentous fungi are especially useful for the conversion of pectin, since they are efficient producers of pectinases. Identification of the fungal D-galacturonate pathway is of fundamental importance for the utilisation of pectin and its conversion to useful products.

FGFR1 regulated gene-expression, cell proliferation and differentiation in the developing midbrain and hindbrain

The neuroectodermal tissue close to the midbrain hindbrain boundary (MHB) is an important secondary organizer in the developing neural tube. This so-called isthmic organizer (IsO) regulates cellular survival, patterning and proliferation in the midbrain (Mb) and rhombomere 1 (R1) of the hindbrain. Signaling molecules of the IsO, such as fibroblast growth factor 8 (FGF8) and WNT1 are expressed in distinct bands of cells around the MHB. It has been previously shown that FGF-receptor 1 (FGFR1) is required for the normal development of this brain region in the mouse embryo. In the present study, we have compared the gene expression profiles of wild-type and Fgfr1 mutant embryos. We show that the loss of Fgfr1 results in the downregulation of several genes expressed close to the MHB and in the disappearance of gene expression gradients in the midbrain and R1. Our microarray screen identified several previously uncharacterized genes which may participate in the development of midbrain R1 region. Our results also show altered neurogenesis in the midbrain and R1 of the Fgfr1 mutants. Interestingly, the neuronal progenitors in midbrain and R1 show different responses to the loss of signaling through FGFR1. As Wnt1 expression at the MHB region requires the FGF signaling pathway, WNT target genes, including Drapc1, were also identified in our screen. The microarray data analysis also suggested that the cells next to the midbrain hindbrain boundary express distinct cell cycle regulators. We showed that the cells close to the border appeared to have unique features. These cells proliferate less rapidly than the surrounding cells. Unlike the cells further away from the boundary, these cells express Fgfr1 but not the other FGF receptors. The slowly proliferating boundary cells are necessary for development of the characteristic isthmic constriction. They may also contribute to compartmentalization of this brain region.

Multifunctional RNA silencing pathway polymerase QDE-1 of Neurospora crassa

Double-stranded RNA and associated proteins are known to regulate the gene expression of most eukaryotic organisms. These regulation pathways have different components, outcomes and distinct nomenclature depending on the model system, and often they are referred to collectively as RNA silencing. In many cases, RNA-dependent RNA polymerases (RdRPs) are found to be involved in the RNA silencing, but their targets, activities, interaction partners and reaction products remain enigmatic. In the filamentous fungus Neurospora crassa, the RdRP QDE-1 is critical for silencing of transgenes a phenomenon known as quelling. In this thesis the structure, biochemical activities and biological functions of QDE-1 were extensively studied. This dimeric RdRP was shown to possess five distinct catalytic in vitro activities that could be dissected by mutagenesis and by altering reaction conditions. The biochemical characterization implied that QDE-1 is actually an active DNA-dependent RNA polymerase that has additional RdRP activity. It also provided a structural explanation for the dimerization and suggested a biological framework for the functions of QDE-1 in vivo. (I) QDE-1 was also studied in a broader context along with the other components of the quelling pathway. It was shown that DNA damage in Neurospora causes a dramatic increase in the expression level of the Argonaute protein QDE-2 as well as the synthesis of a novel class of small RNAs known as qiRNAs. The accumulation of qiRNAs was shown to be dependent on several quelling components, and particularly to be derived from an aberrant ssRNA (aRNA) molecule that is synthesized by QDE-1 in the nucleus. The genomic distribution of qiRNA targets was analyzed and the possible biological significance of qiRNAs was studied. Importantly, qiRNAs are the first class of small RNAs that are induced by DNA damage. (II) After establishing that QDE-1 is a multifunctional RNA polymerase with several activities, template specificities and subcellular locations, the focus was turned onto its interaction partners. It had been previously known that QDE-1 associates with Replication Protein A (RPA), but the RecQ helicase QDE-3 was now shown to regulate this interaction. RPA was also observed to promote QDE-1 dependent dsRNA synthesis in vitro. By characterizing the interplay between QDE-1, QDE-3 and RPA, a working model of quelling and qiRNA pathways in Neurospora was presented. (III) This work sheds light on the complexity of the various RNA silencing pathways of a fungal model system. It shows how an RdRP can regulate gene expression on many levels, and suggests novel lines of research in other eukaryotic organisms.

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.