Determination of the base composition of deoxyribonucleic acid from its buoyant density in CsCl

Structural Capsidomics of Single-Stranded DNA Viruses

Single-stranded DNA (ssDNA) viruses are a diverse group of pathogens with broad host range, including bacteria, archaea, protists, fungi, plants, invertebrates, and vertebrates. Their small compact genomes have evolved to encode multiple proteins. This review focuses on the structure and functional diversity of the icosahedral capsids across the ssDNA viruses. To date, X-ray crystallography and cryo-electron microscopy structural studies have provided detailed capsid architectures for 8 of the 35 ssDNA virus families, illustrating variations in assembly mechanisms, symmetry, and structural adaptations of the capsid. However, common features include the conserved jelly-roll motif of the capsid protein and strategies for genome packaging, also showing evolutionary convergence. The ever-increasing availability of genomic sequences of ssDNA viruses and predictive protein modeling programs, such as using AlphaFold, allows for the extension of structural insights to the less-characterized families. Therefore, this review is a comparative analysis of the icosahedral ssDNA virus families and how the capsid proteins are arranged with different tessellations to form icosahedral spheres. It summarizes the current knowledge, emphasizing gaps in the structural characterization of the ssDNA capsidome, and it underscores the importance of continued exploration to understand the molecular underpinnings of capsid function and evolution. These insights have implications for virology, molecular biology, and therapeutic applications.

Charging of DNA Complexes in Positive-Mode Native Electrospray Ionization Mass Spectrometry

Native mass spectrometry (nMS) provides insights into the structures and dynamics of biomacromolecules in their native-like states by preserving noncovalent interactions through "soft" electrospray ionization (ESI). For native proteins, the number of charges that are acquired scales with the surface area and mass. Here, we explore the effect of highly negatively charged DNA on the ESI charge of protein complexes and find a reduction of the mass-to-charge ratio as well as a greater variation. The charge state distributions of pure DNA assemblies show a lower mass-to-charge ratio than proteins due to their greater density in the gas phase, whereas the charge of protein-DNA complexes can additionally be influenced by the distribution of the ESI charges, ion pairing events, and collapse of the DNA components. Our findings suggest that structural features of protein-DNA complexes can result in lower charge states than expected for proteins.

Establishing Host-Virus Link Through Host Metabolism: Viral DNA SIP Validation Using T4 Bacteriophage and E. coli

DNA Stable Isotope Probing is emerging as a potent methodology for investigating host-virus interactions, based on the essential reliance of viruses on host organisms for the production of virions. Despite the anticipated link between host isotopic compositions and the generated virions, the application of stable isotope probing to viral DNA has never been evaluated on simple biological models. In this study, we assessed the efficacy of this method on the bacteriophage T4 and its host, Escherichia coli. Through the cultivation of E. coli cells on a 13C-enriched substrate and subsequent propagation of T4 bacteriophage, we examine the degree of isotopic enrichment in viral DNA. Our investigation reveals a strong correlation between the proportion of 13C6-D-glucose in the growth substrate and the buoyant density in CsCl gradient of T4 DNA, confirming the validity of DNA SIP in viral ecology. These findings underscore the potential of DNA SIP as a robust tool for characterizing the diversity of viruses infecting hosts with specific metabolic activities and provide then a foundation for further exploration in viral ecology research.

Temperature-induced DNA density transition in phage λ capsid revealed with contrast-matching SANS

Structural details of a genome packaged in a viral capsid are essential for understanding how the structural arrangement of a viral genome in a capsid controls its release dynamics during infection, which critically affects viral replication. We previously found a temperature-induced, solid-like to fluid-like mechanical transition of packaged λ-genome that leads to rapid DNA ejection. However, an understanding of the structural origin of this transition was lacking. Here, we use small-angle neutron scattering (SANS) to reveal the scattering form factor of dsDNA packaged in phage λ capsid by contrast matching the scattering signal from the viral capsid with deuterated buffer. We used small-angle X-ray scattering and cryoelectron microscopy reconstructions to determine the initial structural input parameters for intracapsid DNA, which allows accurate modeling of our SANS data. As result, we show a temperature-dependent density transition of intracapsid DNA occurring between two coexisting phases-a hexagonally ordered high-density DNA phase in the capsid periphery and a low-density, less-ordered DNA phase in the core. As the temperature is increased from 20 °C to 40 °C, we found that the core-DNA phase undergoes a density and volume transition close to the physiological temperature of infection (~37 °C). The transition yields a lower energy state of DNA in the capsid core due to lower density and reduced packing defects. This increases DNA mobility, which is required to initiate rapid genome ejection from the virus capsid into a host cell, causing infection. These data reconcile our earlier findings of mechanical DNA transition in phage.

Genetic and phenotypic assessments for the safety of probiotic Bacillus clausii 088AE

In this study, we report on whole genome sequence analysis of clinically documented, commercial probiotic Bacillus clausii 088AE and genome features contributing to probiotic properties. The whole genome sequence of B. clausii 088AE generated a single scaffold of 4,598,457 bp with 44.74 mol% G + C. This assembled genome sequence annotated by the RAST resulted in 4371 coding genes, 75 tRNAs, and 22 rRNAs. Gene ontology classification indicated 39.5% proteins with molecular function, 44.24% cellular component, and 16.25% proteins involved in biological processes. In taxonomic analysis, B. clausii 088AE shared 99% identity with B. clausii DSM 8716. The gene sequences related to safety and genome stability such as antibiotic resistance (840), virulence factors (706), biogenic amines (1), enterotoxin (0), emetic toxin (0), lanthipeptides (4), prophage (4) and clustered regularly interspaced short palindromic repeats (CRISPR) sequences (11), were identified and evaluated for safety and functions. The absence of functional prophage sequences and the presence of CRISPR indicated an advantage in genome stability. Moreover, the presence of genome features contributing to probiotic characteristics such as acid, and bile salt tolerance, adhesion to the gut mucosa, and environmental resistance ensure the strains survivability when consumed as a probiotic. In conclusion, the absence of risks associated with sequences/genes in the B. clausii 088AE genome and the presence of essential probiotic traits confirm the strain to be safe for use as a probiotic.

Plant-associated fungi support bacterial resilience following water limitation

