Assessment of Species Diversity and Distribution of an Ancient Diatom Lineage Using a DNA Metabarcoding Approach

Background

Continuous efforts to estimate actual diversity and to trace the species distribution and ranges in the natural environments have gone in equal pace with advancements of the technologies in the study of microbial species diversity from microscopic observations to DNA-based barcoding. DNA metabarcoding based on Next Generation Sequencing (NGS) constitutes the latest advancement in these efforts. Here we use NGS data from different sites to investigate the geographic range of six species of the diatom family Leptocylindraceae and to identify possible new taxa within the family.

Methodology/Principal Findings

We analysed the V4 and V9 regions of the nuclear-encoded SSU rDNA gene region in the NGS database of the European ERA-Biodiversa project BioMarKs, collected in plankton and sediments at six coastal sites in European coastal waters, as well as environmental sequences from the NCBI database. All species known in the family Leptocylindraceae were detected in both datasets, but the much larger Illumina V9 dataset showed a higher species coverage at the various sites than the 454 V4 dataset. Sequences identical or similar to the references of Leptocylindrus aporus, L. convexus, L. danicus/hargravesii and Tenuicylindrus belgicus were found in the Mediterranean Sea, North Atlantic Ocean and Black Sea as well as at locations outside Europe. Instead, sequences identical or close to that of L. minimus were found in the North Atlantic Ocean and the Black Sea but not in the Mediterranean Sea, while sequences belonging to a yet undescribed taxon were encountered only in Oslo Fjord and Baffin Bay.

Conclusions/Significance

Identification of Leptocylindraceae species in NGS datasets has expanded our knowledge of the species biogeographic distribution and of the overall diversity of this diatom family. Individual species appear to be widespread, but not all of them are found everywhere. Despite the sequencing depth allowed by NGS and the wide geographic area covered by this study, the diversity of this ancient diatom family appears to be low, at least at the level of the marker used in this study.     

Natural Selection Constrains Neutral Diversity across A Wide Range of Species

The neutral theory of molecular evolution predicts that the amount of neutral polymorphisms within a species will increase proportionally with the census population size (Nc). However, this prediction has not been borne out in practice: while the range of Nc spans many orders of magnitude, levels of genetic diversity within species fall in a comparatively narrow range. Although theoretical arguments have invoked the increased efficacy of natural selection in larger populations to explain this discrepancy, few direct empirical tests of this hypothesis have been conducted. In this work, we provide a direct test of this hypothesis using population genomic data from a wide range of taxonomically diverse species. To do this, we relied on the fact that the impact of natural selection on linked neutral diversity depends on the local recombinational environment. In regions of relatively low recombination, selected variants affect more neutral sites through linkage, and the resulting correlation between recombination and polymorphism allows a quantitative assessment of the magnitude of the impact of selection on linked neutral diversity. By comparing whole genome polymorphism data and genetic maps using a coalescent modeling framework, we estimate the degree to which natural selection reduces linked neutral diversity for 40 species of obligately sexual eukaryotes. We then show that the magnitude of the impact of natural selection is positively correlated with Nc, based on body size and species range as proxies for census population size. These results demonstrate that natural selection removes more variation at linked neutral sites in species with large Nc than those with small Nc and provides direct empirical evidence that natural selection constrains levels of neutral genetic diversity across many species. This implies that natural selection may provide an explanation for this longstanding paradox of population genetics.     

Ebola Virus Genome Plasticity as a Marker of Its Passaging History: A Comparison of In Vitro Passaging to Non-Human Primate Infection

To identify polymorphic sites that could be used as biomarkers of Ebola virus passage history, we repeatedly amplified Ebola virus (Kikwit variant) in vitro and in vivo and performed deep sequencing analysis of the complete genomes of the viral subpopulations. We then determined the sites undergoing selection during passage in Vero E6 cells. Four locations within the Ebola virus Kikwit genome were identified that together segregate cell culture-passaged virus and virus obtained from infected non-human primates. Three of the identified sites are located within the glycoprotein gene (GP) sequence: the poly-U (RNA editing) site at position 6925, as well as positions 6677, and 6179. One site was found in the VP24 gene at position 10833. In all cases, in vitro and in vivo, both populations (majority and minority variants) were maintained in the viral swarm, with rapid selections occurring after a few passages or infections. This analysis approach will be useful to differentiate whether filovirus stocks with unknown history have been passaged in cell culture and may support filovirus stock standardization for medical countermeasure development.     

Plasticity in Dendroclimatic Response across the Distribution Range of Aleppo Pine (Pinus halepensis)

