Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.

Availability of sequences of aquatic and terrestrial taxa in genetic databases (GenBank and BOLD - the Barcode of Life Database) in 2010, 2012 and 2016

Correct species identifications are of tremendous importance for invasion ecology, as mistakes could lead to misdirecting limited resources against harmless species or inaction against problematic ones. DNA barcoding is becoming a promising and reliable tool for species identifications, however the efficacy of such molecular taxonomy depends on gene region(s) that provide a unique sequence to differentiate among species and on availability of reference sequences in existing genetic databases. Here, we assembled a list of aquatic and terrestrial non-indigenous species (NIS) and checked two leading genetic databases for corresponding sequences of six genome regions used for DNA barcoding. The genetic databases were checked in 2010, 2012, and 2016. All four aquatic kingdoms (Animalia, Chromista, Plantae and Protozoa) were initially equally represented in the genetic databases, with 64, 65, 69, and 61% of NIS included, respectively. Sequences for terrestrial NIS were present at rates of 58 and 78% for Animalia and Plantae, respectively. Six years later, the number of sequences for aquatic NIS increased to 75, 75, 74, and 63% respectively, while those for terrestrial NIS increased to 74 and 88% respectively. Genetic databases are marginally better populated with sequences of terrestrial NIS of plants compared to aquatic NIS and terrestrial NIS of animals. The rate at which sequences are added to databases is not equal among taxa. Though some groups of NIS are not detectable at all based on available data - mostly aquatic ones - encouragingly, current availability of sequences of taxa with environmental and/or economic impact is relatively good and continues to increase with time.

PubDNA Finder in a Nutshell. Searching the Life Sciences Literature with Sequences of Nucleic Acids

Biomedical researchers and clinicians working with molecular technologies in routine clinical practice often need to review the available literature to gather information regarding specific sequences of nucleic acids. This includes, for instance, finding articles related to a concrete DNA sequence, or identifying empirically-validated primer/probe sequences to evaluate the presence of different micro-organisms. Unfortunately, these hard and time-consuming tasks often need to be manually performed by researchers themselves since no publicly available biomedical literature search engine, e.g. PubMed, PubMed Central (PMC), etc., provides the required search functionalities. In this article, we describe PubDNA Finder, a web service that enables users to perform advanced searches on PubMed Central-indexed full text articles with sequences of nucleic acids

Mitochondrial DNA reveals hidden diversity for tardigrades from across the Antarctic realm

Antarctica contains some of the most challenging environmental conditions on the planet due to freezing temperatures, prolonged winters and lack of liquid water. Whereas 99.7% of Antarctica is permanently covered by ice and snow, some coastal areas and mountain ridges have remained ice-free and are able to sustain populations of microinvertebrates. Tardigrades are one of the more dominant groups of microfauna in soil and limno-terrestrial habitats, but little is known of their diversity and distribution across Antarctica. Here, we examine tardigrades sampled from across an extensive region of continental Antarctica, and analyse and compare their partial mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences with those from the Antarctic Peninsula, maritime and sub-Antarctica, Tierra del Fuego and other worldwide locations in order to recognise operational taxonomic units (OTUs). From 439 new tardigrade COI sequences, we identified 98 unique haplotypes (85 from Antarctica) belonging to Acutuncus, Diphascon, Echiniscus, Macrobiotus, Milnesium and unidentified Parachela. Operational taxonomic units were delimited by Poisson tree processes and general mixed Yule coalescent methods, resulting in 58 and 55 putative species, respectively. Most tardigrades appear to be locally endemic (i.e. restricted to a single geographic region), but some (e.g. Acutuncus antarcticus (Richters, 1904)) are widespread across continental Antarctica. Our molecular results reveal: (i) greater diversity than has previously been appreciated with distinct OTUs that potentially represent undescribed species, and (ii) a lack of connectivity between most OTUs from continental Antarctica and those from other Antarctic geographical zones.

The theoretical limits of DNA sequence discrimination by linked polyamides

