A MULTILEVEL MODEL TO ADDRESS BATCH EFFECTS IN COPY NUMBER ESTIMATION USING SNP ARRAYS

Submicroscopic changes in chromosomal DNA copy number dosage are common and have been implicated in many heritable diseases and cancers. Recent high-throughput technologies have a resolution that permits the detection of segmental changes in DNA copy number that span thousands of basepairs across the genome. Genome-wide association studies (GWAS) may simultaneously screen for copy number-phenotype and SNP-phenotype associations as part of the analytic strategy. However, genome-wide array analyses are particularly susceptible to batch effects as the logistics of preparing DNA and processing thousands of arrays often involves multiple laboratories and technicians, or changes over calendar time to the reagents and laboratory equipment. Failure to adjust for batch effects can lead to incorrect inference and requires inefficient post-hoc quality control procedures that exclude regions that are associated with batch. Our work extends previous model-based approaches for copy number estimation by explicitly modeling batch effects and using shrinkage to improve locus-specific estimates of copy number uncertainty. Key features of this approach include the use of diallelic genotype calls from experimental data to estimate batch- and locus-specific parameters of background and signal without the requirement of training data. We illustrate these ideas using a study of bipolar disease and a study of chromosome 21 trisomy. The former has batch effects that dominate much of the observed variation in quantile-normalized intensities, while the latter illustrates the robustness of our approach to datasets where as many as 25% of the samples have altered copy number. Locus-specific estimates of copy number can be plotted on the copy-number scale to investigate mosaicism and guide the choice of appropriate downstream approaches for smoothing the copy number as a function of physical position. The software is open source and implemented in the R package CRLMM available at Bioconductor (http:www.bioconductor.org).

EXPLORATION, NORMALIZATION, AND GENOTYPE CALLS OF HIGH DENSITY OLIGONUCLEOTIDE SNP ARRAY DATA

In most microarray technologies, a number of critical steps are required to convert raw intensity measurements into the data relied upon by data analysts, biologists and clinicians. These data manipulations, referred to as preprocessing, can influence the quality of the ultimate measurements. In the last few years, the high-throughput measurement of gene expression is the most popular application of microarray technology. For this application, various groups have demonstrated that the use of modern statistical methodology can substantially improve accuracy and precision of gene expression measurements, relative to ad-hoc procedures introduced by designers and manufacturers of the technology. Currently, other applications of microarrays are becoming more and more popular. In this paper we describe a preprocessing methodology for a technology designed for the identification of DNA sequence variants in specific genes or regions of the human genome that are associated with phenotypes of interest such as disease. In particular we describe methodology useful for preprocessing Affymetrix SNP chips and obtaining genotype calls with the preprocessed data. We demonstrate how our procedure improves existing approaches using data from three relatively large studies including one in which large number independent calls are available. Software implementing these ideas are avialble from the Bioconductor oligo package.

EFFICIENT EVALUATION OF RANKING PROCEDURES WHEN THE NUMBER OF UNITS IS LARGE WITH APPLICATION TO SNP IDENTIFICATION

Simulation-based assessment is a popular and frequently necessary approach to evaluation of statistical procedures. Sometimes overlooked is the ability to take advantage of underlying mathematical relations and we focus on this aspect. We show how to take advantage of large-sample theory when conducting a simulation using the analysis of genomic data as a motivating example. The approach uses convergence results to provide an approximation to smaller-sample results, results that are available only by simulation. We consider evaluating and comparing a variety of ranking-based methods for identifying the most highly associated SNPs in a genome-wide association study, derive integral equation representations of the pre-posterior distribution of percentiles produced by three ranking methods, and provide examples comparing performance. These results are of interest in their own right and set the framework for a more extensive set of comparisons.

Assessment of use of microsatellite polymorphism analysis for improving spatial distribution tracking of echinococcus multilocularis