Drought disrupts soil microbial activity and many biogeochemical processes. Although plant-associated fungi can support plant performance and nutrient cycling during drought, their effects on nearby drought-exposed soil microbial communities are not well resolved. We used H218O quantitative stable isotope probing (qSIP) and 16S rRNA gene profiling to investigate bacterial community dynamics following water limitation in the hyphospheres of two distinct fungal lineages (Rhizophagus irregularis and Serendipita bescii) grown with the bioenergy model grass Panicum hallii. In uninoculated soil, a history of water limitation resulted in significantly lower bacterial growth potential and growth efficiency, as well as lower diversity in the actively growing bacterial community. In contrast, both fungal lineages had a protective effect on hyphosphere bacterial communities exposed to water limitation: bacterial growth potential, growth efficiency, and the diversity of the actively growing bacterial community were not suppressed by a history of water limitation in soils inoculated with either fungus. Despite their similar effects at the community level, the two fungal lineages did elicit different taxon-specific responses, and bacterial growth potential was greater in R. irregularis compared to S. bescii-inoculated soils. Several of the bacterial taxa that responded positively to fungal inocula belong to lineages that are considered drought susceptible. Overall, H218O qSIP highlighted treatment effects on bacterial community structure that were less pronounced using traditional 16S rRNA gene profiling. Together, these results indicate that fungal-bacterial synergies may support bacterial resilience to moisture limitation.

Ectomycorrhizal fungi mediate belowground carbon transfer between pines and oaks

Inter-kingdom belowground carbon (C) transfer is a significant, yet hidden, biological phenomenon, due to the complexity and highly dynamic nature of soil ecology. Among key biotic agents influencing C allocation belowground are ectomycorrhizal fungi (EMF). EMF symbiosis can extend beyond the single tree-fungus partnership to form common mycorrhizal networks (CMNs). Despite the high prevalence of CMNs in forests, little is known about the identity of the EMF transferring the C and how these in turn affect the dynamics of C transfer. Here, Pinus halepensis and Quercus calliprinos saplings growing in forest soil were labeled using a 13CO2 labeling system. Repeated samplings were applied during 36 days to trace how 13C was distributed along the tree-fungus-tree pathway. To identify the fungal species active in the transfer, mycorrhizal fine root tips were used for DNA-stable isotope probing (SIP) with 13CO2 followed by sequencing of labeled DNA. Assimilated 13CO2 reached tree roots within four days and was then transferred to various EMF species. C was transferred across all four tree species combinations. While Tomentella ellisii was the primary fungal mediator between pines and oaks, Terfezia pini, Pustularia spp., and Tuber oligospermum controlled C transfer among pines. We demonstrate at a high temporal, quantitative, and taxonomic resolution, that C from EMF host trees moved into EMF and that C was transferred further to neighboring trees of similar and distinct phylogenies.

Establishment of a Method for Investigating Direct and Indirect Actions of Ionizing Radiation Using Scavenger-free Plasmid DNA

In this study, an improved method using scavenger-free plasmid DNA was established to accurately determine yields of DNA damage induced by direct and indirect actions of ionizing radiation. The scavenger-free plasmid DNA was obtained by dialysis over 5-7 days, and the DNA solvent was replaced with phosphate buffer to completely remove impurities, which could be scavengers of radicals produced as a result of water radiolysis. DNA samples of films and dilute aqueous solutions were used to separately evaluate contributions of the direct and indirect actions of X rays (150-160 kVp). The irradiated DNA was analyzed by agarose gel electrophoresis to quantify strand-break yields. The yields of single-strand breaks (SSBs), n(SSB), were determined to be (6.5 ± 2.0) × 10-10 and (3.1 ± 0.9) × 10-7 SSBs/Gy/Da for the film and solution samples, respectively, showing a significant contribution of hydroxyl radicals (•OH) compared with direct energy depositions from ionizing radiation to DNA. As observed in SSBs, the yields of double-strand breaks (DSBs), n(DSB), were (5.6 ± 1.1) × 10-11 and (1.3 ± 0.2) × 10-8 DSBs/Gy/Da for the film and solution samples, respectively. The yield ratio of DSBs to SSBs, that is, n(DSB)/n(SSB), was 0.091 ± 0.026 for the film samples, while it was much lower for the solution samples (0.045 ± 0.010), indicating that direct actions result in more localized strand breaks relative to indirect actions. Base excision repair enzymes, namely, endonuclease III (Nth) and formamidopyrimidine-DNA glycosylase (Fpg), were utilized after irradiations to convert base lesions and apurinic/apyrimidinic (AP) sites into strand breaks. The amounts of Nth and Fpg for the conversion were optimized to a few units per μg of DNA, although the optimal concentrations can differ among conditions.

Developing an Updated Strategy for Estimating the Free-Energy Parameters in RNA Duplexes

For the last 20 years, it has been common lore that the free energy of RNA duplexes formed from canonical Watson-Crick base pairs (bps) can be largely approximated with dinucleotide bp parameters and a few simple corrective constants that are duplex independent. Additionally, the standard benchmark set of duplexes used to generate the parameters were GC-rich in the shorter duplexes and AU-rich in the longer duplexes, and the length of the majority of the duplexes ranged between 6 and 8 bps. We were curious if other models would generate similar results and whether adding longer duplexes of 17 bps would affect the conclusions. We developed a gradient-descent fitting program for obtaining free-energy parameters-the changes in Gibbs free energy (ΔG), enthalpy (ΔH), and entropy (ΔS), and the melting temperature (Tm)-directly from the experimental melting curves. Using gradient descent and a genetic algorithm, the duplex melting results were combined with the standard benchmark data to obtain bp parameters. Both the standard (Turner) model and a new model that includes length-dependent terms were tested. Both models could fit the standard benchmark data; however, the new model could handle longer sequences better. We developed an updated strategy for fitting the duplex melting data.

Direct observation of ion emission from charged aqueous nanodrops: effects on gaseous macromolecular charging

