The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights.

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Initiation of Parental Genome Reprogramming in Fertilized Oocyte by Splicing Kinase SRPK1-Catalyzed Protamine Phosphorylation.

The paternal genome undergoes a massive exchange of histone with protamine for compaction into sperm during spermiogenesis. Upon fertilization, this process is potently reversed, which is essential for parental genome reprogramming and subsequent activation; however, it remains poorly understood how this fundamental process is initiated and regulated. Here, we report that the previously characterized splicing kinase SRPK1 initiates this life-beginning event by catalyzing site-specific phosphorylation of protamine, thereby triggering protamine-to-histone exchange in the fertilized oocyte. Interestingly, protamine undergoes a DNA-dependent phase transition to gel-like condensates and SRPK1-mediated phosphorylation likely helps open up such structures to enhance protamine dismissal by nucleoplasmin (NPM2) and enable the recruitment of HIRA for H3.3 deposition. Remarkably, genome-wide assay for transposase-accessible chromatin sequencing (ATAC-seq) analysis reveals that selective chromatin accessibility in both sperm and MII oocytes is largely erased in early pronuclei in a protamine phosphorylation-dependent manner, suggesting that SRPK1-catalyzed phosphorylation initiates a highly synchronized reorganization program in both parental genomes.

LncRNA CCTT-mediated RNA-DNA and RNA-protein interactions facilitate the recruitment of CENP-C to centromeric DNA during kinetochore assembly.

Kinetochore assembly on centromeres is central for chromosome segregation, and defects in this process cause mitotic errors and aneuploidy. Besides the well-established protein network, emerging evidence suggests the involvement of regulatory RNA in kinetochore assembly; however, it has remained elusive about the identity of such RNA, let alone its mechanism of action in this critical process. Here, we report CCTT, a previously uncharacterized long non-coding RNA (lncRNA) transcribed from the arm of human chromosome 17, which plays a vital role in kinetochore assembly. We show that CCTT highly localizes to all centromeres via the formation of RNA-DNA triplex and specifically interacts with CENP-C to help engage this blueprint protein in centromeres, and consequently, CCTT loss triggers extensive mitotic errors and aneuploidy. These findings uncover a non-centromere-derived lncRNA that recruits CENP-C to centromeres and shed critical lights on the function of centromeric DNA sequences as anchor points for kinetochore assembly.

Evolution of the germline mutation rate across vertebrates.

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself1. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies2. Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis3. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

The Augmented R-Loop Is a Unifying Mechanism for Myelodysplastic Syndromes Induced by High-Risk Splicing Factor Mutations.

Mutations in several general pre-mRNA splicing factors have been linked to myelodysplastic syndromes (MDSs) and solid tumors. These mutations have generally been assumed to cause disease by the resultant splicing defects, but different mutations appear to induce distinct splicing defects, raising the possibility that an alternative common mechanism is involved. Here we report a chain of events triggered by multiple splicing factor mutations, especially high-risk alleles in SRSF2 and U2AF1, including elevated R-loops, replication stress, and activation of the ataxia telangiectasia and Rad3-related protein (ATR)-Chk1 pathway. We further demonstrate that enhanced R-loops, opposite to the expectation from gained RNA binding with mutant SRSF2, result from impaired transcription pause release because the mutant protein loses its ability to extract the RNA polymerase II (Pol II) C-terminal domain (CTD) kinase-the positive transcription elongation factor complex (P-TEFb)-from the 7SK complex. Enhanced R-loops are linked to compromised proliferation of bone-marrow-derived blood progenitors, which can be partially rescued by RNase H overexpression, suggesting a direct contribution of augmented R-loops to the MDS phenotype.

ceRNA crosstalk mediated by ncRNAs is a novel regulatory mechanism in fish sex determination and differentiation.

Competing endogenous RNAs (ceRNAs) are vital regulators of gene networks in mammals. The involvement of noncoding RNAs (ncRNAs) as ceRNA in genotypic sex determination (GSD) and environmental sex determination (ESD) in fish is unknown. The Chinese tongue sole, which has both GSD and ESD mechanisms, was used to map the dynamic expression pattern of ncRNAs and mRNA in gonads during sex determination and differentiation. Transcript expression patterns shift during the sex differentiation phase, and ceRNA modulation occurs through crosstalk of differentially expressed long ncRNAs (lncRNAs), circular RNAs (circRNAs), microRNAs (miRNAs), and sex-related genes in fish. Of note was the significant up-regulation of a circRNA from the sex-determining gene dmrt1 (circular RNA dmrt1) and a lncRNA, called AMSDT (which stands for associated with male sex differentiation of tongue sole) in Chinese tongue sole testis. These two ncRNAs both share the same miRNA response elements with gsdf, which has an up-regulated expression when they bind to miRNA cse-miR-196 and concurrent down-regulated female sex-related genes to facilitate testis differentiation. This is the first demonstration in fish that ceRNA crosstalk mediated by ncRNAs modulates sexual development and unveils a novel regulatory mechanism for sex determination and differentiation.

