Development and validation of a real-time quantitative PCR assay for rapid identification of Bacillus anthracis in environmental samples

A real-time polymerase chain reaction (PCR) assay was developed for rapid identification of Bacillus anthracis in environmental samples. These samples often harbor Bacillus cereus bacteria closely related to B. anthracis, which may hinder its specific identification by resulting in false positive signals. The assay consists of two duplex real-time PCR: the first PCR allows amplification of a sequence specific of the B. cereus group (B. anthracis, B. cereus, Bacillus thuringiensis, Bacillus weihenstephanensis, Bacillus pseudomycoides, and Bacillus mycoides) within the phosphoenolpyruvate/sugar phosphotransferase system I gene and a B. anthracis specific single nucleotide polymorphism within the adenylosuccinate synthetase gene. The second real-time PCR assay targets the lethal factor gene from virulence plasmid pXO1 and the capsule synthesis gene from virulence plasmid pXO2. Specificity of the assay is enhanced by the use of minor groove binding probes and/or locked nucleic acids probes. The assay was validated on 304 bacterial strains including 37 B. anthracis, 67 B. cereus group, 54 strains of non-cereus group Bacillus, and 146 Gram-positive and Gram-negative bacteria strains. The assay was performed on various environmental samples spiked with B. anthracis or B. cereus spores. The assay allowed an accurate identification of B. anthracis in environmental samples. This study provides a rapid and reliable method for improving rapid identification of B. anthracis in field operational conditions.

ENU-induced phenovariance in mice: inferences from 587 mutations

Background We present a compendium of N-ethyl-N-nitrosourea (ENU)-induced mouse mutations, identified in our laboratory over a period of 10 years either on the basis of phenotype or whole genome and/or whole exome sequencing, and archived in the Mutagenetix database. Our purpose is threefold: 1) to formally describe many point mutations, including those that were not previously disclosed in peer-reviewed publications; 2) to assess the characteristics of these mutations; and 3) to estimate the likelihood that a missense mutation induced by ENU will create a detectable phenotype. Findings In the context of an ENU mutagenesis program for C57BL/6J mice, a total of 185 phenotypes were tracked to mutations in 129 genes. In addition, 402 incidental mutations were identified and predicted to affect 390 genes. As previously reported, ENU shows strand asymmetry in its induction of mutations, particularly favoring T to A rather than A to T in the sense strand of coding regions and splice junctions. Some amino acid substitutions are far more likely to be damaging than others, and some are far more likely to be observed. Indeed, from among a total of 494 non-synonymous coding mutations, ENU was observed to create only 114 of the 182 possible amino acid substitutions that single base changes can achieve. Based on differences in overt null allele frequencies observed in phenotypic vs. non-phenotypic mutation sets, we infer that ENU-induced missense mutations create detectable phenotype only about 1 in 4.7 times. While the remaining mutations may not be functionally neutral, they are, on average, beneath the limits of detection of the phenotypic assays we applied. Conclusions Collectively, these mutations add to our understanding of the chemical specificity of ENU, the types of amino acid substitutions it creates, and its efficiency in causing phenovariance. Our data support the validity of computational algorithms for the prediction of damage caused by amino acid substitutions, and may lead to refined predictions as to whether specific amino acid changes are responsible for observed phenotypes. These data form the basis for closer in silico estimations of the number of genes mutated to a state of phenovariance by ENU within a population of G3 mice.

Identification of the bovine Arachnomelia mutation by massively parallel sequencing implicates sulfite oxidase (SUOX) in bone development

Arachnomelia is a monogenic recessive defect of skeletal development in cattle. The causative mutation was previously mapped to a approximately 7 Mb interval on chromosome 5. Here we show that array-based sequence capture and massively parallel sequencing technology, combined with the typical family structure in livestock populations, facilitates the identification of the causative mutation. We re-sequenced the entire critical interval in a healthy partially inbred cow carrying one copy of the critical chromosome segment in its ancestral state and one copy of the same segment with the arachnomelia mutation, and we detected a single heterozygous position. The genetic makeup of several partially inbred cattle provides extremely strong support for the causality of this mutation. The mutation represents a single base insertion leading to a premature stop codon in the coding sequence of the SUOX gene and is perfectly associated with the arachnomelia phenotype. Our findings suggest an important role for sulfite oxidase in bone development.

