Human Adaptive Immunity Rescues an Inborn Error of Innate Immunity.

The molecular basis of the incomplete penetrance of monogenic disorders is unclear. We describe here eight related individuals with autosomal recessive TIRAP deficiency. Life-threatening staphylococcal disease occurred during childhood in the proband, but not in the other seven homozygotes. Responses to all Toll-like receptor 1/2 (TLR1/2), TLR2/6, and TLR4 agonists were impaired in the fibroblasts and leukocytes of all TIRAP-deficient individuals. However, the whole-blood response to the TLR2/6 agonist staphylococcal lipoteichoic acid (LTA) was abolished only in the index case individual, the only family member lacking LTA-specific antibodies (Abs). This defective response was reversed in the patient, but not in interleukin-1 receptor-associated kinase 4 (IRAK-4)-deficient individuals, by anti-LTA monoclonal antibody (mAb). Anti-LTA mAb also rescued the macrophage response in mice lacking TIRAP, but not TLR2 or MyD88. Thus, acquired anti-LTA Abs rescue TLR2-dependent immunity to staphylococcal LTA in individuals with inherited TIRAP deficiency, accounting for incomplete penetrance. Combined TIRAP and anti-LTA Ab deficiencies underlie staphylococcal disease in this patient.

Hormone-induced enhancer assembly requires an optimal level of hormone receptor multivalent interactions.

Transcription factors (TFs) activate enhancers to drive cell-specific gene programs in response to signals, but our understanding of enhancer assembly during signaling events is incomplete. Here, we show that androgen receptor (AR) forms condensates through multivalent interactions mediated by its N-terminal intrinsically disordered region (IDR) to orchestrate enhancer assembly in response to androgen signaling. AR IDR can be substituted by IDRs from selective proteins for AR condensation capacity and its function on enhancers. Expansion of the poly(Q) track within AR IDR results in a higher AR condensation propensity as measured by multiple methods, including live-cell single-molecule microscopy. Either weakening or strengthening AR condensation propensity impairs its heterotypic multivalent interactions with other enhancer components and diminishes its transcriptional activity. Our work reveals the requirement of an optimal level of AR condensation in mediating enhancer assembly and suggests that alteration of the fine-tuned multivalent IDR-IDR interactions might underlie AR-related human pathologies.

Label-free detection and profiling of individual solution-phase molecules.

Most chemistry and biology occurs in solution, in which conformational dynamics and complexation underlie behaviour and function. Single-molecule techniques1 are uniquely suited to resolving molecular diversity and new label-free approaches are reshaping the power of single-molecule measurements. A label-free single-molecule method2-16 capable of revealing details of molecular conformation in solution17,18 would allow a new microscopic perspective of unprecedented detail. Here we use the enhanced light-molecule interactions in high-finesse fibre-based Fabry-Pérot microcavities19-21 to detect individual biomolecules as small as 1.2 kDa, a ten-amino-acid peptide, with signal-to-noise ratios (SNRs) >100, even as the molecules are unlabelled and freely diffusing in solution. Our method delivers 2D intensity and temporal profiles, enabling the distinction of subpopulations in mixed samples. Notably, we observe a linear relationship between passage time and molecular radius, unlocking the potential to gather crucial information about diffusion and solution-phase conformation. Furthermore, mixtures of biomolecule isomers of the same molecular weight and composition but different conformation can also be resolved. Detection is based on the creation of a new molecular velocity filter window and a dynamic thermal priming mechanism that make use of the interplay between optical and thermal dynamics22,23 and Pound-Drever-Hall (PDH) cavity locking24 to reveal molecular motion even while suppressing environmental noise. New in vitro ways of revealing molecular conformation, diversity and dynamics can find broad potential for applications in the life and chemical sciences.

The HP1 homolog rhino anchors a nuclear complex that suppresses piRNA precursor splicing.

piRNAs guide an adaptive genome defense system that silences transposons during germline development. The Drosophila HP1 homolog Rhino is required for germline piRNA production. We show that Rhino binds specifically to the heterochromatic clusters that produce piRNA precursors, and that binding directly correlates with piRNA production. Rhino colocalizes to germline nuclear foci with Rai1/DXO-related protein Cuff and the DEAD box protein UAP56, which are also required for germline piRNA production. RNA sequencing indicates that most cluster transcripts are not spliced and that rhino, cuff, and uap56 mutations increase expression of spliced cluster transcripts over 100-fold. LacI::Rhino fusion protein binding suppresses splicing of a reporter transgene and is sufficient to trigger piRNA production from a trans combination of sense and antisense reporters. We therefore propose that Rhino anchors a nuclear complex that suppresses cluster transcript splicing and speculate that stalled splicing differentiates piRNA precursors from mRNAs.

