A proof of concept for a targeted enrichment approach to the simultaneous detection and characterization of rickettsial pathogens from clinical specimens

Infection with either Rickettsia prowazekii or Orientia tsutsugamushi is common, yet diagnostic capabilities are limited due to the short window for positive identification. Until now, although targeted enrichment had been applied to increase sensitivity of sequencing-based detection for various microorganisms, it had not been applied to sequencing of R. prowazekii in clinical samples. Additionally, hybridization-based targeted enrichment strategies had only scarcely been applied to qPCR of any pathogens in clinical samples. Therefore, we tested a targeted enrichment technique as a proof of concept and found that it dramatically reduced the limits of detection of these organisms by both qPCR and high throughput sequencing. The enrichment methodology was first tested in contrived clinical samples with known spiked-in concentrations of R. prowazekii and O. tsutsugamushi DNA. This method was also evaluated using clinical samples, resulting in the simultaneous identification and characterization of O. tsutsugamushi directly from clinical specimens taken from sepsis patients. We demonstrated that the targeted enrichment technique is helpful by lowering the limit of detection, not only when applied to sequencing, but also when applied to qPCR, suggesting the technique could be applied more broadly to include other assays and/or microbes for which there are limited diagnostic or detection modalities.

Enricherator: A Bayesian Method for Inferring Regularized Genome-wide Enrichments from Sequencing Count Data