R-ChIP Using Inactive RNase H Reveals Dynamic Coupling of R-loops with Transcriptional Pausing at Gene Promoters.

R-loop, a three-stranded RNA/DNA structure, has been linked to induced genome instability and regulated gene expression. To enable precision analysis of R-loops in vivo, we develop an RNase-H-based approach; this reveals predominant R-loop formation near gene promoters with strong G/C skew and propensity to form G-quadruplex in non-template DNA, corroborating with all biochemically established properties of R-loops. Transcription perturbation experiments further indicate that R-loop induction correlates to transcriptional pausing. Interestingly, we note that most mapped R-loops are each linked to a nearby free RNA end; by using a ribozyme to co-transcriptionally cleave nascent RNA, we demonstrate that such a free RNA end coupled with a G/C-skewed sequence is necessary and sufficient to induce R-loop. These findings provide a topological solution for RNA invasion into duplex DNA and suggest an order for R-loop initiation and elongation in an opposite direction to that previously proposed.

Single-cell transcriptome landscape of zebrafish liver reveals hepatocytes and immune cell interactions in understanding nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is becoming the most common chronic liver disease in the world. Immunity is the major contributing factor in NAFLD; however, the interaction of immune cells and hepatocytes in disease progression has not been fully elucidated. As a popular species for studying NAFLD, zebrafish, whose liver is a complex immune system mediated by immune cells and non-immune cells in maintaining immune tolerance and homeostasis. Understanding the cellular composition and immune environment of zebrafish liver is of great significance for its application in NAFLD. Here, we established a liver atlas that consists of 10 cell types using single-cell RNA sequencing (scRNA-seq). By examining the heterogeneity of hepatocytes and analyzing the expression of NAFLD-associated genes in the specific cluster, we provide a potential target cell model to study NAFLD. Additionally, our analysis identified two subtypes of distinct resident macrophages with inflammatory and non-inflammatory functions and characterized the successive stepwise development of T cell subclusters in the liver. Importantly, we uncovered the possible regulation of macrophages and T cells on target cells of fatty liver by analyzing the cellular interaction between hepatocytes and immune cells. Our data provide valuable information for an in-depth study of immune cells targeting hepatocytes to regulate the immune balance in NAFLD.

Examining the Impact of Phased Nursing within the Chronic Disease Trajectory Model on Glioma Patient Outcomes.

Background: Glioma (GL) , a primary brain tumor, presents significant challenges in patient care due to its complex disease trajectory and psychological impact. Phased nursing interventions, grounded in the Chronic Illness Trajectory Model (CITM), offer a holistic approach to addressing these multifaceted needs. Objective: The objective of this study was to assess the impact of phased nursing within the CITM on the psychological well-being, quality of life, and cancer-related fatigue (CRF) of glioma patients. Methods: A total of 100 GL patients undergoing treatment at our hospital between February 2020 and February 2021 were enrolled in this randomized controlled trial. Patients were randomly assigned to either the control group, which received standard routine care, or the observation group, which received phased nursing interventions based on the CITM framework. The mental state, quality of life, and CRF scores of the patients were assessed using validated measures at baseline and following the intervention period. Statistical analyses were conducted to compare the outcomes between the two groups. Results: The findings revealed that patients in the observation group exhibited significantly higher scores in mental state and quality of life domains compared to those in the control group (P < .05). Additionally, patients receiving phased nursing showed a significant reduction in CRF scores post-intervention. These results indicate that phased nursing within the CITM framework has a beneficial effect on the psychological well-being and overall quality of life of GL patients while also mitigating CRF. Conclusions: Our findings suggest that incorporating phased nursing interventions into the care of GL patients can lead to improvements in psychological outcomes, CRF, and quality of life. These findings underscore the importance of adopting holistic approaches to patient care, particularly in chronic disease management.

