Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA.

DNA that survives in museum specimens, bones and other tissues recovered by archaeologists is invariably fragmented and chemically modified. The extent to which such modifications accumulate over time is largely unknown but could potentially be used to differentiate between endogenous old DNA and present-day DNA contaminating specimens and experiments. Here we examine mitochondrial DNA sequences from tissue remains that vary in age between 18 and 60,000 years with respect to three molecular features: fragment length, base composition at strand breaks, and apparent C to T substitutions. We find that fragment length does not decrease consistently over time and that strand breaks occur preferentially before purine residues by what may be at least two different molecular mechanisms that are not yet understood. In contrast, the frequency of apparent C to T substitutions towards the 5'-ends of molecules tends to increase over time. These nucleotide misincorporations are thus a useful tool to distinguish recent from ancient DNA sources in specimens that have not been subjected to unusual or harsh treatments.

Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705.

BACKGROUND: Analysis of the first reported complete genome sequence of Bifidobacterium longum NCC2705, an actinobacterium colonizing the gastrointestinal tract, uncovered its proteomic relatedness to Streptomyces coelicolor and Mycobacterium tuberculosis. However, a rapid scrutiny by genometric methods revealed a genome organization totally different from all so far sequenced high-GC Gram-positive chromosomes. RESULTS: Generally, the cumulative GC- and ORF orientation skew curves of prokaryotic genomes consist of two linear segments of opposite slope: the minimum and the maximum of the curves correspond to the origin and the terminus of chromosome replication, respectively. However, analyses of the B. longum NCC2705 chromosome yielded six, instead of two, linear segments, while its dnaA locus, usually associated with the origin of replication, was not located at the minimum of the curves. Furthermore, the coorientation of gene transcription with replication was very low. Comparison with closely related actinobacteria strongly suggested that the chromosome of B. longum was misassembled, and the identification of two pairs of relatively long homologous DNA sequences offers the possibility for an alternative genome assembly proposed here below. By genometric criteria, this configuration displays all of the characters common to bacteria, in particular to related high-GC Gram-positives. In addition, it is compatible with the partially sequenced genome of DJO10A B. longum strain. Recently, a corrected sequence of B. longum NCC2705, with a configuration similar to the one proposed here below, has been deposited in GenBank, confirming our predictions. CONCLUSION: Genometric analyses, in conjunction with standard bioinformatic tools and knowledge of bacterial chromosome architecture, represent fast and straightforward methods for the evaluation of chromosome assembly.

Phylogenetic structures of the Holarctic Sorex araneus group and its relationships with S-samniticus, as inferred from mtDNA sequences

The shrews of the Sorex araneus group, characterized by the sexual chromosome complex XY1, Y2 have been intensively studied by morphological, karyotypical, and biochemical analyses. Nevertheless, the phylogenetic relationships among the species belonging to the araneus complex are still under debate, as different approaches gave often contradictory results. In this paper, partial nucleotide sequences of the mitochondrial DNA cytochrome b gene (1011 bp) were determined for 6 species of the araneus group from Eurasia and North America. We also included in the data set the sequences of Sorex samniticus, whose relationships with the araneus group remain controversial. Three other species representing two major karyological groups were also examined. Both parsimony and distance trees strongly support the monophyly of the araneus group. Sorex sumniticus is significantly more closely related to the araneus complex than to the other species included in the analysis. Based on the branching pattern within the araneus group, an attempt has been made to reconstruct the colonization history of the Holarctic region.

A single-base substitution within an intronic repetitive element causes dominant retinitis pigmentosa with reduced penetrance.

We report the study of a large American family displaying autosomal dominant retinitis pigmentosa with reduced penetrance, a form of hereditary retinal degeneration. Although the inheritance pattern and previous linkage mapping pointed to the involvement of the PRPF31 gene, extensive screening of all its exons and their boundaries failed in the past to reveal any mutation. In this work, we sequenced the entire PRPF31 genomic region by both the classical Sanger method and ultrahigh throughput (UHT) sequencing. Among the many variants identified, a single-base substitution (c.1374+654C>G) located deep within intron 13 and inside a repetitive DNA element was common to all patients and obligate asymptomatic carriers. This change created a new splice donor site leading to the synthesis of two mutant PRPF31 isoforms, degraded by nonsense-mediated mRNA decay. As a consequence, amounts of PRPF31 mRNA derived from the mutant allele were very reduced, with no evidence of mutant proteins being synthesized. Our results indicate that c.1374+654C>G causes retinitis pigmentosa via haploinsufficiency, similar to the vast majority of PRPF31 mutations described so far. We discuss the potential of UHT sequencing technologies in mutation screening and the continued identification of pathogenic splicing mutations buried deep within intronic regions.

