Assessing the genetic composition of cotton-top tamarins (Saguinus oedipus) before sweeping anthropogenic impact.

During the last century, the critically endangered cotton-top tamarin (Saguinus oedipus) has been threatened by multiple anthropogenic factors that drastically affected their habitat and population size. As the genetic impact of these pressures is largely unknown, this study aimed to establish a genetic baseline with the use of temporal sampling to determine the genetic makeup before detrimental anthropogenic impact. Genomes were resequenced from a combination of historical museum samples and modern wild samples at low-medium coverage, to unravel how the cotton-top tamarin population structure and genomic diversity may have changed during this period. Our data suggest two populations can be differentiated, probably separated historically by the mountain ranges of the Paramillo Massif in Colombia. Although this population structure persists in the current populations, modern samples exhibit genomic signals consistent with recent inbreeding, such as long runs of homozygosity and a reduction in genome-wide heterozygosity especially in the greater northeast population. This loss is likely the consequence of the population reduction following the mass exportation of cotton-top tamarins for biomedical research in the 1960s, coupled with the habitat loss this species continues to experience. However, current populations have not experienced an increase in genetic load. We propose that the historical genetic baseline established in this study can be used to provide insight into alteration in the modern population influenced by a drastic reduction in population size as well as providing background information to be used for future conservation decision-making for the species.

A simplified bacterial community found within the epidermis than at the epidermal surface of atopic dermatitis patients and healthy controls.

There has been considerable research into the understanding of the healthy skin microbiome. Similarly, there is also a considerable body of research into whether specific microbes contribute to skin disorders, with atopic dermatitis (AD) routinely linked to increased Staphylococcus aureus (S. aureus) colonisation. In this study, the epidermal surface of participants was sampled using swabs, while serial tape-stripping (35 tapes) was performed to sample through the stratum corneum. Samples were taken from AD patients and healthy controls, and the bacterial communities were profiled by metabarcoding the universal V3-V4 16S rRNA region. Results show that the majority of bacterial richness is located within the outermost layers of the stratum corneum, however there were many taxa that were found almost exclusively at the very outermost layer of the epidermis. We therefore hypothesise that tape-stripping can be performed to investigate the 'core microbiome' of participants by removing environmental contaminants. Interestingly, significant community variation between AD patients and healthy controls was only observable at the epidermal surface, yet a number of individual taxa were found to consistently differ with AD status across the entire epidermis (i.e. both the epidermal surface and within the epidermis). Sampling strategy could therefore be tailored dependent on the hypothesis, with sampling for forensic applications best performed using surface swabs and outer tapes, while profiling sub-surface communities may better reflect host genome and immunological status.

Temporal and Spatial Variation of the Skin-Associated Bacteria from Healthy Participants and Atopic Dermatitis Patients.

Several factors have been shown to influence the composition of the bacterial communities inhabiting healthy skin, with variation between different individuals, differing skin depths, and body locations (spatial-temporal variation). Atopic dermatitis (AD) is a chronic skin disease also affecting the skin-associated bacterial communities. While the effects of AD have been studied on these processes individually, few have considered how AD disrupts the spatial-temporal variation of the skin bacteria as a whole (i.e., considered these processes simultaneously). Here, we characterized the skin-associated bacterial communities of healthy volunteers and lesional and nonlesional skin of AD patients by metabarcoding the universal V3-V4 16S rRNA region from tape strip skin samples. We quantified the spatial-temporal variation (interindividual variation, differing skin depths, multiple time points) of the skin-associated bacteria within healthy controls and AD patients, including the relative change induced by AD in each. Interindividual variation correlated with the bacterial community far more strongly than any other factors followed by skin depth and then AD status. There was no significant temporal variation found within either AD patients or healthy controls. The bacterial community was found to vary markedly according to AD severity, and between patients without and with filaggrin mutations. Therefore, future studies may benefit from sampling subsurface epidermal communities and considering AD severity and the host genome in understanding the role of the skin bacterial community within AD pathogenesis rather than considering AD as a presence-absence disorder. IMPORTANCE The bacteria associated with human skin may influence skin barrier function and the immune response. Previous studies have attempted to understand the factors that regulate the skin bacteria, characterizing the spatial-temporal variation of the skin bacteria within unaffected skin. Here, we quantified the effect of AD on the skin bacteria on multiple spatial-temporal factors simultaneously. Although significant community variation between healthy controls and AD patients was observed, the effects of AD on the overall bacterial community were relatively low compared to other measured factors. Results here suggest that changes in specific taxa rather than wholesale changes in the skin bacteria are associated with mild to moderate AD. Further studies would benefit from incorporating the complexity of AD into models to better understand the condition, including AD severity and the host genome, alongside microbial composition.

