Solving genomic puzzles: computational methods for metagenomic binning.

Metagenomics involves the study of genetic material obtained directly from communities of microorganisms living in natural environments. The field of metagenomics has provided valuable insights into the structure, diversity and ecology of microbial communities. Once an environmental sample is sequenced and processed, metagenomic binning clusters the sequences into bins representing different taxonomic groups such as species, genera, or higher levels. Several computational tools have been developed to automate the process of metagenomic binning. These tools have enabled the recovery of novel draft genomes of microorganisms allowing us to study their behaviors and functions within microbial communities. This review classifies and analyzes different approaches of metagenomic binning and different refinement, visualization, and evaluation techniques used by these methods. Furthermore, the review highlights the current challenges and areas of improvement present within the field of research.

ACR: metagenome-assembled prokaryotic and eukaryotic genome refinement tool.

Microbial genome recovery from metagenomes can further explain microbial ecosystem structures, functions and dynamics. Thus, this study developed the Additional Clustering Refiner (ACR) to enhance high-purity prokaryotic and eukaryotic metagenome-assembled genome (MAGs) recovery. ACR refines low-quality MAGs by subjecting them to iterative k-means clustering predicated on contig abundance and increasing bin purity through validated universal marker genes. Synthetic and real-world metagenomic datasets, including short- and long-read sequences, evaluated ACR's effectiveness. The results demonstrated improved MAG purity and a significant increase in high- and medium-quality MAG recovery rates. In addition, ACR seamlessly integrates with various binning algorithms, augmenting their strengths without modifying core features. Furthermore, its multiple sequencing technology compatibilities expand its applicability. By efficiently recovering high-quality prokaryotic and eukaryotic genomes, ACR is a promising tool for deepening our understanding of microbial communities through genome-centric metagenomics.

AFITbin: a metagenomic contig binning method using aggregate l-mer frequency based on initial and terminal nucleotides.

BACKGROUND: Using next-generation sequencing technologies, scientists can sequence complex microbial communities directly from the environment. Significant insights into the structure, diversity, and ecology of microbial communities have resulted from the study of metagenomics. The assembly of reads into longer contigs, which are then binned into groups of contigs that correspond to different species in the metagenomic sample, is a crucial step in the analysis of metagenomics. It is necessary to organize these contigs into operational taxonomic units (OTUs) for further taxonomic profiling and functional analysis. For binning, which is synonymous with the clustering of OTUs, the tetra-nucleotide frequency (TNF) is typically utilized as a compositional feature for each OTU. RESULTS: In this paper, we present AFIT, a new l-mer statistic vector for each contig, and AFITBin, a novel method for metagenomic binning based on AFIT and a matrix factorization method. To evaluate the performance of the AFIT vector, the t-SNE algorithm is used to compare species clustering based on AFIT and TNF information. In addition, the efficacy of AFITBin is demonstrated on both simulated and real datasets in comparison to state-of-the-art binning methods such as MetaBAT 2, MaxBin 2.0, CONCOT, MetaCon, SolidBin, BusyBee Web, and MetaBinner. To further analyze the performance of the purposed AFIT vector, we compare the barcodes of the AFIT vector and the TNF vector. CONCLUSION: The results demonstrate that AFITBin shows superior performance in taxonomic identification compared to existing methods, leveraging the AFIT vector for improved results in metagenomic binning. This approach holds promise for advancing the analysis of metagenomic data, providing more reliable insights into microbial community composition and function. AVAILABILITY: A python package is available at: https://github.com/SayehSobhani/AFITBin .

Expanding and improving nanobody repertoires using a yeast display method: Targeting SARS-CoV-2.

