Comparative performance of the BGISEQ-500 vs Illumina HiSeq2500 sequencing platforms for palaeogenomic sequencing

Field and laboratory guidelines for reliable bioinformatic and statistical analysis of bacterial shotgun metagenomic data

Shotgun metagenomics is an increasingly cost-effective approach for profiling environmental and host-associated microbial communities. However, due to the complexity of both microbiomes and the molecular techniques required to analyze them, the reliability and representativeness of the results are contingent upon the field, laboratory, and bioinformatic procedures employed. Here, we consider 15 field and laboratory issues that critically impact downstream bioinformatic and statistical data processing, as well as result interpretation, in bacterial shotgun metagenomic studies. The issues we consider encompass intrinsic properties of samples, study design, and laboratory-processing strategies. We identify the links of field and laboratory steps with downstream analytical procedures, explain the means for detecting potential pitfalls, and propose mitigation measures to overcome or minimize their impact in metagenomic studies. We anticipate that our guidelines will assist data scientists in appropriately processing and interpreting their data, while aiding field and laboratory researchers to implement strategies for improving the quality of the generated results.

Host-gut microbiota interactions shape parasite infections in farmed Atlantic salmon

Animals and their associated microbiota share long evolutionary histories. However, it is not always clear how host genotype and microbiota interact to affect phenotype. We applied a hologenomic approach to explore how host-microbiota interactions shape lifetime growth and parasite infection in farmed Atlantic salmon (Salmo salar). Multi-omics data sets were generated from the guts of 460 salmon, 82% of which were naturally infected with an intestinal cestode. A single Mycoplasma bacterial strain, MAG01, dominated the gut metagenome of large, non-parasitized fish, consistent with previous studies showing high levels of Mycoplasma in the gut microbiota of healthy salmon. While small and/or parasitized salmon also had high abundance of MAG01, we observed increased alpha diversity in these individuals, driven by increased frequency of low-abundance Vibrionaceae and other Mycoplasma species that carried known virulence genes. Colonization by one of these cestode-associated Mycoplasma strains was associated with host individual genomic variation in long non-coding RNAs. Integrating the multi-omic data sets revealed coordinated changes in the salmon gut mRNA transcriptome and metabolome that correlated with shifts in the microbiota of smaller, parasitized fish. Our results suggest that the gut microbiota of small and/or parasitized fish is in a state of dysbiosis that partly depends on the host genotype, highlighting the value of using a hologenomic approach to incorporate the microbiota into the study of host-parasite dynamics.IMPORTANCEStudying host-microbiota interactions through the perspective of the hologenome is gaining interest across all life sciences. Intestinal parasite infections are a huge burden on human and animal health; however, there are few studies investigating the role of the hologenome during parasite infections. We address this gap in the largest multi-omics fish microbiota study to date using natural cestode infection of farmed Atlantic salmon. We find a clear association between cestode infection, salmon lifetime growth, and perturbation of the salmon gut microbiota. Furthermore, we provide the first evidence that the genetic background of the host may partly determine how the gut microbiota changes during parasite-associated dysbiosis. Our study therefore highlights the value of a hologenomic approach for gaining a more in-depth understanding of parasitism.

Ancient mitochondrial genomes unveil the origins and evolutionary history of New Zealand's enigmatic takahē and moho

Many avian species endemic to Aotearoa New Zealand were driven to extinction or reduced to relict populations following successive waves of human arrival, due to hunting, habitat destruction and the introduction of mammalian predators. Among the affected species were the large flightless South Island takahē (Porphyrio hochstetteri) and the moho (North Island takahē; P. mantelli), with the latter rendered extinct and the former reduced to a single relictual population. Little is known about the evolutionary history of these species prior to their decline and/or extinction. Here we sequenced mitochondrial genomes from takahē and moho subfossils (12 takahē and 4 moho) and retrieved comparable sequence data from takahē museum skins (n = 5) and contemporary individuals (n = 17) to examine the phylogeny and recent evolutionary history of these species. Our analyses suggest that prehistoric takahē populations lacked deep phylogeographic structure, in contrast to moho, which exhibited significant spatial genetic structure, albeit based on limited sample sizes (n = 4). Temporal genetic comparisons show that takahē have lost much of their mitochondrial genetic diversity, likely due to a sudden demographic decline soon after human arrival (~750 years ago). Time-calibrated phylogenetic analyses strongly support a sister species relationship between takahē and moho, suggesting these flightless taxa diverged around 1.5 million years ago, following a single colonisation of New Zealand by a flighted Porphyrio ancestor approximately 4 million years ago. This study highlights the utility of palaeogenetic approaches for informing the conservation and systematic understanding of endangered species whose ranges have been severely restricted by anthropogenic impacts.

COL3A1 is a potential diagnostic biomarker for synovial chondromatosis and affects the cell cycle and migration of chondrocytes

BACKGROUND

Synovial chondromatosis (SC) primarily affects the major joints and is characterized by the formation of benign cartilaginous nodules. In the present study, we evaluated the differences in the histology and gene expression of SC and normal cartilages and further elucidated the function of hub genes in SC.

METHODS

Histological staining and biochemical analysis were performed to measure collagen and glycosaminoglycan (GAG) contents in SC and normal cartilage samples. Then, microarray analysis was performed using knee joint samples (three normal and three SC samples) to identify the differentially expressed genes (DEGs). Subsequently, bioinformatics analysis was performed to identify the hub genes and explore the mechanisms underlying SC. The intersection of the top 10 upregulated DEGs, top 10 downregulated DEGs, and hub genes was validated in SC tissues. Lastly, in vitro experiments and our clinical cohort were used to determine the potential biological functions and diagnostic value, respectively, of the most significant gene.

RESULTS

The GAG and collagen contents were comparable to or higher in SC tissues than in normal tissues. Microarray analysis revealed 143 upregulated and 107 downregulated DEGs in SC. Furthermore, functional enrichment analysis revealed an association between immunity and metabolism-related pathways and SC development. Among 20 hub genes, two intersection genes, namely, collagen type III alpha 1 chain (COL3A1) and HSPA8, were notably expressed in SC tissues, with COL3A1 exhibiting a more significant difference in mRNA expression. Furthermore, COL3A1 can promote chondrocyte migration and cell cycle progression. Additionally, clinical data revealed COL3A1 can be a diagnostic marker for primary SC (AUC = 0.82) and be a positive correlation with neutrophil-to-lymphocyte ratio.

CONCLUSIONS

These results suggest that SC tissues contained the abundant GAG and collagen. COL3A1 can affect the function of chondrocytes and be a diagnostic marker of primary SC patients. These findings provide a novel approach and a fundamental contribution for diagnosis and treatment in SC.

Comparison of the DNBSEQ platform and Illumina HiSeq 2000 for bacterial genome assembly

The Illumina HiSeq platform has been a commonly used option for bacterial genome sequencing. Now the BGI DNA nanoball (DNB) nanoarrays platform may provide an alternative platform for sequencing of bacterial genomes. To explore the impact of sequencing platforms on bacterial genome assembly, quality assessment, sequence alignment, functional annotation, mutation detection, and metagenome mapping, we compared genome assemblies based on sequencing of cultured bacterial species using the HiSeq 2000 and BGISEQ-500 platforms. In addition, simulated reads were used to evaluate the impact of insert size on genome assembly. Genome assemblies based on BGISEQ-500 sequencing exhibited higher completeness and fewer N bases in high GC genomes, whereas HiSeq 2000 assemblies exhibited higher N50. The majority of assembly assessment parameters, sequences of 16S rRNA genes and genomes, numbers of single nucleotide variants (SNV), and mapping to metagenome data did not differ significantly between platforms. More insertions were detected in HiSeq 2000 genome assemblies, whereas more deletions were detected in BGISEQ-500 genome assemblies. Insert size had no significant impact on genome assembly. Taken together, our results suggest that DNBSEQ platforms would be a valid substitute for HiSeq 2000 for bacterial genome sequencing.

