Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Genetic diversity among Mycoplasma species bovine group 7: clonal isolates from an outbreak of polyarthritis, mastitis, and abortion in dairy cattle

A comprehensive genetic analysis of 60 Mycoplasma sp. bovine group 7 isolates from different geographic origins and epidemiological settings is presented. Twenty-four isolates were recovered from the joints of calves during sporadic episodes of polyarthritis in geographically distinct regions of Queensland and New South Wales, Australia, including two clones of the type strain PG5O. A further three Australian isolates were also recovered from the tympanic bulla, retropharyngeal lymph node and the lung and another three isolates had unconfirmed histories. Six isolates originated from Germany, Portugal, Nigeria, and France. Twenty-four epidemiologically related isolates of Mycoplasma sp. bovine group 7 were recovered from multiple tissue sites and body fluids of infected calves with polyarthritis, mastitic milk, and from the stomach contents, lung and liver from aborted foetuses in three large, centrally managed dairy herds in New South Wales, Australia. Restriction endonuclease analysis (REA) of genomic DNA differentiated 29 Cfol profiles among these 60 isolates and grouped all 24 epidemiologically related isolates in a defined pattern showing a clonal origin. Three isolates of this clonal cluster were recovered from mastitic milk and the synovial exudate of clinically-affected calves and appeared sporadically for periods up to 18 months after the initial outbreak of polyarthritis indicating a persistent, close association of the organism with cattle in these herds. The Cfol profile representative of the clonal cluster was distinguishable from profiles of isolates recovered from multiple, unrelated cases of polyarthritis in Queensland and New South Wales and from other countries. All 24 isolates from the clonal cluster possessed a plasmid (pBG7AU) with a molecular size of 1022 bp. DNA sequence analysis of pBG7AU identified two open reading frames sharing 81 and 99% DNA sequence similarity with hypothetical replication control proteins A and B respectively, previously described in plasmid pADB201 isolated from M. mycoides subspecies mycoides. Other isolates of bovine group 7, epidemiologically unrelated to the clonal cluster, including two clones of the type strain PG5O, possessed a similar-sized plasmid. These data confirm that Mycoplasma sp. bovine group 7 is capable of migrating to, and multiplying within, different tissue sites within a single animal and among different animals within a herd.

Genetic testing for glucokinase mutations in clinically selected patients with MODY: a worthwhile investment

The differential diagnosis for children with diabetes includes a group of monogenic diabetic disorders known as maturity-onset diabetes of the young (MODY). So far, six underlying gene defects have been identified. The most common subtypes are caused by mutations in the genes encoding the transcription factor HNF-1a (MODY 3) and the glycolytic enzyme glucokinase (GCK) (MODY 2). MODY 2 is the most benign form of diabetes as the threshold for glucose sensing is elevated resulting in mild, regulated hyperglycemia. MODY 2 may usually be treated with diet alone without risk of microvascular complications. Patients with MODY usually present as children or young adults. Genetic testing for MODY in diabetic subjects is often not performed because of the costs and its unavailability in Switzerland. We describe the impact of the genetic analysis for MODY 2 on diabetes management and treatment costs in a five-year-old girl. The patient and her diabetic mother were both found to have a heterozygous missense mutation (V203A) in the glucokinase gene. The five-year-old girl was started on insulin therapy for her diabetes but because her HbA1c remained between 5.8-6.4% (reference 4.1-5.7%) and her clinical presentation suggested MODY insulin was discontinued. She is now well controlled on a carbohydrate controlled diet regimen only. Omission of insulin treatment made regular blood glucose monitoring unnecessary and removed her risk of hypoglycemia. Costs for the genetic analysis were 500 Euro. At our centre costs for diabetes care of a patient with type 1 diabetes are approximately 2050 Euro/year compared to 410 Euro/year for the care of a patient with MODY 2. In addition, a diagnosis of MODY 2 may reassure patients and their families, as microvascular complications are uncommon. Thus there are both health and financial benefits in diagnosing MODY 2. We recommend genetic testing for MODY 2 in clinically selected patients even though this analysis is currently not available in Switzerland and costs are not necessarily covered by the health insurances.

Clinical, genetic, and functional characterization of adrenocorticotropin receptor mutations using a novel receptor assay.

