Single-nucleotide polymorphism in two representative multidrug-resistant Mycobacterium bovis isolates collected from patients in a Spanish hospital harboring a human infection outbreak

Mycobacterium bovis is the etiological agent of tuberculosis in domestic and wild animals. Its involvement as a human pathogen has been highlighted again with the recent descriptions of transmission through dairy products (18), reactivation or primary infection in human immunodeficiency virus-infected patients (5), and association with meat industry workers, animal keepers, or hunters (3). Strains resistant to antituberculous drugs (M. bovis is naturally resistant to pyrazinamide) pose an additional risk (2). Several studies have demonstrated that mutations in target genes are associated with resistance to antituberculous drugs (4, 7, 10, 11, 16). However, most of them have been developed in Mycobacterium tuberculosis strains and limited data are available regarding M. bovis isolates. The aim of this study was to characterize by sequencing the main genes involved in antibiotic resistance in two multidrug-resistant (MDR) M. bovis isolates in a human outbreak detected in a hospital in Madrid that subsequently spread to several countries (5, 6, 15). The isolates were resistant to 11 drugs, but only their rpoB and katG genes have been analyzed so far (1, 14). We studied the first (93/R1) and last (95/R4) M. bovis isolates of this nosocomial outbreak, characterized by spoligotyping as SB0426 (hexacode 63-5F-5E-7F-FF-60 in the database at www.mbovis.org) (1, 13). Several genes involved in resistance to isoniazid (katG, ahpC, inhA, and the oxyR-ahpC intergenic region), rifampin (rpoB), streptomycin (rrs, rpsL), ethambutol (embB), and quinolones (gyrA) were studied. These genes, or fragments of genes, were amplified and sequenced as previously described (12). The sequence analysis revealed polymorphisms in five (ahpC, rpoB, rpsL, embB, and gyrA) out of nine analyzed genes (Table 1). Nucleotide substitutions in four genes cause a change in the encoded amino acid. Two additional synonymous mutations in ahpC and rpsL differentiated the first and last isolates from the outbreak.

A single nucleotide polymorphism in NEUROD1 is associated with production traits in nelore beef cattle.

Feed efficiency and carcass characteristics are late-measured traits. The detection of molecular markers associated with them can help breeding programs to select animals early in life, and to predict breeding values with high accuracy. The objective of this study was to identify polymorphisms in the functional and positional candidate gene NEUROD1 (neurogenic differentiation 1), and investigate their associations with production traits in reference families of Nelore cattle. A total of 585 steers were used, from 34 sires chosen to represent the variability of this breed. By sequencing 14 animals with extreme residual feed intake (RFI) values, seven single nucleotide polymorphisms (SNPs) in NEUROD1 were identified. The investigation of marker effects on the target traits RFI, backfat thickness (BFT), ribeye area (REA), average body weight (ABW), and metabolic body weight (MBW) was performed with a mixed model using the restricted maximum likelihood method. SNP1062, which changes cytosine for guanine, had no significant association with RFI or REA. However, we found an additive effect on ABW (P ≤ 0.05) and MBW (P ≤ 0.05), with an estimated allele substitution effect of -1.59 and -0.93 kg0.75, respectively. A dominant effect of this SNP for BFT was also found (P ≤ 0.010). Our results are the first that identify NEUROD1 as a candidate that affects BFT, ABW, and MBW. Once confirmed, the inclusion of this SNP in dense panels may improve the accuracy of genomic selection for these traits in Nelore beef cattle as this SNP is not currently represented on SNP chips.

A three-nucleotide mutation altering the Maize streak virus Rep pRBR-interaction motif reduces symptom severity in maize and partially reverts at high frequency without restoring pRBR-Rep binding

Geminivirus infectivity is thought to depend on interactions between the virus replication-associated proteins Rep or RepA and host retinoblastoma-related proteins (pRBR), which control cell-cycle progression. It was determined that the substitution of two amino acids in the Maize streak virus (MSV) RepA pRBR-interaction motif (LLCNE to LLCLK) abolished detectable RepA-pRBR interaction in yeast without abolishing infectivity in maize. Although the mutant virus was infectious in maize, it induced less severe symptoms than the wild-type virus. Sequence analysis of progeny viral DNA isolated from infected maize enabled detection of a high-frequency single-nucleotide reversion of C(601)A in the 3 nt mutated sequence of the Rep gene. Although it did not restore RepA-pRBR interaction in yeast, sequence-specific PCR showed that, in five out of eight plants, the C(601)A reversion appeared by day 10 post-inoculation. In all plants, the C(601)A revertant eventually completely replaced the original mutant population, indicating a high selection pressure for the single-nucleotide reversion. Apart from potentially revealing an alternative or possibly additional function for the stretch of DNA that encodes the apparently non-essential pRBR-interaction motif of MSV Rep, the consistent emergence and eventual dominance of the C(601)A revertant population might provide a useful tool for investigating aspects of MSV biology, such as replication, mutation and evolution rates, and complex population phenomena, such as competition between quasispecies and population turnover. © 2005 SGM.

