Correcting PCR amplification errors in unique molecular identifiers to generate accurate numbers of sequencing molecules

Unique molecular identifiers are random oligonucleotide sequences that remove PCR amplification biases. However, the impact that PCR associated sequencing errors have on the accuracy of generating absolute counts of RNA molecules is underappreciated. We show that PCR errors are a source of inaccuracy in both bulk and single-cell sequencing data, and synthesizing unique molecular identifiers using homotrimeric nucleotide blocks provides an error-correcting solution that allows absolute counting of sequenced molecules.

A rapid multiplex loop-mediated isothermal amplification (mLAMP) assay for detection of Entamoeba histolytica and Giardia duodenalis

We developed a rapid multiplex loop-mediated isothermal amplification (mLAMP) assay for two common intestinal parasites-Entamoeba histolytica and Giardia duodenalis, where early detection may be helpful. The mLAMP assay was optimized for the detection of DNA of E. histolytica (18S rRNA gene) and G. duodenalis (Elongation factor 1 alpha gene) from standard strains by using six specific primers FIP (forward inner primer), BIP (backward inner primer), F3 (forward outer primer), B3 (backward outer primer), loopF (forward loop primer), and loopB (backward loop primer) for each gene target. The amplification time was 16-26 min for E. histolytica and 10-15 min for G. duodenalis, and the parasites could be distinguished based on melting-curve analysis for specific annealing temperatures (Tm) of 84°C-86°C and 88°C-90°C for E. histolytica and G. duodenalis, respectively. The analytical sensitivity was one fg, and no cross-reactivity with other intestinal pathogens was observed. Thus, the mLAMP assay could detect and clearly distinguish E. histolytica and G. duodenalis with a rapid turnaround time and excellent analytical sensitivity and specificity.

Whole Genome Amplification in Preimplantation Genetic Testing in the Era of Massively Parallel Sequencing

Successful whole genome amplification (WGA) is a cornerstone of contemporary preimplantation genetic testing (PGT). Choosing the most suitable WGA technique for PGT can be particularly challenging because each WGA technique performs differently in combination with different downstream processing and detection methods. The aim of this review is to provide insight into the performance and drawbacks of DOP-PCR, MDA and MALBAC, as well as the hybrid WGA techniques most widely used in PGT. As the field of PGT is moving towards a wide adaptation of comprehensive massively parallel sequencing (MPS)-based approaches, we especially focus our review on MPS parameters and detection opportunities of WGA-amplified material, i.e., mappability of reads, uniformity of coverage and its influence on copy number variation analysis, and genomic coverage and its influence on single nucleotide variation calling. The ability of MDA-based WGA solutions to better cover the targeted genome and the ability of PCR-based solutions to provide better uniformity of coverage are highlighted. While numerous comprehensive PGT solutions exploiting different WGA types and adjusted bioinformatic pipelines to detect copy number and single nucleotide changes are available, the ones exploiting MDA appear more advantageous. The opportunity to fully analyse the targeted genome is influenced by the MPS parameters themselves rather than the solely chosen WGA.

Paleo-diatom composition from Santa Barbara Basin deep-sea sediments: a comparison of 18S-V9 and diat-rbcL metabarcoding vs shotgun metagenomics

Sedimentary ancient DNA (sedaDNA) analyses are increasingly used to reconstruct marine ecosystems. The majority of marine sedaDNA studies use a metabarcoding approach (extraction and analysis of specific DNA fragments of a defined length), targeting short taxonomic marker genes. Promising examples are 18S-V9 rRNA (~121-130 base pairs, bp) and diat-rbcL (76 bp), targeting eukaryotes and diatoms, respectively. However, it remains unknown how 18S-V9 and diat-rbcL derived compositional profiles compare to metagenomic shotgun data, the preferred method for ancient DNA analyses as amplification biases are minimised. We extracted DNA from five Santa Barbara Basin sediment samples (up to ~11 000 years old) and applied both a metabarcoding (18S-V9 rRNA, diat-rbcL) and a metagenomic shotgun approach to (i) compare eukaryote, especially diatom, composition, and (ii) assess sequence length and database related biases. Eukaryote composition differed considerably between shotgun and metabarcoding data, which was related to differences in read lengths (~112 and ~161 bp, respectively), and overamplification of short reads in metabarcoding data. Diatom composition was influenced by reference bias that was exacerbated in metabarcoding data and characterised by increased representation of Chaetoceros, Thalassiosira and Pseudo-nitzschia. Our results are relevant to sedaDNA studies aiming to accurately characterise paleo-ecosystems from either metabarcoding or metagenomic data.

Directionally biased habitat shifts and biogeographically informative cytonuclear discordance in the Hawaiian silversword alliance (Compositae)

PREMISE

Expanded phylogenetic analyses of the Hawaiian silversword alliance (Argyroxiphium, Dubautia, Wilkesia; Compositae) were undertaken to assess evolutionary and biogeographic informativeness of cytonuclear discordance and any biases in evolutionary directionality of ecological transitions within this prominent example of adaptive radiation.

METHODS

Samples spanning the geographic and ecological distributions of all recognized taxa were included in phylogenetic and biogeographic analyses of nuclear ribosomal DNA (nrDNA) and cpDNA sequences. Bayesian model testing approaches were used to model ecological evolution and the evolution of nuclear chromosomal arrangements while accounting for phylogenetic uncertainty.

RESULTS

Cytonuclear discordance detected previously appears to reflect chloroplast capture, at least in part, with nrDNA trees being largely congruent with nuclear chromosomal structural data and fine-scale taxonomy. Comparison of biogeographic histories estimated from the posterior distributions of nrDNA and cpDNA trees, including inferred chloroplast-capture events, provides additional resolution of dispersal history, including a back-dispersal to Maui Nui from Hawai'i. A newly resolved major nrDNA clade of endemic Kaua'i taxa that mostly were described as new-to-science since the 1980s strengthens the earlier hypothesis that diversification on Kaua'i has not waned since the island began to decline in area through subsidence and erosion. Bias in habitat shifts was estimated, with transitions from dry-to-mesic or -wet and from wet-to-mesic or -bog habitats dominating diversification of the silversword alliance from a dry-adapted tarweed ancestor.

CONCLUSIONS

The habitat-transition biases estimated here may indicate more limited pathways of ecological evolution than proposed previously for an adaptive radiation involving such major ecological shifts.

An optimized method for the extraction of ancient eukaryote DNA from marine sediments

