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.

The HCP5 single-nucleotide polymorphism: a simple screening tool for prediction of hypersensitivity reaction to abacavir

The HLA-B 5701 allele is predictive of hypersensitivity reaction to abacavir, a response herein termed "ABC-HSR." This study of 1,103 individuals infected with human immunodeficiency virus assessed the usefulness of genotyping a HCP5 single-nucleotide polymorphism (SNP), rs2395029, in relation to ABC-HSR. In populations with European ancestry, rs2395029 is in linkage disequilibrium with HLA-B 5701. The HCP5 SNP was present in all 98 HLA-B 5701-positive individuals and was absent in 999 of 1005 HLA-B 5701-negative individuals. rs2395029 was overrepresented in 25 individuals with clinically likely ABC-HSR, compared with its frequency in 175 ABC-tolerant individuals (80% vs. 2%, respectively; P < .0001). Therefore, HCP5 genotyping could serve as a simple screening tool for ABC-HSR, particularly in settings where sequence-based HLA typing is not available.

Mapping Individual Variation in Male Mating Preference Space: Multiple Choice in a Color Polymorphic Cichlid Fish

Sexual selection theory largely rests on the assumption that populations contain individual variation in mating preferences and that individuals are consistent in their preferences. However, there are few empirical studies of within-population variation and even fewer have examined individual male mating preferences. Here, we studied a color polymorphic population of the Lake Victoria cichlid fish Neochromis omnicaeruleus, a species in which color morphs are associated with different sex-determining factors. Wild-caught males were tested in three-way choice trials with multiple combinations of different females belonging to the three color morphs. Compositional log-ratio techniques were applied to analyze individual male mating preferences. Large individual variation in consistency, strength, and direction of male mating preferences for female color morphs was found and hierarchical clustering of the compositional data revealed the presence of four distinct preference groups corresponding to the three color morphs in addition to a no-preference class. Consistency of individual male mating preferences was higher in males with strongest preferences. We discuss the implications of these findings for our understanding of the mechanisms underlying polymorphism in mating preferences.

Molecular cytogenetics and gene mapping in sheep (Ovis aries, 2n = 54)

The development of a completely annotated sheep genome sequence is a key need for understanding the phylogenetic relationships and genetic diversity among the many different sheep breeds worldwide and for identifying genes controlling economically and physiologically important traits. The ovine genome sequence assembly will be crucial for developing optimized breeding programs based on highly productive, healthy sheep phenotypes that are adapted to modern breeding and production conditions. Scientists and breeders around the globe have been contributing to this goal by generating genomic and cDNA libraries, performing genome-wide and trait-associated analyses of polymorphism, expression analysis, genome sequencing, and by developing virtual and physical comparative maps. The International Sheep Genomics Consortium (ISGC), an informal network of sheep genomics researchers, is playing a major role in coordinating many of these activities. In addition to serving as an essential tool for monitoring chromosome abnormalities in specific sheep populations, ovine molecular cytogenetics provides physical anchors which link and order genome regions, such as sequence contigs, genes and polymorphic DNA markers to ovine chromosomes. Likewise, molecular cytogenetics can contribute to the process of defining evolutionary breakpoints between related species. The selective expansion of the sheep cytogenetic map, using loci to connect maps and identify chromosome bands, can substantially contribute to improving the quality of the annotated sheep genome sequence and will also accelerate its assembly. Furthermore, identifying major morphological chromosome anomalies and micro-rearrangements, such as gene duplications or deletions, that might occur between different sheep breeds and other Ovis species will also be important to understand the diversity of sheep chromosome structure and its implications for cross-breeding. To date, 566 loci have been assigned to specific chromosome regions in sheep and the new cytogenetic map is presented as part of this review. This review will also summarize the current cytogenomic status of the sheep genome, describe current activities in the sheep cytogenomics research sector, and will discuss the cytogenomics data in context with other major sheep genomics projects.

Unravelling the distinct strains of Tharu ancestry

