Miracle-dz/UniFAR-FTUs huggingface dataset:
Deconvolving multiplexed histone modifications in single cells

Recent advances have enabled mapping of histone modifications in single cells1–12, but current methods are constrained to profile only one histone modification per cell. Here we present an integrated experimental and computational framework, scChIX (single-cell chromatin immunocleavage and unmixing), to map multiple histone modifications in single cells. scChIX multiplexes two histone modification profiles together in single cells, then computationally deconvolves the signal using training data from respective histone modification profiles. We first validate this method using purified blood cells and show that although the two repressive marks, H3K27me3 and H3K9me3, are generally mutually exclusive, the transitions between the two regions can vary between cell types. Next we apply scChIX to a heterogenous cell population from mouse bone marrow to generate linked maps of active (H3K4me1) and repressive (H3K27me3) chromatin landscapes in single cells, where coordinates in the active modification map correspond to coordinates in the repressive map. Linked analysis reveals that immunoglobulin genes in the Igkv region are in a repressed chromatin state in pro-B cells, but become activated in B cells. Overall, scChIX unlocks systematic interrogation of the interplay between histone modifications in single cells.
Endosomes deliver ceramide phosphoethanolamine with unique acyl chain anchors to the cleavage furrow during male meiotic cytokinesis

Cell division, wherein one cell divides into two daughter cells, is fundamental to all living organisms. Cytokinesis, the final step in cell division, begins with the formation of an actomyosin contractile ring, positioned midway between the segregated chromosomes. Constriction of the ring with concomitant membrane deposition in a spatiotemporal manner generates a cleavage furrow that physically separates the cytoplasm. Unique lipids with specific biophysical properties have been shown to localize to intercellular bridges (also called midbody) connecting the two dividing cells; however, their biological roles and delivery mechanisms remain largely unknown. In this study, we show that Ceramide phosphoethanolamine (CPE), the structural analog of sphingomyelin, has unique acyl chain anchors in spermatocytes and is essential for meiotic cytokinesis. The head group of CPE is also important for spermatogenesis. We find that aberrant central spindle and contractile ring behavior but not mislocalization of phosphatidylinositol phosphates (PIPs) at the plasma membrane is responsible for the male meiotic cytokinesis defect in CPE deficient animals. Further, we demonstrate the enrichment of CPE in multivesicular bodies marked by Rab7, which in turn localize to cleavage furrow. Volume electron microscopy analysis using correlative light and focused ion beam scanning electron microscopy shows that CPE enriched Rab7 positive endosomes are juxtaposed on contractile ring material. Correlative light and transmission electron microscopy reveal Rab7 positive endosomes as a multivesicular body-like organelle that releases its intraluminal vesicles in the vicinity of ingressing furrows. Genetic ablation of Rab7 or expression of dominant negative Rab11 results in significant meiotic cytokinesis defects. Our results imply that endosomal delivery of CPE to ingressing membranes is crucial for meiotic cytokinesis.
Hippocampal gamma-band oscillopathy in a mouse model of Fragile X Syndrome

Fragile X syndrome (FXS) is the most common inherited form of intellectual disability arising from the loss of fragile X mental retardation protein (FMRP), a protein that plays a central role in neuronal function and plasticity. FXS patients show sensory hypersensitivity, hyperarousal and hippocampus-dependent learning deficits that can be recapitulated in the FMR1 KO mice. Enhanced metabotropic glutamate receptor (mGluR) signaling and muscarinic acetylcholine receptor (mAChR) signaling in the FMR1 KO mouse are implicated as the primary causes of the disease pathogenesis. Furthermore, glutamatergic kainate receptor (KAR) function is reduced in the cortex of the FMR1 KO mice. Of note, activation of these signaling pathways leads to slow gamma-range oscillations in the hippocampus in vitro and abnormal gamma oscillations have been reported in FMR1 KO mice and patients with FXS. Thus, we hypothesized that aberrant activation of these receptors leads to the observed gamma oscillopathy. We recorded gamma oscillations induced by either cholinergic agonist carbachol (CCh), mGluR1/5 agonist Dihydroxyphenylglycine (DHPG) or ionotropic glutamatergic agonist KA from the hippocampal CA3 in WT and FMR1 KO mice in vitro. We show a specific increase in the power of DHPG and CCh-induced gamma oscillations and reduction in the synchronicity of gamma oscillations induced by KA. We further elucidate an aberrant spiking activity during CCh-induced and kainate-induced gamma oscillations which may underlie the altered gamma oscillation synchronization in the FMR1 KO mice. Last, we also noted a reduced incidence of spontaneously-occurring hippocampal sharp wave-ripple events. Our study provides further evidence for aberrant hippocampal rhythms in the FMR1 KO mice and identifies potential signaling pathways underlying gamma band oscillopathy in FXS.
Ancient Gene Duplications, Rather Than Polyploidization, Facilitate Diversification of Petal Pigmentation Patterns in Clarkia gracilis (Onagraceae)