Detecting flying insects using car nets and DNA metabarcoding.

Monitoring insects across space and time is challenging, due to their vast taxonomic and functional diversity. This study demonstrates how nets mounted on rooftops of cars (car nets) and DNA metabarcoding can be applied to sample flying insect richness and diversity across large spatial scales within a limited time period. During June 2018, 365 car net samples were collected by 151 volunteers during two daily time intervals on 218 routes in Denmark. Insect bulk samples were processed with a DNA metabarcoding protocol to estimate taxonomic composition, and the results were compared to known flying insect richness and occurrence data. Insect and hoverfly richness and diversity were assessed across biogeographic regions and dominant land cover types. We detected 15 out of 19 flying insect orders present in Denmark, with high proportions of especially Diptera compared to Danish estimates, and lower insect richness and diversity in urbanized areas. We detected 319 species not known for Denmark and 174 species assessed in the Danish Red List. Our results indicate that the methodology can assess the flying insect fauna at large spatial scales to a wide extent, but may be, like other methods, biased towards certain insect orders.

Diverse variola virus (smallpox) strains were widespread in northern Europe in the Viking Age.

Smallpox, one of the most devastating human diseases, killed between 300 million and 500 million people in the 20th century alone. We recovered viral sequences from 13 northern European individuals, including 11 dated to ~600-1050 CE, overlapping the Viking Age, and reconstructed near-complete variola virus genomes for four of them. The samples predate the earliest confirmed smallpox cases by ~1000 years, and the sequences reveal a now-extinct sister clade of the modern variola viruses that were in circulation before the eradication of smallpox. We date the most recent common ancestor of variola virus to ~1700 years ago. Distinct patterns of gene inactivation in the four near-complete sequences show that different evolutionary paths of genotypic host adaptation resulted in variola viruses that circulated widely among humans.

Universal Dermal Microbiome in Human Skin.

Human skin microbiota has been described as a "microbial fingerprint" due to observed differences between individuals. Current understanding of the cutaneous microbiota is based on sampling the outermost layers of the epidermis, while the microbiota in the remaining skin layers has not yet been fully characterized. Environmental conditions can vary drastically between the cutaneous compartments and give rise to unique communities. We demonstrate that the dermal microbiota is surprisingly similar among individuals and contains a specific subset of the epidermal microbiota. Variability in bacterial community composition decreased significantly from the epidermal to the dermal compartment but was similar among anatomic locations (hip and knee). The composition of the epidermal microbiota was more strongly affected by environmental factors than that of the dermal community. These results indicate a well-conserved dermal community that is functionally distinct from the epidermal community, challenging the current dogma. Future studies in cutaneous disorders and chronic infections may benefit by focusing on the dermal microbiota as a persistent microbial community.IMPORTANCE Human skin microbiota is thought to be unique according to the individual's lifestyle and genetic predisposition. This is true for the epidermal microbiota, while our findings demonstrate that the dermal microbiota is universal between healthy individuals. The preserved dermal microbial community is compositionally unique and functionally distinct to the specific environment in the depth of human skin. It is expected to have direct contact with the immune response of the human host, and research in the communication between host and microbiota should be targeted to this cutaneous compartment. This novel insight into specific microbial adaptation can be used advantageously in the research of chronic disorders and infections of the skin. It can enlighten the alteration between health and disease to the benefit of patients suffering from long-lasting socioeconomic illnesses.