Marine sedimentary ancient DNA (sedaDNA) provides a powerful means to reconstruct marine palaeo-communities across the food web. However, currently there are few optimized sedaDNA extraction protocols available to maximize the yield of small DNA fragments typical of ancient DNA (aDNA) across a broad diversity of eukaryotes. We compared seven combinations of sedaDNA extraction treatments and sequencing library preparations using marine sediments collected at a water depth of 104 m off Maria Island, Tasmania, in 2018. These seven methods contrasted frozen versus refrigerated sediment, bead-beating induced cell lysis versus ethylenediaminetetraacetic acid (EDTA) incubation, DNA binding in silica spin columns versus in silica-solution, diluted versus undiluted DNA in shotgun library preparations to test potential inhibition issues during amplification steps, and size-selection of low molecular-weight (LMW) DNA to increase the extraction efficiency of sedaDNA. Maximum efficiency was obtained from frozen sediments subjected to a combination of EDTA incubation and bead-beating, DNA binding in silica-solution, and undiluted DNA in shotgun libraries, across 45 marine eukaryotic taxa. We present an optimized extraction protocol integrating these steps, with an optional post-library LMW size-selection step to retain DNA fragments of ≤500 base pairs. We also describe a stringent bioinformatic filtering approach for metagenomic data and provide a comprehensive list of contaminants as a reference for future sedaDNA studies. The new extraction and data-processing protocol should improve quantitative paleo-monitoring of eukaryotes from marine sediments, as well as other studies relying on the detection of highly fragmented and degraded eukaryote DNA in sediments.

Molecular evolution of the internal transcribed spacers in red oaks (Quercus sect. Lobatae)

Previous studies of the Internal Transcribed Spacers of the nuclear ribosomal DNA (ITS) in sections Quercus (white oaks), Protobalanus (intermediate or golden cup oaks), Cerris (Cerris oaks), and Ilex (Ilex oaks) suggest that ITS regions undergo full concerted evolution in oaks; however, ITS evolution patterns in red oaks (section Lobatae) are unknown due to scant representation in published work. To determine whether full concerted evolution occurs in red oaks, the purpose of this study was to examine ITS sequences from 40 red oak species. The results show incomplete concerted evolution and the presence of three ITS ribotypes of lengths 505, 609, 601 bp, hereafter referred to as ITS-S (small), I ITS-M (medium), and ITS-L (large), respectively. Thirty species had only one ribotype (ITS-M), nine species had two ribotypes (different combinations of ITS-L, ITS-M, and ITS-S), and only one species had all three ribotypes. Furthermore, examination of these three ribotypes showed that only ITS-M is putatively functional and ITS-L and ITS-S are pseudogenes. Bayesian analysis strongly supported (100%) two pseudogenes clades but provided weak support for the monophyly of a putative functional clade (ITS-M); moreover, within the "functional" clade, species relationships were uncertain and, in most cases, sequences from the same species failed to group together. The results of the current study suggest that ITS may not be appropriate for phylogeny reconstruction of red oaks due to low levels of interspecific variation and incomplete concerted evolution.

Consistent and correctable bias in metagenomic sequencing experiments

Marker-gene and metagenomic sequencing have profoundly expanded our ability to measure biological communities. But the measurements they provide differ from the truth, often dramatically, because these experiments are biased toward detecting some taxa over others. This experimental bias makes the taxon or gene abundances measured by different protocols quantitatively incomparable and can lead to spurious biological conclusions. We propose a mathematical model for how bias distorts community measurements based on the properties of real experiments. We validate this model with 16S rRNA gene and shotgun metagenomics data from defined bacterial communities. Our model better fits the experimental data despite being simpler than previous models. We illustrate how our model can be used to evaluate protocols, to understand the effect of bias on downstream statistical analyses, and to measure and correct bias given suitable calibration controls. These results illuminate new avenues toward truly quantitative and reproducible metagenomics measurements.

Consistent and correctable bias in metagenomic sequencing experiments

Marker-gene and metagenomic sequencing have profoundly expanded our ability to measure biological communities. But the measurements they provide differ from the truth, often dramatically, because these experiments are biased toward detecting some taxa over others. This experimental bias makes the taxon or gene abundances measured by different protocols quantitatively incomparable and can lead to spurious biological conclusions. We propose a mathematical model for how bias distorts community measurements based on the properties of real experiments. We validate this model with 16S rRNA gene and shotgun metagenomics data from defined bacterial communities. Our model better fits the experimental data despite being simpler than previous models. We illustrate how our model can be used to evaluate protocols, to understand the effect of bias on downstream statistical analyses, and to measure and correct bias given suitable calibration controls. These results illuminate new avenues toward truly quantitative and reproducible metagenomics measurements.

Evaluation of multiplex PCR assay for detection of Babesia spp, Ehrlichia canis and Trypanosoma evansi in dogs

A multiplex PCR test was evaluated to detect the DNA of three important dog haemoparasites by comparing with singular PCR counterpart on clinical blood samples of dogs in and around Chennai, Tamil Nadu, India. Initial screening of samples was done by microscopic examination of peripheral blood smear and singular PCR and those found exclusively positive for Babesia spp, Ehrlichia canis and Trypanosoma evansi and concurrent infections were used to standardize multiplex PCR. Amplicons of 619 bp, 377 bp and 227 bp corresponding to Babesia spp (18S rRNA gene), E. canis (VirB9 gene), and T.evansi (VSG gene) respectively were amplified, without any non-specific amplification. The laboratory sensitivity (91.7% to 100%) and specificity (100%) of the multiplex PCR were calculated using 'true positive' and 'true negative' dog blood samples obtained in the initial screening process. Clinical blood samples from 287 dogs were screened using singular PCR and multiplex PCR tests for the presence of genome of Babesia spp, E. canis and T. evansi. The multiplex PCR was found to have high level of diagnostic specificity (97.5%-100%) in the detection of all three dog blood parasites and high level of diagnostic sensitivity (95%) in the detection of T. evansi from field level clinical blood samples compared to the singular PCR. However, the diagnostic sensitivity of the multiplex PCR was found to be low to moderate (40.45%-66.7%) in detection of Babesia spp and E. canis from field level clinical blood samples. The strength of agreement between singular and multiplex PCR assays was 'moderate' (0.445), 'good' (0.708) and 'very good' (0.968) in detection of DNA of Babesia spp, E. canis and T. evansi. The multiplex PCR was found to be 10 fold less sensitive in comparison with the singular PCR counterpart.

Debugging diversity - a pan-continental exploration of the potential of terrestrial blood-feeding leeches as a vertebrate monitoring tool

