High-throughput sequencing on preservative ethanol is effective at jointly examining infraspecific and taxonomic diversity, although bioinformatics pipelines do not perform equally.

While recent sequencing technologies (third generation sequencing) can successfully sequence all copies of nuclear ribosomal DNA (rDNA) markers present within a genome and offer insights into the intragenomic variation of these markers, high intragenomic variation can be a source of confusion for high-throughput species identification using such technologies. High-throughput (HT) amplicon sequencing via PacBio SEQUEL I was used to evaluate the intragenomic variation of the ITS region and D1–D2 LSU domains in nine Cordyceps species, and the accuracy of such technology to identify these species based on molecular phylogenies was also assessed. PacBio sequences within strains showed variable level of intragenomic variation among the studied Cordyceps species with C. blackwelliae showing greater variation than the others. Some variants from a mix of species clustered together outside their respective species of origin, indicative of intragenomic variation that escaped concerted evolution shared between species. Proper selection of consensus sequences from HT amplicon sequencing is a challenge for interpretation of correct species identification. PacBio consensus sequences with the highest number of reads represent the major variants within a genome and gave the best results in terms of species identification.

Assessing phytoplankton diversity is of primary importance for both basic and applied ecological studies. Following the advances in molecular methods, phytoplankton studies are switching from using classical microscopy to high throughput sequencing approaches. However, methodological comparisons of these approaches have rarely been reported. In this study, we compared the two methods, using a unique dataset of multiple water samples taken from a natural freshwater environment. Environmental DNA was extracted from 300 water samples collected weekly during 20 years, followed by high throughput sequencing of amplicons from the 16S and 18S rRNA hypervariable regions. For each water sample, phytoplankton diversity was also estimated using light microscopy. Our study indicates that species compositions detected by light microscopy and 454 high throughput sequencing do not always match. High throughput sequencing detected more rare species and picoplankton than light microscopy, and thus gave a better assessment of phytoplankton diversity. However, when compared to light microscopy, high throughput sequencing of 16S and 18S rRNA amplicons did not adequately identify phytoplankton at the species level. In summary, our study recommends a combined strategy using both morphological and molecular techniques.

Community genomics: a community-wide perspective on within-species genetic diversity.

Gene encoding the large subunit of As(III) oxidase (AioA), an important component of the microbial As(III) oxidation system, is a widely used biomarker to characterize As(III)-oxidizing communities in the environment. However, many studies were restricted to a few sequences generated by clone libraries and Sanger sequencing, which may have underestimated the diversity of As(III)-oxidizers in natural environments. In this study, we designed a primer pair, 1109F (5′-ATC TGG GGB AAY RAC AAY TA−3′) and 1548R (5′-TTC ATB GAS GTS AGR TTC AT−3′), targeting gene sequence encoding for the conserved molybdopterin center of the AioA protein, yielding amplicons approximately 450 bp in size that are feasible for highly parallel amplicon sequencing. By utilizing in silico analyses and the experimental construction of clone libraries using Sanger sequencing, the specificity and resolution of 1109F/1548R are approximated with two other previously published and commonly used primers, i.e., M1-2F/M3-2R and deg1F/deg1R. With the use of the 1109F/1548R primer pair, the taxonomic composition of the aioA genes was similar both according to the Sanger and next-generation sequencing (NGS) platforms. Furthermore, high-throughput amplicon sequencing using the primer pair, 1109F/1548R, successfully identified the well-known As(III)-oxidizers in paddy soils and sediments, and they also revealed the differences in the community structure and composition of As(III)-oxidizers in above two biotopes. The random forest analysis showed that the dissolved As(III) had the highest relative influence on the Chao1 index of the aioA genes. These observations demonstrate that the newly designed PCR primers enhanced the ability to detect the diversity of aioA-encoding microorganisms in environments using highly parallel short amplicon sequencing.

DNA analysis of predator faeces using high-throughput amplicon sequencing (HTS) enhances our understanding of predator-prey interactions. However, conclusions drawn from this technique are constrained by biases that occur in multiple steps of the HTS workflow. To better characterize insectivorous animal diets, we used DNA from a diverse set of arthropods to assess PCR biases of commonly used and novel primer pairs for the mitochondrial gene, cytochrome oxidase C subunit 1 (COI). We compared diversity recovered from HTS of bat guano samples using a commonly used primer pair "ZBJ" to results using the novel primer pair "ANML." To parameterize our bioinformatics pipeline, we created an arthropod mock community consisting of single-copy (cloned) COI sequences. To examine biases associated with both PCR and HTS, mock community members were combined in equimolar amounts both pre- and post-PCR. We validated our system using guano from bats fed known diets and using composite samples of morphologically identified insects collected in pitfall traps. In PCR tests, the ANML primer pair amplified 58 of 59 arthropod taxa (98%), whereas ZBJ amplified 24-40 of 59 taxa (41%-68%). Furthermore, in an HTS comparison of field-collected samples, the ANML primers detected nearly fourfold more arthropod taxa than the ZBJ primers. The additional arthropods detected include medically and economically relevant insect groups such as mosquitoes. Results revealed biases at both the PCR and sequencing levels, demonstrating the pitfalls associated with using HTS read numbers as proxies for abundance. The use of an arthropod mock community allowed for improved bioinformatics pipeline parameterization.