Alveolar echinococcosis (AE)--caused by the cestode Echinococcus multilocularis--is a severe zoonotic disease found in temperate and arctic regions of the northern hemisphere. Even though the transmission patterns observed in different geographical areas are heterogeneous, the nuclear and mitochondrial targets usually used for the genotyping of E. multilocularis have shown only a marked genetic homogeneity in this species. We used microsatellite sequences, because of their high typing resolution, to explore the genetic diversity of E. multilocularis. Four microsatellite targets (EmsJ, EmsK, and EmsB, which were designed in our laboratory, and NAK1, selected from the literature) were tested on a panel of 76 E. multilocularis samples (larval and adult stages) obtained from Alaska, Canada, Europe, and Asia. Genetic diversity for each target was assessed by size polymorphism analysis. With the EmsJ and EmsK targets, two alleles were found for each locus, yielding two and three genotypes, respectively, discriminating European isolates from the other groups. With NAK1, five alleles were found, yielding seven genotypes, including those specific to Tibetan and Alaskan isolates. The EmsB target, a tandem repeated multilocus microsatellite, found 17 alleles showing a complex pattern. Hierarchical clustering analyses were performed with the EmsB findings, and 29 genotypes were identified. Due to its higher genetic polymorphism, EmsB exhibited a higher discriminatory power than the other targets. The complex EmsB pattern was able to discriminate isolates on a regional and sectoral level, while avoiding overdistinction. EmsB will be used to assess the putative emergence of E. multilocularis in Europe.

Beta-D-glucoside utilization by Mycoplasma mycoides subsp. mycoides SC: possible involvement in the control of cytotoxicity towards bovine lung cells

BACKGROUND: Contagious bovine pleuropneumonia (CBPP) caused by Mycoplasma mycoides subsp. mycoides small-colony type (SC) is among the most serious threats for livestock producers in Africa. Glycerol metabolism-associated H2O2 production seems to play a crucial role in virulence of this mycoplasma. A wide number of attenuated strains of M. mycoides subsp. mycoides SC are currently used in Africa as live vaccines. Glycerol metabolism is not affected in these vaccine strains and therefore it does not seem to be the determinant of their attenuation. A non-synonymous single nucleotide polymorphism (SNP) in the bgl gene coding for the 6-phospho-beta-glucosidase (Bgl) has been described recently. The SNP differentiates virulent African strains isolated from outbreaks with severe CBPP, which express the Bgl isoform Val204, from strains to be considered less virulent isolated from CBPP outbreaks with low mortality and vaccine strains, which express the Bgl isoform Ala204. RESULTS: Strains of M. mycoides subsp. mycoides SC considered virulent and possessing the Bgl isoform Val204, but not strains with the Bgl isoform Ala204, do trigger elevated levels of damage to embryonic bovine lung (EBL) cells upon incubation with the disaccharides (i.e., beta-D-glucosides) sucrose and lactose. However, strains expressing the Bgl isoform Val204 show a lower hydrolysing activity on the chromogenic substrate p-nitrophenyl-beta-D-glucopyranoside (pNPbG) when compared to strains that possess the Bgl isoform Ala204. Defective activity of Bgl in M. mycoides subsp. mycoides SC does not lead to H2O2 production. Rather, the viability during addition of beta-D-glucosides in medium-free buffers is higher for strains harbouring the Bgl isoform Val204 than for those with the isoform Ala204. CONCLUSION: Our results indicate that the studied SNP in the bgl gene is one possible cause of the difference in bacterial virulence among strains of M. mycoides subsp. mycoides SC. Bgl does not act as a direct virulence factor, but strains possessing the Bgl isoform Val204 with low hydrolysing activity are more prone to survive in environments that contain high levels of beta-D-glucosides, thus contributing in some extent to mycoplasmaemia.

Functional polymorphism in ABCA1 influences age of symptom onset in coronary artery disease patients

