A fast and accurate method for SARS-CoV-2 genomic tracing.

To contain infectious diseases, it is crucial to determine the origin and transmission routes of the pathogen, as well as how the virus evolves. With the development of genome sequencing technology, genome epidemiology has emerged as a powerful approach for investigating the source and transmission of pathogens. In this study, we first presented the rationale for genomic tracing of SARS-CoV-2 and the challenges we currently face. Identifying the most genetically similar reference sequence to the query sequence is a critical step in genome tracing, typically achieved using either a phylogenetic tree or a sequence similarity search. However, these methods become inefficient or computationally prohibitive when dealing with tens of millions of sequences in the reference database, as we encountered during the COVID-19 pandemic. To address this challenge, we developed a novel genomic tracing algorithm capable of processing 6 million SARS-CoV-2 sequences in less than a minute. Instead of constructing a giant phylogenetic tree, we devised a weighted scoring system based on mutation characteristics to quantify sequences similarity. The developed method demonstrated superior performance compared to previous methods. Additionally, an online platform was developed to facilitate genomic tracing and visualization of the spatiotemporal distribution of sequences. The method will be a valuable addition to standard epidemiological investigations, enabling more efficient genomic tracing. Furthermore, the computational framework can be easily adapted to other pathogens, paving the way for routine genomic tracing of infectious diseases.

PySNV for complex intra-host variation detection.

MOTIVATION: Intra-host variants refer to genetic variations or mutations that occur within an individual host organism. These variants are typically studied in the context of viruses, bacteria, or other pathogens to understand the evolution of pathogens. Moreover, intra-host variants are also explored in the field of tumor biology and mitochondrial biology to characterize somatic mutations and inherited heteroplasmic mutations. Intra-host variants can involve long insertions, deletions, and combinations of different mutation types, which poses challenges in their identification. The performance of current methods in detecting of complex intra-host variants is unknown. RESULTS: First, we simulated a dataset comprising 10 samples with 1869 intra-host variants involving various mutation patterns and benchmarked current variant detection software. The results indicated that though current software can detect most variants with F1-scores between 0.76 and 0.97, their performance in detecting long indels and low frequency variants was limited. Thus, we developed a new software, PySNV, for the detection of complex intra-host variations. On the simulated dataset, PySNV successfully detected 1863 variant cases (F1-score: 0.99) and exhibited the highest Pearson correlation coefficient (PCC: 0.99) to the ground truth in predicting variant frequencies. The results demonstrated that PySNV delivered promising performance even for long indels and low frequency variants, while maintaining computational speed comparable to other methods. Finally, we tested its performance on SARS-CoV-2 replicate sequencing data and found that it reported 21% more variants compared to LoFreq, the best-performing benchmarked software, while showing higher consistency (62% over 54%) within replicates. The discrepancies mostly exist in low-depth regions and low frequency variants. AVAILABILITY AND IMPLEMENTATION: https://github.com/bnuLyndon/PySNV/.

Highly diverse sputum microbiota correlates with the disease severity in patients with community-acquired pneumonia: a longitudinal cohort study.

BACKGROUND: Community-acquired pneumonia (CAP) is a common and serious condition that can be caused by a variety of pathogens. However, much remains unknown about how these pathogens interact with the lower respiratory commensals, and whether any correlation exists between the dysbiosis of the lower respiratory microbiota and disease severity and prognosis. METHODS: We conducted a retrospective cohort study to investigate the composition and dynamics of sputum microbiota in patients diagnosed with CAP. In total, 917 sputum specimens were collected consecutively from 350 CAP inpatients enrolled in six hospitals following admission. The V3-V4 region of the 16 S rRNA gene was then sequenced. RESULTS: The sputum microbiota in 71% of the samples were predominately composed of respiratory commensals. Conversely, 15% of the samples demonstrated dominance by five opportunistic pathogens. Additionally, 5% of the samples exhibited sterility, resembling the composition of negative controls. Compared to non-severe CAP patients, severe cases exhibited a more disrupted sputum microbiota, characterized by the highly dominant presence of potential pathogens, greater deviation from a healthy state, more significant alterations during hospitalization, and sparser bacterial interactions. The sputum microbiota on admission demonstrated a moderate prediction of disease severity (AUC = 0.74). Furthermore, different pathogenic infections were associated with specific microbiota alterations. Acinetobacter and Pseudomonas were more abundant in influenza A infections, with Acinetobacter was also enriched in Klebsiella pneumoniae infections. CONCLUSION: Collectively, our study demonstrated that pneumonia may not consistently correlate with severe dysbiosis of the respiratory microbiota. Instead, the degree of microbiota dysbiosis was correlated with disease severity in CAP patients.