Comparison of FLT3-ITD Detection By High-Throughput Amplicon Sequencing to Routine Diagnostics - a Retrospective Analysis of AMLCG Study Patients

BackgroundHemophagocytic lymphohistocytosis (HLH) is a rare but fatal hyperinflammatory syndrome caused by uncontrolled proliferation of activated macrophages and T lymphocytes secreting high amounts of inflammatory cytokines. Genetic defect is a common cause of HLH. HLH is complicated to be diagnosed as there are many common symptoms with other disorders.Case presentationHere we report on an HLH case caused by 1 bp deletion in gene SH2D1A. Patient was a 3-years-old boy and had fever for more than 8 days. Splenomegaly and hemophagocytosis in bone marrow were observed in examination. The results of the blood analysis suggested the diagnosis of HLH. Genetic test based on high throughput amplicon sequencing was then conducted by targeting all six known HLH-causing genes simultaneously. It took only one single day to accomplish the amplicon sequencing library preparation, sequencing and data analysis. Finally, a novel 1 bp deletion in gene SH2D1A was discovered. The result was also confirmed by Sanger sequencing. The result of the genetic test served as a good basis for further diagnosis of HLH.ConclusionThis is the first case that the disease-causing genetic defect of HLH was quickly determined by high throughput amplicon sequencing. This diagnosis was also confirmed by Sanger sequencing and cross-validated by blood analysis and other clinical criteria. This case suggests that genetic test based on amplicon sequencing is a powerful tool for diagnosis of HLH and other diseases caused by genetic defect.

DNA analysis of predator feces using high-throughput amplicon sequencing (HTS) enhances our understanding of predator-prey interactions. However, conclusions drawn from this technique are constrained by biases that occur in multiple steps of the HTS workflow. To better characterize insectivorous animal diets, we used DNA from a diverse set of arthropods to assess PCR biases of commonly used and novel primer pairs for the mitochondrial gene, cytochrome oxidase C subunit 1 (CO1). We compared diversity recovered from HTS of bat guano samples using a commonly used primer pair “ZBJ” to results using the novel primer pair “ANML”. To parameterize our bioinformatics pipeline, we created an arthropod mock community consisting of single-copy (cloned) CO1 sequences. To examine biases associated with both PCR and HTS, mock community members were combined in equimolar amounts both pre- and post-PCR. We validated our system using guano from bats fed known diets and using composite samples of morphologically identified insects collected in pitfall traps. In PCR tests, the ANML primer pair amplified 58 of 59 arthropod taxa (98%) whereas ZBJ amplified 24 of 59 taxa (41%). Furthermore, in an HTS comparison of field-collected samples, the ANML primers detected nearly four-fold more arthropod taxa than the ZBJ primers. The additional arthropods detected include medically and economically relevant insect groups such as mosquitoes. Results revealed biases at both the PCR and sequencing levels, demonstrating the pitfalls associated with using HTS read numbers as proxies for abundance. The use of an arthropod mock community allowed for improved bioinformatics pipeline parameterization.

Spatial patterns of biological diversity have been extensively studied in ecology and population genetics, because they reflect the forces acting on biodiversity. A growing number of studies have found that genetic (within-species) and species diversity can be correlated in space (the so-called species-gene diversity correlation [SGDC]), which suggests that they are controlled by nonindependent processes. Positive SGDCs are generally assumed to arise from parallel responses of genetic and species diversity to variation in site size and connectivity. However, this argument implicitly assumes a neutral model that has yet to be developed. Here, we build such a model to predict SGDC in a metacommunity. We describe how SGDC emerges from competition within sites and variation in connectivity and carrying capacity among sites. We then introduce the formerly ignored mutation process, which affects genetic but not species diversity. When mutation rate is low, our model confirms that variation in the number of migrants among sites creates positive SGDCs. However, when considering high mutation rates, interactions between mutation, migration, and competition can produce negative SGDCs. Neutral processes thus do not always contribute positively to SGDCs. Our approach provides empirical guidelines for interpreting these novel patterns in natura with respect to evolutionary and ecological forces shaping metacommunities.