The use of environmental DNA (eDNA) has become an applicable noninvasive tool with which to obtain information about biodiversity. A subdiscipline of eDNA is iDNA (invertebrate-derived DNA), where genetic material ingested by invertebrates is used to characterize the biodiversity of the species that served as hosts. While promising, these techniques are still in their infancy, as they have only been explored on limited numbers of samples from only a single or a few different locations. In this study, we investigate the suitability of iDNA extracted from more than 3,000 haematophagous terrestrial leeches as a tool for detecting a wide range of terrestrial vertebrates across five different geographical regions on three different continents. These regions cover almost the full geographical range of haematophagous terrestrial leeches, thus representing all parts of the world where this method might apply. We identify host taxa through metabarcoding coupled with high-throughput sequencing on Illumina and IonTorrent sequencing platforms to decrease economic costs and workload and thereby make the approach attractive for practitioners in conservation management. We identified hosts in four different taxonomic vertebrate classes: mammals, birds, reptiles and amphibians, belonging to at least 42 different taxonomic families. We find that vertebrate blood ingested by haematophagous terrestrial leeches throughout their distribution is a viable source of DNA with which to examine a wide range of vertebrates. Thus, this study provides encouraging support for the potential of haematophagous terrestrial leeches as a tool for detecting and monitoring terrestrial vertebrate biodiversity.

Genomic GC-Content Affects the Accuracy of 16S rRNA Gene Sequencing Based Microbial Profiling due to PCR Bias

Profiling of microbial community composition is frequently performed by partial 16S rRNA gene sequencing on benchtop platforms following PCR amplification of specific hypervariable regions within this gene. Accuracy and reproducibility of this strategy are two key parameters to consider, which may be influenced during all processes from sample collection and storage, through DNA extraction and PCR based library preparation to the final sequencing. In order to evaluate both the reproducibility and accuracy of 16S rRNA gene based microbial profiling using the Ion Torrent PGM platform, we prepared libraries and performed sequencing of a well-defined and validated 20-member bacterial DNA mock community on five separate occasions and compared results with the expected even distribution. In general the applied method had a median coefficient of variance of 11.8% (range 5.5–73.7%) for all 20 included strains in the mock community across five separate sequencing runs, with underrepresented strains generally showing the largest degree of variation. In terms of accuracy, mock community species belonging to Proteobacteria were underestimated, whereas those belonging to Firmicutes were mostly overestimated. This could be explained partly by premature read truncation, but to larger degree their genomic GC-content, which correlated negatively with the observed relative abundances, suggesting a PCR bias against GC-rich species during library preparation. Increasing the initial denaturation time during the PCR amplification from 30 to 120 s resulted in an increased average relative abundance of the three mock community members with the highest genomic GC%, but did not significantly change the overall evenness of the community distribution. Therefore, efforts should be made to optimize the PCR conditions prior to sequencing in order to maximize accuracy.

The Use of Degenerate Primers in qPCR Analysis of Functional Genes Can Cause Dramatic Quantification Bias as Revealed by Investigation of nifH Primer Performance

The measurement of functional gene abundance in diverse microbial communities often employs quantitative PCR (qPCR) with highly degenerate oligonucleotide primers. While degenerate PCR primers have been demonstrated to cause template-specific bias in PCR applications, the effect of such bias on qPCR has been less well explored. We used a set of diverse, full-length nifH gene standards to test the performance of several universal nifH primer sets in qPCR. We found significant template-specific bias in all but the PolF/PolR primer set. Template-specific bias caused more than 1000-fold mis-estimation of nifH gene copy number for three of the primer sets and one primer set resulted in more than 10,000-fold mis-estimation. Furthermore, such template-specific bias will cause qPCR estimates to vary in response to beta-diversity, thereby causing mis-estimation of changes in gene copy number. A reduction in bias was achieved by increasing the primer concentration. We conclude that degenerate primers should be evaluated across a range of templates, annealing temperatures, and primer concentrations to evaluate the potential for template-specific bias prior to their use in qPCR.

Transcriptional activity of PIF and Pong-like Class II transposable elements in Triticeae

Background

Transposable elements are major contributors to genome size and variability, accounting for approximately 70–80% of the maize, barley, and wheat genomes. PIF and Pong-like elements belong to two closely-related element families within the PIF/Harbinger superfamily of Class II (DNA) transposons. Both elements contain two open reading frames; one encodes a transposase (ORF2) that catalyzes transposition of the functional elements and their related non-autonomous elements, while the function of the second is still debated. In this work, we surveyed for PIF- and Pong-related transcriptional activity in 13 diploid Triticeae species, all of which have been previously shown to harbor extensive within-genome diversity of both groups of elements.

Results

The results revealed that PIF elements have considerable transcriptional activity in Triticeae, suggesting that they can escape the initial levels of plant cell control and are regulated at the post-transcriptional level. Phylogenetic analysis of 156 PIF cDNA transposase fragments along with 240 genomic partial transposase sequences showed that most, if not all, PIF clades are transcriptionally competent, and that multiple transposases coexisting within a single genome have the potential to act simultaneously. In contrast, we did not detect any transcriptional activity of Pong elements in any sample.

Conclusions

The lack of Pong element transcription shows that even closely related transposon families can exhibit wide variation in their transposase transcriptional activity within the same genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1028-6) contains supplementary material, which is available to authorized users.

Sensitivity and accuracy of high-throughput metabarcoding methods for early detection of invasive fish species

High-throughput DNA metabarcoding has gained recognition as a potentially powerful tool for biomonitoring, including early detection of aquatic invasive species (AIS). DNA based techniques are advancing, but our understanding of the limits to detection for metabarcoding complex samples is inadequate. For detecting AIS at an early stage of invasion when the species is rare, accuracy at low detection limits is key. To evaluate the utility of metabarcoding in future fish community monitoring programs, we conducted several experiments to determine the sensitivity and accuracy of routine metabarcoding methods. Experimental mixes used larval fish tissue from multiple “common” species spiked with varying proportions of tissue from an additional “rare” species. Pyrosequencing of genetic marker, COI (cytochrome c oxidase subunit I) and subsequent sequence data analysis provided experimental evidence of low-level detection of the target “rare” species at biomass percentages as low as 0.02% of total sample biomass. Limits to detection varied interspecifically and were susceptible to amplification bias. Moreover, results showed some data processing methods can skew sequence-based biodiversity measurements from corresponding relative biomass abundances and increase false absences. We suggest caution in interpreting presence/absence and relative abundance in larval fish assemblages until metabarcoding methods are optimized for accuracy and precision.

Development of a multiplex loop-mediated isothermal amplification method for the simultaneous detection of Salmonella spp. and Vibrio parahaemolyticus

Rapid detection of food-borne pathogens is important in the food industry, to monitor and prevent the spread of these pathogens through contaminated food products. We therefore established a multiplex real-time loop-mediated isothermal amplification (LAMP) assay to simultaneously detect and distinguish Salmonella spp. and Vibrio parahaemolyticus DNA in a single reaction. Two target sequences, one specific for Salmonella and the other specific for Vibrio parahaemolyticus, were amplified by specific LAMP primers in the same reaction tube. After amplification at 65 °C for 60 min, the amplified products were subjected to melting curve analysis and thus could be distinguished based on the different melting temperatures (T_m values) of the two specifically amplified products. The specificity of the multiplex LAMP assay was evaluated using 19 known bacterial strains, including one V. parahaemolyticus and seven Salmonella spp. strains. The multiplex LAMP showed 100% inclusivity and exclusivity, and a detection limit similar to that of multiplex PCR. In addition, we observed and corrected preferential amplification induced by what we call LAMP selection in the multiplex LAMP reaction. In conclusion, our assay was rapid, specific, and quantitative, making it a useful tool for the food industry.