A pervasive question in biological research studying gene regulation, chromatin structure, or genomics is where, and to what extent, does a signal of interest arise genome-wide? This question is addressed using a variety of methods relying on high-throughput sequencing data as their final output, including ChIP-seq for protein-DNA interactions,1 GapR-seq for measuring supercoiling,2 and HBD-seq or DRIP-seq for R-loop positioning.3,4 Current computational methods to calculate genome-wide enrichment of the signal of interest usually do not properly handle the count-based nature of sequencing data, they often do not make use of the local correlation structure of sequencing data, and they do not apply any regularization of enrichment estimates. This can result in unrealistic estimates of the true underlying biological enrichment of interest, unrealistically low estimates of confidence in point estimates of enrichment (or no estimates of confidence at all), unrealistic gyrations in enrichment estimates at very close (<10 bp) genomic loci due to noise inherent in sequencing data, and in a multiple-hypothesis testing problem during interpretation of genome-wide enrichment estimates. We developed a tool called Enricherator to infer genome-wide enrichments from sequencing count data. Enricherator uses the variational Bayes algorithm to fit a generalized linear model to sequencing count data and to sample from the approximate posterior distribution of enrichment estimates (https://github.com/jwschroeder3/enricherator). Enrichments inferred by Enricherator more precisely identify known binding sites in cases where low coverage between binding sites leads to false-positive peak calls in these noisy regions of the genome; these benefits extend to published datasets.

Biogeographical and biodiversity patterns of planktonic microeukaryotes along the tropical western to eastern Pacific Ocean transect revealed by metabarcoding

Microeukaryotic plankton (0.2-200 µm), which are morphologically and genetically highly diverse, play a crucial role in ocean productivity and carbon consumption. The Pacific Ocean (PO), one of the world's largest oligotrophic regions, remains largely unexplored in terms of the biogeography and biodiversity of microeukaryotes based on large-scale sampling. We investigated the horizontal distribution of microeukaryotes along a 16,000 km transect from the west to the east of the PO. The alpha diversity indices showed a distinct decreasing trend from west to east, which was highly correlated with water temperature. The microeukaryotic community, which was clustered into the western, central, and eastern PO groups, displayed a significant distance-decay relationship. Syndiniales, a lineage of parasitic dinoflagellates, was ubiquitously distributed along the transect and dominated the community in terms of both sequence and zero-radius operational taxonomic unit (ZOTU) proportions. The prevailing dominance of Syndiniales-affiliated ZOTUs and their close associations with dinoflagellates, diatoms, and radiolarians, as revealed by SparCC correlation analysis, suggested that parasitism may be an important trophic strategy in the surface waters of the PO. Geographical distance and temperature were the most important environmental factors that significantly correlated with community structure. Overall, our study sheds more light on the distribution pattern of both alpha and beta diversities of microeukaryotic communities and highlighted the importance of parasitisms by Syndiniales across the tropical PO.IMPORTANCEUnderstanding the biogeographical and biodiversity patterns of microeukaryotic communities is essential to comprehending their roles in biogeochemical cycling. In this study, planktonic microeukaryotes were collected along a west-to-east Pacific Ocean transect (ca. 16,000 km). Our study revealed that the alpha diversity indices were highly correlated with water temperature, and the microeukaryotic communities displayed a distinct geographical distance-driven pattern. The predominance of the parasitic dinoflagellate lineage Syndiniales and their close relationship with other microeukaryotic groups suggest that parasitism may be a crucial survival strategy for microeukaryotes in the surface waters of the Pacific Ocean. Our findings expand our understanding of the biodiversity and biogeographical pattern of microeukaryotes and highlight the significance of parasitic Syndiniales in the surface ocean.

Unraveling the shift in bacterial communities profile grown in sediments co-contaminated with chlorolignin waste of pulp-paper mill by metagenomics approach

Pulp-paper mills (PPMs) are known for consistently generating a wide variety of pollutants, that are often unidentified and highly resistant to environmental degradation. The current study aims to investigate the changes in the indigenous bacterial communities profile grown in the sediment co-contaminated with organic and inorganic pollutants discharged from the PPMs. The two sediment samples, designated PPS-1 and PPS-2, were collected from two different sites. Physico-chemical characterization of PPS-1 and PPS-2 revealed the presence of heavy metals (mg kg-1) like Cu (0.009-0.01), Ni (0.005-0.002), Mn (0.078-0.056), Cr (0.015-0.009), Pb (0.008-0.006), Zn (0.225-0.086), Fe (2.124-0.764), Al (3.477-22.277), and Ti (99.792-45.012) along with high content of chlorophenol, and lignin. The comparative analysis of organic pollutants in sediment samples using gas chromatography-mass spectrometry (GC-MS) revealed the presence of major highly refractory compounds, such as stigmasterol, β-sitosterol, hexadecanoic acid, octadecanoic acid; 2,4-di-tert-butylphenol; heptacosane; dimethyl phthalate; hexachlorobenzene; 1-decanol,2-hexyl; furane 2,5-dimethyl, etc in sediment samples which are reported as a potential toxic compounds. Simultaneously, high-throughput sequencing targeting the V3-V4 hypervariable region of the 16S rRNA genes, resulted in the identification of 1,249 and 1,345 operational taxonomic units (OTUs) derived from a total of 115,665 and 119,386 sequences read, in PPS-1 and PPS-2, respectively. Analysis of rarefaction curves indicated a diversity in OTU abundance between PPS-1 (1,249 OTUs) and PPS-2 (1,345 OTUs). Furthermore, taxonomic assignment of metagenomics sequence data showed that Proteobacteria (55.40%; 56.30%), Bacteoidetes (11.30%; 12.20%), and Planctomycetes (5.40%; 4.70%) were the most abundant phyla; Alphproteobacteria (20.50%; 23.50%), Betaproteobacteria (16.00%; 12.30%), and Gammaproteobacteria were the most recorded classes in PPS-1 and PPS-2, respectively. At the genus level, Thiobacillus (7.60%; 4.50%) was the most abundant genera grown in sediment samples. The results indicate significant differences in both the diversity and relative abundance of taxa in the bacterial communities associated with PPS-2 when compared to PPS-1. This study unveils key insights into contaminant characteristics and shifts in bacterial communities within contaminated environments. It highlights the potential for developing efficient bioremediation techniques to restore ecological balance in pulp-paper mill waste-polluted areas, stressing the importance of identifying a significant percentage of unclassified genera and species to explore novel genes.

Comparative analysis of the diversity of symbionts in fat body of long- and short-winged brown planthoppers

The microbial community structure plays an important role in the internal environment of brown planthopper (BPH), Nilaparvata lugens (Hemiptera: Delphacidae), which is an indispensable part to reflect the internal environment of BPH. Wing dimorphism is a strategy for balancing flight and reproduction of insects. Here, quantitative fluorescence PCR was used to analyse the number and changes of the symbionts in the fat body of long- and short-winged BPHs at different developmental stages. A metagenomic library was constructed based on the 16 S rRNA sequence and internal transcribed spacer sequence for high-throughput sequencing, to analyze the community structure and population number of the symbionts of long- and short-winged BPHs, and to make functional prediction. This study enriches the connotation of BPH symbionts, and laid a theoretical foundation for the subsequent study of BPH-symbionts interaction and the function of symbionts in the host.

Unraveling the Physicochemical Properties and Bacterial Communities in Rabbit Meat during Chilled Storage

The freshness and bacterial communities of fresh and salted rabbit meat during 8 days of refrigerated storage at 4 °C were evaluated. The results showed that the addition of 2% salt significantly changed the color of meat, of which the lightness (L*), redness (a*), and yellowness (b*) were lower than that of fresh meat over time. The pH of all samples increased during storage, and meat with salt addition had lower values in comparison to fresh samples over time. The total volatile base nitrogen (TVB-N) concentration increased rapidly in salt-treated meat but was significantly (p < 0.05) lower than that in meat without salt added before 6 days. Over time, the content of thiobarbituric acid reactive substances (TBARS) showed a progressive trend, but a rapid increase occurred in salted meat. High-throughput sequencing showed that the microflora of each sample had a positive trend in alpha diversity and a negative trend in beta diversity. Bacterial taxonomic analysis indicated that the initial microbial flora for chilled rabbit meat was dominated by Shigaella, Bacteroides, and Lactococcus, and the population of Brochothrix and Psychrobacter increased over time and became the dominant spoilage bacterium. In particular, the addition of salt significantly reduced the abundance of Psychrobacter and Brochothrix. These findings might provide valuable information regarding the quality monitoring of rabbit meat during chilled storage.

Supplemented Infant Formula and Human Breast Milk Show Similar Patterns in Modulating Infant Microbiota Composition and Function In Vitro

The gut microbiota of healthy breastfed infants is often dominated by bifidobacteria. In an effort to mimic the microbiota of breastfed infants, modern formulas are fortified with bioactive and bifidogenic ingredients. These ingredients promote the optimal health and development of infants as well as the development of the infant microbiota. Here, we used INFOGEST and an in vitro batch fermentation model to investigate the gut health-promoting effects of a commercial infant formula supplemented with a blend containing docosahexaenoic acid (DHA) (20 mg/100 kcal), polydextrose and galactooligosaccharides (PDX/GOS) (4 g/L, 1:1 ratio), milk fat globule membrane (MFGM) (5 g/L), lactoferrin (0.6 g/L), and Bifidobacterium animalis subsp. lactis, BB-12 (BB-12) (106 CFU/g). Using fecal inoculates from three healthy infants, we assessed microbiota changes, the bifidogenic effect, and the short-chain fatty acid (SCFA) production of the supplemented test formula and compared those with data obtained from an unsupplemented base formula and from the breast milk control. Our results show that even after INFOGEST digestion of the formula, the supplemented formula can still maintain its bioactivity and modulate infants' microbiota composition, promote faster bifidobacterial growth, and stimulate production of SCFAs. Thus, it may be concluded that the test formula containing a bioactive blend promotes infant gut microbiota and SCFA profile to something similar, but not identical to those of breastfed infants.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Genotype Combinations Drive Variability in the Microbiome Configuration of the Rhizosphere of Maize/Bean Intercropping System

In an intercropping system, the interplay between cereals and legumes, which is strongly driven by the complementarity of below-ground structures and their interactions with the soil microbiome, raises a fundamental query: Can different genotypes alter the configuration of the rhizosphere microbial communities? To address this issue, we conducted a field study, probing the effects of intercropping and diverse maize (Zea mays L.) and bean (Phaseolus vulgaris L., Phaseolus coccineus L.) genotype combinations. Through amplicon sequencing of bacterial 16S rRNA genes from rhizosphere samples, our results unveil that the intercropping condition alters the rhizosphere bacterial communities, but that the degree of this impact is substantially affected by specific genotype combinations. Overall, intercropping allows the recruitment of exclusive bacterial species and enhances community complexity. Nevertheless, combinations of maize and bean genotypes determine two distinct groups characterized by higher or lower bacterial community diversity and complexity, which are influenced by the specific bean line associated. Moreover, intercropped maize lines exhibit varying propensities in recruiting bacterial members with more responsive lines showing preferential interactions with specific microorganisms. Our study conclusively shows that genotype has an impact on the rhizosphere microbiome and that a careful selection of genotype combinations for both species involved is essential to achieve compatibility optimization in intercropping.

Gut microbiota changes in patients with hypertension: A systematic review and meta-analysis

Hypertension is a major public health issue worldwide. The imbalance of gut microbiota is thought to play an important role in the pathogenesis of hypertension. The authors conducted the systematic review and meta-analysis to clarify the relationship between gut microbiota and hypertension through conducting an electronic search in six databases. Our meta-analysis included 19 studies and the results showed that compared with healthy controls, Shannon significantly decreased in hypertension [SMD = -0.13, 95%CI (-0.22, -0.04), p = .007]; however, Simpson [SMD = -0.01, 95%CI (-0.14, 0.12), p = .87], ACE [SMD = 0.18, 95%CI (-0.06, 0.43), p = .14], and Chao1 [SMD = 0.11, 95%CI (-0.01, 0.23), p = .08] did not differ significantly between hypertension and healthy controls. The F/B ratio significantly increased in hypertension [SMD = 0.84, 95%CI (0.10, 1.58), p = .03]. In addition, Shannon index was negatively correlated with hypertension [r = -0.12, 95%CI (-0.19, -0.05)], but had no significant correlation with SBP [r = 0.10, 95%CI (-0.19, 0.37)] and DBP [r = -0.39, 95%CI (-0.73, 0.12)]. At the phylum level, the relative abundance of Firmicutes [SMD = -0.01, 95%CI (-0.37, 0.34), p = .94], Bacteroidetes [SMD = -0.15, 95%CI (-0.44, 0.14), p = .30], Proteobacteria [SMD = 0.25, 95%CI (-0.01, 0.51), p = .06], and Actinobacteria [SMD = 0.21, 95%CI (-0.11, 0.53), p = .21] did not differ significantly between hypertension and healthy controls. At the genus level, compared with healthy controls, the relative abundance of Faecalibacterium decreased significantly [SMD = -0.16, 95%CI (-0.28, -0.04), p = .01], while the Streptococcus [SMD = 0.20, 95%CI (0.08, 0.32), p = .001] and Enterococcus [SMD = 0.20, 95%CI (0.08, 0.33), p = .002] significantly increased in hypertension. Available evidence suggests that hypertensive patients may have an imbalance of gut microbiota. However, it still needs further validation by large sample size studies of high quality.

Advances in transposable elements: from mechanisms to applications in mammalian genomics

It has been 70 years since Barbara McClintock discovered transposable elements (TE), and the mechanistic studies and functional applications of transposable elements have been at the forefront of life science research. As an essential part of the genome, TEs have been discovered in most species of prokaryotes and eukaryotes, and the relative proportion of the total genetic sequence they comprise gradually increases with the expansion of the genome. In humans, TEs account for about 40% of the genome and are deeply involved in gene regulation, chromosome structure maintenance, inflammatory response, and the etiology of genetic and non-genetic diseases. In-depth functional studies of TEs in mammalian cells and the human body have led to a greater understanding of these fundamental biological processes. At the same time, as a potent mutagen and efficient genome editing tool, TEs have been transformed into biological tools critical for developing new techniques. By controlling the random insertion of TEs into the genome to change the phenotype in cells and model organisms, critical proteins of many diseases have been systematically identified. Exploiting the TE's highly efficient in vitro insertion activity has driven the development of cutting-edge sequencing technologies. Recently, a new technology combining CRISPR with TEs was reported, which provides a novel targeted insertion system to both academia and industry. We suggest that interrogating biological processes that generally depend on the actions of TEs with TEs-derived genetic tools is a very efficient strategy. For example, excessive activation of TEs is an essential factor in the occurrence of cancer in humans. As potent mutagens, TEs have also been used to unravel the key regulatory elements and mechanisms of carcinogenesis. Through this review, we aim to effectively combine the traditional views of TEs with recent research progress, systematically link the mechanistic discoveries of TEs with the technological developments of TE-based tools, and provide a comprehensive approach and understanding for researchers in different fields.

Five new F10 variants in hereditary factor x deficiency detected by high-throughput sequencing

INTRODUCTION

Factor X deficiency is a rare inherited bleeding disorder. To date, 181 variants are reported in the recently updated F10-gene variant database.

AIM

This study aimed to describe new F10 variants.

METHOD

The F10 gene was analysed in 16 consecutive families with FX deficiency by a targeted high-throughput sequencing approach, including F10, F9, F8 genes, and 78 genes dedicated to haematological malignancies.

RESULTS

We identified 19 variants (17 missense, one nonsense and one frameshift) and two copy number variations. Two patients presenting a combined FVII-FX deficiency showed a loss of one F10 gene copy (del13q34) associated with a missense variant on the remaining allele, leading to a FX:C significantly lower than the FVII:C level and explaining their unusual bleeding history. We reported five novel variants. Three missense variants (p.Glu22Val affecting the signal peptide cleavage site, p.Cys342Tyr removing the disulphide bond between the FX heavy and light chains, and p.Val385Met located in FX peptidase S1 domain) were detected at compound heterozygosis status in three patients with severe bleeding symptoms and FX:C level below 10 IU/dL. Two truncating variants p.Tyr279* and p.Thr434Aspfs*13 leading to an altered FX protein were found at heterozygous state in two patients with mild bleeding history.

CONCLUSION

This study showed the feasibility and the interest of high-throughput sequencing approach for rare bleeding disorders, enabling the report of F10 gene screening in a 3-weeks delay, suitable for clinical use. The description of five new variants may contribute to a better understanding of the phenotype-genotype correlation in FX deficiency.

Process of Forensic Medicine in DNA Identification of Aged Human Remains

Skeleton and teeth are important biological samples. Due to their special structure and strong ability to resist degradation, they are ideal biological materials to retain DNA under natural condition. In many cases, such as historical figure identification, aged skeleton and teeth are usually the only biological samples. However, their DNA is in a state of trace, damage and degradation to different degrees, which requires special experimental treatment to achieve identification. This paper reviews the sample selection, DNA extraction, DNA enrichment and analysis approaches based on relevant research reports in recent years, aiming to promote the further development and improvement of the aged skeleton and teeth identification system.

[Effects of Different Biochar and Effective Microorganism Agent Improvement Approaches on the Nutrient Release Characteristics and Potential of Compost]

The nutrient release characteristics of four types of composts, pure municipal sewage sludge compost, corn straw biochar (CSB) improved compost, effective microorganism agent (EM) improved compost, and CSB+EM improved compost, in coastal wetland soil were examined through a soil incubation experiment. The effects of different composts on the spectral characteristics of soil dissolved organic matter (DOM) and microbial community were also investigated. The results demonstrated that the compost additions could significantly reduce soil pH, while increasing soil electrical conductivity and contents of plant available nutrients (e.g., dissolved organic carbon, NH4+-N, NO3--N, available phosphorus, and available potassium). By comparing the nutrient release potential among the improved composts, the CSB+EM-improved compost (CSB+EM-C) evidently had the highest nutrient release potential. Furthermore, the DOM in CSB+EM-C amended soil exhibited a higher humification degree than that of the other composts. The high-throughput sequencing results indicated that the compost additions increased the relative abundances of dominant bacteria at the phylum level, such as the Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria. CSB+EM-C exhibited a greater potential to improve the relative abundance of these dominant bacteria phyla than other improved composts. Overall, among all the improvement approaches, the combined use of CSB and EM agent was the optimal composting strategy owing to its highest potentials of nutrient supply and soil quality improvement. The present findings can provide a solid scientific theoretical basis for establishing an effective technology strategy involving the combination of municipal sewage sludge utilization and degraded coastal wetland soil remediation.

High throughput miRNA sequencing and bioinformatics analysis identify the mesenchymal cell proliferation and apoptosis related miRNAs during fetal mice palate development

BACKGROUND

Palatogenesis requires a precise spatiotemporal regulation of gene expression. Recent studies indicate that microRNAs (miRNAs) are key factors in normal palatogenesis. The present study aimed to explain the regulatory mechanisms of miRNAs during palate development.

METHODS

Pregnant ICR mice were choose at embryonic day 10.5 (E10.5). Hemotoxylin and eosin (H&E) staining was used to observe the morphological changes during the development of palatal process at embryonic day (E)13.5, E14.0, E14.5, E15.0 and E15.5. The fetal palatal tissues were collected at E13.5, E14.0, E14.5 and E15.0 to explore miRNA expression and function by high throughput sequencing and bioinformatic analysis. Mfuzz cluster analysis was used to look for miRNAs related to the fetal mice palate formation. The target genes of miRNAs were predicted by miRWalk. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed base on target genes. The mesenchymal cell proliferation and apoptosis related miRNAs-genes networks were predicted and constructed using miRWalk and Cytoscape software. The expression of mesenchymal cell proliferation and apoptosis related miRNAs at the E13.5, E14.0, E14.5, and E15.0 was detected by a quantitative real-time PCR (RT-qPCR) assay.

RESULTS

H&E staining found that the palatal process grows vertically along the sides of the tongue at E13.5, the position of the tongue begins to descend and the bilateral palatal processes rise above the tongue at E14.0, the palatal process grows horizontally at E14.5, there is palatal contact fusion at E15.0, and the palatal suture disappeared at E15.5. Nine clusters of miRNA expression changes were identified in the fetal mice palate formation progression, including two reducing trends, two rising trends and five disordered trends. Next, the heatmap showed the miRNA expression from Clusters 4, 6, 9, 12 in the E13.5, E14.0, E14.5 and E15.0 groups. GO functional and KEGG pathway enrichment analysis found target genes of miRNAs in clusters involved in regulation of mesenchymal phenotype and the mitogen-activated protein kinase (MAPK) signaling pathway. Next, mesenchymal phenotype related miRNA-genes networks were constructed. The heatmap showing that the mesenchymal phenotype related miRNA expression of Clusters 4, 6, 9 and 12 at E13.5, E14.0, E14.5 and E15.0. Furthermore, the mesenchymal cell proliferation and apoptosis related miRNA-gene networks were identified in Clusters 6 and 12, including mmu-miR-504-3p-Hnf1b, etc. The expression level of mesenchymal cell proliferation and apoptosis related miRNAs at the E13.5, E14.0, E14.5, and E15.0 was verified by a RT-qPCR assay.

CONCLUSIONS

For the first time, we identified that clear dynamic miRNA expression during palate development. Furthermore, we demonstrated that mesenchymal cell proliferation and apoptosis related miRNAs, genes and the MAPK signaling pathway are important during fetal mice palate development.

Going beyond consensus genome sequences: An innovative SNP-based methodology reconstructs different Ugandan cassava brown streak virus haplotypes at a nationwide scale in Rwanda

Cassava Brown Streak Disease (CBSD), which is caused by cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV), represents one of the most devastating threats to cassava production in Africa, including in Rwanda where a dramatic epidemic in 2014 dropped cassava yield from 3.3 million to 900,000 tonnes (1). Studying viral genetic diversity at the genome level is essential in disease management, as it can provide valuable information on the origin and dynamics of epidemic events. To fill the current lack of genome-based diversity studies of UCBSV, we performed a nationwide survey of cassava ipomovirus genomic sequences in Rwanda by high-throughput sequencing (HTS) of pools of plants sampled from 130 cassava fields in thirteen cassava-producing districts, spanning seven agro-ecological zones with contrasting climatic conditions and different cassava cultivars. HTS allowed the assembly of a nearly complete consensus genome of UCBSV in twelve districts. The phylogenetic analysis revealed high homology between UCBSV genome sequences, with a maximum of 0.8 per cent divergence between genomes at the nucleotide level. An in-depth investigation based on Single Nucleotide Polymorphisms (SNPs) was conducted to explore the genome diversity beyond the consensus sequences. First, to ensure the validity of the result, a panel of SNPs was confirmed by independent reverse transcription polymerase chain reaction (RT-PCR) and Sanger sequencing. Furthermore, the combination of fixation index (FST) calculation and Principal Component Analysis (PCA) based on SNP patterns identified three different UCBSV haplotypes geographically clustered. The haplotype 2 (H2) was restricted to the central regions, where the NAROCAS 1 cultivar is predominantly farmed. RT-PCR and Sanger sequencing of individual NAROCAS1 plants confirmed their association with H2. Haplotype 1 was widely spread, with a 100 per cent occurrence in the Eastern region, while Haplotype 3 was only found in the Western region. These haplotypes' associations with specific cultivars or regions would need further confirmation. Our results prove that a much more complex picture of genetic diversity can be deciphered beyond the consensus sequences, with practical implications on virus epidemiology, evolution, and disease management. Our methodology proposes a high-resolution analysis of genome diversity beyond the consensus between and within samples. It can be used at various scales, from individual plants to pooled samples of virus-infected plants. Our findings also showed how subtle genetic differences could be informative on the potential impact of agricultural practices, as the presence and frequency of a virus haplotype could be correlated with the dissemination and adoption of improved cultivars.