The Role of Z Chromosome Localization Gene psmd9 in Spermatogenesis of Cynoglossus semilaevis.

Proteasome 26S Subunit, Non-ATPase 9 (psmd9) plays an important role in the balance of protamine and the stability of the nucleolar structure during spermatogenesis. In this study, we cloned the psmd9 of Cynoglossus semilaevis and analyzed its expression pattern. psmd9 was identified on the Z chromosome of C. semilaevis, which is considered an interesting candidate gene for spermatogenesis. qRT-PCR and FISH experiments showed that the psmd9 gene was significantly highly expressed in the testes. It is worth noting that the expression level of psmd9 in male fish testes is significantly higher than that in pseudomales. In order to further explore the role of psmd9 in spermatogenesis, a male testicular cell line was used as the experimental material. The results of the psmd9-RNAi and overexpression experiments showed that psmd9 had a synergistic effect with spermatogenesis-related genes dnd1, cfap69, dnah3 and dnajb13, but had an antagonistic effect with ccne2. Our findings offer a scientific foundation for comprehending the role of psmd9 in the spermatogenesis regulatory network of C. semilaevis.

Active retrotransposons help maintain pericentromeric heterochromatin required for faithful cell division.

Retrotransposons are populated in vertebrate genomes, and when active, are thought to cause genome instability with potential benefit to genome evolution. Retrotransposon-derived RNAs are also known to give rise to small endo-siRNAs to help maintain heterochromatin at their sites of transcription; however, as not all heterochromatic regions are equally active in transcription, it remains unclear how heterochromatin is maintained across the genome. Here, we address these problems by defining the origins of repeat-derived RNAs and their specific chromatin locations in Drosophila S2 cells. We demonstrate that repeat RNAs are predominantly derived from active gypsy elements and processed by Dcr-2 into small RNAs to help maintain pericentromeric heterochromatin. We also show in cultured S2 cells that synthetic repeat-derived endo-siRNA mimics are sufficient to rescue Dcr-2-deficiency-induced defects in heterochromatin formation in interphase and chromosome segregation during mitosis, demonstrating that active retrotransposons are required for stable genetic inheritance.

The Role of the Heat Shock Cognate Protein 70 Genes in Sex Determination and Differentiation of Chinese Tongue Sole (Cynoglossus semilaevis).

Fish sex determination can be affected by environmental temperature. This process relies on temperature-sensitive proteins such as heat shock proteins (HSPs). Our previous work found that heat shock cognate proteins (HSCs) may participate in high-temperature associated sex reversal of Chinese tongue sole (Cynoglossus semilaevis). However, the role of hsc genes in responding to high temperature and affecting sex determination/differentiation remains unclear. Here, by using C. semilaevis as model, we identified hsc70 and hsc70-like. hsc70 was abundant in the gonads with a testicular-higher expression at all gonadal development stages except for 6 months post fertilization (mpf). Intriguingly, hsc70-like showed higher expression in testes from 6 mpf on. Both long-term heat treatment during the temperature-sensitive sex-determining period and short-term heat stress at the end of this period caused different expression of hsc70/hsc70-like between sexes. The dual-luciferase assay results also suggested that these genes can respond to high temperature rapidly in vitro. Heat treatment of C. semilaevis testis cells overexpressed with hsc70/hsc70-like could affect the expression of sex-related genes sox9a and cyp19a1a. Our results indicated that hsc70 and hsc70-like were key regulators linking external high-temperature signals with sex differentiation in vivo and provide a new idea for understanding the mechanism by which high temperature affects sex determination/differentiation in teleosts.

RBFox2-miR-34a-Jph2 axis contributes to cardiac decompensation during heart failure.

Heart performance relies on highly coordinated excitation-contraction (EC) coupling, and defects in this critical process may be exacerbated by additional genetic defects and/or environmental insults to cause eventual heart failure. Here we report a regulatory pathway consisting of the RNA binding protein RBFox2, a stress-induced microRNA miR-34a, and the essential EC coupler JPH2. In this pathway, initial cardiac defects diminish RBFox2 expression, which induces transcriptional repression of miR-34a, and elevated miR-34a targets Jph2 to impair EC coupling, which further manifests heart dysfunction, leading to progressive heart failure. The key contribution of miR-34a to this process is further established by administrating its mimic, which is sufficient to induce cardiac defects, and by using its antagomir to alleviate RBFox2 depletion-induced heart dysfunction. These findings elucidate a potential feed-forward mechanism to account for a critical transition to cardiac decompensation and suggest a potential therapeutic avenue against heart failure.