A systematic survey of regional multi-taxon biodiversity: evaluating strategies and coverage.

BACKGROUND: In light of the biodiversity crisis and our limited ability to explain variation in biodiversity, tools to quantify spatial and temporal variation in biodiversity and its underlying drivers are critically needed. Inspired by the recently published ecospace framework, we developed and tested a sampling design for environmental and biotic mapping. We selected 130 study sites (40 × 40 m) across Denmark using stratified random sampling along the major environmental gradients underlying biotic variation. Using standardized methods, we collected site species data on vascular plants, bryophytes, macrofungi, lichens, gastropods and arthropods. To evaluate sampling efficiency, we calculated regional coverage (relative to the known species number per taxonomic group), and site scale coverage (i.e., sample completeness per taxonomic group at each site). To extend taxonomic coverage to organisms that are difficult to sample by classical inventories (e.g., nematodes and non-fruiting fungi), we collected soil for metabarcoding. Finally, to assess site conditions, we mapped abiotic conditions, biotic resources and habitat continuity. RESULTS: Despite the 130 study sites only covering a minute fraction (0.0005%) of the total Danish terrestrial area, we found 1774 species of macrofungi (54% of the Danish fungal species pool), 663 vascular plant species (42%), 254 bryophyte species (41%) and 200 lichen species (19%). For arthropods, we observed 330 spider species (58%), 123 carabid beetle species (37%) and 99 hoverfly species (33%). Overall, sample coverage was remarkably high across taxonomic groups and sufficient to capture substantial spatial variation in biodiversity across Denmark. This inventory is nationally unprecedented in detail and resulted in the discovery of 143 species with no previous record for Denmark. Comparison between plant OTUs detected in soil DNA and observed plant species confirmed the usefulness of carefully curated environmental DNA-data. Correlations among species richness for taxonomic groups were predominantly positive, but did not correlate well among all taxa suggesting differential and complex biotic responses to environmental variation. CONCLUSIONS: We successfully and adequately sampled a wide range of diverse taxa along key environmental gradients across Denmark using an approach that includes multi-taxon biodiversity assessment and ecospace mapping. Our approach is applicable to assessments of biodiversity in other regions and biomes where species are structured along environmental gradient.

Predicting provenance of forensic soil samples: Linking soil to ecological habitats by metabarcoding and supervised classification.

Environmental DNA (eDNA) is increasingly applied in ecological studies, including studies with the primary purpose of criminal investigation, in which eDNA from soil can be used to pair samples or reveal sample provenance. We collected soil eDNA samples as part of a large national biodiversity research project across 130 sites in Denmark. We investigated the potential for soil eDNA metabarcoding in predicting provenance in terms of environmental conditions, habitat type and geographic regions. We used linear regression for predicting environmental gradients of light, soil moisture, pH and nutrient status (represented by Ellenberg Indicator Values, EIVs) and Quadratic Discriminant Analysis (QDA) to predict habitat type and geographic region. eDNA data performed relatively well as a predictor of environmental gradients (R2 > 0.81). Its ability to discriminate between habitat types was variable, with high accuracy for certain forest types and low accuracy for heathland, which was poorly predicted. Geographic region was also less accurately predicted by eDNA. We demonstrated the application of provenance prediction in forensic science by evaluating and discussing two mock crime scenes. Here, we listed the plant species from annotated sequences, which can further aid in identifying the likely habitat or, in case of rare species, a geographic region. Predictions of environmental gradients and habitat types together give an overall accurate description of a crime scene, but care should be taken when interpreting annotated sequences, e.g. due to erroneous assignments in GenBank. Our approach demonstrates that important habitat properties can be derived from soil eDNA, and exemplifies a range of potential applications of eDNA in forensic ecology.

High-Throughput Sequencing-Based Investigation of Viruses in Human Cancers by Multienrichment Approach.