Molecular techniques drive cutting edge advancements in management of cutaneous T cell lymphoma

Cutaneous 5T cell lymphoma (CTCL), characterized by malignant T cells infiltrating the skin with potential for dissemination, remains a challenging disease to diagnose and treat due to disease heterogeneity, treatment resistance, and lack of effective and standardized diagnostic and prognostic clinical tools. Currently, diagnosis of CTCL practically relies on clinical presentation, histopathology, and immunohistochemistry. These methods are collectively fraught with limitations in sensitivity and specificity. Fortunately, recent advances in flow cytometry, polymerase chain reaction, high throughput sequencing, and other molecular techniques have shown promise in improving diagnosis and treatment of CTCL. Examples of these advances include T cell receptor clonotyping via sequencing to detect CTCL earlier in the disease course and single-cell RNA sequencing to identify gene expression patterns that commonly drive CTCL pathogenesis. Experience with these techniques has afforded novel insights which may translate into enhanced diagnostic and therapeutic approaches for CTCL.

TRcaller: a novel tool for precise and ultrafast tandem repeat variant genotyping in massively parallel sequencing reads

Calling tandem repeat (TR) variants from DNA sequences is of both theoretical and practical significance. Some bioinformatics tools have been developed for detecting or genotyping TRs. However, little study has been done to genotyping TR alleles from long-read sequencing data, and the accuracy of genotyping TR alleles from next-generation sequencing data still needs to be improved. Herein, a novel algorithm is described to retrieve TR regions from sequence alignment, and a software program TRcaller has been developed and integrated into a web portal to call TR alleles from both short- and long-read sequences, both whole genome and targeted sequences generated from multiple sequencing platforms. All TR alleles are genotyped as haplotypes and the robust alleles will be reported, even multiple alleles in a DNA mixture. TRcaller could provide substantially higher accuracy (>99% in 289 human individuals) in detecting TR alleles with magnitudes faster (e.g., ∼2 s for 300x human sequence data) than the mainstream software tools. The web portal preselected 119 TR loci from forensics, genealogy, and disease related TR loci. TRcaller is validated to be scalable in various applications, such as DNA forensics and disease diagnosis, which can be expanded into other fields like breeding programs. Availability: TRcaller is available at https://www.trcaller.com/SignIn.aspx.

