TopHat: discovering splice junctions with RNA-Seq

Human-specific epigenomic states in spermatogenesis

Infertility is becoming increasingly common, affecting one in six people globally. Half of these cases can be attributed to male factors, many driven by abnormalities in the process of sperm development. Emerging evidence from genome-wide association studies, genetic screening of patient cohorts, and animal models highlights an important genetic contribution to spermatogenic defects, but comprehensive identification and characterization of the genes critical for male fertility remain lacking. High divergence of gene regulation in spermatogenic cells across species poses challenges for delineating the genetic pathways required for human spermatogenesis using common model organisms. In this study, we leveraged post-translational histone modification and gene transcription data for 15,491 genes in four mammalian species (human, rhesus macaque, mouse, and opossum), to identify human-specific patterns of gene regulation during spermatogenesis. We combined H3K27me3 ChIP-seq, H3K4me3 ChIP-seq, and RNA-seq data to define epigenetic states for each gene at two stages of spermatogenesis, pachytene spermatocytes and round spermatids, in each species. We identified 239 genes that are uniquely active, poised, or dynamically regulated in human spermatogenic cells distinct from the other three species. While some of these genes have been implicated in reproductive functions, many more have not yet been associated with human infertility and may be candidates for further molecular and epidemiologic studies. Our analysis offers an example of the opportunities provided by evolutionary and epigenomic data for broadly screening candidate genes implicated in reproduction, which might lead to discoveries of novel genetic targets for diagnosis and management of male infertility and male contraception.

Early-life milk replacer feeding mediates lipid metabolism disorders induced by colonic microbiota and bile acid profiles to reduce body weight in goat model

BACKGROUND

Dysregulation of lipid metabolism and its consequences on growth performance in young ruminants have attracted attention, especially in the context of alternative feeding strategies. This study aims to elucidate the effects of milk replacer (MR) feeding on growth, lipid metabolism, colonic epithelial gene expression, colonic microbiota composition and systemic metabolism in goat kids compared to breast milk (BM) feeding, addressing a critical knowledge gap in early life nutrition.

METHODS

Ten female goat kids were divided into 2 groups: those fed breast milk (BM group) and those fed a milk replacer (MR group). Over a period of 28 d, body weight was monitored and blood and tissue samples were collected for biochemical, transcriptomic and metabolomic analyses. Profiling of the colonial microbiota was performed using 16S rRNA gene sequencing. Intestinal microbiota transplantation (IMT) experiments in gnotobiotic mice were performed to validate causality.

RESULTS

MR-fed pups exhibited reduced daily body-weight gain due to impaired lipid metabolism as evidenced by lower serum and liver total cholesterol (TC) and non-esterified fatty acid (NEFA) concentrations. Transcriptomic analysis of the colonic epithelium revealed upregulated genes involved in negative regulation of lipid metabolism, concomitant with microbiota shifts characterized by a decrease in Firmicutes and an increase in Actinobacteria. Specifically, genera such as Bifidobacterium and Prevotella were enriched in the MR group, while Clostridium and Faecalibacterium were depleted. Metabolomics analyses confirmed alterations in bile acid and fatty acid metabolic pathways. IMT experiments in mice recapitulated the metabolic phenotype observed in MR-fed goats, confirming the role of the microbiota in modulating host lipid metabolism.

CONCLUSIONS

Milk replacer feeding in goat kids disrupts lipid metabolism and gut microbiota dynamics, resulting in reduced growth rates and metabolic alterations. These findings highlight the importance of early nutritional intervention on metabolic programming and suggest that modulation of the gut microbiota may be a target for improving growth and metabolic health in ruminants. This study contributes to the understanding of nutritional management strategies in livestock and their impact on animal health and productivity.

Mis-splicing of Mitotic Regulators Sensitizes SF3B1-Mutated Human HSCs to CHK1 Inhibition

Splicing factor SF3B1 mutations are frequent somatic lesions in myeloid neoplasms that transform hematopoietic stem cells (HSCs) by inducing mis-splicing of target genes. However, the molecular and functional consequences of SF3B1 mutations in human HSCs and progenitors (HSPCs) remain unclear. Here, we identify the mis-splicing program in human HSPCs as a targetable vulnerability by precise gene editing of SF3B1 K700E mutations in primary CD34+ cells. Mutant SF3B1 induced pervasive mis-splicing and reduced expression of genes regulating mitosis and genome maintenance leading to altered differentiation, delayed G2/M progression, and profound sensitivity to CHK1 inhibition (CHK1i). Mis-splicing or reduced expression of mitotic regulators BUBR1 and CDC27 delayed G2/M transit and promoted CHK1i sensitivity. Clinical CHK1i prexasertib selectively targeted SF3B1-mutant immunophenotypic HSCs and abrogated engraftment in vivo. These findings identify mis-splicing of mitotic regulators in SF3B1-mutant HSPCs as a targetable vulnerability engaged by pharmacological CHK1 inhibition. Significance: In this study, we engineer precise SF3B1 mutations in human HSPCs and identify CHK1 inhibition as a selective vulnerability promoted by mis-splicing of mitotic regulators. These findings uncover the mis-splicing program induced by mutant SF3B1 in human HSPCs and show that it can be therapeutically targeted by clinical CHK1 inhibitors.

Transcription factor ZmEREB97 regulates nitrate uptake in maize (Zea mays) roots

Maize (Zea mays L.) has very strong requirements for nitrogen. However, the molecular mechanisms underlying the regulations of nitrogen uptake and translocation in this species are not fully understood. Here, we report that an APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) transcription factor ZmEREB97 functions as an important regulator in the N signaling network in maize. Predominantly expressed and accumulated in main root and lateral root primordia, ZmEREB97 rapidly responded to nitrate treatment. By overlapping the analyses of differentially expressed genes and conducting a DAP-seq assay, we identified 1,446 potential target genes of ZmEREB97. Among these, 764 genes were coregulated in 2 lines of zmereb97 mutants. Loss of function of ZmEREB97 substantially weakened plant growth under both hydroponic and soil conditions. Physiological characterization of zmereb97 mutant plants demonstrated that reduced biomass and grain yield were both associated with reduced nitrate influx, decreased nitrate content, and less N accumulation. We further demonstrated that ZmEREB97 directly targets and regulates the expression of 6 ZmNRT genes by binding to the GCC-box-related sequences in gene promoters. Collectively, these data suggest that ZmEREB97 is a major positive regulator of the nitrate response and that it plays an important role in optimizing nitrate uptake, offering a target for improvement of nitrogen use efficiency in crops.

Architecture and dynamics of the abscisic acid gene regulatory network

The plant hormone abscisic acid (ABA) regulates essential processes in plant development and responsiveness to abiotic and biotic stresses. ABA perception triggers a post-translational signaling cascade that elicits the ABA gene regulatory network (GRN), encompassing hundreds of transcription factors (TFs) and thousands of transcribed genes. To further our knowledge of this GRN, we performed an RNA-seq time series experiment consisting of 14 time points in the 16 h following a one-time ABA treatment of 5-week-old Arabidopsis rosettes. During this time course, ABA rapidly changed transcription levels of 7151 genes, which were partitioned into 44 coexpressed modules that carry out diverse biological functions. We integrated our time-series data with publicly available TF-binding site data, motif data, and RNA-seq data of plants inhibited in translation, and predicted (i) which TFs regulate the different coexpression clusters, (ii) which TFs contribute the most to target gene amplitude, (iii) timing of engagement of different TFs in the ABA GRN, and (iv) hierarchical position of TFs and their targets in the multi-tiered ABA GRN. The ABA GRN was found to be highly interconnected and regulated at different amplitudes and timing by a wide variety of TFs, of which the bZIP family was most prominent, and upregulation of genes encompassed more TFs than downregulation. We validated our network models in silico with additional public TF-binding site data and transcription data of selected TF mutants. Finally, using a drought assay we found that the Trihelix TF GT3a is likely an ABA-induced positive regulator of drought tolerance.

Functional characterization of maize phytochrome-interacting factor 3 (ZmPIF3) in enhancing salt tolerance in arabidopsis

Soil salinization, a prevalent form of environmental stress, leads to significant soil desertification and impacts agricultural productivity by altering the internal soil environment, slowing cellular metabolism, and modifying cellular architecture. This results in a marked reduction in both the yield and diversity of crops. Maize, which is particularly susceptible to salt stress, serves as a critical model for studying these effects, making the elucidation of its molecular responses essential for crop improvement strategies. This study focuses on the phytochrome-interacting factor 3 (PIF3), previously known for its role in freezing tolerance, to assess its function in salt stress tolerance. Utilizing two transcript variants of maize ZmPIF3 (ZmPIF3.1 and ZmPIF3.2), we engineered Arabidopsis transgenic lines to overexpress these variants and analyzed their phenotypic, physiological, biochemical, and transcriptomic responses to salt stress. Our findings reveal that these transgenic lines displayed not only enhanced salt tolerance but also improved peroxide decomposition and reduced cellular membrane damage. Transcriptome analysis indicated significant roles of hormonal and Ca2+ signaling pathways, along with key transcription factors, in mediating the enhanced salt stress response. This research underscores a novel role for ZmPIF3 in plant salt stress tolerance, offering potential avenues for breeding salt-resistant crop varieties.