Energy status dictates PD-L1 protein abundance and anti-tumor immunity to enable checkpoint blockade.

Aberrant energy status contributes to multiple metabolic diseases, including obesity, diabetes, and cancer, but the underlying mechanism remains elusive. Here, we report that ketogenic-diet-induced changes in energy status enhance the efficacy of anti-CTLA-4 immunotherapy by decreasing PD-L1 protein levels and increasing expression of type-I interferon (IFN) and antigen presentation genes. Mechanistically, energy deprivation activates AMP-activated protein kinase (AMPK), which in turn, phosphorylates PD-L1 on Ser283, thereby disrupting its interaction with CMTM4 and subsequently triggering PD-L1 degradation. In addition, AMPK phosphorylates EZH2, which disrupts PRC2 function, leading to enhanced IFNs and antigen presentation gene expression. Through these mechanisms, AMPK agonists or ketogenic diets enhance the efficacy of anti-CTLA-4 immunotherapy and improve the overall survival rate in syngeneic mouse tumor models. Our findings reveal a pivotal role for AMPK in regulating the immune response to immune-checkpoint blockade and advocate for combining ketogenic diets or AMPK agonists with anti-CTLA4 immunotherapy to combat cancer.

Large-scale migration into Britain during the Middle to Late Bronze Age.

Present-day people from England and Wales have more ancestry derived from early European farmers (EEF) than did people of the Early Bronze Age1. To understand this, here we generated genome-wide data from 793 individuals, increasing data from the Middle to the Late Bronze Age and Iron Age in Britain by 12-fold, and western and central Europe by 3.5-fold. Between 1000 and 875 BC, EEF ancestry increased in southern Britain (England and Wales) but not northern Britain (Scotland) due to incorporation of migrants who arrived at this time and over previous centuries, and who were genetically most similar to ancient individuals from France. These migrants contributed about half the ancestry of people of England and Wales from the Iron Age, thereby creating a plausible vector for the spread of early Celtic languages into Britain. These patterns are part of a broader trend of EEF ancestry becoming more similar across central and western Europe in the Middle to the Late Bronze Age, coincident with archaeological evidence of intensified cultural exchange2-6. There was comparatively less gene flow from continental Europe during the Iron Age, and the independent genetic trajectory in Britain is also reflected in the rise of the allele conferring lactase persistence to approximately 50% by this time compared to approximately 7% in central Europe where it rose rapidly in frequency only a millennium later. This suggests that dairy products were used in qualitatively different ways in Britain and in central Europe over this period.

OpenFold: retraining AlphaFold2 yields new insights into its learning mechanisms and capacity for generalization.

AlphaFold2 revolutionized structural biology with the ability to predict protein structures with exceptionally high accuracy. Its implementation, however, lacks the code and data required to train new models. These are necessary to (1) tackle new tasks, like protein-ligand complex structure prediction, (2) investigate the process by which the model learns and (3) assess the model's capacity to generalize to unseen regions of fold space. Here we report OpenFold, a fast, memory efficient and trainable implementation of AlphaFold2. We train OpenFold from scratch, matching the accuracy of AlphaFold2. Having established parity, we find that OpenFold is remarkably robust at generalizing even when the size and diversity of its training set is deliberately limited, including near-complete elisions of classes of secondary structure elements. By analyzing intermediate structures produced during training, we also gain insights into the hierarchical manner in which OpenFold learns to fold. In sum, our studies demonstrate the power and utility of OpenFold, which we believe will prove to be a crucial resource for the protein modeling community.

UAP56 couples piRNA clusters to the perinuclear transposon silencing machinery.