It has been suggested that gene duplication and polyploidization create opportunities for the evolution of novel characters. However, the connections between the effects of polyploidization and morphological novelties have rarely been examined. In this study, we investigated whether petal pigmentation patterning in an allotetraploid Clarkia gracilis has evolved as a result of polyploidization. C. gracilis is thought to be derived through a recent polyploidization event with two diploid species, C. amoena huntiana and an extinct species that is closely related to C. lassenensis. We reconstructed phylogenetic relationships of the R2R3-MYBs (the regulators of petal pigmentation) from two subspecies of C. gracilis and the two purported progenitors, C. a. huntiana and C. lassenensis. The gene tree reveals that these R2R3-MYB genes have arisen through duplications that occurred before the divergence of the two progenitor species, i.e., before polyploidization. After polyploidization and subsequent gene loss, only one of the two orthologous copies inherited from the progenitors was retained in the polyploid, turning it to diploid inheritance. We examined evolutionary changes in these R2R3-MYBs and in their expression, which reveals that the changes affecting patterning (including expression domain contraction, loss-of-function mutation, cis-regulatory mutation) occurred after polyploidization within the C. gracilis lineages. Our results thus suggest that polyploidization itself is not necessary in producing novel petal color patterns. By contrast, duplications of R2R3-MYB genes in the common ancestor of the two progenitors have apparently facilitated diversification of petal pigmentation patterns.
Resibufogenin Suppresses Triple-Negative Breast Cancer Angiogenesis by Blocking VEGFR2-Mediated Signaling Pathway

Resibufogenin (RBF), an active compound from Bufo bufonis, has been used for the treatment of multiple malignant cancers, including pancreatic cancer, colorectal cancer, and breast cancer. However, whether RBF could exert its antitumor effect by inhibiting angiogenesis remains unknown. Here, we aimed to explore the antiangiogenic activity of RBF and its underlying mechanism on human umbilical vein endothelial cell (HUVEC), and the therapeutic efficacy with regard to antiangiogenesis in vivo using two triple-negative breast cancer (TNBC) models. Our results demonstrated that RBF can inhibit the proliferation, migration, and tube formation of HUVECs in a dose-dependent manner. Spheroid sprouts were thinner and shorter after RBF treatment in vitro 3D spheroid sprouting assay. RBF also significantly suppressed VEGF-mediated vascular network formation in vivo Matrigel plug assay. In addition, Western blot analysis was used to reveal that RBF inhibited the phosphorylation of VEGFR2 and its downstream protein kinases FAK and Src in endothelial cells (ECs). Molecular docking simulations showed that RBF affected the phosphorylation of VEGFR2 by competitively binding to the ATP-bound VEGFR2 kinase domain, thus preventing ATP from providing phosphate groups. Finally, we found that RBF exhibited promising antitumor effect through antiangiogenesis in vivo without obvious toxicity. The present study first revealed the high antiangiogenic activity and the underlying molecular basis of RBF, suggesting that RBF could be a potential antiangiogenic agent for angiogenesis-related diseases.
Transcriptome and Co-expression Network Analyses Reveal Differential Gene Expression and Pathways in Response to Severe Drought Stress in Peanut (Arachis hypogaea L.)