COVID-19, caused by the coronavirus SARS-CoV-2, represents a serious worldwide health issue, with continually emerging new variants challenging current therapeutics. One promising alternate therapeutic avenue is represented by nanobodies, small single-chain antibodies derived from camelids with numerous advantageous properties and the potential to neutralize the virus. For identification and characterization of a broad spectrum of anti-SARS-CoV-2 Spike nanobodies, we further optimized a yeast display method, leveraging a previously published mass spectrometry-based method, using B-cell complementary DNA from the same immunized animals as a source of VHH sequences. Yeast display captured many of the sequences identified by the previous approach, as well as many additional sequences that proved to encode a large new repertoire of nanobodies with high affinities and neutralization activities against different SARS-CoV-2 variants. We evaluated DNA shuffling applied to the three complementarity-determining regions of antiviral nanobodies. The results suggested a surprising degree of modularity to complementarity-determining region function. Importantly, the yeast display approach applied to nanobody libraries from immunized animals allows parallel interrogation of a vast number of nanobodies. For example, we employed a modified yeast display to carry out massively parallel epitope binning. The current yeast display approach proved comparable in efficiency and specificity to the mass spectrometry-based approach, while requiring none of the infrastructure and expertise required for that approach, making these highly complementary approaches that together appear to comprehensively explore the paratope space. The larger repertoires produced maximize the likelihood of discovering broadly specific reagents and those that powerfully synergize in mixtures.

Cryptic-site-specific antibodies to the SARS-CoV-2 receptor binding domain can retain functional binding affinity to spike variants.

IMPORTANCE: Multiple SARS-CoV-2 variants of concern have emerged and caused a significant number of infections and deaths worldwide. These variants of concern contain mutations that might significantly affect antigen-targeting by antibodies. It is therefore important to further understand how antibody binding and neutralization are affected by the mutations in SARS-CoV-2 variants. We highlighted how antibody epitope specificity can influence antibody binding to SARS-CoV-2 spike protein variants and neutralization of SARS-CoV-2 variants. We showed that weakened spike binding and neutralization of Beta (B.1.351) and Omicron (BA.1) variants compared to wildtype are not universal among the panel of antibodies and identified antibodies of a specific binding footprint exhibiting consistent enhancement of spike binding and retained neutralization to Beta variant. These data and analysis can inform how antigen-targeting by antibodies might evolve during a pandemic and prepare for potential future sarbecovirus outbreaks.

Maximum likelihood outperforms binning methods for detecting differences in abundance size spectra across environmental gradients.

Individual body size distributions (ISD) within communities are remarkably consistent across habitats and spatiotemporal scales and can be represented by size spectra, which are described by a power law. The focus of size spectra analysis is to estimate the exponent ( λ ) of the power law. A common application of size spectra studies is to detect anthropogenic pressures. Many methods have been proposed for estimating λ most of which involve binning the data, counting the abundance within bins, and then fitting an ordinary least squares regression in log-log space. However, recent work has shown that binning procedures return biased estimates of λ compared to procedures that directly estimate λ using maximum likelihood estimation (MLE). While it is clear that MLE produces less biased estimates of site-specific λ's, it is less clear how this bias affects the ability to test for changes in λ across space and time, a common question in the ecological literature. Here, we used simulation to compare the ability of two normalised binning methods (equal logarithmic and log2 bins) and MLE to (1) recapture known values of λ , and (2) recapture parameters in a linear regression measuring the change in λ across a hypothetical environmental gradient. We also compared the methods using two previously published body size datasets across a natural temperature gradient and an anthropogenic pollution gradient. Maximum likelihood methods always performed better than common binning methods, which demonstrated consistent bias depending on the simulated values of λ . This bias carried over to the regressions, which were more accurate when λ was estimated using MLE compared to the binning procedures. Additionally, the variance in estimates using MLE methods is markedly reduced when compared to binning methods. The error induced by binning methods can be of similar magnitudes as the variation previously published in experimental and observational studies, bringing into question the effect sizes of previously published results. However, while the methods produced different regression slope estimates, they were in qualitative agreement on the sign of those slopes (i.e. all negative or all positive). Our results provide further support for the direct estimation of λ and its relative variation across environmental gradients using MLE over the more common methods of binning.