piRNAs silence transposons during germline development. In Drosophila, transcripts from heterochromatic clusters are processed into primary piRNAs in the perinuclear nuage. The nuclear DEAD box protein UAP56 has been previously implicated in mRNA splicing and export, whereas the DEAD box protein Vasa has an established role in piRNA production and localizes to nuage with the piRNA binding PIWI proteins Ago3 and Aub. We show that UAP56 colocalizes with the cluster-associated HP1 variant Rhino, that nuage granules containing Vasa localize directly across the nuclear envelope from cluster foci containing UAP56 and Rhino, and that cluster transcripts immunoprecipitate with both Vasa and UAP56. Significantly, a charge-substitution mutation that alters a conserved surface residue in UAP56 disrupts colocalization with Rhino, germline piRNA production, transposon silencing, and perinuclear localization of Vasa. We therefore propose that UAP56 and Vasa function in a piRNA-processing compartment that spans the nuclear envelope.

Experimental observation of non-Abelian topological charges and edge states.

In the last few decades, topological phase1-11 has emerged as a new classification of matter states beyond the Ginzburg-Landau symmetry-breaking paradigm. The underlying global invariant is usually well characterized by integers, such as Chern numbers or winding numbers-the Abelian charges12-15. Very recently, researchers proposed the notion of non-Abelian topological charges16-19, which possess non-commutative and fruitful braiding structures with multiple (more than one) bandgaps tangled together. Here we experimentally observe the non-Abelian topological charges in a time-reversal and inversion-symmetric transmission line network. The quaternion-valued non-Abelian topological charges are clearly mapped onto an eigenstate-frame sphere. Moreover, we find a non-Abelian quotient relation that provides a global perspective on the distribution of edge/domain-wall states. Our work opens the door towards characterization and manipulation of non-Abelian topological charges, which may lead to interesting observables such as trajectory-dependent Dirac/Weyl node collisions in two-dimensional systems16,17,20, admissible nodal line configurations in three dimensions16,19,20, and may provide insight into certain strongly correlated phases of twisted bilayer graphene21.

Genomic insights into the formation of human populations in East Asia.

The deep population history of East Asia remains poorly understood owing to a lack of ancient DNA data and sparse sampling of present-day people1,2. Here we report genome-wide data from 166 East Asian individuals dating to between 6000 BC and AD 1000 and 46 present-day groups. Hunter-gatherers from Japan, the Amur River Basin, and people of Neolithic and Iron Age Taiwan and the Tibetan Plateau are linked by a deeply splitting lineage that probably reflects a coastal migration during the Late Pleistocene epoch. We also follow expansions during the subsequent Holocene epoch from four regions. First, hunter-gatherers from Mongolia and the Amur River Basin have ancestry shared by individuals who speak Mongolic and Tungusic languages, but do not carry ancestry characteristic of farmers from the West Liao River region (around 3000 BC), which contradicts theories that the expansion of these farmers spread the Mongolic and Tungusic proto-languages. Second, farmers from the Yellow River Basin (around 3000 BC) probably spread Sino-Tibetan languages, as their ancestry dispersed both to Tibet-where it forms approximately 84% of the gene pool in some groups-and to the Central Plain, where it has contributed around 59-84% to modern Han Chinese groups. Third, people from Taiwan from around 1300 BC to AD 800 derived approximately 75% of their ancestry from a lineage that is widespread in modern individuals who speak Austronesian, Tai-Kadai and Austroasiatic languages, and that we hypothesize derives from farmers of the Yangtze River Valley. Ancient people from Taiwan also derived about 25% of their ancestry from a northern lineage that is related to, but different from, farmers of the Yellow River Basin, which suggests an additional north-to-south expansion. Fourth, ancestry from Yamnaya Steppe pastoralists arrived in western Mongolia after around 3000 BC but was displaced by previously established lineages even while it persisted in western China, as would be expected if this ancestry was associated with the spread of proto-Tocharian Indo-European languages. Two later gene flows affected western Mongolia: migrants after around 2000 BC with Yamnaya and European farmer ancestry, and episodic influences of later groups with ancestry from Turan.

A Non-canonical Role of YAP/TEAD Is Required for Activation of Estrogen-Regulated Enhancers in Breast Cancer.

YAP/TEAD are nuclear effectors of the Hippo pathway, regulating organ size and tumorigenesis largely through promoter-associated function. However, their function as enhancer regulators remains poorly understood. Through an in vivo proximity-dependent labeling (BioID) technique, we identified YAP1 and TEAD4 protein as co-regulators of ERα on enhancers. The binding of YAP1/TEAD4 to ERα-bound enhancers is augmented upon E2 stimulation and is required for the induction of E2/ERα target genes and E2-induced oncogenic cell growth. Furthermore, their enhancer binding is a prerequisite for enhancer activation marked by eRNA transcription and for the recruitment of the enhancer activation machinery component MED1. The binding of TEAD4 on active ERE-containing enhancers is independent of its DNA-binding behavior, and instead, occurs through protein-tethering trans-binding. Our data reveal a non-canonical function of YAP1 and TEAD4 as ERα cofactors in regulating cancer growth, highlighting the potential of YAP/TEAD as possible actionable drug targets for ERα+ breast cancer.

