Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data.

MOTIVATION: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts. Results: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays.

Simultaneous identification of multiple mammalian species from mixed forensic samples based on mtDNA control region length polymorphism

Molecular species identification in mixed or contaminated biological material has always been problematic. We developed a simple and accurate method for mammal DNA identification in mixtures, based on interspecific mitochondrial DNA control region length polymorphism. Contrary to other published methods dealing with species mixtures, our protocol requires a single universal primer pair and amplification step, and is not based on a pre-defined panel of species. This protocol has been routinely employed by our laboratory for species identification in dozens of human and animal forensic caseworks. Six representative forensic caseworks involving the specific identification of mixed animal samples are reported in this paper, in order to demonstrate the applicability and usefulness of the method.

Dramatic influence of V beta gene polymorphism on an antigen-specific CD8+ T cell response in vivo.

According to recent crystallographic studies, the TCR-alpha beta contacts MHC class I-bound antigenic peptides via the polymorphic V gene-encoded complementarity-determining region 1 beta (CDR1 beta) and the hypervariable (D)J-encoded CDR3 beta and CDR3 alpha domains. To evaluate directly the relative importance of CDR1 beta polymorphism on the fine specificity of T cell responses in vivo, we have taken advantage of congenic V beta a and V beta b mouse strains that differ by a CDR1 polymorphism in the V beta 10 gene segment. The V beta 10-restricted CD8+ T cell response to a defined immunodominant epitope was dramatically reduced in V beta a compared with V beta b mice, as measured either by the expansion of V beta 10+ cells or by the binding of MHC-peptide tetramers. These data indicate that V beta polymorphism has an important impact on TCR-ligand binding in vivo, presumably by modifying the affinity of CDR1 beta-peptide interactions.

Extra-pair paternity, testes size and testosterone level in relation to colour polymorphism in the barn owl Tyto alba

In many bird populations, individuals display one of several genetically inherited colour morphs. Colour polymorphism can be maintained by several mechanisms one of which being frequency-dependent selection with colour morphs signalling alternative mating strategies. One morph may be dominant and territorial, and another one adopt a sneaky behaviour to gain access to fertile females. We tested this hypothesis in the barn owl Tyto alba in which coloration varies from reddish-brown to white. This trait is heritable and neither sensitive to the environment in which individuals live nor to body condition. In Switzerland, reddish-brown males were observed to feed their brood at a higher rate and to produce more offspring than white males. This observation lead us to hypothesize that white males may equalise fitness by investing more effort in extra-pair copulations. This hypothesis predicts that lighter Coloured males produce more extra-pair young, have larger testes and higher levels of circulating testosterone. However, our results are not consistent with these three predictions. First, paternity analyses of 54 broods with a total of 211 offspring revealed that only one young was not sired by the male that was feeding it. Second, testes size was not correlated with male plumage coloration suggesting that white males are not sexually more active. Finally, in nestlings at the time of feather growth testosterone level was not related to plumage coloration suggesting that this androgen is not required for the expression of this plumage trait. Our study therefore indicates that in the barn owl colour polymorphism plays no role in the probability of producing extra-pair young.

SSCP-polymorphism in intron 12 of the CFTR gene recognized by BclI

Source/Description: SSCP analysis of intron 12 of the CFTR gene from PCR products showed an extra band in several DNA samples. Sequencing of the additional fragment extra band revealed a T- A change in the position 1898 + 152 of CFTR (Fig. 1). The change is a polymorphism which can be identified by SSCP or by BclI digestion...

Mspl restriction fragment length polymorphism near exon 10 of cystic fibrosis (CFTR) gene.

Source/Description: The probe used is a 98 bp fragment amplified by PCR from a cDNA clone of the CFTR gene or from genomic DNA corresponding to exon 10, using two primers from this exon (1)...

Malaria haplotype frequency estimation.