Effective population size of an indigenous Swiss cattle breed estimated from linkage disequilibrium

Effective population size is an important parameter for the assessment of genetic diversity within a livestock population and its development over time. If pedigree information is not available, linkage disequilibrium (LD) analysis might offer an alternative perspective for the estimation of effective population size. In this study, 128 individuals of the Swiss Eringer breed were genotyped using the Illumina BovineSNP50 beadchip. We set bin size at 50 kb for LD analysis, assuming that LD for proximal single nucleotide polymorphism (SNP)-pairs reflects distant breeding history while LD from distal SNP-pairs would reflect near history. Recombination rates varied among different regions of the genome. The use of physical distances as an approximation of genetic distances (e.g. setting 1 Mb = 0.01 Morgan) led to an upward bias in LD-based estimates of effective population size for generations beyond 50, while estimates for recent history were unaffected. Correction for restricted sample size did not substantially affect these results. LD-based actual effective population size was estimated in the range of 87-149, whereas pedigree-based effective population size resulted in 321 individuals. For conservation purposes, requiring knowledge of recent history (<50 generations), approximation assuming constant recombination rate seemed adequate.

Genetic variations in IL28B and allergic disease in children

Environmental changes affecting the relationship between the developing immune system and microbial exposure have been implicated in the epidemic rise of allergic disease in developed countries. While early developmental differences in T cell function are well-recognised, there is now emerging evidence that this is related to developmental differences in innate immune function. In this study we sought to examine if differences associated with innate immunity contribute to the altered immune programming recognised in allergic children. Here, we describe for the first time, the association of carriage of the T allele of the tagging single nucleotide polymorphism rs12979860 3 kb upstream of IL28B, encoding the potent innate immune modulator type III interferon lambda (IFN-λ3), and allergy in children (p = 0.004; OR 4.56). Strikingly, the association between rs12979860 genotype and allergic disease is enhanced in girls. Furthermore, carriage of the T allele at rs12979860 correlates with differences in the pro-inflammatory profile during the first five years of life suggesting this contributes to the key differences in subsequent innate immune development in children who develop allergic disease. In the context of rising rates of disease, these immunologic differences already present at birth imply very early interaction between genetic predisposition and prenatal environmental influences.

Genetic variations in bile acid homeostasis are not overrepresented in alcoholic cirrhosis compared to patients with heavy alcohol abuse and absent liver disease

Increased serum bile salt levels have been associated to a single-nucleotide polymorphism in the bile salt export pump (BSEP; ABCB11) in several acquired cholestatic liver diseases but there is little evidence in alcoholic liver disease (ALD). Furthermore, a crosstalk between vitamin D and bile acid synthesis has recently been discovered. Whether this crosstalk has an influence on the course of ALD is unclear to date. Our aim was to analyse the role of genetic polymorphisms in BSEP and the vitamin D receptor gene (NR1I1) on the emergence of cirrhosis in patients with ALD. Therefore, 511 alcoholic patients (131 with cirrhosis and 380 without cirrhosis) underwent ABCB11 genotyping (rs2287622). Of these, 321 (131 with cirrhosis and 190 without cirrhosis) were also tested for NR1I1 polymorphisms (bat-haplotype: BsmI rs1544410, ApaI rs7975232 and TaqI rs731236). Frequencies of ABCB11 and NR1I1 genotypes and haplotypes were compared between alcoholic patients with and without cirrhosis and correlated to serum bile salt, bilirubin and aspartate aminotransferase levels in those with cirrhosis. Frequencies of ABCB11 and NR1I1 genotypes and haplotypes did not differ between the two subgroups and no significant association between genotypes/haplotypes and liver function tests could be determined for neither polymorphism. We conclude that ABCB11 and NR1I1 polymorphisms are obviously not associated with development of cirrhosis in patients with ALD.

Age-modulated association between prefrontal NAA and the BDNF gene