Mechanistic information about how gaseous ions are formed from charged droplets has been difficult to establish because direct observation of nanodrops in a size range relevant to gaseous macromolecular ion formation by optical or traditional mass spectrometry methods is challenging owing to their small size and heterogeneity. Here, the mass and charge of individual aqueous nanodrops between 1-10 MDa (15-32 nm diameter) with ∼50-300 charges are dynamically monitored for 1 s using charge detection mass spectrometry. Discrete losses of minimally solvated singly charged ions occur, marking the first direct observation of ion emission from aqueous nanodrops in late stages of droplet evaporation relevant to macromolecular ion formation in native mass spectrometry. Nanodrop charge depends on the identity of constituent ions, with pure water nanodrops charged slightly above the Rayleigh limit and aqueous solutions containing alkali metal ions charged progressively below the Rayleigh limit with increasing cation size. MS2 capsid ions (∼3.5 MDa; ∼27 nm diameter) are more highly charged from aqueous ammonium acetate than from its biochemically preferred, 100 mM NaCl/10 mM Na phosphate solution, consistent with ion emission reducing the nanodrop and resulting capsid charge. The extent of charging indicates that the capsid partially collapses inside the nanodrops prior to the charging and formation of the dehydrated gaseous ions. These results demonstrate that ion emission can affect macromolecular charging and that conformational changes to macromolecular structure can occur in nanodrops prior to the formation of naked gaseous ions.

Measurement Error and Resolution in Quantitative Stable Isotope Probing: Implications for Experimental Design

Quantitative stable isotope probing (qSIP) estimates isotope tracer incorporation into DNA of individual microbes and can link microbial biodiversity and biogeochemistry in complex communities. As with any quantitative estimation technique, qSIP involves measurement error, and a fuller understanding of error, precision, and statistical power benefits qSIP experimental design and data interpretation. We used several qSIP data sets-from soil and seawater microbiomes-to evaluate how variance in isotope incorporation estimates depends on organism abundance and resolution of the density fractionation scheme. We assessed statistical power for replicated qSIP studies, plus sensitivity and specificity for unreplicated designs. As a taxon's abundance increases, the variance of its weighted mean density declines. Nine fractions appear to be a reasonable trade-off between cost and precision for most qSIP applications. Increasing the number of density fractions beyond that reduces variance, although the magnitude of this benefit declines with additional fractions. Our analysis suggests that, if a taxon has an isotope enrichment of 10 atom% excess, there is a 60% chance that this will be detected as significantly different from zero (with alpha 0.1). With five replicates, isotope enrichment of 5 atom% could be detected with power (0.6) and alpha (0.1). Finally, we illustrate the importance of internal standards, which can help to calibrate per sample conversions of %GC to mean weighted density. These results should benefit researchers designing future SIP experiments and provide a useful reference for metagenomic SIP applications where both financial and computational limitations constrain experimental scope.IMPORTANCE One of the biggest challenges in microbial ecology is correlating the identity of microorganisms with the roles they fulfill in natural environmental systems. Studies of microbes in pure culture reveal much about their genomic content and potential functions but may not reflect an organism's activity within its natural community. Culture-independent studies supply a community-wide view of composition and function in the context of community interactions but often fail to link the two. Quantitative stable isotope probing (qSIP) is a method that can link the identity and functional activity of specific microbes within a naturally occurring community. Here, we explore how the resolution of density gradient fractionation affects the error and precision of qSIP results, how they may be improved via additional experimental replication, and discuss cost-benefit balanced scenarios for SIP experimental design.

Microbial Taxon-Specific Isotope Incorporation with DNA Quantitative Stable Isotope Probing

Quantitative stable isotope probing (qSIP) measures rates of taxon-specific element assimilation in intact microbial communities, utilizing substrates labeled with a heavy isotope.The laboratory protocol for qSIP is nearly identical to that for conventional stable isotope probing, with two key additions: (1) in qSIP, qPCR measurements are conducted on each density fraction recovered after isopycnic separation, and (2) in qSIP, multiple density fractions are sequenced spanning the entire range of densities over which nucleic acids were recovered. qSIP goes beyond identifying taxa assimilating a substrate, as it also allows for measuring that assimilation for each taxon within a given microbial community. Here, we describe an analysis process necessary to determine atom fraction excess of a heavy stable isotope added to an environmental sample for a given taxon's DNA.

Taxon-specific microbial growth and mortality patterns reveal distinct temporal population responses to rewetting in a California grassland soil

Microbial activity increases after rewetting dry soil, resulting in a pulse of carbon mineralization and nutrient availability. The biogeochemical responses to wet-up are reasonably well understood and known to be microbially mediated. Yet, the population level dynamics, and the resulting changes in microbial community patterns, are not well understood as ecological phenomena. Here, we used sequencing of 16S rRNA genes coupled with heavy water (H218O) DNA quantitative stable isotope probing to estimate population-specific rates of growth and mortality in response to a simulated wet-up event in a California annual grassland soil. Bacterial growth and mortality responded rapidly to wet-up, within 3 h, and continued throughout the 168 h incubation, with patterns of sequential growth observed at the phylum level. Of the 37 phyla detected in the prewet community, growth was found in 18 phyla while mortality was measured in 26 phyla. Rapid growth and mortality rates were measurable within 3 h of wet-up but had contrasting characteristics; growth at 3 h was dominated by select taxa within the Proteobacteria and Firmicutes, whereas mortality was taxonomically widespread. Furthermore, across the community, mortality exhibited density-independence, consistent with the indiscriminate shock resulting from dry-down and wet-up, whereas growth was density-dependent, consistent with control by competition or predation. Total aggregated growth across the community was highly correlated with total soil CO2 production. Together, these results illustrate how previously "invisible" population responses can translate quantitatively to emergent observations of ecosystem-scale biogeochemistry.

Optimizing ultracentrifugation conditions for DNA-based stable isotope probing (DNA-SIP)

DNA-SIP (DNA-based stable isotope probing) is increasingly being employed in soil microbial ecology to identify those microbes assimilating the ¹³C/¹⁵N labelled substrate. Isopycnic gradient centrifugation is the primary experimental process for conducting DNA-SIP. However, diverse centrifugal conditions have been used in various recent studies. In order to get the optimum conditions of centrifugation for DNA-SIP, centrifugation time (36, 42, 48, 60 h), speed (45,000, 55,000 rpm) and the initial buoyant density (1.69, 1.71, 1.725 g ml^-1), as were used extensively in related studies, were tested in this experiment with the Vti 65.2 rotor. DNA with either ¹³C-labelling or unlabelled was extracted from a paddy soil pre-incubated with either ¹³C-labelled or natural abundance glucose. After ultracentrifugation, the gene abundance of bacterial 16S rRNA, fungal 18S rRNA, bacterial and archaeal amoA within the fractioned DNA was detected. The results showed that centrifugation for 48 h was enough for the DNA to reach stabilization in the CsCl solution. The initial density of the mixed solution was best adjusted to 1.71 g ml^-1 to ensure that most of the genes were concentrated on the middle fractions of the density gradient. Increasing the centrifugation speed would increase the density gradient of fractions; therefore, 45,000 rpm (184,000 g) was recommended so as to obtain the more widespread pattern of DNA in the centrifugal tube. We hope these findings will assist future researchers to conduct optimum ultracentrifugation for DNA-SIP.

Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding

BACKGROUND

In late December, 2019, patients presenting with viral pneumonia due to an unidentified microbial agent were reported in Wuhan, China. A novel coronavirus was subsequently identified as the causative pathogen, provisionally named 2019 novel coronavirus (2019-nCoV). As of Jan 26, 2020, more than 2000 cases of 2019-nCoV infection have been confirmed, most of which involved people living in or visiting Wuhan, and human-to-human transmission has been confirmed.

METHODS

We did next-generation sequencing of samples from bronchoalveolar lavage fluid and cultured isolates from nine inpatients, eight of whom had visited the Huanan seafood market in Wuhan. Complete and partial 2019-nCoV genome sequences were obtained from these individuals. Viral contigs were connected using Sanger sequencing to obtain the full-length genomes, with the terminal regions determined by rapid amplification of cDNA ends. Phylogenetic analysis of these 2019-nCoV genomes and those of other coronaviruses was used to determine the evolutionary history of the virus and help infer its likely origin. Homology modelling was done to explore the likely receptor-binding properties of the virus.

FINDINGS

The ten genome sequences of 2019-nCoV obtained from the nine patients were extremely similar, exhibiting more than 99·98% sequence identity. Notably, 2019-nCoV was closely related (with 88% identity) to two bat-derived severe acute respiratory syndrome (SARS)-like coronaviruses, bat-SL-CoVZC45 and bat-SL-CoVZXC21, collected in 2018 in Zhoushan, eastern China, but were more distant from SARS-CoV (about 79%) and MERS-CoV (about 50%). Phylogenetic analysis revealed that 2019-nCoV fell within the subgenus Sarbecovirus of the genus Betacoronavirus, with a relatively long branch length to its closest relatives bat-SL-CoVZC45 and bat-SL-CoVZXC21, and was genetically distinct from SARS-CoV. Notably, homology modelling revealed that 2019-nCoV had a similar receptor-binding domain structure to that of SARS-CoV, despite amino acid variation at some key residues.

INTERPRETATION

2019-nCoV is sufficiently divergent from SARS-CoV to be considered a new human-infecting betacoronavirus. Although our phylogenetic analysis suggests that bats might be the original host of this virus, an animal sold at the seafood market in Wuhan might represent an intermediate host facilitating the emergence of the virus in humans. Importantly, structural analysis suggests that 2019-nCoV might be able to bind to the angiotensin-converting enzyme 2 receptor in humans. The future evolution, adaptation, and spread of this virus warrant urgent investigation.

FUNDING

National Key Research and Development Program of China, National Major Project for Control and Prevention of Infectious Disease in China, Chinese Academy of Sciences, Shandong First Medical University.

Simulating metagenomic stable isotope probing datasets with MetaSIPSim

Background

DNA-stable isotope probing (DNA-SIP) links microorganisms to their in-situ function in diverse environmental samples. Combining DNA-SIP and metagenomics (metagenomic-SIP) allows us to link genomes from complex communities to their specific functions and improves the assembly and binning of these targeted genomes. However, empirical development of metagenomic-SIP methods is hindered by the complexity and cost of these studies. We developed a toolkit, ‘MetaSIPSim,’ to simulate sequencing read libraries for metagenomic-SIP experiments. MetaSIPSim is intended to generate datasets for method development and testing. To this end, we used MetaSIPSim generated data to demonstrate the advantages of metagenomic-SIP over a conventional shotgun metagenomic sequencing experiment.

Results

Through simulation we show that metagenomic-SIP improves the assembly and binning of isotopically labeled genomes relative to a conventional metagenomic approach. Improvements were dependent on experimental parameters and on sequencing depth. Community level G + C content impacted the assembly of labeled genomes and subsequent binning, where high community G + C generally reduced the benefits of metagenomic-SIP. Furthermore, when a high proportion of the community is isotopically labeled, the benefits of metagenomic-SIP decline. Finally, the choice of gradient fractions to sequence greatly influences method performance.

Conclusions

Metagenomic-SIP is a valuable method for recovering isotopically labeled genomes from complex communities. We show that metagenomic-SIP performance depends on optimization of experimental parameters. MetaSIPSim allows for simulation of metagenomic-SIP datasets which facilitates the optimization and development of metagenomic-SIP experiments and analytical approaches for dealing with these data.

Genomics of the Pathogenic Clostridia

Whole-genome sequences are now available for all the clinically important clostridia and many of the lesser or opportunistically pathogenic clostridia. The complex clade structures of C. difficile, C. perfringens, and the species that produce botulinum toxins have been delineated by whole-genome sequence analysis. The true clostridia of cluster I show relatively low levels of gross genomic rearrangements within species, in contrast to the species of cluster XI, notably C. difficile, which have been found to have very plastic genomes with significant levels of chromosomal rearrangement. Throughout the clostridial phylotypes, a large proportion of the strain diversity is driven by the acquisition and loss of mobile elements, including phages, plasmids, insertion sequences, and transposons. Genomic analysis has been used to investigate the diversity and spread of C. difficile within hospital settings, the zoonotic transfer of isolates, and the emergence, origins, and geographic spread of epidemic ribotypes. In C. perfringens the clades defined by chromosomal sequence analysis show no indications of clustering based on host species or geographical location. Whole-genome sequence analysis helps to define the different survival and pathogenesis strategies that the clostridia use. Some, such as C. botulinum, produce toxins which rapidly act to kill the host, whereas others, such as C. perfringens and C. difficile, produce less lethal toxins which can damage tissue but do not rapidly kill the host. The genomes provide a resource that can be mined to identify potential vaccine antigens and targets for other forms of therapeutic intervention.