A dual regulatory role for the arbuscular mycorrhizal master regulator RAM1 in tomato

The REQUIRED FOR ARBUSCULAR MYCORRHIZATION1 (RAM1) transcription factor from the GRAS family is well known for its role as a master regulator of the arbuscular mycorrhizal (AM) symbiosis in dicotyledonous and monocotyledonous species, being essential in transcriptional reprogramming for the development and functionality of the arbuscules. In tomato, SlGRAS27 is the putative orthologue of RAM1 (here named SlRAM1), but has not yet been characterized. A reduced colonization of the root and impaired arbuscule formation were observed in SlRAM1-silenced plants, confirming the functional conservation of the RAM1 orthologue in tomato. However, unexpectedly, SlRAM1-overexpressing (UBIL:SlRAM1) plants also showed decreased mycorrhizal colonization. Analysis of non-mycorrhizal UBIL:SlRAM1 roots revealed an overall regulation of AM-related genes and a reduction of strigolactone biosynthesis. Moreover, external application of the strigolactone analogue GR244DO almost completely reversed the negative effects of SlRAM1 overexpression on the frequency of mycorrhization. However, it only partially recovered the pattern of arbuscule distribution observed in control plants. Our results strongly suggest that SlRAM1 has a dual regulatory role during mycorrhization and, in addition to its recognized action as a positive regulator of arbuscule development, it is also involved in different mechanisms for the negative regulation of mycorrhization, including the repression of strigolactone biosynthesis.

Epigenetic and transcriptional consequences in the endosperm of chemically induced transposon mobilization in Arabidopsis

Genomic imprinting, an epigenetic phenomenon leading to parent-of-origin-specific gene expression, has independently evolved in the endosperm of flowering plants and the placenta of mammals-tissues crucial for nurturing embryos. While transposable elements (TEs) frequently colocalize with imprinted genes and are implicated in imprinting establishment, direct investigations of the impact of de novo TE transposition on genomic imprinting remain scarce. In this study, we explored the effects of chemically induced transposition of the Copia element ONSEN on genomic imprinting in Arabidopsis thaliana. Through the combination of chemical TE mobilization and doubled haploid induction, we generated a line with 40 new ONSEN copies. Our findings reveal a preferential targeting of maternally expressed genes (MEGs) for transposition, aligning with the colocalization of H2A.Z and H3K27me3 in MEGs-both previously identified as promoters of ONSEN insertions. Additionally, we demonstrate that chemically-induced DNA hypomethylation induces global transcriptional deregulation in the endosperm, leading to the breakdown of MEG imprinting. This study provides insights into the consequences of chemically induced TE remobilization in the endosperm, revealing that chemically-induced epigenome changes can have long-term consequences on imprinted gene expression.

Heterochromatin formation and remodeling by IRTKS condensates counteract cellular senescence

Heterochromatin, a key component of the eukaryotic nucleus, is fundamental to the regulation of genome stability, gene expression and cellular functions. However, the factors and mechanisms involved in heterochromatin formation and maintenance still remain largely unknown. Here, we show that insulin receptor tyrosine kinase substrate (IRTKS), an I-BAR domain protein, is indispensable for constitutive heterochromatin formation via liquid‒liquid phase separation (LLPS). In particular, IRTKS droplets can infiltrate heterochromatin condensates composed of HP1α and diverse DNA-bound nucleosomes. IRTKS can stabilize HP1α by recruiting the E2 ligase Ubc9 to SUMOylate HP1α, which enables it to form larger phase-separated droplets than unmodified HP1α. Furthermore, IRTKS deficiency leads to loss of heterochromatin, resulting in genome-wide changes in chromatin accessibility and aberrant transcription of repetitive DNA elements. This leads to activation of cGAS-STING pathway and type-I interferon (IFN-I) signaling, as well as to the induction of cellular senescence and senescence-associated secretory phenotype (SASP) responses. Collectively, our findings establish a mechanism by which IRTKS condensates consolidate constitutive heterochromatin, revealing an unexpected role of IRTKS as an epigenetic mediator of cellular senescence.

Fitness adaptations of Japanese encephalitis virus in pigs following vector-free serial passaging

Japanese encephalitis virus (JEV) is a zoonotic mosquito-transmitted Flavivirus circulating in birds and pigs. In humans, JEV can cause severe viral encephalitis with high mortality. Considering that vector-free direct virus transmission was observed in experimentally infected pigs, JEV introduction into an immunologically naïve pig population could result in a series of direct transmissions disrupting the alternating host cycling between vertebrates and mosquitoes. To assess the potential consequences of such a realistic scenario, we passaged JEV ten times in pigs. This resulted in higher in vivo viral replication, increased shedding, and stronger innate immune responses in pigs. Nevertheless, the viral tissue tropism remained similar, and frequency of direct transmission was not enhanced. Next generation sequencing showed single nucleotide deviations in 10% of the genome during passaging. In total, 25 point mutations were selected to reach a frequency of at least 35% in one of the passages. From these, six mutations resulted in amino acid changes located in the precursor of membrane, the envelope, the non-structural 3 and the non-structural 5 proteins. In a competition experiment with two lines of passaging, the mutation M374L in the envelope protein and N275D in the non-structural protein 5 showed a fitness advantage in pigs. Altogether, the interruption of the alternating host cycle of JEV caused a prominent selection of viral quasispecies as well as selection of de novo mutations associated with fitness gains in pigs, albeit without enhancing direct transmission frequency.

Unexpected nuclear hormone receptor and chromatin dynamics regulate estrous cycle dependent gene expression

Chromatin changes in response to estrogen and progesterone are well established in cultured cells, but how they control gene expression under physiological conditions is largely unknown. To address this question, we examined in vivo estrous cycle dynamics of mouse uterus hormone receptor occupancy, chromatin accessibility and chromatin structure by combining RNA-seq, ATAC-seq, HiC-seq and ChIP-seq. Two estrous cycle stages were chosen for these analyses, diestrus (highest estrogen) and estrus (highest progesterone). Unexpectedly, rather than alternating with each other, estrogen receptor alpha (ERα) and progesterone receptor (PGR) were co-bound during diestrus and lost during estrus. Motif analysis of open chromatin followed by hypoxia inducible factor 2A (HIF2A) ChIP-seq and conditional uterine deletion of this transcription factor revealed a novel role for HIF2A in regulating diestrus gene expression patterns that were independent of either ERα or PGR binding. Proteins in complex with ERα included PGR and cohesin, only during diestrus. Combined with HiC-seq analyses, we demonstrate that complex chromatin architecture changes including enhancer switching are coordinated with ERα and PGR co-binding during diestrus and non-hormone receptor transcription factors such as HIF2A during estrus to regulate most differential gene expression across the estrous cycle.

The Personalized Inherited Signature Predisposing to Non-Small-Cell Lung Cancer in Non-Smokers

Lung cancer (LC) continues to be an important public health problem, being the most common form of cancer and a major cause of cancer deaths worldwide. Despite the great bulk of research to identify genetic susceptibility genes by genome-wide association studies, only few loci associated to nicotine dependence have been consistently replicated. Our previously published study in few phenotypically discordant sib-pairs identified a combination of germline truncating mutations in known cancer susceptibility genes in never-smoker early-onset LC patients, which does not present in their healthy sib. These results firstly demonstrated the presence of an oligogenic combination of disrupted cancer-predisposing genes in non-smokers patients, giving experimental support to a model of a "private genetic epidemiology". Here, we used a combination of whole-exome and RNA sequencing coupled with a discordant sib's model in a novel cohort of pairs of never-smokers early-onset LC patients and in their healthy sibs used as controls. We selected rare germline variants predicted as deleterious by CADD and SVM bioinformatics tools and absent in the healthy sib. Overall, we identified an average of 200 variants per patient, about 10 of which in cancer-predisposing genes. In most of them, RNA sequencing data reinforced the pathogenic role of the identified variants showing: (i) downregulation in LC tissue (indicating a "second hit" in tumor suppressor genes); (ii) upregulation in cancer tissue (likely oncogene); and (iii) downregulation in both normal and cancer tissue (indicating transcript instability). The combination of the two techniques demonstrates that each patient has an average of six (with a range from four to eight) private mutations with a functional effect in tumor-predisposing genes. The presence of a unique combination of disrupting events in the affected subjects may explain the absence of the familial clustering of non-small-cell lung cancer. In conclusion, these findings indicate that each patient has his/her own "predisposing signature" to cancer development and suggest the use of personalized therapeutic strategies in lung cancer.

