Bacillus anthracis: Molecular taxonomy, population genetics, phylogeny and patho-evolution

Bacillus anthracis, the etiological agent of anthrax, manifests a particular bimodal lifestyle. This bacterial species alternates between short replication phases of 20-40 generations that strictly require infection of the host, normally causing death, interrupted by relatively long, mostly dormant phases as spores in the environment. Hence, the B. anthracis genome is highly homogeneous. This feature and the fact that strains from nearly all parts of the world have been analysed for canonical single nucleotide polymorphisms (canSNPs) and variable number tandem repeats (VNTRs) has allowed the development of molecular epidemiological and molecular clock models to estimate the age of major diversifications in the evolution of B. anthracis and to trace the global spread of this pathogen, which was mostly promoted by movement of domestic cattle with settlers and by international trade of contaminated animal products. From a taxonomic and phylogenetic point of view, B. anthracis is a member of the Bacillus cereus group. The differentiation of B. anthracis from B. cereus sensu strict, solely based on chromosomal markers, is difficult. However, differences in pathogenicity clearly differentiate B. anthracis from B. cereus and are marked by the strict presence of virulence genes located on the two virulence plasmids pXO1 and pXO2, which both are required by the bacterium to cause anthrax. Conversely, anthrax-like symptoms can also be caused by organisms with chromosomal features that are more closely related to B. cereus, but which carry these virulence genes on two plasmids that largely resemble the B. anthracis virulence plasmids. (C) 2011 Elsevier B.V. All rights reserved.

Novel ADAMTS13 mutations in an obstetric patient with upshaw-schulman syndrome

Upshaw-Schulman syndrome (USS) is a rarely reported congenital form of thrombotic thrombocytopenic purpura (TTP) that results from mutations in the ADAMTS13 gene. Many USS patients are diagnosed during the second or third trimester of their first pregnancy. We present a patient diagnosed with USS following retinal detachments and intrauterine fetal demise at 34 weeks of gestation. The patient's plasma was tested for ADAMTS13 activity, inhibitor, and antibody. Subsequently, she and her first-degree relatives had ADAMTS13 gene sequencing. Initially, the patient was found to have an ADAMTS13 activity of <5% in the absence of an ADAMTS13 inhibitor (FRETS assay) or antibody (immunoassay). Repeat studies in the months following hospital discharge showed persistent, undetectable ADAMTS13 activity and she was given a clinical diagnosis of USS. Molecular sequencing demonstrated two novel missense mutations in the ADAMTS13 gene: one in the maternal exon 17 (p.Ala690Thr due to nucleotide substitution c.2068 G>A) and another in the paternal exon 22 (p.Arg915Cys due to nucleotide substitution c.2746 C>T). In addition to being compound heterozygous for two ADAMTS13 mutations, the patient also had two maternally inherited single nucleotide polymorphisms: p.P618A (exon 16) and p.A732V (exon 18). Her parents and only sister had normal or near-normal ADAMTS13 activity. Each was heterozygous for one of the novel missense mutations. This case highlights the importance of molecular analysis of the ADAMTS13 gene in patients and family members when the severe ADAMTS13 deficiency does not appear to be autoimmune in nature. J. Clin. Apheresis, 2012. © 2012 Wiley Periodicals, Inc.

Molecular characterization of five porcine candidate genes for drip loss in pork

Drip loss is the loss of fluid from a piece of meat without mechanical force and represents an important meat quality trait. Previous work revealed a quantitative trait locus (QTL) for drip loss in pork in an experimental Duroc x Pietrain (DUPI) F2 family on SSC 5. Based on functional data indicating their possible involvement in water holding capacity and their expression in skeletal muscle, we selected five positional candidates (ACO2, ADSL, CBY1, KCNJ4, PLA2AG6) out of 130 predicted genes in the QTL interval for further analysis. We performed a mutation analysis of all coding exons and discovered 204 polymorphisms. We genotyped 39 single nucleotide polymorphisms (SNPs) in 192 Pietrain pigs with extreme drip loss phenotypes and detected a possible association with drip loss for one non-coding SNP in the ADSL gene (ss107793818, p(raw) = 0.021). Correspondingly, ADSL diplotypes were associated with drip loss and pH1 of M. longissimus dorsi. However, after correction for multiple testing, none of the tested SNPs were significantly associated with drip loss. One possible explanation for these results is that one of the QTL-alleles from the experimental DUPI family may be fixed or nearly fixed in the tested Pietrain population.

Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster

Cytochrome P450 enzymes (CYP450s) represent a superfamily of haem-thiolate proteins. CYP450s are most abundant in the liver, a major site of drug metabolism, and play key roles in the metabolism of a variety of substrates, including drugs and environmental contaminants. Interaction of two or more different drugs with the same enzyme can account for adverse effects and failure of therapy. Human CYP3A4 metabolizes about 50% of all known drugs, but little is known about the orthologous CYP450s in horses. We report here the genomic organization of the equine CYP3A gene cluster as well as a comparative analysis with the human CYP3A gene cluster. The equine CYP450 genes of the 3A family are located on ECA 13 between 6.97-7.53 Mb, in a region syntenic to HSA 7 99.05-99.35 Mb. Seven potential, closely linked equine CYP3A genes were found, in contrast to only four genes in the human genome. RNA was isolated from an equine liver sample, and the approximately 1.5-kb coding sequence of six CYP3A genes could be amplified by RT-PCR. Sequencing of the RT-PCR products revealed numerous hitherto unknown single nucleotide polymorphisms (SNPs) in these six CYP3A genes, and one 6-bp deletion compared to the reference sequence (EquCab2.0). The presence of the variants was confirmed in a sample of genomic DNA from the same horse. In conclusion, orthologous genes for the CYP3A family exist in horses, but their number differs from those of the human CYP3A gene family. CYP450 genes of the same family show high homology within and between mammalian species, but can be highly polymorphic.

Design, conduct, and analysis of a multicenter, pharmacogenomic, biomarker study in matched patients with severe sepsis treated with or without drotrecogin Alfa (activated)

A genomic biomarker identifying patients likely to benefit from drotrecogin alfa (activated) (DAA) may be clinically useful as a companion diagnostic. This trial was designed to validate biomarkers (improved response polymorphisms (IRPs)). Each IRP (A and B) contains two single nucleotide polymorphisms that were associated with a differential DAA treatment effect.

First insights into the phylogenetic diversity of Mycobacterium tuberculosis in Nepal

Background Tuberculosis (TB) is a major public health problem in Nepal. Strain variation in Mycobacterium tuberculosis may influence the outcome of TB infection and disease. To date, the phylogenetic diversity of M. tuberculosis in Nepal is unknown. Methods and Findings We analyzed 261 M. tuberculosis isolates recovered from pulmonary TB patients recruited between August 2009 and August 2010 in Nepal. M. tuberculosis lineages were determined by single nucleotide polymorphisms (SNP) typing and spoligotyping. Drug resistance was determined by sequencing the hot spot regions of the relevant target genes. Overall, 164 (62.8%) TB patients were new, and 97 (37.2%) were previously treated. Any drug resistance was detected in 50 (19.2%) isolates, and 16 (6.1%) were multidrug-resistant. The most frequent M. tuberculosis lineage was Lineage 3 (CAS/Delhi) with 106 isolates (40.6%), followed by Lineage 2 (East-Asian lineage, includes Beijing genotype) with 84 isolates (32.2%), Lineage 4 (Euro-American lineage) with 41 (15.7%) isolates, and Lineage 1 (Indo-Oceanic lineage) with 30 isolates (11.5%). Based on spoligotyping, we found 45 different spoligotyping patterns that were previously described. The Beijing (83 isolates, 31.8%) and CAS spoligotype (52, 19.9%) were the dominant spoligotypes. A total of 36 (13.8%) isolates could not be assigned to any known spoligotyping pattern. Lineage 2 was associated with female sex (adjusted odds ratio [aOR] 2.58, 95% confidence interval [95% CI] 1.42–4.67, p = 0.002), and any drug resistance (aOR 2.79; 95% CI 1.43–5.45; p = 0.002). We found no evidence for an association of Lineage 2 with age or BCG vaccination status. Conclusions We found a large genetic diversity of M. tuberculosis in Nepal with representation of all four major lineages. Lineages 3 and 2 were dominating. Lineage 2 was associated with clinical characteristics. This study fills an important gap on the map of the M. tuberculosis genetic diversity in the Asian region.

Two new rapid SNP-typing methods for classifying Mycobacterium tuberculosis complex into the main phylogenetic lineages