Multimorbidity in middle-aged women and COVID-19: binary data clustering for unsupervised binning of rare multimorbidity features and predictive modeling.

BACKGROUND: Multimorbidity is typically associated with deficient health-related quality of life in mid-life, and the likelihood of developing multimorbidity in women is elevated. We address the issue of data sparsity in non-prevalent features by clustering the binary data of various rare medical conditions in a cohort of middle-aged women. This study aims to enhance understanding of how multimorbidity affects COVID-19 severity by clustering rare medical conditions and combining them with prevalent features for predictive modeling. The insights gained can guide the development of targeted interventions and improved management strategies for individuals with multiple health conditions. METHODS: The study focuses on a cohort of 4477 female patients, (aged 45-60) in Piedmont, Italy, and utilizes their multimorbidity data prior to the COVID-19 pandemic from their medical history from 2015 to 2019. The COVID-19 severity is determined by the hospitalization status of the patients from February to May 2020. Each patient profile in the dataset is depicted as a binary vector, where each feature denotes the presence or absence of a specific multimorbidity condition. By clustering the sparse medical data, newly engineered features are generated as a bin of features, and they are combined with the prevalent features for COVID-19 severity predictive modeling. RESULTS: From sparse data consisting of 174 input features, we have created a low-dimensional feature matrix of 17 features. Machine Learning algorithms are applied to the reduced sparsity-free data to predict the Covid-19 hospital admission outcome. The performance obtained for the corresponding models are as follows: Logistic Regression (accuracy 0.72, AUC 0.77, F1-score 0.69), Linear Discriminant Analysis (accuracy 0.7, AUC 0.77, F1-score 0.67), and Ada Boost (accuracy 0.7, AUC 0.77, F1-score 0.68). CONCLUSION: Mapping higher-dimensional data to a low-dimensional space can result in information loss, but reducing sparsity can be beneficial for Machine Learning modeling due to improved predictive ability. In this study, we addressed the issue of data sparsity in electronic health records and created a model that incorporates both prevalent and rare medical conditions, leading to more accurate and effective predictive modeling. The identification of complex associations between multimorbidity and the severity of COVID-19 highlights potential areas of focus for future research, including long COVID and intervention efforts.

Genome-Centric Metagenomics Insights into the Plastisphere-Driven Natural Degradation Characteristics and Mechanism of Biodegradable Plastics in Aquatic Environments.

Biodegradable plastics (BPs) are pervasively available as alternatives to traditional plastics, but their natural degradation characteristics and microbial-driven degradation mechanisms are poorly understood, especially in aquatic environments, the primary sink of plastic debris. Herein, the three-month dynamic degradation process of BPs (the copolymer of poly(butylene adipate-co-terephthalate) and polylactic acid (PLA) (PBAT/PLA) and single PLA) in a natural aquatic environment was investigated, with nonbiodegradable plastics polyvinyl chloride, polypropylene, and polystyrene as controls. PBAT/PLA showed the weight loss of 47.4% at 50 days and severe fragmentation within two months, but no significant decay for other plastics. The significant increase in the specific surface area and roughness and the weakening of hydrophobicity within the first month promoted microbial attachment to the PBAT/PLA surface. Then, a complete microbial succession occurred, including biofilm formation, maturation, and dispersion. Metagenomic analysis indicated that plastispheres selectively enriched degraders. Based on the functional genes involved in BPs degradation, a total of 16 high-quality metagenome-assembled genomes of degraders (mainly Burkholderiaceae) were recovered from the PBAT/PLA plastisphere. These microbes showed the greatest degrading potential at the biofilm maturation stage and executed the functions by PLA_depolymerase, polyesterase, hydrolase, and esterase. These findings will enhance understanding of BPs' environmental behavior and microbial roles on plastic degradation.

The value of hybrid genomes: Building two highly contiguous reference genome assemblies to advance Canis genomic studies.