Phenotypic and transcriptomics characterization uncovers genes underlying tuber yield traits and gene expression marker development in potato under aeroponics

MAIN CONCLUSION

Transcriptome analysis in potato varieties revealed genes associated with tuber yield-related traits and developed gene expression markers. This study aimed to identify genes involved in high tuber yield and its component traits in test potato varieties (Kufri Frysona, Kufri Khyati, and Kufri Mohan) compared to control (Kufri Sutlej). The aeroponic evaluation showed significant differences in yield-related traits in the varieties. Total RNA sequencing was performed using tuber and leaf tissues on the Illumina platform. The high-quality reads (QV > 25) mapping with the reference potato genomes revealed statistically significant (P < 0.05) differentially expressed genes (DEGs) into two categories: up-regulated (> 2 Log2 fold change) and down-regulated (< -2 Log2 fold change). DEGs were characterized by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Collectively, we identified genes participating in sugar metabolism, stress response, transcription factors, phytohormones, kinase proteins, and other genes greatly affecting tuber yield and its related traits. A few selected genes were UDP-glucose glucosyltransferase, glutathion S-transferase, GDSL esterase/lipase, transcription factors (MYB, WRKY, bHLH63, and BURP), phytohormones (auxin-induced protein X10A, and GA20 oxidase), kinase proteins (Kunitz-type tuber invertase inhibitor, BRASSINOSTEROID INSENSITIVE 1-associated receptor kinase 1) and laccase. Based on the selected 17 peptide sequences representing 13 genes, a phylogeny tree and motifs were analyzed. Real time-quantitative polymerase chain reaction (RT-qPCR) analysis was used to validate the RNA-seq results. RT-qPCR based gene expression markers were developed for the genes such as 101 kDa heat shock protein, catechol oxidase B chloroplastic, cysteine protease inhibitor 1, Kunitz-type tuber invertase inhibitor, and laccase to identify high yielding potato genotypes. Thus, our study paved the path for potential genes associated with tuber yield traits in potato under aeroponics.

Transcriptional Modulation during Photomorphogenesis in Rice Seedlings

Light is one of the most important factors regulating plant gene expression patterns, metabolism, physiology, growth, and development. To explore how light may induce or alter transcript splicing, we conducted RNA-Seq-based transcriptome analyses by comparing the samples harvested as etiolated seedlings grown under continuous dark conditions vs. the light-treated green seedlings. The study aims to reveal differentially regulated protein-coding genes and novel long noncoding RNAs (lncRNAs), their light-induced alternative splicing, and their association with biological pathways. We identified 14,766 differentially expressed genes, of which 4369 genes showed alternative splicing. We observed that genes mapped to the plastid-localized methyl-erythritol-phosphate (MEP) pathway were light-upregulated compared to the cytosolic mevalonate (MVA) pathway genes. Many of these genes also undergo splicing. These pathways provide crucial metabolite precursors for the biosynthesis of secondary metabolic compounds needed for chloroplast biogenesis, the establishment of a successful photosynthetic apparatus, and photomorphogenesis. In the chromosome-wide survey of the light-induced transcriptome, we observed intron retention as the most predominant splicing event. In addition, we identified 1709 novel lncRNA transcripts in our transcriptome data. This study provides insights on light-regulated gene expression and alternative splicing in rice.

The dopamine analogue CA140 alleviates AD pathology, neuroinflammation, and rescues synaptic/cognitive functions by modulating DRD1 signaling or directly binding to Abeta

BACKGROUND

We recently reported that the dopamine (DA) analogue CA140 modulates neuroinflammatory responses in lipopolysaccharide-injected wild-type (WT) mice and in 3-month-old 5xFAD mice, a model of Alzheimer's disease (AD). However, the effects of CA140 on Aβ/tau pathology and synaptic/cognitive function and its molecular mechanisms of action are unknown.

METHODS

To investigate the effects of CA140 on cognitive and synaptic function and AD pathology, 3-month-old WT mice or 8-month-old (aged) 5xFAD mice were injected with vehicle (10% DMSO) or CA140 (30 mg/kg, i.p.) daily for 10, 14, or 17 days. Behavioral tests, ELISA, electrophysiology, RNA sequencing, real-time PCR, Golgi staining, immunofluorescence staining, and western blotting were conducted.

RESULTS

In aged 5xFAD mice, a model of AD pathology, CA140 treatment significantly reduced Aβ/tau fibrillation, Aβ plaque number, tau hyperphosphorylation, and neuroinflammation by inhibiting NLRP3 activation. In addition, CA140 treatment downregulated the expression of cxcl10, a marker of AD-associated reactive astrocytes (RAs), and c1qa, a marker of the interaction of RAs with disease-associated microglia (DAMs) in 5xFAD mice. CA140 treatment also suppressed the mRNA levels of s100β and cxcl10, markers of AD-associated RAs, in primary astrocytes from 5xFAD mice. In primary microglial cells from 5xFAD mice, CA140 treatment increased the mRNA levels of markers of homeostatic microglia (cx3cr1 and p2ry12) and decreased the mRNA levels of a marker of proliferative region-associated microglia (gpnmb) and a marker of lipid-droplet-accumulating microglia (cln3). Importantly, CA140 treatment rescued scopolamine (SCO)-mediated deficits in long-term memory, dendritic spine number, and LTP impairment. In aged 5xFAD mice, these effects of CA140 treatment on cognitive/synaptic function and AD pathology were regulated by dopamine D1 receptor (DRD1)/Elk1 signaling. In primary hippocampal neurons and WT mice, CA140 treatment promoted long-term memory and dendritic spine formation via effects on DRD1/CaMKIIα and/or ERK signaling.

CONCLUSIONS

Our results indicate that CA140 improves neuronal/synaptic/cognitive function and ameliorates Aβ/tau pathology and neuroinflammation by modulating DRD1 signaling in primary hippocampal neurons, primary astrocytes/microglia, WT mice, and aged 5xFAD mice.

Telomere-to-telomere genome assembly of Eleocharis dulcis and expression profiles during corm development

Eleocharis dulcis (Burm. f.) Trin. ex Hensch., commonly known as Chinese water chestnut, is a traditional aquatic vegetable in China, and now is widely cultivated throughout the world because of its high nutritional value and unique tastes. Here, we report the assembly of a 493.24 Mb telomere-to-telomere (T2T) genome of E. dulcis accomplished by integrating ONT ultra-long reads, PacBio long reads and Hi-C data. The reference genome was anchored onto 111 gap-free chromosomes, containing 48.31% repeat elements and 33,493 predicted protein-coding genes. Whole genome duplication (WGD) and inter-genomic synteny analyses indicated that chromosome breakage and genome duplication in E. dulcis possibly occurred multiple times during genome evolution after its divergence from a common ancestor with Rhynchospora breviuscula at ca. 35.6 Mya. A comparative time-course transcriptome analysis of corm development revealed different patterns of gene expression between cultivated and wild accessions with the highest number of differentially expressed genes (DEGs, 15,870) at the middle swelling stage and some of the DEGs were significantly enriched for starch metabolic process.

Sporadic inclusion body myositis-derived myotube culture revealed muscle cell-autonomous expression profiles

Sporadic inclusion body myositis (sIBM) is a muscle disease in older people and is characterized by inflammatory cell invasion into intact muscle fibers and rimmed vacuoles. The pathomechanism of sIBM is not fully elucidated yet, and controversy exists as to whether sIBM is a primary autoimmune disease or a degenerative muscle disease with secondary inflammation. Previously, we established a method of collecting CD56-positive myoblasts from human skeletal muscle biopsy samples. We hypothesized that the myoblasts derived from these patients are useful to see the cell-autonomous pathomechanism of sIBM. With these resources, myoblasts were differentiated into myotubes, and the expression profiles of cell-autonomous pathology of sIBM were analyzed. Myoblasts from three sIBM cases and six controls were differentiated into myotubes. In the RNA-sequencing analysis of these "myotube" samples, 104 differentially expressed genes (DEGs) were found to be significantly upregulated by more than twofold in sIBM, and 13 DEGs were downregulated by less than twofold. For muscle biopsy samples, a comparative analysis was conducted to determine the extent to which "biopsy" and "myotube" samples differed. Fifty-three DEGs were extracted of which 32 (60%) had opposite directions of expression change (e.g., increased in biopsy vs decreased in myotube). Apolipoprotein E (apoE) and transmembrane protein 8C (TMEM8C or MYMK) were commonly upregulated in muscle biopsies and myotubes from sIBM. ApoE and myogenin protein levels were upregulated in sIBM. Given that enrichment analysis also captured changes in muscle contraction and development, the triggering of muscle atrophy signaling and abnormal muscle differentiation via MYMK or myogenin may be involved in the pathogenesis of sIBM. The presence of DEGs in sIBM suggests that the myotubes formed from sIBM-derived myoblasts revealed the existence of muscle cell-autonomous degeneration in sIBM. The catalog of DEGs will be an important resource for future studies on the pathogenesis of sIBM focusing on primary muscle degeneration.

