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.

A whole-genome scan for recurrent airway obstruction in Warmblood sport horses indicates two positional candidate regions

Recurrent airway obstruction (RAO), or heaves, is a naturally occurring asthma-like disease that is related to sensitisation and exposure to mouldy hay and has a familial basis with a complex mode of inheritance. A genome-wide scanning approach using two half-sibling families was taken in order to locate the chromosome regions that contribute to the inherited component of this condition in these families. Initially, a panel of 250 microsatellite markers, which were chosen as a well-spaced, polymorphic selection covering the 31 equine autosomes, was used to genotype the two half-sibling families, which comprised in total 239 Warmblood horses. Subsequently, supplementary markers were added for a total of 315 genotyped markers. Each half-sibling family is focused around a severely RAO-affected stallion, and the phenotype of each individual was assessed for RAO and related signs, namely, breathing effort at rest, breathing effort at work, coughing, and nasal discharge, using an owner-based questionnaire. Analysis using a regression method for half-sibling family structures was performed using RAO and each of the composite clinical signs separately; two chromosome regions (on ECA13 and ECA15) showed a genome-wide significant association with RAO at P < 0.05. An additional 11 chromosome regions showed a more modest association. This is the first publication that describes the mapping of genetic loci involved in RAO. Several candidate genes are located in these regions, a number of which are interleukins. These are important signalling molecules that are intricately involved in the control of the immune response and are therefore good positional candidates.

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.

Genetic mapping of the belt pattern in Brown Swiss cattle to BTA3

The white belt pattern of Brown Swiss cattle is characterized by a lack of melanocytes in a stretch of skin around the midsection. This pattern is of variable width and sometimes the belt does not fully circle the body. To identify the gene responsible for this colour variation, we performed linkage mapping of the belted locus using six segregating half-sib families including 104 informative meioses for the belted character. The pedigree confirmed a monogenic autosomal dominant inheritance of the belted phenotype in Brown Swiss cattle. We performed a genome scan using 186 microsatellite markers in a subset of 88 animals of the six families. Linkage with the belt phenotype was detected at the telomeric region of BTA3. Fine-mapping and haplotype analysis using 19 additional markers in this region refined the critical region of the belted locus to a 922-kb interval on BTA3. As the corresponding human and mouse chromosome segments contain no obvious candidate gene for this coat colour trait, the mutation causing the belt pattern in the Brown Swiss cattle might help to identify an unknown gene influencing skin pigmentation.

Patentability and Scope of Protection for DNA Sequence

Since the mapping of the human genome and the technical innovations in the field of biotechnology, patent law has gone through great controversies. Protection is required for an investor to make an investment but how broad should the given protection be? Whether the invention is a mi- cro-organism capable of dissolving crude oil, or the gene of a soya plant, the genetic engineering required for their production entails vast amounts of capi- tal. The policy in that respect is tailored by legislative acts and judicial decisions, ensuring a fair balance be- tween the interests of patent right holders and third parties. However, the policy differs from jurisdiction to jurisdiction, thus creating inconsistencies with re- gards to the given protection to the same invention, and as a result this could deter innovation and pro- mote stagnation. The most active actors shaping the patent policy on an international level are the patent offices of the United States of America, Japan and the European Patent Organization. These three patent offices have set up a cooperation programme in order to promote and improve efficiency with regards to their patent policies on a global scale. However, recent judicial de- velopments have shown that the policy in respect to the field of biotechnology differs between the patent regimes of the United States of America and the two- layer system of the European Patent Organisation/ the European Union.

HLA-DPB1 and DPB2 are genetic loci for systemic sclerosis: a genome-wide association study in Koreans with replication in North Americans.

OBJECTIVE: To identify systemic sclerosis (SSc) susceptibility loci via a genome-wide association study. METHODS: A genome-wide association study was performed in 137 patients with SSc and 564 controls from Korea using the Affymetrix Human SNP Array 5.0. After fine-mapping studies, the results were replicated in 1,107 SSc patients and 2,747 controls from a US Caucasian population. RESULTS: The single-nucleotide polymorphisms (SNPs) (rs3128930, rs7763822, rs7764491, rs3117230, and rs3128965) of HLA-DPB1 and DPB2 on chromosome 6 formed a distinctive peak with log P values for association with SSc susceptibility (P=8.16x10(-13)). Subtyping analysis of HLA-DPB1 showed that DPB1*1301 (P=7.61x10(-8)) and DPB1*0901 (P=2.55x10(-5)) were the subtypes most susceptible to SSc in Korean subjects. In US Caucasians, 2 pairs of SNPs, rs7763822/rs7764491 and rs3117230/rs3128965, showed strong association with SSc patients who had either circulating anti-DNA topoisomerase I (P=7.58x10(-17)/4.84x10(-16)) or anticentromere autoantibodies (P=1.12x10(-3)/3.2x10(-5)), respectively. CONCLUSION: The results of our genome-wide association study in Korean subjects indicate that the region of HLA-DPB1 and DPB2 contains the loci most susceptible to SSc in a Korean population. The confirmatory studies in US Caucasians indicate that specific SNPs of HLA-DPB1 and/or DPB2 are strongly associated with US Caucasian patients with SSc who are positive for anti-DNA topoisomerase I or anticentromere autoantibodies.