CircleBase: an integrated resource and analysis platform for human eccDNAs.

Rapid advances in high-throughput sequencing technologies have led to the discovery of thousands of extrachromosomal circular DNAs (eccDNAs) in the human genome. Loss-of-function experiments are difficult to conduct on circular and linear chromosomes, as they usually overlap. Hence, it is challenging to interpret the molecular functions of eccDNAs. Here, we present CircleBase (http://circlebase.maolab.org), an integrated resource and analysis platform used to curate and interpret eccDNAs in multiple cell types. CircleBase identifies putative functional eccDNAs by incorporating sequencing datasets, computational predictions, and manual annotations. It classifies them into six sections including targeting genes, epigenetic regulations, regulatory elements, chromatin accessibility, chromatin interactions, and genetic variants. The eccDNA targeting and regulatory networks are displayed by informative visualization tools and then prioritized. Functional enrichment analyses revealed that the top-ranked cancer cell eccDNAs were enriched in oncogenic pathways such as the Ras and PI3K-Akt signaling pathways. In contrast, eccDNAs from healthy individuals were not significantly enriched. CircleBase provides a user-friendly interface for searching, browsing, and analyzing eccDNAs in various cell/tissue types. Thus, it is useful to screen for potential functional eccDNAs and interpret their molecular mechanisms in human cancers and other diseases.

Water-Stable Hybrid Lead-Free Perovskite for Negative Temperature Coefficient Thermistors.

Negative temperature coefficient (NTC) thermistors feature higher sensitivities and faster response speeds and thereby have particular applications in many fields. However, current NTC thermistors are mostly based on inorganic ceramic materials, which show obvious drawbacks in material synthesis, property modulation, and flexible film fabrication. Herein, we report, for the first time, the promising application of an inorganic-organic hybrid NTC thermistor. A new lead-free hybrid iodo bismuthate [1,1',1″-(benzene-1,3,5-triyl)tris(3-methyl-1H-imidazol-3-ium)]Bi2I9 [denoted as (Me3TMP)Bi2I9] was synthesized by a "double-free" strategy. (Me3TMP)Bi2I9 features a lead-free binuclear bismuth iodine anion charge compensated by a "classic hydrogen-bond-free" cation. (Me3TMP)Bi2I9 exhibits remarkable stability in water and UV light irradiation and shows the largest temperature sensitivity coefficient among all reported NTC materials. Theoretical calculation and detailed structural analysis disclose that the seriously distorted (BiI6) octahedra are responsible for the intriguing NTC effect for (Me3TMP)Bi2I9.

Toward efficient and high-fidelity metagenomic data from sub-nanogram DNA: evaluation of library preparation and decontamination methods.

BACKGROUND: Shotgun metagenomic sequencing has greatly expanded the understanding of microbial communities in various biological niches. However, it is still challenging to efficiently convert sub-nanogram DNA to high-quality metagenomic libraries and obtain high-fidelity data, hindering the exploration of niches with low microbial biomass. RESULTS: To cope with this challenge comprehensively, we evaluated the performance of various library preparation methods on 0.5 pg-5 ng synthetic microbial community DNA, characterized contaminants, and further applied different in silico decontamination methods. First, we discovered that whole genome amplification prior to library construction led to worse outcomes than preparing libraries directly. Among different non-WGA-based library preparation methods, we found the endonuclease-based method being generally good for different amounts of template and the tagmentation-based method showing specific advantages with 0.5 pg template, based on evaluation metrics including fidelity, proportion of designated reads, and reproducibility. The load of contaminating DNA introduced by library preparation varied from 0.01 to 15.59 pg for different kits and accounted for 0.05 to 45.97% of total reads. A considerable fraction of the contaminating reads were mapped to human commensal and pathogenic microbes, thus potentially leading to erroneous conclusions in human microbiome studies. Furthermore, the best performing in silico decontamination method in our evaluation, Decontam-either, was capable of recovering the real microbial community from libraries where contaminants accounted for less than 10% of total reads, but not from libraries with heavy and highly varied contaminants. CONCLUSIONS: This study demonstrates that high-quality metagenomic data can be obtained from samples with sub-nanogram microbial DNA by combining appropriate library preparation and in silico decontamination methods and provides a general reference for method selection for samples with varying microbial biomass.

The transmission of human mitochondrial DNA in four-generation pedigrees.

Most of the pathogenic variants in mitochondrial DNA (mtDNA) exist in a heteroplasmic state (coexistence of mutant and wild-type mtDNA). Understanding how mtDNA is transmitted is crucial for predicting mitochondrial disease risk. Previous studies were based mainly on two-generation pedigree data, which are limited by the randomness in a single transmission. In this study, we analyzed the transmission of heteroplasmies in 16 four-generation families. First, we found that 57.8% of the variants in the great grandmother were transmitted to the fourth generation. The direction and magnitude of the frequency change during transmission appeared to be random. Moreover, no consistent correlation was identified between the frequency changes among the continuous transmissions, suggesting that most variants were functionally neutral or mildly deleterious and thus not subject to strong natural selection. Additionally, we found that the frequency of one nonsynonymous variant (m.15773G>A) showed a consistent increase in one family, suggesting that this variant may confer a fitness advantage to the mitochondrion/cell. We also estimated the effective bottleneck size during transmission to be 21-71. In summary, our study demonstrates the advantages of multigeneration data for studying the transmission of mtDNA for shedding new light on the dynamics of the mutation frequency in successive generations.

False Coronavirus Disease 2019 Cases due to Contamination by Inactivated Virus Vaccine.

A false-positive severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reverse-transcription polymerase chain reaction result can lead to unnecessary public health measures. We report 2 individuals whose respiratory specimens were contaminated by an inactivated SARS-CoV-2 vaccine strain (CoronaVac), likely at vaccination premises. Incidentally, whole genome sequencing of CoronaVac showed adaptive deletions on the spike protein, which do not result in observable changes of antigenicity.

Dynamics of the Upper Respiratory Tract Microbiota and Its Association with Mortality in COVID-19.

Rationale: Alteration of human respiratory microbiota had been observed in coronavirus disease (COVID-19). How the microbiota is associated with the prognosis in COVID-19 is unclear. Objectives: To characterize the feature and dynamics of the respiratory microbiota and its associations with clinical features in patients with COVID-19. Methods: We conducted metatranscriptome sequencing on 588 longitudinal oropharyngeal swab specimens collected from 192 patients with COVID-19 (including 39 deceased patients) and 95 healthy controls from the same geographic area. Meanwhile, the concentration of 27 cytokines and chemokines in plasma was measured for patients with COVID-19. Measurements and Main Results: The upper respiratory tract (URT) microbiota in patients with COVID-19 differed from that in healthy controls, whereas deceased patients possessed a more distinct microbiota, both on admission and before discharge/death. The alteration of URT microbiota showed a significant correlation with the concentration of proinflammatory cytokines and mortality. Specifically, Streptococcus-dominated microbiota was enriched in recovered patients, and showed high temporal stability and resistance against pathogens. In contrast, the microbiota in deceased patients was more susceptible to secondary infections and became more deviated from the norm after admission. Moreover, the abundance of S. parasanguinis on admission was significantly correlated with prognosis in nonsevere patients (lower vs. higher abundance, odds ratio, 7.80; 95% CI, 1.70-42.05). Conclusions: URT microbiota dysbiosis is a remarkable manifestation of COVID-19; its association with mortality suggests it may reflect the interplay between pathogens, symbionts, and the host immune status. Whether URT microbiota could be used as a biomarker for diagnosis and prognosis of respiratory diseases merits further investigation.

OncoBase: a platform for decoding regulatory somatic mutations in human cancers.

Whole-exome and whole-genome sequencing have revealed millions of somatic mutations associated with different human cancers, and the vast majority of them are located outside of coding sequences, making it challenging to directly interpret their functional effects. With the rapid advances in high-throughput sequencing technologies, genome-scale long-range chromatin interactions were detected, and distal target genes of regulatory elements were determined using three-dimensional (3D) chromatin looping. Herein, we present OncoBase (http://www.oncobase.biols.ac.cn/), an integrated database for annotating 81 385 242 somatic mutations in 68 cancer types from more than 120 cancer projects by exploring their roles in distal interactions between target genes and regulatory elements. OncoBase integrates local chromatin signatures, 3D chromatin interactions in different cell types and reconstruction of enhancer-target networks using state-of-the-art algorithms. It employs informative visualization tools to display the integrated local and 3D chromatin signatures and effects of somatic mutations on regulatory elements. Enhancer-promoter interactions estimated from chromatin interactions are integrated into a network diffusion system that quantitatively prioritizes somatic mutations and target genes from a large pool. Thus, OncoBase is a useful resource for the functional annotation of regulatory noncoding regions and systematically benchmarking the regulatory effects of embedded noncoding somatic mutations in human carcinogenesis.

Genomic Diversity of Severe Acute Respiratory Syndrome-Coronavirus 2 in Patients With Coronavirus Disease 2019.

BACKGROUND: A novel coronavirus (CoV), severe acute respiratory syndrome (SARS)-CoV-2, has infected >75 000 individuals and spread to >20 countries. It is still unclear how fast the virus evolved and how it interacts with other microorganisms in the lung. METHODS: We have conducted metatranscriptome sequencing for bronchoalveolar lavage fluid samples from 8 patients with SARS-CoV-2, and also analyzed data from 25 patients with community-acquired pneumonia (CAP), and 20 healthy controls for comparison. RESULTS: The median number of intrahost variants was 1-4 in SARS-CoV-2-infected patients, ranged from 0 to 51 in different samples. The distribution of variants on genes was similar to those observed in the population data. However, very few intrahost variants were observed in the population as polymorphisms, implying either a bottleneck or purifying selection involved in the transmission of the virus, or a consequence of the limited diversity represented in the current polymorphism data. Although current evidence did not support the transmission of intrahost variants in a possible person-to-person spread, the risk should not be overlooked. Microbiotas in SARS-CoV-2-infected patients were similar to those in CAP, either dominated by the pathogens or with elevated levels of oral and upper respiratory commensal bacteria. CONCLUSION: SARS-CoV-2 evolves in vivo after infection, which may affect its virulence, infectivity, and transmissibility. Although how the intrahost variant spreads in the population is still elusive, it is necessary to strengthen the surveillance of the viral evolution in the population and associated clinical changes.

Transmission of human mtDNA heteroplasmy in the Genome of the Netherlands families: support for a variable-size bottleneck.

Although previous studies have documented a bottleneck in the transmission of mtDNA genomes from mothers to offspring, several aspects remain unclear, including the size and nature of the bottleneck. Here, we analyze the dynamics of mtDNA heteroplasmy transmission in the Genomes of the Netherlands (GoNL) data, which consists of complete mtDNA genome sequences from 228 trios, eight dizygotic (DZ) twin quartets, and 10 monozygotic (MZ) twin quartets. Using a minor allele frequency (MAF) threshold of 2%, we identified 189 heteroplasmies in the trio mothers, of which 59% were transmitted to offspring, and 159 heteroplasmies in the trio offspring, of which 70% were inherited from the mothers. MZ twin pairs exhibited greater similarity in MAF at heteroplasmic sites than DZ twin pairs, suggesting that the heteroplasmy MAF in the oocyte is the major determinant of the heteroplasmy MAF in the offspring. We used a likelihood method to estimate the effective number of mtDNA genomes transmitted to offspring under different bottleneck models; a variable bottleneck size model provided the best fit to the data, with an estimated mean of nine individual mtDNA genomes transmitted. We also found evidence for negative selection during transmission against novel heteroplasmies (in which the minor allele has never been observed in polymorphism data). These novel heteroplasmies are enhanced for tRNA and rRNA genes, and mutations associated with mtDNA diseases frequently occur in these genes. Our results thus suggest that the female germ line is able to recognize and select against deleterious heteroplasmies.

Age-Related and Heteroplasmy-Related Variation in Human mtDNA Copy Number.

The mitochondrial (mt) genome is present in many copies in human cells, and intra-individual variation in mtDNA sequences is known as heteroplasmy. Recent studies found that heteroplasmies are highly tissue-specific, site-specific, and allele-specific, however the functional implications have not been explored. This study investigates variation in mtDNA copy numbers (mtCN) in 12 different tissues obtained at autopsy from 152 individuals (ranging in age from 3 days to 96 years). Three different methods to estimate mtCN were compared: shotgun sequencing (in 4 tissues), capture-enriched sequencing (in 12 tissues) and droplet digital PCR (ddPCR, in 2 tissues). The highest precision in mtCN estimation was achieved using shotgun sequencing data. However, capture-enrichment data provide reliable estimates of relative (albeit not absolute) mtCNs. Comparisons of mtCN from different tissues of the same individual revealed that mtCNs in different tissues are, with few exceptions, uncorrelated. Hence, each tissue of an individual seems to regulate mtCN in a tissue-related rather than an individual-dependent manner. Skeletal muscle (SM) samples showed an age-related decrease in mtCN that was especially pronounced in males, while there was an age-related increase in mtCN for liver (LIV) samples. MtCN in SM samples was significantly negatively correlated with both the total number of heteroplasmic sites and with minor allele frequency (MAF) at two heteroplasmic sites, 408 and 16327. Heteroplasmies at both sites are highly specific for SM, accumulate with aging and are part of functional elements that regulate mtDNA replication. These data support the hypothesis that selection acting on these heteroplasmic sites is reducing mtCN in SM of older individuals.

Potential effects of exploiting the Yunfu pyrite mine (southern China) on soil: evidence from analyzing trace elements in surface soil.

Trace element contamination caused by mining is a serious environmental problem. The potential effects of exploiting the Yunfu pyrite mine (southern China) on soil were investigated by determining trace elements in 56 surface soil samples from the vicinity of the Yunfu pyrite mine. The samples were acid dissolved and measured by an inductively coupled plasma mass spectrometry (ICP-MS). Principal component analysis and hierarchical cluster analysis were used to identify factors influencing the trace element contents and possible sources of the trace elements. The degree of trace element pollution was determined using the geological accumulation index Igeo. Monte Carlo simulations were used to assess the health risks posed. The results show that (1) six factors (parent material, mining activities, ore composition, rainfall, terrain, and other inputs) strongly affected the trace element contents of the soil samples. (2) There were three groups of trace elements, according to their possible sources. One group (Cs, Ga, Ge, Hf, Nb, Rb, Ta, Th, Ti, U, and Zr) mainly originated in parent rocks. Another group (Cr, Ni, Sr, and V) was mainly supplied by industrial plants and traffic emissions. The third group (Ba, Co, Cu, Mn, Pb, and Zn) was mainly supplied through pyrite ore exploitation processes. (3) Some samples were slightly to moderately polluted with Cs, Ga, Ge, Nb, Rb, Ta, and Ti. Most samples were moderately to highly polluted with Ba, Co, Cu, Mn, Pb, and Zn. (4) Trace elements in soil pose strong non-carcinogenic and carcinogenic health risks to people (particularly children) living near the Yunfu pyrite mine.

All-fiber radially/azimuthally polarized lasers based on mode coupling of tapered fibers.

We demonstrate a mode converter with an insertion loss of 0.36 dB based on mode coupling of tapered single-mode and two-mode fibers, and realize all-fiber flexible cylindrical vector lasers at 1550 nm. Attributing to the continuous distribution of a tangential electric field at taper boundaries, the laser is switchable between the radially and azimuthally polarized states by adjusting the input polarization. In the temporal domain, the operation is controllable among continuous-wave, Q-switched, and mode-locked statuses by changing the saturable absorber or pump strength. The duration of Q-switched radially/azimuthally polarized laser spans from 10.4/10.8 to 6/6.4 µs at the pump range of 38 to 58 mW, while that of the mode-locked pulse varies from 39.2/31.9 to 5.6/5.2 ps by controlling the laser bandwidth. The proposed laser combines the features of a cylindrical vector beam, a fiber laser, and an ultrafast pulse, providing a special and cost-effective source for practical applications.

Extensive tissue-related and allele-related mtDNA heteroplasmy suggests positive selection for somatic mutations.

Heteroplasmy in human mtDNA may play a role in cancer, other diseases, and aging, but patterns of heteroplasmy variation across different tissues have not been thoroughly investigated. Here, we analyzed complete mtDNA genome sequences at ∼3,500× average coverage from each of 12 tissues obtained at autopsy from each of 152 individuals. We identified 4,577 heteroplasmies (with an alternative allele frequency of at least 0.5%) at 393 positions across the mtDNA genome. Surprisingly, different nucleotide positions (nps) exhibit high frequencies of heteroplasmy in different tissues, and, moreover, heteroplasmy is strongly dependent on the specific consensus allele at an np. All of these tissue-related and allele-related heteroplasmies show a significant age-related accumulation, suggesting positive selection for specific alleles at specific positions in specific tissues. We also find a highly significant excess of liver-specific heteroplasmies involving nonsynonymous changes, most of which are predicted to have an impact on protein function. This apparent positive selection for reduced mitochondrial function in the liver may reflect selection to decrease damaging byproducts of liver mitochondrial metabolism (i.e., "survival of the slowest"). Overall, our results provide compelling evidence for positive selection acting on some somatic mtDNA mutations.

Early Austronesians: into and out of Taiwan.

A Taiwan origin for the expansion of the Austronesian languages and their speakers is well supported by linguistic and archaeological evidence. However, human genetic evidence is more controversial. Until now, there had been no ancient skeletal evidence of a potential Austronesian-speaking ancestor prior to the Taiwan Neolithic ~6,000 years ago, and genetic studies have largely ignored the role of genetic diversity within Taiwan as well as the origins of Formosans. We address these issues via analysis of a complete mitochondrial DNA genome sequence of an ~8,000-year-old skeleton from Liang Island (located between China and Taiwan) and 550 mtDNA genome sequences from 8 aboriginal (highland) Formosan and 4 other Taiwanese groups. We show that the Liangdao Man mtDNA sequence is closest to Formosans, provides a link to southern China, and has the most ancestral haplogroup E sequence found among extant Austronesian speakers. Bayesian phylogenetic analysis allows us to reconstruct a history of early Austronesians arriving in Taiwan in the north ~6,000 years ago, spreading rapidly to the south, and leaving Taiwan ~4,000 years ago to spread throughout Island Southeast Asia, Madagascar, and Oceania.

Improved calibration of the human mitochondrial clock using ancient genomes.

Reliable estimates of the rate at which DNA accumulates mutations (the substitution rate) are crucial for our understanding of the evolution and past demography of virtually any species. In humans, there are considerable uncertainties around these rates, with substantial variation among recent published estimates. Substitution rates have traditionally been estimated by associating dated events to the root (e.g., the divergence between humans and chimpanzees) or to internal nodes in a phylogenetic tree (e.g., first entry into the Americas). The recent availability of ancient mitochondrial DNA sequences allows for a more direct calibration by assigning the age of the sequenced samples to the tips within the human phylogenetic tree. But studies also vary greatly in the methodology employed and in the sequence panels analyzed, making it difficult to tease apart the causes for the differences between previous estimates. To clarify this issue, we compiled a comprehensive data set of 350 ancient and modern human complete mitochondrial DNA genomes, among which 146 were generated for the purpose of this study and estimated substitution rates using calibrations based both on dated nodes and tips. Our results demonstrate that, for the same data set, estimates based on individual dated tips are far more consistent with each other than those based on nodes and should thus be considered as more reliable.