DnaSP, DNA polymorphism analyses by the coalescent and other methods

SUMMARY: 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 (approximately 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     

DNAlive: a tool for the physical analysis of DNA at the genomic scale

SUMMARY: DNAlive is a tool for the analysis and graphical display of structural and physical characteristics of genomic DNA. The web server implements a wide repertoire of metrics to derive physical information from DNA sequences with a powerful interface to derive 3D information on large sequences of both naked and protein-bound DNAs. Furthermore, it implements a mesoscopic Metropolis code which allows the inexpensive study of the dynamic properties of chromatin fibers. In addition, our server also surveys other protein and genomic databases allowing the user to combine and explore the physical properties of selected DNA in the context of functional features annotated on those regions. AVAILABILITY: http://mmb.pcb.ub.es/DNAlive/ ; http://www.inab.org/     

生物学软件在线粒体DNA序列多态性分析中的应用

对线粒体DNA(mtDNA)序列进行多态性分析,将有利于从母系遗传的角度研究各人群的源流、迁移以及亲缘关系等情况。探讨利用Internet共享生物学软件对各种生物样本mtDNA高变区进行序列比对分析及多态性分析,结果显示DNAstar、ClustalX、DnaSP以及Mega等软件操作简便、易于判读、结果准确可靠,可广泛应用于人类学、系统发育学以及法医鉴定等领域的研究。     

Monitoring the long-term molecular epidemiology of the pneumococcus and detection of potential 'vaccine escape' strains

Background

While the pneumococcal protein conjugate vaccines reduce the incidence in invasive pneumococcal disease (IPD), serotype replacement remains a major concern. Thus, serotype-independent protection with vaccines targeting virulence genes, such as PspA, have been pursued. PspA is comprised of diverse clades that arose through recombination. Therefore, multi-locus sequence typing (MLST)-defined clones could conceivably include strains from multiple PspA clades. As a result, a method is needed which can both monitor the long-term epidemiology of the pneumococcus among a large number of isolates, and analyze vaccine-candidate genes, such as pspA, for mutations and recombination events that could result in ‘vaccine escape’ strains.

Methodology

We developed a resequencing array consisting of five conserved and six variable genes to characterize 72 pneumococcal strains. The phylogenetic analysis of the 11 concatenated genes was performed with the MrBayes program, the single nucleotide polymorphism (SNP) analysis with the DNA Sequence Polymorphism program (DnaSP), and the recombination event analysis with the recombination detection package (RDP).

Results

The phylogenetic analysis correlated with MLST, and identified clonal strains with unique PspA clades. The DnaSP analysis correlated with the serotype-specific diversity detected using MLST. Serotypes associated with more than one ST complex had a larger degree of sequence polymorphism than a serotype associated with one ST complex. The RDP analysis confirmed the high frequency of recombination events in the pspA gene.

Conclusions

The phylogenetic tree correlated with MLST, and detected multiple PspA clades among clonal strains. The genetic diversity of the strains and the frequency of recombination events in the mosaic gene, pspA were accurately assessed using the DnaSP and RDP programs, respectively. These data provide proof-of-concept that resequencing arrays could play an important role within research and clinical laboratories in both monitoring the molecular epidemiology of the pneumococcus and detecting ‘vaccine escape’ strains among vaccine-candidate genes.     

nucleR: a package for non-parametric nucleosome positioning

SUMMARY: nucleR is an R/Bioconductor package for a flexible and fast recognition of nucleosome positioning from next generation sequencing and tiling arrays experiments. The software is integrated with standard high-throughput genomics R packages and allows for in situ visualization as well as to export results to common genome browser formats. AVAILABILITY: Additional information and methodological details can be found at http://mmb.pcb.ub.es/nucleR     

The antitumor molecular mechanism of Alisma orientalis with c-myc DNA: multi-spectroscopic analysis and molecular simulation

Abstract

We prepared extracts of Alisma orientalis from Sichuan and Fujian Province, China. Based on the ratio of alisol B 23-acetate (23B) to alisol A 24-acetate (24A) in two Alisma orientalis extracts, we prepared two mixtures of 24A and 23B (24A:23B?=?1:3 or 1:10). The antitumor molecular mechanism of the monomers 24A and 23B, the two mixtures and the effective components of Alisma orientalis from different habitats were studied. The MTT assay suggested that the difference in the antitumor activity of Alisma orientalis from different habitats was correlated to the ratio of 24A to 23B. The multi-spectroscopic analysis suggested that the effective components, the monomers and mixtures interacted with c-myc DNA in a partial intercalation manner. The binding strength of the alisol acetates to c-myc DNA was consistent with the anticancer activity, indicating that c-myc DNA was the anticancer target. The molecular simulation indicated that the mixtures were all directly bound to different base pairs of c-myc DNA for a superimposed effect, which led to the binding strength of the mixtures to c-myc DNA was stronger than that of the monomers. The molecules in the 1:3 mixture were all bound to different base pairs of c-myc DNA. However, for the 1:10 mixture, seven molecules of 23B bound to the side chain of 24A, resulting in the mixture with a long chain structure which increased the steric hindrance of 24A. As a result, affinity between 24A and c-myc DNA in the 1:10 mixture was weaker than that in the 1:3 mixture.

The antitumor molecular mechanism of the alisol monomers 24A and 23B, the mixtures with different proportions and the effective components of Alisma orientalis from different habitats were studied. The order of the antitumor activity was as follows: Sichuan?>?Fujian, 24A-23B (1:3) > 24A-23B (1:10) > 23B?>?24A. The antitumor activity of Alisma orientalis from different habitats was consistent with the mixtures which were designed according to the contents of the active ingredients of the medicinal materials, indicating that the antitumor activity of Alisma orientalis from Sichuan is better than that from Fujian which is related to the contents of 24A and 23B and the proportion of 1:3 is better than 1:10. The binding strength of the mixtures to c-myc DNA was consistent with the anticancer activity. The mixtures were all directly bound to different base pairs of c-myc DNA for a superimposed effect, which led to the strength of the interaction of the mixtures to c-myc DNA was stronger than that of the monomers. For the 24A-23B (1:3) mixture, the four small molecules bound to c-myc DNA directly and interacted with different base pairs of c-myc DNA. While for the 24A-23B (1:10) mixture, 24A and three 23B molecules interacted with c-myc DNA, the remaining seven 23B molecules bound to the side chain of 24A, which increased the steric hindrance. The binding of the mixture to c-myc DNA was decreased.

Communicated by Ramaswamy H. Sarma     

Polymorphism of the angiotensin-converting enzyme gene in cardiovascular pathology

The aim of the study was to investigate influence of polymorphism of angiotensin-converting enzyme (ACE) gene on peculiarities of clinical process of such cardiovascular pathology as hypertrophic cardiomyopathy, coronary arterial disease and arterial hypertension. The polymorphism of ACE gene was studied in 98 patients: 38 with hypertrophic cardiomyopathy, 35 with coronary arterial disease and 25 with arterial hypertension. Nuclear DNA was extracted from blood leukocytes by phenol-chloroform method. Genotypes of ACE gene were determined by polymeraze chain reaction, followed with electrophoresis in agarose gel. It has been established, that I/D polymorphism of ACE gene has important modificative significance in clinical process at the mentioned diseases.     

DNA分子标记技术在濒危物种保护中的应用

近20年来,随着分子生物学技术的迅猛发展,涌现出一批高效、可靠的DNA分子标记技术.本文论述了限制性片段长度多态性、微卫星DNA、随机扩增多态性DNA、扩增片段长度多态性等DNA分子标记技术的基本原理及技术特点;同时,介绍了DNA分子标记在濒危物种种群遗传学研究、致危因素分析及保护策略的制定等保护生物学方面的应用.     

Genetic diversity and natural selection at the domain I of apical membrane antigen-1 (AMA-1) of Plasmodium falciparum in isolates from Iran

The apical membrane antigen-1 (AMA-1) of Plasmodium falciparum is a prime malaria asexual blood-stage vaccine candidate. Antigenic variation is one of the main obstacles in the development of a universal effective malaria vaccine. The extracellular region of P. falciparum AMA-1 (PfAMA-1) consists of three domains (I-III), of which the domain I is the most diverse region of this antigen. The objective of our study was to investigate and analyze the extent of genetic diversity and the effectiveness of natural selection at the AMA-1 domain I of P. falciparum in isolates from Iran. A fragment of ama-1 gene spanning domain I was amplified by nested PCR from 48 P. falciparum isolates collected from two major malaria endemic areas of Iran during 2009 to August 2010 and sequenced. Genetic polymorphism and statistical analyses were performed using DnaSP and MEGA software packages. Analysis of intrapopulation diversity revealed relatively high nucleotide and haplotype diversity at the PfAMA-1 domain I of Iranian isolates. Neutrality tests provided strong evidence of positive natural selection acting on the sequenced gene region. The findings also demonstrated that, in addition to natural selection, intragenic recombination may contribute to the diversity observed at the domain I. The results obtained will have significant implications in the design and the development of an AMA-1-based vaccine against falciparum malaria.     

Calign: aligning sequences with restricted affine gap penalties.

MOTIVATION: Given a genomic DNA sequence, it is still an open problem to determine its coding regions, i.e. the region consisting of exons and introns. The comparison of cDNA and genomic DNA helps the understanding of coding regions. For such an application, it might be adequate to use the restricted affine gap penalties which penalize long gaps with a constant penalty. RESULTS: Several techniques developed for solving the approximate string-matching problem are employed to yield efficient algorithms for computing the optimal alignment with restricted affine gap penalties. In particular, efficient algorithms can be derived based on the suffix automaton with failure transitions and on the diagonalwise monotonicity of the cost tables. We have implemented the above methods in C on Sun workstations running SunOS Unix. Preliminary experiments show that these approaches are very promising for aligning a cDNA sequence with a genomic DNA sequence. AVAILABILITY: Calign is available free of charge by anonymous ftp at: iubio.bio. indiana.edu, directory: molbio/align, files: calign.driver.c calign. c. Another URL reference for the files is http://iubio.bio.indiana.edu/soft/molbio/align/+ ++calign.c.     

Application of kDNA as a molecular marker to analyse Leishmania infantum diversity in Portugal

Around the Mediterranean basin Leishmania infantum is an important parasite causing canine leishmaniasis and visceral and cutaneous clinical forms in both immunocompetent and immunocompromised humans. Efficient monitoring and evaluation of epidemiology with discriminatory molecular markers are required. We investigated the genetic diversity of L. infantum in Portugal by polymerase chain amplification and restriction fragment length polymorphism analysis of kinetoplastid DNA, as molecular marker. We analysed 120 Portuguese isolates of L. infantum plus 16 other non-Portuguese isolates (as a reference group) from humans, dogs and sand flies. The Portuguese population showed a high degree of polymorphism with a total of 13 profiles identified. The predominant profile was A, which was only detected in the Portuguese samples.

The kinetoplastid DNA PCR-RFLP assay described here was suitable for use directly with biological samples and the profiles obtained were stable during long-term growth in vitro and in laboratory animals.     

A molecular basis for human hypersensitivity to aminoglycoside antibiotics.

We have investigated the distribution of mitochondrial DNA polymorphisms in a rare maternally transmitted genetic trait that causes hypersensitivity to aminoglycoside antibiotics, in the hope that a characterization of its molecular basis might provide a molecular and cellular understanding of aminoglycoside-induced deafness (AGD). Here we report that the frequency of a particular mitochondrial DNA polymorphism, 1555G, is associated nonrandomly with aminoglycoside-induced deafness in two Japanese pedigrees, bringing the frequency of this polymorphism to 5 occurrences in 5 pedigrees of AGD, and in 4 of 78 sporadic cases in which deafness was thought to be the result of aminoglycoside exposure; both frequencies are significantly different from the occurrence of this mutation in the hearing population, which was 0 in 414 individuals surveyed. The 1555G polymorphism occurred in none of 34 aminoglycoside-resistant individuals. We propose a specific molecular mechanism for aminoglycoside hypersensitivity in individuals carrying the 1555G polymorphism, based on the three-dimensional structure of the ribosome, in which the 1555G polymorphism favors aminoglycoside binding sterically, by increasing access to the the ribosome cleft.     

A new progressive-iterative algorithm for multiple structure alignment

MOTIVATION: Multiple structure alignments are becoming important tools in many aspects of structural bioinformatics. The current explosion in the number of available protein structures demands multiple structural alignment algorithms with an adequate balance of accuracy and speed, for large scale applications in structural genomics, protein structure prediction and protein classification. RESULTS: A new multiple structural alignment program, MAMMOTH-mult, is described. It is demonstrated that the alignments obtained with the new method are an improvement over previous manual or automatic alignments available in several widely used databases at all structural levels. Detailed analysis of the structural alignments for a few representative cases indicates that MAMMOTH-mult delivers biologically meaningful trees and conservation at the sequence and structural levels of functional motifs in the alignments. An important improvement over previous methods is the reduction in computational cost. Typical alignments take only a median time of 5 CPU seconds in a single R12000 processor. MAMMOTH-mult is particularly useful for large scale applications. AVAILABILITY: http://ub.cbm.uam.es/mammoth/mult.     

ITS/CHS1靶位序列分析结合SRAP分子标记技术鉴定Nannizzia incurvata临床分离株

以转录间隔区（internal transcribed spacer,ITS）和几丁质合成酶1（chitinase synthase 1,CHS1）为鉴定靶位,结合种系进化分析将2007年至2016年分离自头癣或面癣患儿的12株石膏样小孢子菌复合体成员（现被归入Nannizzia菌属）鉴定到种水平,其中5株为Nannizzia gypsea（Microsporum gypseum）,另7株为Nannizzia incurvata（Microsporum incurvatum）。在此基础上,通过序列相关扩增多态性（SRAP）分子标记进一步证实二者在DNA多位点存在扩增多态性。本研究证实在我国湖北地区存在Nannizzia incurvata,为我国开展Nannizzia菌属的分子流行病学研究提供了实验室依据和支持。     

An optimized algorithm for flux estimation from isotopomer distribution in glucose metabolites

MOTIVATION: Analysis of the conversion of (13)C glucose within the metabolic network allows the evaluation of the biochemical fluxes in interconnecting metabolic pathways. Such analyses require solving hundreds of equations with respect to individual isotopomer concentrations, and this assumes applying special software even for constructing the equations. The algorithm, proposed by others could be improved. METHOD: A C-code linked to the program written in Mathematica (Wolfram Research Inc.), constructs and solves differential equations for all isotopomer concentrations, using the general enzyme characteristics (K(m), equilibrium constant, etc.). This code uses innovative algorithm of determination for the isotopomers-products, thus essentially decreasing the computation time. Feasible metabolic fluxes are provided by the parameters of enzyme kinetics found from the data fitting. RESULTS: The software effectively evaluates metabolic fluxes based on the measured isotopomer distribution, as was illustrated by the analysis of glycolysis and pentose phosphate cycle. The mechanism of transketolase and transaldolase catalysis was shown to induce a specific kind of isotopomer re-distribution, which, despite the significance of its effect, usually is not taken into account. AVAILABILITY: The software could be freely downloaded from the site: http://bq.ub.es/bioqint/label_distribution/.     

Pushing Structural Information into the Yeast Interactome by High-Throughput Protein Docking Experiments

The last several years have seen the consolidation of high-throughput proteomics initiatives to identify and characterize protein interactions and macromolecular complexes in model organisms. In particular, more that 10,000 high-confidence protein-protein interactions have been described between the roughly 6,000 proteins encoded in the budding yeast genome (Saccharomyces cerevisiae). However, unfortunately, high-resolution three-dimensional structures are only available for less than one hundred of these interacting pairs. Here, we expand this structural information on yeast protein interactions by running the first-ever high-throughput docking experiment with some of the best state-of-the-art methodologies, according to our benchmarks. To increase the coverage of the interaction space, we also explore the possibility of using homology models of varying quality in the docking experiments, instead of experimental structures, and assess how it would affect the global performance of the methods. In total, we have applied the docking procedure to 217 experimental structures and 1,023 homology models, providing putative structural models for over 3,000 protein-protein interactions in the yeast interactome. Finally, we analyze in detail the structural models obtained for the interaction between SAM1-anthranilate synthase complex and the MET30-RNA polymerase III to illustrate how our predictions can be straightforwardly used by the scientific community. The results of our experiment will be integrated into the general 3D-Repertoire pipeline, a European initiative to solve the structures of as many as possible protein complexes in yeast at the best possible resolution. All docking results are available at http://gatealoy.pcb.ub.es/HT_docking/.     

Patterns of inheritance with RAPD molecular markers reveal novel types of polymorphism in the honey bee

Summary The polymerase chain reaction (PCR) was used to generate random amplified polymorphic DNA (RAPD) from honey bee DNA samples in order to follow the patterns of inheritance of RAPD markers in a haplodiploid insect. The genomic DNA samples from two parental bees, a haploid drone and a diploid queen, were screened for polymorphism with 68 different tennucleotide primers of random sequence. Parents were scored for the presence or absence of individual bands. An average of 6.3 bands and 1.3 polymorphisms for presence/absence were observed per primer between the parents. Thirteen of these primers were used to determine the inheritance of RAPD marker alleles in the resulting progeny and in haploid drones from a daughter queen. Four types of polymorphisms were observed. Polymorphisms for band presence/absence as well as for band brightness were inherited as dominant markers, meeting Mendelian expectations in haploid and diploid progeny. Polymorphisms for fragment-length were also observed. These segregated in a near 11 ratio in drone progeny. The last type of polymorphism was manifested as a diploid-specific band. Mixing of amplification products after PCR showed that the diploid-specific band was the result of heteroduplex formation from the DNA of alternate alleles in heterozygotes. In two of the four cases of heteroduplex formation, the alternative alleles were manifested as small fragment-length polymorphisms, resulting in co-dominant markers. This is the first demonstration that a proportion of RAPD markers are not inherited in a dominant fashion.     

Isolation, characterization, and molecular identification of bacteria from the red imported fire ant (Solenopsis invicta) midgut

Bacteria were isolated and cultured from the red imported fire ant (Solenopsis invicta) midgut. The small-subunit ribosomal RNA gene, (16s rRNA gene, approximately 1500 bp) was amplified from bacterial genomic DNA using the polymerase chain reaction and consensus sequence primers. Restriction fragment length polymorphism analysis revealed 10 unique profiles, indicating that at least 10 different bacteria are present in red imported fire ant midguts. The 16s rRNA gene sequence was determined for these isolates and queried against the NCBI genetic database. The results identified all isolates to at least the genus level. Antibiotic resistance profiles and biochemical activities were also determined for these species. This work provides the basis for a wider characterization of bacterial distributions in fire ant colonies and provides strains suitable for genetic manipulation to develop novel methods of fire ant control.     

PCR-based molecular markers for the fragrance gene in rice (Oryza sativa. L.)

The genomic DNA clone RG28, linked to the major fragrance gene of rice (fgr), was assessed for polymorphism in order to produce a PCR-based marker for fragrance. A small mono-nucleotide repeat, that was polymorphic between a pair of fragrant and non-fragrant cultivars, was identified and developed into a co-dominant PCR-based marker. The polymorphism-information-content determinations for three microsatellite markers, that have been genetically mapped near RG28, are also presented. These PCR-based markers will be highly useful in distinguishing fragrance-producing alleles from non-fragrance-producing alleles at the fgr locus. Received: 19 October 1999 / Accepted: 16 December 1999     

The molecular basis of high-altitude adaptation in deer mice

Elucidating genetic mechanisms of adaptation is a goal of central importance in evolutionary biology, yet few empirical studies have succeeded in documenting causal links between molecular variation and organismal fitness in natural populations. Here we report a population genetic analysis of a two-locus α-globin polymorphism that underlies physiological adaptation to high-altitude hypoxia in natural populations of deer mice, Peromyscus maniculatus. This system provides a rare opportunity to examine the molecular underpinnings of fitness-related variation in protein function that can be related to a well-defined selection pressure. We surveyed DNA sequence variation in the duplicated α-globin genes of P. maniculatus from high- and low-altitude localities (i) to identify the specific mutations that may be responsible for the divergent fine-tuning of hemoglobin function and (ii) to test whether the genes exhibit the expected signature of diversifying selection between populations that inhabit different elevational zones. Results demonstrate that functionally distinct protein alleles are maintained as a long-term balanced polymorphism and that adaptive modifications of hemoglobin function are produced by the independent or joint effects of five amino acid mutations that modulate oxygen-binding affinity.