Vitellogenesis and the vitellogenin gene family.

Vitellogenin is synthesized under estrogen control in the liver, extensively modified, transported to the ovary, and there processed to the yolk proteins lipovitellin and phosvitin. In the frog Xenopus laevis there are at least four distinct but related vitellogenin genes. The two genes A1 and A2 have a 95 percent sequence homology in their messenger RNA coding regions, and contain 33 introns that interrupt the coding region (exons) at homologous positions. Sequences and lengths of analogous introns differ, and many introns contain repetitive DNA elements. The introns in these two genes that have apparently arisen by duplication have diverged extensively by events that include deletions, insertions, and probably duplications. Rapid evolutionary change involving rearrangements and the presence of repeated DNA suggests that the bulk of the sequences within introns may not have any specific function.

Protein-protein interactions between adenovirus DNA polymerase and nuclear factor I mediate formation of the DNA replication preinitiation complex.

The in vitro adenovirus (Ad) DNA replication system provides an assay to study the interaction of viral and host replication proteins with the DNA template in the formation of the preinitiation complex. This initiation system requires in addition to the origin DNA sequences 1) Ad DNA polymerase (Pol), 2) Ad preterminal protein (pTP), the covalent acceptor for protein-primed DNA replication, and 3) nuclear factor I (NFI), a host cell protein identical to the CCAAT box-binding transcription factor. The interactions of these proteins were studied by coimmunoprecipitation and Ad origin DNA binding assays. The Ad Pol can bind to origin sequences only in the presence of another protein which can be either pTP or NFI. While NFI alone can bind to its origin recognition sequence, pTP does not specifically recognize DNA unless Ad Pol is present. Thus, protein-protein interactions are necessary for the targetting of either Ad Pol or pTP to the preinitiation complex. DNA footprinting demonstrated that the Ad DNA site recognized by the pTP.Pol complex was within the first 18 bases at the end of the template which constitutes the minimal origin of replication. Mutagenesis studies have defined the Ad Pol interaction site on NFI between amino acids 68-150, which overlaps the DNA binding and replication activation domain of this factor. A putative zinc finger on the Ad Pol has been mutated to a product that fails to bind the Ad origin sequences but still interacts with pTP. These results indicate that both protein-protein and protein-DNA interactions mediate specific recognition of the replication origin by Ad DNA polymerase.

Proof of concept for microarray-based detection of DNA-binding oncogenes in cell extracts.

The function of DNA-binding proteins is controlled not just by their abundance, but mainly at the level of their activity in terms of their interactions with DNA and protein targets. Moreover, the affinity of such transcription factors to their target sequences is often controlled by co-factors and/or modifications that are not easily assessed from biological samples. Here, we describe a scalable method for monitoring protein-DNA interactions on a microarray surface. This approach was designed to determine the DNA-binding activity of proteins in crude cell extracts, complementing conventional expression profiling arrays. Enzymatic labeling of DNA enables direct normalization of the protein binding to the microarray, allowing the estimation of relative binding affinities. Using DNA sequences covering a range of affinities, we show that the new microarray-based method yields binding strength estimates similar to low-throughput gel mobility-shift assays. The microarray is also of high sensitivity, as it allows the detection of a rare DNA-binding protein from breast cancer cells, the human tumor suppressor AP-2. This approach thus mediates precise and robust assessment of the activity of DNA-binding proteins and takes present DNA-binding assays to a high throughput level.

Experimental analysis and computer prediction of CTF/NFI transcription factor DNA binding sites.

Accurate prediction of transcription factor binding sites is needed to unravel the function and regulation of genes discovered in genome sequencing projects. To evaluate current computer prediction tools, we have begun a systematic study of the sequence-specific DNA-binding of a transcription factor belonging to the CTF/NFI family. Using a systematic collection of rationally designed oligonucleotides combined with an in vitro DNA binding assay, we found that the sequence specificity of this protein cannot be represented by a simple consensus sequence or weight matrix. For instance, CTF/NFI uses a flexible DNA binding mode that allows for variations of the binding site length. From the experimental data, we derived a novel prediction method using a generalised profile as a binding site predictor. Experimental evaluation of the generalised profile indicated that it accurately predicts the binding affinity of the transcription factor to natural or synthetic DNA sequences. Furthermore, the in vitro measured binding affinities of a subset of oligonucleotides were found to correlate with their transcriptional activities in transfected cells. The combined computational-experimental approach exemplified in this work thus resulted in an accurate prediction method for CTF/NFI binding sites potentially functioning as regulatory regions in vivo.