The ACTH receptor (MC2R) is expressed predominantly in the adrenal cortex, but is one of five G protein-coupled, seven-transmembrane melanocortin receptors (MCRs), all of which bind ACTH to some degree. Testing of MC2R activity is difficult because most cells express endogenous MCRs; hence, ACTH will elicit background activation of assayable reporter systems. Inactivating mutations of MC2R lead to hereditary unresponsiveness to ACTH, also known as familial glucocorticoid deficiency (FGD). These patients are usually seen in early childhood with very low cortisol concentrations, normal mineralocorticoids, hyperpigmentation, and increased bodily growth. Several MC2R mutations have been reported in FGD, but assays of the activities of these mutants are cumbersome. We saw two patients with typical clinical findings of FGD. Genetic analysis showed that patient 1 was homozygous for the mutation R137W, and patient 2 was a compound heterozygote for S74I and Y254C. We tested the activity of these mutations in OS-3 cells, which are unresponsive to ACTH but have intact downstream cAMP signal transduction. OS-3 cells transfected with a cAMP-responsive luciferase reporter plasmid (pCREluc) were unresponsive to ACTH, but cotransfection with a vector expressing human MC2R increased luciferase activity more than 40-fold. Addition of ACTH to cells cotransfected with the pCREluc reporter and wild-type MC2R activated luciferase expression with a 50% effective concentration of 5.5 x 10(-9) M ACTH, which is similar to previously reported values. By contrast, the MC2R mutant R137W had low activity, and the S74I or Y254C mutants elicited no measurable response. This assay provides excellent sensitivity in an easily assayed transient transfection system, providing a more rapid and efficient measurement of ACTH receptor activity.

{\it Hox\/} genes and specification of regional identity in the axial skeleton: Analysis of the individual and combined mutant phenotypes of the paralogous group 4 murine {\it Hox\/} genes; {\it hoxa\/}-4, {\it hoxb\/}-4 and {\it hoxd\/}-4

The Hox gene products are transcription factors involved in specifying regional identity along the anteroposterior body axis. In Drosophila, where these genes are known as HOM-C (Homeotic-complex) genes and where they have been most extensively studied, they are expressed in restricted domains along the anteroposterior axis with different anterior limits. Genetic analysis of a large number of gain- and loss-of-function alleles of these genes has revealed that these genes are important in specifying segmental identity at their anterior limits of expression. Furthermore, there is a functional dominance of posterior genes over anterior genes, such that posterior genes can dominantly specify their developmental programs in spite of the expression of more anterior genes in the same segment. In the mouse, there are four clusters of HOM-C genes, called Hox genes. Thus, there may be up to four genes, called paralogs, that are more highly homologous to each other and to their Drosophila homolog than they are to the other mouse Hox genes. The single mutants for two paralogous genes, hoxa-4 and hoxd-4, presented in this dissertation, are similar to several other mouse Hox mutants in that they show partial, incompletely penetrant homeotic transformations of vertebrae at their anterior limit of expression. These mutants were then bred with hoxb-4 mutants (Ramirez-Solis, et al. 1993) to generate the three possible double mutant combinations as well as the triple mutant. The skeletal phenotypes of these group 4 Hox compound mutants displayed clear alterations in regional identity, such that a nearly complete transformation towards the morphology of the first cervical vertebra occurs. These results suggest a certain degree of functional redundancy among paralogous genes in specifying regional identity. Furthermore, there was a remarkable dose-dependent increase in the number of vertebrae transformed to a first cervical vertebra identity, including the second through the fifth cervical vertebrae in the triple mutant. Thus, these genes are required in a larger anteroposterior domain than is revealed by the single mutant phenotypes alone, such that multiple mutations in these genes result in transformations of vertebrae that are not at their anterior limit of expression. ^

Microsurgical and laser ablation analysis of leaf positioning and dorsoventral patterning in tomato

Leaves are arranged according to regular patterns, a phenomenon referred to as phyllotaxis. Important determinants of phyllotaxis are the divergence angle between successive leaves, and the size of the leaves relative to the shoot axis. Young leaf primordia are thought to provide positional information to the meristem, thereby influencing the positioning of new primordia and hence the divergence angle. On the contrary, the meristem signals to the primordia to establish their dorsoventral polarity, which is a prerequisite for the formation of a leaf blade. These concepts originate from classical microsurgical studies carried out between the 1920s and the 1970s. Even though these techniques have been abandoned in favor of genetic analysis, the resulting insights remain a cornerstone of plant developmental biology. Here, we employ new microsurgical techniques to reassess and extend the classical studies on phyllotaxis and leaf polarity. Previous experiments have indicated that the isolation of an incipient primordium by a tangential incision caused a change of divergence angle between the two subsequent primordia, indicating that pre-existing primordia influence further phyllotaxis. Here.. we repeat these experiments and compare them with the results of laser ablation of incipient primordia. Furthermore. we explore to what extent the different pre-existing primordia influence the size and position of new organs. and hence phyllotaxis. We propose that the two youngest primordia (P-1 and P-2) are sufficient for the approximate positioning of the incipient primordium (I-1), and therefore for the perpetuation of the generative spiral, whereas the direct contact neighbours of I-1 (P-2 and P-3) control its delimitation and hence its exact size and position. Finally. we report L I specific cell ablation experiments suggesting that the meristem L-1 layer is essential for the dorsoventral patterning of leaf primordia.

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