Comparative transcriptomic analysis of Illumina and MGI next-generation sequencing platforms using RUNX3- and ZBTB46-instructed embryonic stem cells

Introduction: We have previously observed phenotypic and developmental changes upon the ectopic expression of the RUNX3 or the ZBTB46 transcription factors in mouse embryonic stem cell (ESC) derived progenitors. In this study, we evaluated the gene expression profiles of the RUNX3- and the ZBTB46-instructed murine ESCs with RNA-seq testing two next-generation sequencing technologies. Methods: We compared the DNA nanoball-based DNBSEQ G400 sequencer (MGI) with the bridge-PCR-based NextSeq 500 instrument (Illumina) for RNA sequencing. Moreover, we also compared two types of MGI sequencing reagents (Standard versus Hot-massive parallel sequencing (MPS)) with the DNBSEQ G400. Results: We observed that both sequencing platforms showed comparable levels of quality, sequencing uniformity, and gene expression profiles. For example, highly overlapping RUNX3- and ZBTB46-regulated gene lists were obtained from both sequencing datasets. Moreover, we observed that the Standard and the Hot-MPS-derived RUNX3- and ZBTB46-regulated gene lists were also considerably overlapped. This transcriptome analysis also helped us to identify differently expressed genes in the presence of the transgenic RUNX3 or ZBTB46. For example, we found that Gzmb, Gzmd, Gzme, Gdf6, and Ccr7 genes were robustly upregulated upon the forced expression of Runx3; on the other hand, Gpx2, Tdpoz4, and Arg2 were induced alongside the ectopic expression of Zbtb46. Discussion: Similar gene expression profile and greatly overlapping RUNX3- and ZBTB46-regulated gene sets were detected with both DNA sequencing platforms. Our analyses demonstrate that both sequencing technologies are suitable for transcriptome profiling and target gene selection. These findings suggest that DNBSEQ G400 represents a cost-effective alternative sequencing platform for gene expression monitoring. Moreover, this analysis provides a resource for exploration of the RUNX3- and ZBTB46-dependent gene regulatory networks.

The history of Coast Salish "woolly dogs" revealed by ancient genomics and Indigenous Knowledge

Ancestral Coast Salish societies in the Pacific Northwest kept long-haired "woolly dogs" that were bred and cared for over millennia. However, the dog wool-weaving tradition declined during the 19th century, and the population was lost. In this study, we analyzed genomic and isotopic data from a preserved woolly dog pelt from "Mutton," collected in 1859. Mutton is the only known example of an Indigenous North American dog with dominant precolonial ancestry postdating the onset of settler colonialism. We identified candidate genetic variants potentially linked with their distinct woolly phenotype. We integrated these data with interviews from Coast Salish Elders, Knowledge Keepers, and weavers about shared traditional knowledge and memories surrounding woolly dogs, their importance within Coast Salish societies, and how colonial policies led directly to their disappearance.

Redefining the Evolutionary History of the Rock Dove, Columba livia, Using Whole Genome Sequences

The domestic pigeon's exceptional phenotypic diversity was key in developing Darwin's Theory of Evolution and establishing the concept of artificial selection. However, unlike its domestic counterpart, its wild progenitor, the rock dove Columba livia has received considerably less attention. Therefore, questions regarding its domestication, evolution, taxonomy, and conservation status remain unresolved. We generated whole-genome sequencing data from 65 historical rock doves that represent all currently recognized subspecies and span the species' original geographic distribution. Our dataset includes 3 specimens from Darwin's collection, and the type specimens of 5 different taxa. We characterized their population structure, genomic diversity, and gene-flow patterns. Our results show the West African subspecies C. l. gymnocyclus is basal to rock doves and domestic pigeons, and suggests gene-flow between the rock dove's sister species C. rupestris, and the ancestor of rock doves after its split from West African populations. These genomes allowed us to propose a model for the evolution of the rock dove in light of the refugia theory. We propose that rock dove genetic diversity and introgression patterns derive from a history of allopatric cycles and dispersion waves during the Quaternary glacial and interglacial periods. To explore the rock dove domestication history, we combined our new dataset with available genomes from domestic pigeons. Our results point to at least 1 domestication event in the Levant that gave rise to all domestic breeds analysed in this study. Finally, we propose a species-level taxonomic arrangement to reflect the evolutionary history of the West African rock dove populations.

BCL11B and the NuRD complex cooperatively guard T-cell fate and inhibit OPA1-mediated mitochondrial fusion in T cells

The nucleosome remodeling and histone deacetylase (NuRD) complex physically associates with BCL11B to regulate murine T-cell development. However, the function of NuRD complex in mature T cells remains unclear. Here, we characterize the fate and metabolism of human T cells in which key subunits of the NuRD complex or BCL11B are ablated. BCL11B and the NuRD complex bind to each other and repress natural killer (NK)-cell fate in T cells. In addition, T cells upregulate the NK cell-associated receptors and transcription factors, lyse NK-cell targets, and are reprogrammed into NK-like cells (ITNKs) upon deletion of MTA2, MBD2, CHD4, or BCL11B. ITNKs increase OPA1 expression and exhibit characteristically elongated mitochondria with augmented oxidative phosphorylation (OXPHOS) activity. OPA1-mediated elevated OXPHOS enhances cellular acetyl-CoA levels, thereby promoting the reprogramming efficiency and antitumor effects of ITNKs via regulating H3K27 acetylation at specific targets. In conclusion, our findings demonstrate that the NuRD complex and BCL11B cooperatively maintain T-cell fate directly by repressing NK cell-associated transcription and indirectly through a metabolic-epigenetic axis, providing strategies to improve the reprogramming efficiency and antitumor effects of ITNKs.

Historic Sampling of a Vanishing Beast: Population Structure and Diversity in the Black Rhinoceros

The black rhinoceros (Diceros bicornis L.) is a critically endangered species historically distributed across sub-Saharan Africa. Hunting and habitat disturbance have diminished both its numbers and distribution since the 19th century, but a poaching crisis in the late 20th century drove them to the brink of extinction. Genetic and genomic assessments can greatly increase our knowledge of the species and inform management strategies. However, when a species has been severely reduced, with the extirpation and artificial admixture of several populations, it is extremely challenging to obtain an accurate understanding of historic population structure and evolutionary history from extant samples. Therefore, we generated and analyzed whole genomes from 63 black rhinoceros museum specimens collected between 1775 and 1981. Results showed that the black rhinoceros could be genetically structured into six major historic populations (Central Africa, East Africa, Northwestern Africa, Northeastern Africa, Ruvuma, and Southern Africa) within which were nested four further subpopulations (Maasailand, southwestern, eastern rift, and northern rift), largely mirroring geography, with a punctuated north-south cline. However, we detected varying degrees of admixture among groups and found that several geographical barriers, most prominently the Zambezi River, drove population discontinuities. Genomic diversity was high in the middle of the range and decayed toward the periphery. This comprehensive historic portrait also allowed us to ascertain the ancestry of 20 resequenced genomes from extant populations. Lastly, using insights gained from this unique temporal data set, we suggest management strategies, some of which require urgent implementation, for the conservation of the remaining black rhinoceros diversity.

The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

Conservation implications of elucidating the Korean wolf taxonomic ambiguity through whole-genome sequencing

The taxonomic status of the now likely extirpated Korean Peninsula wolf has been extensively debated, with some arguing it represents an independent wolf lineage, Canis coreanus. To investigate the Korean wolf's genetic affiliations and taxonomic status, we sequenced and analysed the genomes of a Korean wolf dated to the beginning of the 20th century, and a captive wolf originally from the Pyongyang Central Zoo. Our results indicated that the Korean wolf bears similar genetic ancestry to other regional East Asian populations, therefore suggesting it is not a distinct taxonomic lineage. We identified regional patterns of wolf population structure and admixture in East Asia with potential conservation consequences in the Korean Peninsula and on a regional scale. We find that the Korean wolf has similar genomic diversity and inbreeding to other East Asian wolves. Finally, we show that, in contrast to the historical sample, the captive wolf is genetically more similar to wolves from the Tibetan Plateau; hence, Korean wolf conservation programmes might not benefit from the inclusion of this specimen.