Thermodynamics of force-induced B-DNA melting: Single-strand discreteness matters

Overstretching of B-DNA is currently understood as force-induced melting. Depending on the geometry of the stretching experiment, the force threshold for the overstretching transition is around 65 or 110 pN. Although the mechanisms behind force-induced melting have been correctly described by Rouzina and Bloomfield [Biophys. J. 80, 882 (2001)BIOJAU0006-349510.1016/S0006-3495(01)76067-5], neither force threshold has been exactly calculated by theory. In this work, a detailed analysis of the force-extension curve is presented, based on a description of single-stranded (ss) DNA in terms of the discrete Kratky-Porod model, consistent with (i) the contour length expected from the crystallographically determined monomer distance and (ii) a high value of the elastic stretch modulus arising from covalent bonding. The value estimated for the ss-DNA persistence length, λ=1.0 nm, is at the low end of currently known estimates and reflects the intrinsic stiffness of the partially, or fully stretched state, where electrostatic repulsion effects are expected to be minimal. A detailed analysis of single- and double-stranded DNA free energies provides estimates of the overstretching force thresholds. In the unconstrained geometry, the predicted threshold is 64 pN. In the constrained geometry, after allowing for the entropic penalty of the plectonemic topology of the molten state, the predicted threshold is 111 pN.

On the length, weight and GC content of the human genome

Objective

Basic parameters commonly used to describe genomes including length, weight and relative guanine-cytosine (GC) content are widely cited in absence of a primary source. By using updated data and original software we determined these values to the best of our knowledge as standard reference for the whole human nuclear genome, for each chromosome and for mitochondrial DNA. We also devised a method to calculate the relative GC content in the whole messenger RNA sequence set and in transcriptomes by multiplying the GC content of each gene by its mean expression level.

Results

The male nuclear diploid genome extends for 6.27 Gigabase pairs (Gbp), is 205.00 cm (cm) long and weighs 6.41 picograms (pg). Female values are 6.37 Gbp, 208.23 cm, 6.51 pg. The individual variability and the implication for the DNA informational density in terms of bits/volume were discussed. The genomic GC content is 40.9%. Following analysis in different transcriptomes and species, we showed that the greatest deviation was observed in the pathological condition analysed (trisomy 21 leukaemic cells) and in Caenorhabditis elegans. Our results may represent a solid basis for further investigation on human structural and functional genomics while also providing a framework for other genome comparative analysis.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4137-z) contains supplementary material, which is available to authorized users.

DNA-Based Stable Isotope Probing

Microbiomes on Earth are often considered the most heterogeneous biological entities, but their vital roles in driving global biogeochemical cycles often remain elusive. DNA-based stable isotope probing (DNA-SIP) provides a powerful means to establish a direct link between biogeochemical processes and the taxonomic identities of active microorganisms involved in the processes. Combined with high-throughput sequencing, it significantly aids in deciphering ecophysiological functions of active microorganisms at the level of microbial communities. DNA-SIP relies solely on the propagation of targeted microbial communities, during which the entire genomes of daughter cells are synthesized and increasingly ¹³C-labeled. This growth on ¹³C-labeled substrate in association with cell division provides solid evidence for the functional importance and metabolic potential of targeted microorganisms. The essential prerequisite for a successful DNA-SIP experiment is the identification, with confidence, of isotopically enriched ¹³C-DNA, of which the amount is generally too low to allow for the direct measurement of ¹³C atomic percent of nucleic acid. The ¹³C labeling can be readily identified in the fractionated DNA by quantification of functional genes specific to the known targeted microorganisms, and by high-throughput sequencing of the total microbial communities via 16S rRNA genes without prior knowledge of which microorganisms are ¹³C-labeled (i.e., highly enriched in the heavy fractions relative to ¹²C (natural isotope abundance) control treatments). In this chapter, the protocol for obtaining DNA highly enriched in heavy isotope is presented using diazotrophic methanotrophs in a paddy soil as a case study.

Quantifying local randomness in human DNA and RNA sequences using Erdös motifs

In 1932, Paul Erdös asked whether a random walk constructed from a binary sequence can achieve the lowest possible deviation (lowest discrepancy), for the sequence itself and for all its subsequences formed by homogeneous arithmetic progressions. Although avoiding low discrepancy is impossible for infinite sequences, as recently proven by Terence Tao, attempts were made to construct such sequences with finite lengths. We recognize that such constructed sequences (we call these "Erdös sequences") exhibit certain hallmarks of randomness at the local level: they show roughly equal frequencies of short subsequences, and at the same time exclude trivial periodic patterns. For the human DNA we examine the frequency of a set of Erdös motifs of length-10 using three nucleotides-to-binary mappings. The particular length-10 Erdös sequence is derived from the length-11 Mathias sequence and is identical with the first 10 digits of the Thue-Morse sequence, underscoring the fact that both are deficient in periodicities. Our calculations indicate that: (1) the purine(A and G)/pyridimine(C and T) based Erdös motifs are greatly underrepresented in the human genome, (2) the strong(G and C)/weak(A and T) based Erdös motifs are slightly overrepresented, (3) the densities of the two are negatively correlated, (4) the Erdös motifs based on all three mappings being combined are slightly underrepresented, and (5) the strong/weak based Erdös motifs are greatly overrepresented in the human messenger RNA sequences.

Stable isotope informed genome-resolved metagenomics reveals that Saccharibacteria utilize microbially-processed plant-derived carbon

BACKGROUND

The transformation of plant photosynthate into soil organic carbon and its recycling to CO₂ by soil microorganisms is one of the central components of the terrestrial carbon cycle. There are currently large knowledge gaps related to which soil-associated microorganisms take up plant carbon in the rhizosphere and the fate of that carbon.

RESULTS