Comparative genomics sheds new light on the convergent evolution of infrared vision in snakes

Infrared vision is a highly specialized sensory system that evolved independently in three clades of snakes. Apparently, convergent evolution occurred in the transient receptor potential ankyrin 1 (TRPA1) proteins of infrared-sensing snakes. However, this gene can only explain how infrared signals are received, and not the transduction and processing of those signals. We sequenced the genome of Xenopeltis unicolor, a key outgroup species of pythons, and performed a genome-wide analysis of convergence between two clades of infrared-sensing snakes. Our results revealed pervasive molecular adaptation in pathways associated with neural development and other functions, with parallel selection on loci associated with trigeminal nerve structural organization. In addition, we found evidence of convergent amino acid substitutions in a set of genes, including TRPA1 and TRPM2. The analysis also identified convergent accelerated evolution in non-coding elements near 12 genes involved in facial nerve structural organization and optic nerve development. Thus, convergent evolution occurred across multiple dimensions of infrared vision in vipers and pythons, as well as amino acid substitutions, non-coding elements, genes and functions. These changes enabled independent groups of snakes to develop and use infrared vision.

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

An ancient whole-genome duplication in barnacles contributes to their diversification and intertidal sessile life adaptation

INTRODUCTION

Whole-genome duplication (WGD) is one of the most sudden and dramatic events rarely reported in invertebrates, but its occurrence can lead to physiological, morphological, and behavioral diversification. WGD has also never been reported in barnacles, which is one of the most unique groups of crustaceans with extremely speciallized morphology (calcareous shells) and habits (intertidal sessile lifestyle).

OBJECTIVES

To investigate whether WGD has occurred in barnacles and examine its potential role in driving the adaptive evolution and diversification of barnacles.

METHODS

Based on a newly sequenced and assembled chromosome-level barnacle genome, a novel WGD event has been identified in barnacles through a comprehensive analysis of interchromosomal synteny, the Hox gene cluster, and synonymous substitution distribution.

RESULTS

We provide ample evidences for WGD in the barnacle genomes. Comparative genomic analysis indicates that this WGD event predates the divergence of Thoracicalcarea, occurring more than 247 million years ago. The retained ohnologs from the WGD are primarily enriched in various pathways related to environmental information processing, shedding light on the adaptive evolution and diversification of intertidal sessile lifestyle. In addition, transcriptomic analyses show that most of these ohnologs were differentially expressed following the ebb of tide. And the cytochrome P450 ohnologs with differential expression patterns are subject to subfunctionalization and/or neofunctionalization for intertidal adaptation. Besides WGD, parallel evolution underlying intertidal adaptation has also occurred in barnacles.

CONCLUSION

This study revealed an ancient WGD event in the barnacle genomes, which is potentially associated with the origin and diversification of thoracican barnacles, and may have contributed to the adaptive evolution of their intertidal sessile lifestyle.

An autoregulatory poison exon in Smndc1 is conserved across kingdoms and influences organism growth

Many of the most highly conserved elements in the human genome are "poison exons," alternatively spliced exons that contain premature termination codons and permit post-transcriptional regulation of mRNA abundance through induction of nonsense-mediated mRNA decay (NMD). Poison exons are widely assumed to be highly conserved due to their presumed importance for organismal fitness, but this functional importance has never been tested in the context of a whole organism. Here, we report that a poison exon in Smndc1 is conserved across mammals and plants and plays a molecular autoregulatory function in both kingdoms. We generated mouse and A. thaliana models lacking this poison exon to find its loss leads to deregulation of SMNDC1 protein levels, pervasive alterations in mRNA processing, and organismal size restriction. Together, these models demonstrate the importance of poison exons for both molecular and organismal phenotypes that likely explain their extraordinary conservation.

Angiopoietin-like growth factor-derived peptides as biological activators of adipose-derived mesenchymal stromal cells

Adipose-derived mesenchymal stromal cells (AD-MSCs) are an essential issue in modern medicine. Extensive preclinical and clinical studies have shown that mesenchymal stromal/stem cells, including AD-MSCs, have specific properties (ability to differentiate into other cells, recruitment to the site of injury) of particular importance in the regenerative process. Ongoing research aims to elucidate factors supporting AD-MSC culture and differentiation in vitro. Angiopoietin-like proteins (ANGPTLs), known for their pleiotropic effects in lipid and glucose metabolism, may play a significant role in this context. Regeneration is a complex and dynamic process controlled by many factors. ANGPTL6 (Angiopoietin-related growth factor, AGF), among many activities modulated the biological activity of stem cells. This study examined the influence of synthesized AGF-derived peptides, designated as AGF9 and AGF27, on AD-MSCs. AGF9 and AGF27 enhanced the viability and migration of AD-MSCs and acted as a chemotactic factor for these cells. AGF9 stimulated chondrogenesis and lipid synthesis during AD-MSCs differentiation, influenced AD-MSCs cytokine secretion and modulated transcriptome for such basic cell activities as migration, transport of molecules, and apoptosis. The ability of AGF9 to modulate the biological activity of AD-MSCs warrants the consideration of this peptide a noteworthy therapeutic agent that deserves further investigation for applications in regenerative medicine.

Transcription factors HB21/40/53 trigger inflorescence arrest through abscisic acid accumulation at the end of flowering

Flowers, and hence, fruits and seeds, are produced by the activity of the inflorescence meristem after the floral transition. In plants with indeterminate inflorescences, the final number of flowers produced by the inflorescence meristem is determined by the length of the flowering period, which ends with inflorescence arrest. Inflorescence arrest depends on many different factors, such as the presence of seeds, the influence of the environment, or endogenous factors such as phytohormone levels and age, which modulate inflorescence meristem activity. The FRUITFULL-APETALA2 (FUL-AP2) pathway plays a major role in regulating the end of flowering, likely integrating both endogenous cues and those related to seed formation. Among AP2 targets, HOMEOBOX PROTEIN21 (HB21) has been identified as a putative mediator of AP2 function in the control of inflorescence arrest. HB21 is a homeodomain leucine zipper transcription factor involved in establishing axillary bud dormancy. Here, we characterized the role of HB21 in the control of the inflorescence arrest at the end of flowering in Arabidopsis (Arabidopsis thaliana). HB21, together with HB40 and HB53, are upregulated in the inflorescence apex at the end of flowering, promoting floral bud arrest. We also show that abscisic acid (ABA) accumulation occurs in the inflorescence apex in an HB-dependent manner. Our work suggests a physiological role of ABA in floral bud arrest at the end of flowering, pointing to ABA as a regulator of inflorescence arrest downstream of the HB21/40/53 genes.

The first complete mitochondrial genome of Grossulariaceae: Molecular features, structure recombination, and genetic evolution

BACKGROUND

Mitochondria play crucial roles in the growth, development, and adaptation of plants. Blackcurrant (Ribes nigrum L.) stands out as a significant berry species due to its rich nutritional profile, medicinal properties, and health benefits. Despite its importance, the mitochondrial genome of blackcurrant remains unassembled.

RESULTS

This study presents the first assembly of the mitochondrial genome of R. nigrum in the Grossulariaceae family. The genome spans 450,227 base pairs (bp) and encompasses 39 protein-coding genes (PCGs), 19 transfer RNAs (tRNAs), and three ribosomal RNAs (rRNAs). Protein-coding regions constitute 8.88% of the entire genome. Additionally, we identified 180 simple sequence repeats, 12 tandem repeats, and 432 pairs of dispersed repeats. Notably, the dispersed sequence R1 (cotig3, 1,129 bp) mediated genome recombination, resulting in the formation of two major conformations, namely master and double circles. Furthermore, we identified 731 C-to-U RNA editing sites within the PCGs. Among these, cox1-2, nad1-2, and nad4L-2 were associated with the creation of start codons, whereas atp6-718 and rps10-391 were linked to termination codons. We also detected fourteen plastome fragments within the mitogenome, constituting 1.11% of the total length. Phylogenetic analysis suggests that R. nigrum might have undergone multiple genomic reorganization and/or gene transfer events, resulting in the loss of two PCGs (rps2 and rps11) during its evolutionary history.

CONCLUSIONS

This investigation unveils the molecular characteristics of the R. nigrum mitogenome, shedding light on its evolutionary trajectory and phylogenetic implications. Furthermore, it serves as a valuable reference for evolutionary research and germplasm identification within the genus.

