Search

Inferring gene regulatory networks from single-cell RNA-sequencing trajectories has been an active area of research yet methods are still needed to identify regulators governing cell transitions. We developed DREAMIT (Dynamic Regulation of Expression Across Modules in Inferred Trajectories) to annotate transcription-factor activity along single-cell trajectory branches, using ensembles of relations to target genes. Using a benchmark representing several different tissues, as well as external validation with ATAC-Seq and Perturb-Seq data on hematopoietic cells, the method was found to have higher tissue-specific sensitivity and specificity over competing approaches.

OBJECTIVES: Anal cancer risk is elevated in MSM with HIV (MSMWH). Anal high-risk human papillomavirus (hr-HPV) infection is necessary but insufficient to develop high-grade squamous intraepithelial lesion (HSIL), the anal cancer precursor, suggesting additional factors. We sought to determine whether the microbiome of the anal canal is distinct by comparing it with the microbiome of stool. We also sought to determine whether changes in the anal microbiome are associated with HSIL among MSMWH. DESIGN: Cross-sectional comparison of the microbiome of the anal canal with the microbiome of stool in MSMWH and cross-sectional comparison of the anal microbiome of MSMWH with anal HSIL with the anal microbiome of MSMWH without anal HSIL. METHODS: Sterile swabs were used to sample the anus of MSMWH for microbiome and HPV testing, followed by high-resolution anoscopy. Stool samples were mailed from home. 16S sequencing was used for bacterial identification. Measures of alpha diversity, beta diversity, and differential abundance analysis were used to compare samples. RESULTS: One hundred sixty-six anal samples and 103 matching stool samples were sequenced. Beta diversity showed clustering of stool and anal samples. Of hr-HPV-positive MSMWH, 31 had HSIL and 13 had no SIL. Comparison of the microbiome between these revealed 28 different species. The highest-fold enrichment among MSMWH/hr-HPV/HSIL included pro-inflammatory and carcinogenic Prevotella, Parasuterella, Hungatella, Sneathia, and Fusobacterium species. The anti-inflammatory Anaerostipes caccae showed the greatest reduction among MSMWH/hr-HPV/HSIL. CONCLUSION: The anal microbiome is distinct from stool. A pro-inflammatory and carcinogenic environment may be associated with anal HSIL.

The jaguar (Panthera onca) is the largest living cat species native to the Americas and one of few large American carnivorans to have survived into the Holocene. However, the extent to which jaguar diversity declined during the end-Pleistocene extinction event remains unclear. For example, Pleistocene jaguar fossils from North America are notably larger than the average extant jaguar, leading to hypotheses that jaguars from this continent represent a now-extinct subspecies (Panthera onca augusta) or species (Panthera augusta). Here, we used a hybridization capture approach to recover an ancient mitochondrial genome from a large, late Pleistocene jaguar from Kingston Saltpeter Cave, Georgia, United States, which we sequenced to 26-fold coverage. We then estimated the evolutionary relationship between the ancient jaguar mitogenome and those from other extinct and living large felids, including multiple jaguars sampled across the species' current range. The ancient mitogenome falls within the diversity of living jaguars. All sampled jaguar mitogenomes share a common mitochondrial ancestor ~400 thousand years ago, indicating that the lineage represented by the ancient specimen dispersed into North America from the south at least once during the late Pleistocene. While genomic data from additional and older specimens will continue to improve understanding of Pleistocene jaguar diversity in the Americas, our results suggest that this specimen falls within the variation of extant jaguars despite the relatively larger size and geographic location and does not represent a distinct taxon.

Tandem repeats are frequent across the human genome, and variation in repeat length has been linked to a variety of traits. Recent improvements in long read sequencing technologies have the potential to greatly improve tandem repeat analysis, especially for long or complex repeats. Here, we introduce LongTR, which accurately genotypes tandem repeats from high-fidelity long reads available from both PacBio and Oxford Nanopore Technologies. LongTR is freely available at https://github.com/gymrek-lab/longtr and https://zenodo.org/doi/10.5281/zenodo.11403979 .