DNA methylation by CcrM activates the transcription of two genes required for the division of Caulobacter crescentus.

DNA methylation regulates many processes, including gene expression, by superimposing secondary information on DNA sequences. The conserved CcrM enzyme, which methylates adenines in GANTC sequences, is essential to the viability of several Alphaproteobacteria. In this study, we find that Caulobacter crescentus cells lacking the CcrM enzyme accumulate low levels of the two conserved FtsZ and MipZ proteins, leading to a severe defect in cell division. This defect can be compensated by the expression of the ftsZ gene from an inducible promoter or by spontaneous suppressor mutations that promote FtsZ accumulation. We show that CcrM promotes the transcription of the ftsZ and mipZ genes and that the ftsZ and mipZ promoter regions contain a conserved CGACTC motif that is critical to their activities and to their regulation by CcrM. In addition, our results suggest that the ftsZ promoter has the lowest activity when the CGACTC motif is non-methylated, an intermediate activity when it is hemi-methylated and the highest activity when it is fully methylated. The regulation of ftsZ expression by DNA methylation may explain why CcrM is essential in a subset of Alphaproteobacteria.

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.

Photocontrolled DNA Binding of a Receptor-Targeted Organometallic Ruthenium(II) Complex

A photoactivated ruthenium(II) arene complex has been conjugated to two receptor-binding peptides, a dicarba analogue of octreotide and the Arg-Gly-Asp (RGD) tripeptide. These peptides can act as"tumor-targeting devices" since their receptors are overexpressed on the membranes of tumor cells. Both ruthenium-peptide conjugates are stable in aqueous solution in the dark, but upon irradiation with visible light, the pyridyl-derivatized peptides were selectively photodissociated from the ruthenium complex, as inferred by UV-vis and NMR spectroscopy. Importantly, the reactive aqua species generated from the conjugates, [(η6-p-cym)Ru(bpm)(H2O)]2+, reacted with the model DNA nucleobase 9-ethylguanine as well as with guanines of two DNA sequences, 5′dCATGGCT and 5′dAGCCATG. Interestingly, when irradiation was performed in the presence of the oligonucleotides, a new ruthenium adduct involving both guanines was formed as a consequence of the photodriven loss of p-cymene from the two monofunctional adducts. The release of the arene ligand and the formation of a ruthenated product with a multidentate binding mode might have important implications for the biological activity of such photoactivated ruthenium(II) arene complexes. Finally, photoreactions with the peptide-oligonucleotide hybrid, Phac-His-Gly-Met-linker-p5′dCATGGCT, also led to arene release and to guanine adducts, including a GG chelate. The lack of interaction with the peptide fragment confirms the preference of such organometallic ruthenium(II) complexes for guanine over other potential biological ligands, such as histidine or methionine amino acids.

The role of MAR elements in DNA recombination

The untargeted integration of foreign DNA into the mammalian cell genome, extensively used in gene therapy and biotechnology, remains an incompletely understood process. It is believed to be based on cellular DNA double strand break (DSB) repair machinery and to involve two major steps: i) the formation of long gene arrays (concatemers), and ii) recombination of the resulting concatemer with the genome. The main DSB repair pathways in eukaryotes include non-homologous end-joining (NHEJ), homologous recombination (HR), and microhomology-mediated end-joining (MMEJ). However, it is still not clear, which of these pathways are responsible for transgene integration. Here, we show that NHEJ is not the primary pathway used by mammalian cells in the transgene integration process, while the components of the HR pathway seem to be important for genomic integration but not concatemerization. Instead, concatemer formation appears to be mediated by a subset of the MMEJ pathway, termed synthesis-dependent MMEJ (SD-MMEJ). This mechanism also seems to be preferentially used for plasmid integration into the genome, as confirmed by the analysis of plasmid-to-genome junction sequences, which were found to display an SD-MMEJ pattern. Therefore, we propose the existence of two distinct SD-MMEJ subpathways, relying on different subsets of enzymes. One of these mechanisms appears to be responsible for concatemerization, while the other mechanism, partially dependent in HR enzymes, seems to mediate recombination with the genome. Previous studies performed by our group suggested that matrix attachment regions (MARs), which are epigenetic regulatory DNA elements that participate in the formation of chromatin boundaries and augment transcription, may mediate increased plasmid integration into the genome of CHO cells by stimulating DNA recombination. In the present work, we demonstrate that MAR-mediated plasmid integration results from the enhanced SD-MMEJ pathway. Analysis of transgene integration loci and junction DNA sequences validated the prevalent use of this pathway by the MAR elements to target plasmid DNA into gene-rich areas of the CHO genome. We propose that this finding should in the future help to engineer cells for improved recombinant protein production. In addition to investigating the process of transgene integration, we designed recombination assays to better characterize the components of the MMEJ and SD-MMEJ pathways. We also used CHO cells expressing cycle-sensitive reporter genes to demonstrate a potential role of HR proteins in the cell cycle regulation.

