A Gly98Val mutation in the N-Myc downstream regulated gene 1 (NDRG1) in Alaskan Malamutes with polyneuropathy.

The first cases of early-onset progressive polyneuropathy appeared in the Alaskan Malamute population in Norway in the late 1970s. Affected dogs were of both sexes and were ambulatory paraparetic, progressing to non-ambulatory tetraparesis. On neurologic examination, affected dogs displayed predominantly laryngeal paresis, decreased postural reactions, decreased spinal reflexes and muscle atrophy. The disease was considered eradicated through breeding programmes but recently new cases have occurred in the Nordic countries and the USA. The N-myc downstream-regulated gene (NDRG1) is implicated in neuropathies with comparable symptoms or clinical signs both in humans and in Greyhound dogs. This gene was therefore considered a candidate gene for the polyneuropathy in Alaskan Malamutes. The coding sequence of the NDRG1 gene derived from one healthy and one affected Alaskan Malamute revealed a non-synonymous G>T mutation in exon 4 in the affected dog that causes a Gly98Val amino acid substitution. This substitution was categorized to be "probably damaging" to the protein function by PolyPhen2 (score: 1.000). Subsequently, 102 Alaskan Malamutes from the Nordic countries and the USA known to be either affected (n = 22), obligate carriers (n = 7) or healthy (n = 73) were genotyped for the SNP using TaqMan. All affected dogs had the T/T genotype, the obligate carriers had the G/T genotype and the healthy dogs had the G/G genotype except for 13 who had the G/T genotype. A protein alignment showed that residue 98 is conserved in mammals and also that the entire NDRG1 protein is highly conserved (94.7%) in mammals. We conclude that the G>T substitution is most likely the mutation that causes polyneuropathy in Alaskan Malamutes. Our characterization of a novel candidate causative mutation for polyneuropathy offers a new canine model that can provide further insight into pathobiology and therapy of human polyneuropathy. Furthermore, selection against this mutation can now be used to eliminate the disease in Alaskan Malamutes.

Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Inducible caspase 9 suicide gene to improve the safety of allodepleted T cells after haploidentical stem cell transplantation.

Addback of donor T cells following T cell-depleted stem cell transplantation (SCT) can accelerate immune reconstitution and be effective against relapsed malignancy. After haploidentical SCT, a high risk of graft-versus-host disease (GVHD) essentially precludes this option, unless the T cells are first depleted of alloreactive precursor cells. Even then, the risks of severe GVHD remain significant. To increase the safety of the approach and thereby permit administration of larger T cell doses, we used a suicide gene, inducible caspase 9 (iCasp9), to transduce allodepleted T cells, permitting their destruction should administration have adverse effects. We made a retroviral vector encoding iCasp9 and a selectable marker (truncated CD19). Even after allodepletion (using anti-CD25 immunotoxin), donor T cells could be efficiently transduced, expanded, and subsequently enriched by CD19 immunomagnetic selection to >90% purity. These engineered cells retained antiviral specificity and functionality, and contained a subset with regulatory phenotype and function. Activating iCasp9 with a small-molecule dimerizer rapidly produced >90% apoptosis. Although transgene expression was downregulated in quiescent T cells, iCasp9 remained an efficient suicide gene, as expression was rapidly upregulated in activated (alloreactive) T cells. We have demonstrated the clinical feasibility of this approach after haploidentical transplantation by scaling up production using clinical grade materials.

Gene transfer by viral vectors

