Molecular dating of caprines using ancient DNA sequences of Myotragus balearicus, an extinct endemic Balearic mammal

Background: Myotragus balearicus was an endemic bovid from the Balearic Islands (Western Mediterranean) that became extinct around 6,000-4,000 years ago. The Myotragus evolutionary lineage became isolated in the islands most probably at the end of the Messinian crisis, when the desiccation of the Mediterranean ended, in a geological date established at 5.35 Mya. Thus, the sequences of Myotragus could be very valuable for calibrating the mammalian mitochondrial DNA clock and, in particular, the tree of the Caprinae subfamily, to which Myotragus belongs. Results: We have retrieved the complete mitochondrial cytochrome b gene (1,143 base pairs), plus fragments of the mitochondrial 12S gene and the nuclear 28S rDNA multi-copy gene from a well preserved Myotragus subfossil bone. The best resolved phylogenetic trees, obtained with the cytochrome b gene, placed Myotragus in a position basal to the Ovis group. Using the calibration provided by the isolation of Balearic Islands, we calculated that the initial radiation of caprines can be dated at 6.2 ± 0.4 Mya. In addition, alpine and southern chamois, considered until recently the same species, split around 1.6 ± 0.3 Mya, indicating that the two chamois species have been separated much longer than previously thought. Conclusion: Since there are almost no extant endemic mammals in Mediterranean islands, the sequence of the extinct Balearic endemic Myotragus has been crucial for allowing us to use the Messinian crisis calibration point for dating the caprines phylogenetic tree.

Genetic characterization of the ABO blood group in Neandertals

Background: The high polymorphism rate in the human ABO blood group gene seems to be related to susceptibility to different pathogens. It has been estimated that all genetic variation underlying the human ABO alleles appeared along the human lineage, after the divergence from the chimpanzee lineage. A paleogenetic analysis of the ABO blood group gene in Neandertals allows us to directly test for the presence of the ABO alleles in these extinct humans. Results: We have analysed two male Neandertals that were retrieved under controlled conditions at the El Sidron site in Asturias (Spain) and that appeared to be almost free of modern human DNA contamination. We find a human specific diagnostic deletion for blood group O (O01 haplotype) in both Neandertal individuals. Conclusion: These results suggest that the genetic change responsible for the O blood group in humans predates the human and Neandertal divergence. A potential selective event associated with the emergence of the O allele may have therefore occurred after humans separated from their common ancestor with chimpanzees and before the human-Neandertal population divergence.

Statistical properties of new neutrality tests against population growth.

A number of statistical tests for detecting population growth are described. We compared the statistical power of these tests with that of others available in the literature. The tests evaluated fall into three categories: those tests based on the distribution of the mutation frequencies, on the haplotype distribution, and on the mismatch distribution. We found that, for an extensive variety of cases, the most powerful tests for detecting population growth are Fu"s FS test and the newly developed R2 test. The behavior of the R2 test is superior for small sample sizes, whereas FS is better for large sample sizes. We also show that some popular statistics based on the mismatch distribution are very conservative. Key words: population growth, population expansion, coalescent simulations, neutrality tests

Análisis del ADN mitocondrial de dos muestras del yacimiento paleolítico de El Pirulejo.

En la presente publicación se presentan los resultados del análisis del ADN mitocondrial (ADNmt) de dos individuos del nivel magdaleniense de El Pirulejo (Córdoba) (Asquerino et al., 1991). Las secuencias obtenidas se encuadran dentro de la variabilidad genética del hombre moderno actual, sin presentar similitud alguna con las secuencias de ADN mitocondrial de neandertal publicadas hasta la fecha. La distribución actual de las variantes mitocondriales encontradas en las dos muestras de El Pirulejo es similar a la encontrada en otros dos individuos Paleolíticos (Caramelli et al., 2003), y coherente con el modelo de expansión del hombre anatómicamente moderno desde África.

Ancient DNA analysis of 8000 B.C. near eastern farmers supports an early neolithic pioneer maritime colonization of Mainland Europe through Cyprus and the Aegean Islands

The genetic impact associated to the Neolithic spread in Europe has been widely debated over the last 20 years. Within this context, ancient DNA studies have provided a more reliable picture by directly analyzing the protagonist populations at different regions in Europe. However, the lack of available data from the original Near Eastern farmers has limited the achieved conclusions, preventing the formulation of continental models of Neolithic expansion. Here we address this issue by presenting mitochondrial DNA data of the original Near-Eastern Neolithic communities with the aim of providing the adequate background for the interpretation of Neolithic genetic data from European samples. Sixty-three skeletons from the Pre Pottery Neolithic B (PPNB) sites of Tell Halula, Tell Ramad and Dja'de El Mughara dating between 8,700-6,600 cal. B.C. were analyzed, and 15 validated mitochondrial DNA profiles were recovered. In order to estimate the demographic contribution of the first farmers to both Central European and Western Mediterranean Neolithic cultures, haplotype and haplogroup diversities in the PPNB sample were compared using phylogeographic and population genetic analyses to available ancient DNA data from human remains belonging to the Linearbandkeramik-Alföldi Vonaldiszes Kerámia and Cardial/Epicardial cultures. We also searched for possible signatures of the original Neolithic expansion over the modern Near Eastern and South European genetic pools, and tried to infer possible routes of expansion by comparing the obtained results to a database of 60 modern populations from both regions. Comparisons performed among the 3 ancient datasets allowed us to identify K and N-derived mitochondrial DNA haplogroups as potential markers of the Neolithic expansion, whose genetic signature would have reached both the Iberian coasts and the Central European plain. Moreover, the observed genetic affinities between the PPNB samples and the modern populations of Cyprus and Crete seem to suggest that the Neolithic was first introduced into Europe through pioneer seafaring colonization.