The northern region of the Indian subcontinent is a vast landscape interlaced by diverse ecologies, e.g. the Gangetic plain and the Himalayas. A great number of ethnic groups are found there, displayed as a multitude of languages and cultures. The Tharu represent one of the largest and linguistically most diversified such groups, scattered across the Tarai region of Nepal and bordering Indian states. Their origins are uncertain. Hypotheses have been advanced about an Austroasiatic affinity, Tibeto-Burman origins, as well as aboriginal roots in the Tarai. Several Tharu groups speak a variety of Indo-Aryan languages, but have traditionally been described by ethnographers as representing an East Asian phenotype. Their ancestry and intra-population diversity had previously been tested only for haploid (mtDNA and Y-chromosome) markers in a small portion of the population. This study presents the first systematic genetic survey of the Tharu from both Nepal and the Indian states of Uttaranchal and Uttar Pradesh, using genome-wide SNPs and haploid (mtDNA and Y-chromosome) markers. The results suggest that the 'ethnic' construct of Tharu is likely to have lain in the Tarai region, with a reconstructible radiation to Uttaranchal and Uttar Pradesh. Despite extensive admixture with other local communities, Tharu sub-populations preserve strong genetic signatures that indicate a common ancestry.

Polymorphisms, linkage and mapping of four enzyme loci in the fish genus Xiphophorus (Poeciliidae).

Electrophoretic variants at four additional enzyme loci--two esterases (Est-2, Est-3), retinal lactate dehydrogenase (LDH-1) and mannose phosphate isomerase (MPI)--among three species and four subspecies of fish of the genus Xiphophorus were observed. Electrophoretic patterns in F1 hybrid heterozygotes confirmed the monomeric structures of MPI and the esterase and the tetrametric structure of LDH in these fishes. Variant alleles of all four loci displayed normal Mendelian segregation in backcross and F2 hybrids. Recombination data from backcross hybrids mapped with Haldane's mapping function indicate the four loci to be linked as Est-2--0.43--Est3--0.26--LDH-1--0.19--MPI. Significant interference was detected and apparently concentrated in the Est-3 to MPI region. No significant sex-specific differences in recombination were observed. This group (designated linkage group II) was shown to assort independently from the three loci of linkage group I (adenosine deaminase, glucose-6-phosphate dehydrogenase, and 6-phosphogluconate dehydrogenase) and from glyceraldehyde-3-phosphate dehydrogenase and two isocitrate dehydrogenase loci. Evidence for conservation of the linkage group, at least in part, in other vertebrate species is presented.

Comparative mapping of Andropogoneae: Saccharum L. (sugarcane) and its relation to sorghum and maize

Comparative genetic maps of Papuan Saccharum officinarum L. (2n = 80) and S. robustum (2n = 80) were constructed by using single-dose DNA markers (SDMs). SDM-framework maps of S. officinarum and S. robustum were compared with genetic maps of sorghum and maize by way of anchor restriction fragment length polymorphism probes. The resulting comparisons showed striking colinearity between the sorghum and Saccharum genomes. There were no differences in marker order between S. officinarum and sorghum. Furthermore, there were no alterations in SDM order between S. officinarum and S. robustum. The S. officinarum and S. robustum maps also were compared with the map of the polysomic octoploid S. spontaneum ‘SES 208’ (2n = 64, x = 8), thus permitting relations to homology groups (“chromosomes”) of S. spontaneum to be studied. Investigation of transmission genetics in S. officinarum and S. robustum confirmed preliminary results that showed incomplete polysomy in these species. Because of incomplete polysomy, multiple-dose markers could not be mapped for lack of a genetic model for their segregation. To coalesce S. officinarum and S. robustum linkage groups into homology groups (composed of homologous pairing partners), they were compared with sorghum (2n = 20), which functioned as a synthetic diploid. Groupings suggested by comparative mapping were found to be highly concordant with groupings based on highly polymorphic restriction fragment length polymorphism probes detecting multiple SDMs. The resulting comparative maps serve as bridges to allow information from one Andropogoneae to be used by another, for breeding, ecology, evolution, and molecular biology.

High-resolution mapping and isolation of a yeast artificial chromosome contig containing fw2.2: A major fruit weight quantitative trait locus in tomato