In Silico Probiogenomic Characterization of Lactobacillus delbrueckii subsp. lactis A4 Strain Isolated from an Armenian Honeybee Gut

A Lactobacillus delbrueckii ssp. lactis strain named A4, isolated from the gut of an Armenian honeybee, was subjected to a probiogenomic characterization because of its unusual origin. A whole-genome sequencing was performed, and the bioinformatic analysis of its genome revealed a reduction in the genome size and the number of the genes-a process typical for the adaptation to endosymbiotic conditions. Further analysis of the genome revealed that Lactobacillus delbrueckii ssp. lactis strain named A4 could play the role of a probiotic endosymbiont because of the presence of intact genetic sequences determining antioxidant properties, exopolysaccharides synthesis, adhesion properties, and biofilm formation, as well as an antagonistic activity against some pathogens which is not due to pH or bacteriocins production. Additionally, the genomic analysis revealed significant potential for stress tolerance, such as extreme pH, osmotic stress, and high temperature. To our knowledge, this is the first report of a potentially endosymbiotic Lactobacillus delbrueckii ssp. lactis strain adapted to and playing beneficial roles for its host.

Identification of miRNA-mediated gene regulatory networks in L-methionine exposure counteracts cocaine-conditioned place preference in mice

Background and Aims: Methionine has been proven to inhibit addictive behaviors of cocaine dependence. This study aimed to identify the potential mechanisms of MET relating to its inhibitory effects on cocaine induced cellular and behavioral changes. Methods: MRNA and miRNA high-throughput sequencing of the prefrontal cortex in a mouse model of cocaine conditioned place preference (CPP) combined with L-methionine was performed. Differentially expressed miRNAs (DE-miRNAs) and differentially expressed genes (DEGs) regulated by cocaine and inhibited by L-methionine were identified. DEGs were mapped to STRING database to construct a protein-protein interaction (PPI) network. Then, the identified DEGs were subjected to the DAVID webserver for functional annotation. Finally, miRNA-mRNA regulatory network and miRNA-mRNA-TF regulatory networks were established to screen key DE-miRNAs and coregulation network in Cytoscape. Results: Sequencing data analysis showed that L-methionine reversely regulated genes and miRNAs affected by cocaine. Pathways associated with drug addiction only enriched in CS-down with MC-up genes targeted by DE-miRNAs including GABAergic synapse, Glutamatergic synapse, Circadian entrainment, Axon guidance and Calcium signaling pathway. Drug addiction associated network was formed of 22 DEGs including calcium channel (Cacna1c, Cacna1e, Cacna1g and Cacng8), ephrin receptor genes (Ephb6 and Epha8) and ryanodine receptor genes (Ryr1 and Ryr2). Calcium channel gene network were identified as a core gene network modulated by L-methionine in response to cocaine dependence. Moreover, it was predicted that Grin1 and Fosb presented in TF-miRNA-mRNA coregulation network with a high degree of interaction as hub genes and interacted calcium channels. Conclusion: These identified key genes, miRNA and coregulation network demonstrated the efficacy of L-methionine in counteracting the effects of cocaine CPP. To a certain degree, it may provide some hints to better understand the underlying mechanism on L-methionine in response to cocaine abuse.

From Genomics to Metagenomics in the Era of Recent Sequencing Technologies

Metagenomics, also known as environmental genomics, is the study of the genomic content of a sample of organisms obtained from a common habitat. Metagenomics and other "omics" disciplines have captured the attention of researchers for several decades. The effect of microbes in our body is a relevant concern for health studies. Through sampling the sequences of microbial genomes within a certain environment, metagenomics allows study of the functional metabolic capacity of a community as well as its structure based upon distribution and richness of species. Exponentially increasing number of microbiome literatures illustrate the importance of sequencing techniques which have allowed the expansion of microbial research into areas, including the human gut, antibiotics, enzymes, and more. This chapter illustrates how metagenomics field has evolved with the progress of sequencing technologies.Further, from this chapter, researchers will be able to learn about all current options for sequencing techniques and comparison of their cost and read statistics, which will be helpful for planning their own studies.

Loss of Mitochondrial Genetic Diversity despite Population Growth: The Legacy of Past Wolf Population Declines

Gray wolves (Canis lupus) in the Iberian Peninsula declined substantially in both range and population size in the last few centuries due to human persecution and habitat fragmentation. However, unlike many other western European populations, gray wolves never went extinct in Iberia. Since the minimum number was recorded around 1970, their numbers have significantly increased and then stabilized in recent decades. We analyzed mitochondrial genomes from 54 historical specimens of Iberian wolves from across their historical range using ancient DNA methods. We compared historical and current mitochondrial diversity in Iberian wolves at the 5' end of the control region (n = 17 and 27) and the whole mitochondrial genome excluding the control region (n = 19 and 29). Despite an increase in population size since the 1970s, genetic diversity declined. We identified 10 whole mitochondrial DNA haplotypes in 19 historical specimens, whereas only six of them were observed in 29 modern Iberian wolves. Moreover, a haplotype that was restricted to the southern part of the distribution has gone extinct. Our results illustrate a lag between demographic and genetic diversity changes, and show that after severe population declines, genetic diversity can continue to be lost in stable or even expanding populations. This suggests that such populations may be of conservation concern even after their demographic trajectory has been reversed.

The draft genome and multi-omics analyses reveal new insights into geo-herbalism properties of Citrus grandis 'Tomentosa'

Citrus grandis 'Tomentosa' (CGT) (Huajuhong, HJH) is a widely used medicinal plant, which is mainly produced in Guangdong and Guangxi provinces of South China. Particularly, HJH from Huazhou (HZ) county of Guangdong province has been well-regarded as the best national product for geo-herbalism. But the reasons for geo-herbalism property in HJH from HZ county remains a mystery. Therefore, a multi-omics approach was applied to identify the nature of the geo-herbalism in CGT from three different regions. The comprehensive screening of differential metabolites revealed that the Nobiletin content was significantly different in HZ region compared to other regions, and could be employed as a key indicator to determine the geo-herbalism. Furthermore, the high-quality genome (N50 of 9.12 Mb), coupled with genomics and transcriptomics analyses indicated that CGT and Citrus grandis are closely related, with a predicted divergence time of 19.1 million years ago (MYA), and no recent WGD occurred in the CGT, and the bioactive ingredients of CGT were more abundant than that of Citrus grandis. Interestingly, Nobiletin (Polymethoxyflavones) content was identified as a potential indicator of geo-herbalism, and O-methyltransferase (OMT) genes are involved in the synthesis of Polymethoxyflavones. Further multi-omics analysis led to the identification of a novel OMT gene (CtgOMT1) whose transient overexpression displayed significantly higher Nobiletin content, suggesting that CtgOMT1 was involved in the synthesis of Nobiletin. Overall, our findings provide new data resources for geo-herbalism evaluation, germplasm conservation and insights into Nobiletin biosynthesis pathways for the medicinal plant C. grandis 'Tomentosa'.

Detoxification Gene Families at the Genome-Wide Level of Rhus Gall Aphid Schlechtendalia chinensis

