Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

Forensic identification of urine samples: a comparison between nuclear and mitochondrial DNA markers.

Urine samples from 20 male volunteers of European Caucasian origin were stored at 4 degrees C over a 4-month period in order to compare the identification potential of nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) markers. The amount of nDNA recovered from urines dramatically declined over time. Consequently, nDNA likelihood ratios (LRs) greater than 1,000 were obtained for 100, 70 and 55% of the urines analysed after 6, 60 and 120 days, respectively. For the mtDNA, HVI and HVII sequences were obtained for all samples tested, whatever the period considered. Nevertheless, the highest mtDNA LR of 435 was relatively low compared to its nDNA equivalent. Indeed, LRs obtained with only three nDNA loci could easily exceed this value and are quite easier to obtain. Overall, the joint use of nDNA and mtDNA markers enabled the 20 urine samples to be identified, even after the 4-month period.

Identifying protein-coding genes in genomic sequences.

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

Identification of active oxalotrophic bacteria by Bromodeoxyuridine DNA labeling in a microcosm soil experiments

The oxalate-carbonate pathway (OCP) leads to a potential carbon sink in terrestrial environments. This process is linked to the activity of oxalotrophic bacteria. Although isolation and molecular characterizations are used to study oxalotrophic bacteria, these approaches do not give information on the active oxalotrophs present in soil undergoing the OCP. The aim of this study was to assess the diversity of active oxalotrophic bacteria in soil microcosms using the Bromodeoxyuridine (BrdU) DNA labeling technique. Soil was collected near an oxalogenic tree (Milicia excelsa). Different concentrations of calcium oxalate (0.5%, 1%, and 4% w/w) were added to the soil microcosms and compared with an untreated control. After 12days of incubation, a maximal pH of 7.7 was measured for microcosms with oxalate (initial pH 6.4). At this time point, a DGGE profile of the frc gene was performed from BrdU-labeled soil DNA and unlabeled soil DNA. Actinobacteria (Streptomyces- and Kribbella-like sequences), Gammaproteobacteria and Betaproteobacteria were found as the main active oxalotrophic bacterial groups. This study highlights the relevance of Actinobacteria as members of the active bacterial community and the identification of novel uncultured oxalotrophic groups (i.e. Kribbella) active in soils.

Comparison of human chromosome 21 conserved nongenic sequences (CNGs) with the mouse and dog genomes shows that their selective constraint is independent of their genic environment.

The analysis of conservation between the human and mouse genomes resulted in the identification of a large number of conserved nongenic sequences (CNGs). The functional significance of this nongenic conservation remains unknown, however. The availability of the sequence of a third mammalian genome, the dog, allows for a large-scale analysis of evolutionary attributes of CNGs in mammals. We have aligned 1638 previously identified CNGs and 976 conserved exons (CODs) from human chromosome 21 (Hsa21) with their orthologous sequences in mouse and dog. Attributes of selective constraint, such as sequence conservation, clustering, and direction of substitutions were compared between CNGs and CODs, showing a clear distinction between the two classes. We subsequently performed a chromosome-wide analysis of CNGs by correlating selective constraint metrics with their position on the chromosome and relative to their distance from genes. We found that CNGs appear to be randomly arranged in intergenic regions, with no bias to be closer or farther from genes. Moreover, conservation and clustering of substitutions of CNGs appear to be completely independent of their distance from genes. These results suggest that the majority of CNGs are not typical of previously described regulatory elements in terms of their location. We propose models for a global role of CNGs in genome function and regulation, through long-distance cis or trans chromosomal interactions.

Phylogeography of Mabuya maculilabris (Reptilia) from São Tomé Island (Gulf of Guinea) inferred from mtDNA sequences

The pattern of genetic variation of the lizard Mabuya maculilabris from São Tomé Island (Gulf of Guinea) was investigated using a combination of three mitochondrial DNA gene fragments. Forty-eight haplotypes were recovered among 66 individuals covering the whole island. The genealogy inferred from the most parsimonious network of haplotypes allows us to detect two main and long branches departing from the putative group of oldest haplotypes. The tips of these branches exhibit star-like phylogenies, which may indicate of recently expanded populations, most probably from a small number of founders. A nested clade analysis suggests a complex pattern of past events that gave rise to the extant geographical pattern found in the haplotype distribution: past and allopatric fragmentation, range expansion, restricted gene Xow and long-distance dispersal. These results are consistent with the complex geological history of the island where important volcanic activity with extensive lava Xows has occurred during several periods. Mismatch- distribution analysis and AMOVA also support these conclusions. Substantial genetic structuring among these lizards was detected as well as high levels of diVerentiation between the southern edge populations (particularly those from the Rolas Islet) and the remaining ones. However, variation is low relative to the geological age of the island. Our results indicate that patterns of variation observed in reptiles in other oceanic islands are not indicative of those observed in the islands of the Gulf of Guinea.