A genetic characterization of {\it runt\/} activity during {\it Drosophila\/} embryogenesis

The Drosophila melanogaster gene runt encodes a novel transcriptional regulator that was originally identified on the basis of its key role in embryonic pattern formation. For my thesis I undertook a genetic analysis of runt activity to identify loci that interact with this unique transcriptional regulator. Specifically, I screened the genome with deficiencies for loci that interact with runt in a dose-dependent fashion during early embryogenesis. From this screen I discovered a vital dose-dependent interaction between runt and the achaete-scute complex (AS-C). The characterization of this interaction led to the exciting discovery of important roles for runt in sex determination and neurogenesis (Duffy and Gergen 1991, Duffy et al. 1991). I demonstrated that in sex determination runt is necessary for the normal transcriptional activation of the master sex-determining gene Sx1 and has all the properties of an X:A numerator element. I also showed that runt is required during the early stages of neurogenesis for the normal development of a subset of CNS ganglion mother cells and neurons. In addition, the screen, which focused on the identification and characterization of maternal loci that influence the activity of runt during segmentation, identified several new maternal loci, one of which affects the activity of the maternal posterior group genes on embryonic pattern formation. ^

A genetic screen for mutations that disrupt an auditory response in Drosophila melanogaster

Hearing is one of the last sensory modalities to be subjected to genetic analysis in Drosophila melanogaster. We describe a behavioral assay for auditory function involving courtship among groups of males triggered by the pulse component of the courtship song. In a mutagenesis screen for mutations that disrupt the auditory response, we have recovered 15 mutations that either reduce or abolish this response. Mutant audiograms indicate that seven mutants reduced the amplitude of the response at all intensities. Another seven abolished the response altogether. The other mutant, 5L3, responded only at high sound intensities, indicating that the threshold was shifted in this mutant. Six mutants were characterized in greater detail. 5L3 had a general courtship defect; courtship of females by 5L3 males also was affected strongly. 5P1 males courted females normally but had reduced success at copulation. 5P1 and 5N18 showed a significant decrement in olfactory response, indicating that the defects in these mutations are not specific to the auditory pathway. Two other mutants, 5M8 and 5N30, produced amotile sperm although in 5N30 this phenotype was genetically separable from the auditory phenotype. Finally, a new adult circling behavior phenotype, the pirouette phenotype, associated with massive neurodegeneration in the brain, was discovered in two mutants, 5G10 and 5N18. This study provides the basis for a genetic and molecular dissection of auditory mechanosensation and auditory behavior.

Genetic Dissection of Endocytic Trafficking in Drosophila Using a Horseradish Peroxidase-Bride of Sevenless Chimera: hook Is Required for Normal Maturation of Multivesicular Endosomes

Mutations in the hook gene alter intracellular trafficking of internalized ligands in Drosophila. To dissect this defect in more detail, we developed a new approach to visualize the pathway taken by the Bride of Sevenless (Boss) ligand after its internalization into R7 cells. A chimeric protein consisting of HRP fused to Boss (HRP-Boss) was expressed in R8 cells. This chimera was fully functional: it rescued the boss mutant phenotype, and its trafficking was indistinguishable from that of the wild-type Boss protein. The HRP activity of the chimera was used to follow HRP-Boss trafficking on the ultrastructural level through early and late endosomes in R7 cells. In both wild-type and hook mutant eye disks, HRP-Boss was internalized into R7 cells. In wild-type tissue, Boss accumulated in mature multivesicular bodies (MVBs) within R7 cells; such accumulation was not observed in hook eye disks, however. Quantitative electron microscopy revealed a loss of mature MVBs in hook mutant tissue compared with wild type, whereas more than twice as many multilammelar late endosomes were detected. Our genetic analysis indicates that Hook is required late in endocytic trafficking to negatively regulate delivery from mature MVBs to multilammelar late endosomes and lysosomes.

Genomic and genetic dissection of an archaeal regulon