BACKGROUND: Viruses and other infectious agents cause more than 15% of human cancer cases. High-throughput sequencing-based studies of virus-cancer associations have mainly focused on cancer transcriptome data. METHODS: In this study, we applied a diverse selection of presequencing enrichment methods targeting all major viral groups, to characterize the viruses present in 197 samples from 18 sample types of cancerous origin. Using high-throughput sequencing, we generated 710 datasets constituting 57 billion sequencing reads. RESULTS: Detailed in silico investigation of the viral content, including exclusion of viral artefacts, from de novo assembled contigs and individual sequencing reads yielded a map of the viruses detected. Our data reveal a virome dominated by papillomaviruses, anelloviruses, herpesviruses, and parvoviruses. More than half of the included samples contained 1 or more viruses; however, no link between specific viruses and cancer types were found. CONCLUSIONS: Our study sheds light on viral presence in cancers and provides highly relevant virome data for future reference.

High-throughput sequencing reveals no viral pathogens in eight cases of ocular adnexal extranodal marginal zone B-cell lymphoma.

The purpose of the present study was to perform a next-generation sequencing (NGS) based analysis of viruses in ocular adnexal extranodal marginal zone B-cell lymphoma (EMZL). Eight patients with extraocular EMZL were identified in the archives of Department of Ophthalmology, Copenhagen University Hospital. All cases were validated according to the World Health Organization classification. We subjected samples to enrichment of virion-associated (encapsidated) nucleic acids which included sample homogenization, filtration, and nuclease treatment. Both DNA and RNA were sequenced, and we analyzed the sequencing data for the presence of viral sequences. We detected no pathogenic viruses likely to be associated to development of EMZL. In one case, we detected human polyomavirus 7 and traces of Epstein-Barr virus (EBV) (human herpesvirus 4 (HHV4)) and a human papillomavirus. In conclusion, no viral pathogens were consistently detected in the extraocular EMZL samples when applying NGS-based methods.

Metagenomic analysis of viruses in toilet waste from long distance flights-A new procedure for global infectious disease surveillance.

Human viral pathogens are a major public health threat. Reliable information that accurately describes and characterizes the global occurrence and transmission of human viruses is essential to support national and global priority setting, public health actions, and treatment decisions. However, large areas of the globe are currently without surveillance due to limited health care infrastructure and lack of international cooperation. We propose a novel surveillance strategy, using metagenomic analysis of toilet material from international air flights as a method for worldwide viral disease surveillance. The aim of this study was to design, implement, and evaluate a method for viral analysis of airplane toilet waste enabling simultaneous detection and quantification of a wide range of human viral pathogens. Toilet waste from 19 international airplanes was analyzed for viral content, using viral capture probes followed by high-throughput sequencing. Numerous human pathogens were detected including enteric and respiratory viruses. Several geographic trends were observed with samples originating from South Asia having significantly higher viral species richness as well as higher abundances of salivirus A, aichivirus A and enterovirus B, compared to samples originating from North Asia and North America. In addition, certain city specific trends were observed, including high numbers of rotaviruses in airplanes departing from Islamabad. Based on this study we believe that central sampling and analysis at international airports could be a useful supplement for global viral surveillance, valuable for outbreak detection and for guiding public health resources.

MobiSeq: De novo SNP discovery in model and non-model species through sequencing the flanking region of transposable elements.

In recent years, the availability of reduced representation library (RRL) methods has catalysed an expansion of genome-scale studies to characterize both model and non-model organisms. Most of these methods rely on the use of restriction enzymes to obtain DNA sequences at a genome-wide level. These approaches have been widely used to sequence thousands of markers across individuals for many organisms at a reasonable cost, revolutionizing the field of population genomics. However, there are still some limitations associated with these methods, in particular the high molecular weight DNA required as starting material, the reduced number of common loci among investigated samples, and the short length of the sequenced site-associated DNA. Here, we present MobiSeq, a RRL protocol exploiting simple laboratory techniques, that generates genomic data based on PCR targeted enrichment of transposable elements and the sequencing of the associated flanking region. We validate its performance across 103 DNA extracts derived from three mammalian species: grey wolf (Canis lupus), red deer complex (Cervus sp.) and brown rat (Rattus norvegicus). MobiSeq enables the sequencing of hundreds of thousands loci across the genome and performs SNP discovery with relatively low rates of clonality. Given the ease and flexibility of MobiSeq protocol, the method has the potential to be implemented for marker discovery and population genomics across a wide range of organisms-enabling the exploration of diverse evolutionary and conservation questions.