E-prescription and invisible work in genomics in France

This article aims to analyze the transformations in medical prescription work and infrastructures brought by digitalization. Our fieldwork takes place in the context of precision medicine development based on genomics High Throughput Sequencing (HTS) in France, through the Plan France Médecine Génomique (PFMG 2025). The Plan aims at industrializing the production of genomic testing in clinical context at a national scale, particularly in oncology. To ensure the intensified flow of information between hospitals and HTS platforms required, a centralized process has been organized around two sequencing platforms and the introduction of a new e-prescription software (E-PRES). We start by analyzing how the e-prescription software changes the practices of health professionals by imposing new technological and professional standards. We show that, more than a mere prescription tool, this software is also a monitoring tool for the platforms and prescribers' work, and a support tool for the logistical and work organization. Secondly, we question the division of labor among the different professionals involved in the organizational or technical tasks required. We show that the feasibility of this new form of digitalized prescription relies on an important datawork performed by "small hands" to select, translate and process a vast amount of heterogeneous data.

Conjunctive Analyses of BSA-Seq and BSR-Seq to Identify Candidate Genes Controlling the Black Lemma and Pericarp Trait in Barley

Black barley seeds are a health-beneficial diet resource because of their special chemical composition and antioxidant properties. The black lemma and pericarp (BLP) locus was mapped in a genetic interval of 0.807 Mb on chromosome 1H, but its genetic basis remains unknown. In this study, targeted metabolomics and conjunctive analyses of BSA-seq and BSR-seq were used to identify candidate genes of BLP and the precursors of black pigments. The results revealed that five candidate genes, purple acid phosphatase, 3-ketoacyl-CoA synthase 11, coiled-coil domain-containing protein 167, subtilisin-like protease, and caffeic acid-O-methyltransferase, of the BLP locus were identified in the 10.12 Mb location region on the 1H chromosome after differential expression analysis, and 17 differential metabolites, including the precursor and repeating unit of allomelanin, were accumulated in the late mike stage of black barley. Phenol nitrogen-free precursors such as catechol (protocatechuic aldehyde) or catecholic acids (caffeic, protocatechuic, and gallic acids) may promote black pigmentation. BLP can manipulate the accumulation of benzoic acid derivatives (salicylic acid, 2,4-dihydroxybenzoic acid, gallic acid, gentisic acid, protocatechuic acid, syringic acid, vanillic acid, protocatechuic aldehyde, and syringaldehyde) through the shikimate/chorismite pathway other than the phenylalanine pathway and alter the metabolism of the phenylpropanoid-monolignol branch. Collectively, it is reasonable to infer that black pigmentation in barley is due to allomelanin biosynthesis in the lemma and pericarp, and BLP regulates melanogenesis by manipulating the biosynthesis of its precursors.

First Identification of Fig Virus A and Fig Virus B in Ficus carica in Italy

Extracts of double-stranded RNA from three fig cultivars from an Apulian (Southern Italy) germplasm collection were used for high-throughput sequencing and revealed the presence of two distinct, recently described closteroviruses. Sequences obtained from these Apulian isolates belong to fig virus A and fig virus B and cover 38 and 25% of their RNA genome, respectively. Primer sets designed on selected contigs confirmed the presence of each virus in infected plants. A close phylogenetic relationship, investigated in a fragment of HSP70h protein, occurs among these isolates and the reference genomes. A nucleotide divergence (ranging from 10 to 30% along the different genes) was observed among our isolates and the reference genomes. This is the first finding of these virus species in autochthonous fig accessions in Europe.

Distinct Bacterial Communities in São Jorge Cheese with Protected Designation of Origin (PDO)

São Jorge cheese is an iconic product of the Azores, produced from raw cow's milk and natural whey starter (NWS). Although it is produced according to Protected Designation of Origin (PDO) specifications, the granting of the PDO label depends crucially on sensory evaluation by trained tasters. The aim of this work was to characterize the bacterial diversity of this cheese using next-generation sequencing (NGS) and to identify the specific microbiota that contributes most to its uniqueness as a PDO by distinguishing the bacterial communities of PDO and non-PDO cheeses. The NWS and curd microbiota was dominated by Streptococcus and Lactococcus, whereas Lactobacillus and Leuconostoc were also present in the core microbiota of the cheese along with these genera. Significant differences (p < 0.05) in bacterial community composition were found between PDO cheese and non-certified cheese; Leuconostoc was found to play the chief role in this regard. Certified cheeses were richer in Leuconostoc, Lactobacillus and Enterococcus, but had fewer Streptococcus (p < 0.05). A negative correlation was found between contaminating bacteria, e.g., Staphylococcus and Acinetobacter, and the development of PDO-associated bacteria such as Leuconostoc, Lactobacillus and Enterococcus. A reduction in contaminating bacteria was found to be crucial for the development of a bacterial community rich in Leuconostoc and Lactobacillus, thus justifying the PDO seal of quality. This study has helped to clearly distinguish between cheeses with and without PDO based on the composition of the bacterial community. The characterization of the NWS and the cheese microbiota can contribute to a better understanding of the microbial dynamics of this traditional PDO cheese and can help producers interested in maintaining the identity and quality of São Jorge PDO cheese.

Atypical familial diabetes associated with a novel NEUROD1 nonsense variant

OBJECTIVES

We aimed to identify the origin of atypical diabetes in a family with four generations of diabetes from South Asia. The family members showed different clinical phenotypes. Members of generation one to three were presumed to have type 2 diabetes and generation four to have type 1 diabetes.

CASE PRESENTATION

We performed a genetic analysis of the family using targeted high throughput sequencing.