MOTIVATION: The alignment of sequencing reads is a critical step in the characterization of ancient genomes. However, reference bias and spurious mappings pose a significant challenge, particularly as cutting-edge wet lab methods generate datasets that push the boundaries of alignment tools. Reference bias occurs when reference alleles are favoured over alternative alleles during mapping, whereas spurious mappings stem from either contamination or when endogenous reads fail to align to their correct position. Previous work has shown that these phenomena are correlated with read length but a more thorough investigation of reference bias and spurious mappings for ancient DNA has been lacking. Here, we use a range of empirical and simulated palaeogenomic datasets to investigate the impacts of mapping tools, quality thresholds, and reference genome on mismatch rates across read lengths. RESULTS: For these analyses, we introduce AMBER, a new bioinformatics tool for assessing the quality of ancient DNA mapping directly from BAM-files and informing on reference bias, read length cut-offs and reference selection. AMBER rapidly and simultaneously computes the sequence read mapping bias in the form of the mismatch rates per read length, cytosine deamination profiles at both CpG and non-CpG sites, fragment length distributions, and genomic breadth and depth of coverage. Using AMBER, we find that mapping algorithms and quality threshold choices dictate reference bias and rates of spurious alignment at different read lengths in a predictable manner, suggesting that optimized mapping parameters for each read length will be a key step in alleviating reference bias and spurious mappings. AVAILABILITY AND IMPLEMENTATION: AMBER is available for noncommercial use on GitHub (https://github.com/tvandervalk/AMBER.git). Scripts used to generate and analyse simulated datasets are available on Github (https://github.com/sdolenz/refbias_scripts).

Recently, we reported the discovery of a novel endoglucanase of the glycoside hydrolase family 12 (GH12), designated IfCelS12A, from the haloalkaliphilic anaerobic bacterium Iocasia fonsfrigidae strain SP3-1, which was isolated from a hypersaline pond in the Samut Sakhon province of Thailand (ca. 2017). IfCelS12A exhibits high substrate specificity on carboxymethyl cellulose and amorphous cellulose but low substrate specificity on b-1,3;1,4-glucan. Unlike some endoglucanases of the GH12 family, IfCelS12A does not exhibit hydrolytic activity on crystalline cellulose (i.e., Avicel™). High-Pressure Liquid Chromatography (HPLC) and Thin Layer Chromatography (TLC) analyses of products resulting from IfCelS12-mediated hydrolysis indicate mode of action for this enzyme. Notably, IfCelS12A preferentially hydrolyzes cellotetraoses, cellopentaoses, and cellohexaoses with negligible activity on cellobiose or cellotriose. Kinetic analysis with cellopentaose and barely b-D-glucan as cellulosic substrates were conducted. On cellopentaose, IfCelS12A demonstrates a 16-fold increase in activity (KM = 0.27 mM; kcat = 0.36 s-1; kcat/KM = 1.34 mM-1 s-1) compared to the enzymatic hydrolysis of barley b-D-glucan (KM: 0.04 mM, kcat: 0.51 s-1, kcat/KM = 0.08 mM-1 s-1). Moreover, IfCelS12A enzymatic efficacy is stable in hypersaline sodium chlorids (NaCl) solutions (up to 10% NaCl). Specifically, IfCel12A retains notable activity after 24 h at 2M NaCl (10% saline solution). IfCelS12A used as a cocktail component with other cellulolytic enzymes and in conjunction with mobile sequestration platform technology offers additional options for deconstruction of ionic liquid-pretreated cellulosic feedstock. KEY POINTS: • IfCelS12A from an anaerobic alkaliphile Iocasia fronsfrigidae shows salt tolerance • IfCelS12A in cocktails with other enzymes efficiently degrades cellulosic biomass • IfCelS12A used with mobile enzyme sequestration platforms enhances hydrolysis.

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

The cellular components of tumors and their microenvironment play pivotal roles in tumor progression, patient survival, and the response to cancer treatments. Unveiling a comprehensive cellular profile within bulk tumors via single-cell RNA sequencing (scRNA-seq) data is crucial, as it unveils intrinsic tumor cellular traits that elude identification through conventional cancer subtyping methods. Our contribution, scBeacon, is a tool that derives cell-type signatures by integrating and clustering multiple scRNA-seq datasets to extract signatures for deconvolving unrelated tumor datasets on bulk samples. Through the employment of scBeacon on the The Cancer Genome Atlas (TCGA) cohort, we find cellular and molecular attributes within specific tumor categories, many with patient outcome relevance. We developed a tumor cell-type map to visually depict the relationships among TCGA samples based on the cell-type inferences.

Tuberculosis (TB) remains a leading cause of death from an infectious disease worldwide, partly due to a lack of effective strategies to screen and triage individuals with potential TB. Whole blood RNA signatures have been tested as biomarkers for TB, but have failed to meet the World Health Organizations (WHO) optimal target product profiles (TPP). Here, we use RNA sequencing and machine-learning to investigate the utility of plasma cell-free RNA (cfRNA) as a host-response biomarker for TB in cohorts from Uganda, Vietnam and Philippines. We report a 6-gene cfRNA signature, which differentiates TB-positive and TB-negative individuals with AUC = 0.95, 0.92, and 0.95 in test, training and validation, respectively. This signature meets WHO TPPs (sensitivity: 97.1% [95% CI: 80.9-100%], specificity: 85.2% [95% CI: 72.4-100%]) regardless of geographic location, sample collection method and HIV status. Overall, our results identify plasma cfRNA as a promising host response biomarker to diagnose TB.

MOTIVATION: Phylogenetic placement of a query sequence on a backbone tree is increasingly used across biomedical sciences to identify the content of a sample from its DNA content. The accuracy of such analyses depends on the density of the backbone tree, making it crucial that placement methods scale to very large trees. Moreover, a new paradigm has been recently proposed to place sequences on the species tree using single-gene data. The goal is to better characterize the samples and to enable combined analyses of marker-gene (e.g., 16S rRNA gene amplicon) and genome-wide data. The recent method DEPP enables performing such analyses using metric learning. However, metric learning is hampered by a need to compute and save a quadratically growing matrix of pairwise distances during training. Thus, the training phase of DEPP does not scale to more than roughly 10 000 backbone species, a problem that we faced when trying to use our recently released Greengenes2 (GG2) reference tree containing 331 270 species. RESULTS: This paper explores divide-and-conquer for training ensembles of DEPP models, culminating in a method called C-DEPP. While divide-and-conquer has been extensively used in phylogenetics, applying divide-and-conquer to data-hungry machine-learning methods needs nuance. C-DEPP uses carefully crafted techniques to enable quasi-linear scaling while maintaining accuracy. C-DEPP enables placing 20 million 16S fragments on the GG2 reference tree in 41 h of computation. AVAILABILITY AND IMPLEMENTATION: The dataset and C-DEPP software are freely available at https://github.com/yueyujiang/dataset_cdepp/.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Cell atlases serve as vital references for automating cell labeling in new samples, yet existing classification algorithms struggle with accuracy. Here we introduce SIMS (scalable, interpretable machine learning for single cell), a low-code data-efficient pipeline for single-cell RNA classification. We benchmark SIMS against datasets from different tissues and species. We demonstrate SIMS's efficacy in classifying cells in the brain, achieving high accuracy even with small training sets (

BACKGROUND: Single-cell transcriptome sequencing (scRNA-Seq) has allowed new types of investigations at unprecedented levels of resolution. Among the primary goals of scRNA-Seq is the classification of cells into distinct types. Many approaches build on existing clustering literature to develop tools specific to single-cell. However, almost all of these methods rely on heuristics or user-supplied parameters to control the number of clusters. This affects both the resolution of the clusters within the original dataset as well as their replicability across datasets. While many recommendations exist, in general, there is little assurance that any given set of parameters will represent an optimal choice in the trade-off between cluster resolution and replicability. For instance, another set of parameters may result in more clusters that are also more replicable. RESULTS: Here, we propose Dune, a new method for optimizing the trade-off between the resolution of the clusters and their replicability. Our method takes as input a set of clustering results-or partitions-on a single dataset and iteratively merges clusters within each partitions in order to maximize their concordance between partitions. As demonstrated on multiple datasets from different platforms, Dune outperforms existing techniques, that rely on hierarchical merging for reducing the number of clusters, in terms of replicability of the resultant merged clusters as well as concordance with ground truth. Dune is available as an R package on Bioconductor: https://www.bioconductor.org/packages/release/bioc/html/Dune.html . CONCLUSIONS: Cluster refinement by Dune helps improve the robustness of any clustering analysis and reduces the reliance on tuning parameters. This method provides an objective approach for borrowing information across multiple clusterings to generate replicable clusters most likely to represent common biological features across multiple datasets.

Intermuscular adipose tissue (IMAT) is a relatively understudied adipose depot located between muscle fibers. IMAT content increases with age and BMI and is associated with metabolic and muscle degenerative diseases; however, an understanding of the biological properties of IMAT and its interplay with the surrounding muscle fibers is severely lacking. In recent years, single-cell and nuclei RNA sequencing have provided us with cell type-specific atlases of several human tissues. However, the cellular composition of human IMAT remains largely unexplored due to the inherent challenges of its accessibility from biopsy collection in humans. In addition to the limited amount of tissue collected, the processing of human IMAT is complicated due to its proximity to skeletal muscle tissue and fascia. The lipid-laden nature of the adipocytes makes it incompatible with single-cell isolation. Hence, single nuclei RNA sequencing is optimal for obtaining high-dimensional transcriptomics at single-cell resolution and provides the potential to uncover the biology of this depot, including the exact cellular composition of IMAT. Here, we present a detailed protocol for nuclei isolation and library preparation of frozen human IMAT for single nuclei RNA sequencing. This protocol allows for the profiling of thousands of nuclei using a droplet-based approach, thus providing the capacity to detect rare and low-abundant cell types.

In droplet-based single-cell and single-nucleus RNA-seq assays, systematic contamination of ambient RNA molecules biases the quantification of gene expression levels. Existing methods correct the contamination for all genes globally. However, there lacks specific evaluation of correction efficacy for varying contamination levels. Here, we show that DecontX and CellBender under-correct highly contaminating genes, while SoupX and scAR over-correct lowly/non-contaminating genes. Here, we develop scCDC as the first method to detect the contamination-causing genes and only correct expression levels of these genes, some of which are cell-type markers. Compared with existing decontamination methods, scCDC excels in decontaminating highly contaminating genes while avoiding over-correction of other genes.

Since the onset of the coronavirus disease 2019 (COVID-19) pandemic, SARS-CoV-2 variants capable of breakthrough infections have attracted global attention. These variants have significant mutations in the receptor-binding domain (RBD) of the spike protein and the membrane (M) protein, which may imply an enhanced ability to evade immune responses. In this study, an examination of co-mutations within the spike RBD and their potential correlation with mutations in the M protein was conducted. The EVmutation method was utilized to analyze the distribution of the mutations to elucidate the relationship between the mutations in the spike RBD and the alterations in the M protein. Additionally, the Sequence-to-Sequence Transformer Model (S2STM) was employed to establish mapping between the amino acid sequences of the spike RBD and M proteins, offering a novel and efficient approach for streamlined sequence analysis and the exploration of their interrelationship. Certain mutations in the spike RBD, G339D-S373P-S375F and Q493R-Q498R-Y505, are associated with a heightened propensity for inducing mutations at specific sites within the M protein, especially sites 3 and 19/63. These results shed light on the concept of mutational synergy between the spike RBD and M proteins, illuminating a potential mechanism that could be driving the evolution of SARS-CoV-2.

UNLABELLED: A comprehensive microbial surveillance was conducted at NASAs Mars 2020 spacecraft assembly facility (SAF), where whole-genome sequencing (WGS) of 110 bacterial strains was performed. One isolate, designated 179-BFC-A-HST, exhibited less than 80% average nucleotide identity (ANI) to known species, suggesting a novel organism. This strain demonstrated high-level resistance [minimum inhibitory concentration (MIC) >256 mg/L] to third-generation cephalosporins, including ceftazidime, cefpodoxime, combination ceftazidime/avibactam, and the fourth-generation cephalosporin cefepime. The results of a comparative genomic analysis revealed that 179-BFC-A-HST is most closely related to Virgibacillus halophilus 5B73CT, sharing an ANI of 78.7% and a digital DNA-DNA hybridization (dDDH) value of 23.5%, while their 16S rRNA gene sequences shared 97.7% nucleotide identity. Based on these results and the recent recognition that the genus Virgibacillus is polyphyletic, strain 179-BFC-A-HST is proposed as a novel species of a novel genus, Tigheibacillus jepli gen. nov., sp. nov (type strain 179-BFC-A-HST = DSM 115946T = NRRL B-65666T), and its closest neighbor, V. halophilus, is proposed to be reassigned to this genus as Tigheibacillus halophilus comb. nov. (type strain 5B73CT = DSM 21623T = JCM 21758T = KCTC 13935T). It was also necessary to reclassify its second closest neighbor Virgibacillus soli, as a member of a novel genus Paracerasibacillus, reflecting its phylogenetic position relative to the genus Cerasibacillus, for which we propose Paracerasibacillus soli comb. nov. (type strain CC-YMP-6T = DSM 22952T = CCM 7714T). Within Amphibacillaceae (n = 64), P. soli exhibited 11 antibiotic resistance genes (ARG), while T. jepli encoded for 3, lacking any known β-lactamases, suggesting resistance from variant penicillin-binding proteins, disrupting cephalosporin efficacy. P. soli was highly resistant to azithromycin (MIC >64 mg/L) yet susceptible to cephalosporins and penicillins. IMPORTANCE: The significance of this research extends to understanding microbial survival and adaptation in oligotrophic environments, such as those found in SAF. Whole-genome sequencing of several strains isolated from Mars 2020 mission assembly cleanroom facilities, including the discovery of the novel species Tigheibacillus jepli, highlights the resilience and antimicrobial resistance (AMR) in clinically relevant antibiotic classes of microbes in nutrient-scarce settings. The study also redefines the taxonomic classifications within the Amphibacillaceae family, aligning genetic identities with phylogenetic data. Investigating ARG and virulence factors (VF) across these strains illuminates the microbial capability for resistance under resource-limited conditions while emphasizing the role of human-associated VF in microbial survival, informing sterilization practices and microbial management in similar oligotrophic settings beyond spacecraft assembly cleanrooms such as pharmaceutical and medical industry cleanrooms.

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

IntroductionLate-onset Alzheimer's disease (LOAD) has a strong genetic component. Participants in Long-Life Family Study (LLFS) exhibit delayed onset of dementia, offering a unique opportunity to investigate LOAD genetics.MethodsWe conducted a whole genome sequence analysis of 3475 LLFS members. Genetic associations were examined in six independent studies (N = 14,260) with a wide range of LOAD risk. Association analysis in a sub-sample of the LLFS cohort (N = 1739) evaluated the association of LOAD variants with beta amyloid (Aβ) levels.ResultsWe identified several single nucleotide polymorphisms (SNPs) in tight linkage disequilibrium within the MTUS2 gene associated with LOAD (rs73154407, p = 7.6 × 10-9). Association of MTUS2 variants with LOAD was observed in the five independent studies and was significantly stronger within high levels of Aβ42/40 ratio compared to lower amyloid.DiscussionMTUS2 encodes a microtubule associated protein implicated in the development and function of the nervous system, making it a plausible candidate to investigate LOAD biology.HighlightsLong-Life Family Study (LLFS) families may harbor late onset Alzheimer's dementia (LOAD) variants. LLFS whole genome sequence analysis identified MTUS2 gene variants associated with LOAD. The observed LLFS variants generalized to cohorts with wide range of LOAD risk. The association of MTUS2 with LOAD was stronger within high levels of beta amyloid. Our results provide evidence for MTUS2 gene as a novel LOAD candidate locus.

In the human genome, heterozygous sites refer to genomic positions with a different allele or nucleotide variant on the maternal and paternal chromosomes. Resolving these allelic differences by chromosomal copy, also known as phasing, is achievable on a short-read sequencer when using a library preparation method that captures long-range genomic information. TELL-Seq is a library preparation that captures long-range genomic information with the aid of molecular identifiers (barcodes). The same barcode is used to tag the reads derived from the same long DNA fragment within a range of up to 200 kilobases (kb), generating linked-reads. This strategy can be used to phase an entire genome. Here, we introduce a TELL-Seq protocol developed for targeted applications, enabling the phasing of enriched loci of varying sizes, purity levels, and heterozygosity. To validate this protocol, we phased 2-200 kb loci enriched with different methods: CRISPR/Cas9-mediated excision coupled with pulse-field electrophoresis for the longest fragments, CRISPR/Cas9-mediated protection from exonuclease digestion for mid-size fragments, and long PCR for the shortest fragments. All selected loci have known clinical relevance: BRCA1, BRCA2, MLH1, MSH2, MSH6, APC, PMS2, SCN5A-SCN10A, and PKI3CA. Collectively, the analyses show that TELL-Seq can accurately phase 2-200 kb targets using a short-read sequencer.

ObjectivesIndividuals increasingly experience delays or interruptions in schooling; we evaluate the association between these non-traditional education trajectories and mental health.MethodsUsing year-by-year education data for 7,501 National Longitudinal Survey of Youth 1979 participants, ages 14-48 (262,535 person-years of education data), we applied sequence analysis and a clustering algorithm to identify educational trajectory groups, incorporating both type and timing to credential. Linear regression models, adjusted for early-life confounders, evaluated relationships between educational trajectories and mental health component scores (MCS) from the 12-item short form instrument at age 50. We evaluated effect modification by race, gender, and race by gender.ResultsWe identified 24 distinct educational trajectories based on highest credential and educational timing. Compared to high school (HS) diplomas, < HS (beta=-3.41, 95%CI:-4.74,-2.07) and general educational development credentials (GEDs) predicted poorer MCS (beta=-2.07,95%CI:-3.16,-0.98). The following educational trajectories predicted better MCS: some college immediately after High School (beta=1.52, 95%CI:0.68,2.37), Associate degrees after long interruptions (beta=1.73, 95%CI:0.27,3.19), and graduate school soon after Bachelor's completion (beta=1.13, 95%CI:0.21,2.06). Compared to White men, Black women especially benefited from educational credentials higher than HS in predicting MCS.ConclusionsBoth type and timing of educational credential predicted mental health. Black women's mental higher especially benefited from higher educational credentials.

Pangenome references address biases of reference genomes by storing a representative set of diverse haplotypes and their alignment, usually as a graph. Alternate alleles determined by variant callers can be used to construct pangenome graphs, but advances in long-read sequencing are leading to widely available, high-quality phased assemblies. Constructing a pangenome graph directly from assemblies, as opposed to variant calls, leverages the graph's ability to represent variation at different scales. Here we present the Minigraph-Cactus pangenome pipeline, which creates pangenomes directly from whole-genome alignments, and demonstrate its ability to scale to 90 human haplotypes from the Human Pangenome Reference Consortium. The method builds graphs containing all forms of genetic variation while allowing use of current mapping and genotyping tools. We measure the effect of the quality and completeness of reference genomes used for analysis within the pangenomes and show that using the CHM13 reference from the Telomere-to-Telomere Consortium improves the accuracy of our methods. We also demonstrate construction of a Drosophila melanogaster pangenome.

MOTIVATION: Taxonomic classification of short reads and taxonomic profiling of metagenomic samples are well-studied yet challenging problems. The presence of species belonging to groups without close representation in a reference dataset is particularly challenging. While k-mer-based methods have performed well in terms of running time and accuracy, they tend to have reduced accuracy for such novel species. Thus, there is a growing need for methods that combine the scalability of k-mers with increased sensitivity. RESULTS: Here, we show that using locality-sensitive hashing (LSH) can increase the sensitivity of the k-mer-based search. Our method, which combines LSH with several heuristics techniques including soft lowest common ancestor labeling and voting, is more accurate than alternatives in both taxonomic classification of individual reads and abundance profiling. AVAILABILITY AND IMPLEMENTATION: CONSULT-II is implemented in C++, and the software, together with reference libraries, is publicly available on GitHub https://github.com/bo1929/CONSULT-II.

BACKGROUND: Pneumonia is a common cause of morbidity and mortality, yet a causative pathogen is identified in a minority of cases. Plasma microbial cell-free DNA sequencing may improve diagnostic yield in immunocompromised patients with pneumonia. METHODS: In this prospective, multicenter, observational study of immunocompromised adults undergoing bronchoscopy to establish a pneumonia etiology, plasma microbial cell-free DNA sequencing was compared to standardized usual care testing. Pneumonia etiology was adjudicated by a blinded independent committee. The primary outcome, additive diagnostic value, was assessed in the Per Protocol population (patients with complete testing results and no major protocol deviations) and defined as the percent of patients with an etiology of pneumonia exclusively identified by plasma microbial cell-free DNA sequencing. Clinical additive diagnostic value was assessed in the Per Protocol subgroup with negative usual care testing. RESULTS: Of 257 patients, 173 met Per Protocol criteria. A pneumonia etiology was identified by usual care in 52/173 (30.1%), plasma microbial cell-free DNA sequencing in 49/173 (28.3%) and the combination of both in 73/173 (42.2%) patients. Plasma microbial cell-free DNA sequencing exclusively identified an etiology of pneumonia in 21/173 patients (additive diagnostic value 12.1%, 95% confidence interval [CI], 7.7% to 18.0%, P < .001). In the Per Protocol subgroup with negative usual care testing, plasma microbial cell-free DNA sequencing identified a pneumonia etiology in 21/121 patients (clinical additive diagnostic value 17.4%, 95% CI, 11.1% to 25.3%). CONCLUSIONS: Non-invasive plasma microbial cell-free DNA sequencing significantly increased diagnostic yield in immunocompromised patients with pneumonia undergoing bronchoscopy and extensive microbiologic and molecular testing. CLINICAL TRIALS REGISTRATION: NCT04047719.

Ancient environmental DNA (aeDNA) from lake sediments has yielded remarkable insights for the reconstruction of past ecosystems, including suggestions of late survival of extinct species. However, translocation and lateral inflow of DNA in sediments can potentially distort the stratigraphic signal of the DNA. Using three different approaches on two short lake sediment cores of the Yamal peninsula, West Siberia, with ages spanning only the past hundreds of years, we detect DNA and identified mitochondrial genomes of multiple mammoth and woolly rhinoceros individuals-both species that have been extinct for thousands of years on the mainland. The occurrence of clearly identifiable aeDNA of extinct Pleistocene megafauna (e.g. >400 K reads in one core) throughout these two short subsurface cores, along with specificities of sedimentology and dating, confirm that processes acting on regional scales, such as extensive permafrost thawing, can influence the aeDNA record and should be accounted for in aeDNA paleoecology.

Marine macroalgae produce abundant and diverse polysaccharides, which contribute substantially to the organic matter exported to the deep ocean. Microbial degradation of these polysaccharides plays an important role in the turnover of macroalgal biomass. Various members of the Planctomycetes-Verrucomicrobia-Chlamydia (PVC) superphylum are degraders of polysaccharides in widespread anoxic environments. In this study, we isolated a novel anaerobic bacterial strain NLcol2T from microbial mats on the surface of marine sediments offshore Santa Barbara, CA, USA. Based on 16S ribosomal RNA (rRNA) gene and phylogenomic analyses, strain NLcol2T represents a novel species within the Pontiella genus in the Kiritimatiellota phylum (within the PVC superphylum). Strain NLcol2T is able to utilize various monosaccharides, disaccharides, and macroalgal polysaccharides such as agar and ɩ-carrageenan. A near-complete genome also revealed an extensive metabolic capacity for anaerobic degradation of sulfated polysaccharides, as evidenced by 202 carbohydrate-active enzymes (CAZymes) and 165 sulfatases. Additionally, its ability of nitrogen fixation was confirmed by nitrogenase activity detected during growth on nitrogen-free medium, and the presence of nitrogenases (nifDKH) encoded in the genome. Based on the physiological and genomic analyses, this strain represents a new species of bacteria that may play an important role in the degradation of macroalgal polysaccharides and with relevance to the biogeochemical cycling of carbon, sulfur, and nitrogen in marine environments. Strain NLcol2T (= DSM 113125T = MCCC 1K08672T) is proposed to be the type strain of a novel species in the Pontiella genus, and the name Pontiella agarivorans sp. nov. is proposed.IMPORTANCEGrowth and intentional burial of marine macroalgae is being considered as a carbon dioxide reduction strategy but elicits concerns as to the fate and impacts of this macroalgal carbon in the ocean. Diverse heterotrophic microbial communities in the ocean specialize in these complex polymers such as carrageenan and fucoidan, for example, members of the Kiritimatiellota phylum. However, only four type strains within the phylum have been cultivated and characterized to date, and there is limited knowledge about the metabolic capabilities and functional roles of related organisms in the environment. The new isolate strain NLcol2T expands the known substrate range of this phylum and further reveals the ability to fix nitrogen during anaerobic growth on macroalgal polysaccharides, thereby informing the issue of macroalgal carbon disposal.

Opinion 129 addresses the status of Firmicutes corrig. Gibbons and Murray 1978 (Approved Lists 1980). The name has the category 'division' and was included in the Approved Lists of Bacterial Names, although that category had previously been removed from the International Code of Nomenclature of Bacteria (1975 revision onwards). When the category 'phylum' was introduced into the International Code of Nomenclature of Prokaryotes (ICNP) in 2021, equivalence between 'phylum' and 'division' was not stipulated. Since the definition of the taxonomic categories and their relative order is one of the principal tasks of every code of nomenclature, the inclusion of Firmicutes corrig. Gibbons and Murray 1978 in the Approved Lists was an error. The name is either not validly published or illegitimate because its category is not covered by the ICNP. If Firmicutes corrig. Gibbons and Murray 1978 (Approved Lists 1980) was a validly published phylum name, it would be illegitimate because it would contravene Rule 8, which does not permit any deviation from the requirement to derive a phylum name from the name of the type genus. Since Firmicutes corrig. Gibbons and Murray 1978 is also part of a 'misfitting megaclassification' recognized in Opinion 128, the name is rejected, without any pre-emption regarding a hypothetically validly published name Firmicutes at the rank of phylum. Gracilicutes Gibbons and Murray 1978 (Approved Lists 1980) and Anoxyphotobacteriae Gibbons and Murray 1978 (Approved Lists 1980) are also rejected. The validly published phylum names have a variety of advantages over their not validly published counterparts and cannot be replaced with ad hoc names suggested in the literature. To ease the transition, it is recommended to mention the not validly published phylum names which strongly deviate in spelling from their validly published counterparts along with the latter in publications during the next years.

In single-cell RNA sequencing (scRNA-Seq), gene expression is assessed individually for each cell, allowing the investigation of developmental processes, such as embryogenesis and cellular differentiation and regeneration, at unprecedented resolution. In such dynamic biological systems, cellular states form a continuum, e.g., for the differentiation of stem cells into mature cell types. This process is often represented via a trajectory in a reduced-dimensional representation of the scRNA-Seq dataset. While many methods have been suggested for trajectory inference, it is often unclear how to handle multiple biological groups or conditions, e.g., inferring and comparing the differentiation trajectories of wild-type and knock-out stem cell populations. In this manuscript, we present condiments, a method for the inference and downstream interpretation of cell trajectories across multiple conditions. Our framework allows the interpretation of differences between conditions at the trajectory, cell population, and gene expression levels. We start by integrating datasets from multiple conditions into a single trajectory. By comparing the cells conditions along the trajectorys path, we can detect large-scale changes, indicative of differential progression or fate selection. We also demonstrate how to detect subtler changes by finding genes that exhibit different behaviors between these conditions along a differentiation path.

BackgroundThe accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.ResultsA total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.ConclusionsThese validated results show significant improvement over current bovine genome annotations.

BackgroundRNA-seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants (SNVs) in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme that mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung adenocarcinoma cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, the use of short-read RNA-seq has limited the community's ability to interrogate both forms of RNA variation simultaneously.ResultsWe employ long-read sequencing technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We develop a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generate nanopore data with high sequence accuracy from H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We apply our workflow to identify key inosine isoform associations to help clarify the prominence of ADAR in tumorigenesis.ConclusionsUltimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns.

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1 Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher's performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.

Sample multiplexing enables pooled analysis during single-cell RNA sequencing workflows, thereby increasing throughput and reducing batch effects. A challenge for all multiplexing techniques is to link sample-specific barcodes with cell-specific barcodes, then demultiplex sample identity post-sequencing. However, existing demultiplexing tools fail under many real-world conditions where barcode cross-contamination is an issue. We therefore developed deMULTIplex2, an algorithm inspired by a mechanistic model of barcode cross-contamination. deMULTIplex2 employs generalized linear models and expectation-maximization to probabilistically determine the sample identity of each cell. Benchmarking reveals superior performance across various experimental conditions, particularly on large or noisy datasets with unbalanced sample compositions.

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

OBJECTIVES: Viral lower respiratory tract infection (vLRTI) contributes to substantial morbidity and mortality in children. Diagnosis is typically confirmed by reverse transcriptase polymerase chain reaction (RT-PCR) of nasopharyngeal specimens in hospitalized patients; however, it is unknown whether nasopharyngeal detection accurately reflects presence of virus in the lower respiratory tract (LRT). This study evaluates agreement between viral detection from nasopharyngeal specimens by RT-PCR compared with metagenomic next-generation RNA sequencing (RNA-Seq) from tracheal aspirates (TAs). DESIGN: This is an analysis of of a seven-center prospective cohort study. SETTING: Seven PICUs within academic childrens hospitals in the United States. PATIENTS: Critically ill children (from 1 mo to 18 yr) who required mechanical ventilation via endotracheal tube for greater than or equal to 72 hours. INTERVENTIONS: We evaluated agreement in viral detection between paired upper and LRT samples. Results of clinical nasopharyngeal RT-PCR were compared with TA RNA-Seq. Positive and negative predictive agreement and Cohens Kappa were used to assess agreement. MEASUREMENTS AND MAIN RESULTS: Of 295 subjects with paired testing available, 200 (68%) and 210 (71%) had positive viral testing by RT-PCR from nasopharyngeal and RNA-Seq from TA samples, respectively; 184 (62%) were positive by both nasopharyngeal RT-PCR and TA RNA-Seq for a virus, and 69 (23%) were negative by both methods. Nasopharyngeal RT-PCR detected the most abundant virus identified by RNA-Seq in 92.4% of subjects. Among the most frequent viruses detected, respiratory syncytial virus demonstrated the highest degree of concordance (κ = 0.89; 95% CI, 0.83-0.94), whereas rhinovirus/enterovirus demonstrated lower concordance (κ = 0.55; 95% CI, 0.44-0.66). Nasopharyngeal PCR was more likely to detect multiple viruses than TA RNA-Seq (54 [18.3%] vs 24 [8.1%], p ≤ 0.001). CONCLUSIONS: Viral nucleic acid detection in the upper versus LRT reveals good overall agreement, but concordance depends on the virus. Further studies are indicated to determine the utility of LRT sampling or the use of RNA-Seq to determine LRTI etiology.

T(8;21)(q22;q22), which generates the AML1-ETO fusion oncoprotein, is a common chromosomal abnormality in acute myeloid leukemia (AML) patients. Despite having favorable prognosis, 40% of patients will relapse, highlighting the need for innovative models and application of the newest technologies to study t(8;21) leukemogenesis. Currently, available AML1-ETO mouse models have limited utility for studying the pre-leukemic stage because AML1-ETO produces mild hematopoietic phenotypes and no leukemic transformation. Conversely, overexpression of a truncated variant, AML1-ETO9a (AE9a), promotes fully penetrant leukemia and is too potent for studying pre-leukemic changes. To overcome these limitations, we devised a germline-transmitted Rosa26 locus AE9a knock-in mouse model that moderately overexpressed AE9a and developed leukemia with long latency and low penetrance. We observed pre-leukemic alterations in AE9a mice, including skewing of progenitors towards granulocyte/monocyte lineages and replating of stem and progenitor cells. Next, we performed single-cell RNA sequencing to identify specific cell populations that contribute to these pre-leukemic phenotypes. We discovered a subset of common myeloid progenitors that have heightened granulocyte/monocyte bias in AE9a mice. We also observed dysregulation of key hematopoietic transcription factor target gene networks, blocking cellular differentiation. Finally, we identified Sox4 activation as a potential contributor to stem cell self-renewal during the pre-leukemic stage.

Objective: The objective is to characterize transcriptomic profiles and immune cell composition and distribution in juvenile idiopathic arthritis (JIA) synovial biopsies, assess for associations of these features with clinical parameters, and compare JIA and rheumatoid arthritis (RA) synovial features. Methods: RNA sequencing (RNASeq) was performed on 24 samples, with pathway analysis and inference of relative abundance of immune cell subsets based on gene expression data. Two multiplex fluorescence immunohistochemistry (IHC) panels were performed on 28 samples (including 13 on which RNASeq was performed), staining for CD206− classical and CD206+ nonclassical macrophages, and CD8+ and CD4+ T and B lymphocytes. Data were compared to a published series of early RA synovial biopsies. Results: Pathway analysis of the most variably expressed genes (n = 339) identified a B and plasma cell signature as the main driver of heterogeneity in JIA synovia, with strong overlap between JIA and RA synovitis. Multiplex IHC confirmed heterogeneity of immune cell infiltration. M1‐like macrophage–rich synovial lining was associated with greater lining hypertrophy and higher (CD45+) pan–immune cell and CD8+ T cell infiltration. Conclusion: Our study indicates significant similarities between JIA and RA synovitis. Similar to RA, JIA synovia may be broadly categorized into two groups: (1) those with an inflammatory/adaptive immune transcriptomic signature, M1‐like macrophage and CD8+ T cell infiltration, and thicker, M1‐like macrophage–rich synovial lining, and (2) those with an M2‐like macrophage transcriptomic signature, greater M2/M1‐like macrophage ratios, and thinner, M2‐like macrophage–rich synovial lining. Synovial features were not significantly associated with clinical parameters, likely because of group size and heterogeneity. [ABSTRACT FROM AUTHOR]

Objective: Glucocorticoid‐induced tumor necrosis factor receptor superfamily–related protein (GITR), with its ligand (GITRL), plays an important role in CD4+ T cell–mediated autoimmunity. This study aimed to investigate the underlying mechanisms of GITRL in primary Sjögren syndrome (pSS). Methods: Patients with pSS and healthy controls were recruited. Serum GITRL and Th17‐related cytokines were determined. RNA sequencing was performed to decipher key signal pathways. Nonobese diabetes (NOD) mice were adopted as experimental Sjögren models and recombinant adeno‐associated virus (rAAV) transduction was conducted to verify the therapeutic potentials of targeting GITRL in vivo. Results: Serum GITRL was significantly higher in patients with pSS and showed a positive correlation with leukopenia, thrombocytopenia, autoantibodies, lung involvement, and disease activity. Serum GITRL was correlated with Th17‐related cytokines. GITRL promoted the expansion of Th17 and Th17.1 cells. Expansion of granulocyte‐macrophage colony‐stimulating factor positive (GM‐CSF+) CD4+ T cells induced by GITRL could be inhibited by blockade of GITRL. Moreover, GM‐CSF could stimulate GITRL expression on monocytes. RNA sequencing revealed mammalian target of rapamycin complexes 1 (mTORC1) might be the key modulator. The increased phosphorylation of S6 and STAT3 and the expansion of Th17 and Th17.1 cells induced by GITRL were effectively inhibited by rapamycin, suggesting a GITRL–mTORC1–GM‐CSF positive loop in pathogenic Th17 response in pSS. Administration of an rAAV vector expressing short hairpin RNA targeting GITRL alleviated disease progression in NOD mice. Conclusion: Our results identified the pathogenic role of GITRL in exacerbating disease activity and promoting pathogenic Th17 response in pSS through a GITRL–mTORC1–GM‐CSF loop. These findings suggest GITRL might be a promising therapeutic target in the treatment of pSS. [ABSTRACT FROM AUTHOR]

Purpose: With the rapid development of China's social work sector, the increasing job pressures, and risks of professional burnout among social workers have become more prevalent. This study examined the relationship between Chinese social workers' attitudes toward evidence-based practice (EBP) and burnout, exploring the mediating mechanisms of evidence-based knowledge (EBK) and service quality perception (SQP). Materials and methods: We applied PROCESS 4.2 macro in SPSS to analyze the data from 5,931 social workers, testing the sequential mediation effects of EBK and SQP between their attitudes toward EBP and burnout. Results: The findings revealed: (1) Attitudes toward EBP had significant indirect positive effects on burnout; (2) EBK partially mediated the relationship between EBP attitude and burnout; (3) SQP partially mediated the relationship between attitudes toward EBP and burnout; (4) Attitudes toward EBP had a sequential mediated effect on burnout through EBK and SQP. Discussion: The findings emphasize the need to implement targeted interventions and training programs to foster positive attitudes toward EBP, promote continuous professional development, and provide access to EBP resources. Moreover, nurturing EBK and SQP could help alleviate burnout by improving social workers' ability to address client issues and enhance their sense of confidence and accomplishment. Conclusion: This study fills a research gap by providing empirical evidence on the negative correlation between Chinese social workers' attitudes toward EBP and burnout, while demonstrating the mediating roles of EBK and SQP. [ABSTRACT FROM AUTHOR]

Aloesone is a bioactive natural product and biosynthetic precursor of rare glucosides found in rhubarb and some aloe plants including Aloe vera. This study aimed to investigate biocatalytic aloesone glycosylation and more than 400 uridine diphosphate-dependent glycosyltransferase (UGT) candidates, including multifunctional and promiscuous enzymes from a variety of plant species were assayed. As a result, 137 selective aloesone UGTs were discovered, including four from the natural producer rhubarb. Rhubarb UGT72B49 was further studied and its catalytic constants (k cat = 0.00092 ± 0.00003 s−1, K M = 30 ± 2.5 μM) as well as temperature and pH optima (50 °C and pH 7, respectively) were determined. We further aimed to find an efficient aloesone glycosylating enzyme with potential application for biocatalytic production of the glucoside. We discovered UGT71C1 from Arabidopsis thaliana as an efficient aloesone UGT showing a 167-fold higher catalytic efficiency compared to that of UGT72B49. Interestingly, sequence analysis of all the 137 newly identified aloesone UGTs showed that they belong to different phylogenetic groups, with the highest representation in groups B, D, E, F and L. Finally, our study indicates that aloesone C -glycosylation is highly specific and rare, since it was not possible to achieve in an efficient manner with any of the 422 UGTs assayed, including multifunctional GTs and 28 known C -UGTs. [ABSTRACT FROM AUTHOR]

Background: Previous research has shown a significant link between gut microbiota in children with autism spectrum disorder (ASD) and attention‐deficit hyperactivity disorder (ADHD). However, much remains unknown because of the heterogeneity of disorders and the potential confounders such as dietary patterns and control group variations. Methods: Children aged 6–12 years who had been clinically diagnosed with ASD and/or ADHD, their unaffected neurotypical siblings, and non‐related neurotypical volunteers were recruited cross‐sectionally. The ASD diagnosis was confirmed using the Autism Diagnostic Observation Schedule‐2 (ADOS‐2) in all patients, including those with ADHD. Standardized DNA extraction and sequencing methods were used to compare gut microbial alpha‐diversity among the groups. Dietary diversity was calculated from a standardized dietary questionnaire form. We compared the difference in gut microbiome between patients with ASD and/or ADHD with neurotypical siblings and non‐related neurotypical controls. Results: Ninety‐eight subjects were included in the study (18 with ASD, 19 with ADHD, 20 with both ASD and ADHD, 13 neurotypical siblings, and 28 non‐related neurotypical controls). The alpha‐diversity indices, such as Chao 1 and Shannon index, showed a significant difference between the groups in a Linear mixed‐effect model (F(4, 93) = 4.539, p =.02), (F(4, 93) = 3.185, p =.017), respectively. In a post‐hoc pairwise comparison, patients with ASD had lower alpha‐diversity compared with non‐related controls after Bonferroni correction. Dietary diversity shown in Shannon index did not differ among the groups (F(4, 84) = 1.494, p =.211). Conclusions: Our study indicates disorder‐specific microbiome differences in patients with ASD. In future research on gut microbiota in neurodevelopmental disorders, it is necessary to consider the impact of ASD and ADHD co‐occurrence, and strictly control for background information such as diet, to elucidate the gut–microbiota interaction in ASD and ADHD for exploring the potential of therapeutic interventions. [ABSTRACT FROM AUTHOR]

Objectives: Pendred syndrome, an autosomal recessive disorder, is often associated with pathogenic variants of the SLC26A4 gene that encodes the pendrin protein. Given its autosomal recessive inheritance, tracing the family history and screening siblings become crucial once a diagnosis of Pendred syndrome is confirmed. This case report aims to underscore the variability in inner ear morphology within a family diagnosed with Pendred syndrome, all carrying the same SLC26A4 gene mutation. Methods: A chart review and a review of the literature. Results: We present a family of 4, all of whom possess sensorineural hearing loss due to the same homozygous SLC26A4 variant c.919-2A>G. Intriguingly, clinical manifestations, especially inner ear deformities, displayed variability among family members. Notably, 1 family member exhibited a normal cochleovestibular structure morphology, which was rarely reported in the literature. Conclusions: This report highlights the significance of genetic testing and familial consultation when a proband exhibits typical Pendred syndrome symptoms. It also underscores that the inner ear morphology can exhibit variability among family members, even with the same homozygous SLC26A4 variant. [ABSTRACT FROM AUTHOR]

Aim: The transcriptional factor HIF‐1α is recognized for its contribution to cardioprotection against acute ischemia/reperfusion injury. Adaptation to chronic hypoxia (CH) is known to stabilize HIF‐1α and increase myocardial ischemic tolerance. However, the precise role of HIF‐1α in mediating the protective effect remains incompletely understood. Methods: Male wild‐type (WT) mice and mice with partial Hif1a deficiency (hif1a+/−) were exposed to CH for 4 weeks, while their respective controls were kept under normoxic conditions. Subsequently, their isolated perfused hearts were subjected to ischemia/reperfusion to determine infarct size, while RNA‐sequencing of isolated cardiomyocytes was performed. Mitochondrial respiration was measured to evaluate mitochondrial function, and western blots were performed to assess mitophagy. Results: We demonstrated enhanced ischemic tolerance in WT mice induced by adaptation to CH compared with their normoxic controls and chronically hypoxic hif1a+/− mice. Through cardiomyocyte bulk mRNA sequencing analysis, we unveiled significant reprogramming of cardiomyocytes induced by CH emphasizing mitochondrial processes. CH reduced mitochondrial content and respiration and altered mitochondrial ultrastructure. Notably, the reduced mitochondrial content correlated with enhanced autophagosome formation exclusively in chronically hypoxic WT mice, supported by an increase in the LC3‐II/LC3‐I ratio, expression of PINK1, and degradation of SQSTM1/p62. Furthermore, pretreatment with the mitochondrial division inhibitor (mdivi‐1) abolished the infarct size‐limiting effect of CH in WT mice, highlighting the key role of mitophagy in CH‐induced cardioprotection. Conclusion: These findings provide new insights into the contribution of HIF‐1α to cardiomyocyte survival during acute ischemia/reperfusion injury by activating the selective autophagy pathway. [ABSTRACT FROM AUTHOR]

Background: Non‐obstructive azoospermia is the most severe form of male infertility. A testicular biopsy is required for the diagnosis of non‐obstructive azoospermia, and the causal factors for non‐obstructive azoospermia remain unknown. Objectives: To reduce the risk of multiple biopsies and identify factors that contribute to non‐obstructive azoospermia, we proposed an integrated approach for the preoperative diagnosis and clinical management of non‐obstructive azoospermia by applying the chromosome‐spreading technique and whole‐exome sequencing. Materials and Methods: Between July 2020 and December 2022, after ruling out definitive obstructive azoospermia and non‐obstructive azoospermia patients with testicular volume < 6 mL, 20 patients with non‐obstructive azoospermia who underwent preoperative testicular diagnostic biopsy using testicular sperm aspiration were subjected to retrospective analysis. Results: Microscopic examination identified four patients with sperm cells, and 16 without sperm cells. Routine pathological analysis classified one patient as normal spermatogenesis, three as hypospermatogenesis, five as maturation arrest, nine as Sertoli cell‐only, and two as unable to judge. With chromosome‐spreading technology using routine cell suspension samples for microscopic examination, 18 patient diagnoses were validated, and two patients without a definitive diagnosis were supplemented. Detection of the Y chromosome and a well‐organized whole‐exome sequencing analysis revealed potential genetic factors. Discussion and Conclusion: The full use of testicular biopsy is beneficial for the diagnosis of azoospermia, as it avoids the risk of multiple biopsies. Moreover, in combination with whole‐exome sequencing, clinicians can obtain more information regarding the pathogenesis of non‐obstructive azoospermia, which may guide treatment. [ABSTRACT FROM AUTHOR]

This study uses data from a retrospective survey of the life history of participants from five cohorts born between 1930 and 1979, to analyze the role of adverse childhood experiences (ACEs). We found that ACEs have significant effects on the transition of young Chinese people into adulthood. Young people who have experienced more ACEs were more likely to transition to adulthood earlier, and accelerate their acceptance of adult social roles, and young adults with a high incidence of ACEs also have an increasing probability of falling into the disadvantaged or least successful trajectory to adulthood. More ACEs are associated with a decreased probability of being classified into the delayed and advantaged trajectories into adulthood characterized by high education and stable work, and an increased probability of being classified into the early and disadvantaged trajectories into adulthood characterized by unemployment, early marriage, and more children. [ABSTRACT FROM AUTHOR]

Lycium Barbarum Polysaccharides (LBP) can benefit lipid parameters such as total cholesterol, triglyceride, and high-density lipoprotein levels and upregulate the level of Firmicutes, increase the diversity of gut microbiota and reduce metabolic disorders, finally relieving weight gain of obese rats. But it cannot reverse the outcome of obesity. Over 30 differential metabolites and four pathways are altered by LBP. [ABSTRACT FROM AUTHOR]

In this study, a novel mitovirus, tentatively designated as "Alternaria alternata mitovirus 2" (AaMV2), was isolated from the fungus Alternaria alternata f. sp. mali causing apple leaf blotch disease. The complete genome of AaMV2 is 3,157 nucleotides in length, with an A+U content of 68.10%. The genome has a single large open reading frame (ORF) encoding an RNA-dependent RNA polymerase (RdRp) protein with a molecular mass of 98.10 kDa. BLAST analysis revealed that AaMV2 has the highest sequence identity to Leptosphaeria biglobosa mitovirus 6, with 79.76% and 82.86% identity at the amino acid and nucleotide level, respectively. Phylogenetic analysis suggested that AaMV2 is a new member of the genus Duamitovirus within the family Mitoviridae. This is the first report of the complete genome sequence analysis of a mitovirus in A. alternata. [ABSTRACT FROM AUTHOR]

• LGB-IPs is a new LGB-based optimal feature fusion model for the accurate prediction of STY phosphorylation sites. • LGB-IPs explores ten different feature descriptors and assesses its discriminative capability using five different classifiers. • Extensive cross-validation and independent assessment shows that LGB-IPs outperformed the single feature models consistently. SARS-CoV-2′s global spread has instigated a critical health and economic emergency, impacting countless individuals. Understanding the virus's phosphorylation sites is vital to unravel the molecular intricacies of the infection and subsequent changes in host cellular processes. Several computational methods have been proposed to identify phosphorylation sites, typically focusing on specific residue (S/T) or Y phosphorylation sites. Unfortunately, current predictive tools perform best on these specific residues and may not extend their efficacy to other residues, emphasizing the urgent need for enhanced methodologies. In this study, we developed a novel predictor that integrated all the residues (STY) phosphorylation sites information. We extracted ten different feature descriptors, primarily derived from composition, evolutionary, and position-specific information, and assessed their discriminative power through five classifiers. Our results indicated that Light Gradient Boosting (LGB) showed superior performance, and five descriptors displayed excellent discriminative capabilities. Subsequently, we identified the top two integrated features have high discriminative capability and trained with LGB to develop the final prediction model, LGB-IPs. The proposed approach shows an excellent performance on 10-fold cross-validation with an ACC, MCC, and AUC values of 0.831, 0.662, 0.907, respectively. Notably, these performances are replicated in the independent evaluation. Consequently, our approach may provide valuable insights into the phosphorylation mechanisms in SARS-CoV-2 infection for biomedical researchers. [ABSTRACT FROM AUTHOR]

Background: Congenital anomalies of the kidney and urinary tract (CAKUT) represent a frequent cause of pediatric kidney failure. CNVs, as a major class of genomic variations, can also affect miRNA regions. Common CNV corresponding miRNAs (cCNV-miRNAs) are functional variants regulating crucial processes which could affect urinary system development. Thus, we hypothesize that cCNV-miRNAs are associated with CAKUT occurrence and its expressivity. Methods: The extraction and filtering of common CNVs, identified in control samples deposited in publicly available databases gnomAD v2.1 and dbVar, were coupled with mapping of miRNA sequences using UCSC Genome Browser. After verification of the mapped miRNAs using referent miRBase V22.1, prioritization of cCNV-miRNA candidates has been performed using bioinformatic annotation and literature research. Genotyping of miRNA gene copy numbers for MIR9-3, MIR511, and MIR1299, was conducted on 221 CAKUT patients and 192 controls using TaqMan™ technology. Results: We observed significantly different MIR9-3 and MIR1299 gene copy number distribution between CAKUT patients and controls (Chi-square, P = 0.006 and P = 0.0002, respectively), while difference of MIR511 copy number distribution showed nominal significance (Chi-square, P = 0.027). The counts of less and more than two of MIR1299 copy numbers were more frequent within CAKUT patients compared to controls (P = 0.01 and P = 0.008, respectively) and also in cohort of patients with anomalies of the urinary tract compared to controls (P = 0.016 and P = 0.003, respectively). Conclusions: Copy number variations of miRNA genes represent a novel avenue in clarification of the inheritance complexity in CAKUT and provide potential evidence about the association of common genetic variation with CAKUT phenotypes. [ABSTRACT FROM AUTHOR]