We present a Bayesian approach for estimating the relative frequencies of multi-single nucleotide polymorphism (SNP) haplotypes in populations of the malaria parasite Plasmodium falciparum by using microarray SNP data from human blood samples. Each sample comes from a malaria patient and contains one or several parasite clones that may genetically differ. Samples containing multiple parasite clones with different genetic markers pose a special challenge. The situation is comparable with a polyploid organism. The data from each blood sample indicates whether the parasites in the blood carry a mutant or a wildtype allele at various selected genomic positions. If both mutant and wildtype alleles are detected at a given position in a multiply infected sample, the data indicates the presence of both alleles, but the ratio is unknown. Thus, the data only partially reveals which specific combinations of genetic markers (i.e. haplotypes across the examined SNPs) occur in distinct parasite clones. In addition, SNP data may contain errors at non-negligible rates. We use a multinomial mixture model with partially missing observations to represent this data and a Markov chain Monte Carlo method to estimate the haplotype frequencies in a population. Our approach addresses both challenges, multiple infections and data errors.

Current status of a model system: the gene Gp-9 and its association with social organization in fire ants.

The Gp-9 gene in fire ants represents an important model system for studying the evolution of social organization in insects as well as a rich source of information relevant to other major evolutionary topics. An important feature of this system is that polymorphism in social organization is completely associated with allelic variation at Gp-9, such that single-queen colonies (monogyne form) include only inhabitants bearing B-like alleles while multiple-queen colonies (polygyne form) additionally include inhabitants bearing b-like alleles. A recent study of this system by Leal and Ishida (2008) made two major claims, the validity and significance of which we examine here. After reviewing existing literature, analyzing the methods and results of Leal and Ishida (2008), and generating new data from one of their study sites, we conclude that their claim that polygyny can occur in Solenopsis invicta in the U.S.A. in the absence of expression of the b-like allele Gp-9(b) is unfounded. Moreover, we argue that available information on insect OBPs (the family of proteins to which GP-9 belongs), on the evolutionary/population genetics of Gp-9, and on pheromonal/behavioral control of fire ant colony queen number fails to support their view that GP-9 plays no role in the chemosensory-mediated communication that underpins regulation of social organization. Our analyses lead us to conclude that there are no new reasons to question the existing consensus view of the Gp-9 system outlined in Gotzek and Ross (2007).

Geographical and altitudinal population genetic structure of two dung fly species with contrasting mobility and temperature preference.

Local adaptation of populations requires some degree of spatio-temporal isolation. Previous studies of the two dung fly species Scathophaga stercoraria and Sepsis cynipsea have revealed low levels of geographic and altitudinal genetic differentiation in quantitative life history and morphological traits, but instead high degrees of phenotypic plasticity. These patterns suggest that gene flow is extensive despite considerable geographic barriers and large spatio-temporal variation in selection on body size and related traits. In this study we addressed this hypothesis by investigating genetic differentiation of dung fly populations throughout Switzerland based on the same 10 electrophoretic loci in each species. Overall, we found no significant geographic differentiation of populations for either species. This is inconsistent with the higher rates of gene flow expected due to better flying capacity of the larger S. stercoraria. However, heterozygote deficiencies within populations indicated structuring on a finer scale, seen for several loci in S. cynipsea, and for the locus PGM (Phosphoglucomutase) in S. stercoraria. Additionally, S. cynipsea showed a tendency towards a greater gene diversity at higher altitudes, mediated primarily by the locus MDH (malate dehydrogenase), at which a second allele was only present in populations above 1000 m. This may be caused by increased environmental stress at higher altitudes in this warm-adapted species. MDH might thus be a candidate locus subject to thermal selection in this species, but this remains to be corroborated by direct evidence. In S. stercoraria, no altitudinal variation was found.

Genetic regulation of colony social organization in fire ants: an integrative overview