Avaluació de l’estructuració genètica humana a la Península Ibèrica i a d’altres regions del sudoest europeu: implicacions poblacionals i epidemiològiques

Els avenços en tècniques de genotipat de polimorfismes genètics a gran escala estan liderant una revolució en el camp de l’epidemiologia genètica i la genètica de poblacions humanes. La informació aportada per aquestes tècniques ha evidenciat l’existència d’estructuracions poblacionals que poden augmentar l’error en els estudis d’associació a escala genòmica (GWAS, genome-wide association studies). Estudis recents han demostrat la presència d’aquestes estructuracions a nivell interregional i intrarregional a Europa. El present projecte ha avaluat el grau d’estructuració genètica en poblacions de la Península Ibèrica i altres regions del sudoest europeu (Itàlia i França) per quantificar l’impacte que aquesta potencial estructuració pot tenir en el disseny d’estudis d’associació GWAS i reconstruir la història demogràfica de les poblacions de la Mediterrània. Per aconseguir aquests objectius, s’han analitzat mostres de DNA de 770 individus de 26 poblacions de la Península Ibèrica, França, Itàlia i d’altres països de la Mediterrània. Aquestes mostres van ser genotipades per 240000 SNPs utilitzant l’array 250K StyI d’Affymetrix en el marc d’aquest projecte o mitjançant altres arrays d’Affymetrix en els projectes internacionals HapMap i POPRES. S’han realitzat anàlisis estadístiques incloent anàlisis de components principals, Fst, identitat per descendència, desequilibri de lligament, barreres genètiques, etc. Aquests resultats han permés construir un marc de referència de la variabilitat en aquesta regió, avaluar el seu impacte en estudis d’associació i proposar mesures per evitar l’increment de qualsevol tipus d’error (tipus I i II) en estudis nacionals i internacionals. A més, també han permés reconstruir la història de les poblacions humanes de la Mediterrània així com analitzar les seves relacions demogràfiques. Donada la duració limitada d’aquesta acció (24 mesos, d’octubre de 2010 a setembre de 2012), els resultats d’aquest projecte es troben actualment en fase de redacció i conduiran a diverses publicacions en revistes internacionals i a la preparació de comunicacions a congressos.

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Reassessing the role of mitochondrial DNA mutations in autism spectrum disorder

Background: There is increasing evidence that impairment of mitochondrial energy metabolism plays an important role in the pathophysiology of autism spectrum disorders (ASD; OMIM number: 209850). A significant proportion of ASD cases display biochemical alterations suggestive of mitochondrial dysfunction and several studies have reported that mutations in the mitochondrial DNA (mtDNA) molecule could be involved in the disease phenotype. Methods: We analysed a cohort of 148 patients with idiopathic ASD for a number of mutations proposed in the literature as pathogenic in ASD. We also carried out a case control association study for the most common European haplogroups (hgs) and their diagnostic single nucleotide polymorphisms (SNPs) by comparing cases with 753 healthy and ethnically matched controls.Results: We did not find statistical support for an association between mtDNA mutations or polymorphisms and ASD.Conclusions: Our results are compatible with the idea that mtDNA mutations are not a relevant cause of ASD and the frequent observation of concomitant mitochondrial dysfunction and ASD could be due to nuclear factors influencing mitochondrion functions or to a more complex interplay between the nucleus and the mitochondrion/mtDNA.

Assessment of microbial community structure changes by amplified ribosomal DNA restriction analysis (ARDRA)

Amplified ribosomal DNA restriction analysis (ARDRA) is a simple method based on restriction endonuclease digestion of the amplified bacterial 16S rDNA. In this study we have evaluated the suitability of this method to detect differences in activated sludge bacterial communities fed on domestic or industrial wastewater, and subject to different operational conditions. The ability of ARDRA to detect these differences has been tested in modified Ludzack-Ettinger (MLE) configurations. Samples from three activated sludge wastewater treatment plants (WWTPs) with the MLE configuration were collected for both oxic and anoxic reactors, and ARDRA patterns using double enzyme digestions AluI+MspI were obtained. A matrix of Dice similarity coefficients was calculated and used to compare these restriction patterns. Differences in the community structure due to influent characteristics and temperature could be observed, but not between the oxic and anoxic reactors of each of the three MLE configurations. Other possible applications of ARDRA for detecting and monitoring changes in activated sludge systems are also discussed