Mutations in two pollen self-incompatibility factors in geographically marginal populations of Solanum habrochaites impact mating system transitions and reproductive isolation

PREMISE OF THE STUDY

Self-incompatibility (SI) is a mechanism that prevents inbreeding in many plant species. The mutational breakdown of SI occurs frequently, yet relatively little is known about the evolutionary steps involved in the progressive loss of pistil and pollen SI function.

METHODS

In Solanaceae, SI is the S-RNase-based gametophytic type. We used SI and SC populations of the wild tomato species Solanum habrochaites to study natural variation for two pollen SI factors: a Cullin1 (CUL1) protein and an S-locus F-box protein (SLF-23). Pollen compatibility was assessed on an allotriploid tester line encoding an S-RNase recognized by SLF-23. Both pollen factors are required for compatibility on this tester line. Complementation tests and gene sequencing were used to identify mutations in CUL1 or SLF-23.

KEY RESULTS

We detected loss-of-function mutations in CUL1 and/or SLF-23 in SC populations collected near the northern and southern geographic margins of this taxon's natural range. Nonmarginal SC and all SI accessions expressed mostly functional alleles of these pollen factors. Comparison of the CUL1 sequences identified several shared deletion mutations present in both northern and southern margin SC accessions.

CONCLUSIONS

Loss-of-function mutations in CUL1 and SLF-23 likely became fixed relatively late during SI to SC transitions, after loss of pistil SI function. Mutations in CUL1 establish unilateral incompatibility with SI populations and strengthen reproductive isolation. Point mutations common to northern and southern SC biotypes likely derive from shared ancestral variants found in more central SI populations.

Across-shelf distribution of blue mussel larvae in the northern Gulf of Maine: consequences for population connectivity and a species range boundary

Studies of population connectivity have largely focused on along-shelf, as opposed to across-shelf, processes. We hypothesized that a discontinuity in across-shelf mixing caused by the divergence of the Eastern Maine Coastal Current (EMCC) from shore acts as an ecological barrier to the supply of mussel larvae to the coast. Existing data on the relative abundance of two congeneric blue mussels, Mytilus edulis and M. trossulus, were analysed to quantify the association of M. trossulus with the colder temperature signal of the EMCC and generate larval distribution predictions. We then sampled the across-shelf distribution of larvae along two transects during 2011. Larvae were identified using restriction digests of PCR amplicons from the mitochondrial 16S rDNA. Mytilus edulis larvae were consistently abundant on either the inshore and offshore transect ends, but not homogeneously distributed across the shelf, while M. trossulus larvae were less common throughout the study area. The divergence of the EMCC from shore appears to create a break in the connectivity of M. edulis populations by isolating those inshore of the EMCC from upstream larval sources. Across-shelf transport processes can thus produce connectivity patterns that would not be predicted solely on the basis of along-shelf processes.

Diversity, abundance, and evolutionary dynamics of Pong-like transposable elements in Triticeae

Pong-like elements are members of the PIF/Harbinger superfamily of DNA transposons that has been described in many plants, animals, and fungi. Most Pong elements contain two open reading frames (ORFs). One encodes a transposase (ORF2) that catalyzes transposition of Pong and related non-autonomous elements, while the function of the second is unknown. Little is known about the evolutionary history of Pong elements in flowering plants. In this work, we present the first comprehensive analysis of the diversity, abundance, and evolution of the Pong-like transposase gene in the genomes of 21 diploid species from the wheat tribe, Triticeae, and we present the first convincing evidence of horizontal transfer of nuclear-encoded Pong elements in any organism. A phylogenetic analysis of nearly 300 Pong sequences based on a conserved region of the transposase domain revealed a complex evolutionary history of Pong elements that can be best explained by ancestral polymorphism, followed by differential evolutionary success of some transposase lineages, and by occasional horizontal transfer between phylogenetically distant genera. In addition, we used transposon display to estimate the abundance of the transposase gene within Triticeae genomes, and our results revealed varying levels of Pong proliferation, with numbers of transposase copies ranging from 22 to 92. Comparisons of Pong transposase abundance to flow cytometry estimates of genome size revealed that larger Triticeae genome size was not correlated with transposase abundance.

Enzymological description of multitemplate PCR-Shrinking amplification bias by optimizing the polymerase-template ratio

Multitemplate polymerase chain reaction (PCR) is used for preparative and analytical applications in diagnostics and research. Classical PCR and qPCR are two basic setups with many possible experimental modifications. Classical PCR is a method of choice to obtain enough material for subsequent sophisticated applications such as construction of libraries for next-generation sequencing or high-throughput screening. Sequencing and Single Nucleotide Primer Extension (SNuPE) employ one-strand synthesis and represent a distinct variant of analytical DNA synthesis. In all these applications, maintaining the initial ratio of templates and avoiding underestimation of minority templates is desired. Here, we demonstrate that different templates can amplify independently at low template concentrations (typical in qPCR setups, in which the polymerase concentration is usually several orders of magnitude higher than the template concentration). However, rare templates can be diluted in an effort to keep DNA amplification in the exponential phase, or template concentration can be biased by differences in amplification efficiency. Moreover, amplification of templates present in low concentrations is more vulnerable to stochastic events that lead to proportional changes in the product ratio, as well as by incomplete amplification leading to chimera formation. These undesired effects can be compensated for by using highly processive polymerases with high and equal affinity to different primer-template complexes. Novel enhanced polymerases are desired. With increasing concentration of a primer-template of interest, the system becomes more deterministic. Nevertheless, marked deviation from independent exponential amplification occurs when the total template concentration starts to approach the polymerase concentration. The primer-template complexes compete for enzyme molecules, and the amount of products grows arithmetically-the system starts to obey Michaelis-Menten kinetics. Synthesis of rare products in a multitemplate mixture can run more easily under the detection limit in such conditions, although it would be unequivocally detectable in a single template assay. When fishing out rare template variants, the best processive polymerases should be used to decrease both amplification and detection limits. The possibility of stochastic events, should be taken into account to correctly interpret the obtained data.

Effect of the enzyme and PCR conditions on the quality of high-throughput DNA sequencing results