There is increasing evidence that strain variation in Mycobacterium tuberculosis complex (MTBC) might influence the outcome of tuberculosis infection and disease. To assess genotype-phenotype associations, phylogenetically robust molecular markers and appropriate genotyping tools are required. Most current genotyping methods for MTBC are based on mobile or repetitive DNA elements. Because these elements are prone to convergent evolution, the corresponding genotyping techniques are suboptimal for phylogenetic studies and strain classification. By contrast, single nucleotide polymorphisms (SNP) are ideal markers for classifying MTBC into phylogenetic lineages, as they exhibit very low degrees of homoplasy. In this study, we developed two complementary SNP-based genotyping methods to classify strains into the six main human-associated lineages of MTBC, the "Beijing" sublineage, and the clade comprising Mycobacterium bovis and Mycobacterium caprae. Phylogenetically informative SNPs were obtained from 22 MTBC whole-genome sequences. The first assay, referred to as MOL-PCR, is a ligation-dependent PCR with signal detection by fluorescent microspheres and a Luminex flow cytometer, which simultaneously interrogates eight SNPs. The second assay is based on six individual TaqMan real-time PCR assays for singleplex SNP-typing. We compared MOL-PCR and TaqMan results in two panels of clinical MTBC isolates. Both methods agreed fully when assigning 36 well-characterized strains into the main phylogenetic lineages. The sensitivity in allele-calling was 98.6% and 98.8% for MOL-PCR and TaqMan, respectively. Typing of an additional panel of 78 unknown clinical isolates revealed 99.2% and 100% sensitivity in allele-calling, respectively, and 100% agreement in lineage assignment between both methods. While MOL-PCR and TaqMan are both highly sensitive and specific, MOL-PCR is ideal for classification of isolates with no previous information, whereas TaqMan is faster for confirmation. Furthermore, both methods are rapid, flexible and comparably inexpensive.

Pediatric and adult eosinophilic esophagitis: similarities and differences

Early in the 1990s, several case series described adults suffering from dysphagia and children with refractory reflux symptoms, both accompanied by an eosinophil-predominant infiltration, thereby conclusively distinguishing it from gastroesophageal reflux disease. Eosinophilic esophagitis (EoE) was recognized as its own entity in the adult and in the pediatric literature. In the last decade, evidence has accumulated that EoE represents a T-helper (Th)2-type inflammatory disease. Remodeling of the esophagus is a hallmark of EoE, leading to esophageal dysfunction and bolus impaction. Familial occurrence and disease association with single-nucleotide polymorphisms underscore the influence of genetics in this disease. Eosinophilic esophagitis may affect individuals at any age, although the clinical presentation is highly age dependent. There is a significant allergic bias in the EoE population, with the majority of patients having concurrent allergic rhinitis, asthma, eczema, and/or a history of atopy. One noteworthy difference is that in children, EoE seems to be primarily a food antigen-driven disease, whereas in adults, mainly aeroallergen sensitization has been observed. Treatment modalities for EoE include the 3Ds: drugs, diet, and dilation. The crucial question of whether adult and pediatric EoE are different phenotypes of one single entity or whether we are confronted with two different diseases is still open. Here, we review similarities and differences between EoE in adults and children.

Common obesity risk alleles in childhood attention-deficit/hyperactivity disorder

Children with attention-deficit/hyperactivity disorder (ADHD) have a higher rate of obesity than children without ADHD. Obesity risk alleles may overlap with those relevant for ADHD. We examined whether risk alleles for an increased body mass index (BMI) are associated with ADHD and related quantitative traits (inattention and hyperactivity/impulsivity). We screened 32 obesity risk alleles of single nucleotide polymorphisms (SNPs) in a genome-wide association study (GWAS) for ADHD based on 495 patients and 1,300 population-based controls and performed in silico analyses of the SNPs in an ADHD meta-analysis comprising 2,064 trios, 896 independent cases, and 2,455 controls. In the German sample rs206936 in the NUDT3 gene (nudix; nucleoside diphosphate linked moiety X-type motif 3) was associated with ADHD risk (OR: 1.39; P = 3.4 × 10(-4) ; Pcorr  = 0.01). In the meta-analysis data we found rs6497416 in the intronic region of the GPRC5B gene (G protein-coupled receptor, family C, group 5, member B; P = 7.2 × 10(-4) ; Pcorr  = 0.02) as a risk allele for ADHD. GPRC5B belongs to the metabotropic glutamate receptor family, which has been implicated in the etiology of ADHD. In the German sample rs206936 (NUDT3) and rs10938397 in the glucosamine-6-phosphate deaminase 2 gene (GNPDA2) were associated with inattention, whereas markers in the mitogen-activated protein kinase 5 gene (MAP2K5) and in the cell adhesion molecule 2 gene (CADM2) were associated with hyperactivity. In the meta-analysis data, MAP2K5 was associated with inattention, GPRC5B with hyperactivity/impulsivity and inattention and CADM2 with hyperactivity/impulsivity. Our results justify further research on the elucidation of the common genetic background of ADHD and obesity.