Identification and core gene-mining of Weighted Gene Co-expression Network Analysis-based co-expression modules related to flood resistance in quinoa seedlings

BACKGROUND

As an emerging food crop with high nutritional value, quinoa has been favored by consumers in recent years; however, flooding, as an abiotic stress, seriously affects its growth and development. Currently, reports on the molecular mechanisms related to quinoa waterlogging stress responses are lacking; accordingly, the core genes related to these processes were explored via Weighted Gene Co-expression Network Analysis (WGCNA).

RESULTS

Based on the transcriptome data, WGCNA was used to construct a co-expression network of weighted genes associated with flooding resistance-associated physiological traits and metabolites. Here, 16 closely related co-expression modules were obtained, and 10 core genes with the highest association with the target traits were mined from the two modules. Functional annotations revealed the biological processes and metabolic pathways involved in waterlogging stress, and four candidates related to flooding resistance, specifically AP2/ERF, MYB, bHLH, and WRKY-family TFs, were also identified.

CONCLUSIONS

These results provide clues to the identification of core genes for quinoa underlying quinoa waterlogging stress responses. This could ultimately provide a theoretical foundation for breeding new quinoa varieties with flooding tolerance.

Integrated analysis of transcriptomics and defense-related phytohormones to discover hub genes conferring maize Gibberella ear rot caused by Fusarium Graminearum

BACKGROUND

Gibberella ear rot (GER) is one of the most devastating diseases in maize growing areas, which directly reduces grain yield and quality. However, the underlying defense response of maize to pathogens infection is largely unknown.

RESULTS

To gain a comprehensive understanding of the defense response in GER resistance, two contrasting inbred lines 'Nov-82' and 'H10' were used to explore transcriptomic profiles and defense-related phytohormonal alterations during Fusarium graminearum infection. Transcriptomic analysis revealed 4,417 and 4,313 differentially expressed genes (DEGs) from the Nov-82 and H10, respectively, and 647 common DEGs between the two lines. More DEGs were obviously enriched in phenylpropanoid biosynthesis, secondary metabolites biosynthesis, metabolic process and defense-related pathways. In addition, the concentration of the defense-related phytohormones, jasmonates (JAs) and salicylates (SAs), was greatly induced after the pathogen infection. The level of JAs in H10 was more higher than in Nov-82, whereas an opposite pattern for the SA between the both lines. Integrated analysis of the DEGs and the phytohormones revealed five vital modules based on co-expression network analysis according to their correlation. A total of 12 hub genes encoding fatty acid desaturase, subtilisin-like protease, ethylene-responsive transcription factor, 1-aminocyclopropane-1-carboxylate oxidase, and sugar transport protein were captured from the key modules, indicating that these genes might play unique roles in response to pathogen infection, CONCLUSIONS: Overall, our results indicate that large number DEGs related to plant disease resistance and different alteration of defensive phytohormones were activated during F. graminearum infection, providing new insight into the defense response against pathogen invasion, in addition to the identified hub genes that can be further investigated for enhancing maize GER resistance.

Vestigial like 4 regulates the adipogenesis of classical brown adipose tissue

Brown adipose tissue (BAT) is mammals' primary non-shivering thermogenesis organ, and the molecular mechanisms regulating BAT growth and adipogenesis are largely unknown. The Hippo-YAP pathway has been well-known for controlling organ size, and Vestigial like 4 (VGLL4) is a transcriptional regulator that modulates the Hippo-YAP pathway by competing against YAP for binding to TEAD proteins. In this study, we dissected the function of VGLL4 in regulating BAT development. We generated a conventional Vgll4 mutant mouse line, in which the two Tondu (TDU) domains of VGLL4 were disrupted. We found that deletion of the TDU domains of VGLL4 resulted in perinatal lethality and paucity of the interscapular BAT. Histological and magnetic resonance imaging studies confirmed that the adipogenesis of BAT was impaired in Vgll4 mutants. Adeno-associated virus (AAV) mediated, brown adipocyte-specific overexpression of VGLL4 increased BAT volume and protected the adult male mice from acute cold stress. Genomic studies suggest that VGLL4/TEAD1 complex directly regulates the myogenic and adipogenic gene expression programs of BAT. In conclusion, our data identify VGLL4 as a previously unrecognized adipogenesis factor that regulates classical BAT development.

Early Developmental Characteristics and Features of a Three-Dimensional Retinal Organoid Model of X-Linked Juvenile Retinoschisis

X-linked juvenile retinoschisis (XLRS) is a hereditary retinal degeneration affecting young males caused by mutations in the retinoschisin (RS1) gene. We generated human induced pluripotent stem cells (hiPSCs) from XLRS patients and established three-dimensional retinal organoids (ROs) for disease investigation. This disease model recapitulates the characteristics of XLRS, exhibiting defects in RS1 protein production and photoreceptor cell development. XLRS ROs also revealed dysregulation of Na/K-ATPase due to RS1 deficiency and increased ERK signaling pathway activity. Transcriptomic analyses of XLRS ROs showed decreased expression of retinal cells, particularly photoreceptor cells. Furthermore, relevant recovery of the XLRS phenotype was observed when co-cultured with control ROs derived from healthy subject during the early stages of differentiation. In conclusion, our in vitro XLRS RO model presents a valuable tool for elucidating the pathophysiological mechanisms underlying XLRS, offering insights into disease progression. Additionally, this model serves as a robust platform for the development and optimization of targeted therapeutic strategies, potentially improving treatment outcomes for patients with XLRS.

The single RRM domain-containing protein SARP1 is required for establishment of the separation zone in Arabidopsis

Abscission is the shedding of plant organs in response to developmental and environmental cues. Abscission involves cell separation between two neighboring cell types, residuum cells (RECs) and secession cells (SECs) in the floral abscission zone (AZ) in Arabidopsis thaliana. However, the regulatory mechanisms behind the spatial determination that governs cell separation are largely unknown. The class I KNOTTED-like homeobox (KNOX) transcription factor BREVIPEDICELLUS (BP) negatively regulates AZ cell size and number in Arabidopsis. To identify new players participating in abscission, we performed a genetic screen by activation tagging a weak complementation line of bp-3. We identified the mutant ebp1 (enhancer of BP1) displaying delayed floral organ abscission. The ebp1 mutant showed a concaved surface in SECs and abnormally stacked cells on the top of RECs, in contrast to the precisely separated surface in the wild-type. Molecular and histological analyses revealed that the transcriptional programming during cell differentiation in the AZ is compromised in ebp1. The SECs of ebp1 have acquired REC-like properties, including cuticle formation and superoxide production. We show that SEPARATION AFFECTING RNA-BINDING PROTEIN1 (SARP1) is upregulated in ebp1 and plays a role in the establishment of the cell separation layer during floral organ abscission in Arabidopsis.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Transcriptome Analysis Reveals the Mechanism of Quinoa Polysaccharides Inhibiting 3T3-L1 Preadipocyte Proliferation

Quinoa is a highly nutritious and biologically active crop. Prior studies have demonstrated that quinoa polysaccharides exhibit anti-obesity activity. This investigation confirmed that quinoa polysaccharides have the ability to inhibit the growth of 3T3-L1 preadipocytes. The objective of transcriptome research was to investigate the mechanism of quinoa water-extracted polysaccharides and quinoa alkaline-extracted polysaccharides that hinder the growth of 3T3-L1 preadipocytes. There were 2194 genes that showed differential expression between untreated cells and those treated with high concentrations of quinoa water-extracted polysaccharides (QWPHs). There were 1774 genes that showed differential expression between untreated cells and those treated with high concentrations of quinoa alkaline-extracted polysaccharides (QAPHs). Through gene ontology and KEGG pathway analysis, 20 characteristic pathways are found significantly enriched between the untreated group and the QAPH and QWPH groups. These pathways include the NOD-like receptor, Hepatitis C, and the PI3K-Akt signaling pathway. Atp13A4 and Gbgt1 have been identified as genes that are upregulated and downregulated in both the untreated group and the QWPH group, as well as in the untreated group and the QAPH group. These findings establish a theoretical foundation for exploring quinoa polysaccharides as an anti-obesity agent.

Establishing a model composed of immune-related gene-modules to predict tumor immunotherapy response

