<title>Introduction</title> <p>Asthma involves genetic and environmental factors in its development, chronicity and severity ##REF##11818486##[1]##, ##REF##11005790##[2]##. Although some of its underlying mechanisms have been elucidated in recent years, more work is needed to gain a clearer understanding of genetic determinants. The mapping of asthma has been one of the most important areas of human genetics in the last two decades. According to an overview by Blumenthal (2005), twelve complete and two incomplete genome scans for asthma have been published, identifying a total of twenty chromosomal linked regions to asthma ##REF##15764904##[3]##. Discrepancies often appeared between linkage studies ##REF##15764904##[3]##–##REF##12413520##[5]##, leading to the use of standard phenotype definitions and founder populations as a way to decrease phenotypic and genetic heterogeneity ##REF##12083968##[6]##–##REF##11326283##[8]##.</p> <p>To date, common strategies have been employed to identify genes involved in asthma predisposition. Linkage studies followed by positional cloning identified six genes while association studies identified over a hundred genes, the majority of these having quite small effects on asthma susceptibility (see ##REF##14748924##[9]##, ##REF##16395390##[10]## for a review). Here, we describe a novel methodological approach that combines classical genetic approaches with a computational data-mining tool. In this regard, a genome-wide scan for asthma and atopy in families originating from the Saguenay–Lac-St-Jean (SLSJ) founder population (Northeastern Quebec, Canada) ##REF##7717408##[11]##–##REF##8091226##[16]## has been combined with the Genes to Diseases (G2D) new computational tool in the prioritization of asthma candidate genes ##REF##16115313##[17]##, ##UREF##0##[18]##. G2D performs the selection of the candidates in the chromosomal regions genetically linked to a disease by highlighting genes whose functions are related to the phenotype of the disease according to a data mining analysis of the bibliographic and protein public databases, or according to the genes already known to be associated with the same or a similar phenotype.</p> <p>Single nucleotide polymorphisms (SNPs) within the candidate genes prioritized by the G2D tool have been subsequently used to conduct an association study with asthma, atopy and allergic asthma phenotypes in the SLSJ asthma familial collection. The chromosomal regions around associated SNPs have then been sequenced in order to find causal mutations, and a replication of the positive association findings has been assessed in the Childhood Asthma Management Program (CAMP) independent cohort. The goal was to apply the G2D tool in the search of genetic determinants associated with a complex trait, using asthma as a model. This approach led to the hypothesis-driven identification of ten genes that may not have been selected otherwise (<italic>LPA</italic>, <italic>NOX3</italic>, <italic>SNX9</italic>, <italic>VIL2</italic>, <italic>VIP</italic>, <italic>ADAM8</italic>, <italic>DOCK1</italic>, <italic>FANK1</italic>, <italic>GPR123</italic> and <italic>PTPRE</italic>). Of these, one positive genetic association, resisting to corrections for multiple testing, has been found between the protein tyrosine phosphatase receptor type E gene (<italic>PTPRE</italic>) and allergic asthma in the SLSJ sample.</p>
<title>Methods</title> <title>Subjects</title> <p>Clinical evaluation and phenotyping criteria of the Saguenay-Lac-St-Jean (SLSJ) subjects have been described in recent reports ##REF##15282199##[19]##–##REF##17363771##[21]## and summarized in ##TAB##0##Table 1##. This familial sample is predominantly composed of probands that reported an onset age of asthma below 12 years old (81.6% of the probands). The mean age of onset for the probands is 7 years and the mean age of onset for the asthmatic family members is 22 years. The entire sample has been used for the association study while for the genome scan, the first 79 recruited families, that were available at the time when the genome scan was performed have been used (n = 609 individuals–see supplementary ##SUPPL##1##Table S1## for subjects characteristics and studied phenotypes). The Chicoutimi Hospital local ethics committee approved the study and all subjects provided informed consent.</p> <title>Genome Scan</title> <p>DNA was extracted for all SLSJ participants from whole blood by using the QIAGEN genomic purification procedure (QIAGEN Inc., Valencia, CA). Genotyping was completed on 367 autosomal and 21 X-chromosome microsatellite markers evenly spaced throughout the genome (average spacing of 9.2 cM). The marker set is a modification of the Cooperative Human Linkage Centre Screening Set (<ext-link ext-link-type="uri" xlink:href="http://gai.nci.nih.gov/CHLC/">http://gai.nci.nih.gov/CHLC/</ext-link>, version 6.0), showing an average heterozygosity of 0.72 in our data set. Each primer was amplified separately and then pooled into panels of eight markers and products were interrogated using ABI 3700 sequencers (Applied Biosystems, Foster City, CA) with a size standard ladder. Duplicate of two CEPH control DNAs and one water were included in each genotyping plate. In our data set of 195939 autosomal genotypes, 91.4% of the alleles were called and the proportion of observed Mendelian error was 1.2%. Linkage of chromosomal regions with putative genetic risk factors for a given trait was assessed by evaluating the extent of excess sharing of alleles identical by descent in affected relatives within families. Test statistics are reported on the LOD scale. Briefly, a multipoint, one-parameter likelihood ratio test that is robust against incompleteness of marker data (when the descent of alleles in a pedigree is not fully known) was used ##REF##9345087##[22]##, ##REF##8651312##[23]##. We moreover evaluated the specific contribution of mothers and fathers to the test of linkage, to look for parent-of-origin effects. The above tests of linkage, linkage through mothers and linkage through fathers are described in details in ##REF##16451620##[24]##.</p> <title>Genes to Diseases (G2D)</title> <p>The G2D tool has been applied in the two best genome scan susceptibility regions: 6q26 between markers D6S476 and D6S305, and 10q26.3 between markers D10S1223 and D10S1248. We used two of the approaches considered in G2D ##UREF##0##[18]## to pre-select gene lists in these regions. The first approach uses a description of the phenotype to point to genes in a region. This method works with automatically derived relationships between the disease symptoms (as MeSH C terms) and gene features (as Gene Ontology or GO terms ##REF##10802651##[25]##) that are obtained from the literature and Entrez gene database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/">http://www.ncbi.nlm.nih.gov/</ext-link>). We call this procedure the “PHENOTYPE” method. The second approach consists on automatically finding genes in a region that are similar to other genes previously associated to asthma. To do this, G2D measures the semantic distance between the annotations of the “known genes” and the annotations of genes in the problem region assigned by homology searches. We call this procedure the “KNOWN GENES” method. For details about the algorithm see the G2D web site at <ext-link ext-link-type="uri" xlink:href="http://www.ogic.ca/projects/g2d_2/">http://www.ogic.ca/projects/g2d_2/</ext-link> and ##REF##16115313##[17]##, ##UREF##0##[18]##.</p> <title>Association Study</title> <title>SNP selection</title> <p>The HapMap database (<ext-link ext-link-type="uri" xlink:href="http://www.hapmap.org/">http://www.hapmap.org/</ext-link>) has been used to identify SNPs assumed to be polymorphic in the SLSJ population. TagSNPs were then selected with the tagger program implemented in the Haploview software (version 3.32)##REF##15297300##[26]## using an r<sup>2</sup> cutoff of 0.8 and a minor allele frequency (MAF) over 0.10 to cover each whole gene. SNPs were also prioritized on their localization (coding or untranslated regions-see supplementary ##SUPPL##2##Table S2##). All SNPs are referred using their reference sequence number (rs#).</p> <title>SNP genotyping</title> <p>Eighty SNPs have been genotyped by the Sequenom® matrix-assisted laser desorption/ionization time-of-flight mass array spectrometer (Sequenom Inc., San Diego, CA) (##SUPPL##2##Table S2##). Sequenom primers were designed using the Sequenom SNP Assay Design software version 3.0 for iPLEX reactions. A total of 74 assays were designed for a single multiplex reaction. The assay group file containing the PCR primers and the iPLEX extension probes can be supplied on request to the corresponding author. The protocol and reaction conditions are in accordance with the manufacturer ##UREF##1##[27]##. The genotypes were viewed and analyzed using the MassARRAY Typer software version 3.4 (Sequenom Inc., San Diego, CA). The ten remaining SNPs (##SUPPL##2##Table S2##) have been genotyped by the TaqMan® SNP Genotyping Assays (Applied Biosystems, Foster City, CA) using the Rotor-Gene™ real-time PCR (Corbett Research Ltd, Sydney, Australia). Protocol and method were supplied by the manufacturer and PCR conditions were optimized to get a good cluster separation between different genotypes (see supplementary ##SUPPL##3##Table S3## for PCR conditions). Genotypes were attributed by the Rotor Gene software using the scatter graph analysis option.</p> <title>Statistical analysis</title> <p>Family-based association testing has been performed with the FBAT software (version 1.7) using an empirical estimate of the variance ##REF##11058432##[28]##–##REF##11055368##[30]## to correctly account for linkage. A Sidak correction for multiple testing has been applied on the p-values accounting for the effective number of independent phenotypes and SNPs, according to the definition of Li and Ji ##REF##16077740##[31]##, as implemented in the SNPSpD program ##REF##14997420##[32]## (see the online supporting ##SUPPL##0##Text S1## file for supplementary details). Parent-specific transmission disequilibrium tests were performed using sib_tdt from the ASPEX package (<ext-link ext-link-type="uri" xlink:href="http://aspex.sourceforge.net/">http://aspex.sourceforge.net/</ext-link>). Mendelian errors have been assessed by FBAT and Hardy-Weinberg equilibrium has been assessed with Haploview software (version 3.32)##REF##15297300##[26]##.</p> <title>
<italic>PTPRE</italic> Sequencing</title> <p>Forty unrelated SLSJ subjects (validated with the BALSAC database ##REF##12282730##[33]##) presenting full-fit allergic asthmatic criteria and that have contributed to the <italic>PTPRE</italic> association were selected. <italic>PTPRE</italic> sequence information was obtained from Ensembl database (<ext-link ext-link-type="uri" xlink:href="http://www.ensembl.org/index.html">http://www.ensembl.org/index.html</ext-link>, release 46). The sequencing was divided in six regions that spanned 2.7 kb, starting from the exon 2 to the exon 4. Oligonucleotides and PCR conditions are listed in ##SUPPL##4##Table S4##. Amplification products were purified with multiscreen PCR plates (Millipore Corporation, Billerica, MA), sequenced with BigDye terminator v3.1 chemistry following instructions of the manufacturer and analyzed on a 3100 Genetic analyzer (Applied Biosystems, Foster City, CA). Sequence analysis was performed with Codoncode Aligner software (CodonCode Corporation, <ext-link ext-link-type="uri" xlink:href="http://www.codoncode.com/">http://www.codoncode.com/</ext-link>). The identified SNPs presenting MAF over 0.05 have been genotyped in the SLSJ sample using a Sequenom panel and analyzed with FBAT, as described above. Newly described SNPs have been submitted to NCBI (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/">http://www.ncbi.nlm.nih.gov/</ext-link>) SNP database.</p> <title>Replication Study</title> <p>The <italic>PTPRE</italic> rs7081735 association has been assessed in the Childhood Asthma Management Program (CAMP) study ##REF##10027502##[34]##, ##REF##11027739##[35]##. This analysis includes the 497 non-Hispanic white children and their parents for whom adequate DNA was available. The genotyping of the <italic>PTPRE</italic> rs7081735 was performed using a TaqMan® SNP Genotyping Assays (Applied Biosystems, Foster City, CA) (see ##SUPPL##3##Table S3## for PCR conditions). Plates were scanned using the 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City, CA) and genotypes were assigned by the SDS 2.2 software using the scatter graph analysis option. The Institutional Review Board of the Brigham and Women's Hospital (BWH), as well as those of the other CAMP study centers, approved this study. Informed assent and consent were obtained from the study participants and their parents to collect DNA for genetic studies.</p>
<title>Results</title> <title>Genome Scan</title> <p>The genome-wide linkage scan analysis revealed at least two regions showing suggestive evidence for linkage, as well as differential maternal and paternal contribution (##FIG##0##Figure 1##). The two regions are 6q26 for asthma (LOD = 1.54, p = 0.0038) and 10q26.3 for atopy (LOD = 2.82, p = 0.00016). In these two regions, affected sibs tend to share more alleles inherited from the mothers than from the fathers. The linkage tests through the mothers reach a LOD of 2.19 (p = 0.00074) in 6q26 for asthma and a LOD of 2.96 (p = 0.00011) in 10q26.3 for atopy. The respective LODs obtained through the fathers at the same loci are only 0.01 (p = 0.40) and 0.58 (p = 0.05). Moreover, even though the 10q region does not show a great strength of linkage with asthma (LOD = 0.57, p = 0.051), the asthmatic sibs tend to share these alleles when received from their mothers (LOD = 2.82, p = 0.00016). Each of the two chromosomal regions show LOD values (either for tests of linkage, or parent-specific LODs, or both) that are, in order of magnitude, consistent with what has been defined as «suggestive» for linkage by Lander and Kruglyak ##REF##7581446##[36]##, which are expected to occur once per whole genome scan on average. Because of their great hypothesis generating potential, the 6q26 and 10q26 regions as well as asthma and atopy phenotypes were selected for the following G2D and association studies.</p> <title>G2D</title> <p>We applied the two algorithms PHENOTYPE and KNOWN GENES to the 6q26 and 10q26.3 regions. For the PHENOTYPE method, we used the OMIM ##REF##17170002##[37]## record 600807 as input, particularly the MeSH C terms from the MEDLINE references in OMIM entry 600807 which deals with susceptibility to asthma and asthma related traits. The most prevailing MeSH C terms are “Asthma” and “Bronchial Hyperreactivity”, but the list also includes “Hypersensitivity”, “Respiratory Hypersensitivity” and “Eosinophilia”. We thus considered that the terms that refer to asthma were “Asthma” and “Bronchial Hyperreactivity” and that the terms that refer to atopy were “Hypersensitivity”, “Respiratory Hypersensitivity” and “Eosinophilia”. The highest scoring GO terms, associated to these MeSH C terms, describe a variety of molecular functions and processes that include leukotriene and interleukin signaling, glutathione metabolism, etc (see header of supplementary ##SUPPL##5##Table S5##). We applied this method to the 6q26 region between D6S476 and D6S305, the two markers directly flanking the linkage peak seen in the region, which corresponds to the 10.48 MB band between positions 151,685,574 and 162,165,587 of chromosome 6. After discarding candidates that did not overlap with any known or hypothetical Entrez Gene sequence, 16 genes remained (see supplementary ##SUPPL##5##Table S5-A##). A similar analysis was carried between D10S1223 and D10S1248, the two markers directly flanking the linkage peak in 10q26.3, between positions 129,150,822 and 130,982,363 in chromosome 10. In that case, 12 candidates were obtained (see supplementary ##SUPPL##5##Table S5-B##). For the KNOWN GENES method, we compiled a list of genes reported to be associated with asthma and atopy from the literature ##REF##15143320##[38]## and from the Genetic Association Database GAD ##REF##15118671##[39]## (supplementary ##SUPPL##6##Table S6##). We then extracted the GO annotation of those genes in Entrez Gene ##REF##17170002##[37]##. We derived a scoring system for the candidates according to the minimal semantic distance between their GO annotation and the ones from the compiled known-gene list. For example, genes annotated with GO terms such as “dipeptidyl-peptidase IV activity” or “chemokine receptor binding” would score high as candidates. The method takes into account the hierarchical structure of GO as well as the specificity of GO terms. In that sense, similarity with more infrequent terms receive higher scores. For the complete list of GO terms see the header of supplementary ##SUPPL##7##Table S7##. We applied the KNOWN GENES method to both regions 6q26 and 10q26.3, obtaining 15 and 10 genes, respectively, after filtering out those candidates that did not overlap with either known or hypothetical genes (see supplementary ##SUPPL##7##Table S7##).</p> <title>Use of a complementary method based on genomic sequence</title> <p>We applied the G2D complementary Disease Gene Prediction (DGP) tool ##REF##15181176##[40]## to both genetically linked regions in order to predict the involvement of genes in inherited disease by their sequence features. This method analyses the probability of a gene to be associated with any disease phenotype. Genes with no associated phenotype and a probability greater than 0.7 were retained. With this criterion the DGP method identified one candidate in chromosome 10 (MMP21) and three candidates in the 6q26 region: <italic>TFB1M</italic> [MIM:607033], <italic>RGS17</italic> [MIM:607191] and <italic>VIP</italic> [MIM:192320].</p> <title>Final list of candidate genes</title> <p>Genes that received the higher scores in the pre-selected lists, and those that were pointed by the two different analyses were preferred, allowing the construction of a list of 17 candidates (displayed in the ##TAB##1##Table 2##). We then applied a candidate gene approach to select a final list of ten genes with the best biological potential related to asthma pathophysiology for genotyping (five in each chromosomal region): <italic>LPA</italic>, <italic>NOX3</italic>, <italic>SNX9</italic>, <italic>VIL2</italic>, <italic>VIP</italic>, <italic>ADAM8</italic>, <italic>DOCK1</italic>, <italic>FANK1</italic>, <italic>GPR123</italic> and <italic>PTPRE</italic> (marked in bold in ##TAB##1##Table 2##).</p> <title>Association Study</title> <p>A final panel of 91 SNPs (75 tagSNPs and 16 non-tagSNPs) was selected among the ten candidate genes (supplementary ##SUPPL##2##Table S2##). Of these, four were non-polymorphic, six failed Sequenom genotyping assays, none presented deviation from Hardy-Weinberg equilibrium (p-values &gt;0.001) and one presented more than two Mendelian errors (##SUPPL##2##Table S2##). For the remaining 80 SNPs, genotypes of individuals with Mendelian errors were considered as missing data in the FBAT analyses. Genotyping presented a mean success rate of 99.0% for the Sequenom assays and a mean success rate of 98.3% for the TaqMan assays. Accounting for the residual correlation between the tagSNPs, it is estimated that the 80 partially correlated SNPs correspond to an effective number of 59 independent ones ##REF##16077740##[31]##. As for the effective number of phenotypes, simulation indicates that the three studied phenotypes (asthma, atopy and allergic asthma) correspond to 1.75 effective independent ones (see supplementary ##SUPPL##0##Text S1## file). Accordingly, the estimated total effective number of independent tests is 103.25 (59 effective independent SNPs×1.75 effective independent phenotypes). Thus, applying Sidak correction, the p-value threshold of significance is estimated to be 0.000483.</p> <p>The FBAT single marker analyses were performed under an additive genetic model for each SNP and the three studied phenotypes. For the sake of brevity, only results showing a p-value under 0.05 before correction for multiple testing for one or more phenotypes are presented in ##TAB##2##Table 3##. For the asthma phenotype, minor alleles LPA_rs12175867C, GPR123_rs11101913T, GPR123_rs11101932T and GPR123_rs12257731A were overtransmitted to the asthmatic probands, suggesting a susceptibility effect of these alleles for asthma (0.0085&lt;p&lt;0.047). Inversely, minor alleles DOCK1_rs1051039G, PTPRE_rs4369314A and PTPRE_rs7081735G were undertransmitted to the asthmatic probands, suggesting a protective effect of these alleles for asthma (0.010&lt;p&lt;0.049). For the atopy phenotype, only ADAM8_rs11101672G and PTPRE_rs7081735G minor alleles have been undertransmitted to the atopic probands, suggesting a protective effect (p = 0.039 and 0.037, respectively). Finally, for the allergic asthma phenotype, minor alleles LPA_rs12175867C, GPR123_rs11101913T, GPR123_rs11101932T and GPR123_rs12257731A were overtransmitted to the allergic asthmatic probands, suggesting a susceptibility effect of these alleles for allergic asthma (0.035&lt;p&lt;0.041). Inversely, minor alleles ADAM8_rs11101672G, GPR123_rs11101916A, GPR123_rs761777G, PTPRE_rs11016002A, PTPRE_rs4002572C and PTPRE_rs7081735G were undertransmitted to the allergic asthmatic probands, suggesting a protective effect of these alleles for allergic asthma (0.000463&lt;p&lt;0.037). None of the SNPs reported above showed a significantly greater extent of transmission distortion from mothers than from fathers, thus not providing insights to the observed parental distortions seen in the linkage results (not shown).</p> <title>
<italic>PTPRE</italic> Sequencing</title> <p>
<italic>PTPRE</italic> rs7081735 shows the strongest association to allergic asthma (p = 0.000463) and is the only SNP shown to be significant after multiple testing correction (corrected p = 0.0478). We thus sequenced strategic genomic regions around the rs7081735, including coding regions in order to identify causal mutation. ##FIG##1##Figure 2## shows the <italic>PTPRE</italic> sequenced regions, the rs7081735 localization and the nine identified variants (diamonds), including four novel ones (c.86172A&gt;G, c.140901G&gt;A, c.140903G&gt;A and c.141102C&gt;T) (##TAB##3##Table 4##). It is worth to note that any of these identified polymorphisms may affect either <italic>PTPRE</italic> isoforms. A family-based association analysis has been performed between the three studied phenotypes and variants presenting a MAF over 0.05 that were not included in the first genotyping panel, which was the case for four of the nine identified variants (rs7911506, c.86172A&gt;G, rs7895103 and c.140901G&gt;A). Only rs7895103 showed a modest association with asthma, at a level that does not provide additional insights (p = 0.013, compared to p = 0.000463 for rs7081735).</p> <title>Replication Study</title> <p>Considering that the <italic>PTPRE</italic> rs7081735 association to allergic asthma is the strongest in the SLSJ sample, we evaluated it in an independent familial cohort, the Childhood Asthma Management Program (CAMP) ##REF##10027502##[34]##. The genotyping completion rate was 96% and no discordance was observed upon repeat genotyping of two random plates. The minor allele frequency was 0.33 and was in Hardy-Weinberg equilibrium (p = 0.90). Family-based association showed no significant association for asthma or atopy phenotypes (data not shown). However, to ensure that the selection of the CAMP study for the replication is appropriate and to demonstrate that the association found could result from a childhood subset of asthma instead of an adulthood one, we stratified the SLSJ association analyses considering only the probands that reported an asthma age of onset below 12 years old. Thus, the positive association found for <italic>PTPRE</italic> rs7081735 and allergic asthma (p = 0.000463, ##TAB##2##Table 3##) remains significant when considering only those probands (Family number = 74; Z for the minor allele = −3.166; p = 0.001546). Even with the loss of statistical power due to the stratified analysis, this comparison allowed to assume that the association found for <italic>PTPRE</italic> is probably more related to a childhood asthma, comforting the choice of the CAMP cohort as a replication study.</p>
<title>Discussion</title> <p>This study proposes a novel approach for the selection of candidate genes for asthma association studies using the computational G2D tool to find genetic determinants for this disease. Based on a genome-wide scan performed in an asthmatic familial sample from the SLSJ founder population, we selected the two best-linked regions (6q26 and 10q26.3) and applied the G2D data mining approach ##REF##16115313##[17]##, ##UREF##0##[18]## to identify ten candidate genes for an association in the SLSJ sample. Among these, five (<italic>LPA</italic>, <italic>ADAM8</italic>, <italic>DOCK1</italic>, <italic>GPR123</italic> and <italic>PTPRE</italic>) presented modest associations with asthma, atopy or allergic asthma. After corrections for multiple testing, only the <italic>PTPRE</italic> rs7081735 association to allergic asthma remained significant. These findings demonstrate that the G2D tool can be useful in the selection of candidate genes for asthma genetic studies.</p> <p>Because the <italic>PTPRE</italic> association to allergic asthma remained significant after correction for multiple testing, we sequenced strategic regions around the rs7081735, aiming to find the causal mutation. Sequencing allowed the identification of five known and four novel variants. Testing four SNPs with MAF &gt;0.05 did not identify additional associated SNPs with asthma or atopy related phenotypes in the SLSJ sample. Thus, the rs7081735, located in the 5′ untranslated region, could be the causal variant, or be in linkage disequilibrium with an unknown causal mutation.</p> <p>
<italic>PTPRE</italic> is a member of the protein tyrosine phosphatase (PTPs) family, which includes genes that are important regulators of signal transduction pathways involved in various cellular processes such as control of metabolic pathways, cellular adhesion, cell cycle progression and immune response ##REF##15186772##[41]##, ##REF##10723800##[42]##. <italic>PTPRE</italic> encodes two different isoforms, cytoplasmic and transmembrane ##REF##8579581##[43]##, based on its different promoters ##REF##9914474##[44]##. PTPs receptors participate in transmembrane signaling and cellular adhesion processes, whereas intracellular PTPs take part in signal transduction within the cell ##REF##9795134##[45]##. For example in mice, PTPε-deficient macrophages present abnormalities in the regulation of the respiratory burst and the production of cytokines in response to bacterial lipopolysaccharide, suggesting a role of the PTPε isoform in inflammation as well as in host defense ##REF##11485326##[46]##. However, <italic>PTPRE</italic> has been shown to be highly expressed in peripheral human monocytes and granulocytes, and antigen-receptor stimulation induces the expression of <italic>PTPRE</italic> in activated lymphocytes ##REF##12121439##[47]##. These observations and our finding suggest a protective effect of <italic>PTPRE</italic> in allergic asthma, and lead us to hypothesize that this potential protective role could involve the leukocyte cellular processes in the limitation of lung inflammation following an allergen sensitization. Further work is needed to define the <italic>PTPRE</italic> possible role in asthma pathophysiology.</p> <p>The <italic>PTPRE</italic> rs7081735 association has been evaluated in the independent CAMP cohort. Results showed no positive association for asthma and atopic phenotypes. Taking into account that the SLSJ and the CAMP studies are both family-based designed, well powered ##REF##15322985##[48]##, ##REF##17702965##[49]## and presenting a childhood onset asthma, the <italic>PTPRE</italic> association lack of replication in the CAMP study can be explained by other reasons: natural variability of asthma history ##REF##16950268##[50]##, ##REF##17627712##[51]##, differences of proband mean age between samples (SLSJ mean age of 18, and CAMP mean age of 8, respectively), differences in the genetic background of the two populations ##REF##11882781##[52]## (SLSJ individuals descend predominantly from French European founders ##REF##7717408##[11]##–##REF##8091226##[16]## whereas CAMP individuals come from a white North-American admixed population), or population specific gene-gene and gene-environment interactions ##REF##11882781##[52]##. Based on these observations, we conclude that the <italic>PTPRE</italic> rs7081735 association to allergic asthma is more penetrant in the SLSJ population, possibly resulting from an interaction with others genes and/or environmental factors that are more common in this founder population.</p> <p>In summary, this study demonstrates that the G2D tool can be useful in the prioritization of candidate genes for a complex disease as it allowed us to find a novel asthma genetic association with the <italic>PTPRE</italic> gene in the SLSJ familial asthma sample. This association represents a potential protective factor for asthma pathogenesis, as it is more likely related to a childhood onset asthma. The present genetic study is an example of how the combination of different methodological approaches can be relevant to target asthma genetic determinants and to motivate further genetic and functional investigations.</p>
<p>Conceived and designed the experiments: CPI MAAN CL. Performed the experiments: CP AT CPI CL. Analyzed the data: KT ML GMH CPI CL. Contributed reagents/materials/analysis tools: KT ML CP. Wrote the paper: KT.</p> <p>Asthma is a complex trait for which different strategies have been used to identify its environmental and genetic predisposing factors. Here, we describe a novel methodological approach to select candidate genes for asthma genetic association studies. In this regard, the Genes to Diseases (G2D) computational tool has been used in combination with a genome-wide scan performed in a sub-sample of the Saguenay−Lac-St-Jean (SLSJ) asthmatic familial collection (n = 609) to identify candidate genes located in two suggestive loci shown to be linked with asthma (6q26) and atopy (10q26.3), and presenting differential parent-of-origin effects. This approach combined gene selection based on the G2D data mining analysis of the bibliographic and protein public databases, or according to the genes already known to be associated with the same or a similar phenotype. Ten genes (<italic>LPA</italic>, <italic>NOX3</italic>, <italic>SNX9</italic>, <italic>VIL2</italic>, <italic>VIP</italic>, <italic>ADAM8</italic>, <italic>DOCK1</italic>, <italic>FANK1</italic>, <italic>GPR123</italic> and <italic>PTPRE</italic>) were selected for a subsequent association study performed in a large SLSJ sample (n = 1167) of individuals tested for asthma and atopy related phenotypes. Single nucleotide polymorphisms (n = 91) within the candidate genes were genotyped and analysed using a family-based association test. The results suggest a protective association to allergic asthma for <italic>PTPRE</italic> rs7081735 in the SLSJ sample (p = 0.000463; corrected p = 0.0478). This association has not been replicated in the Childhood Asthma Management Program (CAMP) cohort. Sequencing of the regions around rs7081735 revealed additional polymorphisms, but additional genotyping did not yield new associations. These results demonstrate that the G2D tool can be useful in the selection of candidate genes located in chromosomal regions linked to a complex trait.</p>
<title>Supporting Information</title>
<p>We thank families of the Saguenay–Lac-St-Jean for their participation in this study. We also thank Janet Murphy for her participation in the genome scan as well as Alexandre Belisle, Pierre Lepage and Charleen Salesse for the genotyping assays. K. Tremblay is an AllerGen PhD trainee and is supported by the Fondation de l'Université Laval studentship. T.J. Hudson received an Investigator Award from the Canadian Institutes of Health Research (CIHR) and a Clinician-scientist Award in Translational Research from the Burroughs Wellcome Fund. M.A. Andrade-Navarro holds a Canada Research Chair in Bioinformatics. C. Laprise is the chairholder of the Canada Research Chair (<ext-link ext-link-type="uri" xlink:href="http://www.chairs.gc.ca">http://www.chairs.gc.ca</ext-link>) on genetic determinants in asthma and the director of the Genetics platform of the Respiratory Health Network of the Fonds de la recherche en santé du Québec. We also thank all subjects of the CAMP study for their ongoing participation. We acknowledge the CAMP investigators and research team, supported by NHLBI, for collection of CAMP Genetic Ancillary Study data. All work on data collected from the CAMP Genetic Ancillary Study was conducted at the Channing Laboratory of the Brigham and Women's Hospital under appropriate CAMP policies and human subject's protections.</p>
<fig id="pone-0002907-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002907.g001</object-id><label>Figure 1</label><caption><title>Genome scan results summary.</title><p>Results from the tests of linkage with atopy (top three panels) and asthma (bottom three panels) reported on the LOD scale. For each phenotype, the top panel shows the results of the test of linkage (excess allele sharing), the middle panel shows results from the tests of linkage through mothers (excess allele sharing transmitted from mothers) and the bottom panel shows results from the tests of linkage through fathers (excess allele sharing transmitted from fathers).</p></caption></fig> <fig id="pone-0002907-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002907.g002</object-id><label>Figure 2</label><caption><title>
<italic>PTPRE</italic> gene sequenced regions and identified SNPs scaled location.</title><p>The black thick boxes above the gene define its sequenced parts, which are identified by the same numbers used in the ##TAB##3##Table 4##, in which exact chromosomal positions are available. Studied SNPs are represented below the gene. The TagSNPs correspond to an asterisk (*) and the SNPs identified by sequencing correspond to a diamond (◊). All <italic>PTPRE</italic> numbers for the discovered SNPs are based on the mRNA sequence NM_006504 (variant 1, receptor form) and from the NCBI (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/">http://www.ncbi.nlm.nih.gov/</ext-link>) SNP database (build 127). Image source: HapMap (<ext-link ext-link-type="uri" xlink:href="http://www.hapmap.org/">www.hapmap.org/</ext-link>) October 2007 (Genome Browser, <ext-link ext-link-type="uri" xlink:href="http://genome.ucsc.edu/">http://genome.ucsc.edu/</ext-link>, version 1.69); modified according to our study design.</p></caption></fig>
<table-wrap id="pone-0002907-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002907.t001</object-id><label>Table 1</label><caption><title>Clinical and phenotypic characteristics of the Saguenay–Lac-St-Jean association study sample subjects.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Probands (n = 226)</td><td align="left" rowspan="1" colspan="1">Family members (n = 941)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Male: Female ratio</td><td align="left" rowspan="1" colspan="1">1: 1.2</td><td align="left" rowspan="1" colspan="1">1: 1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">Mean age in years (range)</td><td align="left" rowspan="1" colspan="1">18 (3–46)</td><td align="left" rowspan="1" colspan="1">44 (2–96)</td></tr><tr><td align="left" rowspan="1" colspan="1">Smoking status (n (%))</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Never</td><td align="left" rowspan="1" colspan="1">185 (83.3)</td><td align="left" rowspan="1" colspan="1">424 (45.8)</td></tr><tr><td align="left" rowspan="1" colspan="1">Ex-smoker</td><td align="left" rowspan="1" colspan="1">12 (5.4)</td><td align="left" rowspan="1" colspan="1">307 (33.2)</td></tr><tr><td align="left" rowspan="1" colspan="1">Smoker</td><td align="left" rowspan="1" colspan="1">25 (11.3)</td><td align="left" rowspan="1" colspan="1">195 (21.0)</td></tr><tr><td align="left" rowspan="1" colspan="1">FEV<sub>1</sub> % of predicted value (SD) <xref ref-type="table-fn" rid="nt101">*</xref>
</td><td align="left" rowspan="1" colspan="1">92.3 (16.3)</td><td align="left" rowspan="1" colspan="1">94.4 (20.5)</td></tr><tr><td align="left" rowspan="1" colspan="1">PC<sub>20</sub> in mg/ml (SD) <xref ref-type="table-fn" rid="nt102">†</xref>
</td><td align="left" rowspan="1" colspan="1">2.68 (3.59)</td><td align="left" rowspan="1" colspan="1">10.87 (5.23)</td></tr><tr><td align="left" rowspan="1" colspan="1">Serum IgE in µg/l (SD) <xref ref-type="table-fn" rid="nt103">‡</xref>
</td><td align="left" rowspan="1" colspan="1">233.2 (4.6)</td><td align="left" rowspan="1" colspan="1">108.0 (4.2)</td></tr><tr><td colspan="3" align="left" rowspan="1">Number of persons with subphenotypes (%)</td></tr><tr><td align="left" rowspan="1" colspan="1">Asthma <xref ref-type="table-fn" rid="nt104">§</xref>
</td><td align="left" rowspan="1" colspan="1">226/226 (100)</td><td align="left" rowspan="1" colspan="1">353/935 (37.8)</td></tr><tr><td align="left" rowspan="1" colspan="1">Atopy <xref ref-type="table-fn" rid="nt105">∏</xref>
</td><td align="left" rowspan="1" colspan="1">185/224 (82.6)</td><td align="left" rowspan="1" colspan="1">468/920 (50.9)</td></tr><tr><td align="left" rowspan="1" colspan="1">Allergic asthma (Asthma+Atopy)</td><td align="left" rowspan="1" colspan="1">185/224 (82.6)</td><td align="left" rowspan="1" colspan="1">239/344 (69.5)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002907-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002907.t002</object-id><label>Table 2</label><caption><title>Genes identified by G2D data mining analysis and those selected for the association study based on their number of appearance in G2D analyses and on their biological function.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gene <xref ref-type="table-fn" rid="nt107">*</xref>
</td><td align="left" rowspan="1" colspan="1">GeneID/[MIM]</td><td align="left" rowspan="1" colspan="1">Analysis <xref ref-type="table-fn" rid="nt108">†</xref>
</td><td align="left" rowspan="1" colspan="1">Gene definition <xref ref-type="table-fn" rid="nt109">‡</xref>
</td><td align="left" rowspan="1" colspan="1">Gene function <xref ref-type="table-fn" rid="nt109">‡</xref>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>ARID1B</italic> 6q26</td><td align="left" rowspan="1" colspan="1">57492/[MIM:NA]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">AT rich interactive domain 1B (SWI1-like)</td><td align="left" rowspan="1" colspan="1">Involved in transcriptional activation and repression of select genes by chromatin remodeling (alteration of DNA-nucleosome topology). Binds DNA non-specifically.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ESR1</italic> 6q26</td><td align="left" rowspan="1" colspan="1">2099/[MIM:133430]</td><td align="left" rowspan="1" colspan="1">Yes</td><td align="left" rowspan="1" colspan="1">Estrogen receptor 1</td><td align="left" rowspan="1" colspan="1">Nuclear hormone receptor. The steroid hormones and their receptors are involved in the regulation of eukaryotic gene expression and affect cellular proliferation and differentiation in target tissues.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>IGF2R</italic> 6q26</td><td align="left" rowspan="1" colspan="1">3482/[MIM:147280]</td><td align="left" rowspan="1" colspan="1">Yes</td><td align="left" rowspan="1" colspan="1">Insulin-like growth factor 2 receptor</td><td align="left" rowspan="1" colspan="1">Transmembrane protein with a short cytoplasmic tail containing an internalization signal. This receptor binds insulin-like growth factor IGF2.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>LPA</italic></bold> 6q26</td><td align="left" rowspan="1" colspan="1">4018/[MIM:152200]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">Lipoprotein, Lp(a)</td><td align="left" rowspan="1" colspan="1">Apo(a) is the main constituent of lipoprotein(a). It has a serine proteinase activity and is able of autoproteolysis. Inhibits tissue-type plasminogen activator 1.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>NOX3</italic></bold> 6q26</td><td align="left" rowspan="1" colspan="1">50508 [MIM:607105]</td><td align="left" rowspan="1" colspan="1">Yes</td><td align="left" rowspan="1" colspan="1">NADPH oxidase 3</td><td align="left" rowspan="1" colspan="1">NADPH oxidases, such as NOX3, are plasma membrane-associated enzymes found in many cell types. They catalyze the production of superoxide by a 1-electron reduction of oxygen, using NADPH as the electron donor.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>SLC22A2</italic> (OCT2) 6q26</td><td align="left" rowspan="1" colspan="1">6582/[MIM:602608]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">Solute carrier family 22 (organic cation transporter), member 2</td><td align="left" rowspan="1" colspan="1">Polyspecific transporter of organic cations, mainly expressed in kidney (luminal membrane of distal tube).</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>SNX9</italic></bold> 6q26</td><td align="left" rowspan="1" colspan="1">51429/[MIM:605952]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">Sorting nexin 9</td><td align="left" rowspan="1" colspan="1">This gene encodes a member of the sorting nexin family. Members of this family contain a phox (PX) domain, which is a phosphoinositide binding domain, and are involved in several stages of intracellular trafficking.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>TFB1M</italic> 6q26</td><td align="left" rowspan="1" colspan="1">51106/[MIM:607033]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">Transcription factor B1, mitochondrial</td><td align="left" rowspan="1" colspan="1">The transcription of genes from mitochondrial DNA requires a mitochondrial RNA polymerase and a DNA-binding transcription factor. Transcription factor B1 (TFB1M) is a part of this transcription complex and is implied in the rRNA processing.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>VIL2</italic></bold> 6q26</td><td align="left" rowspan="1" colspan="1">7430/[MIM:123900]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">Villin 2 (ezrin)</td><td align="left" rowspan="1" colspan="1">The cytoplasmic peripheral membrane protein that functions as a protein-tyrosine kinase substrate in microvilli. Serves as an intermediate between the plasma membrane and the actin cytoskeleton. Plays a key role in cell surface structure adhesion, migration, and organization.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>VIP</italic></bold> 6q26</td><td align="left" rowspan="1" colspan="1">7432/[MIM:192320]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">Vasoactive intestinal peptid</td><td align="left" rowspan="1" colspan="1">Secreted protein belongs to the glucagon family. It stimulates myocardial contractility, causes vasodilatation, increases glycogenolysis, lowers arterial blood pressure and relaxes the smooth muscle of trachea, stomach and gall bladder.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>ADAM8</italic></bold> 10q26.3</td><td align="left" rowspan="1" colspan="1">101/[MIM:602267]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">ADAM metallopeptidase domain 8</td><td align="left" rowspan="1" colspan="1">Membrane-anchored proteins structurally related to snake venom disintegrins expressed in granulocyte, monocyte, and macrophage. Implicated in a variety of biological processes involving cell-cell and cell-matrix interactions. May be involved in cell adhesion during neurodegeneration.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>CPXM2</italic> 10q26.3</td><td align="left" rowspan="1" colspan="1">119587/[MIM:NA]</td><td align="left" rowspan="1" colspan="1">Yes</td><td align="left" rowspan="1" colspan="1">Carboxypeptidase X, member 2</td><td align="left" rowspan="1" colspan="1">Implicated in cell adhesion and proteolysis.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>DOCK1</italic></bold> 10q26.3</td><td align="left" rowspan="1" colspan="1">1793/[MIM:601403]</td><td align="left" rowspan="1" colspan="1">Yes</td><td align="left" rowspan="1" colspan="1">Dedicator of cytokinesis 1</td><td align="left" rowspan="1" colspan="1">Role in signaling from focal adhesions. This gene product binds to the SH3 domain of CRK protein. It may regulate cell surface extension and may have a role in the cell surface extension of an engulfing cell around a dying cell during apoptosis.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>FANK1</italic></bold> 10q26.3</td><td align="left" rowspan="1" colspan="1">92565/[MIM:611640]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">Fibronectin type III and ankyrin repeat domains 1</td><td align="left" rowspan="1" colspan="1">NA</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>GPR123</italic></bold> 10q26.3</td><td align="left" rowspan="1" colspan="1">84435/[MIM:NA]</td><td align="left" rowspan="1" colspan="1">No</td><td align="left" rowspan="1" colspan="1">G protein-coupled receptor 123</td><td align="left" rowspan="1" colspan="1">Orphan receptor : multi-pass membrane protein.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>MKI67</italic> 10q26.3</td><td align="left" rowspan="1" colspan="1">4288/[MIM:176741]</td><td align="left" rowspan="1" colspan="1">Yes</td><td align="left" rowspan="1" colspan="1">Antigen identified by monoclonal antibody Ki-67</td><td align="left" rowspan="1" colspan="1">Predominantly localized in the G1 phase in the perinucleolar region and in the nuclear matrix. In mitosis, present on all chromosomes. Thought to be required for maintaining cell proliferation. Regulation of progression through cell cycle.</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>PTPRE</italic></bold> 10q26.3</td><td align="left" rowspan="1" colspan="1">5791/[MIM:600926]</td><td align="left" rowspan="1" colspan="1">Yes</td><td align="left" rowspan="1" colspan="1">Protein tyrosine phosphatase, receptor type E</td><td align="left" rowspan="1" colspan="1">PTPs are known to be signaling molecules that regulate a variety of cellular processes including cell growth, differentiation, mitotic cycle, and oncogenic transformation.</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002907-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002907.t003</object-id><label>Table 3</label><caption><title>Significant Family-Based Association Test (FBAT) results between the ten G2D candidate genes studied SNPs and asthma, atopy and allergic asthma phenotypes under an additive genetic model.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gene</td><td align="left" rowspan="1" colspan="1">SNP</td><td align="left" rowspan="1" colspan="1">Base change</td><td align="left" rowspan="1" colspan="1">MAF</td><td colspan="3" align="left" rowspan="1">Asthma</td><td colspan="3" align="left" rowspan="1">Atopy</td><td colspan="3" align="left" rowspan="1">Allergic Asthma</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">Z</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">Z</td><td align="left" rowspan="1" colspan="1">P</td><td align="left" rowspan="1" colspan="1">N</td><td align="left" rowspan="1" colspan="1">Z</td><td align="left" rowspan="1" colspan="1">P</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>LPA</italic>
</td><td align="left" rowspan="1" colspan="1">rs12175867</td><td align="left" rowspan="1" colspan="1">T&gt;C</td><td align="left" rowspan="1" colspan="1">0.16</td><td align="left" rowspan="1" colspan="1">81</td><td align="left" rowspan="1" colspan="1">1.98</td><td align="left" rowspan="1" colspan="1">0.0474</td><td align="left" rowspan="1" colspan="1">68</td><td align="left" rowspan="1" colspan="1">1.94</td><td align="left" rowspan="1" colspan="1">0.0530</td><td align="left" rowspan="1" colspan="1">71</td><td align="left" rowspan="1" colspan="1">2.04</td><td align="left" rowspan="1" colspan="1">0.0412</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ADAM8</italic>
</td><td align="left" rowspan="1" colspan="1">rs11101672</td><td align="left" rowspan="1" colspan="1">C&gt;G</td><td align="left" rowspan="1" colspan="1">0.23</td><td align="left" rowspan="1" colspan="1">95</td><td align="left" rowspan="1" colspan="1">−1.76</td><td align="left" rowspan="1" colspan="1">0.0784</td><td align="left" rowspan="1" colspan="1">79</td><td align="left" rowspan="1" colspan="1">−2.06</td><td align="left" rowspan="1" colspan="1">0.0390</td><td align="left" rowspan="1" colspan="1">82</td><td align="left" rowspan="1" colspan="1">−2.12</td><td align="left" rowspan="1" colspan="1">0.0340</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>DOCK1</italic>
</td><td align="left" rowspan="1" colspan="1">rs1051039</td><td align="left" rowspan="1" colspan="1">C&gt;G</td><td align="left" rowspan="1" colspan="1">0.41</td><td align="left" rowspan="1" colspan="1">117</td><td align="left" rowspan="1" colspan="1">−1.97</td><td align="left" rowspan="1" colspan="1">0.0486</td><td align="left" rowspan="1" colspan="1">101</td><td align="left" rowspan="1" colspan="1">−0.72</td><td align="left" rowspan="1" colspan="1">0.4739</td><td align="left" rowspan="1" colspan="1">103</td><td align="left" rowspan="1" colspan="1">−0.48</td><td align="left" rowspan="1" colspan="1">0.6336</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>GPR123</italic>
</td><td align="left" rowspan="1" colspan="1">rs11101913</td><td align="left" rowspan="1" colspan="1">C&gt;T</td><td align="left" rowspan="1" colspan="1">0.29</td><td align="left" rowspan="1" colspan="1">112</td><td align="left" rowspan="1" colspan="1">2.63</td><td align="left" rowspan="1" colspan="1">0.0085</td><td align="left" rowspan="1" colspan="1">94</td><td align="left" rowspan="1" colspan="1">1.81</td><td align="left" rowspan="1" colspan="1">0.0707</td><td align="left" rowspan="1" colspan="1">97</td><td align="left" rowspan="1" colspan="1">2.04</td><td align="left" rowspan="1" colspan="1">0.0412</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">rs11101916</td><td align="left" rowspan="1" colspan="1">G&gt;A</td><td align="left" rowspan="1" colspan="1">0.18</td><td align="left" rowspan="1" colspan="1">95</td><td align="left" rowspan="1" colspan="1">−1.79</td><td align="left" rowspan="1" colspan="1">0.0731</td><td align="left" rowspan="1" colspan="1">73</td><td align="left" rowspan="1" colspan="1">−1.31</td><td align="left" rowspan="1" colspan="1">0.1916</td><td align="left" rowspan="1" colspan="1">80</td><td align="left" rowspan="1" colspan="1">−2.16</td><td align="left" rowspan="1" colspan="1">0.0308</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">rs11101932</td><td align="left" rowspan="1" colspan="1">C&gt;T</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">80</td><td align="left" rowspan="1" colspan="1">2.26</td><td align="left" rowspan="1" colspan="1">0.0240</td><td align="left" rowspan="1" colspan="1">60</td><td align="left" rowspan="1" colspan="1">1.74</td><td align="left" rowspan="1" colspan="1">0.0816</td><td align="left" rowspan="1" colspan="1">66</td><td align="left" rowspan="1" colspan="1">2.26</td><td align="left" rowspan="1" colspan="1">0.0397</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">rs761777</td><td align="left" rowspan="1" colspan="1">A&gt;G</td><td align="left" rowspan="1" colspan="1">0.29</td><td align="left" rowspan="1" colspan="1">113</td><td align="left" rowspan="1" colspan="1">−1.30</td><td align="left" rowspan="1" colspan="1">0.1948</td><td align="left" rowspan="1" colspan="1">92</td><td align="left" rowspan="1" colspan="1">−1.00</td><td align="left" rowspan="1" colspan="1">0.3152</td><td align="left" rowspan="1" colspan="1">99</td><td align="left" rowspan="1" colspan="1">−2.31</td><td align="left" rowspan="1" colspan="1">0.0209</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">rs12257731</td><td align="left" rowspan="1" colspan="1">G&gt;A</td><td align="left" rowspan="1" colspan="1">0.12</td><td align="left" rowspan="1" colspan="1">77</td><td align="left" rowspan="1" colspan="1">2.15</td><td align="left" rowspan="1" colspan="1">0.0314</td><td align="left" rowspan="1" colspan="1">60</td><td align="left" rowspan="1" colspan="1">1.65</td><td align="left" rowspan="1" colspan="1">0.0995</td><td align="left" rowspan="1" colspan="1">65</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">0.0353</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>PTPRE</italic>
</td><td align="left" rowspan="1" colspan="1">rs11016002</td><td align="left" rowspan="1" colspan="1">A&gt;T</td><td align="left" rowspan="1" colspan="1">0.35</td><td align="left" rowspan="1" colspan="1">113</td><td align="left" rowspan="1" colspan="1">−1.04</td><td align="left" rowspan="1" colspan="1">0.2983</td><td align="left" rowspan="1" colspan="1">89</td><td align="left" rowspan="1" colspan="1">−1.33</td><td align="left" rowspan="1" colspan="1">0.1836</td><td align="left" rowspan="1" colspan="1">93</td><td align="left" rowspan="1" colspan="1">−2.47</td><td align="left" rowspan="1" colspan="1">0.0136</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">rs4369314</td><td align="left" rowspan="1" colspan="1">G&gt;A</td><td align="left" rowspan="1" colspan="1">0.21</td><td align="left" rowspan="1" colspan="1">96</td><td align="left" rowspan="1" colspan="1">−2.57</td><td align="left" rowspan="1" colspan="1">0.0101</td><td align="left" rowspan="1" colspan="1">81</td><td align="left" rowspan="1" colspan="1">−1.65</td><td align="left" rowspan="1" colspan="1">0.0999</td><td align="left" rowspan="1" colspan="1">83</td><td align="left" rowspan="1" colspan="1">−1.67</td><td align="left" rowspan="1" colspan="1">0.0958</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">rs4002572</td><td align="left" rowspan="1" colspan="1">C&gt;T</td><td align="left" rowspan="1" colspan="1">0.45</td><td align="left" rowspan="1" colspan="1">113</td><td align="left" rowspan="1" colspan="1">−1.61</td><td align="left" rowspan="1" colspan="1">0.1067</td><td align="left" rowspan="1" colspan="1">86</td><td align="left" rowspan="1" colspan="1">−1.03</td><td align="left" rowspan="1" colspan="1">0.3047</td><td align="left" rowspan="1" colspan="1">93</td><td align="left" rowspan="1" colspan="1">−2.09</td><td align="left" rowspan="1" colspan="1">0.0367</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">rs7081735</td><td align="left" rowspan="1" colspan="1">A&gt;G</td><td align="left" rowspan="1" colspan="1">0.32</td><td align="left" rowspan="1" colspan="1">111</td><td align="left" rowspan="1" colspan="1">−2.23</td><td align="left" rowspan="1" colspan="1">0.0256</td><td align="left" rowspan="1" colspan="1">88</td><td align="left" rowspan="1" colspan="1">−2.09</td><td align="left" rowspan="1" colspan="1">0.0371</td><td align="left" rowspan="1" colspan="1">91</td><td align="left" rowspan="1" colspan="1">−3.50</td><td align="left" rowspan="1" colspan="1">
<bold>0.000463</bold>
</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002907-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002907.t004</object-id><label>Table 4</label><caption><title>
<italic>PTPRE</italic> sequenced regions and characteristics of identified SNPs.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Sequenced regions [bp] <xref ref-type="table-fn" rid="nt113">*</xref>
</td><td align="left" rowspan="1" colspan="1">Identified SNP <xref ref-type="table-fn" rid="nt114">†</xref>
</td><td align="left" rowspan="1" colspan="1">Base change</td><td align="left" rowspan="1" colspan="1">MAF (SLSJ) <xref ref-type="table-fn" rid="nt115">‡</xref>
</td><td align="left" rowspan="1" colspan="1">MAF (CEU) <xref ref-type="table-fn" rid="nt116">§</xref>
</td><td align="left" rowspan="1" colspan="1">SNP location <xref ref-type="table-fn" rid="nt113">*</xref>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1 [75542–75977]</td><td colspan="5" align="left" rowspan="1">
<italic>No SNPs found</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">2 [85761–86205]</td><td align="left" rowspan="1" colspan="1">
<bold>rs7911506</bold>
</td><td align="left" rowspan="1" colspan="1">C&gt;T</td><td align="left" rowspan="1" colspan="1">0.49</td><td align="left" rowspan="1" colspan="1">0.296</td><td align="left" rowspan="1" colspan="1">Intron 2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>c.86172 A&gt;G</bold>
</td><td align="left" rowspan="1" colspan="1">A&gt;G</td><td align="left" rowspan="1" colspan="1">0.31</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">Intron 2</td></tr><tr><td align="left" rowspan="1" colspan="1">3 [88011–88457]</td><td align="left" rowspan="1" colspan="1">rs4369314</td><td align="left" rowspan="1" colspan="1">G&gt;A</td><td align="left" rowspan="1" colspan="1">0.17</td><td align="left" rowspan="1" colspan="1">0.192</td><td align="left" rowspan="1" colspan="1">Intron 2</td></tr><tr><td align="left" rowspan="1" colspan="1">4 [92004–92507]</td><td align="left" rowspan="1" colspan="1">rs7081735</td><td align="left" rowspan="1" colspan="1">A&gt;G</td><td align="left" rowspan="1" colspan="1">0.47</td><td align="left" rowspan="1" colspan="1">0.342</td><td align="left" rowspan="1" colspan="1">Intron 2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>rs7895103</bold>
</td><td align="left" rowspan="1" colspan="1">T&gt;C</td><td align="left" rowspan="1" colspan="1">0.05</td><td align="left" rowspan="1" colspan="1">0.217</td><td align="left" rowspan="1" colspan="1">Intron 2</td></tr><tr><td align="left" rowspan="1" colspan="1">5 [134156–134591]</td><td colspan="5" align="left" rowspan="1">
<italic>No SNPs found</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">6 [140684–141152]</td><td align="left" rowspan="1" colspan="1">rs4262642</td><td align="left" rowspan="1" colspan="1">C&gt;A</td><td align="left" rowspan="1" colspan="1">0.19</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">Intron 3</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>c.140901 G&gt;A</bold>
</td><td align="left" rowspan="1" colspan="1">G&gt;A</td><td align="left" rowspan="1" colspan="1">0.40</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">Intron 3</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">c.140903 G&gt;A</td><td align="left" rowspan="1" colspan="1">G&gt;A</td><td align="left" rowspan="1" colspan="1">0.03</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">Intron 3</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">c.141102 C&gt;T</td><td align="left" rowspan="1" colspan="1">C&gt;T</td><td align="left" rowspan="1" colspan="1">0.03</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">Intron 4</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0002907.s001"><label>Text S1</label><caption><p>(0.03 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002907.s002"><label>Table S1</label><caption><p>Genome-wide scan Saguenay-Lac-St-Jean subjects clinical characteristics and studied phenotypes</p><p>(0.06 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002907.s003"><label>Table S2</label><caption><p>Characteristics of the 91 selected SNPs</p><p>(0.25 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002907.s004"><label>Table S3</label><caption><p>Taqman PCR conditions</p><p>(0.06 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002907.s005"><label>Table S4</label><caption><p>Oligonucleotides used for PTPRE sequencing</p><p>(0.05 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002907.s006"><label>Table S5</label><caption><p>Genes selected by G2D « PHENOTYPE » analysis</p><p>(0.07 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002907.s007"><label>Table S6</label><caption><p>Genes known or suspected to be associated with asthma that were considered for the G2D «KNOWN GENES» analysis</p><p>(0.07 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002907.s008"><label>Table S7</label><caption><p>Genes selected by G2D « KNOWN GENES » analysis</p><p>(0.06 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>*</label><p>FEV<sub>1</sub> = Geometric mean of the force expiratory volume in one second evaluated for 208 probands and for 683 family members.</p></fn><fn id="nt102"><label>†</label><p>PC<sub>20</sub> = Concentration of methacholine inducing a 20% fall in FEV1. Geometric mean and SD were obtained from the log transformed PC20 values. Evaluated for 186 probands and for 614 family members.</p></fn><fn id="nt103"><label>‡</label><p>IgE = Immunoglobulin E serum concentration. Geometric mean and SD were obtained from the log transformed IgE values. Evaluated for 211 probands and 704 family members.</p></fn><fn id="nt104"><label>§</label><p>Present asthma or past documented clinical history of asthma. The reported mean age of onset is 7 years among the asthmatic probands and 22 years among the asthmatic family members.</p></fn><fn id="nt105"><label>∏</label><p>Defined as at least one positive response on skin prick testing (wheal diameter &gt;3 mm at 10 minutes).</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt106"><label/><p>Genes marked in bold indicate that they have been selected for the association study. NA = Not available.</p></fn><fn id="nt107"><label>*</label><p>Gene symbol, chromosomal location and GeneID number (Entrez Gene identifier in the NCBI database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/">http://www.ncbi.nlm.nih.gov/</ext-link>)).</p></fn><fn id="nt108"><label>†</label><p>Appears in at least two G2D analyses.</p></fn><fn id="nt109"><label>‡</label><p>From NCBI (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/">http://www.ncbi.nlm.nih.gov/</ext-link>) and GeneCards (<ext-link ext-link-type="uri" xlink:href="http://www.genecards.org/">http://www.genecards.org/</ext-link>) databases.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt110"><label/><p>Abbreviations used: SNP = Single nucleotide polymorphism, MAF = Minor allele frequency, N = Number of families contributing to the statistic, Z = Z score, P = p-value (significance threshold of 0.000483; p-values under are marked in bold).</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt111"><label/><p>SNPs marked in bold have been genotyped and analyzed.</p></fn><fn id="nt112"><label/><p>Abbreviations used: bp = base pairs, SNP = single nucleotide polymorphism, MAF = Minor allele frequency, SLSJ = Saguenay-Lac-St-Jean, CEU = CEPH Utah residents with ancestry from northern and western Europe, NA = Not available.</p></fn><fn id="nt113"><label>*</label><p>Base pair numbering according to NCBI (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/">http://www.ncbi.nlm.nih.gov/</ext-link>) genome assembly (Build 36.2).</p></fn><fn id="nt114"><label>†</label><p>If known from NCBI SNP and Ensembl (<ext-link ext-link-type="uri" xlink:href="http://www.ensembl.org/index.html">http://www.ensembl.org/index.html</ext-link>) databases, reported by the NCBI rs#.</p></fn><fn id="nt115"><label>‡</label><p>Calculated on 39 individuals originating from the SLSJ population.</p></fn><fn id="nt116"><label>§</label><p>Obtained from the HapMap international project database (<ext-link ext-link-type="uri" xlink:href="http://www.hapmap.org/">http://www.hapmap.org/</ext-link>).</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>The Canadian Institutes of Health Research and the Respiratory Health Network of the Fonds de la recherche en santé du Québec financially supported this work. The Childhood Asthma Management Program Genetics Ancillary Study is supported by U01 HL075419, U01 HL65899, P01 HL083069, R01 HL 086601, and T32 HL07427 from the National Heart, Lung and Blood Institute, National Institutes of Health.</p></fn></fn-group>
<graphic id="pone-0002907-t001-1" xlink:href="pone.0002907.t001"/> <graphic xlink:href="pone.0002907.g001"/> <graphic id="pone-0002907-t002-2" xlink:href="pone.0002907.t002"/> <graphic id="pone-0002907-t003-3" xlink:href="pone.0002907.t003"/> <graphic xlink:href="pone.0002907.g002"/> <graphic id="pone-0002907-t004-4" xlink:href="pone.0002907.t004"/>
<media xlink:href="pone.0002907.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002907.s002.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002907.s003.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002907.s004.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002907.s005.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002907.s006.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002907.s007.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002907.s008.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["18"], "element-citation": ["\n"], "surname": ["Perez-Iratxeta", "Bork", "Andrade-Navarro"], "given-names": ["C", "P", "MA"], "year": ["2007"], "article-title": ["Update of the G2D tool for prioritization of gene candidates to inherited diseases."], "source": ["Nucleic Acids Res Jul"], "volume": ["1"], "fpage": ["W212"], "lpage": ["W216"]}, {"label": ["27"], "element-citation": ["\n"], "surname": ["Oeth", "Beaulieu", "Park", "Kosman", "del Mistro"], "given-names": ["P", "M", "C", "D", "G"], "year": ["2005"], "article-title": ["iPLEX assay: Increased plexing efficiency and flexibility for MassArray system through single base primer extension with mass-modified terminators."], "source": ["SEQUENOM Application Note"]}]
PLoS One. 2008 Aug 6; 3(8):e2907