CONCLUSIONS

We identified a novel nonsense variant in the neurogenic differentiation 1 (NEUROD1) gene, co-segregating with diabetes. The variant was located in the DNA-binding domain, altering a protein residue that was very well conserved among different species.

Exploration of the core gene signatures and mechanisms between NAFLD and sarcopenia through transcriptomic level

INTRODUCTION

The increased prevalence of non-alcoholic fatty liver disease (NAFLD) and sarcopenia among the elderly are facing a significant challenge to the world's health systems. Our study aims to identify the coexpressed genes in NAFLD and sarcopenia patients.

METHODS

We downloaded the transcriptome data of NAFLD tissue from patients, as well as muscle tissues from sarcopenia patients, from the GEO database in order to investigate the shared transcriptional regulation mechanisms between these two diseases. Then, focusing on the genes that were frequently expressed in these diseases, together with GSVA and WGCNA, we utilized a range of analysis methods to identify the main co-expressed genes in both diseases by taking intersections. We investigated these changes after learning that they mostly affected lipid metabolism and oxidative stress injury pathways.

RESULTS

By analyzing these genes and their interactions with transcription factors and proteins, we were able to identify 8 genes that share common patterns. From these 8 genes, we were possible to forecast potential future medicines. Our research raises the possibility of NAFLD and sarcopenia transcriptome regulatory pathways in aging populations.

DISCUSSION

In conclusion, a complete transcription pattern mapping was carried out in order to identify the core genes, underlying biological mechanisms, and possible therapeutic targets that regulate aging in NAFLD and sarcopenia patients. It provides novel insights and proof in favor of decreasing the increased prevalence of sarcopenia in the elderly caused by NAFLD.

Changes in bacterial diversity and community structure in drinking water distribution system revealed by high throughput sequencing

For microbiological management of water quality, it is important to identify bacteria and to understand the community structure. To analyze the community structure during water purification and distribution, we selected a distribution system in which water from other water treatment facilities was not mixed with the target water. Changes in the bacterial community structure during treatment and distribution processes in a slow filtration water treatment facility were analyzed using 16S rRNA gene amplicon sequencing with a portable sequencer MinION. The microbial diversity was reduced by chlorination. The genus level diversity increased during distribution and this diversity was maintained through to the terminal tap water. Yersinia and Aeromonas were dominant in the intake water, and Legionella was dominant in the slow sand filtered water. Chlorination greatly reduced the relative abundance of Yersinia, Aeromonas, and Legionella, and these bacteria were not detected in the terminal tap water. Sphingomonas, Starkeya and Methylobacterium became dominant in the water after chlorination. These bacteria could be used as important indicator bacteria to provide useful information for microbiological control in drinking water distribution systems.

Corrigendum: Artemisinin relieves myocardial ischemia-reperfusion injury via modulating miR-29b-3p and hemicentin 1

[This corrects the article DOI: 10.3389/fphar.2022.918966.].

Comparative Analysis of Tomato Brown Rugose Fruit Virus Isolates Shows Limited Genetic Diversity

Tomato is an important vegetable in the United States and around the world. Recently, tomato brown rugose fruit virus (ToBRFV), an emerging tobamovirus, has impacted tomato crops worldwide and can result in fruit loss. ToBRFV causes severe symptoms, such as mosaic, puckering, and necrotic lesions on leaves; other symptoms include brown rugose and marbling on fruits. More importantly, ToBRFV can overcome resistance in tomato cultivars carrying the Tm-2² locus. In this study, we recovered ToBRFV sequences from tomato seeds, leaves, and fruits from the U.S., Mexico, and Peru. Samples were pre-screened using a real-time RT-PCR assay prior to high-throughput sequencing. Virus draft genomes from 22 samples were assembled and analyzed against more than 120 publicly available genomes. Overall, most sequenced isolates were similar to each other and did not form a distinct population. Phylogenetic analysis revealed three clades within the ToBRFV population. Most of the isolates (95%) clustered in clade 3. Genetic analysis revealed differentiation between the three clades indicating minor divergence occurring. Overall, pairwise identity showed limited genetic diversity among the isolates in this study with worldwide isolates, with a pairwise identity ranging from 99.36% and 99.97%. The overall population is undergoing high gene flow and population expansion with strong negative selection pressure at all ToBRFV genes. Based on the results of this study, it is likely that the limited ToBRFV diversity is associated with the rapid movement and eradication of ToBRFV-infected material between countries.

Effects of Wheat Tempering with Slightly Acidic Electrolyzed Water on the Microbiota and Flour Characteristics

Slightly acidic electrolyzed water (SAEW) was prepared and used as wheat tempering water. This study explored the impacts of tempering with SAEW on microbial load and diversity and quality properties of wheat flour. As SAEW volume ratio increased, the residual level of total plate counts (TPC) and mould/yeast counts (MYC) decreased dramatically (p < 0.05). Based on genomics analysis, bacterial 16S rRNA gene and fungal ITS1 gene region were performed to characterize the changes in microbial communities’ composition and diversity in response to SAEW treatment. SAEW optimal volume ratio (6.5:10, v/v) of SAEW with distilled water influenced wheat microbiome composition, with a higher microbial diversity and abundance discovered on the control grains. Bacteroidetes of predominant bacterial phylum and Ascomycota of the most abundant fungal phylum were reduced after SAEW optimal volume ratio tempering. The flour yield is higher and ash content is lower than the control samples. Falling number and “b*” in terms of colour markedly increased. DSC (Differential Scanning Calorimetry) test showed that To (onset temperature), Tp (peak temperature), and Tc (conclusion temperature) were significantly decreased in thermal characteristics of flour. Gluten content, protein content, ΔH and pasting properties tests showed no significant change. It can be concluded that SAEW should be applied on wheat tempering for producing clean wheat flour. ANOVA and Tukey’s honestly significant difference (HSD) test were used for the analysis of variance and differences between the experimental and control groups, with p < 0.05.

Diversity, structure, and distribution of bacterioplankton and diazotroph communities in the Bay of Bengal during the winter monsoon

The Bay of Bengal (BoB) is conventionally believed to be a low productive, oligotrophic marine ecosystem, where the diazotroph communities presumed to play a vital role in adding "new" nitrogen through the nitrogen fixation process. However, the diazotroph communities in the oceanic region of the BoB are still poorly understood though it represents most of the seawater volume. The present study investigated a detailed account of the bacterioplankton community structure and distribution in the oceanic BoB during the winter monsoon using high throughput sequencing targeting the 16S rRNA and nifH genes. Our study observed diverse groups of bacterioplankton communities in the BoB including both cyanobacterial and non-cyanobacterial phylotypes. Cyanobacteria (Prochlorococcus spp. and Synechococcus spp.) and Proteobacteria (mainly α-, γ-, and δ-Proteobacteria) were the most abundant groups within the bacterial communities, possessing differential vertical distribution patterns. Cyanobacteria were more abundant in the surface waters, whereas Proteobacteria dominated the deeper layers (75 m). However, within the diazotroph communities, Proteobacteria (mainly γ-Proteobacteria) were the most dominant groups than Cyanobacteria. Function prediction based on PICRUSt revealed that nitrogen fixation might more active to add fixed nitrogen in the surface waters, while nitrogen removal pathways (denitrification and anammox) might stronger in deeper layers. Canonical correspondence analysis (CCA) indicated that temperature, salinity, and silicate were major environmental factors driving the distribution of bacterial communities. Additionally, phosphate was also an important factor in regulating the diazotroph communities in the surface water. Overall, this study provided detailed information on bacterial communities and their vital role in the nitrogen cycles in oligotrophic ecosystems.

Research on the Gut Microbiota of Hainan Black Goat

The intestine of animals is a complex micro-ecosystem containing a large number of microbiomes, which is essential for the host's health development. The Hainan black goat with good resistance and adaptability is a unique species in Hainan, China. These unique physiological characteristics are inseparable from their intestinal microbiota. In this study, high-throughput sequencing was used to investigate bacterial communities in different segments of the intestinal tract of Hainan black goat. The results showed that the indices of Chao1 and ACE in the cecum and colon were significantly greater than those in the ileum (p = 0.007, 0.018). According to PCoA, the intestinal flora composition of the cecum and colon is almost equivalent. In contexts of the phylum, Firmicutes, Bacteroidota, and Pseudomonadota were the dominant phyla in the gut of the Hainan black goat. While in context of the genus, the dominant groups in the gut of black goats mainly include Ruminococcaceae_UCG-005, Bacteroides, Paeniclostridium, Christensenellaceae_R-7_group, Rikenellaceae_RC9_gut_group, and Eubacterium coprostanoligenes _group, Prevotella_1, they have different proportions in different intestinal segments. The gut microbiota of Hainan black goat is mainly Firmicutes, Bacteroidota, and Pseudomonadota. Influenced by the intestinal location where they colonize, the large intestine has a more complex intestinal flora than the small intestine. In contrast, there are only minor differences between the caecum and the colon in the large intestine.