Library preparation protocols for high-throughput DNA sequencing (HTS) include amplification steps in which errors can build up. In order to have confidence in the sequencing data, it is important to understand the effects of different Taq polymerases and PCR amplification protocols on the DNA molecules sequenced. We compared thirteen enzymes in three different marker systems: simple, single copy nuclear gene and complex multi-gene family. We also tested a modified PCR protocol, which has been suggested to reduce errors associated with amplification steps. We find that enzyme choice has a large impact on the proportion of correct sequences recovered. The most complex marker systems yielded fewer correct reads, and the proportion of correct reads was greatly affected by the enzyme used. Modified cycling conditions did reduce the number of incorrect sequences obtained in some cases, but enzyme had a much greater impact on the number of correct reads. Thus, the coverage required for the safe identification of genotypes using one of the low quality enzymes could be seven times larger than with more efficient enzymes in a biallelic system with equal amplification of the two alleles. Consequently, enzyme selection for downstream HTS has important consequences, especially in complex genetic systems.

Improved pipeline for reducing erroneous identification by 16S rRNA sequences using the Illumina MiSeq platform

The cost of DNA sequencing has decreased due to advancements in Next Generation Sequencing. The number of sequences obtained from the Illumina platform is large, use of this platform can reduce costs more than the 454 pyrosequencer. However, the Illumina platform has other challenges, including bioinformatics analysis of large numbers of sequences and the need to reduce erroneous nucleotides generated at the 3'-ends of the sequences. These erroneous sequences can lead to errors in analysis of microbial communities. Therefore, correction of these erroneous sequences is necessary for accurate taxonomic identification. Several studies that have used the Illumina platform to perform metagenomic analyses proposed curating pipelines to increase accuracy. In this study, we evaluated the likelihood of obtaining an erroneous microbial composition using the MiSeq 250 bp paired sequence platform and improved the pipeline to reduce erroneous identifications. We compared different sequencing conditions by varying the percentage of control phiX added, the concentration of the sequencing library, and the 16S rRNA gene target region using a mock community sample composed of known sequences. Our recommended method corrected erroneous nucleotides and improved identification accuracy. Overall, 99.5% of the total reads shared 95% similarity with the corresponding template sequences and 93.6% of the total reads shared over 97% similarity. This indicated that the MiSeq platform can be used to analyze microbial communities at the genus level with high accuracy. The improved analysis method recommended in this study can be applied to amplicon studies in various environments using high-throughput reads generated on the MiSeq platform.

The Single-Cell Lab or How to Perform Single-Cell Molecular Analysis

The increasing interest towards cellular heterogeneity within cell populations has pushed the development of new protocols to isolate and analyze single cells. PCR-based amplification techniques are widely used in this field. However, setting up an experiment and analyzing the results can sometimes be challenging. The aim of this chapter is to provide a general overview on single-cell PCR analysis focusing on the potential pitfalls and on the possible solutions to successfully perform the analysis.

Bias in Whole Genome Amplification: Causes and Considerations

Whole genome amplification (WGA) is a widely used molecular technique that is becoming increasingly necessary in genetic research on a range of sample types including individual cells, fossilized remains and entire ecosystems. Multiple methods of WGA have been developed, each with specific strengths and weaknesses, but with a common defect in that each method distorts the initial template DNA during the course of amplification. The type, extent, and circumstance of the bias vary with the WGA method and particulars of the template DNA. In this review, we endeavor to discuss the types of bias introduced, the susceptibility of common WGA techniques to these bias types, and the interdependence between bias and characteristics of the template DNA. Finally, we attempt to illustrate some of the criteria specific to the analytical platform and research application that should be considered to enable combination of the appropriate WGA method, template DNA, sequencing platform, and intended use for optimal results.

Multi-template polymerase chain reaction

PCR is a formidable and potent technology that serves as an indispensable tool in a wide range of biological disciplines. However, due to the ease of use and often lack of rigorous standards many PCR applications can lead to highly variable, inaccurate, and ultimately meaningless results. Thus, rigorous method validation must precede its broad adoption to any new application. Multi-template samples possess particular features, which make their PCR analysis prone to artifacts and biases: multiple homologous templates present in copy numbers that vary within several orders of magnitude. Such conditions are a breeding ground for chimeras and heteroduplexes. Differences in template amplification efficiencies and template competition for reaction compounds undermine correct preservation of the original template ratio. In addition, the presence of inhibitors aggravates all of the above-mentioned problems. Inhibitors might also have ambivalent effects on the different templates within the same sample. Yet, no standard approaches exist for monitoring inhibitory effects in multitemplate PCR, which is crucial for establishing compatibility between samples.

Incomplete sequence homogenization in 45S rDNA multigene families: intermixed IGS heterogeneity within the single NOR locus of the polyploid species Medicago arborea (Fabaceae)

BACKGROUND AND AIMS

Ribosomal sequences have become the classical example of the genomic homogenization of nuclear multigene families. Despite theoretical advantages and modelling predictions that support concerted evolution of the 45S rDNA, several reports have found intragenomic polymorphisms. However, the origins and causes of these rDNA polymorphisms are difficult to assess because seed plants show a wide range of 45S rDNA loci number variation, especially in polyploids. Medicago arborea is a tetraploid species that has a single 45S rDNA locus. This feature makes this species a suitable case study to assess the fate of ribosomal IGS homogenization in polyploid species showing nucleolus organizer region (NOR) reduction.

METHODS

The intergenic spacer (IGS) region was amplified by long PCR and the fragments were cloned and sequenced by a primer-walking strategy. The physical mapping of the whole and partial IGS variants was assessed by fluorescent in situ hybridization (FISH) and fibre-FISH methods on mitotic chromosomes and extended DNA fibres, respectively.

KEY RESULTS

Two IGS fragments of 4·8 and 3·5 kb were obtained showing structural features of functional sequences. The shorter variant appears to be a truncated copy of the 4·8 kb fragment that lacks the duplication of the transcription initiation site region and the entire D region. The physical localization of the two IGS variants on metaphase chromosomes and extended DNA fibres using FISH corroborated their joint presence within the same locus. In addition, no spatial structure of the two variants was detected within the NOR.

CONCLUSIONS

The results suggest that full sequence homogenization is not operating within the NOR locus of M. arborea. The structure of the NOR locus reported here departs from the models of IGS heterogeneity present in plants and caution against assuming the widespread belief that intragenomic ribosomal heterogeneity is mainly due to sequence variation between paralogous loci.

Evaluating bias of illumina-based bacterial 16S rRNA gene profiles