At present, tumor immunotherapy has been widely applied to treat various cancers. However, the accuracy of predicting treatment efficacy has not yet achieved a significant breakthrough. This study aimed to construct a prediction model based on the modified WGCNA algorithm to precisely judge the anti-tumor immune response. First, we used a murine colon cancer model to screen corresponding DEGs according to different groups. GSEA was used to analyze the potential mechanisms of the immune-related DEGs (irDEGs) in each group. Subsequently, the intersection of the irDEGs in every group was acquired, and 7 gene-modules were mapped. Finally, 4 gene-modules including cogenes, antiPD-1 immu-genes, chemo immu-genes and comb immu-genes, were selected for subsequent study. Furthermore, a clinical dataset of gastric cancer patients receiving immunotherapy was enrolled, and the irDEGs were identified. A total of 34 vital irDEGs were obtained from the intersections of the vital irDEGs and the four gene-modules. Next, the vital irDEGs were analyzed by the modified WGCNA algorithm, and the correlation coefficients between the 4 gene-modules and the response status to immunotherapy were calculated. Thus, a prediction model based on correlation coefficients was built, and the corresponding model scores were acquired. The AUC calculated according to the model score was 0.727, which was non-inferior to that of the ESTIMATE score and the TIDE score. Meanwhile, the AUC calculated according to the classification of the model scores was 0.705, which was non-inferior to that of the ESTIMATE classification and the TIDE classification. The prediction accuracy of the model was validated in clinical datasets of other cancers.

Using clusterProfiler to characterize multiomics data

With the advent of multiomics, software capable of multidimensional enrichment analysis has become increasingly crucial for uncovering gene set variations in biological processes and disease pathways. This is essential for elucidating disease mechanisms and identifying potential therapeutic targets. clusterProfiler stands out for its comprehensive utilization of databases and advanced visualization features. Importantly, clusterProfiler supports various biological knowledge, including Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, through performing over-representation and gene set enrichment analyses. A key feature is that clusterProfiler allows users to choose from various graphical outputs to visualize results, enhancing interpretability. This protocol describes innovative ways in which clusterProfiler has been used for integrating metabolomics and metagenomics analyses, identifying and characterizing transcription factors under stress conditions, and annotating cells in single-cell studies. In all cases, the computational steps can be completed within ~2 min. clusterProfiler is released through the Bioconductor project and can be accessed via https://bioconductor.org/packages/clusterProfiler/ .

RBM10 loss induces aberrant splicing of cytoskeletal and extracellular matrix mRNAs and promotes metastatic fitness

RBM10 modulates transcriptome-wide cassette exon splicing. Loss-of-function RBM10 mutations are enriched in thyroid cancers with distant metastases. Analysis of transcriptomes and genes mis-spliced by RBM10 loss showed pro-migratory and RHO/RAC signaling signatures. RBM10 loss increases cell velocity. Cytoskeletal and ECM transcripts subject to exon-inclusion events included vinculin (VCL), tenascin C (TNC) and CD44. Knockdown of the VCL exon inclusion transcript in RBM10-null cells reduced cell velocity, whereas knockdown of TNC and CD44 exon-inclusion isoforms reduced invasiveness. RAC1-GTP levels were increased in RBM10-null cells. Mouse Hras G12V /Rbm1O KO thyrocytes develop metastases that are reversed by RBM10 or by combined knockdown of VCL, CD44 and TNC inclusion isoforms. Thus, RBM10 loss generates exon inclusions in transcripts regulating ECM-cytoskeletal interactions, leading to RAC1 activation and metastatic competency. Moreover, a CRISPR-Cas9 screen for synthetic lethality with RBM10 loss identified NFkB effectors as central to viability, providing a therapeutic target for these lethal thyroid cancers.

Transcriptome analysis reveals genes associated with late blight resistance in potato

Late blight is a serious disease of potato worldwide. Our study aimed to unveil genes involved in late blight resistance in potato by RNA-seq analysis after artificial inoculation under controlled conditions. In this study, two potato somatic hybrids (P7 and Crd6) and three varieties such as Kufri Girdhari, Kufri Jyoti and Kufri Bahar (control) were used. Transcriptiome analysis revealed statistically significant (p < 0.05) differentially expressed genes (DEGs), which were analysed into up-regulated and down-regulated genes. Further, DEGs were functionally characterized by the Gene Ontology annotations and the Kyoto Encyclopedia of Genes and Genomes pathways. Overall, some of the up-regulated genes in resistant genotypes were disease resistance proteins such as CC-NBS-LRR resistance protein, ankyrin repeat family protein, cytochrome P450, leucine-rich repeat family protein/protein kinase family, and MYB transcription factor. Sequence diversity analysis based on 38 peptide sequences representing 18 genes showed distinct variation and the presence of three motifs in 15 amino acid sequences. Selected genes were also validated by real-time quantitative polymerase chain reaction analysis. Interestingly, gene expression markers were developed for late blight resistant genotypes. Our study elucidates genes involved in imparting late blight resistance in potato, which will be beneficial for its management strategies in the future.

Secreted clusterin inhibits tumorigenesis by modulating tumor cells and macrophages in human meningioma

BACKGROUND

Meningioma is the most common primary intracranial tumor with a high frequency of postoperative recurrence, yet the biology of the meningioma malignancy process is still obscure.

METHODS

To identify potential therapeutic targets and tumor suppressors, we performed single-cell transcriptome analysis through meningioma malignancy, which included 18 samples spanning normal meninges, benign and high-grade in situ tumors, and lung metastases, for extensive transcriptome characterization. Tumor suppressor candidate gene and molecular mechanism were functionally validated at the animal model and cellular levels.

RESULTS

Comprehensive analysis and validation in mice and clinical cohorts indicated clusterin (CLU) had suppressive function for meningioma tumorigenesis and malignancy by inducing mitochondria damage and triggering type 1 interferon pathway dependent on its secreted isoform, and the inhibition effect was enhanced by TNFα as TNFα also induced type 1 interferon pathway. Meanwhile, both intra- and extracellular CLU overexpression enhanced macrophage polarization towards M1 phenotype and TNFα production, thus promoting tumor killing and phagocytosis.

CONCLUSIONS

CLU might be a key brake of meningioma malignance by synchronously modulating tumor cells and their microenvironment. Our work provides comprehensive insights into meningioma malignancy and a potential therapeutic strategy.

Abscisic Acid Induces DNA Methylation Alteration in Genes Related to Berry Ripening and Stress Response in Grape (Vitis vinifera L)

Abscisic acid (ABA) is a major regulator of nonclimacteric fruit ripening, with its processes involving epigenetic mechanisms. It remains unclear whether DNA methylation is associated with ABA-regulated ripening. In this study, we investigated the patterns of DNA methylation and gene expression following ABA treatment in grape berries by using whole-genome bisulfite sequencing and RNA-sequencing. ABA application changed global DNA methylation in grapes. The hyper-/hypo-differently methylated regions were enriched in defense-related metabolism, degreening processes, or ripening-related metabolic pathways. Many differentially expressed genes showed an alteration in DNA methylation after ABA treatment. Specifically, ten downregulated genes with hypermethylation in promoters were involved in the ripening process, ABA homeostasis/signaling, and stress response. Nine upregulated genes exhibiting hypo-methylation in promoters were related to the ripening process and stress response. These findings demonstrated ABA-induced DNA alteration of ripening related and stress-responsive genes during grape ripening, which provides new insights of the epigenetic regulation of ABA on fruit ripening.

Identifying Key Genes Involved in Axillary Lymph Node Metastasis in Breast Cancer Using Advanced RNA-Seq Analysis: A Methodological Approach with GLMQL and MAS

Our study aims to address the methodological challenges frequently encountered in RNA-Seq data analysis within cancer studies. Specifically, it enhances the identification of key genes involved in axillary lymph node metastasis (ALNM) in breast cancer. We employ Generalized Linear Models with Quasi-Likelihood (GLMQLs) to manage the inherently discrete and overdispersed nature of RNA-Seq data, marking a significant improvement over conventional methods such as the t-test, which assumes a normal distribution and equal variances across samples. We utilize the Trimmed Mean of M-values (TMMs) method for normalization to address library-specific compositional differences effectively. Our study focuses on a distinct cohort of 104 untreated patients from the TCGA Breast Invasive Carcinoma (BRCA) dataset to maintain an untainted genetic profile, thereby providing more accurate insights into the genetic underpinnings of lymph node metastasis. This strategic selection paves the way for developing early intervention strategies and targeted therapies. Our analysis is exclusively dedicated to protein-coding genes, enriched by the Magnitude Altitude Scoring (MAS) system, which rigorously identifies key genes that could serve as predictors in developing an ALNM predictive model. Our novel approach has pinpointed several genes significantly linked to ALNM in breast cancer, offering vital insights into the molecular dynamics of cancer development and metastasis. These genes, including ERBB2, CCNA1, FOXC2, LEFTY2, VTN, ACKR3, and PTGS2, are involved in key processes like apoptosis, epithelial-mesenchymal transition, angiogenesis, response to hypoxia, and KRAS signaling pathways, which are crucial for tumor virulence and the spread of metastases. Moreover, the approach has also emphasized the importance of the small proline-rich protein family (SPRR), including SPRR2B, SPRR2E, and SPRR2D, recognized for their significant involvement in cancer-related pathways and their potential as therapeutic targets. Important transcripts such as H3C10, H1-2, PADI4, and others have been highlighted as critical in modulating the chromatin structure and gene expression, fundamental for the progression and spread of cancer.