Brain-derived neurotrophic factor (BDNF) has been implicated in the pathophysiology of psychiatric and neurological disorders and in the mechanisms of antidepressant pharmacotherapy. Psychiatric and neurological conditions have also been associated with reduced brain levels of N-acetyl-aspartate (NAA), which has been used as a putative marker of neural integrity. However, few studies have explored the relationship between BDNF polymorphisms and NAA levels directly. Here, we present data from a single-voxel proton magnetic resonance spectroscopy study of 64 individuals and explore the relationship between BDNF polymorphisms and prefrontal NAA level. Our results indicate an association between a single nucleotide polymorphism (SNP) within BDNF, known as rs1519480, and reduced NAA level (p = 0.023). NAA levels were further predicted by age and Asian ancestry. There was a significant rs1519480 × age interaction on NAA level (p = 0.031). Specifically, the effect of rs1519480 on NAA level became significant at age ⩾34.17 yr. NAA level decreased with advancing age for genotype TT (p = 0.001) but not for genotype CT (p = 0.82) or CC (p = 0.34). Additional in silico analysis of 142 post-mortem brain samples revealed an association between the same SNP and reduced BDNF mRNA expression in the prefrontal cortex. The rs1519480 SNP influences BDNF mRNA expression and has an impact on prefrontal NAA level over time. This genetic mechanism may contribute to inter-individual variation in cognitive performance seen during normal ageing, as well as contributing to the risk for developing psychiatric and neurological conditions.

Globally dispersed Y chromosomal haplotypes in wild and domestic sheep

To date, investigations of genetic diversity and the origins of domestication in sheep have utilised autosomal microsatellites and variation in the mitochondrial genome. We present the first analysis of both domestic and wild sheep using genetic markers residing on the ovine Y chromosome. Analysis of a single nucleotide polymorphism (oY1) in the SRY promoter region revealed that allele A-oY1 was present in all wild bighorn sheep (Ovis canadensis), two subspecies of thinhorn sheep (Ovis dalli), European Mouflon (Ovis musimon) and the Barbary (Ammontragis lervia). A-oY1 also had the highest frequency (71.4%) within 458 domestic sheep drawn from 65 breeds sampled from Africa, Asia, Australia, the Caribbean, Europe, the Middle East and Central Asia. Sequence analysis of a second locus, microsatellite SRYM18, revealed a compound repeat array displaying fixed differences, which identified bighorn and thinhorn sheep as distinct from the European Mouflon and domestic animals. Combined genotypic data identified 11 male-specific haplotypes that represented at least two separate lineages. Investigation of the geographical distribution of each haplotype revealed that one (H6) was both very common and widespread in the global sample of domestic breeds. The remaining haplotypes each displayed more restricted and informative distributions. For example, H5 was likely founded following the domestication of European breeds and was used to trace the recent transportation of animals to both the Caribbean and Australia. A high rate of Y chromosomal dispersal appears to have taken place during the development of domestic sheep as only 12.9% of the total observed variation was partitioned between major geographical regions.

A HIDDEN MARKOV MODEL FOR JOINT ESTIMATION OF GENOTYPE AND COPY NUMBER IN HIGH-THROUGHPUT SNP CHIPS

Amplifications and deletions of chromosomal DNA, as well as copy-neutral loss of heterozygosity have been associated with diseases processes. High-throughput single nucleotide polymorphism (SNP) arrays are useful for making genome-wide estimates of copy number and genotype calls. Because neighboring SNPs in high throughput SNP arrays are likely to have dependent copy number and genotype due to the underlying haplotype structure and linkage disequilibrium, hidden Markov models (HMM) may be useful for improving genotype calls and copy number estimates that do not incorporate information from nearby SNPs. We improve previous approaches that utilize a HMM framework for inference in high throughput SNP arrays by integrating copy number, genotype calls, and the corresponding confidence scores when available. Using simulated data, we demonstrate how confidence scores control smoothing in a probabilistic framework. Software for fitting HMMs to SNP array data is available in the R package ICE.

A noncoding melanophilin gene (MLPH) SNP at the splice donor of exon 1 represents a candidate causal mutation for coat color dilution in dogs

Coat color dilution in several breeds of dog is characterized by a specific pigmentation phenotype and sometimes accompanied by hair loss and recurrent skin inflammation, the so-called color dilution alopecia or black hair follicular dysplasia. Coat color dilution (d) is inherited as a Mendelian autosomal recessive trait. In a previous study, MLPH polymorphisms showed perfect cosegregation with the dilute phenotype within breeds. However, different dilute haplotypes were found in different breeds, and no single polymorphism was identified in the coding sequence that was likely to be causative for the dilute phenotype. We resequenced the 5'-region of the canine MLPH gene and identified a strong candidate single nucleotide polymorphism within the nontranslated exon 1, which showed perfect association to the dilute phenotype in 65 dilute dogs from 7 different breeds. The A/G polymorphism is located at the last nucleotide of exon 1 and the mutant A-allele is predicted to reduce splicing efficiency 8-fold. An MLPH mRNA expression study using quantitative reverse transcriptase-polymerase chain reaction confirmed that dd animals had only about approximately 25% of the MLPH transcript compared with DD animals. These results provide preliminary evidence that the reported regulatory MLPH mutation might represent a causal mutation for coat color dilution in dogs.

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.