Previous studies of canid population and evolutionary genetics have relied on high-quality domestic dog reference genomes that have been produced primarily for biomedical and trait mapping studies in dog breeds. However, the absence of highly contiguous genomes from other Canis species like the gray wolf and coyote, that represent additional distinct demographic histories, may bias inferences regarding interspecific genetic diversity and phylogenetic relationships. Here, we present single haplotype de novo genome assemblies for the gray wolf and coyote, generated by applying the trio-binning approach to long sequence reads generated from the genome of a female first-generation hybrid produced from a gray wolf and coyote mating. The assemblies were highly contiguous, with contig N50 sizes of 44.6 and 42.0 Mb for the wolf and coyote, respectively. Genome scaffolding and alignments between the two Canis assemblies and published dog reference genomes showed near complete collinearity, with one exception: a coyote-specific chromosome fission of chromosome 13 and fusion of the proximal portion of that chromosome with chromosome 8, retaining the Canis-typical haploid chromosome number of 2n = 78. We evaluated mapping quality for previous RADseq data from 334 canids and found nearly identical mapping quality and patterns among canid species and regional populations regardless of the genome used for alignment (dog, coyote, or gray wolf). These novel wolf and coyote genome reference assemblies will be important resources for proper and accurate inference of Canis demography, taxonomic evaluation, and conservation genetics.

Antimicrobial resistance and its risks evaluation in wetlands on the Qinghai-Tibetan Plateau.

Amidst the global antimicrobial resistance (AMR) crisis, antibiotic resistance has permeated even the most remote environments. To understand the dissemination and evolution of AMR in minimally impacted ecosystems, the resistome and mobilome of wetlands across the Qinghai-Tibetan Plateau and its marginal regions were scrutinized using metagenomic sequencing techniques. The composition of wetland microbiomes exhibits significant variability, with dominant phyla including Proteobacteria, Actinobacteria, Bacteroidetes, and Verrucomicrobia. Notably, a substantial abundance of Antibiotic Resistance Genes (ARGs) and Mobile Genetic Elements (MGEs) was detected, encompassing 17 ARG types, 132 ARG subtypes, and 5 types of MGEs (Insertion Sequences, Insertions Sequences, Genomic Islands, Transposons, and Integrative Conjugative Elements). No significant variance was observed in the prevalence of resistome and mobilome across different wetland types (i.e., the Yellow River, other rivers, lakes, and marshes) (R=-0.5882, P=0.607). The co-occurrence of 74 ARG subtypes and 22 MGEs was identified, underscoring the pivotal role of MGEs in shaping ARG pools within the Qinghai-Tibetan Plateau wetlands. Metagenomic binning and analysis of assembled genomes (MAGs) revealed that 93 out of 206 MAGs harbored ARGs (45.15 %). Predominantly, Burkholderiales, Pseudomonadales, and Enterobacterales were identified as the primary hosts of these ARGs, many of which represent novel species. Notably, a substantial proportion of ARG-carrying MAGs also contained MGEs, reaffirming the significance of MGEs in AMR dissemination. Furthermore, utilizing the arg_ranker framework for risk assessment unveiled severe contamination of high-risk ARGs across most plateau wetlands. Moreover, some prevalent human pathogens were identified as potential hosts for these high-risk ARGs, posing substantial transmission risks. This study aims to investigate the prevalence of resistome and mobilome in wetlands, along with evaluating the risk posed by high-risk ARGs. Such insights are crucial for informing environmental protection strategies and facilitating the management of water resources on the Qinghai-Tibetan Plateau.

Oxygen mediated mobilization and co-occurrence of antibiotic resistance in lab-scale bioreactor using metagenomic binning.