Restoration of native folding of single-stranded DNA sequences through reverse mutations: an indication of a new epigenetic mechanism

We used in vivo (biological), in silico (computational structure prediction), and in vitro (model sequence folding) analyses of single-stranded DNA sequences to show that nucleic acid folding conservation is the selective principle behind a high-frequency single-nucleotide reversion observed in a three-nucleotide mutated motif of the Maize streak virus replication associated protein (Rep) gene. In silico and in vitro studies showed that the three-nucleotide mutation adversely affected Rep nucleic acid folding, and that the single-nucleotide reversion [C(601)A] restored wild-type-like folding. In vivo support came from infecting maize with mutant viruses: those with Rep genes containing nucleotide changes predicted to restore a wild-type-like fold [A(601)/G(601)] preferentially accumulated over those predicted to fold differently [C(601)/T(601)], which frequently reverted to A(601) and displaced the original population. We propose that the selection of native nucleic acid folding is an epigenetic effect, which might have broad implications in the evolution of plants and their viruses.

Facile mutant identification via a single parental backcross method and application of whole genome sequencing based mapping pipelines

Forward genetic screens have identified numerous genes involved in development and metabolism, and remain a cornerstone of biological research. However, to locate a causal mutation, the practice of crossing to a polymorphic background to generate a mapping population can be problematic if the mutant phenotype is difficult to recognize in the hybrid F2 progeny, or dependent on parental specific traits. Here in a screen for leaf hyponasty mutants, we have performed a single backcross of an Ethane Methyl Sulphonate (EMS) generated hyponastic mutant to its parent. Whole genome deep sequencing of a bulked homozygous F2 population and analysis via the Next Generation EMS mutation mapping pipeline (NGM) unambiguously determined the causal mutation to be a single nucleotide polymorphisim (SNP) residing in HASTY, a previously characterized gene involved in microRNA biogenesis. We have evaluated the feasibility of this backcross approach using three additional SNP mapping pipelines; SHOREmap, the GATK pipeline, and the samtools pipeline. Although there was variance in the identification of EMS SNPs, all returned the same outcome in clearly identifying the causal mutation in HASTY. The simplicity of performing a single parental backcross and genome sequencing a small pool of segregating mutants has great promise for identifying mutations that may be difficult to map using conventional approaches.

Extending genetic linkage analysis to diffusion tensor images to map single gene effects on brain fiber architecture

We extended genetic linkage analysis - an analysis widely used in quantitative genetics - to 3D images to analyze single gene effects on brain fiber architecture. We collected 4 Tesla diffusion tensor images (DTI) and genotype data from 258 healthy adult twins and their non-twin siblings. After high-dimensional fluid registration, at each voxel we estimated the genetic linkage between the single nucleotide polymorphism (SNP), Val66Met (dbSNP number rs6265), of the BDNF gene (brain-derived neurotrophic factor) with fractional anisotropy (FA) derived from each subject's DTI scan, by fitting structural equation models (SEM) from quantitative genetics. We also examined how image filtering affects the effect sizes for genetic linkage by examining how the overall significance of voxelwise effects varied with respect to full width at half maximum (FWHM) of the Gaussian smoothing applied to the FA images. Raw FA maps with no smoothing yielded the greatest sensitivity to detect gene effects, when corrected for multiple comparisons using the false discovery rate (FDR) procedure. The BDNF polymorphism significantly contributed to the variation in FA in the posterior cingulate gyrus, where it accounted for around 90-95% of the total variance in FA. Our study generated the first maps to visualize the effect of the BDNF gene on brain fiber integrity, suggesting that common genetic variants may strongly determine white matter integrity.

Further characterization and determination of the single amino acid change in the tsI138 reovirus thermosensitive mutant