Linked polyamides bind in the minor groove of double-stranded DNA in a partially sequence-specific manner. This report analyzes the theoretical limits of DNA sequence discrimination by linked polyamides composed of two to four different types of heterocyclic rings, determining (i) the optimal choice of base-binding specificity for each ring and (ii) the optimal design for a polyamide composed of these rings to target a given DNA sequence and designed to maximize the fraction of the total polyamide binding to the specified target sequence relative to all other sequences. The results show that, fortuitously, polyamides composed of pyrrole, a naturally occurring G-excluding element, and imidazole, a rationally designed G-favoring element, have features similar to the theoretical optimum design for polyamides composed of two different rings. The results also show that, in polyamides composed of two or three types of heterocyclic rings, choosing a nonspecific “placeholder” ring, which binds equally strongly to each of the four bases, along with one or two base-specific rings will often enhance sequence specificity over a polyamide composed entirely of base-specific rings.

Transgene organization in rice engineered through direct DNA transfer supports a two-phase integration mechanism mediated by the establishment of integration hot spots

Organization of transgenes in rice transformed through direct DNA transfer strongly suggests a two-phase integration mechanism. In the “preintegration” phase, transforming plasmid molecules (either intact or partial) are spliced together. This gives rise to rearranged transgenic sequences, which upon integration do not contain any interspersed plant genomic sequences. Subsequently, integration of transgenic DNA into the host genome is initiated. Our experiments suggest that the original site of integration acts as a hot spot, facilitating subsequent integration of successive transgenic molecules at the same locus. The resulting transgenic locus may have plant DNA separating the transgenic sequences. Our data indicate that transformation through direct DNA transfer, specifically particle bombardment, generally results in a single transgenic locus as a result of this two-phase integration mechanism. Transgenic plants generated through such processes may, therefore, be more amenable to breeding programs as the single transgenic locus will be easier to characterize genetically. Results from direct DNA transfer experiments suggest that in the absence of protein factors involved in exogenous DNA transfer through Agrobacterium, the qualitative and/or quantitative efficiency of transformation events is not compromised. Our results cast doubt on the role of Agrobacterium vir genes in the integration process.

Induction of microsatellite instability by oxidative DNA damage

Instability of repetitive sequences, both in intronic sequences and within coding regions, has been demonstrated to be a hallmark of genomic instability in human cancer. Understanding how these mutational events arise may provide an opportunity for prevention or early intervention in cancer development. To study the source of this instability, we have identified a region of the β-lactamase gene that is tolerant to the insertion of fragments of exogenous DNA as large as 1,614 bp with minimal loss of enzyme activity, as determined by antibiotic resistance. Fragments inserted out-of-frame render Escherichia coli sensitive to antibiotic, and compensatory frameshift mutations that restore the reading frame of β-lactamase can be selected on the basis of antibiotic resistance. We have utilized this site to insert a synthetic microsatellite sequence within the β-lactamase gene and selected for mutations yielding frameshifts. This assay provides for detection of one frameshift mutation in a background of 106 wild-type sequences. Mismatch repair deficiency increased the observed frameshift frequency ≈300-fold. Exposure of plasmid containing microsatellite sequences to hydrogen peroxide resulted in frameshift mutations that were localized exclusively to the microsatellite sequences, whereas DNA damage by UV or N-methyl-N′-nitro-N-nitrosoguanidine did not result in enhanced mutagenesis. We postulate that in tumor cells, endogenous production of oxygen free radicals may be a major factor in promoting instability of microsatellite sequences. This β-lactamase assay may provide a sensitive methodology for the detection and quantitation of mutations associated with the development of cancer.

A Z-DNA binding domain present in the human editing enzyme, double-stranded RNA adenosine deaminase

Editing of RNA changes the read-out of information from DNA by altering the nucleotide sequence of a transcript. One type of RNA editing found in all metazoans uses double-stranded RNA (dsRNA) as a substrate and results in the deamination of adenosine to give inosine, which is translated as guanosine. Editing thus allows variant proteins to be produced from a single pre-mRNA. A mechanism by which dsRNA substrates form is through pairing of intronic and exonic sequences before the removal of noncoding sequences by splicing. Here we report that the RNA editing enzyme, human dsRNA adenosine deaminase (DRADA1, or ADAR1) contains a domain (Zα) that binds specifically to the left-handed Z-DNA conformation with high affinity (KD = 4 nM). As formation of Z-DNA in vivo occurs 5′ to, or behind, a moving RNA polymerase during transcription, recognition of Z-DNA by DRADA1 provides a plausible mechanism by which DRADA1 can be targeted to a nascent RNA so that editing occurs before splicing. Analysis of sequences related to Zα has allowed identification of motifs common to this class of nucleic acid binding domain.