Cloning of the Maternal Effector Gene org and Its Regulation by lncRNA ORG-AS in Chinese Tongue Sole (Cynoglossus semilaevis).

Maternal effector genes (MEGs) encode maternal RNA and protein, accumulating in the cytoplasm of oocytes. During oocyte development, MEGs participate in oocyte meiosis and promote oocyte development. And MEGs can also regulate maternal transcriptome stability and promote maternal-zygotic transition (MTZ) in early embryonic development. Long noncoding RNAs (lncRNAs), as new epigenetic regulators, can regulate gene expression at both the transcriptional and post-transcriptional levels through cis- or trans-regulation. The oogenesis-related gene org is a germ-cell-specific gene in fish, but the role of org in embryonic development and oogenesis has rarely been studied, and the knowledge of the lncRNA-mediated regulation of org is limited. In this study, we cloned and identified the org gene of Chinese tongue sole (Cynoglossus semilaevis), and we identified a lncRNA named lncRNA ORG-anti-sequence (ORG-AS), located at the reverse overlapping region of org. The results of qRT-PCR and FISH demonstrated that org was highly expressed during the early stages of embryonic development and oogenesis and was located in the cytoplasm of oocytes. ORG-AS was expressed at low levels in the ovary and colocalized with org in the cytoplasm of oocytes. In vitro experiments showed that overexpression of ORG-AS inhibited org expression. These results suggest that org, as a MEG in C. semilaevis, participates in the MTZ and the oogenesis. The lncRNA ORG-AS negatively regulates the gene expression of org through trans-regulation. These new findings broaden the function of MEGs in embryonic development and the oogenesis of bony fish and prove that lncRNAs are important molecular factors regulating org.

Capturing the interactome of newly transcribed RNA.

We combine the labeling of newly transcribed RNAs with 5-ethynyluridine with the characterization of bound proteins. This approach, named capture of the newly transcribed RNA interactome using click chemistry (RICK), systematically captures proteins bound to a wide range of RNAs, including nascent RNAs and traditionally neglected nonpolyadenylated RNAs. RICK has identified mitotic regulators amongst other novel RNA-binding proteins with preferential affinity for nonpolyadenylated RNAs, revealed a link between metabolic enzymes/factors and nascent RNAs, and expanded the known RNA-bound proteome of mouse embryonic stem cells. RICK will facilitate an in-depth interrogation of the total RNA-bound proteome in different cells and systems.

Footprints of global change in marine life: Inferring past environment based on DNA methylation and gene expression marks.

Ocean global warming affects the distribution, life history and physiology of marine life. Extreme events, like marine heatwaves, are increasing in frequency and intensity. During sensitive stages of early fish development, the consequences may be long-lasting and mediated by epigenetic mechanisms. Here, we used European sea bass as a model to study the possible long-lasting effects of a marine heatwave during early development. We measured DNA methylation and gene expression in four tissues (brain, muscle, liver and testis) and detected differentially methylated regions (DMRs). Six genes were differentially expressed and contained DMRs three years after exposure to increased temperature, indicating direct phenotypic consequences and representing persistent changes. Interestingly, nine genes contained DMRs around the same genomic regions across tissues, therefore consisting of common footprints of developmental temperature in environmentally responsive loci. These loci are, to our knowledge, the first metastable epialleles (MEs) described in fish. MEs may serve as biomarkers to infer past life history events linked with persistent consequences. These results highlight the importance of subtle phenotypic changes mediated by epigenetics to extreme weather events during sensitive life stages. Also, to our knowledge, it is the first time the molecular effects of a marine heatwave during the lifetime of individuals are assessed. MEs could be used in surveillance programs aimed at determining the footprints of climate change on marine life. Our study paves the way for the identification of conserved MEs that respond equally to environmental perturbations across species. Conserved MEs would constitute a tool of assessment of global change effects in marine life at a large scale.

Early response to heat stress in Chinese tongue sole (Cynoglossus semilaevis): performance of different sexes, candidate genes and networks.