To initiate our clinical trial for chemotherapy protection, I established the retroviral vector system for human MDR1 cDNA gene transfer. The human MDR1 cDNA continued to be expressed in the transduced bone marrow cells after four cohorts of serial transplants, 17 months after the initial transduction and transplant. In addition, we used this retroviral vector pVMDR1 to transduce human bone marrow and peripheral blood CD34$\sp+$ cells on stromal monolayer in the presence of hematopoietic growth factors. These data suggest that the retroviral vector pVMDR1 could modify hematopoietic precursor cells with a capacity for long-term self renewal. Thus, it may be possible to use the MDR1 retroviruses to confer chemotherapeutic protection on human normal hematopoietic precursor cells of ovarian and breast cancer patients in whom high doses of MDR drugs may be required to control the diseases.^ Another promising vector system is recombinant adeno-associated virus (rAAV) vector. An impediment to use rAAV vectors is that production of rAAV vectors for clinical use is extremely cumbersome and labor intensive. First I set up the rAAV vector system in our laboratory and then, I focused on studies related to the production of rAAV vectors for clinical use. By using a self-inactivating retroviral vector carrying a selection marker under the control of the CMV immediate early promoter and an AAV genome with the deletion of both ITRs, I have developed either a transient or a stable method to produce rAAV vectors. These methods involve infection only and can generate high-titer rAAV vectors (up to 2 x 10$\sp5$ cfu/ml of CVL) with much less work.^ Although recombinant adenoviral vectors hardly infect early hematopoietic precursor cells lacking $\alpha\sb v\beta\sb5$ or $\alpha\sb v\beta\sb3$ integrin on their surface, but efficiently infect other cells, we can use these properties of adenoviral vectors for bone marrow purging as well as for development of new viral vectors such as pseudotyped retroviral vectors and rAAV vectors. Replacement of self-inactivating retroviral vectors by recombinant adenoviral vectors will facilitate the above strategies for production of new viral vectors. In order to accomplish these goals, I developed a new method which is much more efficient than the current methods to construct adenoviral vectors. This method involves a cosmid vector system which is utilized to construct the full-length recombinant adenoviral vectors in vitro.^ First, I developed an efficient and flexible method for in vitro construction of the full-length recombinant adenoviral vectors in the cosmid vector system by use of a three-DNA fragment ligation. Then, this system was improved by use of a two-DNA fragment ligation. The cloning capacity of recombinant adenoviral vectors constructed by this method to develop recombinant adenoviral vectors depends on the efficiency of transfection only. No homologous recombination is required for development of infectious adenoviral vectors. Thus, the efficiency of generating the recombinant adenoviral vectors by the cosmid method reported here was much higher than that by the in vitro direct ligation method or the in vivo homologous recombination method reported before. This method of the in vitro construction of recombinant adenoviral vectors in the cosmid vector system may facilitate the development of adenoviral vector for human gene therapy. (Abstract shortened by UMI.) ^

Prevalence and prognostic significance of TMPRSS2-ERG gene fusion in lymph node positive prostate cancers.

BACKGROUND TMPRSS2-ERG gene fusion is the most frequent genetic alteration in prostate cancer. However, information about its distribution in lymph node positive prostate cancers and the prognostic significance in these advanced tumors is unknown. METHODS Gene fusion status was determined by fluorescence in situ hybridization on a tissue-microarray constructed from 119 hormone-naïve nodal positive, surgically treated prostate cancers containing samples from the primary tumors and corresponding lymph node metastases. Data were correlated with various tumor features (Gleason score, stage, cancer volume, nodal tumor burden) and biochemical recurrence-free, disease-specific, and overall survival. RESULTS TMPRSS2-ERG fusion was detected in 43.5% of the primary tumors. Conversely, only 29.9% of the metastasizing components showed the fusion. Concordance in TMPRSS2-ERG status between primary tumors and metastases was 70.9% (Kappa 0.39); 20.9% and 8.1% of the patients showed the mutation solely in their primary tumors and metastases, respectively. TMPRSS2-ERG fusion was not correlated with specific histopathological tumor features but predicted favorable biochemical recurrence-free, disease-specific and overall survival independently when present in the primary tumor (P < 0.05 each). CONCLUSION TMPRSS2-ERG fusion is more frequent in primary prostate cancer than in corresponding metastases suggesting no selection of fusion-positive cells in the metastatic process. The gene fusion in primary tumors independently predicts favorable outcome.

Mutations upstream of fabI in triclosan resistant Staphylococcus aureus strains are associated with elevated fabI gene expression

Background The enoyl-acyl carrier protein (ACP) reductase enzyme (FabI) is the target for a series of antimicrobial agents including novel compounds in clinical trial and the biocide triclosan. Mutations in fabI and heterodiploidy for fabI have been shown to confer resistance in S. aureus strains in a previous study. Here we further determined the fabI upstream sequence of a selection of these strains and the gene expression levels in strains with promoter region mutations. Results Mutations in the fabI promoter were found in 18% of triclosan resistant clinical isolates, regardless the previously identified molecular mechanism conferring resistance. Although not significant, a higher rate of promoter mutations were found in strains without previously described mechanisms of resistance. Some of the mutations identified in the clinical isolates were also detected in a series of laboratory mutants. Microarray analysis of selected laboratory mutants with fabI promoter region mutations, grown in the absence of triclosan, revealed increased fabI expression in three out of four tested strains. In two of these strains, only few genes other than fabI were upregulated. Consistently with these data, whole genome sequencing of in vitro selected mutants identified only few mutations except the upstream and coding regions of fabI, with the promoter mutation as the most probable cause of fabI overexpression. Importantly the gene expression profiling of clinical isolates containing similar mutations in the fabI promoter also showed, when compared to unrelated non-mutated isolates, a significant up-regulation of fabI. Conclusions In conclusion, we have demonstrated the presence of C34T, T109G, and A101C mutations in the fabI promoter region of strains with fabI up-regulation, both in clinical isolates and/or laboratory mutants. These data provide further observations linking mutations upstream fabI with up-regulated expression of the fabI gene.