Exon splice enhancer mutation (GH-E32A) causes autosomal dominant growth hormone deficiency

CONTEXT AND OBJECTIVE: Alteration of exon splice enhancers (ESE) may cause autosomal dominant GH deficiency (IGHD II). Disruption analysis of a (GAA) (n) ESE motif within exon 3 by introducing single-base mutations has shown that single nucleotide mutations within ESE1 affect pre-mRNA splicing. DESIGN, SETTING, AND PATIENTS: Confirming the laboratory-derived data, a heterozygous splice enhancer mutation in exon 3 (exon 3 + 2 A-->C) coding for GH-E32A mutation of the GH-1 gene was found in two independent pedigrees, causing familial IGHD II. Because different ESE mutations have a variable impact on splicing of exon 3 of GH and therefore on the expression of the 17.5-kDa GH mutant form, the GH-E32A was studied at the cellular level. INTERVENTIONS AND RESULTS: The splicing of GH-E32A, assessed at the protein level, produced significantly increased amounts of 17.5-kDa GH isoform (55% of total GH protein) when compared with the wt-GH. AtT-20 cells coexpressing both wt-GH and GH-E32A presented a significant reduction in cell proliferation as well as GH production after forskolin stimulation when compared with the cells expressing wt-GH. These results were complemented with confocal microscopy analysis, which revealed a significant reduction of the GH-E32A-derived isoform colocalized with secretory granules, compared with wt-GH. CONCLUSION: GH-E32A mutation found within ESE1 weakens recognition of exon 3 directly, and therefore, an increased production of the exon 3-skipped 17.5-kDa GH isoform in relation to the 22-kDa, wt-GH isoform was found. The GH-E32A mutant altered stimulated GH production as well as cell proliferation, causing IGHD II.

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.

Predictive value of known and novel alleles of CYP2B6 for efavirenz plasma concentrations in HIV-infected individuals

To assess the association of CYP2B6 allelic diversity with efavirenz (EFV) pharmacokinetics, we performed extensive genotyping of 15 relevant single nucleotide polymorphism in 169 study participants, and full resequencing of CYP2B6 in individuals with abnormal EFV plasma levels. Seventy-seven (45.5%) individuals carried a known (CYP2B6*6, *11, *15, or *18) or new loss/diminished-function alleles. Resequencing defined two new loss-of-function alleles: allele *27 (marked by 593T>C [M198T]), that results in 85% decrease in enzyme activity and allele *28 (marked by 1132C>T), that results in protein truncation at arginine 378. Median AUC levels were 188.5 microg h/ml for individuals homozygous for a loss/diminished-function allele, 58.6 microg h/ml for carriers, and 43.7 microg h/ml for noncarriers (P<0.0001). Individuals with a poor metabolizer genotype had a likelihood ratio of 35 (95% CI, 11-110) of presenting very high EFV plasma levels. CYP2B6 poor metabolizer genotypes explain to a large extent EFV pharmacokinetics and identify individuals at risk of extremely elevated EFV plasma levels.

Severe phenotype with cis-acting heterozygous PMP22 mutations

We report on a 20-year-old male with severe Charcot-Marie-Tooth (CMT) disease and a de novo deletion (c.281delG, p.G94AfsX17) on the paternal PMP22 allele harboring c.353C>T (p.T118M). RNA-based sequence analysis confirmed the absence of nonsense-mediated decay and the presence of the mutant transcripts in Epstein-Barr virus-transformed lymphoblastoid cells of our patient. His clinical findings included early onset of polyneuropathy, loss of muscle mass with distal pareses, hammer toes, and progressive scoliosis. There was no neuropsychological alteration. Our results suggest that the deletion c.281delG alone is responsible for the severe CMT phenotype. To the best of our knowledge, this is the second report on a proven paternal origin of a de novo single-base mutation in the PMP22 gene.