We conducted an experiment in which common wild oats (Avena fatua) were grown in a ¹³CO₂ atmosphere and the rhizosphere and non-rhizosphere soil was sampled for genomic analyses. Density gradient centrifugation of DNA extracted from soil samples enabled distinction of microbes that did and did not incorporate the ¹³C into their DNA. A 1.45-Mbp genome of a Saccharibacteria (TM7) was identified and, despite the microbial complexity of rhizosphere soil, curated to completion. The genome lacks many biosynthetic pathways, including genes required to synthesize DNA de novo. Rather, it requires externally derived nucleotides for DNA and RNA synthesis. Given this, we conclude that rhizosphere-associated Saccharibacteria recycle DNA from bacteria that live off plant exudates and/or phage that acquired ¹³C because they preyed upon these bacteria and/or directly from the labeled plant DNA. Isotopic labeling indicates that the population was replicating during the 6-week period of plant growth. Interestingly, the genome is ~ 30% larger than other complete Saccharibacteria genomes from non-soil environments, largely due to more genes for complex carbon utilization and amino acid metabolism. Given the ability to degrade cellulose, hemicellulose, pectin, starch, and 1,3-β-glucan, we predict that this Saccharibacteria generates energy by fermentation of soil necromass and plant root exudates to acetate and lactate. The genome also encodes a linear electron transport chain featuring a terminal oxidase, suggesting that this Saccharibacteria may respire aerobically. The genome encodes a hydrolase that could breakdown salicylic acid, a plant defense signaling molecule, and genes to interconvert a variety of isoprenoids, including the plant hormone zeatin.

CONCLUSIONS

Rhizosphere Saccharibacteria likely depend on other bacteria for basic cellular building blocks. We propose that isotopically labeled CO₂ is incorporated into plant-derived carbon and then into the DNA of rhizosphere organisms capable of nucleotide synthesis, and the nucleotides are recycled into Saccharibacterial genomes.

SIPSim: A Modeling Toolkit to Predict Accuracy and Aid Design of DNA-SIP Experiments

DNA Stable isotope probing (DNA-SIP) is a powerful method that links identity to function within microbial communities. The combination of DNA-SIP with multiplexed high throughput DNA sequencing enables simultaneous mapping of in situ assimilation dynamics for thousands of microbial taxonomic units. Hence, high throughput sequencing enabled SIP has enormous potential to reveal patterns of carbon and nitrogen exchange within microbial food webs. There are several different methods for analyzing DNA-SIP data and despite the power of SIP experiments, it remains difficult to comprehensively evaluate method accuracy across a wide range of experimental parameters. We have developed a toolset (SIPSim) that simulates DNA-SIP data, and we use this toolset to systematically evaluate different methods for analyzing DNA-SIP data. Specifically, we employ SIPSim to evaluate the effects that key experimental parameters (e.g., level of isotopic enrichment, number of labeled taxa, relative abundance of labeled taxa, community richness, community evenness, and beta-diversity) have on the specificity, sensitivity, and balanced accuracy (defined as the product of specificity and sensitivity) of DNA-SIP analyses. Furthermore, SIPSim can predict analytical accuracy and power as a function of experimental design and community characteristics, and thus should be of great use in the design and interpretation of DNA-SIP experiments.

QUANTUM AND CONTINUOUS EVOLUTION OF DNA BASE COMPOSITION IN THE YEAST GENUS PICHIA

This paper investigates the noncontinuous nature and evolution of the base composition of nuclear DNA (expressed as mol% guanine + cytosine) in species of the yeast genus Pichia (sensu Kurtzman, 1984b). The pattern of change in the G + C contents in species of this genus, which range from about 27 to 52 mol%, was evaluated. When specifically those species of Pichia were analyzed that have evolved in necroses of cactus species and associated Drosophila, a periodic change in the G + C contents of approximately 3.0-3.2 mol% was detected by a "bootstrapping" method, Fourier analysis, and a nonlinear trigonometric model. Pichia species occurring in exudates of broad-leaved deciduous trees or associated Drosophila and substrates such as soil and water ("other") showed a periodicity of 2.5-2.6 mol%, whereas species associated with conifers and associated bark beetles showed no significant periodicity. Periodicity in the most recent association (cactus and resident Drosophila) as compared to the lack of periodicity in the oldest association (conifer-beetle) may indicate mixed evolutionary processes. Low mol% G + C values appear more frequently in the relatively recent cactus and Drosophila-associated yeast species. In addition, low mol% G + C species do not display the ancestral bud-meiosis mode of sexual reproduction which occurs frequently in medium to high mol% G + C yeasts. It was found that the mol% G + C content of the Drosophila- and cactus-associated Pichia species is positively correlated with the number of compounds fermented or respired by these yeast species. Possible reasons for the periodic changes in mol% G + C content accompanying speciation include aneuploidy, allopolyploidy, the presence of nuclear plasmids, and regular differences in moderately repetitive portions of DNA. Since significant DNA complementarity is virtually limited to species within a relatively narrow G + C group, this suggests that there are at least two processes which alter the G + C content between species, one saltational and one continuous.

Anaerobic carboxydotrophic bacteria in geothermal springs identified using stable isotope probing

Carbon monoxide (CO) is a potential energy and carbon source for thermophilic bacteria in geothermal environments. Geothermal sites ranging in temperature from 45 to 65°C were investigated for the presence and activity of anaerobic CO-oxidizing bacteria. Anaerobic CO oxidation potentials were measured at up to 48.9 μmoles CO g⁻¹ (wet weight) day⁻¹ within five selected sites. Active anaerobic carboxydotrophic bacteria were identified using ¹³CO DNA stable isotope probing (SIP) combined with pyrosequencing of 16S rRNA genes amplified from labeled DNA. Bacterial communities identified in heavy DNA fractions were predominated by Firmicutes, which comprised up to 95% of all sequences in ¹³CO incubations. The predominant bacteria that assimilated ¹³C derived from CO were closely related (>98% 16S rRNA gene sequence identity) to genera of known carboxydotrophs including Thermincola, Desulfotomaculum, Thermolithobacter, and Carboxydocella, although a few species with lower similarity to known bacteria were also found that may represent previously unconfirmed CO-oxidizers. While the distribution was variable, many of the same OTUs were identified across sample sites from different temperature regimes. These results show that bacteria capable of using CO as a carbon source are common in geothermal springs, and that thermophilic carboxydotrophs are probably already quite well known from cultivation studies.

Quantitative microbial ecology through stable isotope probing