Target-gene amplicon sequencing is the most exploited high-throughput sequencing application in microbial ecology. The targets are taxonomically relevant genes, with 16S rRNA being the gold standard for bacteria. As for fungi, the most commonly used target is the internal transcribed spacer (ITS). However, the uneven ITS length among species may promote preferential amplification and sequencing and incorrect estimation of their abundance. Therefore, the use of different targets is desirable. We evaluated the use of three different target amplicons for the characterization of fungal diversity. After an in silico primer evaluation, we compared three amplicons (the ITS1-ITS2 region [ITS1-2], 18S ribosomal small subunit RNA, and the D1/D2 domain of the 26S ribosomal large subunit RNA), using biological samples and a mock community of common fungal species. All three targets allowed for accurate identification of the species present. Nevertheless, high heterogeneity in ITS1-2 length was found, and this caused an overestimation of the abundance of species with a shorter ITS, while both 18S and 26S amplicons allowed for more reliable quantification. We demonstrated that ITS1-2 amplicon sequencing, although widely used, may lead to an incorrect evaluation of fungal communities, and efforts should be made to promote the use of different targets in sequencing-based microbial ecology studies.IMPORTANCE Amplicon-sequencing approaches for fungi may rely on different targets affecting the diversity and abundance of the fungal species. An increasing number of studies will address fungal diversity by high-throughput amplicon sequencing. The description of the communities must be accurate and reliable in order to draw useful insights and to address both ecological and biological questions. By analyzing a mock community and several biological samples, we demonstrate that using different amplicon targets may change the results of fungal microbiota analysis, and we highlight how a careful choice of the target is fundamental for a thorough description of the fungal communities.

Raw milk microbiota are complex communities with a significant impact on the hygienic, sensory and technological quality of milk products. However, there is a lack of knowledge on factors determining their composition. In the present study, four bulk tank milk samples of two farms at two different time points were analyzed in detail for their microbiota using cultivation and 16S rRNA amplicon sequencing. Diversity in samples from the first time point was assessed via cultivation of 500 aerobic mesophilic bacterial isolates in each sample. A high biodiversity of 70 and 110 species per sample was determined, of which 25–28% corresponded to yet unknown taxa. The isolates were dominated by Gram-positive members of the genera Staphylococcus, Corynebacterium, Streptococcus, or Janibacter, whilst Chryseobacterium and Acinetobacter were most abundant among the Gram-negative taxa. At the second time point, samples of the same farms were analyzed via both cultivation (1,500 individual colonies each) and high-throughput 16S rRNA gene amplicon sequencing. The latter revealed a threefold higher biodiversity at the genus level, as anaerobic or fastidious species were also detected. However, cultivation identified genera not captured by sequencing, indicating that both approaches are complementary. Using amplicon sequencing, the relative abundance of a few genera was distorted, which seems to be an artifact of sample preparation. Therefore, attention needs to be paid to the library preparation procedure with special emphasis on cell lysis and PCR.

Biodiversity is the sum total of all living things on the earth with particular reference to the profound variety in structure,function and genetic constitution. It includes both number and frequency of species or genes in a given assemblage and the variety of resulting ecosystems in a region. It is usually considered at three different levels: genetic, species and ecological diversities. Genus Drosophila belongs to the family Drosophilidae (class Insecta, order Diptera), characterized by rich species diversity at global level and also in India, which is a megadiverse country. At global level, more than 1500 species have been described and several thousands estimated. Hawaiian Islands are particularly rich in species diversity with more than 500 species which provides a unique opportunity to study evolution in genus Drosophila. About 150 species of Drosophila have been reported from India. Certain species of Drosophila found in India have been investigated for genetic diversity within the species. In this regard, Drosophila ananassae is noteworthy. It is a cosmopolitan and domestic species with common occurrence in India and is endowed with many genetic peculiarities. Population genetics and evolutionary studies in this species have revealed as to how genetic diversity within a species play an important role in adaptation of populations to varying environments. In addition, the work carried on D. melanogaster, D. nasuta, D. bipectinata and certain other species in India has shown that these species vary in degree and pattern of genetic diversity, and have evolved different mechanisms for adjusting to their environments. The ecological adaptations to various kinds of stress studied in certain species of Drosophila inhabiting the Indian subcontinent are also discussed.