Senescence-related gene expression profiles of rosette leaves of Arabidopsis thaliana: Leaf age versus plant age

Senescence is a form of programmed cell death (PCD) which leads to the death of whole organs, e.g., leaves or flowers, and eventually to the death of entire plants. Like all forms of PCD, senescence is a highly regulated and energy consuming process. Senescence parameters, like protein content, chlorophyll content, expression of photosynthesis-associated genes or senescence-associated genes (SAGs), reveal that senescence occurs in old leaves derived from young plants (6 week old) as well as in young leaves derived from older plants (8 week old), indicating that it is governed by the actual age of the leaves. in order to analyse the differential gene expression profiles during leaf senescence, hybridizations of high-density genome arrays were performed with: i) individual leaves within the rosette of a 6-week-old plant and ii) leaves of the same position within the rosette but harvested from plants of different ages, ranging from 5 to 8 weeks. Cluster and genetree analyses, according to the expression pattern revealed that genes which are up-regulated with respect to the age of the entire plant, showed completely different expression profiles with respect to the age of the individual leaves within one rosette. This was observed even though the actual difference in leaf age was approximately the same. This indicates that gene expression appears to be governed by different parameters: i) the age of the individual leaf and ii) the age and developmental stage of the entire plant.

Replication and fine-mapping of a QTL for recurrent airway obstruction in European Warmblood horses.

Recurrent airway obstruction (RAO), or 'heaves', is a common performance-limiting allergic respiratory disease of mature horses. It is related to sensitization and exposure to mouldy hay and has a familial basis with a complex mode of inheritance. In a previous study, we detected a QTL for RAO on ECA 13 in a half-sib family of European Warmblood horses. In this study, we genotyped additional markers in the family and narrowed the QTL down to about 1.5 Mb (23.7-25.2 Mb). We detected the strongest association with SNP BIEC2-224511 (24,309,405 bp). We also obtained SNP genotypes in an independent cohort of 646 unrelated Warmblood horses. There was no genome-wide significant association with RAO in these unrelated horses. However, we performed a genotypic association study of the SNPs on ECA 13 in these unrelated horses, and the SNP BIEC2-224511 also showed the strongest association with RAO in the unrelated horses (p(raw) = 0.00037). The T allele at this SNP was associated with RAO both in the family and the unrelated horses. Thus, the association study in the unrelated animals provides independent support for the previously detected QTL. The association study allows further narrowing of the QTL interval to about 0.5 Mb (24.0-24.5 Mb). We sequenced the coding regions of the genes in the critical region but did not find any associated coding variants. Therefore, the causative variant underlying this QTL is likely to be a regulatory mutation.

Molecular cloning, expression analysis and assignment of the porcine tumor necrosis factor superfamily member 10 gene (TNFSF10) to SSC13q34-->q36 by fluorescence in situ hybridization and radiation hybrid mapping

We have cloned the complete coding region of the porcine TNFSF10 gene. The porcine TNFSF10 cDNA has an ORF of 870 nucleotides and shares 85% identity with human TNFSF10, and 75% and 72% identity with rat and mouse Tnfsf10 coding sequences, respectively. The deduced porcine TNFSF10 protein consists of 289 amino acids with the calculated molecular mass of 33.5 kDa and a predicted pI of 8.15. The amino acid sequence similarities correspond to 86, 72 and 70% when compared with human, rat and mouse sequences, respectively. Northern blot analysis detected TNFSF10-specific transcripts (approximately 1.7 kb) in various organs of a 10-week-old pig, suggesting ubiquitous expression. Real-time RT-PCR studies of various organs from fetal (days 73 and 98) and postnatal stages (two weeks, eight months) demonstrated developmental and tissue-specific regulation of TNFSF10 mRNA abundance. The chromosomal location of the porcine TNFSF10 gene was determined by FISH of a specific BAC clone to metaphase chromosomes. This TNFSF10 BAC clone has been assigned to SSC13q34-->q36. Additionally, the localization of the TNFSF10 gene was verified by RH mapping on the porcine IMpRH panel.

Simple measures of climate, soil properties and plant traits predict national-scale grassland soil carbon stocks