Phosphorylation of EZH2 by AMPK Suppresses PRC2 Methyltransferase Activity and Oncogenic Function.

Sustained energy starvation leads to activation of AMP-activated protein kinase (AMPK), which coordinates energy status with numerous cellular processes including metabolism, protein synthesis, and autophagy. Here, we report that AMPK phosphorylates the histone methyltransferase EZH2 at T311 to disrupt the interaction between EZH2 and SUZ12, another core component of the polycomb repressive complex 2 (PRC2), leading to attenuated PRC2-dependent methylation of histone H3 at Lys27. As such, PRC2 target genes, many of which are known tumor suppressors, were upregulated upon T311-EZH2 phosphorylation, which suppressed tumor cell growth both in cell culture and mouse xenografts. Pathologically, immunohistochemical analyses uncovered a positive correlation between AMPK activity and pT311-EZH2, and higher pT311-EZH2 correlates with better survival in both ovarian and breast cancer patients. Our finding suggests that AMPK agonists might be promising sensitizers for EZH2-targeting cancer therapies.

Fluorogenic DNA-PAINT for faster, low-background super-resolution imaging.

DNA-based points accumulation for imaging in nanoscale topography (DNA-PAINT) is a powerful super-resolution microscopy method that can acquire high-fidelity images at nanometer resolution. It suffers, however, from high background and slow imaging speed, both of which can be attributed to the presence of unbound fluorophores in solution. Here we present two-color fluorogenic DNA-PAINT, which uses improved imager probe and docking strand designs to solve these problems. These self-quenching single-stranded DNA probes are conjugated with a fluorophore and quencher at the terminals, which permits an increase in fluorescence by up to 57-fold upon binding and unquenching. In addition, the engineering of base pair mismatches between the fluorogenic imager probes and docking strands allowed us to achieve both high fluorogenicity and the fast binding kinetics required for fast imaging. We demonstrate a 26-fold increase in imaging speed over regular DNA-PAINT and show that our new implementation enables three-dimensional super-resolution DNA-PAINT imaging without optical sectioning.

MATR3-antisense LINE1 RNA meshwork scaffolds higher-order chromatin organization.

Long interspersed nuclear elements (LINEs) play essential roles in shaping chromatin states, while the factors that cooperate with LINEs and their roles in higher-order chromatin organization remain poorly understood. Here, we show that MATR3, a nuclear matrix protein, interplays with antisense LINE1 (AS L1) RNAs to form a meshwork via phase separation, providing a dynamic platform for chromatin spatial organization. MATR3 and AS L1 RNAs affect the nuclear localization of each other. After MATR3 depletion, the chromatin, particularly H3K27me3-modified chromatin, redistributes in the cell nuclei. Topologically associating domains (TADs) that highly transcribe MATR3-associated AS L1 RNAs show decreased intra-TAD interactions in both AML12 and ES cells. MATR3 depletion increases the accessibility of H3K27me3 domains adjacent to MATR3-associated AS L1, without affecting H3K27me3 modifications. Furthermore, amyotrophic lateral sclerosis (ALS)-associated MATR3 mutants alter biophysical features of the MATR3-AS L1 RNA meshwork and cause an abnormal H3K27me3 staining. Collectively, we reveal a role of the meshwork formed by MATR3 and AS L1 RNAs in gathering chromatin in the nucleus.

Colossal barocaloric effects in plastic crystals.