Phylogenetic relationships of Apodemus Kaup, 1829 (Rodentia: Muridae) species in the Eastern Mediterranean inferred from mitochondrial DNA, with emphasis on iranian species

Wood mice of the genus Apodemus are widely distributed in Eurasia, with the Eastern Mediterranean being considered as a hotspot. Indeed, numerous species have been documented in Iran, including A. witherbyi, A. hyrcanicus, A. uralensis, A. avicennicus, A. hermonensis, and A. arianus. In this study, 129 specimens were collected from different Iranian localities and two specimens from Afghanistan. The animals were identified taxonomically and their phylogenetic relationships were investigated using cytochrome b mitochondrial DNA sequences. Five species of the genus Apodemus were identified in Iran, including A. hyrcanicus, A. witherbyi, A. cf. ponticus, A. uralensis, and A. mystacinus, beside, A. pallipes from Afghanistan. This study found no evidence of A. flavicollis or A. sylvaticus in Iran, despite their occurrence in Turkey, shedding doubt on the status of A. flavicollis in Iran, Asia Minor, and the Levant. Phylogenetic analyses imply that A. witherbyi has priority over A. avicennicus, A. hermonensis, and A. iconicus. Estimation of the divergence time for these taxa suggests a separation at around 7.2 Ma for the subgenera Karstomys (including A. mystacinus and A. epimelas) and Sylvaemus (including A. flavicollis, A. sylvaticus, A. uralensis, A. pallipes, A. hyrcanicus, A. witherbyi, and A. cf. ponticus). Within the subgenus Karstomys, the divergence times for A. mystacinus and A. epimelas were between 3.0 and 6.1 Ma, and divergence times within the subgenus Sylvaemus were between 5.2 and 6.9 Ma for A. witherbyi and other species in this subgenus. It is postulated that vicariance events including the uplifting of the Zagros Mountains and Anatolian Plateau in the middle Miocene and climate oscillations during the Messinian Salinity Crisis besides formation of the Hyrcanian tertiary forests during the Neogene probably played substantial roles in the radiation and distribution of the genus Apodemus in the Eastern Mediterranean.

Fragmentation of contaminant and endogenous DNA in ancient samples determined by shotgun sequencing; prospects for human palaeogenomics

Despite the successful retrieval of genomes from past remains, the prospects for human palaeogenomics remain unclear because of the difficulty of distinguishing contaminant from endogenous DNA sequences. Previous sequence data generated on high-throughput sequencing platforms indicate that fragmentation of ancient DNA sequences is a characteristic trait primarily arising due to depurination processes that create abasic sites leading to DNA breaks.

Occurrence of virulence-related sequences and phylogenetic analysis of commensal and pathogenic avian Escherichia coli strains (APEC)

The presence of iron uptake (irp-2, fyuA, sitA, fepC, iucA), adhesion (iha, lpfA O157/O141, lpfA O157/O154, efa, toxB) and invasion (inv, ial-related DNA sequences and assignment to the four main Escherichia coli phylogenetic groups (A, B1, B2 e D) were determined in 30 commensal E. coli strains isolated from healthy chickens and in 49 APEC strains isolated from chickens presenting clinical signs of septicemia (n=24) swollen head syndrome (n=14) and omphalitis (n=11) by PCR. None of the strains presented DNA sequences related to the inv, ial, efa, and toxB genes. DNA sequences related to lpfA O157/O154, iucA, fepC, and irp-2 genes were significantly found among pathogenic strains, where iucA gene was associated with septicemia and swollen head syndrome and fepC and irp-2 genes were associated with swollen head syndrome strains. Phylogenetic typing showed that commensal and omphalitis strains belonged mainly to phylogenetic Group A and swollen head syndrome to phylogenetic Group D. Septicemic strains were assigned in phylogenetic Groups A and D. These data could suggest that clonal lineage of septicemic APEC strains have a multiple ancestor origin; one from a pathogenic bacteria ancestor and other from a non-pathogenic ancestor that evolved by the acquisition of virulence related sequences through horizontal gene transfer. Swollen head syndrome may constitute a pathogenic clonal group. By the other side, omphalitis strains probably constitute a non-pathogenic clonal group, and could cause omphalitis as an opportunistic infection. The sharing of virulence related sequences by human pathogenic E. coli and APEC strains could indicate that APEC strains could be a source of virulence genes to human strains and could represent a zoonotic risk.