Soil carbon (C) storage is a key ecosystem service. Soil C stocks play a vital role in soil fertility and climate regulation, but the factors that control these stocks at regional and national scales are unknown, particularly when their composition and stability are considered. As a result, their mapping relies on either unreliable proxy measures or laborious direct measurements. Using data from an extensive national survey of English grasslands, we show that surface soil (0–7 cm) C stocks in size fractions of varying stability can be predicted at both regional and national scales from plant traits and simple measures of soil and climatic conditions. Soil C stocks in the largest pool, of intermediate particle size (50–250 μm), were best explained by mean annual temperature (MAT), soil pH and soil moisture content. The second largest C pool, highly stable physically and biochemically protected particles (0·45–50 μm), was explained by soil pH and the community abundance-weighted mean (CWM) leaf nitrogen (N) content, with the highest soil C stocks under N-rich vegetation. The C stock in the small active fraction (250–4000 μm) was explained by a wide range of variables: MAT, mean annual precipitation, mean growing season length, soil pH and CWM specific leaf area; stocks were higher under vegetation with thick and/or dense leaves. Testing the models describing these fractions against data from an independent English region indicated moderately strong correlation between predicted and actual values and no systematic bias, with the exception of the active fraction, for which predictions were inaccurate. Synthesis and applications. Validation indicates that readily available climate, soils and plant survey data can be effective in making local- to landscape-scale (1–100 000 km2) soil C stock predictions. Such predictions are a crucial component of effective management strategies to protect C stocks and enhance soil C sequestration.

Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted re-sequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrl änder dogs. To our knowledge, this is the first time this family trio WGS-approach, has successfully been used to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrl änder dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G>C (p.R52P) was found to be concordant in eight additional cases and 16 healthy Kromfohrl änder dogs.

Genome characterization, prevalence and distribution of a Macula-like virus from Apis mellifera and Varroa destructor

Around 14 distinct virus species-complexes have been detected in honeybees, each with one or more strains or sub-species. Here we present the initial characterization of an entirely new virus species-complex discovered in honeybee (Apis mellifera L.) and varroa mite (Varroa destructor) samples from Europe and the USA. The virus has a naturally poly-adenylated RNA genome of about 6500 nucleotides with a genome organization and sequence similar to the Tymoviridae (Tymovirales; Tymoviridae), a predominantly plant-infecting virus family. Literature and laboratory analyses indicated that the virus had not previously been described. The virus is very common in French apiaries, mirroring the results from an extensive Belgian survey, but could not be detected in equally-extensive Swedish and Norwegian bee disease surveys. The virus appears to be closely linked to varroa, with the highest prevalence found in varroa samples and a clear seasonal distribution peaking in autumn, coinciding with the natural varroa population development. Sub-genomic RNA analyses show that bees are definite hosts, while varroa is a possible host and likely vector. The tentative name of Bee Macula-like virus (BeeMLV) is therefore proposed. A second, distantly related Tymoviridae-like virus was also discovered in varroa transcriptomes, tentatively named Varroa Tymo-like virus (VTLV).

Physical and functional mapping of a novel tumor suppressor locus for prostate cancer within chromosome 10p12.31-q11

Prostate cancer remains the second leading cause of male cancer deaths in the United States, yet the molecular mechanisms underlying this disease remain largely unknown. Cytogenetic and molecular analyses of prostate tumors suggest a consistent association with the loss of chromosome 10. Previously, we have defined a novel tumor suppressor locus PAC-1 within chromosome 10pter-q11. Introduction of the short arm of chromosome 10 into a prostatic adenocarcinoma cell line PC-3H resulted in dramatic tumor suppression and restoration of a programmed cell death pathway. Using a combined approach of comparative genomic hybridization and microsatellite analysis of PC-3H, I have identified a region of hemizygosity within 10p12-p15. This region has been shown to be involved in frequent loss of heterozygosity in gliomas and melanoma. To functionally dissect the region within chromosome 10p containing PAC-1, we developed a strategy of serial microcell fusion, a technique that allows the transfer of defined fragments of chromosome 10p into PC-3H. Serial microcell fusion was used to transfer defined 10p fragments into a mouse A9 fibrosarcoma cell line. Once characterized by FISH and microsatellite analyses, the 10p fragments were subsequently transferred into PC-3H to generate a panel of microcell hybrid clones containing overlapping deletions of chromosome 10p. In vivo and microsatellite analyses of these PC hybrids identified a small chromosome 10p fragment (an estimated 31 Mb in size inclusive of the centromere) that when transferred into the PC-3H background, resulted in significant tumor suppression and limited a region of functional tumor suppressor activity to chromosome 10p12.31-q11. This region coincides with a region of LOH demonstrated in prostate cancer. These studies demonstrate the utility of this approach as a powerful tool to limit regions of functional tumor suppressor activity. Furthermore, these data used in conjunction with data generated by the Human Genome Project lent a focused approach to identify candidate tumor suppressor genes involved in prostate cancer. ^