BackgroundSpecies- and genetic diversity can change in parallel, resulting in a species-genetic diversity correlation (SGDC) and raising the question if the same drivers influence both biological levels of diversity. The SGDC can be either positive or negative, depending on whether the species diversity and the genetic diversity of the measured species respond in the same or opposite way to drivers. Using a traditional species diversity approach together with ultra-conserved elements and high throughput sequencing, we evaluated the SGDCs in benthic macrofauna communities in the Baltic Sea, a geologically young brackish water sea characterised by its steep salinity gradient and low species richness. Assessing SGDCs from six focal marine invertebrate species from different taxonomic groups and with differing life histories and ecological functions on both a spatial and temporal scale gives a more comprehensive insight into the community dynamics of this young ecosystem and the extrinsic factors that might drive the SGDCs.ResultsNo significant correlations between species diversity and genetic diversity were found for any of the focal species. However, both negative and positive trends of SGDCs for the individual focal species were observed. When examining the environmental drivers, no common trends between the species were found, even when restricting the analysis to specific taxonomic classes. Additionally, there were no common environmental factors driving the diversity relationships for species sharing the same SGDC trend (positive or negative). Local population dynamics, together with the invasion history of the individual species and their unique adaptation to the distinctive environment of the Baltic Sea, are expected to be of major influence on the outcome of the SGDCs.ConclusionsThe present results highlight the importance of assessing SGDCs using multiple species, not just a single indicator species. This emphasises a need to pay attention to the ecology and life history of the focal species. This study also provides insight into the large differences in both patterns and drivers of genetic diversity, which is important when including genetic biodiversity in conservation plans. We conclude that the effects of environmental and biological factors and processes that affects diversity patterns at both the community and genetic levels are likely species dependent, even in an environment such as the Baltic Sea with strong environmental gradients.

Analysis of the correlation between genetic and species diversity in a temperate forest: Variation in cohorts and effect of disturbance

This monograph presents results of research on genetic diversity of 8 leafy liverwort species differing in reproductive mode. The frst 4 species in Poland are regarded as sterile and reproduce only vegetatively: Bazzania trilobata, Trichoc-olea tomentella, Lophozia hatcheri, and Mylia anomala. The next 4 are fertile, including the monoecious Lepidozia reptans and Calypogeia integristipula as well as the dioecious Mylia taylorii and Tritomaria quinquedentata. For each species, 9-10 populations were sampled. In total, 4744 gametophytes from 73 populations were examined by isozyme analysis. The level of their genetic diversity (total, HT, and within populations, HS) was high, higher than in thallose liverworts, but comparable to the genetic diversity of mosses or even some species of vascular plants. Thus the traditional opinion that the entire group of liverworts has a much lower level of genetic diversity than mosses is erroneous, as it holds true only for thallose liverworts (Metzgeriidae and Marchantiopsida). My results indicate that the effect of reproductive mode on genetic diversity in leafy liverworts is lower than in vascular plants. Sterile and fertile species of liverworts exhibited similar levels of genetic diversity. Moreover, both groups included species that had both high and low levels of HT and HS. In fertile species, monoecious and dioecious species also did not differ signifcantly in genetic diversity, but dioecious liverworts had slightly higher total diversity (HT) than monoecious species. In most of the studied leafy liverworts, the share of genetic diversity within populations in the total genetic diversity of species is greater than between populations. The percentage share of variation among populations (ΦPT) in the total genetic variation was correlated with the total genetic diversity of the species (HT). In species with high HT, differences between populations tended to be rather small. By contrast, in species with lower HT, the percentage share of differentiation among populations in the total diversity of species was much higher. My results confrm theory, based on studies by Kimura, that the main causes of genetic diversity of bryophytes are neutral somatic mutations developing in various vegetative parts of plants. The separation of branches or other plant sections with somatic mutations, followed by the growth of new shoots, can increase the level of genetic diversity. The high level of genetic diversity in sterile liverworts indicates that vegetative reproduction has a greater infuence on the level of genetic diversity than recombination. My results suggest also that mutation rates are similar in closely related species, but species with a wider ecological range exhibit higher genetic diversity because the variability of habitats can infuence the rate and type of somatic mutations. Accordingly, species inhabiting more diverse environments may be more genetically diverse. Patches of the studied species generally consisted of several genotypes (MLGs). Two types of distribution of genotypes in patches were noticed. Patches of species with low total diversity (HT), were often dominated by 1-2 genotypes, which constituted the major part of a patch. In patches of species with higher HT, there was no tendency to form patches with predomination of a single genotype. Different genotypes constituted similar proportions of a patch. In all the studied leafy liverwort species there was a high degree of repeatability of the same genotypes (MLGs) in plants from various patches within the same population or in various populations. Probably the main cause of this is the independent repeatability of the same mutations in different specimens.