Gene flow and the maintenance of species boundaries

Hybrid zones are regions where individuals from genetically differentiated populations meet and mate, resulting in at least some offspring of mixed ancestry. Patterns of gene flow (introgression) in hybrid zones vary across the genome, allowing assessment of the role of individual genes or genome regions in reproductive isolation. Here, we document patterns of introgression between two recently diverged species of field crickets. We sampled at a very fine spatial scale and genotyped crickets for 110 highly differentiated single nucleotide polymorphisms (SNPs) identified through transcriptome scans. Using both genomic and geographic cline analysis, we document remarkably abrupt transitions (<100 m) in allele frequencies for 50 loci, despite high levels of gene flow at other loci. These are among the steepest clines documented for any hybridizing taxa. Furthermore, the cricket hybrid zone provides one of the clearest examples of the semi-permeability of species boundaries. Comparisons between data from the fine-scale transect and data (for the same set of markers) from sampling a much larger area in a different region of the cricket hybrid zone reveal consistent patterns of introgression for individual loci. The consistency in patterns of introgression between these two distant and distinct regions of the hybrid zone suggests that strong selection is acting to maintain abrupt discontinuities within the hybrid zone and that genomic regions with restricted introgression likely include genes that contribute to nonecological prezygotic barriers.

Color polymorphism and sex ratio distortion in a cichlid fish as an incipient stage in sympatric speciation by sexual selection

We investigated a Lake Victoria cichlid with a complex colour polymorphism that apparently represents one original species and two incipient species, all of which are sympatric. In laboratory breeding experiments we observed sex ratio distortion in certain matings between original and incipient species. Mate choice experiments show that males of the incipient species exhibit mating preferences against the original species, and males and females of the original species exhibit strong mating preferences against the incipient species. Mating preferences might evolve by sex ratio selection to avoid matings with distorted progeny sex ratios. Phenotype frequencies in nature suggest that mating preferences translate into mating frequencies, thus restricting gene flow and exerting disruptive sexual selection between the original and incipient species. The incipient species do not differ in morphology or ecology from the original species, implying that colour polymorphism, associated with sex ratio distortion, can be an incipient stage in sympatric speciation, and that disruption of gene flow can precede ecological differentiation

The gene for histone RNA hairpin binding protein is located on human chromosome 4 and encodes a novel type of RNA binding protein.

The hairpin structure at the 3' end of animal histone mRNAs controls histone RNA 3' processing, nucleocytoplasmic transport, translation and stability of histone mRNA. Functionally overlapping, if not identical, proteins binding to the histone RNA hairpin have been identified in nuclear and polysomal extracts. Our own results indicated that these hairpin binding proteins (HBPs) bind their target RNA as monomers and that the resulting ribonucleoprotein complexes are extremely stable. These features prompted us to select for HBP-encoding human cDNAs by RNA-mediated three-hybrid selection in Saccharomyces cerevesiae. Whole cell extract from one selected clone contained a Gal4 fusion protein that interacted with histone hairpin RNA in a sequence- and structure-specific manner similar to a fraction enriched for bovine HBP, indicating that the cDNA encoded HBP. DNA sequence analysis revealed that the coding sequence did not contain any known RNA binding motifs. The HBP gene is composed of eight exons covering 19.5 kb on the short arm of chromosome 4. Translation of the HBP open reading frame in vitro produced a 43 kDa protein with RNA binding specificity identical to murine or bovine HBP. In addition, recombinant HBP expressed in S. cerevisiae was functional in histone pre-mRNA processing, confirming that we have indeed identified the human HBP gene.

Selection signatures in Shetland ponies.