The extremely halophilic archaeon Halobacterium sp. NRC-1 can grow phototrophically by means of light-driven proton pumping by bacteriorhodopsin in the purple membrane. Here, we show by genetic analysis of the wild type, and insertion and double-frame shift mutants of Bat that this transcriptional regulator coordinates synthesis of a structural protein and a chromophore for purple membrane biogenesis in response to both light and oxygen. Analysis of the complete Halobacterium sp. NRC-1 genome sequence showed that the regulatory site, upstream activator sequence (UAS), the putative binding site for Bat upstream of the bacterio-opsin gene (bop), is also present upstream to the other Bat-regulated genes. The transcription regulator Bat contains a photoresponsive cGMP-binding (GAF) domain, and a bacterial AraC type helix–turn–helix DNA binding motif. We also provide evidence for involvement of the PAS/PAC domain of Bat in redox-sensing activity by genetic analysis of a purple membrane overproducer. Five additional Bat-like putative regulatory genes were found, which together are likely to be responsible for orchestrating the complex response of this archaeon to light and oxygen. Similarities of the bop-like UAS and transcription factors in diverse organisms, including a plant and a γ-proteobacterium, suggest an ancient origin for this regulon capable of coordinating light and oxygen responses in the three major branches of the evolutionary tree of life. Finally, sensitivity of four of five regulon genes to DNA supercoiling is demonstrated and correlated to presence of alternating purine–pyrimidine sequences (RY boxes) near the regulated promoters.

The genetic mediation of individual differences in sensitivity to pain and its inhibition

The underlying bases of the considerable interindividual variability in pain-related traits are starting to be revealed. Although the relative importance of genes versus experience in human pain perception remains unclear, rodent populations display large and heritable differences in both nociceptive and analgesic sensitivity. The identification and characterization of particularly divergent populations provides a powerful initial step in the genetic analysis of pain, because these models can be exploited to identify genes contributing to the behavior-level variability. Ultimately, DNA sequence differences representing the differential alleles at pain-relevant genes can be identified. Thus, by using a combination of “top-down” and “bottom-up” strategies, we are now able to genetically dissect even complex biological traits like pain. The present review summarizes the current progress toward these ends in both humans and rodents.

Genetic complexity of pathogen perception by plants: The example of Rcr3, a tomato gene required specifically by Cf-2

Genetic analysis of plant–pathogen interactions has demonstrated that resistance to infection is often determined by the interaction of dominant plant resistance (R) genes and dominant pathogen-encoded avirulence (Avr) genes. It was postulated that R genes encode receptors for Avr determinants. A large number of R genes and their cognate Avr genes have now been analyzed at the molecular level. R gene loci are extremely polymorphic, particularly in sequences encoding amino acids of the leucine-rich repeat motif. A major challenge is to determine how Avr perception by R proteins triggers the plant defense response. Mutational analysis has identified several genes required for the function of specific R proteins. Here we report the identification of Rcr3, a tomato gene required specifically for Cf-2-mediated resistance. We propose that Avr products interact with host proteins to promote disease, and that R proteins “guard” these host components and initiate Avr-dependent plant defense responses.

Multilocus analysis of extracellular putative virulence proteins made by group A Streptococcus: Population genetics, human serologic response, and gene transcription

Species of pathogenic microbes are composed of an array of evolutionarily distinct chromosomal genotypes characterized by diversity in gene content and sequence (allelic variation). The occurrence of substantial genetic diversity has hindered progress in developing a comprehensive understanding of the molecular basis of virulence and new therapeutics such as vaccines. To provide new information that bears on these issues, 11 genes encoding extracellular proteins in the human bacterial pathogen group A Streptococcus identified by analysis of four genomes were studied. Eight of the 11 genes encode proteins with a LPXTG(L) motif that covalently links Gram-positive virulence factors to the bacterial cell surface. Sequence analysis of the 11 genes in 37 geographically and phylogenetically diverse group A Streptococcus strains cultured from patients with different infection types found that recent horizontal gene transfer has contributed substantially to chromosomal diversity. Regions of the inferred proteins likely to interact with the host were identified by molecular population genetic analysis, and Western immunoblot analysis with sera from infected patients confirmed that they were antigenic. Real-time reverse transcriptase–PCR (TaqMan) assays found that transcription of six of the 11 genes was substantially up-regulated in the stationary phase. In addition, transcription of many genes was influenced by the covR and mga trans-acting gene regulatory loci. Multilocus investigation of putative virulence genes by the integrated approach described herein provides an important strategy to aid microbial pathogenesis research and rapidly identify new targets for therapeutics research.