Massively parallel sequencing of 16S rRNA genes enables the comparison of terrestrial, aquatic, and host-associated microbial communities with sufficient sequencing depth for robust assessments of both alpha and beta diversity. Establishing standardized protocols for the analysis of microbial communities is dependent on increasing the reproducibility of PCR-based molecular surveys by minimizing sources of methodological bias. In this study, we tested the effects of template concentration, pooling of PCR amplicons, and sample preparation/interlane sequencing on the reproducibility associated with paired-end Illumina sequencing of bacterial 16S rRNA genes. Using DNA extracts from soil and fecal samples as templates, we sequenced pooled amplicons and individual reactions for both high (5- to 10-ng) and low (0.1-ng) template concentrations. In addition, all experimental manipulations were repeated on two separate days and sequenced on two different Illumina MiSeq lanes. Although within-sample sequence profiles were highly consistent, template concentration had a significant impact on sample profile variability for most samples. Pooling of multiple PCR amplicons, sample preparation, and interlane variability did not influence sample sequence data significantly. This systematic analysis underlines the importance of optimizing template concentration in order to minimize variability in microbial-community surveys and indicates that the practice of pooling multiple PCR amplicons prior to sequencing contributes proportionally less to reducing bias in 16S rRNA gene surveys with next-generation sequencing.

Internal validation of human mitochondrial DNA quantification using real-time PCR

The quantity of mitochondrial DNA (mtDNA) template added for amplification and subsequent dye terminator reactions is critical for obtaining quality sequence data. Validation of a human mtDNA real-time quantitative PCR (qPCR) assay demonstrated its high degree of reproducibility and precision as well as an extremely sensitive threshold of detection (0.0001 pg/μL or approximately six human mtDNA copies/μL). A study of 35 nonprobative bone and teeth evidence samples revealed that 20 pg of mtDNA template is recommended for successful HV1 and HV2 sequence analysis; however, as little as 0.013 pg can generate a full mtDNA profile when using enhanced amplification reactions. The assay can also detect PCR inhibition and is useful for identifying samples that may benefit from re-purification. Overall, the assay is an excellent method to quantify mtDNA and is useful for determining the best analytical approach for successful sequencing.

PCR-based isolation of multigene families: lessons from the avian MHC class IIB

The amount of sequence data available today highly facilitates the access to genes from many gene families. Primers amplifying the desired genes over a range of species are readily obtained by aligning conserved gene regions, and laborious gene isolation procedures can often be replaced by quicker PCR-based approaches. However, in the case of multigene families, PCR-based approaches bear the often ignored risk of incomplete isolation of family members. This problem is most prominent in gene families with highly variable and thus unpredictable number of gene copies among species, such as in the major histocompatibility complex (MHC). In this study, we (i) report new primers for the isolation of the MHC class IIB (MHCIIB) gene family in birds and (ii) share our experience with isolating MHCIIB genes from an unprecedented number of avian species from all over the avian phylogeny. We report important and usually underappreciated problems encountered during PCR-based multigene family isolation and provide a collection of measures to help significantly improving the chance of successfully isolating complete multigene families using PCR-based approaches.

Capturing the population structure of microparasites: using ITS-sequence data and a pooled DNA approach

The internal transcribed spacer (ITS) region of nuclear ribosomal DNA is a common marker not only for the molecular identification of different taxa and strains, but also for the analysis of population structure of wild microparasite communities. Importantly, the multicopy nature of this region allows the amplification of low-quantity samples of the target DNA, a common problem in studies on unicellular, unculturable microparasites. We analysed ITS sequences from the protozoan parasite Caullerya mesnili (class Ichthyosporea) infecting waterflea (Daphnia) hosts, across several host population samples. We showed that analysing representative ITS-types [as identified by statistical parsimony network analysis (SPN)] is a suitable method to address relevant polymorphism. The spatial patterns were consistent regardless of whether parasite DNA was extracted from individual hosts or pooled host samples. Remarkably, the efficiency in detecting different sequence types was even higher after sample pooling. As shown by simulations, an easily manageable number of sequences from pooled DNA samples are sufficient to resolve the spatial population structure in this system. In summary, the ITS region analysed from pooled DNA samples can provide valuable insights into the spatial and temporal dynamics of microparasites. Moreover, the application of SPN analysis is a good alternative to the well-established neighbour-joining method (NJ) for the identification of representative ITS-types. SPN can even outperform NJ by joining most of the singleton sequences to representative sequence clusters.

External and semi-internal controls for PCR amplification of homologous sequences in mixed templates

In a mixed template, the presence of homologous target DNA sequences creates environments that almost inevitably give rise to artifacts and biases during PCR. Heteroduplexes, chimeras, and skewed template-to-product ratios are the exclusive attributes of mixed template PCR and never occur in a single template assay. Yet, multi-template PCR has been used without appropriate attention to quality control and assay validation, in spite of the fact that such practice diminishes the reliability of results. External and internal amplification controls became obligatory elements of good laboratory practice in different PCR assays. We propose the inclusion of an analogous approach as a quality control system for multi-template PCR applications. The amplification controls must take into account the characteristics of multi-template PCR and be able to effectively monitor particular assay performance. This study demonstrated the efficiency of a model mixed template as an adequate external amplification control for a particular PCR application. The conditions of multi-template PCR do not allow implementation of a classic internal control; therefore we developed a convenient semi-internal control as an acceptable alternative. In order to evaluate the effects of inhibitors, a model multi-template mix was amplified in a mixture with DNAse-treated sample. Semi-internal control allowed establishment of intervals for robust PCR performance for different samples, thus enabling correct comparison of the samples. The complexity of the external and semi-internal amplification controls must be comparable with the assumed complexity of the samples. We also emphasize that amplification controls should be applied in multi-template PCR regardless of the post-assay method used to analyze products.

An objective, rapid and reproducible method for scoring AFLP peak-height data that minimizes genotyping error

Amplified fragment length polymorphism (AFLP) fingerprint data are now commonly collected using DNA sequencers. AFLP genotypes are still often scored by eye from such data - a time-consuming, error-prone and subjective process. We present a semi-automated method of genotyping sequencer-collected AFLPs at predefined fragment locations (loci) within the fingerprint. Our method uses thresholds of AFLP-polymerase chain reaction-product fluorescence intensity (peak height) in order to: (i) exclude AFLP loci that are likely to contribute high rates of error to data sets, and (ii) determine the AFLP phenotype (fragment presence or absence) at the retained loci. Error rate analysis is an integral part of this process and is used to determine optimal thresholds that minimize genotyping error, while maximizing the numbers of retained loci. We show that application of this method to a large AFLP data set allows genotype calls that are rapid, objective and repeatable, facilitating the extraction of reliable genotype data for molecular ecological studies.

Intragenomic Profiling Using Multicopy Genes: The rDNA Internal Transcribed Spacer Sequences of the Freshwater Sponge Ephydatia fluviatilis