The complete amino acid sequence for brain beta spectrin (beta fodrin): relationship to globin sequences.

The amino acid sequence of mouse brain beta spectrin (beta fodrin), deduced from the nucleotide sequence of complementary DNA clones, reveals that this non-erythroid beta spectrin comprises 2363 residues, with a molecular weight of 274,449 Da. Brain beta spectrin contains three structural domains and we suggest the position of several functional domains including f-actin, synapsin I, ankyrin and spectrin self association sites. Analysis of deduced amino acid sequences indicated striking homology and similar structural characteristics of brain beta spectrin repeats beta 11 and beta 12 to globins. In vitro analysis has demonstrated that heme is capable of specific attachment to brain spectrin, suggesting possible new functions in electron transfer, oxygen binding, nitric oxide binding or heme scavenging.

DNA binding properties of peroxisome proliferator-activated receptor subtypes on various natural peroxisome proliferator response elements. Importance of the 5'-flanking region.

The three subtypes of the peroxisome proliferator-activated receptors (PPARalpha, beta/delta, and gamma) form heterodimers with the 9-cis-retinoic acid receptor (RXR) and bind to a common consensus response element, which consists of a direct repeat of two hexanucleotides spaced by one nucleotide (DR1). As a first step toward understanding the molecular mechanisms determining PPAR subtype specificity, we evaluated by electrophoretic mobility shift assays the binding properties of the three PPAR subtypes, in association with either RXRalpha or RXRgamma, on 16 natural PPAR response elements (PPREs). The main results are as follows. (i) PPARgamma in combination with either RXRalpha or RXRgamma binds more strongly than PPARalpha or PPARbeta to all natural PPREs tested. (ii) The binding of PPAR to strong elements is reinforced if the heterodimerization partner is RXRgamma. In contrast, weak elements favor RXRalpha as heterodimerization partner. (iii) The ordering of the 16 natural PPREs from strong to weak elements does not depend on the core DR1 sequence, which has a relatively uniform degree of conservation, but correlates with the number of identities of the 5'-flanking nucleotides with respect to a consensus element. This 5'-flanking sequence is essential for PPARalpha binding and thus contributes to subtype specificity. As a demonstration of this, the PPARgamma-specific element ARE6 PPRE is able to bind PPARalpha only if its 5'-flanking region is exchanged with that of the more promiscuous HMG PPRE.

Fluorescence in situ hybridization of potato somatohaploids and their somatic hybrid donors using two Solanum brevidens specific sequences

Selostus: Perunan somaattisten hybridien ja niiden somatohaploidien fluoresenssi in situ -hybridisaatio Solanum brevidens -lajin spesifisten sekvenssien avulla

Helical repeat of DNA in the region of homologous pairing.

The process of DNA strand exchange during general genetic recombination is initiated within protein-stabilized synaptic filaments containing homologous regions of interacting DNA molecules. The RecA protein in bacteria and its analogs in eukaryotic organisms start this process by forming helical filamentous complexes on single-stranded or partially single-stranded DNA molecules. These complexes then progressively bind homologous double-stranded DNA molecules so that homologous regions of single- and double-stranded DNA molecules become aligned in register while presumably winding around common axis. The topological assay presented herein allows us to conclude that in synaptic complexes containing homologous single- and double-stranded DNA molecules, all three DNA strands have a helicity of approximately 19 nt per turn.

The Vi element. A transposon-like repeated DNA sequence interspersed in the vitellogenin locus of Xenopus laevis.

A repeated DNA element in Xenopus laevis is described that is present in about 7500 copies dispersed throughout the genome. It was first identified in the 5' flanking region of one vitellogenin gene and was therefore named the Vi element. Seven copies are present within the vitellogenin gene region, three of them within introns of the genes A1, A2 and B2, and the other four copies in the gene flanking regions. Four of these copies have been sequenced. The Vi element is bounded by a well-conserved 13 base-pair inverted repeat; in addition, it is flanked by a three base-pair direct repeat that appears to be site-specific. The length of these four copies varies from 112 to 469 base-pairs; however, sequence homology between the different copies is very high. Their structural characteristics suggest that length heterogeneity may have arisen by either unequal recombinations, deletions or tandem duplications. Altogether, the characteristics and properties of the Vi element indicate that it might represent a mobile genetic element. One of the four copies sequenced is inserted close (position -535) to the transcription initiation site of the vitellogenin gene B2 in a region otherwise showing considerable homology with the closely related gene B1. Nevertheless, the presence of the Vi element does not seem to influence significantly the estrogen-controlled expression of gene B2. In addition, three alleles of this gene created by length polymorphism in intron 3 and in the Vi element inserted near the transcription initiation site are described.