A high-resolution physical and genetic map of a major fruit weight quantitative trait locus (QTL), fw2.2, has been constructed for a region of tomato chromosome 2. Using an F2 nearly isogenic line mapping population (3472 individuals) derived from Lycopersicon esculentum (domesticated tomato) × Lycopersicon pennellii (wild tomato), fw2.2 has been placed near TG91 and TG167, which have an interval distance of 0.13 ± 0.03 centimorgan. The physical distance between TG91 and TG167 was estimated to be ≤ 150 kb by pulsed-field gel electrophoresis of tomato DNA. A physical contig composed of six yeast artificial chromosomes (YACs) and encompassing fw2.2 was isolated. No rearrangements or chimerisms were detected within the YAC contig based on restriction fragment length polymorphism analysis using YAC-end sequences and anchored molecular markers from the high-resolution map. Based on genetic recombination events, fw2.2 could be narrowed down to a region less than 150 kb between molecular markers TG91 and HSF24 and included within two YACs: YAC264 (210 kb) and YAC355 (300 kb). This marks the first time, to our knowledge, that a QTL has been mapped with such precision and delimited to a segment of cloned DNA. The fact that the phenotypic effect of the fw2.2 QTL can be mapped to a small interval suggests that the action of this QTL is likely due to a single gene. The development of the high-resolution genetic map, in combination with the physical YAC contig, suggests that the gene responsible for this QTL and other QTLs in plants can be isolated using a positional cloning strategy. The cloning of fw2.2 will likely lead to a better understanding of the molecular biology of fruit development and to the genetic engineering of fruit size characteristics.

Centromere mapping and orientation of the molecular linkage map of rice (Oryza sativa L.).

Rice has become a model cereal plant for molecular genetic research. Rice has the most comprehensive molecular linkage maps with more than 2000 DNA markers and shows synteny and colinearity with the maps of other cereal crops. Until now, however, no information was available about the positions of centromeres and arm locations of markers on the molecular linkage map. Secondary and telotrisomics were used to assign restriction fragment length polymorphism markers to specific chromosome arms and thereby to map the positions of centromeres. More than 170 restriction fragment length polymorphism markers were assigned to specific chromosome arms through gene dosage analysis using the secondary and telotrisomics and the centromere positions were mapped on all 12 linkage groups. The orientations of seven linkage groups were reversed to fit the "short arm on top" convention and the corrected map is presented.

Structure of the gene for porcine peptide antibiotic PR-39, a cathelin gene family member: comparative mapping of the locus for the human peptide antibiotic FALL-39.

PR-39 is a porcine 39-aa peptide antibiotic composed of 49% proline and 24% arginine, with an activity against Gram-negative bacteria comparable to that of tetracycline. In Escherichia coli, it inhibits DNA and protein synthesis. PR-39 was originally isolated from pig small intestine, but subsequent cDNA cloning showed that the gene is expressed in the bone marrow. The open reading frame of the clone showed that PR-39 is made as 173-aa precursor whose proregion belongs to the cathelin family. The PR39 gene, which is rather compact and spans only 1784 bp has now been sequenced. The coding information is split into four exons. The first exon contains the signal sequence of 29 residues and the first 37 residues of the cathelin propart. Exons 2 and 3 contain only cathelin information, while exon 4 codes for the four C-terminal cathelin residues and the mature PR-39 peptide extended by three residues. The sequenced upstream region (1183 bp) contains four potential recognition sites for NF-IL6 and three for APRF, transcription factors known to regulate genes for both cytokines and acute phase response factors. Genomic hybridizations revealed a fairly high level of restriction fragment length polymorphism and indicated that there are at least two copies of the PR39 gene in the pig genome. PR39 was mapped to pig chromosome 13 by linkage and in situ hybridization mapping. The gene for the human peptide antibiotic FALL-39 (also a member of the cathelin family) was mapped to human chromosome 3, which is homologous to pig chromosome 13.

Epidermal Growth Factor Gene (EGF) Polymorphism and Risk of Melanocytic Neoplasiay

A common single nucleotide polymorphism (SNP) in the 5' untranslated region (5'UTR) of the epidermal growth factor (EGF) gene modulates the level of transcription of this gene and hence is associated with serum levels of EGF. This variant may be associated with melanoma risk, but conflicting findings have been reported. An Australian melanoma case-control sample was typed for the EGF+61A>G transversion (rs4444903). The sample comprised 753 melanoma cases from 738 families stratified by family history of melanoma and 2387 controls from 645 unselected twin families. Ancestry of the cases and controls was recorded, and the twins had undergone skin examination to assess total body nevus count, degree of freckling and pigmentation phenotype. SNP genotyping was carried out via primer extension followed by matrix-assisted laser desorption time of flight (MALDI-TOF) mass spectroscopy. The EGIF+61 SNP was not found to be significantly associated with melanoma status or with development of nevi or freckles. Among melanoma cases, however, G homozygotes had thicker tumors (p=0.05), in keeping with two previous studies. The EGF polymorphism does not appear to predispose to melanoma or nevus development, but its significant association with tumor thickness implies that it may be a useful marker of prognosis.

Eye colour: portals into pigmentation genes and ancestry