Microbe-induced coordination of plant iron-sulfur metabolism enhances high-light-stress tolerance of Arabidopsis

High-light stress strongly limits agricultural production in subtropical and tropical regions owing to photo-oxidative damage, decreased growth, and decreased yield. Here, we investigated whether beneficial microbes can protect plants under high-light stress. We found that Enterobacter sp. SA187 (SA187) supports the growth of Arabidopsis thaliana under high-light stress by reducing the accumulation of reactive oxygen species and maintaining photosynthesis. Under high-light stress, SA187 triggers dynamic changes in the expression of Arabidopsis genes related to fortified iron metabolism and redox regulation, thereby enhancing the antioxidative glutathione/glutaredoxin redox system of the plant. Genetic analysis showed that the enhancement of iron and sulfur metabolism by SA187 is coordinated by ethylene signaling. In summary, beneficial microbes could be an effective and inexpensive means of enhancing high-light-stress tolerance in plants.

MAST4 regulates stem cell maintenance with DLX3 for epithelial development and amelogenesis

The asymmetric division of stem cells permits the maintenance of the cell population and differentiation for harmonious progress. Developing mouse incisors allows inspection of the role of the stem cell niche to provide specific insights into essential developmental phases. Microtubule-associated serine/threonine kinase family member 4 (Mast4) knockout (KO) mice showed abnormal incisor development with low hardness, as the size of the apical bud was decreased and preameloblasts were shifted to the apical side, resulting in amelogenesis imperfecta. In addition, Mast4 KO incisors showed abnormal enamel maturation, and stem cell maintenance was inhibited as amelogenesis was accelerated with Wnt signal downregulation. Distal-Less Homeobox 3 (DLX3), a critical factor in tooth amelogenesis, is considered to be responsible for the development of amelogenesis imperfecta in humans. MAST4 directly binds to DLX3 and induces phosphorylation at three residues within the nuclear localization site (NLS) that promotes the nuclear translocation of DLX3. MAST4-mediated phosphorylation of DLX3 ultimately controls the transcription of DLX3 target genes, which are carbonic anhydrase and ion transporter genes involved in the pH regulation process during ameloblast maturation. Taken together, our data reveal a novel role for MAST4 as a critical regulator of the entire amelogenesis process through its control of Wnt signaling and DLX3 transcriptional activity.

OsEIL2 balances rice immune responses against (hemi)biotrophic and necrotrophic pathogens via the salicylic acid and jasmonic acid synergism

Plants typically activate distinct defense pathways against various pathogens. Heightened resistance to one pathogen often coincides with increased susceptibility to another pathogen. However, the underlying molecular basis of this antagonistic response remains unclear. Here, we demonstrate that mutants defective in the transcription factor ETHYLENE-INSENSITIVE 3-LIKE 2 (OsEIL2) exhibited enhanced resistance to the biotrophic bacterial pathogen Xanthomonas oryzae pv oryzae and to the hemibiotrophic fungal pathogen Magnaporthe oryzae, but enhanced susceptibility to the necrotrophic fungal pathogen Rhizoctonia solani. Furthermore, necrotroph-induced OsEIL2 binds to the promoter of OsWRKY67 with high affinity, leading to the upregulation of salicylic acid (SA)/jasmonic acid (JA) pathway genes and increased SA/JA levels, ultimately resulting in enhanced resistance. However, biotroph- and hemibiotroph-induced OsEIL2 targets OsERF083, resulting in the inhibition of SA/JA pathway genes and decreased SA/JA levels, ultimately leading to reduced resistance. Our findings unveil a previously uncharacterized defense mechanism wherein two distinct transcriptional regulatory modules differentially mediate immunity against pathogens with different lifestyles through the transcriptional reprogramming of phytohormone pathway genes.

Nuclear receptor Nur77 and Yin-Yang 1 synergistically increase mitochondrial abundance and activity in macrophages

Mitochondrial biogenesis requires precise regulation of both mitochondrial-encoded and nuclear-encoded genes. Nuclear receptor Nur77 is known to regulate mitochondrial metabolism in macrophages and skeletal muscle. Here, we compared genome-wide Nur77 binding site and target gene expression in these two cell types, which revealed conserved regulation of mitochondrial genes and enrichment of motifs for the transcription factor Yin-Yang 1 (YY1). We show that Nur77 and YY1 interact, that YY1 increases Nur77 activity, and that their binding sites are co-enriched at mitochondrial ribosomal protein gene loci in macrophages. Nur77 and YY1 co-expression synergistically increases Mrpl1 expression as well as mitochondrial abundance and activity in macrophages but not skeletal muscle. As such, we identify a macrophage-specific Nur77-YY1 interaction that enhances mitochondrial metabolism.

Temporal gene expression during asexual development of the apicomplexan Sarcocystis neurona

Asexual replication in the apicomplexan Sarcocystis neurona involves two main developmental stages: the motile extracellular merozoite and the sessile intracellular schizont. Merozoites invade host cells and transform into schizonts that undergo replication via endopolygeny to form multiple (64) daughter merozoites that are invasive to new host cells. Given that the capabilities of the merozoite vary significantly from the schizont, the patterns of transcript levels throughout the asexual lifecycle were determined and compared in this study. RNA-Seq data were generated from extracellular merozoites and four intracellular schizont development time points. Of the 6,938 genes annotated in the S. neurona genome, 6,784 were identified in the transcriptome. Of these, 4,111 genes exhibited significant differential expression between the merozoite and at least one schizont development time point. Transcript levels were significantly higher for 2,338 genes in the merozoite and 1,773 genes in the schizont stages. Included in this list were genes encoding the secretory pathogenesis determinants (SPDs), which encompass the surface antigen and SAG-related sequence (SAG/SRS) and the secretory organelle proteins of the invasive zoite stage (micronemes, rhoptries, and dense granules). As anticipated, many of the S. neurona SPD gene transcripts were abundant in merozoites. However, several SPD transcripts were elevated in intracellular schizonts, suggesting roles unrelated to host cell invasion and the initial establishment of the intracellular niche. The hypothetical genes that are potentially unique to the genus Sarcocystis are of particular interest. Their conserved expression patterns are instructive for future investigations into the possible functions of these putative Sarcocystis-unique genes.

IMPORTANCE

The genus Sarcocystis is an expansive clade within the Apicomplexa, with the species S. neurona being an important cause of neurological disease in horses. Research to decipher the biology of S. neurona and its host-pathogen interactions can be enhanced by gene expression data. This study has identified conserved apicomplexan orthologs in S. neurona, putative Sarcocystis-unique genes, and gene transcripts abundant in the merozoite and schizont stages. Importantly, we have identified distinct clusters of genes with transcript levels peaking during different intracellular schizont development time points, reflecting active gene expression changes across endopolygeny. Each cluster also has subsets of transcripts with unknown functions, and investigation of these seemingly Sarcocystis-unique transcripts will provide insights into the interesting biology of this parasite genus.

Transcriptomic Analysis of Pubertal and Adult Virgin Mouse Mammary Epithelial and Stromal Cell Populations

Conflicting data exist as to how mammary epithelial cell proliferation changes during the reproductive cycle. To study the effect of endogenous hormone fluctuations on gene expression in the mouse mammary gland, we performed bulk RNAseq analyses of epithelial and stromal cell populations that were isolated either during puberty or at different stages of the adult virgin estrous cycle. Our data confirm prior findings that proliferative changes do not occur in every mouse in every cycle. We also show that during the estrous cycle the main gene expression changes occur in adipocytes and fibroblasts. Finally, we present a comprehensive overview of the Wnt gene expression landscape in different mammary gland cell types in pubertal and adult mice. This work contributes to understanding the effects of physiological hormone fluctuations and locally produced signaling molecules on gene expression changes in the mammary gland during the reproductive cycle and should be a useful resource for future studies investigating gene expression patterns in different cell types across different developmental timepoints.

A comprehensive workflow for optimizing RNA-seq data analysis

BACKGROUND

Current RNA-seq analysis software for RNA-seq data tends to use similar parameters across different species without considering species-specific differences. However, the suitability and accuracy of these tools may vary when analyzing data from different species, such as humans, animals, plants, fungi, and bacteria. For most laboratory researchers lacking a background in information science, determining how to construct an analysis workflow that meets their specific needs from the array of complex analytical tools available poses a significant challenge.

RESULTS