Density-based hierarchical clustering of pyro-sequences on a large scale--the case of fungal ITS1.

MOTIVATION: Analysis of millions of pyro-sequences is currently playing a crucial role in the advance of environmental microbiology. Taxonomy-independent, i.e. unsupervised, clustering of these sequences is essential for the definition of Operational Taxonomic Units. For this application, reproducibility and robustness should be the most sought after qualities, but have thus far largely been overlooked. RESULTS: More than 1 million hyper-variable internal transcribed spacer 1 (ITS1) sequences of fungal origin have been analyzed. The ITS1 sequences were first properly extracted from 454 reads using generalized profiles. Then, otupipe, cd-hit-454, ESPRIT-Tree and DBC454, a new algorithm presented here, were used to analyze the sequences. A numerical assay was developed to measure the reproducibility and robustness of these algorithms. DBC454 was the most robust, closely followed by ESPRIT-Tree. DBC454 features density-based hierarchical clustering, which complements the other methods by providing insights into the structure of the data. AVAILABILITY: An executable is freely available for non-commercial users at ftp://ftp.vital-it.ch/tools/dbc454. It is designed to run under MPI on a cluster of 64-bit Linux machines running Red Hat 4.x, or on a multi-core OSX system. CONTACT: dbc454@vital-it.ch or nicolas.guex@isb-sib.ch.

Molecular phylogenetics of soricid shrews (Mammalia) based on mitochondrial cytochrome b gene sequences: with special reference to the Soricinae

Molecular phylogeny of soricid shrews (Soricidae, Eulipotyphla, Mammalia) based on 1140 bp mitochondrial cytochrome b gene (cytb) sequences was inferred by the maximum likelihood (ML) method. All 13 genera of extant Soricinae and two genera of Crocidurinae were included in the analyses. Anourosorex was phylogenetically distant from the main groupings within Soricinae and Crocidurinae in the ML tree. Thus, it could not be determined to which subfamily Anourosorex should be assigned: Soricinae, Crocidurinae or a new subfamily. Soricinae (excluding Anourosorex) should be divided into four tribes: Neomyini, Notiosoricini, Soricini and Blarinini. However, monophyly of Blarinini was not robust in the present data set. Also, branching orders among tribes of Soricinae and those among genera of Neomyini could not be determined because of insufficient phylogenetic information of the cytb sequences. For water shrews of Neomyini (Chimarrogale, Nectogale and Neomys), monophyly of Neomys and the Chimarrogale-Nectogale group could not be verified, which implies the possibility of multiple origins for the semi-aquatic mode of living among taxa within Neomyini. Episoriculus may contain several separate genera. Blarinella was included in Blarinini not Soricini, based on the cytb sequences, but the confidence level was rather low; hence more phylogenetic information is needed to determine its phylogenetic position. Furthermore, some specific problems of taxonomy of soricid shrews were clarified, for example phylogeny of local populations of Notiosorex crawfordi, Chimarrogale himalayica and Crocidura attenuata.

Identification and potential origin of invasive clawed frogs Xenopus (Anura: Pipidae) in Sicily based on mitochondrial and nuclear DNA

African clawed frogs of the widespread polytypic species Xenopus laevis Daudin, 1802 (ranging large parts of sub-Saharan Africa) have been spreading since the 1940s, and have established reproductive populations in Europe, Asia and the Americas, where they can have negative impact as competitors of native amphibians and as disease vectors for chytridomycosis or ranaviruses. Here we use two mitochondrial (cytochrome b, 16S rDNA) and one nuclear (RAG 1: Recombination Associated Gene 1) DNA markers to infer the potential origin of invasive clawed frogs from Sicily that represent the largest invasive population in Europe. Identical mtDNA haplotypes match with those of Xenopus laevis, and Sicilian clawed frogs very probably belong to a lineage from the Cape Region of South Africa, most likely originating from a laboratory stock. Nuclear data support this conclusion. Identical mtDNA sequences (cyt b, 16S) of frogs sampled across their range in Sicily suggest the occurrence of a single source population and a potential bottleneck at their release, but faster evolving multilocus nuclear data (microsatellites, SNPs) on the population genetics would be important in the future to better support this hypothesis