Activities and response to DNA damage of latent and active sequence-specific DNA binding forms of mouse p53

The mouse p53 protein generated by alternative splicing (p53as) has amino acid substitutions at its C terminus that result in constitutively active sequence-specific DNA binding (active form), whereas p53 protein itself binds inefficiently (latent form) unless activated by C-terminal modification. Exogenous p53as expression activated transcription of reporter plasmids containing p53 binding sequences and inhibited growth of mouse and human cells lacking functional endogenous p53. Inducible p53as in stably transfected p53 null fibroblasts increased p21WAF1/Cip-1/Sdi and decreased bcl-2 protein steady-state levels. Endogenous p53as and p53 proteins differed in response to cellular DNA damage. p53 protein was induced transiently in normal keratinocytes and fibroblasts whereas p53as protein accumulation was sustained in parallel with induction of p21WAF1/Cip-1/Sdi protein and mRNA, in support of p53as transcriptional activity. Endogenous p53 and p53as proteins in epidermal tumor cells responded to DNA damage with different kinetics of nuclear accumulation and efficiencies of binding to a p53 consensus DNA sequence. A model is proposed in which C-terminally distinct p53 protein forms specialize in functions, with latent p53 forms primarily for rapid non-sequence-specific binding to sites of DNA damage and active p53 forms for sustained regulation of transcription and growth.

The active digestion of uniparental chloroplast DNA in a single zygote of Chlamydomonas reinhardtii is revealed by using the optical tweezer

The non-Mendelian inheritance of organelle genes is a phenomenon common to almost all eukaryotes, and in the isogamous alga Chlamydomonas reinhardtii, chloroplast (cp) genes are transmitted from the mating type positive (mt+) parent. In this study, the preferential disappearance of the fluorescent cp nucleoids of the mating type negative (mt−) parent was observed in living young zygotes. To study the change in cpDNA molecules during the preferential disappearance, the cpDNA of mt+ or mt− origin was labeled separately with bacterial aadA gene sequences. Then, a single zygote with or without cp nucleoids was isolated under direct observation by using optical tweezers and investigated by nested PCR analysis of the aadA sequences. This demonstrated that cpDNA molecules are digested completely during the preferential disappearance of mt− cp nucleoids within 10 min, whereas mt+ cpDNA and mitochondrial DNA are protected from the digestion. These results indicate that the non-Mendelian transmission pattern of organelle genes is determined immediately after zygote formation.

An intercalation inhibitor altering the target selectivity of DNA damaging agents: Synthesis of site-specific aflatoxin B1 adducts in a p53 mutational hotspot

Aflatoxin B1 (AFB1) is a potent human carcinogen implicated in the etiology of hepatocellular carcinoma. Upon metabolic activation to the reactive epoxide, AFB1 forms DNA adducts primarily at the N7 position of guanines. To elucidate more fully the molecular mechanism of AFB1-induced mutagenesis, an intercalation inhibitor was designed to probe the effects of intercalation by AFB1 epoxide on its reaction with DNA. DNA duplexes were prepared consisting of a target strand containing multiple potentially reactive guanines and a nontarget strand containing a cis-syn thymidine-benzofuran photoproduct. Because the covalently linked benzofuran moiety physically occupies an intercalation site, we reasoned that such a site would be rendered inaccessible to AFB1 epoxide. By strategic positioning of this intercalation inhibitor in the intercalation site 5′ to a specific guanine, the adduct yield at that site was greatly diminished, indicating that intercalation by AFB1 epoxide contributes favorably to adduct formation. Using this approach it has been possible to simplify the production of site-specifically modified oligonucleotides containing AFB1 adducts in the sequence context of a p53 mutational hotspot. Moreover, we report herein isolation of site-specifically AFB1-modified oligonucleotides in sequences containing multiple guanines. Use of intercalation inhibitors will facilitate both investigation of the ability of other carcinogens to intercalate into DNA and the synthesis of specific carcinogen-DNA adducts.