Phyllosphere microbial community of cigar tobacco and its corresponding metabolites

Cigar is made of a typical fermented tobacco where the microbiota inhabits within an alkaline environment. Our current understanding on cigar fermentation is far from thorough. This work employed both high-throughput sequencing and chromatography-mass spectrometric technologies to provide new scientific reference for this specific fermented system. Typical cigar samples from different regions (the Caribbeans, South America, East Asia, and Southeast Asia) were investigated. The results show that Firmicutes, Actinobacteria, Proteobacteria, Ascomycota, and Basidiomycota were the predominant phyla in the cigar samples. Rather than the fungal community, it was the bacterial community structures that played vital roles to differentiate the cigar from different regions: Staphylococcus was the dominant genus in the Americas; Bacillus was the dominant genus in Southeast Asia; while in East Asia, there was no dominant genus. Such differences in community structure then affected the microflora metabolism. The correlation between microbiota and metabolites revealed that Aspergillaceae, Cercospora, and Staphylococcus were significantly correlated with sclareolide; Bacillus were positively associated with isophorone. Alcaligenaceae was significantly and positively correlated with L-nicotine and hexadecanoic acid, methyl ester. GRAPHICAL ABSTRACT.

Examining the effects of Salmonella phage on the caecal microbiota and metabolome features in Salmonella-free broilers

Bacteriophages selectively infect and kill their target bacterial host, being a promising approach to controlling zoonotic bacteria in poultry production. To ensure confidence in its use, fundamental questions of safety and toxicity monitoring of phage therapy should be raised. Due to its high specificity, a minimal impact on the gut ecology is expected; however, more in-depth research into key parameters that influence the success of phage interventions has been needed to reach a consensus on the impact of bacteriophage therapy in the gut. In this context, this study aimed to investigate the interaction of phages with animals; more specifically, we compared the caecum microbiome and metabolome after a Salmonella phage challenge in Salmonella-free broilers, evaluating the role of the phage administration route. To this end, we employed 45 caecum content samples from a previous study where Salmonella phages were administered via drinking water or feed for 24 h from 4, 5 to 6-weeks-old broilers. High-throughput 16S rRNA gene sequencing showed a high level of similarity (beta diversity) but revealed a significant change in alpha diversity between broilers with Salmonella-phage administered in the drinking water and control. Our results showed that the phages affected only a few genera of the microbiota's structure, regardless of the administration route. Among these, we found a significant increase in Streptococcus and Sellimonas in the drinking water and Lactobacillus, Anaeroplasma and Clostridia_vadinBB60_group in the feed. Nevertheless, the LC-HRMS-based metabolomics analyses revealed that despite few genera were significantly affected, a substantial number of metabolites, especially in the phage administered in the drinking water were significantly altered (64 and 14 in the drinking water and feed groups, respectively). Overall, our study shows that preventive therapy with bacteriophages minimally alters the caecal microbiota but significantly impacts their metabolites, regardless of the route of administration.

Mixed Infection of Blackcurrant with a Novel Cytorhabdovirus and Black Currant-Associated Nucleorhabdovirus

A virome screen was performed on a new breeding line, KB1, of blackcurrant. Rhabdovirus-like particles were observed by electron microscopy in ultrathin sections of flower stalks, and the complete genome sequence of a novel virus, provisionally named blackcurrant rhabdovirus 2 (BCRV2), was determined and verified using high-throughput sequencing. The genomic organization of BCRV2 was characteristic of cytorhabdoviruses (family Rhabdoviridae) and included seven genes: 3 ́- N-P´-P-P3-M-G-L -5 ́. BLASTP analysis revealed that the putative L protein had the highest amino acid sequence identity (75 %) with strawberry virus 2. BCRV2 was detected in Cryptomyzusgaleopsidis, but efficient transmission by this aphid was not confirmed. Of note, we observed coinfection of the KB1 line with blackcurrant-associated rhabdovirus (BCaRV) by RT-PCR. This is likely the first evidence of the presence of a cyto- and a nucleorhabdovirus in a single host.

Biological and Genetic Characterization of Physostegia Chlorotic Mottle Virus in Europe Based on Host Range, Location, and Time

Application of high throughput sequencing (HTS) technologies enabled the first identification of Physostegia chlorotic mottle virus (PhCMoV) in 2018 in Austria. Subsequently, PhCMoV was detected in Germany and Serbia on tomatoes showing severe fruit mottling and ripening anomalies. We report here how prepublication data-sharing resulted in an international collaboration across eight laboratories in five countries, enabling an in-depth characterization of PhCMoV. The independent studies converged toward its recent identification in eight additional European countries and confirmed its presence in samples collected 20 years ago (2002). The natural plant host range was expanded from two to nine species across seven families, and we confirmed the association of PhCMoV presence with severe fruit symptoms on economically important crops such as tomato, eggplant, and cucumber. Mechanical inoculations of selected isolates in the greenhouse established the causality of the symptoms on a new indexing host range. In addition, phylogenetic analysis showed a low genomic variation across the 29 near-complete genome sequences available. Furthermore, a strong selection pressure within a specific ecosystem was suggested by nearly identical sequences recovered from different host plants through time. Overall, this study describes the European distribution of PhCMoV on multiple plant hosts, including economically important crops on which the virus can cause severe fruit symptoms. This work demonstrates how to efficiently improve knowledge on an emergent pathogen by sharing HTS data and provides a solid knowledge foundation for further studies on plant rhabdoviruses.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Identification, Characterization, and Detection of a Novel Strawberry Cytorhabdovirus

In 2020, a novel agent was discovered in strawberry, a rhabdovirus closely related to lettuce necrotic yellows virus. The new virus, named strawberry virus 2 (StrV-2), was discovered in an accession of the Fragaria virus collection of the National Clonal Germplasm Repository (NCGR), and for this reason, it was studied in-depth. The complete StrV-2 genome was obtained and investigated in silico. Transmission was assessed using two aphid species whereas a multiplex RT-PCR test targeting plant and virus genes was developed and used to screen the NCGR Fragaria virus collection.

Librarian: A quality control tool to analyse sequencing library compositions

BACKGROUND

Robust analysis of DNA sequencing data needs to include a set of quality control steps to ensure that technical bias is kept to a minimum. A metric easily obtained is the frequency of each of the nucleobases for each position across all sequencing reads. Here, we explore the differences in nucleobase compositions of various library types produced by standard experimental methodologies.  Methods: We obtained the compositions of nearly 3000 publicly available datasets and subjected them to Uniform Manifold Approximation and Projection (UMAP) dimensionality reduction for a two-dimensional representation of their composition characteristics.   Results: We find that most library types result in a specific composition profile. We use this to give an estimate of how strongly the composition of a test library resembles the profiles of previously published libraries, and how likely the test sample is to be of a particular type. We introduce Librarian, a user-friendly web application and command line tool which enables checking base compositions of test libraries against known library types.   Conclusions: Library preparation methods strongly influence the per position nucleobase content. By comparing test libraries to a database of previously published library types we can make predictions regarding the library preparation method. Librarian is a user-friendly tool to access this information for quality assurance purposes as discrepancies can flag potential irregularities very early on.

Librarian: A quality control tool to analyse sequencing library compositions

Background

Robust analysis of DNA sequencing data needs to include a set of quality control steps to ensure that technical bias is kept to a minimum. A metric easily obtained is the frequency of each of the nucleobases for each position across all sequencing reads. Here, we explore the differences in nucleobase compositions of various library types produced by standard experimental methodologies.

Methods

We obtained the compositions of nearly 3000 publicly available datasets and subjected them to Uniform Manifold Approximation and Projection (UMAP) dimensionality reduction for a two-dimensional representation of their composition characteristics.

Results