We investigated the variability of the climate-growth relationship of Aleppo pine across its distribution range in the Mediterranean Basin. We constructed a network of tree-ring index chronologies from 63 sites across the region. Correlation function analysis identified the relationships of tree-ring index to climate factors for each site. We also estimated the dominant climatic gradients of the region using principal component analysis of monthly, seasonal, and annual mean temperature and total precipitation from 1,068 climatic gridpoints. Variation in ring width index was primarily related to precipitation and secondarily to temperature. However, we found that the dendroclimatic relationship depended on the position of the site along the climatic gradient. In the southern part of the distribution range, where temperature was generally higher and precipitation lower than the regional average, reduced growth was also associated with warm and dry conditions. In the northern part, where the average temperature was lower and the precipitation more abundant than the regional average, reduced growth was associated with cool conditions. Thus, our study highlights the substantial plasticity of Aleppo pine in response to different climatic conditions. These results do not resolve the source of response variability as being due to either genetic variation in provenance, to phenotypic plasticity, or a combination of factors. However, as current growth responses to inter-annual climate variability vary spatially across existing climate gradients, future climate-growth relationships will also likely be determined by differential adaptation and/or acclimation responses to spatial climatic variation. The contribution of local adaptation and/or phenotypic plasticity across populations to the persistence of species under global warming could be decisive for prediction of climate change impacts across populations. In this sense, a more complex forest dynamics modeling approach that includes the contribution of genetic variation and phenotypic plasticity can improve the reliability of the ecological inferences derived from the climate-growth relationships.     

环境DNA宏条形码在生物多样性研究与监测中的应用

环境DNA宏条形码(eDNA metabarcoding)技术通过提取水体、土壤、空气中的环境DNA,使用引物PCR扩增与高通量测序,进行物种鉴定与生物多样性评估.作为一种新的监测技术,相比于传统监测技术更加快捷、准确以及对自然环境的破坏小,因此在一定程度上改变了我们调查地球生物多样性的方式.本文综述了环境DNA宏条形...     

Multi-Scales Analysis of Primate Diversity and Protected Areas at a Megadiverse Region

In this paper, we address the question of what proportion of biodiversity is represented within protected areas. We assessed the effectiveness of different protected area types at multiple scales in representing primate biodiversity in the Brazilian Legal Amazon. We used point locality data and distribution data for primate species within 1°, 0.5°, and 0.25° spatial resolution grids, and computed the area of reserves within each cell. Four different approaches were used – no reserves (A), exclusively strict use reserves (B), strict and sustainable use reserves (C), and strict and sustainable use reserves and indigenous lands (D). We used the complementarity concept to select reserve networks. The proportions of cells that were classified as reserves at a grid resolution of 1° were 37%, 64%, and 88% for approaches B, C and D, respectively. Our comparison of these approaches clearly showed the effect of an increase in area on species representation. Representation was consistently higher at coarser resolutions, indicating the effect of grain size. The high number of irreplaceable cells for selected networks identified based on approach A could be attributed to the use of point locality occurrence data. Although the limited number of point occurrences for some species may have been due to a Wallacean shortfall, in some cases it may also be the result of an actual restricted geographic distribution. The existing reserve system cannot be ignored, as it has an established structure, legal protection status, and societal recognition, and undoubtedly represents important elements of biodiversity. However, we found that strict use reserves (which are exclusively dedicated to biodiversity conservation) did not effectively represent primate species. This finding may be related to historical criteria for selecting reserves based on political, economic, or social motives.     

Primate evolution of a dispersed human repetitive DNA sequence

A dispersed middle repetitive DNA sequence isolated originally from human chromosome 12 did not show homology with rodent DNA under standard conditions of Southern DNA blot analysis. The evolutionary relationship of this human repetitive DNA to that of other primates was investigated using three hybridization methods: DNA dot blot, Southern DNA blot analysis, and chromosome in situ hybridization. Homology with the human repetitive DNA was found throughout the suborder Anthropoidea, in fourteen ape and New and Old World monkey species. In addition, the human pattern of hybridization to noncentromeric regions of all chromosomes was seen. No hybridization by any of the three techniques was found in five species of the suborder Prosimii. The phenomenon of marked differences in sequence homology and copy number of dispersed repetitive DNA from closely related species has been observed in protozoans (Plasmodia), Drosophila, sea urchins, mice and the great apes (Hominoidea). We report here a similar phenomenon that may have occurred at an early stage in primate evolution.     

Primate evolution at the DNA level and a classification of hominoids