The Rhus gall aphid Schlechtendalia chinensis uses the species Rhus chinensis as its primary host plant, on which galls are produced. The galls have medicinal properties and can be used in various situations due to their high tannin content. Detoxification enzymes play significant roles in the insect lifecycle. In this study, we focused on five detoxification gene families, i.e., glutathione-S-transferase (GST), ABC transporter (ABC), Carboxylesterase (CCE), cyto-chrome P450 (CYP), and UDP-glycosyltransferase (UDP), and manually annotated 144 detoxification genes of S. chinensis using genome-wide techniques. The detoxification genes appeared mostly on chromosome 1, where a total of two pair genes were identified to show tandem duplications. There were 38 gene pairs between genomes of S. chinensis and Acyrthosiphon pisum in the detoxification gene families by collinear comparison. Ka/Ks ratios showed that detoxification genes of S. chinensis were mainly affected by purification selection during evolution. The gene expression numbers of P450s and ABCs by transcriptome sequencing data were greater, while gene expression of CCEs was the highest, suggesting they might be important in the detoxification process. Our study has firstly identified the genes of the different detoxification gene families in the S. chinensis genome, and then analyzed their general features and expression, demonstrating the importance of the detoxification genes in the aphid and providing new information for further research.

Benchmarking of ATAC Sequencing Data From BGI's Low-Cost DNBSEQ-G400 Instrument for Identification of Open and Occupied Chromatin Regions

Background: Chromatin falls into one of two major subtypes: closed heterochromatin and euchromatin which is accessible, transcriptionally active, and occupied by transcription factors (TFs). The most widely used approach to interrogate differences in the chromatin state landscape is the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). While library generation is relatively inexpensive, sequencing depth requirements can make this assay cost-prohibitive for some laboratories. Findings: Here, we benchmark data from Beijing Genomics Institute's (BGI) DNBSEQ-G400 low-cost sequencer against data from a standard Illumina instrument (HiSeqX10). For comparisons, the same bulk ATAC-seq libraries generated from pluripotent stem cells (PSCs) and fibroblasts were sequenced on both platforms. Both instruments generate sequencing reads with comparable mapping rates and genomic context. However, DNBSEQ-G400 data contained a significantly higher number of small, sub-nucleosomal reads (>30% increase) and a reduced number of bi-nucleosomal reads (>75% decrease), which resulted in narrower peak bases and improved peak calling, enabling the identification of 4% more differentially accessible regions between PSCs and fibroblasts. The ability to identify master TFs that underpin the PSC state relative to fibroblasts (via HOMER, HINT-ATAC, TOBIAS), namely, foot-printing capacity, were highly similar between data generated on both platforms. Integrative analysis with transcriptional data equally enabled direct recovery of three published 3-factor combinations that have been shown to induce pluripotency. Conclusion: Other than a small increase in peak calling sensitivity for DNBSEQ-G400 data (BGI), both platforms enable comparable levels of open chromatin identification for ATAC-seq library sequencing, yielding similar analytical outcomes, albeit at low-data generation costs in the case of the BGI instrument.

Grey wolf genomic history reveals a dual ancestry of dogs

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived^1–8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000–30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

DNA from ancient wolves spanning 100,000 years sheds light on wolves’ evolutionary history and the genomic origin of dogs.

Microbiome "Inception": an Intestinal Cestode Shapes a Hierarchy of Microbial Communities Nested within the Host

The concept of a holobiont, a host organism and its associated microbial communities, encapsulates the vital role the microbiome plays in the normal functioning of its host. Parasitic infections can disrupt this relationship, leading to dysbiosis. However, it is increasingly recognized that multicellular parasites are themselves holobionts. Intestinal parasites share space with the host gut microbiome, creating a system of nested microbiomes within the primary host. However, how the parasite, as a holobiont, interacts with the host holobiont remains unclear, as do the consequences of these interactions for host health. Here, we used 16S amplicon and shotgun metagenomics sequencing to characterize the microbiome of the intestinal cestode Eubothrium and its effect on the gut microbiome of its primary host, Atlantic salmon. Our results indicate that cestode infection is associated with salmon gut dysbiosis by acting as a selective force benefiting putative pathogens and potentially introducing novel bacterial species to the host. Our results suggest that parasitic cestodes may themselves be holobionts nested within the microbial community of their holobiont host, emphasizing the importance of also considering microbes associated with parasites when studying intestinal parasitic infections. IMPORTANCE The importance of the parasite microbiome is gaining recognition. Of particular concern is understanding how these parasite microbiomes influence host-parasite interactions and parasite interactions with the vertebrate host microbiome as part of a system of nested holobionts. However, there are still relatively few studies focusing on the microbiome of parasitic helminths in general and almost none on cestodes in particular, despite the significant burden of disease caused by these parasites globally. Our study provides insights into a system of significance to the aquaculture industry, cestode infections of Atlantic salmon and, more broadly, expands our general understanding of parasite-microbiome-host interactions and introduces a new element, the microbiome of the parasite itself, which may play a critical role in modulating the host microbiome, and, therefore, the host response, to parasite infection.

Genetic evidence for post-glacial expansion from a southern refugium in the eastern moa (Emeus crassus)

Cycles of glacial expansion and contraction throughout the Pleistocene drove increases and decreases, respectively, in the geographical range and population size of many animal species. Genetic data have revealed that during glacial maxima the distribution of many Eurasian animals was restricted to small refugial areas, from which species expanded to reoccupy parts of their former range as the climate warmed. It has been suggested that the extinct eastern moa (Emeus crassus)-a large, flightless bird from New Zealand-behaved analogously during glacial maxima, possibly surviving only in a restricted area of lowland habitat in the southern South Island of New Zealand during the Last Glacial Maximum (LGM). However, previous studies have lacked the power and geographical sampling to explicitly test this hypothesis using genetic data. Here we analyse 46 ancient mitochondrial genomes from Late Pleistocene and Holocene bones of the eastern moa from across their post-LGM distribution. Our results are consistent with a post-LGM increase in the population size and genetic diversity of eastern moa. We also demonstrate that genetic diversity was higher in eastern moa from the southern extent of their range, supporting the hypothesis that they expanded from a single glacial refugium following the LGM.

A two-step metagenomics approach for the identification and mitochondrial DNA contig assembly of vertebrate prey from the blood meals of common vampire bats (Desmodus rotundus)

The feeding behaviour of the sanguivorous common vampire bat (Desmodus rotundus) facilitates the transmission of pathogens that can impact both human and animal health. To formulate effective strategies in controlling the spread of diseases, there is a need to obtain information on which animals they feed on. One DNA-based approach, shotgun sequencing, can be used to obtain such information. Even though it is costly, shotgun sequencing can be used to simultaneously retrieve prey and vampire bat mitochondrial DNA for population studies within one round of sequencing. However, due to the challenges of analysing shotgun sequenced metagenomic data such as false negatives/positives and typically low proportion of reads mapped to diet items, shotgun sequencing has not been used for the identification of prey from common vampire bat blood meals. To overcome these challenges and generate longer mitochondrial contigs which could be useful for prey population studies, we shotgun sequenced common vampire bat blood meal samples (n = 8) and utilised a two-step metagenomic approach based on combining existing bioinformatic workflows (alignment and mtDNA contig assembly) to identify prey. After validating our results from detections made through metabarcoding, we accurately identified the common vampire bats’ prey in six out of eight samples without any false positives. We also generated prey mitochondrial contig lengths between 138 bp to 3231 bp (median = 770 bp, Q1 = 262 bp, Q3 = 1766 bp). This opens the potential to conduct phylogenetic and phylogeographic monitoring of elusive prey species in future studies, through the analyses of blood meal metagenomic data.

Echinococcus spp. and genotypes infecting humans in Tibet Autonomous Region of China: a molecular investigation with near-complete/complete mitochondrial sequences

Background

Molecular markers are essential to identify Echinococcus species and genotypes in areas with multiple Echinococcus species to understand their epidemiology and pathology. Tibet Autonomous Region (TAR) is one of the areas worst hit by echinococcosis. However, molecular epidemiology is still missing among echinococcosis patients in TAR. This research explored the Echinococcus species and genotypes infecting humans in TAR and the population diversity and the possible origin of G1 in TAR.