DnaSP, DNA polymorphism analyses by the coalescent and other methods

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 (~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

Antibiotic resistance genes in the bacteriophage DNA fraction of environmental samples

Antibiotic resistance is an increasing global problem resulting from the pressure of antibiotic usage, greater mobility of the population, and industrialization. Many antibiotic resistance genes are believed to have originated in microorganisms in the environment, and to have been transferred to other bacteria through mobile genetic elements. Among others, ß-lactam antibiotics show clinical efficacy and low toxicity, and they are thus widely used as antimicrobials. Resistance to ß-lactam antibiotics is conferred by ß-lactamase genes and penicillin-binding proteins, which are chromosomal- or plasmid-encoded, although there is little information available on the contribution of other mobile genetic elements, such as phages. This study is focused on three genes that confer resistance to ß-lactam antibiotics, namely two ß-lactamase genes (blaTEM and blaCTX-M9) and one encoding a penicillin-binding protein (mecA) in bacteriophage DNA isolated from environmental water samples. The three genes were quantified in the DNA isolated from bacteriophages collected from 30 urban sewage and river water samples, using quantitative PCR amplification. All three genes were detected in the DNA of phages from all the samples tested, in some cases reaching 104 gene copies (GC) of blaTEM or 102 GC of blaCTX-M and mecA. These values are consistent with the amount of fecal pollution in the sample, except for mecA, which showed a higher number of copies in river water samples than in urban sewage. The bla genes from phage DNA were transferred by electroporation to sensitive host bacteria, which became resistant to ampicillin. blaTEM and blaCTX were detected in the DNA of the resistant clones after transfection. This study indicates that phages are reservoirs of resistance genes in the environment.

Oligomerization and DNA binding of Ler, a master regulator of pathogenicity of enterohemorrhagic and enteropathogenic Escherichia coli

Ler is a DNA-binding, oligomerizable protein that regulates pathogenicity islands in enterohemorrhagic and enteropathogenic Escherichia coli strains. Ler counteracts the transcriptional silencing effect of H-NS, another oligomerizable nucleoid-associated protein. We studied the oligomerization of Ler in the absence and presence of DNA by atomic force microscopy. Ler forms compact particles with a multimodal size distribution corresponding to multiples of 35 units of Ler. DNA wraps around Ler particles that contain more than 1516 Ler monomers. The resulting shortening of the DNA contour length is in agreement with previous measurements of the length of DNA protected by Ler in footprinting assays. We propose that the repetition unit corresponds to the number of monomers per turn of a tight helical Ler oligomer. While the repressor (H-NS) and anti-repressor (Ler) have similar DNA-binding domains, their oligomerization domains are unrelated. We suggest that the different oligomerization behavior of the two proteins explains the opposite results of their interaction with the same or proximal regions of DNA.

A nuclear defect in the 4p16 region predisposes to multiple mitochondrial DNA deletions in families with Wolfram syndrome.

Wolfram syndrome is a progressive neurodegenerative disorder transmitted in an autosomal recessive mode. We report two Wolfram syndrome families harboring multiple deletions of mitochondrial DNA. The deletions reached percentages as high as 85-90% in affected tissues such as the central nervous system of one patient, while in other tissues from the same patient and from other members of the family, the percentages of deleted mitochondrial DNA genomes were only 1-10%. Recently, a Wolfram syndrome gene has been linked to markers on 4p16. In both families linkage between the disease locus and 4p16 markers gave a maximum multipoint lod score of 3.79 at theta = 0 (Pi<0.03) with respect to D4S431. In these families, the syndrome was caused by mutations in this nucleus-encoded gene which deleteriously interacts with the mitochondrial genome. This is the first evidence of the implication of both genomes in a recessive disease.

A nuclear defect in the 4p16 region predisposes to multiple mitochondrial DNA deletions in families with Wolfram syndrome.

Wolfram syndrome is a progressive neurodegenerative disorder transmitted in an autosomal recessive mode. We report two Wolfram syndrome families harboring multiple deletions of mitochondrial DNA. The deletions reached percentages as high as 85-90% in affected tissues such as the central nervous system of one patient, while in other tissues from the same patient and from other members of the family, the percentages of deleted mitochondrial DNA genomes were only 1-10%. Recently, a Wolfram syndrome gene has been linked to markers on 4p16. In both families linkage between the disease locus and 4p16 markers gave a maximum multipoint lod score of 3.79 at theta = 0 (Pi<0.03) with respect to D4S431. In these families, the syndrome was caused by mutations in this nucleus-encoded gene which deleteriously interacts with the mitochondrial genome. This is the first evidence of the implication of both genomes in a recessive disease.

DNA sequence of the gene encoding the Zc1 protein from Zea mays W64A

The Zein-2 component named Zc 1 corresponds to a storage protein of an apparent M.W. of 16 kDa present in maize endosperm.