Summary The genetic distances among primate lineages estimated from orthologous noncoding nucleotide sequences of -type globin loci and their flanking and intergenic DNA agree closely with the distances (delta T₅₀H values) estimated by cross hybridization of total genomic single-copy DNAs. These DNA distances and the maximum parsimony tree constructed for the nucleotide sequence orthologues depict a branching pattern of primate lineages that is essentially congruent with the picture from phylogenetic analyses of morphological characters. The molecular evidence, however, resolves ambiguities in the morphological picture and provides an objective view of the cladistic position of humans among the primates. The molecular data group humans with chimpanzees in subtribe Hominina, with gorillas in tribe Hominini, orangutans in subfamily Homininae, gibbons in family Hominidae, Old World monkeys in infraorder Catarrhini, New World monkeys in semisuborder Anthropoidea, tarsiers in suborder Haplorhini, and strepsirhines (lemuriforms and lorisiforms) in order Primates. A seeming incongruency between organismal and molecular levels of evolution, namely that morphological evolution appears to have speeded up in higher primates, especially in the lineage to humans, while molecular evolution has slowed down, may have the trivial explanation that relatively small genetic changes may sometimes result in marked phenotypic changes.     

Benefits and Limitations of DNA Barcoding and Metabarcoding in Herbal Product Authentication

Introduction

Herbal medicines play an important role globally in the health care sector and in industrialised countries they are often considered as an alternative to mono‐substance medicines. Current quality and authentication assessment methods rely mainly on morphology and analytical phytochemistry‐based methods detailed in pharmacopoeias. Herbal products however are often highly processed with numerous ingredients, and even if these analytical methods are accurate for quality control of specific lead or marker compounds, they are of limited suitability for the authentication of biological ingredients.

Objective

To review the benefits and limitations of DNA barcoding and metabarcoding in complementing current herbal product authentication.

Method

Recent literature relating to DNA based authentication of medicinal plants, herbal medicines and products are summarised to provide a basic understanding of how DNA barcoding and metabarcoding can be applied to this field.

Results

Different methods of quality control and authentication have varying resolution and usefulness along the value chain of these products. DNA barcoding can be used for authenticating products based on single herbal ingredients and DNA metabarcoding for assessment of species diversity in processed products, and both methods should be used in combination with appropriate hyphenated chemical methods for quality control.

Conclusions

DNA barcoding and metabarcoding have potential in the context of quality control of both well and poorly regulated supply systems. Standardisation of protocols for DNA barcoding and DNA sequence‐based identification are necessary before DNA‐based biological methods can be implemented as routine analytical approaches and approved by the competent authorities for use in regulated procedures. © 2017 The Authors. Phytochemical Analysis Published by John Wiley & Sons Ltd.     

Broad Anatomical Variation within a Narrow Wood Density Range—A Study of Twig Wood across 69 Australian Angiosperms

ObjectivesJust as people with the same weight can have different body builds, woods with the same wood density can have different anatomies. Here, our aim was to assess the magnitude of anatomical variation within a restricted range of wood density and explore its potential ecological implications.MethodsTwig wood of 69 angiosperm tree and shrub species was analyzed. Species were selected so that wood density varied within a relatively narrow range (0.38–0.62 g cm^-3). Anatomical traits quantified included wood tissue fractions (fibres, axial parenchyma, ray parenchyma, vessels, and conduits with maximum lumen diameter below 15 μm), vessel properties, and pith area. To search for potential ecological correlates of anatomical variation the species were sampled across rainfall and temperature contrasts, and several other ecologically-relevant traits were measured (plant height, leaf area to sapwood area ratio, and modulus of elasticity).ResultsDespite the limited range in wood density, substantial anatomical variation was observed. Total parenchyma fraction varied from 0.12 to 0.66 and fibre fraction from 0.20 to 0.74, and these two traits were strongly inversely correlated (r = -0.86, P < 0.001). Parenchyma was weakly (0.24 ≤|r|≤ 0.35, P < 0.05) or not associated with vessel properties nor with height, leaf area to sapwood area ratio, and modulus of elasticity (0.24 ≤|r|≤ 0.41, P < 0.05). However, vessel traits were fairly well correlated with height and leaf area to sapwood area ratio (0.47 ≤|r|≤ 0.65, all P < 0.001). Modulus of elasticity was mainly driven by fibre wall plus vessel wall fraction rather than by the parenchyma component.ConclusionsOverall, there seem to be at least three axes of variation in xylem, substantially independent of each other: a wood density spectrum, a fibre-parenchyma spectrum, and a vessel area spectrum. The fibre-parenchyma spectrum does not yet have any clear or convincing ecological interpretation.     

Variable Transforming Activity of Nonglucosylated and Fragmented T4 DNA