Methods

Cyst samples were collected in one echinococcosis-designated hospital in TAR. Echinococcus species and genotypes were identified through a maximum-likelihood approach with near-complete/complete mtDNA using IQ-TREE. Phylogenetic networks were built with PopART, and the phylogeographical diffusion pattern was identified using a Bayesian discrete phylogeographic method.

Results

Using phylogenetic trees made with near-complete/complete mtDNA obtained from 92 cysts from TAR patients, the Echinococcus species and genotypes infecting humans in TAR were identified as Echinococcus granulosus (s.s.) G1 (81, 88.04%), accounting for the majority, followed by G6 of the E. canadensis cluster (6, 6.52%), E. granulosus (s.s.) G3 (3, 3.26%), and E. multilocularis (2, 2.17%). An expansion trend and a possible recent bottleneck event were confirmed among the G1 samples in TAR. Adding the other near-complete mtDNA of G1 samples globally from the literature, we identified the possible phylogeographic origin of the G1 samples in TAR as Turkey.

Conclusions

Using near-complete/complete mtDNA sequences of Echinococcus spp. obtained from echinococcosis patients, a variety of Echinococcus species and genotypes infecting humans throughout TAR were identified. As far as we know, this is the first comprehensive molecular investigation of Echinococcus species and genotypes infecting humans throughout TAR. We identified, for the first time to our knowledge, the possible origin of the G1 in TAR. We also enriched the long mtDNA database of Echinococcus spp. and added two complete E. multilocularis mtDNA sequences from human patients. These findings will improve our knowledge of echinococcosis, help to refine the targeted echinococcosis control measures, and serve as a valuable baseline for monitoring the Echinococcus species and genotypes mutations and trends of the Echinococcus spp. population in TAR.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05199-6.

Probing the genomic limits of de-extinction in the Christmas Island rat

Three principal methods are under discussion as possible pathways to “true” de-extinction; i.e., back-breeding, cloning, and genetic engineering.^1,2 Of these, while the latter approach is most likely to apply to the largest number of extinct species, its potential is constrained by the degree to which the extinct species genome can be reconstructed. We explore this question using the extinct Christmas Island rat (Rattus macleari) as a model, an endemic rat species that was driven extinct between 1898 and 1908.3, 4, 5 We first re-sequenced its genome to an average of >60× coverage, then mapped it to the reference genomes of different Rattus species. We then explored how evolutionary divergence from the extant reference genome affected the fraction of the Christmas Island rat genome that could be recovered. Our analyses show that even when the extremely high-quality Norway brown rat (R. norvegicus) is used as a reference, nearly 5% of the genome sequence is unrecoverable, with 1,661 genes recovered at lower than 90% completeness, and 26 completely absent. Furthermore, we find the distribution of regions affected is not random, but for example, if 90% completeness is used as the cutoff, genes related to immune response and olfaction are excessively affected. Ultimately, our approach demonstrates the importance of applying similar analyses to candidates for de-extinction through genome editing in order to provide critical baseline information about how representative the edited form would be of the extinct species.

•

Evolutionary divergence limits the completeness of extinct species genomes

•

The extinct Christmas Island rat was re-sequenced to ca. 60× coverage

•

Nevertheless, 4.85% of the Norway brown rat genome remains absent after mapping

•

Absences are not random; immune response and olfaction are excessively affected

Transcriptome signature changes in the liver of a migratory passerine

The liver plays a principal role in avian migration. Here, we characterised the liver transcriptome of a long-distance migrant, the Northern Wheatear (Oenanthe oenanthe), sampled at different migratory stages, looking for molecular processes linked with adaptations to migration. The analysis of the differentially expressed genes suggested changes in the periods of the circadian rhythm, variation in the proportion of cells in G1/S cell-cycle stages and the putative polyploidization of this cell population. This may explain the dramatic increment in the liver's metabolic capacities towards migration. Additionally, genes involved in anti-oxidative stress, detoxification and innate immune responses, lipid metabolism, inflammation and angiogenesis were regulated. Lipophagy and lipid catabolism were active at all migratory stages and increased towards the fattening and fat periods, explaining the relevance of lipolysis in controlling steatosis and maintaining liver health. Our study clears the way for future functional studies regarding long-distance avian migration.

A chromosome-level genome of Antechinus flavipes provides a reference for an Australian marsupial genus with male death after mating

The 15 species of small carnivorous marsupials that comprise the genus Antechinus exhibit semelparity, a rare life-history strategy in mammals where synchronized death occurs after one breeding season. Antechinus males, but not females, age rapidly (demonstrate organismal senescence) during the breeding season and show promise as new animal models of ageing. Some antechinus species are also threatened or endangered. Here, we report a chromosome-level genome of a male yellow-footed antechinus Antechinus flavipes. The genome assembly has a total length of 3.2 Gb with a contig N50 of 51.8 Mb and a scaffold N50 of 636.7 Mb. We anchored and oriented 99.7% of the assembly on seven pseudochromosomes and found that repetitive DNA sequences occupy 51.8% of the genome. Draft genome assemblies of three related species in the subfamily Phascogalinae, two additional antechinus species (Antechinus argentus and A. arktos) and the iteroparous sister species Murexia melanurus, were also generated. Preliminary demographic analysis supports the hypothesis that climate change during the Pleistocene isolated species in Phascogalinae and shaped their population size. A transcriptomic profile across the A. flavipes breeding season allowed us to identify genes associated with aspects of the male die-off. The chromosome-level A. flavipes genome provides a steppingstone to understanding an enigmatic life-history strategy and a resource to assist the conservation of antechinuses.

The mutationathon highlights the importance of reaching standardization in estimates of pedigree-based germline mutation rates

In the past decade, several studies have estimated the human per-generation germline mutation rate using large pedigrees. More recently, estimates for various nonhuman species have been published. However, methodological differences among studies in detecting germline mutations and estimating mutation rates make direct comparisons difficult. Here, we describe the many different steps involved in estimating pedigree-based mutation rates, including sampling, sequencing, mapping, variant calling, filtering, and appropriately accounting for false-positive and false-negative rates. For each step, we review the different methods and parameter choices that have been used in the recent literature. Additionally, we present the results from a ‘Mutationathon,’ a competition organized among five research labs to compare germline mutation rate estimates for a single pedigree of rhesus macaques. We report almost a twofold variation in the final estimated rate among groups using different post-alignment processing, calling, and filtering criteria, and provide details into the sources of variation across studies. Though the difference among estimates is not statistically significant, this discrepancy emphasizes the need for standardized methods in mutation rate estimations and the difficulty in comparing rates from different studies. Finally, this work aims to provide guidelines for computational and statistical benchmarks for future studies interested in identifying germline mutations from pedigrees.

Transcriptomics de novo sequencing data of Messastrum gracile SE-MC4 under exponential and stationary growth stages

Messastrum gracile SE-MC4 is a non-model microalga exhibiting superior oil-accumulating abilities. However, biomass production in M. gracile SE-MC4 is limited due to low cell proliferation especially after prolonged cultivation under oil-inducing culture conditions. Present data consist of next generation RNA sequencing data of M. gracile SE-MC4 under exponential and stationary growth stages. RNA of six samples were extracted and sequenced with insert size of 100 bp paired-end strategy using BGISEQ-500 platform to produce a total of 59.64 Gb data with 314 million reads. Sequences were filtered and de novo assembled to form 53,307 number of gene sequences. Sequencing data were deposited in National Center for Biotechnology Information (NCBI) and can be accessed via BioProject ID PRJNA552165. This information can be used to enhance biomass production in M. gracile SE-MC4 and other microalgae aimed towards improving biodiesel development.

Gut Microbiota Linked with Reduced Fear of Humans in Red Junglefowl Has Implications for Early Domestication