Many temperature-sensitive mutants have been isolated in early studies of mammalian reovirus. However, the bio- logical properties and nature of the genetic alterations remain incompletely explored for most of these mutants. The mutation harbored by the tsI138 mutant was already assigned to the L3 gene encoding the l1 protein. In the present study, this mu- tant was further studied as a possible tool to establish the role of the putative l1 enzymatic activities in viral multiplication. It was observed that synthesis of viral proteins is only marginally reduced, while it was difficult to recover viral particles at the nonpermissive temperature. A single nucleotide substitution resulting in an amino acid change was found; the position of this amino acid is consistent with a probable defect in assembly of the inner capsid at the nonpermissive temperature.

Heterozygosity and diversity analysis using mapped single nucelotide polymorphisms in a faba bean inbreeding programme

The aim of this study was to convert existing faba bean (Vicia faba L.) single nucleotide polymorphism (SNP) markers from cleaved amplification polymorphic sequences and SNaPshot® formats, which are expensive and time-consuming, to the more convenient KBiosciences competitive allele‐specific PCR (KASP) assay format. Out of 80 assays designed, 75 were validated, though a core set of 67 of the most robust markers is recommended for further use. The 67 best KASP SNP assays were used across two generations of single seed descent to detect unintended outcrossing and to track and quantify loss of heterozygosity, a capability that will significantly increase the efficiency and performance of pure line production and maintenance. This same set of assays was also used to examine genetic relationships between the 67 members of the partly inbred panel, and should prove useful for line identification and diversity studies in the future.

A single multiplex PCR and SNaPshot minisequencing reaction of 42 SNPs to classify admixture populations into mitochondrial DNA haplogroups

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A single multiplex SNaPshot reaction of 42 SNPs to classify admixture populations into mitochondrial DNA haplogroups

SNaPshot minisequencing reaction is in increasing use because of its fast detection of many polymorphisms in a single assay. In this work we described a highly sensitive single nucleotide polymorphisms (SNPs) typing method with detection of 42 mitochondrial DNA (mtDNA) SNPs in a single PCR and SNaPshot multiplex reaction in order to allow haplogroup classification in Latin American admixture population. We validated the panel typing 160 Brazilian individuals. DNA was extracted from blood spotted on filter paper using Chelex protocol. Forty SNPs were selected targeting haplogroup-specific mutations in Europeans, Africans and Asians (only precursors of Native Americans haplogroups A2, B2, C1, and D1) and two non-coding SNPs were chosen to increase the power of discrimination between individuals (SNPs positions 16,519 and 16,362). It was done using a modified version of a previously published multiplex SNaPshot minisequencing reaction established to resolve European haplogroups, adding SNPs targeting Africans (L0, L1, L2, L3, and L*) and Asians (A, B, C, and D) haplogroups based on SNPs described at PhyloTree.org build 2. PCR primers were designed using PerlPrimer software and checked with the Autodimer program. Thirty-three primer-pairs were used to amplify 42 SNPs. Using this panel, we were able to successfully classify 160 individuals into their correct haplogroups. Complete SNP profiles were obtained from 10. pg of total DNA. We conclude that it is possible to build and genotype more than 40 mtDNA SNPs in a single multiplex PCR and SNaPshot reaction, with sensitivity and reliability, resolving haplogroup classification in admixture populations. © 2011.

Identificação de SNPs (Single Nucleiotide Polymorphisms) no gene colina monooxigenase relacionado ao metabolismo da glicina betaína em Eucalyptus

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Identificação de Single Nucleotilde Polymorphisms (SMPs) no gene Nove-cis-epoxicaroteníde dioxigenase (NCED) em Eucalyptus

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A macroscopic kinetic model for DNA polymerase elongation and high-fidelity nucleotide election

The enzymatically catalyzed template-directed extension of ssDNA/primer complex is an impor-tant reaction of extraordinary complexity. The DNA polymerase does not merely facilitate the insertion of dNMP, but it also performs rapid screening of substrates to ensure a high degree of fidelity. Several kinetic studies have determined rate constants and equilibrium constants for the elementary steps that make up the overall pathway. The information is used to develop a macro-scopic kinetic model, using an approach described by Ninio [Ninio J., 1987. Alternative to the steady-state method: derivation of reaction rates from first-passage times and pathway probabili-ties. Proc. Natl. Acad. Sci. U.S.A. 84, 663–667]. The principle idea of the Ninio approach is to track a single template/primer complex over time and to identify the expected behavior. The average time to insert a single nucleotide is a weighted sum of several terms, in-cluding the actual time to insert a nucleotide plus delays due to polymerase detachment from ei-ther the ternary (template-primer-polymerase) or quaternary (+nucleotide) complexes and time delays associated with the identification and ultimate rejection of an incorrect nucleotide from the binding site. The passage times of all events and their probability of occurrence are ex-pressed in terms of the rate constants of the elementary steps of the reaction pathway. The model accounts for variations in the average insertion time with different nucleotides as well as the in-fluence of G+C content of the sequence in the vicinity of the insertion site. Furthermore the model provides estimates of error frequencies. If nucleotide extension is recognized as a compe-tition between successful insertions and time delaying events, it can be described as a binomial process with a probability distribution. The distribution gives the probability to extend a primer/template complex with a certain number of base pairs and in general it maps annealed complexes into extension products.