Drought is one of the major abiotic stress factors limiting peanut production. It causes the loss of pod yield during the pod formation stage. Here, one previously identified drought-tolerant cultivar, “L422” of peanut, was stressed by drought (35 ± 5%) at pod formation stage for 5, 7, and 9 days. To analyze the drought effects on peanut, we conducted physiological and transcriptome analysis in leaves under well-watered (CK1, CK2, and CK3) and drought-stress conditions (T1, T2, and T3). By transcriptome analysis, 3,586, 6,730, and 8,054 differentially expressed genes (DEGs) were identified in “L422” at 5 days (CK1 vs T1), 7 days (CK2 vs T2), and 9 days (CK3 vs T3) of drought stress, respectively, and 2,846 genes were common DEGs among the three-time points. Furthermore, the result of weighted gene co-expression network analysis (WGCNA) revealed one significant module that was closely correlated between drought stress and physiological data. A total of 1,313 significantly up-/down-regulated genes, including 61 transcription factors, were identified in the module at three-time points throughout the drought stress stage. Additionally, six vital metabolic pathways, namely, “MAPK signaling pathway-plant,” “flavonoid biosynthesis,” “starch and sucrose metabolism,” “phenylpropanoid biosynthesis,” “glutathione metabolism,” and “plant hormone signal transduction” were enriched in “L422” under severe drought stress. Nine genes responding to drought tolerance were selected for quantitative real-time PCR (qRT-PCR) verification and the results agreed with transcriptional profile data, which reveals the reliability and accuracy of transcriptome data. Taken together, these findings could lead to a better understanding of drought tolerance and facilitate the breeding of drought-resistant peanut cultivars.
Biphenyl 2,3-Dioxygenase in Pseudomonas alcaliphila JAB1 Is Both Induced by Phenolics and Monoterpenes and Involved in Their Transformation

The involvement of bacterial aromatic ring-hydroxylating dioxygenases (ARHDs) in the degradation of aromatic pollutants, such as polychlorinated biphenyls (PCBs), has been well studied. However, there is considerable speculation as to the origin of this ability. One hypothesis is centered on a connection between the ability to degrade aromatic pollutants and the necessity of soil bacteria to cope with and/or utilize secondary plant metabolites (SPMs). To investigate this connection, we researched the involvement of biphenyl 2,3-dioxygenase (BPDO), an ARHD essential for the degradation of PCBs, in the metabolism of SPMs in the soil bacterium Pseudomonas alcaliphila JAB1, a versatile degrader of PCBs. We demonstrated the ability of the strain JAB1 to transform a variety of SPMs, namely the flavonoids apigenin, flavone, flavanone, naringenin, fisetin, quercetin, morin, and catechin, caffeic acid, trans-cinnamic acid, and the monoterpenes (S)-limonene and (R)-carvone. Of those, the transformation of flavone, flavanone, and (S)-limonene was conditioned by the activity of JAB1-borne BPDO and thus was researched in more detail, and we found evidence for the limonene monooxygenase activity of the BPDO. Furthermore, the bphA gene in the strain JAB1 was demonstrated to be induced by a wide range of SPMs, with monoterpenes being the strongest inducers of the SPMs tested. Thus, our findings contribute to the growing body of evidence that ARHDs not only play a role in the catabolism of aromatic pollutants, but also of natural plant-derived aromatics, and this study supports the hypothesis that ARHDs participate in ecological processes mediated by SPMs.
Alterations of the Human Gut Microbiota in Intrahepatic Cholestasis of Pregnancy

Background and Aims Women with severe intrahepatic cholestasis of pregnancy (ICP) are at higher risks of fetal complications and without effective treatments. Changes in gut microbiota in pregnancy were found to be related to the altered intestinal bile acid composition, so we aimed to explore the alterations of microbiota in the gut of ICP patients. Methods A total of 90 women were recruited, including 45 ICP patients and 45 healthy controls. The gut microbiota communities of ICP group were compared to control group through 16S ribosomal RNA gene sequencing. The results were then confirmed by real-time polymerase chain reaction (PCR) and generalized linear model (GLM). Furthermore, we analyzed the relationships between microbiota and the severity of ICP. Results A total of seven genera and nine taxa with differential abundances between the ICP patients and the controls were identified. All of the seven genera were verified through real-time PCR, and three key genera Parabacteroides, Flavonifractor, and Megamonas were confirmed by using the GLM model. Further analysis found that the genera Escherichia_Shigella, Olsenella, and Turicibacter were enriched in the severe ICP group, the microbial gene function related to biosynthesis of unsaturated fatty acids and propanoate metabolism were also increased in them. Conclusions Overall, our study was the first in Asia to demonstrate an association between gut microbiota and ICP. Our findings would contribute to a better understanding of the occurrence of ICP.
Tissue-Dependent Variation Profiles of Tea Quality-Related Metabolites in New Shoots of Tea Accessions