Multicopy genes, like ribosomal RNA genes (rDNA), are widely used to describe and distinguish individuals. Despite concerted evolution that homogenizes a large number of rDNA gene copies, the presence of different gene variants within a genome has been reported. Characterization of an organism by defining every single variant of tens to thousands of rDNA repeat units present in a eukaryotic genome would be quite unreasonable. Here we provide an alternative approach for the characterization of a set of internal transcribed spacer sequences found within every rDNA repeat unit by implementing direct sequencing methodology. The prominent allelic variants and their relative amounts characterizing an individual can be described by a single sequencing electropherogram of the mixed amplicon containing the variants present within the genome. We propose a method for rational analysis of heterogeneity of multicopy genes by compiling a profile based on quantification of different sequence variants of the internal transcribed spacers of the freshwater sponge Ephydatia fluviatilis as an example. In addition to using conventional substitution analysis, we have developed a mathematical method, the proportion model method, to quantify the relative amounts of allelic variants of different length using data from direct sequencing of the heterogeneous amplicon. This method is based on determining the expected signal intensity values (corresponding to peak heights from the sequencing electropherogram) by sequencing clones from the same or highly similar amplicon and comparing hypothesized combinations against the values obtained by direct sequencing of the heterogeneous amplicon. This method allowed to differentiate between all specimens analysed.

Quality control in multi-tag pyrosequencing of microbial communities

Multi-tag pyrosequencing has become a key method in the analysis of microbial community composition. However, it is well known that kinetic bias during the initial PCR amplification of such microbial communities can dramatically distort amplicon abundance prior to downstream emulsion PCR and pyrosequencing. Here we present a simple protocol combining length-heterogeneity PCR fingerprinting with pyrosequencing to ensure the linearity of microbial community amplification. The method employs a fluorescently labeled reverse primer along with multi-tagged forward primers to initially amplify the microbial community. The resulting labeled amplicons are then fingerprinted, purified, and quantitated prior to emulsion PCR and pyrosequencing. Our data demonstrates: (i) use of this protocol results in a distribution of sequences showing linear amplification following emulsion PCR when compared with the initial length-heterogeneity PCR fingerprints, and (ii) that the added tags and labels do not have a negative effect on overall microbial community profiles.

Effectiveness of annealing blocking primers versus restriction enzymes for characterization of generalist diets: unexpected prey revealed in the gut contents of two coral reef fish species

Characterization of predator-prey interactions is challenging as researchers have to rely on indirect methods that can be costly, biased and too imprecise to elucidate the complexity of food webs. DNA amplification and sequencing techniques of gut and fecal contents are promising approaches, but their success largely depends on the ability to amplify the taxonomic array of prey consumed and then match prey amplicons with reference sequences. When little a priori information on diet is available or a generalist predator is targeted, versatile primer sets (also referred to as universal or general primers) as opposed to group- or species-specific primer sets are the most powerful to unveil the full range of prey consumed. However, versatile primers are likely to preferentially amplify the predominant, less degraded predator DNA if no manipulation is performed to exclude this confounding DNA template. In this study we compare two approaches that eliminate the confounding predator template: restriction digestion and the use of annealing blocking primers. First, we use a preliminary DNA barcode library provided by the Moorea BIOCODE project to 1) evaluate the cutting frequency of commercially available restriction enzymes and 2) design predator specific annealing blocking primers. We then compare the performance of the two predator removal strategies for the detection of prey templates using two versatile primer sets from the gut contents of two generalist coral reef fish species sampled in Moorea. Our study demonstrates that blocking primers should be preferentially used over restriction digestion for predator DNA removal as they recover greater prey diversity. We also emphasize that a combination of versatile primers may be required to best represent the breadth of a generalist's diet.

Assignment of SNP allelic configuration in polyploids using competitive allele-specific PCR: application to citrus triploid progeny

BACKGROUND

Polyploidy is a major component of eukaryote evolution. Estimation of allele copy numbers for molecular markers has long been considered a challenge for polyploid species, while this process is essential for most genetic research. With the increasing availability and whole-genome coverage of single nucleotide polymorphism (SNP) markers, it is essential to implement a versatile SNP genotyping method to assign allelic configuration efficiently in polyploids.

SCOPE

This work evaluates the usefulness of the KASPar method, based on competitive allele-specific PCR, for the assignment of SNP allelic configuration. Citrus was chosen as a model because of its economic importance, the ongoing worldwide polyploidy manipulation projects for cultivar and rootstock breeding, and the increasing availability of SNP markers.

CONCLUSIONS

