OSIRISv1.2: a named entity recognition system for sequence variants of genes in biomedical literature

Background: Single Nucleotide Polymorphisms, among other type of sequence variants, constitute key elements in genetic epidemiology and pharmacogenomics. While sequence data about genetic variation is found at databases such as dbSNP, clues about the functional and phenotypic consequences of the variations are generally found in biomedical literature. The identification of the relevant documents and the extraction of the information from them are hampered by the large size of literature databases and the lack of widely accepted standard notation for biomedical entities. Thus, automatic systems for the identification of citations of allelic variants of genes in biomedical texts are required. Results: Our group has previously reported the development of OSIRIS, a system aimed at the retrieval of literature about allelic variants of genes http://ibi.imim.es/osirisform.html. Here we describe the development of a new version of OSIRIS (OSIRISv1.2, http://ibi.imim.es/OSIRISv1.2.html webcite) which incorporates a new entity recognition module and is built on top of a local mirror of the MEDLINE collection and HgenetInfoDB: a database that collects data on human gene sequence variations. The new entity recognition module is based on a pattern-based search algorithm for the identification of variation terms in the texts and their mapping to dbSNP identifiers. The performance of OSIRISv1.2 was evaluated on a manually annotated corpus, resulting in 99% precision, 82% recall, and an F-score of 0.89. As an example, the application of the system for collecting literature citations for the allelic variants of genes related to the diseases intracranial aneurysm and breast cancer is presented. Conclusion: OSIRISv1.2 can be used to link literature references to dbSNP database entries with high accuracy, and therefore is suitable for collecting current knowledge on gene sequence variations and supporting the functional annotation of variation databases. The application of OSIRISv1.2 in combination with controlled vocabularies like MeSH provides a way to identify associations of biomedical interest, such as those that relate SNPs with diseases.

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

False-negative PCR result due to gene polymorphism: the example of Neisseria meningitidis.

Early treatment of meningococcal meningitis is mandatory but may negate the cerebrospinal fluid culture. Etiological diagnosis then mainly relies on PCR. Here, we report a case of false-negative results for real-time PCR for a Neisseria meningitidis serogroup B isolate with a polymorphism in the ctrA gene.

Draft genome of the red harvester ant Pogonomyrmex barbatus.

We report the draft genome sequence of the red harvester ant, Pogonomyrmex barbatus. The genome was sequenced using 454 pyrosequencing, and the current assembly and annotation were completed in less than 1 y. Analyses of conserved gene groups (more than 1,200 manually annotated genes to date) suggest a high-quality assembly and annotation comparable to recently sequenced insect genomes using Sanger sequencing. The red harvester ant is a model for studying reproductive division of labor, phenotypic plasticity, and sociogenomics. Although the genome of P. barbatus is similar to other sequenced hymenopterans (Apis mellifera and Nasonia vitripennis) in GC content and compositional organization, and possesses a complete CpG methylation toolkit, its predicted genomic CpG content differs markedly from the other hymenopterans. Gene networks involved in generating key differences between the queen and worker castes (e.g., wings and ovaries) show signatures of increased methylation and suggest that ants and bees may have independently co-opted the same gene regulatory mechanisms for reproductive division of labor. Gene family expansions (e.g., 344 functional odorant receptors) and pseudogene accumulation in chemoreception and P450 genes compared with A. mellifera and N. vitripennis are consistent with major life-history changes during the adaptive radiation of Pogonomyrmex spp., perhaps in parallel with the development of the North American deserts.

Heterozygote advantage and the maintenance of polymorphism for multilocus traits.

Explaining how polymorphism is maintained in the face of selection remains a puzzle since selection tends to erode genetic variation. Provided an infinitely large unsubdivided population and no frequency-dependance of selective values, heterozygote advantage is the text book explanation for the maintenance of polymorphism when selection acts at a diallelic locus. Here, we investigate whether this remains true when selection acts at multiple diallelic loci. We use five different definitions of heterozygote advantage that largely cover this concept for multiple loci. Using extensive numerical simulations, we found no clear associations between the presence of any of the five definitions of heterozygote advantage and the maintenance of polymorphism at all loci. The strength of the association decreases as the number of loci increases or as recombination decreases. We conclude that heterozygote advantage cannot be a general mechanism for the maintenance of genetic polymorphism at multiple loci. These findings suggest that a correlation between the number of heterozygote loci and fitness is not warranted on theoretical ground.

A polymorphism in the leptin gene promoter is associated with anemia in patients with HIV disease.