Parasitic infections and resource economy of Danish Iron Age settlement through ancient DNA sequencing.

In this study, we screen archaeological soil samples by microscopy and analyse the samples by next generation sequencing to obtain results with parasites at species level and untargeted findings of plant and animal DNA. Three separate sediment layers of an ancient man-made pond in Hoby, Denmark, ranging from 100 BC to 200 AD, were analysed by microscopy for presence of intestinal worm eggs and DNA analysis were performed to identify intestinal worms and dietary components. Ancient DNA of parasites, domestic animals and edible plants revealed a change in use of the pond over time reflecting the household practice in the adjacent Iron Age settlement. The most abundant parasite found belonged to the Ascaris genus, which was not possible to type at species level. For all sediment layers the presence of eggs of the human whipworm Trichuris trichiura and the beef tapeworm Taenia saginata suggests continuous disposal of human faeces in the pond. Moreover, the continuous findings of T. saginata further imply beef consumption and may suggest that cattle were living in the immediate surrounding of the site throughout the period. Findings of additional host-specific parasites suggest fluctuating presence of other domestic animals over time: Trichuris suis (pig), Parascaris univalens (horse), Taenia hydatigena (dog and sheep). Likewise, alternating occurrence of aDNA of edible plants may suggest changes in agricultural practices. Moreover, the composition of aDNA of parasites, plants and vertebrates suggests a significant change in the use of the ancient pond over a period of three centuries.

Circular DNA elements of chromosomal origin are common in healthy human somatic tissue.

The human genome is generally organized into stable chromosomes, and only tumor cells are known to accumulate kilobase (kb)-sized extrachromosomal circular DNA elements (eccDNAs). However, it must be expected that kb eccDNAs exist in normal cells as a result of mutations. Here, we purify and sequence eccDNAs from muscle and blood samples from 16 healthy men, detecting ~100,000 unique eccDNA types from 16 million nuclei. Half of these structures carry genes or gene fragments and the majority are smaller than 25 kb. Transcription from eccDNAs suggests that eccDNAs reside in nuclei and recurrence of certain eccDNAs in several individuals implies DNA circularization hotspots. Gene-rich chromosomes contribute to more eccDNAs per megabase and the most transcribed protein-coding gene in muscle, TTN (titin), provides the most eccDNAs per gene. Thus, somatic genomes are rich in chromosome-derived eccDNAs that may influence phenotypes through altered gene copy numbers and transcription of full-length or truncated genes.

Correction to: Cancer associated fibroblasts (CAFs) are activated in cutaneous basal cell carcinoma and in the peritumoural skin.

After publication of the original article [1] it was identified that order of the author list had been presented incorrectly. The author Robert Gniadecki's surname was also incorrect in the original article.

Algorithm for post-clustering curation of DNA amplicon data yields reliable biodiversity estimates.

DNA metabarcoding is promising for cost-effective biodiversity monitoring, but reliable diversity estimates are difficult to achieve and validate. Here we present and validate a method, called LULU, for removing erroneous molecular operational taxonomic units (OTUs) from community data derived by high-throughput sequencing of amplified marker genes. LULU identifies errors by combining sequence similarity and co-occurrence patterns. To validate the LULU method, we use a unique data set of high quality survey data of vascular plants paired with plant ITS2 metabarcoding data of DNA extracted from soil from 130 sites in Denmark spanning major environmental gradients. OTU tables are produced with several different OTU definition algorithms and subsequently curated with LULU, and validated against field survey data. LULU curation consistently improves α-diversity estimates and other biodiversity metrics, and does not require a sequence reference database; thus, it represents a promising method for reliable biodiversity estimation.