Examining The Genetic, Social-Cognitive, And Neural Correlates Of Autism Spectrum Disorder Behaviors In Typically Developing Chi

Autism spectrum disorders (ASD) are pervasive developmental disorders that affect approximately 1 in 50 children (Blumberg et al., 2013). Due to the social nature of the deficits that characterize the disorders, many have classified them as disorders of social cognition, which is the process that individuals use in order to successfully interact with members of their own species (Frith & Frith, 2007). Previous research has typically neglected the spectrum nature of ASD in favor of a more categorical approach of ¿autistic¿ versus ¿non-autistic,¿ but the spectrum requires a more continuous approach. Thus, the present study sought to examine the genetic, social-cognitive, and neural correlates of ASD-like traits as well as the relationship between these dimensions in typically developing children. Parents and children completed several quantitative measures examining several areas of social-cognitive functioning, including theory of mind and social functioning, restricted/repetitive behaviors and interests, and adaptive/maladaptive functioning. Children were also asked to undergo an EEG and both parents and children contributed a saliva sample that was used to sequence four single nucleotide polymorphisms (SNPs) of the OXTR gene, rs1042778, rs53576, rs2254298, and rs237897. We successfully demonstrated a significant relationship between behavioral measures of social-cognition and differences in face perception via the N170. However, the directionality of these relationships varied based on the behavioral measure and particular N170 difference scores. We also found support for the associations between the G_G allelic combination of rs1042778 and the A_A and A_G allelic combinations of rs2254298 and increased ASD-like behavior with decreased social-cognitive functioning. In contrast, our results contradict previous findings with rs237897 and imply that individuals with the A_A and A_G genotypes are less similar to those with ASD and have higher social cognitive functioning than those with the G_G genotype. In conclusion, we have demonstrated the existence of ASD-like traits in typically developing children and have shown a link between behavioral, genetic, and neural correlates of social-cognition. These findings demonstrate the importance of considering autism as a spectrum disorder and provide support for the move to a more continuous approach to neurodevelopmental disorders.

β-Arrestin2 influences the response to methadone in opioid-dependent patients

β-Arrestin2 (ARRB2) is a component of the G-protein-coupled receptor complex and is involved in μ-opioid and dopamine D(2) receptor signaling, two central processes in methadone signal transduction. We analyzed 238 patients in methadone maintenance treatment (MMT) and identified a haplotype block (rs34230287, rs3786047, rs1045280 and rs2036657) spanning almost the entire ARRB2 locus. Although none of these single nucleotide polymorphisms (SNPs) leads to a change in amino-acid sequence, we found that for all the SNPs analyzed, with exception of rs34230287, homozygosity for the variant allele confers a nonresponding phenotype (n=73; rs1045280C and rs2036657G: OR=3.1, 95% CI=1.5-6.3, P=0.004; rs3786047A: OR=2.5, 95% CI=1.2-5.1, P=0.02) also illustrated by a 12-fold shorter period of negative urine screening (P=0.01). The ARRB2 genotype may thus contribute to the interindividual variability in the response to MMT and help to predict response to treatment.

Identification of SNPs in the cystic fibrosis interactome influencing pulmonary progression in cystic fibrosis

There is growing evidence that the great phenotypic variability in patients with cystic fibrosis (CF) not only depends on the genotype, but apart from a combination of environmental and stochastic factors predominantly also on modifier gene effects. It has been proposed that genes interacting with CF transmembrane conductance regulator (CFTR) and epithelial sodium channel (ENaC) are potential modifiers. Therefore, we assessed the impact of single-nucleotide polymorphisms (SNPs) of several of these interacters on CF disease outcome. SNPs that potentially alter gene function were genotyped in 95 well-characterized p.Phe508del homozygous CF patients. Linear mixed-effect model analysis was used to assess the relationship between sequence variants and the repeated measurements of lung function parameters. In total, we genotyped 72 SNPs in 10 genes. Twenty-five SNPs were used for statistical analysis, where we found strong associations for one SNP in PPP2R4 with the lung clearance index (P ≤ 0.01), the specific effective airway resistance (P ≤ 0.005) and the forced expiratory volume in 1 s (P ≤ 0.005). In addition, we identified one SNP in SNAP23 to be significantly associated with three lung function parameters as well as one SNP in PPP2R1A and three in KRT19 to show a significant influence on one lung function parameter each. Our findings indicate that direct interacters with CFTR, such as SNAP23, PPP2R4 and PPP2R1A, may modify the residual function of p.Phe508del-CFTR while variants in KRT19 may modulate the amount of p.Phe508del-CFTR at the apical membrane and consequently modify CF disease.