Refrigeration is of vital importance for modern society-for example, for food storage and air conditioning-and 25 to 30 per cent of the world's electricity is consumed for refrigeration1. Current refrigeration technology mostly involves the conventional vapour compression cycle, but the materials used in this technology are of growing environmental concern because of their large global warming potential2. As a promising alternative, refrigeration technologies based on solid-state caloric effects have been attracting attention in recent decades3-5. However, their application is restricted by the limited performance of current caloric materials, owing to small isothermal entropy changes and large driving magnetic fields. Here we report colossal barocaloric effects (CBCEs) (barocaloric effects are cooling effects of pressure-induced phase transitions) in a class of disordered solids called plastic crystals. The obtained entropy changes in a representative plastic crystal, neopentylglycol, are about 389 joules per kilogram per kelvin near room temperature. Pressure-dependent neutron scattering measurements reveal that CBCEs in plastic crystals can be attributed to the combination of extensive molecular orientational disorder, giant compressibility and highly anharmonic lattice dynamics of these materials. Our study establishes the microscopic mechanism of CBCEs in plastic crystals and paves the way to next-generation solid-state refrigeration technologies.

Visual pigment-deficient cones survive and mediate visual signaling despite the lack of outer segments.

Rhodopsin and cone opsins are essential for light detection in vertebrate rods and cones, respectively. It is well established that rhodopsin is required for rod phototransduction, outer segment disk morphogenesis, and rod viability. However, the roles of cone opsins are less well understood. In this study, we adopted a loss-of-function approach to investigate the physiological roles of cone opsins in mice. We showed that cones lacking cone opsins do not form normal outer segments due to the lack of disk morphogenesis. Surprisingly, cone opsin-deficient cones survive for at least 12 mo, which is in stark contrast to the rapid rod degeneration observed in rhodopsin-deficient mice, suggesting that cone opsins are dispensable for cone viability. Although the mutant cones do not respond to light directly, they maintain a normal dark current and continue to mediate visual signaling by relaying the rod signal through rod-cone gap junctions. Our work reveals a striking difference between the role of rhodopsin and cone opsins in photoreceptor viability.

Inhibition of EZH2 transactivation function sensitizes solid tumors to genotoxic stress.

Drugs that block the activity of the methyltransferase EZH2 are in clinical development for the treatment of non-Hodgkin lymphomas harboring EZH2 gain-of-function mutations that enhance its polycomb repressive function. We have previously reported that EZH2 can act as a transcriptional activator in castration-resistant prostate cancer (CRPC). Now we show that EZH2 inhibitors can also block the transactivation activity of EZH2 and inhibit the growth of CRPC cells. Gene expression and epigenomics profiling of cells treated with EZH2 inhibitors demonstrated that in addition to derepressing gene expression, these compounds also robustly down-regulate a set of DNA damage repair (DDR) genes, especially those involved in the base excision repair (BER) pathway. Methylation of the pioneer factor FOXA1 by EZH2 contributes to the activation of these genes, and interaction with the transcriptional coactivator P300 via the transactivation domain on EZH2 directly turns on the transcription. In addition, CRISPR-Cas9-mediated knockout screens in the presence of EZH2 inhibitors identified these BER genes as the determinants that underlie the growth-inhibitory effect of EZH2 inhibitors. Interrogation of public data from diverse types of solid tumors expressing wild-type EZH2 demonstrated that expression of DDR genes is significantly correlated with EZH2 dependency and cellular sensitivity to EZH2 inhibitors. Consistent with these findings, treatment of CRPC cells with EZH2 inhibitors dramatically enhances their sensitivity to genotoxic stress. These studies reveal a previously unappreciated mechanism of action of EZH2 inhibitors and provide a mechanistic basis for potential combination cancer therapies.

Bergamot essential oil improves CUMS-induced depression-like behaviour in rats by protecting the plasticity of hippocampal neurons.

Bergamot essential oil (BEO) is an extract of the bergamot fruit with significant neuroprotective effect. This study was to investigate the effects and the underlying mechanism of BEO in mitigating depression. GC-MS were used to identify its constituents. Antidepressive properties of BEO were evaluated by sucrose preference test (SPT), force swimming test (FST) and open field test (OFT). Nissl staining was used to determine the number of Nissl bodies in hippocampus (HIPP) of rats. Changes in HIPP dendritic length and dendritic spine density were detected by Golgi-Cox staining. Immunohistochemistry and Western blot were used to detect the postsynaptic density protein-95 (PSD-95) and synaptophysin (SYP) in the HIPP of rats. The enzyme-linked immunosorbent assay was used to determine the 5-hydroxytryptamine (5-HT), insulin-like growth factor 1 (IGF-1) and interleukin-1ß (IL-1ß) in the HIPP, serum and cerebrospinal fluid (CSF) of rats. Inhaled BEO significantly improved depressive behaviour in chronic unpredictable mild stress (CUMS) rats. BEO increased Nissl bodies, dendritic length and spine density, PSD-95 and SYP protein in the HIPP. Additionally, BEO upregulated serum 5-HT, serum and CSF IGF-1, while downregulating serum IL-1ß. Collectively, inhaled BEO mitigates depression by protecting the plasticity of hippocampal neurons, hence, providing novel insights into treatment of depression.