Cancer associated fibroblasts (CAFs) are activated in cutaneous basal cell carcinoma and in the peritumoural skin.

BACKGROUND: Cutaneous basal cell carcinoma (BCC) is the commonest cancer worldwide. BCC is locally invasive and the surrounding stromal microenvironment is pivotal for tumourigenesis. Cancer associated fibroblasts (CAFs) in the microenvironment are essential for tumour growth in a variety of neoplasms but their role in BCC is poorly understood. METHODS: Material included facial BCC and control skin from the peritumoural area and from the buttocks. With next-generation sequencing (NGS) we compared mRNA expression between BCC and peritumoural skin. qRT-PCR, immunohistochemical and immunofluorescent staining were performed to validate the NGS results and to investigate CAF-related cyto-and chemokines. RESULTS: NGS revealed upregulation of 65 genes in BCC coding for extracellular matrix components pointing at CAF-related matrix remodeling. qRT-PCR showed increased mRNA expression of CAF markers FAP-α, PDGFR-ß and prolyl-4-hydroxylase in BCC. Peritumoural skin (but not buttock skin) also exhibited high expression of PDGFR-ß and prolyl-4-hydroxylase but not FAP-α. We found a similar pattern for the CAF-associated chemokines CCL17, CCL18, CCL22, CCL25, CXCL12 and IL6 with high expression in BCC and peritumoural skin but absence in buttock skin. Immunofluorescence revealed correlation between FAP-α and PDGFR-ß and CXCL12 and CCL17. CONCLUSION: Matrix remodeling is the most prominent molecular feature of BCC. CAFs are present within BCC stroma and associated with increased expression of chemokines involved in tumour progression and immunosuppression (CXCL12, CCL17). Fibroblasts from chronically sun-exposed skin near tumours show gene expression patterns resembling that of CAFs, indicating that stromal fibroblasts in cancer-free surgical BCC margins exhibit a tumour promoting phenotype.

Mitochondrial DNA of pre-last glacial maximum red deer from NW Spain suggests a more complex phylogeographical history for the species.

The major climatic oscillations that characterized the Quaternary had a great influence on the evolution and distribution of several species. During cold periods, the distribution of temperate-adapted species became fragmented with many surviving in southern refugia (Iberian, Italian, and Balkan Peninsulas). Red deer was one of the species that contracted its original range to southern refugia. Currently, two main lineages have been described for the species: western and eastern. We have analyzed fossils pre-dating the last glacial maximum (LGM) from Liñares cave (NW Spain) that belongs to the peripheral range of the western clade, and fossils from the Danish Holocene belonging to the central part of the same clade. Phylogenetic analyses place our samples in the western clade. However, some specimens from Liñares represent an early split in the tree along with other pre-LGM western samples from previous studies. Despite low bootstrap values in the Bayesian phylogenies, haplotype networks connect these foreign haplotypes to the eastern clade. We suggest a mixed phylogeographical model to explain this pattern with range expansions from the east during the expansion phase after the cold periods in marine isotope stage 3. We find slight isolation by distance in post-LGM populations that could be a consequence of the recolonization from southern refugia after the LGM.

Substitutions of short heterologous DNA segments of intragenomic or extragenomic origins produce clustered genomic polymorphisms.

In a screen for unexplained mutation events we identified a previously unrecognized mechanism generating clustered DNA polymorphisms such as microindels and cumulative SNPs. The mechanism, short-patch double illegitimate recombination (SPDIR), facilitates short single-stranded DNA molecules to invade and replace genomic DNA through two joint illegitimate recombination events. SPDIR is controlled by key components of the cellular genome maintenance machinery in the gram-negative bacterium Acinetobacter baylyi. The source DNA is primarily intragenomic but can also be acquired through horizontal gene transfer. The DNA replacements are nonreciprocal and locus independent. Bioinformatic approaches reveal occurrence of SPDIR events in the gram-positive human pathogen Streptococcus pneumoniae and in the human genome.