Investigating the role of single point mutations in the human proteome: a computational study

In the post genomic era with the massive production of biological data the understanding of factors affecting protein stability is one of the most important and challenging tasks for highlighting the role of mutations in relation to human maladies. The problem is at the basis of what is referred to as molecular medicine with the underlying idea that pathologies can be detailed at a molecular level. To this purpose scientific efforts focus on characterising mutations that hamper protein functions and by these affect biological processes at the basis of cell physiology. New techniques have been developed with the aim of detailing single nucleotide polymorphisms (SNPs) at large in all the human chromosomes and by this information in specific databases are exponentially increasing. Eventually mutations that can be found at the DNA level, when occurring in transcribed regions may then lead to mutated proteins and this can be a serious medical problem, largely affecting the phenotype. Bioinformatics tools are urgently needed to cope with the flood of genomic data stored in database and in order to analyse the role of SNPs at the protein level. In principle several experimental and theoretical observations are suggesting that protein stability in the solvent-protein space is responsible of the correct protein functioning. Then mutations that are found disease related during DNA analysis are often assumed to perturb protein stability as well. However so far no extensive analysis at the proteome level has investigated whether this is the case. Also computationally methods have been developed to infer whether a mutation is disease related and independently whether it affects protein stability. Therefore whether the perturbation of protein stability is related to what it is routinely referred to as a disease is still a big question mark. In this work we have tried for the first time to explore the relation among mutations at the protein level and their relevance to diseases with a large-scale computational study of the data from different databases. To this aim in the first part of the thesis for each mutation type we have derived two probabilistic indices (for 141 out of 150 possible SNPs): the perturbing index (Pp), which indicates the probability that a given mutation effects protein stability considering all the “in vitro” thermodynamic data available and the disease index (Pd), which indicates the probability of a mutation to be disease related, given all the mutations that have been clinically associated so far. We find with a robust statistics that the two indexes correlate with the exception of all the mutations that are somatic cancer related. By this each mutation of the 150 can be coded by two values that allow a direct comparison with data base information. Furthermore we also implement computational methods that starting from the protein structure is suited to predict the effect of a mutation on protein stability and find that overpasses a set of other predictors performing the same task. The predictor is based on support vector machines and takes as input protein tertiary structures. We show that the predicted data well correlate with the data from the databases. All our efforts therefore add to the SNP annotation process and more importantly found the relationship among protein stability perturbation and the human variome leading to the diseasome.

Mapping nucleosome position at single base-pair resolution by using site-directed hydroxyl radicals.

A base-pair resolution method for determining nucleosome position in vitro has been developed to com- plement existing, less accurate methods. Cysteaminyl EDTA was tethered to a recombinant histone octamer via a mutant histone H4 with serine 47 replaced by cysteine. When assembled into nucleosome core particles, the DNA could be cut site specifically by hydroxyl radical-catalyzed chain scission by using the Fenton reaction. Strand cleavage occurs mainly at a single nucleotide close to the dyad axis of the core particle, and assignment of this location via the symmetry of the nucleosome allows base-pair resolution mapping of the histone octamer position on the DNA. The positions of the histone octamer and H3H4 tetramer were mapped on a 146-bp Lytechinus variegatus 5S rRNA sequence and a twofold-symmetric derivative. The weakness of translational determinants of nucleosome positioning relative to the overall affinity of the histone proteins for this DNA is clearly demonstrated. The predominant location of both histone octamer and H3H4 tetramer assembled on the 5S rDNA is off center. Shifting the nucleosome core particle position along DNA within a conserved rotational phase could be induced under physiologically relevant conditions. Since nucleosome shifting has important consequences for chromatin structure and gene regulation, an approach to the thermodynamic characterization of this movement is proposed. This mapping method is potentially adaptable for determining nucleosome position in chromatin in vivo.