Nonglucosylated and fragmented T4 DNA shows a gene-specific variation in transformation efficiency.     

Validation of Next-Generation Sequencing of Entire Mitochondrial Genomes and the Diversity of Mitochondrial DNA Mutations in Oral Squamous Cell Carcinoma

Background

Oral squamous cell carcinoma (OSCC) is mainly caused by smoking and alcohol abuse and shows a five-year survival rate of ~50%. We aimed to explore the variation of somatic mitochondrial DNA (mtDNA) mutations in primary oral tumors, recurrences and metastases.

Methods

We performed an in-depth validation of mtDNA next-generation sequencing (NGS) on an Illumina HiSeq 2500 platform for its application to cancer tissues, with the goal to detect low-level heteroplasmies and to avoid artifacts. Therefore we genotyped the mitochondrial genome (16.6 kb) from 85 tissue samples (tumors, recurrences, resection edges, metastases and blood) collected from 28 prospectively recruited OSCC patients applying both Sanger sequencing and high-coverage NGS (~35,000 reads per base).

Results

We observed a strong correlation between Sanger sequencing and NGS in estimating the mixture ratio of heteroplasmies (r = 0.99; p<0.001). Non-synonymous heteroplasmic variants were enriched among cancerous tissues. The proportions of somatic and inherited variants in a given gene region were strongly correlated (r = 0.85; p<0.001). Half of the patients shared mutations between benign and cancerous tissue samples. Low level heteroplasmies (<10%) were more frequent in benign samples compared to tumor samples, where heteroplasmies >10% were predominant. Four out of six patients who developed a local tumor recurrence showed mutations in the recurrence that had also been observed in the primary tumor. Three out of five patients, who had tumor metastases in the lymph nodes of their necks, shared mtDNA mutations between primary tumors and lymph node metastases. The percentage of mutation heteroplasmy increased from the primary tumor to lymph node metastases.

Conclusions

We conclude that Sanger sequencing is valid for heteroplasmy quantification for heteroplasmies ≥10% and that NGS is capable of reliably detecting and quantifying heteroplasmies down to the 1%-level. The finding of shared mutations between primary tumors, recurrences and metastasis indicates a clonal origin of malignant cells in oral cancer.     

Plasticity of Lythrum salicaria and Phragmites australis growth characteristics across a European geographical gradient

Plants of Lythrum salicaria and Phragmites australis originating from localities across the European north–south geographical gradient were cultivated in parallel in an outdoor tub experiment. A strong correlation was found between growth and morphometric characteristics related to plant size (plant height, basal diameter, aboveground- and belowground plant biomass, etc.) and the position of the respective populations along the north–south gradient. Plants of both L. salicaria and P. australis from the southern localities grew taller and more vigorously and flowered later than plants from relatively more northern localities. From this point of view, the plants originating from south European populations were comparable to invasive North American plants. Our study indicates that explanation of the competitive success of populations invading new geographical areas may involve the role of geographic gradients within the species native range.     

High Genetic Diversity and Structured Populations of the Oriental Fruit Moth in Its Range of Origin

The oriental fruit moth Grapholita ( = Cydia) molesta is a key fruit pest globally. Despite its economic importance, little is known about its population genetics in its putative native range that includes China. We used five polymorphic microsatellite loci and two mitochondrial gene sequences to characterize the population genetic diversity and genetic structure of G. molesta from nine sublocations in three regions of a major fruit growing area of China. Larval samples were collected throughout the season from peach, and in late season, after host switch by the moth to pome fruit, also from apple and pear. We found high numbers of microsatellite alleles and mitochondrial DNA haplotypes in all regions, together with a high number of private alleles and of haplotypes at all sublocations, providing strong evidence that the sampled area belongs to the origin of this species. Samples collected from peach at all sublocations were geographically structured, and a significant albeit weak pattern of isolation-by-distance was found among populations, likely reflecting the low flight capacity of this moth. Interestingly, populations sampled from apple and pear in the late season showed a structure differing from that of populations sampled from peach throughout the season, indicating a selective host switch of a certain part of the population only. The recently detected various olfactory genotypes in G. molesta may underly this selective host switch. These genetic data yield, for the first time, an understanding of population dynamics of G. molesta in its native range, and of a selective host switch from peach to pome fruit, which may have a broad applicability to other global fruit production areas for designing suitable pest management strategies.     

Diversity of hopanoids and squalene-hopene cyclases across a tropical land-sea gradient