Expression of colony social organization in fire ants appears to be under the control of a single Mendelian factor of large effect. Variation in colony queen number in Solenopsis invicta and its relatives is associated with allelic variation at the gene Gp-9, but not with variation at other unlinked genes; workers regulate queen identity and number on the basis of Gp-9 genotypic compatibility. Nongeneticfactors, such as prior social experience, queen reproductive status, and local environment, have negligible effects on queen number which illustrates the nearly complete penetrance of Gp-9. As predicted, queen number can be manipulated experimentally by altering worker Gp-9 genotype frequencies. The Gp-9 allele lineage associated with polygyny in South American fire? ants has been retained across multiple speciation events, which may signal the action of balancing selection to maintain social polymorphism in these species. Moreover positive selection is implicated in driving the molecular evolution of Gp-9 in association with the origin of polygyny. The identity of the product of Gp-9 as an odorant-binding protein suggests plausible scenarios for its direct involvement in the regulation of queen number via a role in chemical communication. While these and other lines of evidence show that Gp-9 represents a legitimate candidate gene of major effect, studies aimed at determining (i) the biochemical pathways in which GP-9 functions; (ii) the phenotypic effects of molecular variation at Gp-9 and other pathway genes; and (iii) the potential involvement of genes in linkage disequilibrium with Gp-9 are needed to elucidate the genetic architecture underlying social organization in fire ants. Information that reveals the links between molecular variation, individual phenotype, and colony-level behaviors, combined with behavioral models that incorporate details of the chemical communication involved in regulating queen number will yield a novel integrated view of the evolutionary changes underlying a key social adaptation.

Genomewide association study of atazanavir pharmacokinetics and hyperbilirubinemia in AIDS Clinical Trials Group protocol A5202.

BACKGROUND: Atazanavir-associated hyperbilirubinemia can cause premature discontinuation of atazanavir and avoidance of its initial prescription. We used genomewide genotyping and clinical data to characterize determinants of atazanavir pharmacokinetics and hyperbilirubinemia in AIDS Clinical Trials Group protocol A5202. METHODS: Plasma atazanavir pharmacokinetics and indirect bilirubin concentrations were characterized in HIV-1-infected patients randomized to atazanavir/ritonavir-containing regimens. A subset had genomewide genotype data available. RESULTS: Genomewide assay data were available from 542 participants, of whom 475 also had data on estimated atazanavir clearance and relevant covariates available. Peak bilirubin concentration and relevant covariates were available for 443 participants. By multivariate analysis, higher peak on-treatment bilirubin levels were found to be associated with the UGT1A1 rs887829 T allele (P=6.4×10), higher baseline hemoglobin levels (P=4.9×10), higher baseline bilirubin levels (P=6.7×10), and slower plasma atazanavir clearance (P=8.6×10). For peak bilirubin levels greater than 3.0 mg/dl, the positive predictive value of a baseline bilirubin level of 0.5 mg/dl or higher with hemoglobin concentrations of 14 g/dl or higher was 0.51, which increased to 0.85 with rs887829 TT homozygosity. For peak bilirubin levels of 3.0 mg/dl or lower, the positive predictive value of a baseline bilirubin level less than 0.5 mg/dl with a hemoglobin concentration less than 14 g/dl was 0.91, which increased to 0.96 with rs887829 CC homozygosity. No polymorphism predicted atazanavir pharmacokinetics at genomewide significance. CONCLUSION: Atazanavir-associated hyperbilirubinemia is best predicted by considering UGT1A1 genotype, baseline bilirubin level, and baseline hemoglobin level in combination. Use of ritonavir as a pharmacokinetic enhancer may have abrogated genetic associations with atazanavir pharmacokinetics.

DnaSP, DNA polymorphism analyses by the coalescent and other methods

DnaSP is a software package for the analysis of DNA polymorphism data. Present version introduces several new modules and features which, among other options allow: (1) handling big data sets (~5 Mb per sequence); (2) conducting a large number of coalescent-based tests by Monte Carlo computer simulations; (3) extensive analyses of the genetic differentiation and gene flow among populations; (4) analysing the evolutionary pattern of preferred and unpreferred codons; (5) generating graphical outputs for an easy visualization of results. Availability: The software package, including complete documentation and examples, is freely available to academic users from: http://www.ub.es/dnasp