Bacteria grow and transform elements at different rates, and as yet, quantifying this variation in the environment is difficult. Determining isotope enrichment with fine taxonomic resolution after exposure to isotope tracers could help, but there are few suitable techniques. We propose a modification to stable isotope probing (SIP) that enables the isotopic composition of DNA from individual bacterial taxa after exposure to isotope tracers to be determined. In our modification, after isopycnic centrifugation, DNA is collected in multiple density fractions, and each fraction is sequenced separately. Taxon-specific density curves are produced for labeled and nonlabeled treatments, from which the shift in density for each individual taxon in response to isotope labeling is calculated. Expressing each taxon's density shift relative to that taxon's density measured without isotope enrichment accounts for the influence of nucleic acid composition on density and isolates the influence of isotope tracer assimilation. The shift in density translates quantitatively to isotopic enrichment. Because this revision to SIP allows quantitative measurements of isotope enrichment, we propose to call it quantitative stable isotope probing (qSIP). We demonstrated qSIP using soil incubations, in which soil bacteria exhibited strong taxonomic variations in (18)O and (13)C composition after exposure to [(18)O]water or [(13)C]glucose. The addition of glucose increased the assimilation of (18)O into DNA from [(18)O]water. However, the increase in (18)O assimilation was greater than expected based on utilization of glucose-derived carbon alone, because the addition of glucose indirectly stimulated bacteria to utilize other substrates for growth. This example illustrates the benefit of a quantitative approach to stable isotope probing.

Large-scale evaluation of experimentally determined DNA G+C contents with whole genome sequences of prokaryotes

Historically, DNA G+C content has played a critical role in the description of bacterial and archaeal species. Despite its importance in prokaryote taxonomy, its accuracy has been questioned due to methodological heterogeneity and measurement errors of conventional methods. Here we investigated the extent of accuracy of experimentally determined DNA G+C contents by comparing the reference values calculated from whole genome sequences. The large-scale comparison revealed that G+C contents determined by high-performance liquid chromatography and buoyant density centrifugation methods were more similar to the genome-derived reference values than those generated by thermal denaturation method. However, there was a substantial degree of discrepancy in DNA G+C contents between values obtained by conventional methods and genome-derived reference values. The majority of the differences between them fell out of the acceptable range (i.e. 1 mol% G+C content difference) for species delimitation of prokaryotes. In contrast, when average nucleotide identity (ANI) was correlated to G+C difference among genomes, most G+C difference was confined to less than 1% within species. Therefore, erroneous conventional methods are not meaningful in the description of bacterial and archaeal species. For taxonomic purposes, DNA G+C content should be determined by calculating directly from high-quality genome sequences with at least 16× or higher sequencing depth of coverage.

Taxonomic use of DNA G+C content and DNA-DNA hybridization in the genomic age

The G+C content of a genome is frequently used in taxonomic descriptions of species and genera. In the past it has been determined using conventional, indirect methods, but it is nowadays reasonable to calculate the DNA G+C content directly from the increasingly available and affordable genome sequences. The expected increase in accuracy, however, might alter the way in which the G+C content is used for drawing taxonomic conclusions. We here re-estimate the literature assumption that the G+C content can vary up to 3-5 % within species using genomic datasets. The resulting G+C content differences are compared with DNA-DNA hybridization (DDH) similarities calculated in silico using the GGDC web server, with 70% similarity as the gold standard threshold for species boundaries. The results indicate that the G+C content, if computed from genome sequences, varies no more than 1% within species. Statistical models based on larger differences alone can reject the hypothesis that two strains belong to the same species. Because DDH similarities between two non-type strains occur in the genomic datasets, we also examine to what extent and under which conditions such a similarity could be <70% even though the similarity of either strain to a type strain was ≥ 70%. In theory, their similarity could be as low as 50%, whereas empirical data suggest a boundary closer (but not identical) to 70%. However, it is shown that using a 50% boundary would not affect the conclusions regarding the DNA G+C content. Hence, we suggest that discrepancies between G+C content data provided in species descriptions on the one hand and those recalculated after genome sequencing on the other hand ≥ 1% are due to significant inaccuracies of the applied conventional methods and accordingly call for emendations of species descriptions.

Characterization of chloroplast DNA in Chlamydomonas eugametos and C. moewusii and its inheritance in hybrid progeny

The density, molecular weight, and cellular repetition of DNA molecules associated with the β-DNA satellite of the interfertile algae Chlamydomonas eugametos and C. moewusii are reported. The similarities between these values and those for the chloroplast DNA (cpDNA) in the related alga Chlamydomonas reinhardtii indicate that these satellites represent cpDNA. The buoyant densities of C. eugametos and C. moewusii cpDNAs are indistinguishable from one another, as are those of their respective nuclear DNAs. These densities differ slightly from the densities of the homologous components of C. reinhardtii whole cell DNA. All three species differ with respect to additional minor satellite DNAs and low molecular weight DNAs of unknown cellular location.Differences in the Aval and Smal restriction endonuclease fragmentation patterns of C. eugametos and C. moewusii cpDNAs were employed to study the inheritance of cpDNA in an F1 hybrid which had inherited a non-Mendelian streptomycin resistance marker (sr-2) from the C. eugametos mating-type plus (mt (+)) parent and in two homoplasmic mitotic segregants from a B 1 hybrid (F1 × C. moewusii) which had been initially heteroplasmic for the resistance marker. Although the cpDNA patterns in the F1 hybrid were similar to those of the C. eugametos ml (1) parent, important differences were noted which suggest that recombination between C. eugametos and C. moewusii cpDNA had occurred. Homoplasmic streptomycin resistant and sensitive mitotic segregants recovered from the B1 hybrid product reveal Aval restriction patterns similar to those of the respective resistant and sensitive parents. These data are consistent with the hypothesis that the sr-2 marker is located in cpDNA and that C. eugametos and C. moewusii cpDNA sequences can coexist in the same chloroplast and, at least sometimes, segregate without extensive recombination. The transmission of low molecular weight DNAs characteristic of C. moewusii but of unknown cellular origin shows no direct correlation with the transmission of the sr-2 marker.

Dispersive labelling of Chlamydomonas chloroplast DNA in (15)N- (14)N density transfer experiments