ATP-binding-cassette-transporter-A1 (ABCA1) plays a pivotal role in intracellular cholesterol removal, exerting a protective effect against atherosclerosis. ABCA1 gene severe mutations underlie Tangier disease, a rare Mendelian disorder that can lead to premature coronary artery disease (CAD), with age of CAD onset being two decades earlier in mutant homozygotes and one decade earlier in heterozygotes than in mutation non-carriers. It is unknown whether common polymorphisms in ABCA1 could influence age of symptom onset of CAD in the general population. We examined common promoter and non-synonymous coding polymorphisms in relation to age of symptom onset in a group of CAD patients (n = 1164), and also carried out in vitro assays to test effects of the promoter variations on ABCA1 promoter transcriptional activity and effects of the coding variations on ABCA1 function in mediating cellular cholesterol efflux. Age of symptom onset was found to be associated with the promoter - 407G > C polymorphism, being 2.82 years higher in C allele homozygotes than in G allele homozygotes and intermediate in heterozygotes (61.54, 59.79 and 58.72 years, respectively; P = 0.002). In agreement, patients carrying ABCA1 haplotypes containing the -407C allele had higher age of symptom onset. Patients of the G/G or G/C genotype of the -407G > C polymorphism had significant coronary artery stenosis (>75%) at a younger age than those of the C/C genotype (P = 0.003). Reporter gene assays showed that ABCA1 haplotypes bearing the -407C allele had higher promoter activity than haplotypes with the -407G allele. Functional analyses of the coding polymorphisms showed an effect of the V825I substitution on ABCA1 function, with the 825I variant having higher activity in mediating cholesterol efflux than the wild-type (825V). A trend towards higher symptom onset age in 825I allele carriers was observed. The data indicate an influence of common ABCA1 functional polymorphisms on age of symptom onset in CAD patients.

Heterosis for biomass-related traits in Arabidopsis investigated by quantitative trait loci analysis of the triple testcross design with recombinant inbred lines

Arabidopsis thaliana has emerged as a leading model species in plant genetics and functional genomics including research on the genetic causes of heterosis. We applied a triple testcross (TTC) design and a novel biometrical approach to identify and characterize quantitative trait loci (QTL) for heterosis of five biomass-related traits by (i) estimating the number, genomic positions, and genetic effects of heterotic QTL, (ii) characterizing their mode of gene action, and (iii) testing for presence of epistatic effects by a genomewide scan and marker x marker interactions. In total, 234 recombinant inbred lines (RILs) of Arabidopsis hybrid C24 x Col-0 were crossed to both parental lines and their F1 and analyzed with 110 single-nucleotide polymorphism (SNP) markers. QTL analyses were conducted using linear transformations Z1, Z2, and Z3 calculated from the adjusted entry means of TTC progenies. With Z1, we detected 12 QTL displaying augmented additive effects. With Z2, we mapped six QTL for augmented dominance effects. A one-dimensional genome scan with Z3 revealed two genomic regions with significantly negative dominance x additive epistatic effects. Two-way analyses of variance between marker pairs revealed nine digenic epistatic interactions: six reflecting dominance x dominance effects with variable sign and three reflecting additive x additive effects with positive sign. We conclude that heterosis for biomass-related traits in Arabidopsis has a polygenic basis with overdominance and/or epistasis being presumably the main types of gene action.

Intrasexual competition among females and the stabilization of a conspicuous colour polymorphism in a Lake Victoria cichlid fish

The maintenance of colour polymorphisms within populations has been a long-standing interest in evolutionary ecology. African cichlid fish contain some of the most striking known cases of this phenomenon. Intrasexual selection can be negative frequency dependent when males bias aggression towards phenotypically similar rivals, stabilizing male colour polymorphisms. We propose that where females are territorial and competitive, aggression biases in females may also promote coexistence of female morphs. We studied a polymorphic population of the cichlid fish Neochromis omnicaeruleus from Lake Victoria, in which three distinct female colour morphs coexist: one plain brown and two blotched morphs. Using simulated intruder choice tests in the laboratory, we show that wild-caught females of each morph bias aggression towards females of their own morph, suggesting that females of all three morphs may have an advantage when their morph is locally the least abundant. This mechanism may contribute to the establishment and stabilization of colour polymorphisms. Next, by crossing the morphs, we generated sisters belonging to different colour morphs. We find no sign of aggression bias in these sisters, making pleiotropy unlikely to explain the association between colour and aggression bias in wild fish, which is maintained in the face of gene flow. We conclude that female-female aggression may be one important force for stabilizing colour polymorphism in cichlid fish.

Preference polymorphism for coloration but no speciation in a population of Lake Victoria cichlids