We find that most library types result in a specific composition profile. We use this to give an estimate of how strongly the composition of a test library resembles the profiles of previously published libraries, and how likely the test sample is to be of a particular type. We introduce Librarian, a user-friendly web application and command line tool which enables checking base compositions of test libraries against known library types.

Conclusions

Library preparation methods strongly influence the per position nucleobase content. By comparing test libraries to a database of previously published library types we can make predictions regarding the library preparation method. Librarian is a user-friendly tool to access this information for quality assurance purposes as discrepancies can flag potential irregularities very early on.

Structural and Dynamic Analysis of Leaf-Associated Fungal Community of Walnut Leaves Infected by Leaf Spot Disease Based Illumina High-Throughput Sequencing Technology

Leaf-associated microbiota is vital in plant-environment interactions and is the basis for micro-ecological regulation. However, there are no studies on the direct differences in microbial community composition between disease-susceptible and healthy walnut leaves. This study collected five samples of healthy and infected leaves (all leaves with abnormal spots were considered diseased leaves) from May to October 2018. Differences in fungal diversity (Chao1 index, Shannon index, and Simpson index) and community structure were observed by sequencing and analyzing diseased and healthy leaf microbial communities by Illumina HiSeq sequencing technology. The main fungal phyla of walnut leaf-associated were Ascomycota, Basidiomycota, and Glomeromycota. Diversity indices (Shannon and Chao1 index values) of healthy leaves differed significantly in the late stages of disease onset. The results showed that the fungal species that differed considerably between the healthy and infected groups differed, and the fungal species that differed significantly between the healthy and infected groups changed with the development of the leaf disease. Critical control time points were determined by analyzing the population dynamics of pathogenic fungi. Leaf-associated microorganisms are abundant and diverse, and fungal identification and diversity studies are helpful for developing more appropriate walnut management strategies.

A multiview model for relative and absolute microbial abundances

The absolute abundance of bacterial taxa in human host-associated environments plays a critical role in reproductive and gastrointestinal health. However, obtaining the absolute abundance of many bacterial species is typically prohibitively expensive. In contrast, relative abundance data for many species are comparatively cheap and easy to collect (e.g., with universal primers for the 16S rRNA gene). In this paper, we propose a method to jointly model relative abundance data for many taxa and absolute abundance data for a subset of taxa. Our method provides point and interval estimates for the absolute abundance of all taxa. Crucially, our proposal accounts for differences in the efficiency of taxon detection in the relative and absolute abundance data. We show that modeling taxon-specific efficiencies substantially reduces the estimation error for absolute abundance, and controls the coverage of interval estimators. We demonstrate the performance of our proposed method via a simulation study, a study of the effect of HIV acquisition on microbial abundances, and a sensitivity study where we jackknife the taxa with observed absolute abundances.

Survey for Virus Diversity in Common Bean (Phaseolus vulgaris) Fields and the Detection of a Novel Strain of Cowpea polerovirus 1 in Zambia

The production of common bean (Phaseolus vulgaris L.) is adversely affected by virus-like diseases globally, but little is known about the occurrence, distribution, and diversity of common bean-infecting viruses in Zambia. Consequently, field surveys were conducted during the 2018 season in 128 fields across six provinces of Zambia and 640 common bean leaf tissue samples were collected with (n = 585) or without (n = 55) symptoms. The prevalence of symptomatic fields was 100%, but incidence of symptomatic plants ranged from 32 to 67.5%. Metagenomic analyses of nine composite samples and a single plant sample of interest revealed the occurrence of isolates of Bean common mosaic necrosis virus, Bean common mosaic virus, Cowpea aphid-borne mosaic virus, Peanut mottle virus, Southern bean mosaic virus (SBMV), Cucumber mosaic virus, Phaseolus vulgaris alphaendornavirus 1 (PvEV-1), PvEV-2, Ethiopian tobacco bushy top virus (ETBTV), and a novel strain of Cowpea polerovirus 1 (CPPV1-Pv) of 5,902 nt in length. While CPPV1-Pv was consistently detected in mixed infection with ETBTV and its satellite RNA molecule, based on results of mechanical transmission assays it does not appear to be involved in disease etiology, suggesting that its role may be limited to being a helper virus for the umbravirus. Screening of the survey samples by real-time PCR for the viruses detected by high-throughput sequencing revealed the prevalence of single (65.2% or 417/640) over mixed (1.9% or 12/640) infections in the samples. SBMV was the most frequently detected virus, occurring in ∼29.4% (188/640) of the samples and at a prevalence rate of 58.6% (75/128) across fields. The results showed that diverse virus species are present in Zambian common bean fields and the information will be useful for the management of common bean viral diseases.

Study on microbial community of "green-covering" Tuqu and the effect of fortified autochthonous Monascus purpureus on the flavor components of light-aroma-type Baijiu

"Green-covering" Tuqu (TQ), as one of Xiaoqu, is a special fermentative starter (also known as Jiuqu in Chinese) that originated in southern China and is characterized by a layer of green mold covering (Aspergillus clavatus) the surface and (sometimes) with a red heart. It plays a vital role in producing light-aroma-type Baijiu (LATB). However, to date, the microbiota that causes red heart of TQ remain largely unexplored, and it is still unclear how these microbiota influence on the quality of LATB. In this study, two types of TQ, one with a red heart (RH) and another with a non-red heart (NRH), were investigated by high throughput sequencing (HTS) and directional screening of culture-dependent methods. The obtained results revealed the differences in the microbial communities of different TQ and led to the isolation of two species of Monascus. Interestingly, the results of high performance liquid chromatography (HPLC) detection showed that citrinin was not detected, indicating that Monascus isolated from TQ was no safety risk, and the contents of gamma-aminobutyric acid in the fermented grains of RH were higher than that of NRH during the fermentation. Selecting the superior autochthonous Monascus (M1) isolated from the TQ to reinoculate into the TQ-making process, established a stable method for producing the experimental "red heart" Tuqu (ERH), which confirmed that the cause of "red heart" was the growth of Monascus strains. After the lab-scale production test, ERH increased ethyl ester production and reduced higher alcohols production. In addition, Monascus had an inhibitory effect on the growth of Saccharomyces and Aspergillus. This study provides the safe, health-beneficial, and superior fermentation strains and strategies for improving the quality of TQ and LATB.

Metataxonomic Mapping of the Microbial Diversity of Irish and Eastern Mediterranean Cheeses

The distinct sensorial characteristics of local cheeses influence consumer preferences, and make an essential contribution to the local economy. Microbial diversity in cheese is among the fundamental contributors to sensorial and qualitative characteristics. However, knowledge regarding the existence of microbial patterns associated with regional production practices in ripened cheeses remains limited. The present research was conducted to test the hypothesis that the background metagenome of cheeses could be used as a marker of their origin. We compared Irish versus Eastern Mediterranean cheeses-namely Greek and Cypriot-using High Throughput Sequencing (HTS). The study identified a significantly distinct separation among cheeses originating from the three different countries, in terms of the total microbial community composition. The use of machine learning and biomarkers discovery algorithms defined key microbes that differentiate each geographic region. Finally, the development of interaction networks revealed that the key species developed mostly negative interactions with the other members of the communities, highlighting their dominance in the community. The findings of the present research demonstrate that metagenome could indeed be used as a biological marker of the origin of mature cheeses, and could provide further insight into the dynamics of microbial community composition in ripened cheeses.

Cap-snatching as a possible contributor to photosynthesis shut-off

Cap-snatching is a mechanism applied by segmented, negative strand (-) RNA viruses (NSVs) to initiate genome transcription. So far, the cap donor source of cytoplasmic-replicating NSVs has remained elusive. Recently, studies pointed to processing body (P body, PB) as the potential source for providing capped RNAs but conclusive evidence is still lacking. To attempt identifying these sources, here the 5' non-viral leader sequences of Tomato spotted wilt virus (TSWV) N mRNAs were analysed by high-throughput sequencing (HTS) from plants subjected to normal and heat-stress conditions, and subsequently mapped on host donor transcripts. The majority of non-viral heterogenous, host-derived leader sequences ranged in size between ~10-20 nt and contained A or AG residues at the cleavage site and the presence of certain sequence motifs. Mapping the capped-leader sequences to the 5' UTR region of genes encoded by the Nicotiana tabacum genome, identified 348 donor genes and which were specifically enriched in cellular photosynthesis pathway. Nineteen of those were clearly expressed differentially at normal condition versus heat-stress conditions. Although the results did not point towards snatching of capped-RNA leader sequences from certain cytoplasmic RNA granules in particular, they indicated photosynthesis downregulation (and development of disease symptoms) partially result from cap-snatching.