Several metabolites define tea quality in new tea shoots composed of leaf and stem. To improve tea quality for breeding, it is important to understand the tissue-dependent genetic mechanisms and metabolic network responsible for the profile of tea quality-related metabolites. We analyzed the volatiles and specialized metabolites as the tea quality-related metabolites in leaves and stems of new shoots in 30 tea accessions to understand the tissue variation and network between tea quality-related metabolites. Our results provided the tissue-dependent variation network in the tea quality-related metabolites, including volatiles in new leaves and stems in tea accessions. Each volatile content in tea accessions showed the coefficient of variation ranging from 58.7 to 221.9% and 54.2 to 318.3% in new leaves and new stems, respectively. The accumulation pattern of tea quality-related metabolites in new leaves and stems varied depending on the accession. When comparing tea genetic populations, the profile of tea quality-related metabolites of new leaves, but not new stems, was the key to distinguishing tea genetic populations by chemical indicators. We described the network between tea quality-related metabolites, especially the dense network in new leaves. These results also will provide the key information for metabolic engineering and the selection of breeding materials in tea plants based on the tea quality-related metabolites and aid in understanding their molecular mechanisms and network of metabolic variation.
Tolerance to Abiotic Factors of Microsclerotia and Mycelial Pellets From Metarhizium robertsii, and Molecular and Ultrastructural Changes During Microsclerotial Differentiation

Metarhizium species fungi are able to produce resistant structures termed microsclerotia, formed by compact and melanized threads of hyphae. These propagules are tolerant to desiccation and produce infective conidia; thus, they are promising candidates to use in biological control programs. In this study, we investigated the tolerance to both ultraviolet B (UV-B) radiation and heat of microsclerotia of Metarhizium robertsii strain ARSEF 2575. We also adapted the liquid medium and culture conditions to obtain mycelial pellets from the same isolate in order to compare these characteristics between both types of propagules. We followed the peroxisome biogenesis and studied the oxidative stress during differentiation from conidia to microsclerotia by transmission electron microscopy after staining with a peroxidase activity marker and by the expression pattern of genes potentially involved in these processes. We found that despite their twice smaller size, microsclerotia exhibited higher dry biomass, yield, and conidial productivity than mycelial pellets, both with and without UV-B and heat stresses. From the 16 genes measured, we found an induction after 96-h differentiation in the oxidative stress marker genes MrcatA, MrcatP, and Mrgpx; the peroxisome biogenesis factors Mrpex5 and Mrpex14/17; and the photoprotection genes Mrlac1 and Mrlac2; and Mrlac3. We concluded that an oxidative stress scenario is induced during microsclerotia differentiation in M. robertsii and confirmed that because of its tolerance to desiccation, heat, and UV-B, this fungal structure could be an excellent candidate for use in biological control of pests under tropical and subtropical climates where heat and UV radiation are detrimental to entomopathogenic fungi survival and persistence.
Welfare Assessment of 30 Dairy Goat Farms in the Midwestern United States

Dairy goat animal welfare assessment protocols have been developed and conducted in Europe and the United Kingdom for dairy goats; however, there are no published reports of large-scale welfare assessment for dairy goats on farms in the Midwestern United States (US). Therefore, the objective of this study was to perform welfare assessment of lactating dairy goats and identify the most prevalent welfare issues on 30 farms across the Midwestern US. Thirty dairy goat farms (self-selected) were enrolled in the study if they shipped milk for human consumption (regardless of herd size). The number of lactating does on each farm ranged from 34 to 6,500 goats, with a median number of 158 lactating does (mean ± SD: 602 ± 1,708 lactating does). The protocol used was developed from available literature on goat welfare assessment but modified for use in the Midwestern US. Observations were made without handling the animals and included 22 animal-based indicators evaluated at the group- and individual-level. The observations were conducted during ~3–5 h during a milking session (either morning or afternoon) and time in the home pen. Principal components analysis (PCA) was carried out on the welfare assessment data from each farm. The first two dimensions of the PCA explained 34.8% of the variation. The PCA biplot indicated correlations between indicators. The most prevalent conditions observed across the 30 farms included any knee calluses (80.9%), any claw overgrowth (51.4%), poor hygiene (14.9%), skin lesions (8.9%), poor hair coat condition (8.3%) and any ear pathology (8.0%). These results are the first to provide the Midwestern US dairy goat industry with information to improve commercial dairy goat welfare.
Evidence for Peripheral Immune Activation in Parkinson’s Disease