BACKGROUND: Temperature is known to affect living organisms and alter the expression of responsive genes, which affects a series of life processes, such as development, reproduction and metabolism. Several genes and gene families have been involved in high temperature responses, such as heat shock protein (hsp) family, Jumonji family and genes related to cortisol synthesis. Gonad is a vital organ related to the existence of a species. However, the comprehensive understanding of gonadal responses to environmental temperature is limited. RESULTS: To explore the effects of environmental temperature on genes and gene networks in gonads, we performed acute heat treatment (48 h) on Chinese tongue sole (Cynoglossus semilaevis). Gonadal transcriptome analysis was conducted on females, pseudomales and males exposed to high (28 °C) and normal (22 °C) temperatures. A total of 1226.24 million clean reads were obtained from 18 libraries. Principal component analysis (PCA) and differentially expressed gene (DEG) analysis revealed different performance of sex responses to heat stress. There were 4565, 790 and 1117 specific genes altered their expression level in females, pseudomales and males, respectively. Of these, genes related to hsp gene family, cortisol synthesis and metabolism and epigenetic regulation were involved in early heat response. Furthermore, a total of 1048 DEGs were shared among females, pesudomales and males, which may represent the inherent difference between high and normal temperatures. Genes, such as eef1akmt3, eef1akmt4, pnmt and hsp family members, were found. CONCLUSIONS: Our results depicted for the first time the gonadal gene expression under acute high temperature treatment in Chinese tongue sole. The findings may provide a clue for understanding the responses of genes and networks to environmental temperature.

Nuclear matrix factor hnRNP U/SAF-A exerts a global control of alternative splicing by regulating U2 snRNP maturation.

The nuclear matrix-associated hnRNP U/SAF-A protein has been implicated in diverse pathways from transcriptional regulation to telomere length control to X inactivation, but the precise mechanism underlying each of these processes has remained elusive. Here, we report hnRNP U as a regulator of SMN2 splicing from a custom RNAi screen. Genome-wide analysis by CLIP-seq reveals that hnRNP U binds virtually to all classes of regulatory noncoding RNAs, including all snRNAs required for splicing of both major and minor classes of introns, leading to the discovery that hnRNP U regulates U2 snRNP maturation and Cajal body morphology in the nucleus. Global analysis of hnRNP U-dependent splicing by RNA-seq coupled with bioinformatic analysis of associated splicing signals suggests a general rule for splice site selection through modulating the core splicing machinery. These findings exemplify hnRNP U/SAF-A as a potent regulator of nuclear ribonucleoprotein particles in diverse gene expression pathways.

SPOP-containing complex regulates SETD2 stability and H3K36me3-coupled alternative splicing.

Trimethylation of histone H3K36 is a chromatin mark associated with active gene expression, which has been implicated in coupling transcription with mRNA splicing and DNA damage response. SETD2 is a major H3K36 trimethyltransferase, which has been implicated as a tumor suppressor in mammals. Here, we report the regulation of SETD2 protein stability by the proteasome system, and the identification of SPOP, a key subunit of the CUL3 ubiquitin E3 ligase complex, as a SETD2-interacting protein. We demonstrate that SPOP is critically involved in SETD2 stability control and that the SPOP/CUL3 complex is responsible for SETD2 polyubiquitination both in vivo and in vitro ChIP-Seq analysis and biochemical experiments demonstrate that modulation of SPOP expression confers differential H3K36me3 on SETD2 target genes, and induce H3K36me3-coupled alternative splicing events. Together, these findings establish a functional connection between oncogenic SPOP and tumor suppressive SETD2 in the dynamic regulation of gene expression on chromatin.

Distinct splicing signatures affect converged pathways in myelodysplastic syndrome patients carrying mutations in different splicing regulators.

Myelodysplastic syndromes (MDS) are heterogeneous myeloid disorders with prevalent mutations in several splicing factors, but the splicing programs linked to specific mutations or MDS in general remain to be systematically defined. We applied RASL-seq, a sensitive and cost-effective platform, to interrogate 5502 annotated splicing events in 169 samples from MDS patients or healthy individuals. We found that splicing signatures associated with normal hematopoietic lineages are largely related to cell signaling and differentiation programs, whereas MDS-linked signatures are primarily involved in cell cycle control and DNA damage responses. Despite the shared roles of affected splicing factors in the 3' splice site definition, mutations in U2AF1, SRSF2, and SF3B1 affect divergent splicing programs, and interestingly, the affected genes fall into converging cancer-related pathways. A risk score derived from 11 splicing events appears to be independently associated with an MDS prognosis and AML transformation, suggesting potential clinical relevance of altered splicing patterns in MDS.