Revealing the Microbiome of Four Different Thermal Springs in Turkey with Environmental DNA Metabarcoding

Simple Summary

The physicochemical conditions of thermal springs are one of the most significant barriers for detecting microbial life. According to the findings of various studies, high-throughput DNA sequencing technology can be utilized to perform more precise and detailed microbiome assessments. The main goal of this paper was to determine the microbiome in a thermal spring by metabarcoding environmental DNA obtained from four different sources and revealing how temperature and chemical composition affect the microbiome. This research also aimed to gather information that will aid in determining the best gene region and bioinformatic pipeline. The findings revealed a link between four different thermal springs’ physicochemical parameters and microbial composition and we found various manipulable steps in this study. This research is also first comprehensive thermal spring metabarcoding study conducted in Turkey.

Abstract

One of the most significant challenges for detecting microbial life in thermal springs by conventional techniques such as culturing is these places’ physicochemical (temperature, heavy metal content, pH, etc.) conditions. Data from several studies suggest that high-throughput DNA sequencing technologies can be used to perform more accurate and detailed microbiome analyses. The primary aim of this paper was to determine the microbiome in the thermal source by metabarcoding environmental DNA isolated from four different sources and reveal the reflection of differences caused by temperature and chemical content on the microbiome. DNA was extracted from water filtered with enclosed filters and using the Illumina high-throughput sequencing platform, V3 and V4 regions of the 16S rRNA gene were sequenced. The results showed a correlation between physicochemical conditions and microorganism composition of four different thermal springs. Springs with extremely high temperature (89–90 °C) were dominated by hyperthermophiles such as Hydrogenobacter and Thermus, while a spring with a high temperature (52 °C) was dominated by thermophiles such as Thermoanaerobaculum and Desulfurispora, and a spring with a low temperature (26 °C) and high salinity was dominated by halophiles and sulfur-oxidizers such as Hydrogenovibrio and Sulfirimonas. With this research, we observed many manipulable steps according to the work of interest. This study sought to obtain data that will help decide the right gene region and choose the optimal bioinformatic pipeline.

Genome-Wide Association Mapping of Hulless Barely Phenotypes in Drought Environment

Drought stress is one of the main factors restricting hulless barley (Hordeum vulgare L. var. nudum Hook. f.) yield. Genome-wide association study was performed using 269 lines of hulless barley to identify single-nucleotide polymorphism (SNP) markers associated with drought-resistance traits. The plants were cultured under either normal or drought conditions, and various quantitative traits including shoot fresh weight, shoot dry weight, root fresh weight, root dry weight, leaf fresh weight, leaf saturated fresh weight, leaf dry weight, ratio of root and shoot fresh weight, ratio of root and shoot dry weight, shoot water loss rate, root water loss rate, leaf water content and leaf relative water content, and field phenotypes including main spike length, grain number per plant, grain weight per plant, thousand grain weight (TGW), main spike number, plant height, and effective spike number of plants were collected. After genotyping the plants, a total of 8,936,130 highly consistent population SNP markers were obtained with integrity > 0.5 and minor allele frequency > 0.05. Eight candidate genes potentially contributed to the hulless barley drought resistance were obtained at loci near significant SNPs. For example, EMB506, DCR, and APD2 genes for effective spike number of plants, ABCG11 gene for main spike number (MEN), CLPR2 gene for main spike length, YIP4B gene for root and shoot dry weight (RSWD), and GLYK and BTS genes for TGW. The SNPs and candidate genes identified in this study will be useful in hulless barley breeding under drought resistance.

Seagrass Colonization Alters Diversity, Abundance, Taxonomic, and Functional Community Structure of Benthic Microbial Eukaryotes

Seagrass form high productive ecosystems in coastal environments. However, the effects of these coastal plants on the structure and function of the belowground eukaryotic microbiome remain elusive. In this study, we characterized the community of microbial eukaryotes (microeukaryotes) in both vegetated and unvegetated sediments using 18S rRNA gene amplicon sequencing and quantitative PCR. Analysis of sequencing data showed that the eelgrass (Zostera marina) colonization decreased the alpha diversity indices of benthic microeukaryotes. Apicomplexa represented an average of 83% of reads across all samples, with a higher proportion at the vegetated sites. The taxonomic community structure was significantly different between these two types of sediments, for which the concentration ofNH 4 + in sediment porewater and salinity could account. Phylogenetic analyses of long 18S rRNA genes (around 1,030 bp) indicated these apicomplexan parasites are closely related to gregarine Lecudina polymorpha. Determination of 18S rRNA gene abundances provided evidence that the eelgrass markedly promoted the biomass of the gregarine and all microeukaryotes in the seagrass-colonized sediments and confirmed that the gregarine was hosted by a polychaete species. Significantly higher gene abundances of heterotrophs and mixotrophs were found at the vegetated sites, which could be explained by the finer sediments and short supply of dissolved inorganic nitrogen, respectively. The pigmented protists were more abundant in 18S rRNA gene copies at the lower and higher pH levels than at the intermediate. Nevertheless, the fractions of heterotrophs and phototrophs in the community were significantly related to porewater N:P ratio. These results indicate that seagrass colonization significantly induces an increase in overall biomass and a decrease in diversity of benthic microeukaryotes, making them more heterotrophic. This study also highlights that the hotspot of eukaryotic parasites could be linked with the high productivity of a natural ecosystem.

Landscape Composition and Soil Physical-Chemical Properties Drive the Assemblages of Bacteria and Fungi in Conventional Vegetable Fields

The soil microbiome is crucial for improving the services and functioning of agroecosystems. Numerous studies have demonstrated the potential of soil physical–chemical properties in driving the belowground microbial assemblages in different agroecosystems. However, not much is known about the assemblage of bacteria and fungi in response to soil physical–chemical properties and the surrounding landscape composition in different vegetable fields of a highly intensive agricultural system. Here, we investigated the effects of soil physical–chemical properties and landscape composition on the community trends of bacteria and fungi in two different soil compartments (bulk and rhizospheric soils) of two different brassica crop types (Chinese cabbage and flower cabbage). The results revealed that bulk soil had a higher alpha diversity of both bacteria and fungi than rhizospheric soil. Each of the soil physical–chemical properties and landscape compositions contributed differently to driving the community structure of distinct bacterial and fungal taxa in both soil compartments and crop types. The higher proportions of forest, grassland, and cultivated land, along with the higher amount of soil calcium in flower cabbage fields, promote the assemblage of Gammaproteobacteria, Actinobacteria, Oxyophotobacteria, Agaricomycetes, and Eurotiomycetes. On the other hand, in Chinese cabbage fields, the increased amounts of iron, zinc, and manganese in the soil together with higher proportions of non-brassica crops in the surrounding landscape strongly support the assemblage of Deltaproteobacteria, Gemmatimonadetes, Bacilli, Clostridia, Alphaproteobacteria, an unknown bacterial species Subgroup-6, Mortierellomycetes, Rhizophlyctidomycetes, and Chytridiomycetes. The findings of this study provide the most comprehensive, comparative, and novel insights related to the bacterial and fungal responses in a highly intensive vegetable growing system for the improvement of the soil fertility and structure. These are important clues for the identification of key bacteria and fungi contributing to the plant–environment interactions and are of a practical significance for landscape-based ecological pest management.

distAngsd: Fast and accurate inference of genetic distances for Next Generation Sequencing data

Commonly used methods for inferring phylogenies were designed before the emergence of high throughput sequencing and can generally not accommodate the challenges associated with noisy, diploid sequencing data. In many applications, diploid genomes are still treated as haploid through the use of ambiguity characters; while the uncertainty in genotype calling - arising as a consequence of the sequencing technology - is ignored. In order to address this problem we describe two new probabilistic approaches for estimating genetic distances: distAngsd-geno and distAngsd-nuc, both implemented in a software suite named distAngsd. These methods are specifically designed for next generation sequencing data, utilize the full information from the data, and take uncertainty in genotype calling into account. Through extensive simulations, we show that these new methods are markedly more accurate and have more stable statistical behaviors than other currently available methods for estimating genetic distances - even for very low depth data with high error rates.