Sub-lethal levels of antibiotic stimulate bacteria to generate reactive oxygen species (ROS) that promotes emergence and spread of antibiotic resistance mediated by mobile genetic elements (MGEs). Nevertheless, the influence of dissolved oxygen (DO) levels on mobility of antibiotic resistance genes (ARGs) in response to ROS-induced stress remains elusive. Thus, the study employs metagenomic assembly and binning approaches to decipher mobility potential and co-occurrence frequency of ARGs and MGEs under hyperoxic (5.5-7 mgL- 1), normoxic (2.5-4 mgL- 1), and hypoxic (0.5-1 mgL- 1) conditions in lab-scale bioreactor for 6 months. Among 163 high-quality metagenome-assembled genomes (MAGs) recovered from 13 metagenomes, 42 MAGs harboured multiple ARGs and were assigned to priority pathogen group. Total ARG count increased by 4.3 and 2.5% in hyperoxic and normoxic, but decreased by 0.53% in hypoxic conditions after 150 days. On contrary, MGE count increased by 7.3-1.3% in all the DO levels, with only two ARGs showed positive correlation with MGEs in hypoxic compared to 20 ARGs under hyperoxic conditions. Opportunistic pathogens (Escherichia, Klebsiella, Clostridium, and Proteus) were detected as potential hosts of ARGs wherein co-localisation of critical ARG gene cassette (sul1, dfr1,adeF, and qacC) were identified in class 1 integron/Tn1 family transposons. Thus, enhanced co-occurrence frequency of ARGs with MGEs in pathogens suggested promotion of ARGs mobility under oxidative stress. The study offers valuable insights into ARG dissemination and hosts dynamics that is essential for controlling oxygen-related stress for mitigating MGEs and ARGs in the environment.

FFT-Based Probability Density Imaging of Euler Solutions.

When using traditional Euler deconvolution optimization strategies, it is difficult to distinguish between anomalies and their corresponding Euler tails (those solutions are often distributed outside the anomaly source, forming "tail"-shaped spurious solutions, i.e., misplaced Euler solutions, which must be removed or marked) with only the structural index. The nonparametric estimation method based on the normalized B-spline probability density (BSS) is used to separate the Euler solution clusters and mark different anomaly sources according to the similarity and density characteristics of the Euler solutions. For display purposes, the BSS needs to map the samples onto the estimation grid at the points where density will be estimated in order to obtain the probability density distribution. However, if the size of the samples or the estimation grid is too large, this process can lead to high levels of memory consumption and excessive computation times. To address this issue, a fast linear binning approximation algorithm is introduced in the BSS to speed up the computation process and save time. Subsequently, the sample data are quickly projected onto the estimation grid to facilitate the discrete convolution between the grid and the density function using a fast Fourier transform. A method involving multivariate B-spline probability density estimation based on the FFT (BSSFFT), in conjunction with fast linear binning appropriation, is proposed in this paper. The results of two random normal distributions show the correctness of the BSS and BSSFFT algorithms, which is verified via a comparison with the true probability density function (pdf) and Gaussian kernel smoothing estimation algorithms. Then, the Euler solutions of the two synthetic models are analyzed using the BSS and BSSFFT algorithms. The results are consistent with their theoretical values, which verify their correctness regarding Euler solutions. Finally, the BSSFFT is applied to Bishop 5X data, and the numerical results show that the comprehensive analysis of the 3D probability density distributions using the BSSFFT algorithm, derived from the Euler solution subset of x0,y0,z0, can effectively separate and locate adjacent anomaly sources, demonstrating strong adaptability.

On the Accurate Estimation of Information-Theoretic Quantities from Multi-Dimensional Sample Data.