Shetland ponies were selected for numerous traits including small stature, strength, hardiness and longevity. Despite the different selection criteria, Shetland ponies are well known for their small stature. We performed a selection signature analysis including genome-wide SNPs of 75 Shetland ponies and 76 large-sized horses. Based upon this dataset, we identified a selection signature on equine chromosome (ECA) 1 between 103.8 Mb and 108.5 Mb. A total of 33 annotated genes are located within this interval including the IGF1R gene at 104.2 Mb and the ADAMTS17 gene at 105.4 Mb. These two genes are well known to have a major impact on body height in numerous species including humans. Homozygosity mapping in the Shetland ponies identified a region with increased homozygosity between 107.4 Mb and 108.5 Mb. None of the annotated genes in this region have so far been associated with height. Thus, we cannot exclude the possibility that the identified selection signature on ECA1 is associated with some trait other than height, for which Shetland ponies were selected.

ARRANGEMENT OF GENE LOCI IN EUPLOID CHINESE HAMSTER CELLS AND GENETIC CONSEQUENCES OF REARRANGEMENT IN CHO CELLS

In both euploid Chinese hamster (Cricetulus griseus) cells and pseudodiploid Chinese hamster ovary (CHO) cells, gene assignments were accomplished by G band chromosome and isozyme analysis (32 isozymes) of interspecific somatic cell hybrids obtained after HAT selection of mouse CL 1D (TK('-)) cells which were PEG-fused with either euploid Chinese hamster cells or HPRT('-) CHO cells. Hybrids slowly segregated hamster chromosomes. Clone panels consisting of independent hybrid clones and subclones containing different combinations of Chinese hamster chromosomes and isozymes were established from each type of fusion.^ These clone panels enabled us to provisionally assign the loci for: nucleoside phosphorylase (NP), glyoxalase (GLO), glutathione reductase (GSR), adenosine kinase (ADK), esterase D (ESD), peptidases B and S (PEPB and -S) and phosphoglucomutase 2 (PGM2, human nomenclature) to chromosome 1; adenylate kinase 1 (AK1), adenosine deaminase (ADA) and inosine triosephosphatase (ITP) to chromosome 6; triosephosphate isomerase (TPI) to chromosome 8; and glucose phosphate isomerse (GPI) and peptidase D (PEPD) to chromosome 9.^ We also confirm the assignments of 6-phosphogluconate dehydrogenase (PGD), PGM1, enolase 1 (ENO1) and diptheria toxin sensitivity (DTS) to chromosome 2 as well as provisionally assign galactose-1-phosphate uridyl transferase (GALT) and AK2 to chromosome 2. Selection in either HAT or BrdU for hybrids that had retained or lost the chromosome carrying the locus for TK enabled us to assign the loci for TK, galactokinase (GALK) and acid phosphatase 1 (ACP1) to Chinese hamster chromosome 7.^ These results are discussed in relation to current theories on the basis for high frequency of drug resistant autosomal recessive mutants in CHO cells and conservation of mammalian autosomal linkage groups. ^

INCREASED EXPRESSION OF ONE OF TWO ADENOSINE DEAMINASE ALLELES IN A HUMAN CHORIOCARCINOMA CELL LINE FOLLOWING SELECTION WITH ADENINE NUCLEOSIDES

The human choriocarcinoma cell line JEG-3 is heterozygous at the adenosine deaminase (ADA) gene locus. Both allelic genes are under strong but incomplete repression causing a very low level expression of the gene locus. Because cytotoxic adenosine analogues such as 9-(beta)-D arabinofuranosyladenine (ara-A) and 9-(beta)-D xylofuranosyladenine (xyl-A) can be specifically detoxified by the action of ADA, these analogues were used to select for JEG-3 derived cells which had increased ADA expression. When JEG-3 cells were subjected to a multi-step, successively increasing dosage of either ara-A or xyl-A, resistant cells with increased ADA expression were generated. This increased ADA expression in the resistant cells was unstable, so that when the selective pressure was removed, cellular ADA expression would decrease. Subclone analysis of xyl-A resistant cells revealed that compared to parental JEG-3 cells, individual resistant cells had either elevated ADA levels or decreased adenosine kinase (ADK) levels or both. This altered ADA and ADK expression in the resistant cells were found to be independent events. Because of high endogenous tissue conversion factor (TCF) expression in the JEG-3 cells, the allelic nature of the increased ADA expression in most of the resistant cells could not be determined. However, several resistant subcloned cells were found to have lost TCF expression. These TCF('-) cells expressed only the ADA*2 allelic gene product. Cell fusion experiments demonstrated that the ADA*1 allelic gene was intact and functional in the A3-1A7 cell line. Chromosomal analysis of the A3-1A7 cells showed that they had no double-minutes or homogeneously staining chromosomal regions, although a pair of new chromosomes were found in these cells. Segregation analysis of the hybrid cells indicated that an ADA*2 allelic gene was probably located on this new chromosome. The analysis of the A3-1A7 cell line suggested that the expression of only ADA 2 in these cells was the result of possibly a cis-deregulation of the ADA gene locus or more probably an amplification of the ADA*2 allelic gene. Two effective positive selection systems for ADA('+) cells were also developed and tested. These selection systems should eventually lead to the isolation of the ADA gene.^

Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^

Conditional expression of Sry, and Wnt4 by gene-targeting: Studies in sexual and skeletal development

Sry and Wnt4 cDNAs were individually introduced into the ubiquitously-expressed Rosa26 ( R26) locus by gene targeting in embryonic stem (ES) cells to create a conditional gene expression system in mice. In the targeted alleles, expression of these cDNAs should be blocked by a neomycin resistance selection cassette that is flanked by loxP sites. Transgene expression should be activated after the blocking cassette is deleted by Cre recombinase. ^ To test this conditional expression system, I have bred R26-stop- Sry and R26-stop-Wnt4 heterozygotes with a MisRII-Cre mouse line that expresses Cre in the gonads of both sexes. Analysis of these two types of bigenic heterozygotes indicated that their gonads developed normally like those of wild types. However, one XX R26-Sry/R26-Sry; MisR2-Cre/+ showed epididymis-like structures resembling those of males. In contrast, only normal phenotypes were observed in XY R26-Wnt4/R26-Wnt4; MisR2-Cre /+ mice. To interpret these results, I have tested for Cre recombinase activity by Southern blot and transcription of the Sry and Wnt4 transgenes by RT-PCR. Results showed that bigenic mutants had insufficient activation of the transgenes in their gonads at E12.5 and E13.5. Therefore, the failure to observe mutant phenotypes may have resulted from low activity of MisR2-Cre recombination at the appropriate time. ^ Col2a1-Cre transgenic mice express Cre in differentiating chondrocytes. R26-Wnt4; Col2a1-Cre bigenic heterozygous mice were found to exhibit a dramatic alteration in growth presumably caused by Wnt4 overexpression during chondrogenesis. R26-Wnt4; Col2a1-Cre mice exhibited dwarfism beginning approximately 10 days after birth. In addition, they also had craniofacial abnormalities, and had delayed ossification of the lumbar vertebrate and pelvic bones. Histological analysis of the growth plates of R26-Wnt4; Col2a1-Cre mice revealed less structural organization and a delay in onset of the primary and secondary ossification centers. Molecular studies confirmed that overexpression of Wnt4 causes decreased proliferation and early maturation of chondrocytes. In addition, R26-Wnt4; Col2a1-Cre mice had decreased expression of vascular endothelial growth factor (VEGF), suggesting that defects in vascularization may contribute to the dwarf phenotype. Finally, 9-month-old R26-Wnt4; Col2a1-Cre mice had significantly more fat cells in the marrow cavities of their metaphysis long bones, implying that long-term overexpression of Wnt4may cause bone marrow pathologies. In conclusion, Wnt4 was activated by Col2a1-Cre recombinase and was overexpressed in the growth plate, resulting in aberrant proliferation and differentiation of chondrocytes, and ultimately leads to dwarfism in mice. ^

Network-based analysis of genome-wide scans of natural selection

The genomic era brought by recent advances in the next-generation sequencing technology makes the genome-wide scans of natural selection a reality. Currently, almost all the statistical tests and analytical methods for identifying genes under selection was performed on the individual gene basis. Although these methods have the power of identifying gene subject to strong selection, they have limited power in discovering genes targeted by moderate or weak selection forces, which are crucial for understanding the molecular mechanisms of complex phenotypes and diseases. Recent availability and rapid completeness of many gene network and protein-protein interaction databases accompanying the genomic era open the avenues of exploring the possibility of enhancing the power of discovering genes under natural selection. The aim of the thesis is to explore and develop normal mixture model based methods for leveraging gene network information to enhance the power of natural selection target gene discovery. The results show that the developed statistical method, which combines the posterior log odds of the standard normal mixture model and the Guilt-By-Association score of the gene network in a naïve Bayes framework, has the power to discover moderate/weak selection gene which bridges the genes under strong selection and it helps our understanding the biology under complex diseases and related natural selection phenotypes.^