Molecular characterization of the equine collagen, type IX, alpha 2 (COL9A2) gene on horse chromosome 2p16-->p15

The mammalian collagen, type IX, alpha 2 gene (COL9A2) encodes the alpha-2 chain of type IX collagen and is located on horse chromosome 2p16-->p14 harbouring a quantitative trait locus for osteochondrosis. We isolated a bacterial artificial chromosome (BAC) clone containing the equine COL9A2 gene and determined the complete genomic sequence of this gene. Cloning and characterization of equine COL9A2 revealed that the equine gene consists of 32 exons spanning approximately 15 kb. The COL9A2 transcript encodes a single protein of 688 amino acids. Thirty two single nucleotide polymorphisms (SNPs) equally distributed in the gene were detected in a mutation scan of eight unrelated Hanoverian warmblood stallions, including one SNP that affects the amino acid sequence of COL9A2. Comparative analyses between horse, human, mouse and rat indicate that the chromosomal location of equine COL9A2 is in agreement with known chromosomal synteny relationships. The comparison of the gene structure and transcript revealed a high degree of conservation towards the other mammalian COL9A2 genes. We chose three informative SNPs for association and linkage disequilibrium tests in three to five paternal half-sib families of Hanoverian warmblood horses consisting of 44 to 75 genotyped animals. The test statistics did not reach the significance threshold of 5% and so we could not show an association of COL9A2 with equine osteochondrosis.

Rapid and reliable genotyping of polymorphic loci modifying correct splicing of CFTR pre-mRNA using mass spectrometry

We describe a fast and unambiguous method for haplotyping the (TG)mTn repeat in IVS8 and determining three other single nucleotide polymorphisms (SNPs) in exons 10, 14a and 24 in the cystic fibrosis transmembrane conductance regulator (CFTR) gene affecting correct splicing of the CFTR pre-mRNA using primer extension and mass spectrometry. The diagnostic products are generated by primer extension (PEX) reactions, which require a single detection primer complementary to a region downstream of a target strand's variable site. On addition of a polymerase and an appropriate mixture of dNTP's and 2', 3'-dideoxynucleotide triphosphates (ddNTP's), the primer is extended through the mutation region until the first ddNTP is incorporated and the mass of the extension products determines the composition of the variable site. Analysis of patient DNA assigned the correct and unambiguous haplotype for the (TG)mTn repeat in intron 8 of the CFTR gene. Additional crucial SNPs influencing correct splicing in exon 10, 14 and 24 can easily be detected by biplexing the assay to genotype allelic variants important for correct splicing of the CFTR pre-mRNA. Different PEX reactions with subsequent mass spectrometry generate sufficient data, to enable unambiguous and easy haplotyping of the (TG)mTn repeat in the CFTR gene. The method can be easily extended to the inclusion of additional SNPs of interest by biplexing some of the PEX reactions. All experimental steps required for PEX are amenable to the high degree of automation desirable for a high-throughput diagnostic setting, facilitating the work of clinicians involved in the diagnosis of non-classic cystic fibrosis.

ASSOCIATON TESTS THAT ACCOMMODATE GENOTYPING ERRORS

High-throughput SNP arrays provide estimates of genotypes for up to one million loci, often used in genome-wide association studies. While these estimates are typically very accurate, genotyping errors do occur, which can influence in particular the most extreme test statistics and p-values. Estimates for the genotype uncertainties are also available, although typically ignored. In this manuscript, we develop a framework to incorporate these genotype uncertainties in case-control studies for any genetic model. We verify that using the assumption of a “local alternative” in the score test is very reasonable for effect sizes typically seen in SNP association studies, and show that the power of the score test is simply a function of the correlation of the genotype probabilities with the true genotypes. We demonstrate that the power to detect a true association can be substantially increased for difficult to call genotypes, resulting in improved inference in association studies.