By utilizing RNA-seq data from plants, animals, and fungi, it was observed that different analytical tools demonstrate some variations in performance when applied to different species. A comprehensive experiment was conducted specifically for analyzing plant pathogenic fungal data, focusing on differential gene analysis as the ultimate goal. In this study, 288 pipelines using different tools were applied to analyze five fungal RNA-seq datasets, and the performance of their results was evaluated based on simulation. This led to the establishment of a relatively universal and superior fungal RNA-seq analysis pipeline that can serve as a reference, and certain standards for selecting analysis tools were derived for reference. Additionally, we compared various tools for alternative splicing analysis. The results based on simulated data indicated that rMATS remained the optimal choice, although consideration could be given to supplementing with tools such as SpliceWiz.

CONCLUSION

The experimental results demonstrate that, in comparison to the default software parameter configurations, the analysis combination results after tuning can provide more accurate biological insights. It is beneficial to carefully select suitable analysis software based on the data, rather than indiscriminately choosing tools, in order to achieve high-quality analysis results more efficiently.

Chromosome-scale genome assemblies of sexually dimorphic male and female Acrossocheilus fasciatus

Acrossocheilus fasciatus is a stream-dwelling fish species of the Barbinae subfamily. It is valued for its colorfully striped appearance and delicious meat. This species is also characterized by apparent sexual dimorphism and toxic ovum. Biology and aquaculture researches of A. fasciatus are hindered by the lack of a high-quality reference genome. Here, we report chromosome-level genome assemblies of the male and female A. fasciatus. The HiFi-only genome assemblies for both female and male individuals were 899.13 Mb (N50 length of 32.58 Mb) and 885.68 Mb (N50 length of 33.06 Mb), respectively. Notably, a substantial proportion of the assembled sequences, accounting for 96.15% and 98.35% for female and male genomes, respectively, were successfully anchored onto 25 chromosomes utilizing Hi-C data. We annotated the female assembly as a reference genome and identified a total of 400.62 Mb (44.56%) repetitive sequences, 27,392 protein-coding genes, and 35,869 ncRNAs. The high-quality male and female reference genomes will provide genomic resources for developing sex-specific molecular markers, inform single-sex breeding, and elucidate genetic mechanisms of sexual dimorphism.

Machine Learning-Based Etiologic Subtyping of Ischemic Stroke Using Circulating Exosomal microRNAs

Ischemic stroke is a major cause of mortality worldwide. Proper etiological subtyping of ischemic stroke is crucial for tailoring treatment strategies. This study explored the utility of circulating microRNAs encapsulated in extracellular vesicles (EV-miRNAs) to distinguish the following ischemic stroke subtypes: large artery atherosclerosis (LAA), cardioembolic stroke (CES), and small artery occlusion (SAO). Using next-generation sequencing (NGS) and machine-learning techniques, we identified differentially expressed miRNAs (DEMs) associated with each subtype. Through patient selection and diagnostic evaluation, a cohort of 70 patients with acute ischemic stroke was classified: 24 in the LAA group, 24 in the SAO group, and 22 in the CES group. Our findings revealed distinct EV-miRNA profiles among the groups, suggesting their potential as diagnostic markers. Machine-learning models, particularly logistic regression models, exhibited a high diagnostic accuracy of 92% for subtype discrimination. The collective influence of multiple miRNAs was more crucial than that of individual miRNAs. Additionally, bioinformatics analyses have elucidated the functional implications of DEMs in stroke pathophysiology, offering insights into the underlying mechanisms. Despite limitations like sample size constraints and retrospective design, our study underscores the promise of EV-miRNAs coupled with machine learning for ischemic stroke subtype classification. Further investigations are warranted to validate the clinical utility of the identified EV-miRNA biomarkers in stroke patients.

The role of cis-zeatin in enhancing high-temperature resistance and fucoxanthin biosynthesis in Phaeodactylum tricornutum

Phaeodactylum tricornutum a prominent source of industrial fucoxanthin production, faces challenges in its application due to its tolerance to high-temperature environments. This study investigates the physiological responses of P. tricornutum to high-temperature stress and its impact on fucoxanthin content, with a specific focus on the role of cis-zeatin. The results reveal that high-temperature stress inhibits P. tricornutum's growth and photosynthetic activity, leading to a decrease in fucoxanthin content. Transcriptome analysis shows that high temperature suppresses the expression of genes related to photosynthesis (e.g., psbO, psbQ, and OEC) and fucoxanthin biosynthesis (e.g., PYS, PDS1, and PSD2), underscoring the negative effects of high temperature on P. tricornutum. Interestingly, genes associated with cis-zeatin biosynthesis and cytokinesis signaling pathways exhibited increased expression under high-temperature conditions, indicating a potential role of cis-zeatin signaling in response to elevated temperatures. Content measurements confirm that high temperature enhances cis-zeatin content. Furthermore, the exogenous addition of cytokinesis mimetics or inhibitors significantly affected P. tricornutum's high-temperature resistance. Overexpression of the cis-zeatin biosynthetic enzyme gene tRNA DMATase enhanced P. tricornutum's resistance to high-temperature stress, while genetic knockout of tRNA DMATase reduced its resistance to high temperatures. Therefore, this research not only uncovers a novel mechanism for high-temperature resistance in P. tricornutum but also offers a possible alga species that can withstand high temperatures for the industrial production of fucoxanthin, offering valuable insights for practical utilization.IMPORTANCEThis study delves into Phaeodactylum tricornutum's response to high-temperature stress, specifically focusing on cis-zeatin. We uncover inhibited growth, reduced fucoxanthin, and significant cis-zeatin-related gene expression under high temperatures, highlighting potential signaling mechanisms. Crucially, genetic engineering and exogenous addition experiments confirm that the change in cis-zeatin levels could influence P. tricornutum's resistance to high-temperature stress. This breakthrough deepens our understanding of microalgae adaptation to high temperatures and offers an innovative angle for industrial fucoxanthin production. This research is a pivotal step toward developing heat-resistant microalgae for industrial use.

A chronic signaling TGFb zebrafish reporter identifies immune response in melanoma

Developmental signaling pathways associated with growth factors such as TGFb are commonly dysregulated in melanoma. Here we identified a human TGFb enhancer specifically activated in melanoma cells treated with TGFB1 ligand. We generated stable transgenic zebrafish with this TGFb Induced Enhancer driving green fluorescent protein (TIE:EGFP). TIE:EGFP was not expressed in normal melanocytes or early melanomas but was expressed in spatially distinct regions of advanced melanomas. Single-cell RNA-sequencing revealed that TIE:EGFP+ melanoma cells down-regulated interferon response while up-regulating a novel set of chronic TGFb target genes. ChIP-sequencing demonstrated that AP-1 factor binding is required for activation of chronic TGFb response. Overexpression of SATB2, a chromatin remodeler associated with tumor spreading, showed activation of TGFb signaling in early melanomas. Confocal imaging and flow cytometric analysis showed that macrophages localize to TIE:EGFP+ regions and preferentially phagocytose TIE:EGFP+ melanoma cells compared to TIE:EGFP- melanoma cells. This work identifies a TGFb induced immune response and demonstrates the need for the development of chronic TGFb biomarkers to predict patient response to TGFb inhibitors.

Investigation of selective glucocorticoid receptor modulation in high-grade serous ovarian cancer PDX models

OBJECTIVE

In ovarian cancer (OvCa), tumor cell high glucocorticoid receptor (GR) has been associated with poor patient prognosis. In vitro, GR activation inhibits chemotherapy-induced OvCa cell death in association with transcriptional upregulation of genes encoding anti-apoptotic proteins. A recent randomized phase II study demonstrated improvement in progression-free survival (PFS) for heavily pre-treated OvCa patients randomized to receive therapy with a selective GR modulator (SGRM) plus chemotherapy compared to chemotherapy alone. We hypothesized that SGRM therapy would improve carboplatin response in OvCa patient-derived xenograft (PDX).

METHODS

Six high-grade serous (HGS) OvCa PDX models expressing GR mRNA (NR3C1) and protein were treated with chemotherapy +/- SGRM. Tumor size was measured longitudinally by peritoneal transcutaneous ultrasonography.

RESULTS

One of the 6 GR-positive PDX models showed a significant improvement in PFS with the addition of a SGRM. Interestingly, the single model with an improved PFS was least carboplatin sensitive. Possible explanations for the modest SGRM activity include the high carboplatin sensitivity of 5 of the PDX tumors and the potential that SGRMs activate the tumor invasive immune cells in patients (absent from immunocompromised mice). The level of tumor GR protein expression alone appears insufficient for predicting SGRM response.

CONCLUSION

The significant improvement in PFS shown in 1 of the 6 models after treatment with a SGRM plus chemotherapy underscores the need to determine predictive biomarkers for SGRM therapy in HGS OvCa and to better identify patient subgroups that are most likely to benefit from adding GR modulation to chemotherapy.