Using information-theoretic quantities in practical applications with continuous data is often hindered by the fact that probability density functions need to be estimated in higher dimensions, which can become unreliable or even computationally unfeasible. To make these useful quantities more accessible, alternative approaches such as binned frequencies using histograms and k-nearest neighbors (k-NN) have been proposed. However, a systematic comparison of the applicability of these methods has been lacking. We wish to fill this gap by comparing kernel-density-based estimation (KDE) with these two alternatives in carefully designed synthetic test cases. Specifically, we wish to estimate the information-theoretic quantities: entropy, Kullback-Leibler divergence, and mutual information, from sample data. As a reference, the results are compared to closed-form solutions or numerical integrals. We generate samples from distributions of various shapes in dimensions ranging from one to ten. We evaluate the estimators' performance as a function of sample size, distribution characteristics, and chosen hyperparameters. We further compare the required computation time and specific implementation challenges. Notably, k-NN estimation tends to outperform other methods, considering algorithmic implementation, computational efficiency, and estimation accuracy, especially with sufficient data. This study provides valuable insights into the strengths and limitations of the different estimation methods for information-theoretic quantities. It also highlights the significance of considering the characteristics of the data, as well as the targeted information-theoretic quantity when selecting an appropriate estimation technique. These findings will assist scientists and practitioners in choosing the most suitable method, considering their specific application and available data. We have collected the compared estimation methods in a ready-to-use open-source Python 3 toolbox and, thereby, hope to promote the use of information-theoretic quantities by researchers and practitioners to evaluate the information in data and models in various disciplines.

A large-scale investigation and identification of methicillin-resistant Staphylococcus aureus based on peaks binning of matrix-assisted laser desorption ionization-time of flight MS spectra.

Recent studies have demonstrated that the matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) could be used to detect superbugs, such as methicillin-resistant Staphylococcus aureus (MRSA). Due to an increasingly clinical need to classify between MRSA and methicillin-sensitive Staphylococcus aureus (MSSA) efficiently and effectively, we were motivated to develop a systematic pipeline based on a large-scale dataset of MS spectra. However, the shifting problem of peaks in MS spectra induced a low effectiveness in the classification between MRSA and MSSA isolates. Unlike previous works emphasizing on specific peaks, this study employs a binning method to cluster MS shifting ions into several representative peaks. A variety of bin sizes were evaluated to coalesce drifted or shifted MS peaks to a well-defined structured data. Then, various machine learning methods were performed to carry out the classification between MRSA and MSSA samples. Totally 4858 MS spectra of unique S. aureus isolates, including 2500 MRSA and 2358 MSSA instances, were collected by Chang Gung Memorial Hospitals, at Linkou and Kaohsiung branches, Taiwan. Based on the evaluation of Pearson correlation coefficients and the strategy of forward feature selection, a total of 200 peaks (with the bin size of 10 Da) were identified as the marker attributes for the construction of predictive models. These selected peaks, such as bins 2410-2419, 2450-2459 and 6590-6599 Da, have indicated remarkable differences between MRSA and MSSA, which were effective in the prediction of MRSA. The independent testing has revealed that the random forest model can provide a promising prediction with the area under the receiver operating characteristic curve (AUC) at 0.8450. When comparing to previous works conducted with hundreds of MS spectra, the proposed scheme demonstrates that incorporating machine learning method with a large-scale dataset of clinical MS spectra may be a feasible means for clinical physicians on the administration of correct antibiotics in shorter turn-around-time, which could reduce mortality, avoid drug resistance and shorten length of stay in hospital in the future.

Genome-resolved metagenomics using environmental and clinical samples.

Recent advances in high-throughput sequencing technologies and computational methods have added a new dimension to metagenomic data analysis i.e. genome-resolved metagenomics. In general terms, it refers to the recovery of draft or high-quality microbial genomes and their taxonomic classification and functional annotation. In recent years, several studies have utilized the genome-resolved metagenome analysis approach and identified previously unknown microbial species from human and environmental metagenomes. In this review, we describe genome-resolved metagenome analysis as a series of four necessary steps: (i) preprocessing of the sequencing reads, (ii) de novo metagenome assembly, (iii) genome binning and (iv) taxonomic and functional analysis of the recovered genomes. For each of these four steps, we discuss the most commonly used tools and the currently available pipelines to guide the scientific community in the recovery and subsequent analyses of genomes from any metagenome sample. Furthermore, we also discuss the tools required for validation of assembly quality as well as for improving quality of the recovered genomes. We also highlight the currently available pipelines that can be used to automate the whole analysis without having advanced bioinformatics knowledge. Finally, we will highlight the most widely adapted and actively maintained tools and pipelines that can be helpful to the scientific community in decision making before they commence the analysis.