Female mating preference based on male nuptial coloration has been suggested to be an important source of diversifying selection in the radiation of Lake Victoria cichlid fish. Initial variation in female preference is a prerequisite for diversifying selection; however, it is rarely studied in natural populations. In clear water areas of Lake Victoria, the sibling species Pundamilia pundamilia with blue males and Pundamilia nyererei with red males coexist, intermediate phenotypes are rare, and most females have species-assortative mating preferences. Here, we study a population of Pundamilia that inhabits turbid water where male coloration is variable from reddish to blue with most males intermediate. We investigated male phenotype distribution and female mating preferences. Male phenotype was unimodally distributed with a mode on intermediate color in 1 year and more blue-shifted in 2 other years. In mate choice experiments with females of the turbid water population and males from a clearer water population, we found females with a significant and consistent preference for P. pundamilia (blue) males, females with such preferences for P. nyererei (red) males, and many females without a preference. Hence, female mating preferences in this population could cause disruptive selection on male coloration that is probably constrained by the low signal transduction of the turbid water environment. We suggest that if environmental signal transduction was improved and the preference/color polymorphism was stabilized by negative frequency-dependent selection, divergent sexual selection might separate the 2 morphs into reproductively isolated species resembling the clear water species P. pundamilia and P. nyererei.

Influence of ABCB1, ABCC1, ABCC2, and ABCG2 haplotypes on the cellular exposure of nelfinavir in vivo

OBJECTIVES: The human immunodeficiency virus protease inhibitor nelfinavir is substrate of polyspecific drug transporters encoded by ABCB1 (P-glycoprotein), ABCC1 (MRP1) and ABCC2 (MRP2), and an inhibitor of BCRP, encoded by ABCG2. Genetic polymorphism in these genes may be associated with changes in transport function. METHODS: A comprehensive evaluation of single nucleotide polymorphisms (39 SNPs in ABCB1, 7 in ABCC1, 27 in ABCC2, and 16 in ABCG2), and inferred haplotypes was done to assess possible associations of genetic variants with cellular exposure of nelfinavir in vivo. Analysis used peripheral mononuclear cells from individuals receiving nelfinavir (n=28). Key results were re-examined in a larger sample size (n=129) contributing data on plasma drug levels. RESULTS AND CONCLUSIONS: There was no significant association between cellular nelfinavir area under the curve (AUC) and SNPs or haplotypes at ABCC1, ABCC2, ABCG2. There was an association with cellular exposure for two loci in strong linkage disequilibrium: ABCB1 3435C>T; AUCTT>AUCCT>AUCCC (ratio 2.1, 1.4, 1, Ptrend=0.01), and intron 26 +80T>C; AUCCC> AUCCT > AUCTT (ratio 2.4, 1.3, 1, Ptrend=0.006). Haplotypic analysis using tagging SNPs did not improve the single SNP association values.

Influence of CYP2B6 polymorphism on plasma and intracellular concentrations and toxicity of efavirenz and nevirapine in HIV-infected patients

BACKGROUND: Efavirenz (EFV) and nevirapine (NVP) are metabolized by cytochrome P450 2B6 (CYP2B6). Allele 516 G>T (Gln172His) is associated with diminished activity of this isoenzyme, and may lead to differences in drug exposure. METHODS: We evaluated this allele as a pharmacogenetic marker of EFV and NVP pharmacokinetics and EFV toxicity in 167 participants receiving EFV and 59 receiving NVP recruited within the genetics project of the Swiss HIV Cohort Study. Drug concentrations were measured in plasma and in peripheral blood mononuclear cells (PBMCs) from the same sample. Neuropsychological toxicity of EFV (sleep disorders, mood disorders, fatigue) was assessed using a standardized questionnaire. RESULTS AND CONCLUSIONS: CYP2B6 516TT was associated with greater plasma and intracellular exposure to EFV, and greater plasma exposure to NVP. Intracellular drug concentration, and CYP2B6 genotype were predictors of EFV neuropsychological toxicity. CYP2B6 genotyping may be useful to complement an individualization strategy based on plasma drug determinations to increase the safety and tolerability of EFV.