To study factors associated with anemia and its effect on survival in HIV-infected persons treated with modern combined antiretroviral therapy (cART), we characterized the prevalence of anemia in the Veterans Aging Cohort Study (VACS) and used a candidate gene approach to identify proinflammatory gene single nucleotide polymorphisms (SNPs) associated with anemia in HIV disease. The study comprised 1597 HIV(+) and 865 HIV(-) VACS subjects with DNA, blood, and annotated clinical data available for analysis. Anemia was defined according to World Health Organization criteria (hemoglobin < 13 g/dL and < 12 g/dL in men and women, respectively). The prevalence of anemia in HIV(+) and HIV(-) subjects was 23.1% and 12.9%, respectively. Independent of HIV status, anemia was present in 23.4% and 8% in blacks and whites, respectively. Analysis of our candidate genes revealed that the leptin -2548 G/A SNP was associated with anemia in HIV(+), but not HIV(-), patients, with the AA and AG genotypes significantly predicting anemia (P < .003 and P < .039, respectively, logistic regression). This association was replicated in an independent cohort of HIV(+) women. Our study provides novel insight into the association between genetic variability in the leptin gene and anemia in HIV(+) individuals.

Transposase and cointegrase: specialized transposition proteins of the bacterial insertion sequence IS21 and related elements.

The bacterial insertion sequence IS21 shares with many insertion sequences a two-step, reactive junction transposition pathway, for which a model is presented in this review: a reactive junction with abutted inverted repeats is first formed and subsequently integrated into the target DNA. The reactive junction occurs in IS21-IS21 tandems and IS21 minicircles. In addition, IS21 shows a unique specialization of transposition functions. By alternative translation initiation, the transposase gene codes for two products: the transposase, capable of promoting both steps of the reactive junction pathway, and the cointegrase, which only promotes the integration of reactive junctions but with higher efficiency. This review also includes a survey of the IS21 family and speculates on the possibility that other members present a similar transpositional specialization.

IL28B expression depends on a novel TT/-G polymorphism which improves HCV clearance prediction.

Approximately 3% of the world population is chronically infected with the hepatitis C virus (HCV), with potential development of cirrhosis and hepatocellular carcinoma. Despite the availability of new antiviral agents, treatment remains suboptimal. Genome-wide association studies (GWAS) identified rs12979860, a polymorphism nearby IL28B, as an important predictor of HCV clearance. We report the identification of a novel TT/-G polymorphism in the CpG region upstream of IL28B, which is a better predictor of HCV clearance than rs12979860. By using peripheral blood mononuclear cells (PBMCs) from individuals carrying different allelic combinations of the TT/-G and rs12979860 polymorphisms, we show that induction of IL28B and IFN-γ-inducible protein 10 (IP-10) mRNA relies on TT/-G, but not rs12979860, making TT/-G the only functional variant identified so far. This novel step in understanding the genetic regulation of IL28B may have important implications for clinical practice, as the use of TT/G genotyping instead of rs12979860 would improve patient management.

Evaluation of adding a second marker to overcome Staphylococcus aureus spa typing homoplasies.

The utility of sequencing a second highly variable locus in addition to the spa gene (e.g., double-locus sequence typing [DLST]) was investigated to overcome limitations of a Staphylococcus aureus single-locus typing method. Although adding a second locus seemed to increase discriminatory power, it was not sufficient to definitively infer evolutionary relationships within a single multilocus sequence type (ST-5).

Resistance to nucleoside analog reverse transcriptase inhibitors mediated by human immunodeficiency virus type 1 p6 protein.

Resistance of human immunodeficiency virus type 1 (HIV-1) to antiretroviral agents results from target gene mutation within the pol gene, which encodes the viral protease, reverse transcriptase (RT), and integrase. We speculated that mutations in genes other that the drug target could lead to drug resistance. For this purpose, the p1-p6(gag)-p6(pol) region of HIV-1, placed immediately upstream of pol, was analyzed. This region has the potential to alter Pol through frameshift regulation (p1), through improved packaging of viral enzymes (p6(Gag)), or by changes in activation of the viral protease (p6(Pol)). Duplication of the proline-rich p6(Gag) PTAP motif, necessary for late viral cycle activities, was identified in plasma virus from 47 of 222 (21.2%) patients treated with nucleoside analog RT inhibitor (NRTI) antiretroviral therapy but was identified very rarely from drug-naïve individuals. Molecular clones carrying a 3-amino-acid duplication, APPAPP (transframe duplication SPTSPT in p6(Pol)), displayed a delay in protein maturation; however, they packaged a 34% excess of RT and exhibited a marked competitive growth advantage in the presence of NRTIs. This phenotype is reminiscent of the inoculum effect described in bacteriology, where a larger input, or a greater infectivity of an organism with a wild-type antimicrobial target, leads to escape from drug pressure and a higher MIC in vitro. Though the mechanism by which the PTAP region participates in viral maturation is not known, duplication of this proline-rich motif could improve assembly and packaging at membrane locations, resulting in the observed phenotype of increased infectivity and drug resistance.