Domestication of animals can lead to profound phenotypic modifications within short evolutionary time periods, and for many species behavioral selection is likely at the forefront of this process. Animal studies have strongly implicated that the gut microbiome plays a major role in host behavior and cognition through the microbiome-gut-brain axis. Consequently, herein, it is hypothesized that host gut microbiota may be one of the earliest phenotypes to change as wild animals were domesticated. Here, the gut microbiome community in two selected lines of red junglefowl that are selected for either high or low fear of humans up to eight generations is examined. Microbiota profiles reveal taxonomic differences in gut bacteria known to produce neuroactive compounds between the two selection lines. Gut-brain module analysis by means of genome-resolved metagenomics identifies enrichment in the microbial synthesis and degradation potential of metabolites associated with fear extinction and reduces anxiety-like behaviors in low fear fowls. In contrast, high fear fowls are enriched in gut-brain modules from the butyrate and glutamate pathways, metabolites associated with fear conditioning. Overall, the results identify differences in the composition and functional potential of the gut microbiota across selection lines that may provide insights into the mechanistic explanations of the domestication process.

Kouprey (Bos sauveli) genomes unveil polytomic origin of wild Asian Bos

The evolution of the genera Bos and Bison, and the nature of gene flow between wild and domestic species, is poorly understood, with genomic data of wild species being limited. We generated two genomes from the likely extinct kouprey (Bos sauveli) and analyzed them alongside other Bos and Bison genomes. We found that B. sauveli possessed genomic signatures characteristic of an independent species closely related to Bos javanicus and Bos gaurus. We found evidence for extensive incomplete lineage sorting across the three species, consistent with a polytomic diversification of the major ancestry in the group, potentially followed by secondary gene flow. Finally, we detected significant gene flow from an unsampled Asian Bos-like source into East Asian zebu cattle, demonstrating both that the full genomic diversity and evolutionary history of the Bos complex has yet to be elucidated and that museum specimens and ancient DNA are valuable resources to do so.

•

We generated two genomes from the likely extinct kouprey (Bos sauveli)

•

Extensive mt and nuclear-genome-wide incomplete lineage sorting across wild Asian Bos

•

Initial polytomic diversification of the wild Asian Bos—kouprey, banteng, and gaur

An efficient method for simultaneous species, individual, and sex identification via in-solution single nucleotide polymorphism capture from low-quality scat samples

Understanding predator population dynamics is important for conservation management because of the critical roles predators play within ecosystems. Noninvasive genetic sampling methods are useful for the study of predators like canids that can be difficult to capture or directly observe. Here, we introduce the FAECES* method (Fast and Accurate Enrichment of Canid Excrement for Species* and other analyses) which expands the toolbox for canid researchers and conservationists by using in-solution hybridization sequence capture to produce single nucleotide polymorphism (SNP) genotypes for multiple canid species from scat-derived DNA using a single enrichment. We designed a set of hybridization probes to genotype both coyotes (Canis latrans) and kit foxes (Vulpes macrotis) at hundreds of polymorphic SNP loci and we tested the probes on both tissues and field-collected scat samples. We enriched and genotyped by sequencing 52 coyote and 70 kit fox scats collected in and around a conservation easement in the Nevada Mojave Desert. We demonstrate that the FAECES* method produces genotypes capable of differentiating coyotes and kit foxes, identifying individuals and their sex, and estimating genetic diversity and effective population sizes, even using highly degraded, low-quantity DNA extracted from scat. We found that the study area harbours a large and diverse population of kit foxes and a relatively smaller population of coyotes. By replicating our methods in the future, conservationists can assess the impacts of management decisions on canid populations. The method can also be adapted and applied more broadly to enrich and sequence multiple loci from any species of interest using scat or other noninvasive genetic samples.

An accurate assignment test for extremely low-coverage whole-genome sequence data

Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.

Comparative Performance of the MGISEQ-2000 and Illumina X-Ten Sequencing Platforms for Paleogenomics

The MGISEQ-2000 sequencer is widely used in various omics studies, but the performance of this platform for paleogenomics has not been evaluated. We here compare the performance of MGISEQ-2000 with the Illumina X-Ten on ancient human DNA using four samples from 1750BCE to 60CE. We found there were only slight differences between the two platforms in most parameters (duplication rate, sequencing bias, θ, δS, and λ). MGISEQ-2000 performed well on endogenous rate and library complexity although X-Ten had a higher average base quality and lower error rate. Our results suggest that MGISEQ-2000 and X-Ten have comparable performance, and MGISEQ-2000 can be an alternative platform for paleogenomics sequencing.

Evolutionary history of the extinct Sardinian dhole

The Sardinian dhole (Cynotherium sardous)¹ was an iconic and unique canid species that was endemic to Sardinia and Corsica until it became extinct at the end of the Late Pleistocene.^2-5 Given its peculiar dental morphology, small body size, and high level of endemism, several extant canids have been proposed as possible relatives of the Sardinian dhole, including the Asian dhole and African hunting dog ancestor.^3,6-9 Morphometric analyses^3,6,8-12 have failed to clarify the evolutionary relationship with other canids.We sequenced the genome of a ca-21,100-year-old Sardinian dhole in order to understand its genomic history and clarify its phylogenetic position. We found that it represents a separate taxon from all other living canids from Eurasia, Africa, and North America, and that the Sardinian dhole lineage diverged from the Asian dhole ca 885 ka. We additionally detected historical gene flow between the Sardinian and Asian dhole lineages, which ended approximately 500-300 ka, when the land bridge between Sardinia and mainland Italy was already broken, severing their population connectivity. Our sample showed low genome-wide diversity compared to other extant canids-probably a result of the long-term isolation-that could have contributed to the subsequent extinction of the Sardinian dhole.

Ancient and modern genomes unravel the evolutionary history of the rhinoceros family

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Nasal Immunization With Small Molecule Mast Cell Activators Enhance Immunity to Co-Administered Subunit Immunogens

Mast cell activators are a novel class of mucosal vaccine adjuvants. The polymeric compound, Compound 48/80 (C48/80), and cationic peptide, Mastoparan 7 (M7) are mast cell activators that provide adjuvant activity when administered by the nasal route. However, small molecule mast cell activators may be a more cost-efficient adjuvant alternative that is easily synthesized with high purity compared to M7 or C48/80. To identify novel mast cell activating compounds that could be evaluated for mucosal vaccine adjuvant activity, we employed high-throughput screening to assess over 55,000 small molecules for mast cell degranulation activity. Fifteen mast cell activating compounds were down-selected to five compounds based on in vitro immune activation activities including cytokine production and cellular cytotoxicity, synthesis feasibility, and selection for functional diversity. These small molecule mast cell activators were evaluated for in vivo adjuvant activity and induction of protective immunity against West Nile Virus infection in BALB/c mice when combined with West Nile Virus envelope domain III (EDIII) protein in a nasal vaccine. We found that three of the five mast cell activators, ST101036, ST048871, and R529877, evoked high levels of EDIII-specific antibody and conferred comparable levels of protection against WNV challenge. The level of protection provided by these small molecule mast cell activators was comparable to the protection evoked by M7 (67%) but markedly higher than the levels seen with mice immunized with EDIII alone (no adjuvant 33%). Thus, novel small molecule mast cell activators identified by high throughput screening are as efficacious as previously described mast cell activators when used as nasal vaccine adjuvants and represent next-generation mast cell activators for evaluation in mucosal vaccine studies.

Historical population declines prompted significant genomic erosion in the northern and southern white rhinoceros (Ceratotherium simum)