Several recent papers have tried to address the genetic determination of eye colour via microsatellite linkage, testing of pigmentation candidate gene polymorphisms and the genome wide analysis of SNP markers that are informative for ancestry. These studies show that the OCA2 gene on chromosome 15 is the major determinant of brown and/or blue eye colour but also indicate that other loci will be involved in the broad range of hues seen in this trait in Europeans.

Fine mapping of the tomato I-3 gene for fusarium wilt resistance and elimination of a co-segregating resistance gene analogue as a candidate for I-3

The I-3 gene from the wild tomato species Lycopersicon pennellii confers resistance to race 3 of the devastating vascular wilt pathogen Fusarium oxysporum f. sp. lycopersici. As an initial step in a positional cloning strategy for the isolation of I-3, we converted restriction fragment length polymorphism and conserved orthologue set markers, known genes and a resistance gene analogue (RGA) mapping to the I-3 region into PCR-based sequence characterised amplified region (SCAR) and cleaved amplified polymorphic sequence (CAPS) markers. Additional PCR-based markers in the I-3 region were generated using the randomly amplified DNA fingerprinting (RAF) technique. SCAR, CAPS and RAF markers were used for high-resolution mapping around the I-3 locus. The I-3 gene was localised to a 0.3-cM region containing a RAF marker, eO6, and an RGA, RGA332. RGA332 was cloned and found to correspond to a putative pseudogene with at least two loss-of-function mutations. The predicted pseudogene belongs to the Toll interleukin-1 receptor-nucleotide-binding site-leucine-rich-repeat sub-class of plant disease resistance genes. Despite the presence of two RGA332 homologues in L. esculentum, DNA gel blot and PCR analysis suggests that no other homologues are present in lines carrying I-3 that could be alternative candidates for the gene.

Mapping genes for resistance to net form of net blotch and strip rust in barley (Hordeum vulgare L.)

A dihaploid mapping population comprising 65 lines was developed between barley parent varieties Tallon and Kaputar and used to construct a genetic linkage map. This map, comprising 195 amplified fragment length polymorphism and 38 simple sequence repeat markers, was used to identify markers linked to the net form of net blotch (Pyrenophora teres f.sp. teres) and to stripe rust (Puccinia striiformis f.sp. hordei) in barley. The population was screened with five pathotypes of net blotch at the seedling stage in the glasshouse and subjected to a natural inoculation in Hermitage, Queensland. Stripe rust screening was conducted at the adult plant stage in Toluca, Mexico. Analyses of the markers were performed using Mapmanager and Qgene software. One region on chromosome 6H was highly significantly associated with resistance to the net blotch (R2 = 79%). This association was consistent for all pathotypes studied. One region on chromosome 5H was found to be highly significantly associated with resistance to stripe rust (R2= 65%). There are a number of very closely linked markers showing strong associations in these regions, and these markers present an opportunity for marker assisted selection of these traits in barley breeding programs.

The promoter polymorphism -232C/G of the PCK1 gene is associated with type 2 diabetes in a UK-resident South Asian population

Background - The PCK1 gene, encoding cytosolic phosphoenolpyruvate carboxykinase (PEPCKC), has previously been implicated as a candidate gene for type 2 diabetes (T2D) susceptibility. Rodent models demonstrate that over-expression of Pck1 can result in T2D development and a single nucleotide polymorphism (SNP) in the promoter region of human PCK1 (-232C/G) has exhibited significant association with the disease in several cohorts. Within the UK-resident South Asian population, T2D is 4 to 6 times more common than in indigenous white Caucasians. Despite this, few studies have reported on the genetic susceptibility to T2D in this ethnic group and none of these has investigated the possible effect of PCK1 variants. We therefore aimed to investigate the association between common variants of the PCK1 gene and T2D in a UK-resident South Asian population of Punjabi ancestry, originating predominantly from the Mirpur area of Azad Kashmir, Pakistan. Methods - We used TaqMan assays to genotype five tagSNPs covering the PCK1 gene, including the -232C/G variant, in 903 subjects with T2D and 471 normoglycaemic controls. Results - Of the variants studied, only the minor allele (G) of the -232C/G SNP demonstrated a significant association with T2D, displaying an OR of 1.21 (95% CI: 1.03 - 1.42, p = 0.019). Conclusion - This study is the first to investigate the association between variants of the PCK1 gene and T2D in South Asians. Our results suggest that the -232C/G promoter polymorphism confers susceptibility to T2D in this ethnic group.