(15)N-(14)N density transfer experiments with synchronized vegetative cultures of Chlamydomonas reinhardtii revealed a dispersive labelling of chloroplast DNA (cpDNA) while the labelling of nuclear DNA was consistent with semiconservative replication. The dispersive labelling of cpDNA was progressive and extensive as after less than two net doublings of this DNA in (14)N-medium no significant amount of fully heavy, (15)N-strands could be detected in denatured cpDNA preparations; the average size of DNA in these preparations corresponded to 6% of the intact chloroplast genome or about 12 kbp. The density shifts of native cpDNA samples were found to be consistent with the net amounts of cpDNA synthesized. This observation indicates that essentially all (15)N atoms incorporated prior to the transfer were conserved and that metabolic turnover of cpDNA was probably absent. Our results are best explained by the exchange of homologous single-stranded segments between cpDNA molecules to form heteroduplex regions and by each DNA molecule undergoing several rounds of heteroduplex formation.

Cesium chloride-bisbenzimide gradients for separation of phytoplasma and plant DNA

CsCl-bisbezimide gradient centrifugation is a method to separate phytoplasma from host plant DNA. Bisbenzimide forms a complex with A + T-rich DNA thereby lowering its relative density. During centrifugation the A + T-rich phytoplasma DNA is spatially separated from the less A + T-rich host plant DNA. The difference in buoyant density between phytoplasma DNA and plant DNA varies according to the host-pathogen combination. The phytoplasma DNA forms a distinct band above the host plant DNA and can be collected. Depending on the phytoplasma titer and the scale of extraction, highly purified DNA is obtained in sufficient quantities for the construction of a genomic library, a sequencing project or hybridization studies.

A liquid chromatography - mass spectrometry method to measure ¹³C-isotope enrichment for DNA stable-isotope probing

DNA stable-isotope probing (DNA-SIP) is a cultivation-independent technique that makes it possible to associate metabolic function and taxonomic identity in a wide range of terrestrial and aquatic environments. In DNA-SIP, DNA is labeled via the assimilation of a labeled growth substrate that is subsequently used to identify microorganisms involved in assimilation of the substrate. However, the labeling time has to be sufficient to obtain labeled DNA but not so long such that cross-feeding of ¹³C-labeled metabolites from the primary consumers to nontarget species can occur. Confirmation that the DNA is isotopically labeled in DNA-SIP assays can be achieved using an isotope ratio mass spectrometer. In this study, we describe the development of a method using liquid chromatography (HPLC) coupled to a quadrupole mass spectrometer (QMS) to measure the ¹³C enrichment of thymine incorporated into DNA in Escherichia coli cultures fed with [¹³C]acetate. The method involved the hydrolysis of DNA extracted from the cultures that released the nucleotides, followed by the separation of the thymine by HPLC on a reverse-phase C₈ column in isocratic elution mode and the detection and quantification of ¹³C-labeled thymine by QMS. To mimic a DNA-SIP assay, a DNA mixture was made using ¹³C-labeled E. coli DNA with DNA extracted from five bacterial species. The HPLC-MS method was able to measure the correct proportion of ¹³C-DNA in the mix. This method can then be used as an alternative to the use of isotope ratio mass spectrometry in DNA-SIP assays.

Mitochondrial DNA: an advance in eukaryotic cell biology in the 1960s

Between 1950 and 1960 mitochondria were recognized as well-characterized organelles of animal and fungal cells. They shared more functional autonomy than other cellular structures. The transmission of some mitochondrial characteristics did not obey Mendelian rules and followed cytoplasmic inheritance patterns. Was this situation a consequence of still unknown complexities? We present a personal account on how approaches were set up to test very different hypotheses. In the end, it was shown that mitochondria had their own DNA, mitochondrial DNA, and that this molecule carried information specific to these organelles.

Isolation and characterization of virus-specific double-stranded DNA from tissues infected by bean golden mosaic virus

A double-stranded (ds) DNA which may be a replication intermediate was isolated from bean (Phaseolus vulgaris L. "Top Crop") leaves systemically infected with bean golden mosaic virus, a whitefly-transmitted plant virus with a genome of circular single-stranded (ss) DNA. The isolation method used phenol/chloroform extraction, hydroxyapatite column chromatography, and rate-zonal centrifugation. The dsDNA had sequences complementary to those of viral DNA. The guanine-plus-cytosine content was 35%, and the sedimentation coefficient in alkaline sucrose density gradients was similar to that of viral ssDNA. Digestion of the dsDNA by Hha I endonuclease produced fragments that corresponded exactly in number and size with those produced by complete digestion of circular viral ssDNA by Hha I, when the fragments were denatured and analyzed on polyacrylamide gels. The dsDNA molecule was a circular structure with one discontinuity in one strand; hybridization results suggest that some of a the dsDNA has a discontinuity in the viral strand and some has a discontinuity in the nonviral strand. On the basis of these structures for the dsDNA, a preliminary model for replication of viral DNA is discussed.

Characterization of the DNA of a Nonoccluded Baculovirus, Hz-1V

The DNA of the nonoccluded baculovirus (Hz-1V) obtained from the IMC-Hz-1 cell line was characterized by physicochemical and restriction endonuclease techniques. Hz-1V DNA isolated from purified virus had buoyant densities of 1.58 and 1.54 g/ml in CsCl-ethidium bromide density gradients, which corresponded to supercoiled and to relaxed circular and linear DNA, respectively. Neutral CsCl equilibrium centrifugation indicated that the Hz-1V DNA had a buoyant density of 1.7024 g/ml, which corresponded to a guanine-plus-cytosine (G+C) content of 43%. Thermal denaturation indicated a high G+C domain(s) in the Hz-1V genomic DNA. The domain(s), which included about 11% of the total genomic DNA, exhibited a T(m) of 97 degrees C. The remaining portion (89%) of the DNA had a T(m) of 86.5 degrees C. The T(m)s corresponded to G+C contents of 42 and 67%, respectively. The mean genetic complexity of Hz-1V DNA determined by DNA reassociation kinetic analysis was found to be 152 x 10(6). A possible rapidly reassociating component comprising approximately 13% of the genome was observed. The mean molecular weights from restriction endonuclease digests were 159 x 10(6) for both HindIII and EcoRI. Genomic heterogeneity was found in both the wild-type Hz-1V stock and in two plaque isolates. Of 12 single-plaque isolates, 3 basic restriction endonuclease DNA fragment patterns were observed. The molecular size estimates from electron microscopic contour lengths of uncloned viral DNA ranged from 70 to 158 megadaltons, and the mode was the 130- to 140-megadalton class.