Cutting edge: a Toll-like receptor 2 polymorphism that is associated with lepromatous leprosy is unable to mediate mycobacterial signaling.

Toll-like receptors (TLRs) are key mediators of the innate immune response to microbial pathogens. We investigated the role of TLRs in the recognition of Mycobacterium leprae and the significance of TLR2Arg(677)Trp, a recently discovered human polymorphism that is associated with lepromatous leprosy. In mice, TNF-alpha production in response to M. leprae was essentially absent in TLR2-deficient macrophages. Similarly, human TLR2 mediated M. leprae-dependent activation of NF-kappaB in transfected Chinese hamster ovary and human embryonic kidney 293 cells, with enhancement of this signaling in the presence of CD14. In contrast, activation of NF-kappaB by human TLR2Arg(677)Trp was abolished in response to M. leprae and Mycobacterium tuberculosis. The impaired function of this TLR2 variant provides a molecular mechanism for the poor cellular immune response associated with lepromatous leprosy and may have important implications for understanding the pathogenesis of other mycobacterial infections.

Covariation between plumage colour polymorphism and diet in the Barn Owl Tyto alba

Several hypotheses might explain the evolution and maintenance of colour morphs within animal populations. The 'alternative foraging strategy' hypothesis states that alternative colour morphs exploit different ecological niches. This hypothesis predicts that morphs differ in diet, either because foraging success on alternative prey species is morph-dependent or because differently coloured individuals exploit alternative habitats. I examined this prediction in the Barn Owl Tyto alba, a bird that varies in plumage coloration continuously from dark reddish-brown to white. On the European continent, Owls are light-coloured (subspecies T. a. alba) in the south and reddish-brown (T. a. guttata) in the north; in central Europe the two subspecies interbreed, generating many colour variants. If plumage coloration indicates alternative foraging strategies, in sympatry dark- and light-coloured owls should consume prey species that are typical of the diets of T. a. guttata and T. a. alba in allopatry, respectively. In line with this prediction, both in allopatry and in sympatry in Switzerland T. a. guttata fed primarily upon Common Voles Microtus arvalis and T. a. alba upon Wood Mice Apodemus spp. Statistical analyses suggest that morph-dependent diet did not arise from a non-random habitat distribution of owls with respect to plumage coloration. This suggests that foraging success upon alternative prey is morph-dependent.

A large deletion in the adRP gene PRPF31: evidence that haploinsufficiency is the cause of disease.

PURPOSE: To report a large deletion that encompasses more than 90% of PRPF31 gene and two other neighboring genes in their entirety in an adRP pedigree that appears to show only the typical clinical features of retinitis pigmentosa. METHODS: To identify PRPF31 mutation in a dominant RP family (ADRP2) previously linked to the RP11 locus, the 14 exons of PRPF31 were screened for mutations by direct sequencing. To investigate the possibility of a large deletion, microsatellite markers near PRPF31 gene were analyzed by non-denaturing PAGE. RESULTS: Initial screening of PRPF31 gene in the ADRP2 family did not reveal an obvious mutation. A large deletion was however suspected due to lack of heterozygosity for nearly all PRPF31 intragenic single nucleotide polymorphysm (SNPs). In order to estimate the size of the deletion, SNPs and microsatellite markers spanning and flanking PRPF31 were analyzed in the entire ADRP2 family. Haplotype analysis with the above markers suggested a deletion of approximately 30 kb that included the putative promoter region of a novel gene OSCAR, the entire genomic content of genes NDUFA3, TFPT and more than 90% of PRPF31 gene. Sequence analysis of the region flanking the potential deletion showed a high presence of Alu elements implicating Alu mediated recombination as the mechanism responsible for this event. CONCLUSIONS: This mutation provides evidence that haploinsufficiency rather than aberrant function of mutated proteins is the cause of disease in these adRP patients with mutations in PRPF31 gene.

Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.