Large vertebrates are extremely sensitive to anthropogenic pressure, and their populations are declining fast. The white rhinoceros (Ceratotherium simum) is a paradigmatic case: this African megaherbivore has suffered a remarkable decline in the last 150 years due to human activities. Its subspecies, the northern (NWR) and the southern white rhinoceros (SWR), however, underwent opposite fates: the NWR vanished quickly, while the SWR recovered after the severe decline. Such demographic events are predicted to have an erosive effect at the genomic level, linked to the extirpation of diversity, and increased genetic drift and inbreeding. However, there is little empirical data available to directly reconstruct the subtleties of such processes in light of distinct demographic histories. Therefore, we generated a whole-genome, temporal data set consisting of 52 resequenced white rhinoceros genomes, representing both subspecies at two time windows: before and during/after the bottleneck. Our data reveal previously unknown population structure within both subspecies, as well as quantifiable genomic erosion. Genome-wide heterozygosity decreased significantly by 10% in the NWR and 36% in the SWR, and inbreeding coefficients rose significantly by 11% and 39%, respectively. Despite the remarkable loss of genomic diversity and recent inbreeding it suffered, the only surviving subspecies, the SWR, does not show a significant accumulation of genetic load compared to its historical counterpart. Our data provide empirical support for predictions about the genomic consequences of shrinking populations, and our findings have the potential to inform the conservation efforts of the remaining white rhinoceroses.

Historical Demographic Processes Dominate Genetic Variation in Ancient Atlantic Cod Mitogenomes

Ancient DNA (aDNA) approaches have been successfully used to infer the long-term impacts of climate change, domestication, and human exploitation in a range of terrestrial species. Nonetheless, studies investigating such impacts using aDNA in marine species are rare. Atlantic cod (Gadus morhua), is an economically important species that has experienced dramatic census population declines during the last century. Here, we investigated 48 ancient mitogenomes from historical specimens obtained from a range of archeological excavations in northern Europe dated up to 6,500 BCE. We compare these mitogenomes to those of 496 modern conspecifics sampled across the North Atlantic Ocean and adjacent seas. Our results confirm earlier observations of high levels of mitogenomic variation and a lack of mutation-drift equilibrium—suggestive of population expansion. Furthermore, our temporal comparison yields no evidence of measurable mitogenomic changes through time. Instead, our results indicate that mitogenomic variation in Atlantic cod reflects past demographic processes driven by major historical events (such as oscillations in sea level) and subsequent gene flow rather than contemporary fluctuations in stock abundance. Our results indicate that historical and contemporaneous anthropogenic pressures such as commercial fisheries have had little impact on mitogenomic diversity in a wide-spread marine species with high gene flow such as Atlantic cod. These observations do not contradict evidence that overfishing has had negative consequences for the abundance of Atlantic cod and the importance of genetic variation in implementing conservation strategies. Instead, these observations imply that any measures toward the demographic recovery of Atlantic cod in the eastern Atlantic, will not be constrained by recent loss of historical mitogenomic variation.

The genomic origin of Zana of Abkhazia

Enigmatic phenomena have sparked the imagination of people around the globe into creating folkloric creatures. One prime example is Zana of Abkhazia (South Caucasus), a well-documented 19th century female who was captured living wild in the forest. Zana's appearance was sufficiently unusual, that she was referred to by locals as an Almasty-the analog of Bigfoot in the Caucasus. Although the exact location of Zana's burial site was unknown, the grave of her son, Khwit, was identified in 1971. The genomes of Khwit and the alleged Zana skeleton were sequenced to an average depth of ca. 3× using ancient DNA techniques. The identical mtDNA and parent-offspring relationship between the two indicated that the unknown woman was indeed Zana. Population genomic analyses demonstrated that Zana's immediate genetic ancestry can likely be traced to present-day East-African populations. We speculate that Zana might have had a genetic disorder such as congenital generalized hypertrichosis which could partially explain her strange behavior, lack of speech, and long body hair. Our findings elucidate Zana's unfortunate story and provide a clear example of how prejudices of the time led to notions of cryptic hominids that are still held and transmitted by some today.

Progress in forensic bone DNA analysis: Lessons learned from ancient DNA

Research on ancient and forensic DNA is related in many ways, and the two fields must deal with similar obstacles. Therefore, communication between these two communities has the potential to improve results in both research fields. Here, we present the insights gained in the ancient DNA community with regard to analyzing DNA from aged skeletal material and the potential use of the developed protocols in forensic work. We discuss the various steps, from choosing samples for DNA extraction to deciding between classical PCR amplification and massively parallel sequencing approaches. Based on the progress made in ancient DNA analyses combined with the requirements of forensic work, we suggest that there is substantial potential for incorporating ancient DNA approaches into forensic protocols, a process that has already begun to a considerable extent. However, taking full advantage of the experiences gained from ancient DNA work will require comparative studies by the forensic DNA community to tailor the methods developed for ancient samples to the specific needs of forensic studies and case work. If successful, in our view, the benefits for both communities would be considerable.

A Fast and Efficient Single-stranded Genomic Library Preparation Method Optimized for Ancient DNA

We present a protocol to prepare extracted DNA for sequencing on the Illumina sequencing platform that has been optimized for ancient and degraded DNA. Our approach, the Santa Cruz Reaction or SCR, uses directional splinted ligation of Illumina’s P5 and P7 adapters to convert natively single-stranded DNA and heat denatured double-stranded DNA into sequencing libraries in a single enzymatic reaction. To demonstrate its efficacy in converting degraded DNA molecules, we prepare 5 ancient DNA extracts into sequencing libraries using the SCR and 2 of the most commonly used approaches for preparing degraded DNA for sequencing: BEST, which targets and converts double-stranded DNA, and ssDNA2.0, which targets and converts single-stranded DNA. We then compare the efficiency with which each approach recovers unique molecules, or library complexity, given a standard amount of DNA input. We find that the SCR consistently outperforms the BEST protocol in recovering unique molecules and, despite its relative simplicity to perform and low cost per library, has similar performance to ssDNA2.0 across a wide range of DNA inputs. The SCR is a cost- and time-efficient approach that minimizes the loss of unique molecules and makes accessible a taxonomically, geographically, and a temporally broader sample of preserved remains for genomic analysis.

Rapid identification and metagenomics analysis of the adenovirus type 55 outbreak in Hubei using real-time and high-throughput sequencing platforms

The rise in human adenovirus (HAdV) infections poses a serious challenge to public health in China. Real-time (RT) sequencing provides solutions for achieving rapid pathogen identification during outbreaks, whereas high-throughput sequencing yields higher sequence accuracy. In the present study, we report the outcomes of applying nanopore and BGI platforms in the identification and genomic analysis of an HAdV outbreak in Hubei province, China in May of 2019. A mixed sample of nine nasopharyngeal swabs and one single sample were submitted to direct nanopore sequencing (MinION device), generating their first HAdV-55 reads within 13 and 20 min, respectively. The sequences were confirmed by RT-polymerase chain reaction (PCR). Ten HAdV-positive samples were further sequenced using next-generation high-throughput sequencing (BGISEQ-500 device). Phylogenetic analysis revealed that the outbreak strain had a close genetic relation to strains isolated in Sichuan province. Metagenomic analysis showed that HAdV-55 was not a dominant species in samples from which the whole HAdV-55 genome could not be assembled. The present results highlight the value of combining sequencing platforms and using mixed samples for nucleic acid enrichment in pathogen detection of infectious disease outbreaks.

Metagenomics: A viable tool for reconstructing herbivore diet