Fifteen SNP markers were successfully designed that produced clear allele signals that were in agreement with previous genotyping results at the diploid level. The analysis of DNA mixes between two haploid lines (Clementine and pummelo) at 13 different ratios revealed a very high correlation (average = 0·9796; s.d. = 0·0094) between the allele ratio and two parameters [θ angle = tan(-1) (y/x) and y' = y/(x + y)] derived from the two normalized allele signals (x and y) provided by KASPar. Separated cluster analysis and analysis of variance (ANOVA) from mixed DNA simulating triploid and tetraploid hybrids provided 99·71 % correct allelic configuration. Moreover, triploid populations arising from 2n gametes and interploid crosses were easily genotyped and provided useful genetic information. This work demonstrates that the KASPar SNP genotyping technique is an efficient way to assign heterozygous allelic configurations within polyploid populations. This method is accurate, simple and cost-effective. Moreover, it may be useful for quantitative studies, such as relative allele-specific expression analysis and bulk segregant analysis.

Ultradeep sequencing detects GNAQ and GNA11 mutations in cell-free DNA from plasma of patients with uveal melanoma

Elevated levels of cell-free DNA (cfDNA) are frequently observed in tumor patients. Activating mutations in exon 4 (R183) and exon 5 (Q209) of GNAQ and GNA11 are almost exclusively found in uveal melanoma, thus providing a highly specific marker for the presence of circulating tumor DNA (ctDNA). To establish a reliable, noninvasive assay that might allow early detection and monitoring of metastatic disease, we determined the proportion of GNAQ or GNA11 mutant reads in cfDNA of uveal melanoma patients by ultradeep sequencing. Cell-free DNA from 28 uveal melanoma patients with metastases or extraocular growth was isolated and quantified by real-time polymerase chain reaction (PCR) (7–1550 ng DNA/mL plasma). GNAQ and GNA11 regions of interest were amplified in 22 of 28 patients and ultradeep sequencing of amplicons was performed to detect even low proportions of mutant reads. We detected Q209 mutations (2–38% mutant reads) in either GNAQ or GNA11 in the plasma of 9 of 22 metastasized patients. No correlation between the proportion of mutant reads and the concentration of cfDNA could be detected. Among the nine ctDNA-positive patients, four had metastases in bone, whereas no metastases were detected in the 13 ctDNA-negative patients at this location (P = 0.025). Furthermore, ctDNA-positive patients tended to be younger at initial diagnosis and show larger metastases. The results show that ultradeep amplicon sequencing can be used to detect tumor DNA in plasma of metastasized uveal melanoma patients. It remains to be shown if this approach can be used for early detection of disseminated tumor disease.

Reconstruction of phylogenetic relationships in a highly reticulate group with deep coalescence and recent speciation (Hieracium, Asteraceae)

Phylogeny reconstruction based on multiple unlinked markers is often hampered by incongruent gene trees, especially in closely related species complexes with high degrees of hybridization and polyploidy. To investigate the particular strengths and limitations of chloroplast DNA (cpDNA), low-copy nuclear and multicopy nuclear markers for elucidating the evolutionary history of such groups, we focus on Hieracium s.str., a predominantly apomictic genus combining the above-mentioned features. Sequences of the trnV-ndhC and trnT-trnL intergenic spacers were combined for phylogenetic analyses of cpDNA. Part of the highly variable gene for squalene synthase (sqs) was applied as a low-copy nuclear marker. Both gene trees were compared with previous results based on the multicopy external transcribed spacer (ETS) of the nuclear ribosomal DNA. The power of the different markers to detect hybridization varied, but they largely agreed on particular hybrid and allopolyploid origins. The same crown groups of species were recognizable in each dataset, but basal relationships were strongly incongruent among cpDNA, sqs and ETS trees. The ETS tree was considered as the best approximation of the species tree. Both cpDNA and sqs trees showed basal polytomies as well as merging or splitting of species groups of non-hybrid taxa. These patterns can be best explained by a rapid diversification of the genus with ancestral polymorphism and incomplete lineage sorting. A hypothetical scenario of Hieracium speciation based on all available (including non-molecular) evidence is depicted. Incorporation of seemingly contradictory information helped to better understand species origins and evolutionary patterns in this notoriously difficult agamic complex.

Reconstructing the complex evolutionary origin of wild allopolyploid tobaccos (Nicotiana section suaveolentes)

Nicotiana (Solanaceae) provides an ideal system for understanding polyploidization, a pervasive and powerful evolutionary force in plants, as this genus contains several groups of allotetraploids that formed at different times from different diploid progenitors. However, the parental lineages of the largest group of allotetraploids, Nicotiana section Suaveolentes, have been problematic to identify. Using data from four regions of three low-copy nuclear genes, nuclear ribosomal DNA, and regions of the plastid genome, we have reconstructed the evolutionary origin of sect. Suaveolentes and identified the most likely diploid progenitors by using a combination of gene trees and network approaches to uncover the most strongly supported evidence of species relationships. Our analyses best support a scenario where a member of the sect. Sylvestres lineage acted as the paternal progenitor and a member of either sect. Petunioides or sect. Noctiflorae that also contained introgressed DNA from the other, or a hypothetical hybrid species between these two sections, was the maternal progenitor. Nicotiana exemplifies many of the factors that can complicate the reconstruction of polyploid evolutionary history and highlights how reticulate evolution at the diploid level can add even greater complexity to allopolyploid genomes.

Ocean viruses: rigorously evaluating the metagenomic sample-to-sequence pipeline

As new environments are studied, viruses consistently emerge as important and prominent players in natural and man-made ecosystems. However, much of what we know is built both upon the foundation of the culturable minority and using methods that are often insufficiently ground-truthed. Here, we review the modern culture-independent viral metagenomic sample-to-sequence pipeline and how next-generation sequencing techniques are drastically altering our ability to systematically and rigorously evaluate them. Together, a series of studies quantitatively evaluate existing and new methods that allow-even for ultra-low DNA samples-the generation of replicable, near-quantitative datasets that maximize inter-comparability and biological inference.

Tracking temporal changes of bacterial community fingerprints during the initial stages of composting

Abstract The initial phase of composting is the most dynamic part of the process and is characterized by rapid increases in temperature, large swings in pH, and the degradation of simple organic compounds. DNA samples were taken from an active compost system to determine the microbial 16S rRNA gene sequences that were present during this phase. We observed two significant shifts in the composition of the microbial community, one between 12 and 24 h and the other between 60 and 72 h into the process using automated 16S-23S rRNA intergenic spacer amplification (ARISA). The 16S rRNA gene sequences adjoining the most common ARISA fragments at each time point were determined. We found that sequences related to lactic acid bacteria were most common during the first 60 h and Bacillus-type sequences were most common between 72 and 96 h. While the temperature increased steadily over the first 96 h, the pH dropped after 12 h and increased after 60 h correlating with the shift from Bacillus to lactic acid sequences and the later return to Bacillus-type sequences.

Extreme multiple reticulate origins of the Pteris cadieri complex (Pteridaceae)

The Pteris cadieri complex displays extensive morphological variation and seems to have originated through hybridization. However, the members of this complex reproduce by apogamy, which usually limits genetic variation. To evaluate the hypotheses of hybrid origins, the pattern of evolution in this species complex is reconstructed. Multiple methodologies were used. Diploids, triploids, and tetraploids were identified by chromosome counts and flow cytometry. Nuclear DNA markers (cytosolic phosphoglucose isomerase gene, PgiC) were used, together with chloroplast DNA markers (atpB-rbcL spacer and rbcL gene) to infer the biparental and maternal lineages of the Pteris cadieri complex. The three cpDNA haplotype groups and five PgiC alleles found in this study indicate that the evolution of the Pteris cadieri complex has been extremely reticulate. Up to 11 taxa belonging to eight morphs were identified. By comparing genetic variation in the Pteris cadieri in two independent areas, Hainan and Taiwan, we inferred that hybridization has occurred independently in different areas. Furthermore, we found evidence for phenological divergence (evergreen and deciduous) within Taiwan. We propose that the Pteris cadieri complex originated from different genetic lineages through multiple hybridizations in different geographical areas, leading to its present morphological diversity.

DNA barcoding of lichenized fungi demonstrates high identification success in a floristic context

• Efforts are currently underway to establish a standard DNA barcode region for fungi; we tested the utility of the internal transcribed spacer (ITS) of nuclear ribosomal DNA for DNA barcoding in lichen-forming fungi by sampling diverse species across eight orders. • Amplification of the ITS region (ITS1-5.8S-ITS2) was conducted for 351 samples, encompassing 107, 55 and 28 species, genera and families, respectively, of lichenized fungi. We assessed the ability of the entire ITS vs the ITS2 alone to discriminate between species in a taxonomic dataset (members of the genus Usnea) and a floristic dataset. • In the floristic dataset, 96.3% of sequenced samples could be assigned to the correct species using ITS or ITS2; a barcode gap for ITS is present in 92.1% of species. Although fewer species have a barcode gap in the taxonomic dataset (73.3% with ITS and 68.8% with ITS2), up to 94.1% of samples were assigned to the correct species using BLAST. • While discrimination between the most closely related species will remain challenging, our results demonstrate the potential to identify a high percentage of specimens to the correct species, and the remainder to the correct genus, when using DNA barcoding in a floristic context.