Single-spoke binning: Reducing motion artifacts in abdominal radial stack-of-stars imaging.

PURPOSE: To increase the effectiveness of respiratory gating in radial stack-of-stars MRI, particularly when imaging at high spatial resolutions or with multiple echoes. METHODS: Free induction decay (FID) navigators were integrated into a three-dimensional gradient echo radial stack-of-stars pulse sequence. These navigators provided a motion signal with a high temporal resolution, which allowed single-spoke binning (SSB): each spoke at each phase encode step was sorted individually to the corresponding motion state of the respiratory signal. SSB was compared with spoke-angle binning (SAB), in which all phase encode steps of one projection angle were sorted without the use of additional navigator data. To illustrate the benefit of SSB over SAB, images of a motion phantom and of six free-breathing volunteers were reconstructed after motion-gating using either method. Image sharpness was quantitatively compared using image gradient entropies. RESULTS: The proposed method resulted in sharper images of the motion phantom and free-breathing volunteers. Differences in gradient entropy were statistically significant (p = 0.03) in favor of SSB. The increased accuracy of motion-gating led to a decrease of streaking artifacts in motion-gated four-dimensional reconstructions. To consistently estimate respiratory signals from the FID-navigator data, specific types of gradient spoiler waveforms were required. CONCLUSION: SSB allowed high-resolution motion-corrected MR imaging, even when acquiring multiple gradient echo signals or large acquisition matrices, without sacrificing accuracy of motion-gating. SSB thus relieves restrictions on the choice of pulse sequence parameters, enabling the use of motion-gated radial stack-of-stars MRI in a broader domain of clinical applications.

Motion-resolved four-dimensional abdominal diffusion-weighted imaging using PROPELLER EPI (4D-DW-PROPELLER-EPI).

PURPOSE: To develop a distortion-free motion-resolved four-dimensional diffusion-weighted PROPELLER EPI (4D-DW-PROPELLER-EPI) technique for benefiting clinical abdominal radiotherapy (RT). METHODS: An improved abdominal 4D-DWI technique based on 2D diffusion-weighted PROPELLER-EPI (2D-DW-PROPELLER-EPI), termed 4D-DW-PROPELLER-EPI, was proposed to improve the frame rate of repeated data acquisition and produce distortion-free 4D-DWI images. Since the radial or PROPELLER sampling with golden-angle rotation can achieve an efficient k-space coverage with a flexible time-resolved acquisition, the golden-angle multi-blade acquisition was used in the proposed 4D-DW-PROPELLER-EPI to improve the performance of data sorting. A new k-space and blade (K-B) amplitude binning method was developed for the proposed 4D-DW-PROPELLER-EPI to optimize the number of blades and the k-space uniformity before performing conventional PROPELLER-EPI reconstruction, by using two metrics to evaluate the adequacy of the acquired data. The proposed 4D-DW-PROPELLER-EPI was preliminarily evaluated in both simulation experiments and in vivo experiments with varying frame rates and different numbers of repeated acquisition. RESULTS: The feasibility of achieving distortion-free 4D-DWI images by using the proposed 4D-DW-PROPELLER-EPI technique was demonstrated in both digital phantom and healthy subjects. Evaluation of the 4D completeness metrics shows that the K-B amplitude binning method could simultaneously improve the acquisition efficiency and data reconstruction performance for 4D-DW-PROPELLER-EPI. CONCLUSION: 4D-DW-PROPELLER-EPI with K-B amplitude binning is an advanced technique that can provide distortion-free 4D-DWI images for resolving respiratory motion, and may benefit the application of image-guided abdominal RT.

Home-based microbial solution to boost crop growth in low-fertility soil.