Integrated metabolome, full-length sequencing, and transcriptome analyses unveil the molecular mechanisms of color formation of the canary yellow and red bracts of Bougainvillea × buttiana 'Chitra'.

Bougainvillea is a typical tropical flower of great ornamental value due to its colorful bracts. The molecular mechanism behind color formation is not well-understood. Therefore, this research conducted metabolome analysis, transcriptome analysis, and multi-flux full-length sequencing in two color bracts of Bougainvillea × buttiana 'Chitra' to investigate the significantly different metabolites (SDMs) and differentially expressed genes (DEGs). Overall, 261 SDMs, including 62 flavonoids and 26 alkaloids, were detected, and flavonols and betalains were significantly differentially accumulated among the two bracts. Furthermore, the complete-length transcriptome of Bougainvillea × buttiana was also developed, which contained 512 493 non-redundant isoforms. Among them, 341 210 (66.58%) displayed multiple annotations in the KOG, GO, NR, KEGG, Pfam, Swissprot, and NT databases. RNA-seq findings revealed that 3610 DEGs were identified between two bracts. Co-expression analysis demonstrated that the DEGs and SDMs involved in flavonol metabolism (such as CHS, CHI, F3H, FLS, CYP75B1, kaempferol, and quercetin) and betacyanin metabolism (DODA, betanidin, and betacyanins) were the main contributors for the canary yellow and red bract formation, respectively. Further investigation revealed that several putative transcription factors (TFs) might interact with the promoters of the genes mentioned above. The expression profiles of the putative TFs displayed that they may positively and negatively regulate the structural genes' expression profiles. The data revealed a potential regulatory network between important genes, putative TFs, and metabolites in the flavonol and betacyanin biosynthesis of Bougainvillea × buttiana 'Chitra' bracts. These findings will serve as a rich genetic resource for future studies that could create new color bracts.

Castration alters the ileum microbiota of Holstein bulls and promotes beef flavor compounds.

BACKGROUND: In the beef industry, bull calves are usually castrated to improve flavor and meat quality; however, this can reduce their growth and slaughter performance. The gut microbiota is known to exert a significant influence on growth and slaughter performance. However, there is a paucity of research investigating the impact of castration on gut microbiota composition and its subsequent effects on slaughter performance and meat flavor. RESULT: The objective of this study was to examine the processes via which castration hinders slaughter productivity and enhances meat quality. Bull and castrated calves were maintained under the same management conditions, and at slaughter, meat quality was assessed, and ileum and epithelial tissue samples were obtained. The research employed metagenomic sequencing and non-targeted metabolomics techniques to investigate the makeup of the microbiota and identify differential metabolites. The findings of this study revealed the Carcass weight and eye muscle area /carcass weight in the bull group were significantly higher than those in the steer group. There were no significant differences in the length, width, and crypt depth of the ileum villi between the two groups. A total of 53 flavor compounds were identified in the two groups of beef, of which 16 were significantly higher in the steer group than in the bull group, and 5 were significantly higher in the bull group than in the steer group. In addition, bacteria, Eukaryota, and virus species were significantly separated between the two groups. The lipid metabolism pathways of α-linolenic acid, linoleic acid, and unsaturated fatty acids were significantly enriched in the Steers group. Compared with the steer group, the organic system pathway is significantly enriched in the bull group. The study also found that five metabolites (LPC (0:0/20:3), LPC (20:3/0:0), LPE (0:0/22:5), LPE (22:5/0:0), D-Mannosamine), and three species (s_Cloning_vector_Hsp70_LexA-HP1, s_Bacteroides_Coprophilus_CAG: 333, and s_Clostridium_nexile-CAG: 348) interfere with each other and collectively have a positive impact on the flavor compounds of beef. CONCLUSIONS: These findings provide a basic understanding that under the same management conditions, castration does indeed reduce the slaughter performance of bulls and improve the flavor of beef. Microorganisms and metabolites contribute to these changes through interactions.