Background Accumulating evidence has revealed that peripheral immunity is involved in Parkinson’s disease (PD). However, the results regarding the percentage of T-cell subsets are inconsistent, and the changes of immunoglobins levels have been seldom studied in PD patients. Methods Serum levels of the percentage of T-cell subsets and immunoglobulins were measured in 761 PD patients and 761 age- and gender-matched healthy controls. The correlations between the variables of peripheral immune activation (PIA) and the clinical characteristics of PD were analyzed using correlation analysis. Results The pooled results showed that PD patients had higher proportional levels of CD3+ T and CD4+ T lymphocytes than healthy controls. CD8+ T cell percentages were similar in PD patients and controls, and the CD4/CD8 ratio was significantly higher in the PD population. No significant differences in IgG, IgA, or IgM levels between these two groups were found. CD4+ T cell percentage was inversely correlated with the H&Y stage, and IgG level was positively correlated with disease duration and UPDRS part III. Subgroup analyses showed that these associations existed in female patients, but not in male patients. Conclusion The enhanced immune activation in the peripheral system is indicated in PD, and dynamic alterations in CD4+ T cell percentage and IgG level suggest an active role for peripheral immunity in the disease progression, especially in female PD patients.
LINC02257, an Enhancer RNA of Prognostic Value in Colon Adenocarcinoma, Correlates With Multi-Omics Immunotherapy-Related Analysis in 33 Cancers

Accumulated evidence supports that long non-coding RNAs (lncRNAs) are involved significantly in the development of human cancers. Enhancer RNAs (eRNAs), a subtype of lncRNAs, have recently attracted much attention about their roles in carcinogenesis. Colon adenocarcinoma is one of the most commonly diagnosed tumors with unfavorable prognosis. It highlights the great significance of screening and identifying novel biomarkers. More importantly, it remains to be elucidated with respect to the function of eRNAs in colon adenocarcinoma, as is in pan-cancers. The expression of LINC02257 was determined based on the data obtained from The Cancer Genome Atlas (TCGA). Further evaluation was performed on the basis of the following analyses: clinicopathology and survival analysis, gene ontology (GO) terms, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, as well as multi-omics immunotherapy-related analysis and co-expression analysis. The statistical analysis was conducted in R software, and immune cell infiltration of LINC02257 expression in cancers was investigated by using the CIBERSORT algorithm. By large-scale data mining, our study highlighted that a total of 39 eRNA genes were associated with colon adenocarcinoma prognosis, among which 25 eRNAs showed significant associations with their predicted target genes. LINC02257 was identified as the most significant survival-associated eRNA, with DUSP10 as its target gene. Besides, the high expression of LINC02257 in colon adenocarcinoma was more vulnerable to unfavorable prognosis and correlated with various clinical characteristics. GO and KEGG analyses revealed that LINC02257 was closely correlated with extracellular matrix organization via the PI3K-Akt signaling pathway. Besides, LINC02257 expression correlated with a multi-omics analysis of 33 cancer types, such as survival analysis [overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI)] and immunotherapy-related analysis [tumor microenvironment (TME), tumor mutational burden (TMB), and microsatellite instability (MSI)]. Finally, we investigated the co-expression genes of LINC02257 and its potential signaling pathways across different cancer types. LINC02257 is screened and can function as an independent prognostic biomarker through the PI3K-Akt signaling pathway for colon adenocarcinoma. Simultaneously, LINC02257 may be a multifaceted and significant immunotherapy-related eRNA in different cancers.