Soil microbial inoculants are expected to boost crop productivity under climate change and soil degradation. However, the efficiency of native vs commercialized microbial inoculants in soils with different fertility and impacts on resident microbial communities remain unclear. We investigated the differential plant growth responses to native synthetic microbial community (SynCom) and commercial plant growth-promoting rhizobacteria (PGPR). We quantified the microbial colonization and dynamic of niche structure to emphasize the home-field advantages for native microbial inoculants. A native SynCom of 21 bacterial strains, originating from three typical agricultural soils, conferred a special advantage in promoting maize growth under low-fertility conditions. The root : shoot ratio of fresh weight increased by 78-121% with SynCom but only 23-86% with PGPRs. This phenotype correlated with the potential robust colonization of SynCom and positive interactions with the resident community. Niche breadth analysis revealed that SynCom inoculation induced a neutral disturbance to the niche structure. However, even PGPRs failed to colonize the natural soil, they decreased niche breadth and increased niche overlap by 59.2-62.4%, exacerbating competition. These results suggest that the home-field advantage of native microbes may serve as a basis for engineering crop microbiomes to support food production in widely distributed poor soils.

The role of structural variants in pest adaptation and genome evolution of the Colorado potato beetle, Leptinotarsa decemlineata (Say).

Structural variation has been associated with genetic diversity and adaptation. Despite these observations, it is not clear what their relative importance is for evolution, especially in rapidly adapting species. Here, we examine the significance of structural polymorphisms in pesticide resistance evolution of the agricultural super-pest, the Colorado potato beetle, Leptinotarsa decemlineata. By employing a parent offspring trio sequencing procedure, we develop highly contiguous reference genomes to characterize structural variation. These updated assemblies represent >100-fold improvement of contiguity and include derived pest and ancestral nonpest individuals. We identify >200,000 structural variations, which appear to be nonrandomly distributed across the genome as they co-occur with transposable elements and genes. Structural variations intersect with exons in a large proportion of gene annotations (~20%) that are associated with insecticide resistance (including cytochrome P450s), development, and transcription. To understand the role structural variations play in adaptation, we measure their allele frequencies among an additional 57 individuals using whole genome resequencing data, which represents pest and nonpest populations of North America. Incorporating multiple independent tests to detect the signature of natural selection using SNP data, we identify 14 genes that are probably under positive selection, include structural variations, and SNPs of elevated frequency within the pest lineages. Among these, three are associated with insecticide resistance based on previous research. One of these genes, CYP4g15, is coinduced during insecticide exposure with glycosyltransferase-13, which is a duplicated gene enclosed within a structural variant adjacent to the CYP4g15 genic region. These results demonstrate the significance of structural variations as a genomic feature to describe species history, genetic diversity, and adaptation.

Identification of novel Carnobacterium maltaromaticum strains in bone marrow samples of patients with acute myeloid leukemia using a metagenomic binning approach.

The aim of this study aimed to examine the existence of a bacterial metagenome in the bone marrow of patients with acute myeloid leukemia (AML). We re-examined whole-genome sequencing data from the bone marrow samples of seven patients with AML, four of whom were remitted after treatment, for metagenomic analysis. After the removal of human reads, unmapped reads were used to profile the species-level composition. We used the metagenomic binning approach to confirm whether the identified taxon was a complete genome of known or novel strains. We observed a unique and novel microbial signature in which Carnobacterium maltaromaticum was the most abundant species in five patients with AML or remission. The complete genome of C. maltaromaticum "BMAML_KR01," which was observed in all samples, was 100% complete with 8.5% contamination and closely clustered with C. maltaromaticum strains DSM20730 and SF668 based on single nucleotide polymorphism variations. We identified five unique proteins that could contribute to cancer progression and 104 virulent factor proteins in the BMAML_KR01 genome. To our knowledge, this is the first report of a new strain of C. maltaromaticum in patients with AML. The presence of C. maltaromaticum and its new strain in patients indicates an urgent need to validate the existence of this bacterium and evaluate its pathophysiological role.