Dual roles for DNA sequence identity and the mismatch repair system in the regulation of mitotic crossing-over in yeast

Sequence divergence acts as a potent barrier to homologous recombination; much of this barrier derives from an antirecombination activity exerted by mismatch repair proteins. An inverted repeat assay system with recombination substrates ranging in identity from 74% to 100% has been used to define the relationship between sequence divergence and the rate of mitotic crossing-over in yeast. To elucidate the role of the mismatch repair machinery in regulating recombination between mismatched substrates, we performed experiments in both wild-type and mismatch repair defective strains. We find that a single mismatch is sufficient to inhibit recombination between otherwise identical sequences, and that this inhibition is dependent on the mismatch repair system. Additional mismatches have a cumulative negative effect on the recombination rate. With sequence divergence of up to approximately 10%, the inhibitory effect of mismatches results mainly from antirecombination activity of the mismatch repair system. With greater levels of divergence, recombination is inefficient even in the absence of mismatch repair activity. In both wild-type and mismatch repair defective strains, an approximate log-linear relationship is observed between the recombination rate and the level of sequence divergence.

Molecular keys to speciation: DNA polymorphism and the control of genetic exchange in enterobacteria

Speciation involves the establishment of genetic barriers between closely related organisms. The extent of genetic recombination is a key determinant and a measure of genetic isolation. The results reported here reveal that genetic barriers can be established, eliminated, or modified by manipulating two systems which control genetic recombination, SOS and mismatch repair. The extent of genetic isolation between enterobacteria is a simple mathematical function of DNA sequence divergence. The function does not depend on hybrid DNA stability, but rather on the number of blocks of sequences identical in the two mating partners and sufficiently large to allow the initiation of recombination. Further, there is no obvious discontinuity in the function that could be used to define a level of divergence for distinguishing species.

Strained DNA is kinked by low concentrations of Zn2+

A novel atomic force microscope with a magnetically oscillated tip has provided unprecedented resolution of small DNA fragments spontaneously adsorbed to mica and imaged in situ in the presence of divalent ions. Kinks (localized bends of average angle 78°) were observed in axially strained minicircles consisting of tandemly repeated d(A)5 and d(GGGCC[C]) sequences. The frequency of kinks in identical minicircles increased 4-fold in the presence of 1 mM Zn2+ compared with 1 mM Mg2+. Kinking persisted in mixed Mg2+/Zn2+ electrolytes until the Zn2+ concentration dropped below 100 μM, indicating that this type of kinking may occur under physiological conditions. Kinking appears to replace intrinsic bending, and statistical analysis shows that kinks are not localized within any single sequence element. A surprisingly small free energy is associated with kink formation.

DNA polymerase δ isolated from Schizosaccharomyces pombe contains five subunits

DNA polymerase δ (pol δ) plays an essential role in DNA replication, repair, and recombination. We have purified pol δ from Schizosaccharomyces pombe more than 103-fold and demonstrated that the polymerase activity of purified S. pombe pol δ is completely dependent on proliferating cell nuclear antigen and replication factor C. SDS/PAGE analysis of the purified fraction indicated that the pol δ complex consists of five subunits that migrate with apparent molecular masses of 125, 55, 54, 42, and 22 kDa. Western blot analysis indicated that the 125, 55, and 54 kDa proteins are the large catalytic subunit (Pol3), Cdc1, and Cdc27, respectively. The identity of the other two subunits, p42 and p22, was determined following proteolytic digestion and sequence analysis of the resulting peptides. The peptide sequences derived from the p22 subunit indicated that this subunit is identical to Cdm1, previously identified as a multicopy suppressor of the temperature-sensitive cdc1-P13 mutant, whereas peptide sequences derived from the p42 subunit were identical to a previously uncharacterized ORF located on S. pombe chromosome 1.