Bacterial hopanoids are ubiquitous in Earth surface environments. They hold promise as environmental and ecological biomarkers, if the phylogeny and physiological drivers of hopanoid biosynthesis can be linked with the distribution of hopanoids observed across a breadth of samples. Here we survey the diversity of hopanoid cyclases from a land‐sea gradient across the island of San Salvador, in the easternmost part of the Bahamas. The distribution of lipids was determined for the same sites, for the first time overlaying quantification of bacteriohopanepolyols with sqhC phylogeny. The results are similar to previous reports: environmental sqhCs average < 65% translated amino acid identity to their closest named relatives, and sequences from putative Proteobacteria dominate. Additionally, a new and apparently ubiquitous group of marine hopanoid producers is identified; it has no identifiable close relatives. The greatest diversity of hopanoid lipids occurs in soil, but hopanoids represent a minor fraction of total soil‐derived lipids. Marine samples contain fewer identifiable hopanoids, but they are more abundant as a fraction of the total extractable lipids. In soil, the dominant compounds are 35‐aminobacteriohopane‐32,33,34‐triol and adenosylhopane. In an upper estuarine sample, bacteriohopanetetrol and 32,35‐anhydrobacteriohopanetetrol dominate; while in lower estuarine and open marine samples, the most abundant are bacteriohopanetetrol and bacteriohopaneribonolactone. Cyclitol ethers are trace components in the soil, absent in the estuary, and of moderate abundance in the open marine setting, suggesting a dominant marine source. Conversely, aminotriol and aminotetrol decrease in abundance or disappear completely from land to ocean, while 2‐methyldiplopterol shows the opposite trend. Small quantities of 2‐methylbacteriohopanepolyols are detectable in all samples. The overall hopanoid distributions may correlate to the major phylogenetic families of hopanoid producers or to the environments in which they are found.     

High Cryptic Diversity across the Global Range of the Migratory Planktonic Copepods Pleuromamma piseki and P. gracilis

Although holoplankton are ocean drifters and exhibit high dispersal potential, a number of studies on single species are finding highly divergent genetic clades. These cryptic species complexes are important to discover and describe, as identification of common marine species is fundamental to understanding ecosystem dynamics. Here we investigate the global diversity within Pleuromamma piseki and P. gracilis, two dominant members of the migratory zooplankton assemblage in subtropical and tropical waters worldwide. Using DNA sequence data from the mitochondrial gene cytochrome c oxidase subunit II (mtCOII) from 522 specimens collected across the Pacific, Atlantic and Indian Oceans, we discover twelve well-resolved genetically distinct clades in this species complex (Bayesian posterior probabilities >0.7; 6.3–17% genetic divergence between clades). The morphologically described species P. piseki and P. gracilis did not form monophyletic groups, rather they were distributed throughout the phylogeny and sometimes co-occurred within well-resolved clades: this result suggests that morphological characters currently used for taxonomic identification of P. gracilis and P. piseki may be inaccurate as indicators of species’ boundaries. Cryptic clades within the species complex ranged from being common to rare, and from cosmopolitan to highly restricted in distribution across the global ocean. These novel lineages appear to be ecologically divergent, with distinct biogeographic distributions across varied pelagic habitats. We hypothesize that these mtDNA lineages are distinct species and suggest that resolving their systematic status is important, given the ecological significance of the genus Pleuromamma in subtropical-tropical waters worldwide.     

Dietary Plasticity of Generalist and Specialist Ungulates in the Namibian Desert: A Stable Isotopes Approach

Desert ungulates live in adverse ecosystems that are particularly sensitive to degradation and global climate change. Here, we asked how two ungulate species with contrasting feeding habits, grazing gemsbok (Oryx g. gazella) and browsing springbok (Antidorcas marsupialis), respond to an increase in food availability during a pronounced rain period. We used a stable isotope approach to delineate the feeding habits of these two ungulates in the arid Kunene Region of Namibia. Our nineteen months field investigation included two time periods of drought when food availability for ungulates was lowest and an intermediate period with extreme, unusual rainfalls. We documented thirteen isotopically distinct food sources in the isotopic space of the study area. Our results indicated a relatively high dietary plasticity of gemsbok, which fed on a mixture of plants, including more than 30% of C3 plants during drought periods, but almost exclusively on C4 and CAM plant types when food was plentiful. During drought periods, the inferred gemsbok diets also consisted of up to 25% of Euphorbia damarana; an endemic CAM plant that is rich in toxic secondary plant compounds. In contrast, springbok were generalists, feeding on a higher proportion of C3 than C4/CAM plants, irrespective of environmental conditions. Our results illustrate two dietary strategies in gemsbok and springbok which enable them to survive and coexist in the hostile Kunene arid ecosystem.