Metagenomics can generate data on the diet of herbivores, without the need for primer selection and PCR enrichment steps as is necessary in metabarcoding. Metagenomic approaches to diet analysis have remained relatively unexplored, requiring validation of bioinformatic steps. Currently, no metagenomic herbivore diet studies have utilized both chloroplast and nuclear markers as reference sequences for plant identification, which would increase the number of reads that could be taxonomically informative. Here, we explore how in silico simulation of metagenomic data sets resembling sequences obtained from faecal samples can be used to validate taxonomic assignment. Using a known list of sequences to create simulated data sets, we derived reliable identification parameters for taxonomic assignments of sequences. We applied these parameters to characterize the diet of western capercaillies (Tetrao urogallus) located in Norway, and compared the results with metabarcoding trnL P6 loop data generated from the same samples. Both methods performed similarly in the number of plant taxa identified (metagenomics 42 taxa, metabarcoding 43 taxa), with no significant difference in species resolution (metagenomics 24%, metabarcoding 23%). We further observed that while metagenomics was strongly affected by the age of faecal samples, with fresh samples outperforming old samples, metabarcoding was not affected by sample age. On the other hand, metagenomics allowed us to simultaneously obtain the mitochondrial genome of the western capercaillies, thereby providing additional ecological information. Our study demonstrates the potential of utilizing metagenomics for diet reconstruction but also highlights key considerations as compared to metabarcoding for future utilization of this technique.

Transcriptional Landscape of Vero E6 Cells during Early Swine Acute Diarrhea Syndrome Coronavirus Infection

Swine acute diarrhea syndrome coronavirus (SADS-CoV) is a newly emerged and highly pathogenic virus that is associated with fatal diarrhea disease in piglets, causing significant economic losses to the pig industry. At present, the research on the pathogenicity and molecular mechanisms of host-virus interactions of SADS-CoV are limited and remain poorly understood. Here, we investigated the global gene expression profiles of SADS-CoV-infected Vero E6 cells at 12, 18, and 24 h post-infection (hpi) using the RNA-sequencing. As a result, a total of 3324 differentially expressed genes (DEG) were identified, most of which showed a down-regulated expression pattern. Functional enrichment analyses indicated that the DEGs are mainly involved in signal transduction, cellular transcription, immune and inflammatory response, and autophagy. Collectively, our results provide insights into the changes in the cellular transcriptome during early infection of SADS-CoV and may provide information for further study of molecular mechanisms.

Efficient and stable metabarcoding sequencing data using a DNBSEQ-G400 sequencer validated by comprehensive community analyses

Metabarcoding is a widely used method for fast characterization of microbial communities in complex environmental samples. However, the selction of sequencing platform can have a noticeable effect on the estimated community composition. Here, we evaluated the metabarcoding performance of a DNBSEQ-G400 sequencer developed by MGI Tech using 16S and internal transcribed spacer (ITS) markers to investigate bacterial and fungal mock communities, as well as the ITS2 marker to investigate the fungal community of 1144 soil samples, with additional technical replicates. We show that highly accurate sequencing of bacterial and fungal communities is achievable using DNBSEQ-G400. Measures of diversity and correlation from soil metabarcoding showed that the results correlated highly with those of different machines of the same model, as well as between different sequencing modes (single-end 400 bp and paired-end 200 bp). Moderate, but significant differences were observed between results produced with different sequencing platforms (DNBSEQ-G400 and MiSeq); however, the highest differences can be caused by selecting different primer pairs for PCR amplification of taxonomic markers. These differences suggested that care is needed while jointly analyzing metabarcoding data from differenet experiments. This study demonstrated the high performance and accuracy of DNBSEQ-G400 for short-read metabarcoding of microbial communities. Our study also produced datasets to allow further investigation of microbial diversity.

TDRD7 participates in lens development and spermiogenesis by mediating autophagosome maturation

In humans, TDRD7 (tudor domain containing 7) mutations lead to a syndrome combining congenital cataracts (CCs) and non-obstructive azoospermia (NOA), characterized by abnormal lens development and spermiogenesis. However, the molecular mechanism underlying TDRD7's functions in eye and testicular development are still largely unknown. Here, we show that the depletion of this gene in mice and humans resulted in the accumulation of autophagosomes and the disruption of macroautophagic/autophagic flux. The disrupted autophagic flux in tdrd7-deficient mouse embryonic fibroblasts (MEFs) was caused by a failure of autophagosome fusion with lysosomes. Furthermore, transcriptome analysis and biochemical assays showed that TDRD7 might directly bind to Tbc1d20 mRNAs and downregulate its expression, which is a key regulator of autophagosome maturation, resulting in the disruption of autophagosome maturation. In addition, we provide evidence to show that TDRD7-mediated autophagosome maturation maintains lens transparency by facilitating the removal of damaged proteins and organelles from lens fiber cells and the biogenesis of acrosome. Altogether, our results showed that TDRD7 plays an essential role in the maturation of autophagosomes and that tdrd7 deletion results in eye defects and testicular abnormalities in mice, implicating disrupted autophagy might be the mechanism that contributes to lens development and spermiogenesis defects in human.Abbreviations: CB: chromatoid bodies; CC: congenital cataract; CTSD: cathepsin D; DMSO: dimethyl sulfoxide; LAMP1: lysosomal-associated membrane protein 1; LECs: lens epithelial cells; MAP1LC3/LC3/Atg8: microtubule-associated protein 1 light chain 3; MEFs: mouse embryonic fibroblasts; NOA: non-obstructive azoospermia; OFZ: organelle-free zone; RG: RNA granules; SQSTM1/p62: sequestosome 1; TBC1D20: TBC1 domain family member 20; TDRD7: tudor domain containing 7; TEM: transmission electron microscopy; WT: wild type.

Single-cell RNA-seq reveals novel mitochondria-related musculoskeletal cell populations during adult axolotl limb regeneration process

While the capacity to regenerate tissues or limbs is limited in mammals, including humans, axolotls are able to regrow entire limbs and major organs after incurring a wound. The wound blastema has been extensively studied in limb regeneration. However, due to the inadequate characterization of ECM and cell subpopulations involved in the regeneration process, the discovery of the key drivers for human limb regeneration remains unknown. In this study, we applied large-scale single-cell RNA sequencing to classify cells throughout the adult axolotl limb regeneration process, uncovering a novel regeneration-specific mitochondria-related cluster supporting regeneration through energy providing and the ECM secretion (COL2+) cluster contributing to regeneration through cell-cell interactions signals. We also discovered the dedifferentiation and re-differentiation of the COL1+/COL2+ cellular subpopulation and exposed a COL2-mitochondria subcluster supporting the musculoskeletal system regeneration. On the basis of these findings, we reconstructed the dynamic single-cell transcriptome of adult axolotl limb regenerative process, and identified the novel regenerative mitochondria-related musculoskeletal populations, which yielded deeper insights into the crucial interactions between cell clusters within the regenerative microenvironment.

A Long-Term Conserved Satellite DNA That Remains Unexpanded in Several Genomes of Characiformes Fish Is Actively Transcribed

Eukaryotic genomes contain large amounts of repetitive DNA sequences, such as tandemly repeated satellite DNAs (satDNAs). These sequences are highly dynamic and tend to be genus- or species-specific due to their particular evolutionary pathways, although there are few unusual cases of conserved satDNAs over long periods of time. Here, we used multiple approaches to reveal that an satDNA named CharSat01-52 originated in the last common ancestor of Characoidei fish, a superfamily within the Characiformes order, ∼140-78 Ma, whereas its nucleotide composition has remained considerably conserved in several taxa. We show that 14 distantly related species within Characoidei share the presence of this satDNA, which is highly amplified and clustered in subtelomeric regions in a single species (Characidium gomesi), while remained organized as small clusters in all the other species. Defying predictions of the molecular drive of satellite evolution, CharSat01-52 shows similar values of intra- and interspecific divergence. Although we did not provide evidence for a specific functional role of CharSat01-52, its transcriptional activity was demonstrated in different species. In addition, we identified short tandem arrays of CharSat01-52 embedded within single-molecule real-time long reads of Astyanax paranae (536 bp-3.1 kb) and A. mexicanus (501 bp-3.9 kb). Such arrays consisted of head-to-tail repeats and could be found interspersed with other sequences, inverted sequences, or neighbored by other satellites. Our results provide a detailed characterization of an old and conserved satDNA, challenging general predictions of satDNA evolution.