Gene expression across mammalian organ development

Comparative analysis of testicular fusion in Spodoptera litura (cutworm) and Bombyx mori (silkworm): Histological and transcriptomic insights

Spodoptera litura commonly known as the cutworm, is among the most destructive lepidopteran pests affecting over 120 plants species. The powerful destructive nature of this lepidopteran is attributable to its high reproductive capacity. The testicular fusion that occurs during metamorphosis from larvae to pupa in S.litura positively influences the reproductive success of the offspring. In contrast, Bombyx mori, the silkworm, retains separate testes throughout its life and does not undergo this fusion process. Microscopic examination reveals that during testicular fusion in S.litura, the peritoneal sheath becomes thinner and more translucent, whereas in B.mori, the analogous region thickens. The outer basement membrane in S.litura exhibits fractures, discontinuity, and uneven thickness accompanied by a significant presence of cellular secretions, large cell size, increased vesicles, liquid droplets, and a proliferation of rough endoplasmic reticulum and mitochondria. In contrast, the testicular peritoneal sheath of B.mori at comparable developmental stage exhibits minimal change. Comparative transcriptomic analysis of the testicular peritoneal sheath reveals a substantial difference in gene expression between the two species. The disparity in differential expressed genes (DEGs) is linked to an enrichment of numerous transcription factors, intracellular signaling pathways involving Ca2+ and GTPase, as well as intracellular protein transport and signaling pathways. Meanwhile, structural proteins including actin, chitin-binding proteins, membrane protein fractions, cell adhesion, extracellular matrix proteins are predominantly identified. Moreover, the study highlights the enrichment of endopeptidases, serine proteases, proteolytic enzymes and matrix metalloproteins, which may play a role in the degradation of the outer membrane. Five transcription factors-Slforkhead, Slproline, Slcyclic, Slsilk, and SlD-ETS were identified, and their expression pattern were confirmed by qRT-PCR. they are candidates for participating in the regulation of testicular fusion. Our findings underscore significant morphological and trancriptomic variation in the testicular peritoneal sheath of S.litura compared to the silkworm, with substantial changes at the transcriptomic level coinciding with testicular fusion. The research provides valuable clues for understanding the complex mechanisms underlying this unique phenomenon in insects.

The hepatokine FGF21 stopped lipogenesis and reduced testosterone production in mLTC-1 Leydig Cell Line

Beyond their link to metabolic issues like type 2 diabetes, factors like lifestyle, environment, and excess weight may also influence fertility. Fibroblast growth factor 21 (FGF21), a liver-derived hormone linked to energy balance, has recently emerged as a potential player in female mammalian reproduction. In male, only two studies have described potential effects of FGF21 on fertility. A recent study has described a negative correlation observed in obese patients presenting a low testosterone level associated with elevated FGF21 plasma levels. To investigate the role of FGF21 in steroidogenesis, we have studied the involvement of FGF21 in lipid and steroid activity by Leydig cells. Leydig cell model expressed all FGF21 receptors and β-Klotho cofactor as determined by RT-qPCR and by western-blot. Cultured mLTC-1 Leydig cell line exposed to increasing FGF21 concentration induced phosphorylation (Ser 473) of Akt and modified the CREB factor activity, suggesting the functionality of the FGF21 pathway. FGF21 consequences on mLTC-1 Leydig cells are inhibition of the lipid synthesis, leading to a reduction in the content of lipid droplets. The drop in lipid synthesis is associated with a reduction in the amount of lipids (mainly PUFA, cholesterol esterified, and triglycerides) as measured by lipidomic approach. The main consequence is to reduce the quantity of cholesterol, the steroid precursor, in mLTC-1 Leydig cells and is associated with a low production in testosterone. The decrease in androgens was also associated with a reduction in the steroid enzyme genes expression, which are under the control of CREB activity, and present a lower activity due to low cAMP intracellular levels. In vivo, steroid production was lowering after FGF21 administration in adult male mice associated to a decrease in progressive motility and velocity of sperm. In addition, these experimental data are reinforced by a data mining analysis focused on "gonad" terms in 1,319,905 article references showing the link already described between FGF21 with the fatty acids pathways, cholesterol storage, and steroid production. In conclusion, we demonstrated that Leydig cells in the testes present a functional FGF21 pathway, which regulates lipid metabolism and steroid function. In mLTC-1 Leydig cells, FGF21 reduced cholesterol, PUFA content, and testosterone production. Finally, this work highlighted that the hepatokine FGF21 could have a negative impact on androgen synthesis and testicular activity.

DOT1-like histone lysine methyltransferase is critical for adult vessel maintenance and functions

OBJECTIVE

Disruptor of telomeric silencing 1-like (DOT1L) is the only known histone H3K79 methyltransferase essential for the development of the embryonic cardiovascular system, including the heart, blood vessels, and lymphatic vessels, through transcriptional regulation. Our previous study demonstrated that Dot1l deletion results in aberrant lymphatic development and function. However, its precise function in the postnatal cardiovascular system remains unknown.

METHODS

Using conditional and inducible Dot1l knockout (KO) mice, along with a reporter strain carrying the Geo gene at the Dot1l locus, DOT1L expression and its function in the vascular system during postnatal life were investigated. To assess vessel morphology and vascular permeability, we administered Latex or Evans blue dye to KO mice. In addition, in vitro tube formation and cell migration assays were performed using DOT1L-depleted human umbilical vein endothelial cells (HUVECs). Changes in the expression of vascular genes in HUVECs were measured by quantitative polymerase chain reaction.

RESULTS

Our findings demonstrate that conditional Dot1l knockout in the Tg (Tie2-cre) strain results in abnormal blood vessel formation and lymphatic anomalies in the intestine. In a mouse model of Rosa26-creER-mediated inducible Dot1l knockout, we observed vascular phenotypes, including increased vascular permeability and brain hemorrhage, when DOT1L was deleted in adulthood. Additionally, DOT1L depletion in cultured HUVECs led to impaired cell migration and tube formation, likely due to altered gene transcription. These findings highlight the essential role of DOT1L in maintaining vascular integrity and function during embryonic development and postnatal life.

CONCLUSION

Our study revealed that DOT1L is required for the maintenance of adult vascular function through the regulation of gene expression.

GEGA (Gallus Enriched Gene Annotation): an online tool providing genomics and functional information across 47 tissues for a chicken gene-enriched atlas gathering Ensembl and Refseq genome annotations

GEGA is a user-friendly tool designed to navigate through various genomic and functional information related to an enriched gene atlas in chicken that integrates the gene catalogues from the two reference databases, NCBI-RefSeq and EMBL-Ensembl/GENCODE, along with four additional rich resources such as FAANG and NONCODE. Using the latest GRCg7b genome assembly, GEGA encompasses a total of 78 323 genes, including 24 102 protein-coding genes (PCGs) and 44 428 long non-coding RNAs (lncRNAs), significantly increasing the number of genes provided by each resource independently. However, GEGA is more than just a gene database. It offers a range of features that allow us to go deeper into the functional aspects of these genes. Users can explore gene expression and co-expression profiles across 47 tissues from 36 datasets and 1400 samples, discover tissue-specific variations and their expression as a function of sex or age and extract orthologous genes or their genomic configuration relative to the closest gene. For the communities interested in a specific gene, a list of genes or a quantitative trait locus region in chicken, GEGA's user-friendly interface facilitates efficient gene analysis, easy downloading of results and a multitude of graphical representations, from genomic information to detailed visualization of expression levels.

Loss of Elp1 in cerebellar granule cell progenitors models ataxia phenotype of Familial Dysautonomia

Familial Dysautonomia (FD) is an autosomal recessive disorder caused by a splice site mutation in the gene ELP1, which disproportionally affects neurons. While classically characterized by deficits in sensory and autonomic neurons, neuronal defects in the central nervous system have also been described. Although ELP1 expression remains high in the normal developing and adult cerebellum, its role in cerebellar development is unknown. To explore the role of Elp1 in the cerebellum, we knocked out Elp1 in cerebellar granule cell progenitors (GCPs) and examined the outcome on animal behavior and cellular composition. We found that GCP-specific conditional knockout of Elp1 (Elp1cKO) resulted in ataxia by 8 weeks of age. Cellular characterization showed that the animals had smaller cerebella with fewer granule cells. This defect was already apparent as early as 7 days after birth, when Elp1cKO animals also had fewer mitotic GCPs and shorter Purkinje dendrites. Through molecular characterization, we found that loss of Elp1 was associated with an increase in apoptotic cell death and cell stress pathways in GCPs. Our study demonstrates the importance of ELP1 in the developing cerebellum, and suggests that loss of Elp1 in the GC lineage may also play a role in the progressive ataxia phenotypes of FD patients.

Lagging Brain Gene Expression Patterns of Drosophila melanogaster Young Adult Males Confound Comparisons Between Sexes

Many species, including fruit flies (Drosophila melanogaster), are sexually dimorphic. Phenotypic variation in morphology, physiology, and behavior can affect development, reproduction, health, and aging. Therefore, designating sex as a variable and sex-blocking should be considered when designing experiments. The brain regulates phenotypes throughout the lifespan by balancing survival and reproduction, and sex-specific development at each life stage is likely. Changes in morphology and physiology are governed by differential gene expression, a quantifiable molecular marker for age- and sex-specific variations. We assessed the fruit fly brain transcriptome at three adult ages for gene expression signatures of sex, age, and sex-by-age: 6698 genes were differentially expressed between sexes, with the most divergence at 3 days. Between ages, 31.1% of 6084 differentially expressed genes (1890 genes) share similar expression patterns from 3 to 7 days in females, and from 7 to 14 days in males. Most of these genes (90.5%, 1712) were upregulated and enriched for chemical stimulus detection and/or cilium regulation. Our data highlight an important delay in male brain gene regulation compared to females. Because significant delays in expression could confound comparisons between sexes, studies of sexual dimorphism at phenotypically comparable life stages rather than chronological age should be more biologically relevant.

A single-cell atlas of spatial and temporal gene expression in the mouse cranial neural plate

The formation of the mammalian brain requires regionalization and morphogenesis of the cranial neural plate, which transforms from an epithelial sheet into a closed tube that provides the structural foundation for neural patterning and circuit formation. Sonic hedgehog (SHH) signaling is important for cranial neural plate patterning and closure, but the transcriptional changes that give rise to the spatially regulated cell fates and behaviors that build the cranial neural tube have not been systematically analyzed. Here we used single-cell RNA sequencing to generate an atlas of gene expression at six consecutive stages of cranial neural tube closure in the mouse embryo. Ordering transcriptional profiles relative to the major axes of gene expression predicted spatially regulated expression of 870 genes along the anterior-posterior and mediolateral axes of the cranial neural plate and reproduced known expression patterns with over 85% accuracy. Single-cell RNA sequencing of embryos with activated SHH signaling revealed distinct SHH-regulated transcriptional programs in the developing forebrain, midbrain, and hindbrain, suggesting a complex interplay between anterior-posterior and mediolateral patterning systems. These results define a spatiotemporally resolved map of gene expression during cranial neural tube closure and provide a resource for investigating the transcriptional events that drive early mammalian brain development.

Adaptive Evolution and Functional Differentiation of Testis-Expressed Genes in Theria

Gene expression patterns differ in different tissues, and the expression pattern of genes in the mammalian testis is known to be extremely variable in different species. To clarify how the testis transcriptomic pattern has evolved in particular species, we examined the evolution of the adult testis transcriptome in Theria using 10 species: two marsupials (opossum and Tasmanian devil), six eutherian (placental) mammals (human, chimpanzee, bonobo, gorilla, rhesus macaque, and mouse), and two outgroup species (platypus and chicken). We show that 22 testis-expressed genes are marsupial-specific, suggesting their acquisition in the stem lineage of marsupials after the divergence from eutherians. Despite the time length of the eutherian stem lineage being similar to that of the marsupial lineage, acquisition of testis-expressed genes was not found in the stem lineage of eutherians; rather, their expression patterns differed by species, suggesting rapid gene evolution in the eutherian ancestors. Fifteen testis-expressed genes are therian-specific, and for three of these genes, the evolutionary tempo is markedly faster in eutherians than in marsupials. Our phylogenetic analysis of Rho GTPase-activating protein 28 (ARHGAP28) suggests the adaptive evolution of this gene in the eutherians, probably together with the expression pattern differentiation.

Loss-of-function of the Zinc Finger Homeobox 4 (ZFHX4) gene underlies a neurodevelopmental disorder

8q21.11 microdeletions encompassing the gene encoding transcription factor ZFHX4, have previously been associated by us with a syndromic form of intellectual disability, hypotonia, decreased balance and hearing loss. Here, we report on 57 individuals, 52 probands and 5 affected family members, with protein truncating variants (n=36), (micro)deletions (n=20) or an inversion (n=1) affecting ZFHX4 with variable developmental delay and intellectual disability, distinctive facial characteristics, morphological abnormalities of the central nervous system, behavioral alterations, short stature, hypotonia, and occasionally cleft palate and anterior segment dysgenesis. The phenotypes associated with 8q21.11 microdeletions and ZFHX4 intragenic loss-of-function variants largely overlap, identifying ZFHX4 as the main driver for the microdeletion syndrome, although leukocyte-derived DNA shows a mild common methylation profile for (micro)deletions only. We identify ZFHX4 as a transcription factor that is increasingly expressed during human brain development and neuronal differentiation. Furthermore, ZFHX4 interacting factors identified via IP-MS in neural progenitor cells, suggest an important role for ZFHX4 in cellular and developmental pathways, especially during histone modifications, cytosolic transport and development. Additionally, using CUT&RUN, we observed that ZFHX4 binds with the promoter regions of genes with crucial roles in embryonic, neuron and axon development. Since loss-of-function variants in ZFHX4 are found with consistent dysmorphic facial features, we investigated whether the disruption of zfhx4 causes craniofacial abnormalities in zebrafish. First-generation (F0) zfhx4 crispant zebrafish, (mosaic) mutant for zfhx4 loss-of-function variants, have significantly shorter Meckel's cartilages and smaller ethmoid plates compared to control zebrafish. Furthermore, behavioral assays show a decreased movement frequency in the zfhx4 crispant zebrafish in comparison with control zebrafish larvae. Although further research is needed, our in vivo work suggests a role for zfhx4 in facial skeleton patterning, palatal development and behavior.

Disentangling cell-intrinsic and extrinsic factors underlying evolution

A key goal of developmental biology is to determine the extent to which cells and organs develop autonomously, as opposed to requiring interactions with other cells or environmental factors. Chimeras have played a foundational role in this by enabling qualitative classification of cell-intrinsically vs. extrinsically driven processes. Here, we extend this framework to precisely decompose evolutionary divergence in any quantitative trait into cell-intrinsic, extrinsic, and intrinsic-extrinsic interaction components. Applying this framework to thousands of gene expression levels in reciprocal rat-mouse chimeras, we found that the majority of their divergence is attributable to cell-intrinsic factors, though extrinsic factors also play an integral role. For example, a rat-like extracellular environment extrinsically up-regulates the expression of a key transcriptional regulator of the endoplasmic reticulum (ER) stress response in some but not all cell types, which in turn strongly predicts extrinsic up-regulation of its target genes and of the ER stress response pathway as a whole. This effect is also seen at the protein level, suggesting propagation through multiple regulatory levels. Applying our framework to a cellular trait, neuronal differentiation, revealed a complex interaction of intrinsic and extrinsic factors. Finally, we show that imprinted genes are dramatically mis-expressed in species-mismatched environments, suggesting that mismatch between rapidly evolving intrinsic and extrinsic mechanisms controlling gene imprinting may contribute to barriers to interspecies chimerism. Overall, our conceptual framework opens new avenues to investigate the mechanistic basis of developmental processes and evolutionary divergence across myriad quantitative traits in any multicellular organism.

Sexual size dimorphism in mammals is associated with changes in the size of gene families related to brain development

In mammals, sexual size dimorphism often reflects the intensity of sexual selection, yet its connection to genomic evolution remains unexplored. Gene family size evolution can reflect shifts in the relative importance of different molecular functions. Here, we investigate the associate between brain development gene repertoire to sexual size dimorphism using 124 mammalian species. We reveal significant changes in gene family size associations with sexual size dimorphism. High levels of dimorphism correlate with an expansion of gene families enriched in olfactory sensory perception and a contraction of gene families associated with brain development functions, many of which exhibited particularly high expression in the human adult brain. These findings suggest a relationship between intense sexual selection and alterations in gene family size. These insights illustrate the complex interplay between sexual dimorphism, gene family size evolution, and their roles in mammalian brain development and function, offering a valuable understanding of mammalian genome evolution.

hnRNPs: roles in neurodevelopment and implication for brain disorders

Heterogeneous nuclear ribonucleoproteins (hnRNPs) constitute a family of multifunctional RNA-binding proteins able to process nuclear pre-mRNAs into mature mRNAs and regulate gene expression in multiple ways. They comprise at least 20 different members in mammals, named from A (HNRNP A1) to U (HNRNP U). Many of these proteins are components of the spliceosome complex and can modulate alternative splicing in a tissue-specific manner. Notably, while genes encoding hnRNPs exhibit ubiquitous expression, increasing evidence associate these proteins to various neurodevelopmental and neurodegenerative disorders, such as intellectual disability, epilepsy, microcephaly, amyotrophic lateral sclerosis, or dementias, highlighting their crucial role in the central nervous system. This review explores the evolution of the hnRNPs family, highlighting the emergence of numerous new members within this family, and sheds light on their implications for brain development.

PRDM6 promotes medulloblastoma by repressing chromatin accessibility and altering gene expression

SNCAIP duplication may promote Group 4 medulloblastoma via induction of PRDM6, a poorly characterized member of the PRDF1 and RIZ1 homology domain-containing (PRDM) family of transcription factors. Here, we investigated the function of PRDM6 in human hindbrain neuroepithelial stem cells and tested PRDM6 as a driver of Group 4 medulloblastoma. We report that human PRDM6 localizes predominantly to the nucleus, where it causes widespread repression of chromatin accessibility and complex alterations of gene expression patterns. Genome-wide mapping of PRDM6 binding reveals that PRDM6 binds to chromatin regions marked by histone H3 lysine 27 trimethylation that are located within, or proximal to, genes. Moreover, we show that PRDM6 expression in neuroepithelial stem cells promotes medulloblastoma. Surprisingly, medulloblastomas derived from PRDM6-expressing neuroepithelial stem cells match human Group 3, but not Group 4, medulloblastoma. We conclude that PRDM6 expression has oncogenic potential but is insufficient to drive Group 4 medulloblastoma from neuroepithelial stem cells. We propose that both PRDM6 and additional factors, such as specific cell-of-origin features, are required for Group 4 medulloblastoma. Given the lack of PRDM6 expression in normal tissues and its oncogenic potential shown here, we suggest that PRDM6 inhibition may have therapeutic value in PRDM6-expressing medulloblastomas.

The mouse multi-organ proteome from infancy to adulthood

The early-life organ development and maturation shape the fundamental blueprint for later-life phenotype. However, a multi-organ proteome atlas from infancy to adulthood is currently not available. Herein, we present a comprehensive proteomic analysis of ten mouse organs (brain, heart, lung, liver, kidney, spleen, stomach, intestine, muscle and skin) at three crucial developmental stages (1-, 4- and 8-weeks after birth) acquired using data-independent acquisition mass spectrometry. We detect and quantify 11,533 protein groups across the ten organs and obtain 115 age-related differentially expressed protein groups that are co-expressed in all organs from infancy to adulthood. We find that spliceosome proteins prevalently play crucial regulatory roles in the early-life development of multiple organs, and detect organ-specific expression patterns and sexual dimorphism. This multi-organ proteome atlas provides a fundamental resource for understanding the molecular mechanisms underlying early-life organ development and maturation.

The Evolution of Ultraconserved Elements in Vertebrates

Ultraconserved elements were discovered two decades ago, arbitrarily defined as sequences that are identical over a length ≥ 200 bp in the human, mouse, and rat genomes. The definition was subsequently extended to sequences ≥ 100 bp identical in at least three of five mammalian genomes (including dog and cow), and shown to have undergone rapid expansion from ancestors in fish and strong negative selection in birds and mammals. Since then, many more genomes have become available, allowing better definition and more thorough examination of ultraconserved element distribution and evolutionary history. We developed a fast and flexible analytical pipeline for identifying ultraconserved elements in multiple genomes, dedUCE, which allows manipulation of minimum length, sequence identity, and number of species with a detectable ultraconserved element according to specified parameters. We suggest an updated definition of ultraconserved elements as sequences ≥ 100 bp and ≥97% sequence identity in ≥50% of placental mammal orders (12,813 ultraconserved elements). By mapping ultraconserved elements to ∼200 species, we find that placental ultraconserved elements appeared early in vertebrate evolution, well before land colonization, suggesting that the evolutionary pressures driving ultraconserved element selection were present in aquatic environments in the Cambrian-Devonian periods. Most (>90%) ultraconserved elements likely appeared after the divergence of gnathostomes from jawless predecessors, were largely established in sequence identity by early Sarcopterygii evolution-before the divergence of lobe-finned fishes from tetrapods-and became near fixed in the amniotes. Ultraconserved elements are mainly located in the introns of protein-coding and noncoding genes involved in neurological and skeletomuscular development, enriched in regulatory elements, and dynamically expressed throughout embryonic development.

ABC Efflux Transporters and Solute Carriers in the Early Developing Brain of a Marsupial Monodelphis domestica (South American Gray Short-Tailed Opossum)

This study used a marsupial Monodelphis domestica, which is born very immature and most of its development is postnatal without placental protection. RNA-sequencing (RNA-Seq) was used to identify the expression of influx and efflux transporters (ATP-binding cassettes [ABCs] and solute carriers [SLCs]) and metabolizing enzymes in brains of newborn to juvenile Monodelphis. Results were compared to published data in the developing eutherian rat. To test the functionality of these transporters at similar ages, the entry of paracetamol (acetaminophen) into the brain and cerebrospinal fluid (CSF) was measured using liquid scintillation counting following a single administration of the drug along with its radiolabelled tracer [3H]. Drug permeability studies found that in Monodelphis, brain entry of paracetamol was already restricted at P5; it decreased further in the first week of life and then remained stable until the oldest age group tested (P110). Transcriptomic analysis of Monodelphis brain showed that expression of transporters and their metabolizing enzymes in early postnatal (P) pups (P0, P5, and P8) was relatively similar, but by P109, many more transcripts were identified. When transcriptomes of newborn Monodelphis brain and E19 rat brain and placenta were compared, several transporters present in the rat placenta were also found in the newborn Monodelphis brain. These were absent from E19 rat brain but were present in the adult rat brain. These data indicate that despite its extreme immaturity, the newborn Monodelphis brain may compensate for the lack of placental protection during early brain development by upregulating protective mechanisms, which in eutherian animals are instead present in the placenta.

In Utero Extrahepatic Bile Duct Damage and Repair: Implications for Biliary Atresia

Biliary atresia (BA) is a cholangiopathy affecting the extrahepatic bile duct (EHBD) of newborns. The etiology and pathophysiology of BA are not fully understood; however, multiple causes of damage and obstruction of the neonatal EHBD have been identified. Initial damage to the EHBD likely occurs before birth. We discuss how different developmental stages in utero and birth itself could influence the susceptibility of the fetal EHBD to damage and a damaging wound-healing response. We propose that a damage-repair response of the fetal and neonatal EHBD involving redox stress and a program of fetal wound healing could-regardless of the cause of the initial damage-lead to either obstruction and BA or repair of the duct and recovery. This overarching concept should guide future research targeted toward identification of factors that contribute to recovery as opposed to progression of injury and fibrosis. Viewing BA through the lens of an in utero damage-repair response could open up new avenues for research and suggests exciting new therapeutic targets.

Disentangling the splicing factor programs underlying complex molecular phenotypes

The regulation of exon inclusion through alternative splicing tunes the cell's behavior by increasing the functional diversity of the transcriptome and the proteome. Splicing factors work in concert to generate gene isoform pools that contribute to cell phenotypes yet their activity is controlled by multiple regulatory and signaling layers. This hinders identification of functional, phenotype-specific splicing factors using traditional single-omic measurements, such as their mutational state or expression. To address this challenge, we propose repurposing the virtual inference of protein activity by enriched regulon analysis (VIPER) to measure splicing factor activity solely from their downstream exon transcriptomic inclusion signatures. This approach is effective in assessing the effect of co-occurring splicing factor perturbations, as well as their post-translational regulation. As proof of concept, we dissect recurrent splicing factor programs underlying tumorigenesis including aberrantly activated factors acting as oncogenes and inactivated ones acting as tumor suppressors, which are undetectable by more conventional methodologies. Activation and inactivation of these cancer splicing programs effectively stratifies overall survival, as well as cancer hallmarks such as proliferation and immune evasion. Altogether, repurposing network-based inference of protein activity for splicing factor networks distills common, functionally relevant splicing programs in otherwise heterogeneous molecular contexts.

Transcriptomic profiling of gastrointestinal tracts in dairy cattle during lactation reveals molecular adaptations for milk synthesis

During lactation, dairy cattle's digestive tract requires significant adaptations to meet the increased nutrient demands for milk production. As we attempt to improve milk-related traits through selective pressure, it is crucial to understand the biological functions of the epithelia of the rumen, small intestine, and colonic tissues in response to changes in physiological state driven by changes in nutrient demands for milk synthesis. In this study, we obtained a total of 108 transcriptome profiles from three tissues (epithelia of the colon, duodenum, and rumen) of five Holstein cows, spanning eight time points from the early, mid, late lactation periods to the dry period. On average 97.06% of reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We analyzed 27,607 gene expression patterns at multiple periods, enabling direct comparisons within and among tissues during different lactation stages, including early and peak lactation. We identified 1645, 813, and 2187 stage-specific genes in the colon, duodenum, and rumen, respectively, which were enriched for common or specific biological functions among different tissues. Time series analysis categorized the expressed genes within each tissue into four clusters. Furthermore, when the three tissues were analyzed collectively, 36 clusters of similarly expressed genes were identified. By integrating other comprehensive approaches such as gene co-expression analyses, functional enrichment, and cell type deconvolution, we gained profound insights into cattle lactation, revealing tissue-specific characteristics of the gastrointestinal tract and shedding light on the intricate molecular adaptations involved in nutrient absorption, immune regulation, and cellular processes for milk synthesis during lactation.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Characterization of dUTPase expression in mouse postnatal development and adult neurogenesis

The enzyme dUTPase has an essential role in maintaining genomic integrity. In mouse, nuclear and mitochondrial isoforms of the enzyme have been described. Here we present the isoform-specific mRNA expression levels in different murine organs during development using RT-qPCR. In this study, we analyzed organs of 14.5-day embryos and of postnatal 2-, 4-, 10-week- and 13-month-old mice. We demonstrate organ-, sex- and developmental stage-specific differences in the mRNA expression levels of both isoforms. We found high mRNA expression level of the nuclear isoform in the embryo brain, and the expression level remained relatively high in the adult brain as well. This was surprising, since dUTPase is known to play an important role in proliferating cells, and mass production of neural cells is completed by adulthood. Thus, we investigated the pattern of the dUTPase protein expression specifically in the adult brain with immunostaining and found that dUTPase is present in the germinative zones, the subventricular and the subgranular zones, where neurogenesis occurs and in the rostral migratory stream where neuroblasts migrate to the olfactory bulb. These novel findings suggest that dUTPase may have a role in cell differentiation and indicate that accurate dTTP biosynthesis can be vital, especially in neurogenesis.

Genetic and developmental divergence in the neural crest programme between cichlid fish species

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic programme in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to inter-specific morphological differences. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared to the teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes - particularly those controlled by sox10s - might be involved in generating morphological diversification between species and lays the groundwork for further investigations into mechanisms underpinning vertebrate NC diversification.

IMMUNOTAR - Integrative prioritization of cell surface targets for cancer immunotherapy

Cancer remains a leading cause of mortality globally. Recent improvements in survival have been facilitated by the development of less toxic immunotherapies; however, identifying targets for immunotherapies remains a challenge in the field. To address this challenge, we developed IMMUNOTAR, a computational tool that systematically prioritizes and identifies candidate immunotherapeutic targets. IMMUNOTAR integrates user-provided RNA-sequencing or proteomics data with quantitative features extracted from publicly available databases based on predefined optimal immunotherapeutic target criteria and quantitatively prioritizes potential surface protein targets. We demonstrate the utility and flexibility of IMMUNOTAR using three distinct datasets, validating its effectiveness in identifying both known and new potential immunotherapeutic targets within the analyzed cancer phenotypes. Overall, IMMUNOTAR enables the compilation of data from multiple sources into a unified platform, allowing users to simultaneously evaluate surface proteins across diverse criteria. By streamlining target identification, IMMUNOTAR empowers researchers to efficiently allocate resources and accelerate immunotherapy development.

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

BACKGROUND

DNA methylation is an important epigenetic modification which has numerous roles in modulating genome function. Its levels are spatially correlated across the genome, typically high in repressed regions but low in transcription factor (TF) binding sites and active regulatory regions. However, the mechanisms establishing genome-wide and TF binding site methylation patterns are still unclear.

RESULTS

Here we use a comparative approach to investigate the association of DNA methylation to TF binding evolution in mammals. Specifically, we experimentally profile DNA methylation and combine this with published occupancy profiles of five distinct TFs (CTCF, CEBPA, HNF4A, ONECUT1, FOXA1) in the liver of five mammalian species (human, macaque, mouse, rat, dog). TF binding sites are lowly methylated, but they often also have intermediate methylation levels. Furthermore, biding sites are influenced by the methylation status of CpGs in their wider binding regions even when CpGs are absent from the core binding motif. Employing a classification and clustering approach, we extract distinct and species-conserved patterns of DNA methylation levels at TF binding regions. CEBPA, HNF4A, ONECUT1, and FOXA1 share the same methylation patterns, while CTCF's differ. These patterns characterize alternative functions and chromatin landscapes of TF-bound regions. Leveraging our phylogenetic framework, we find DNA methylation gain upon evolutionary loss of TF occupancy, indicating coordinated evolution. Furthermore, each methylation pattern has its own evolutionary trajectory reflecting its genomic contexts.

CONCLUSIONS

Our epigenomic analyses indicate a role for DNA methylation in TF binding changes across species including that specific DNA methylation profiles characterize TF binding and are associated with their regulatory activity, chromatin contexts, and evolutionary trajectories.

Neuroinflammation and Lysosomal Abnormalities Characterise the Essential Role for Oxidation Resistance 1 in the Developing and Adult Cerebellum

Loss-of-function mutations in the TLDc family of proteins cause a range of severe childhood-onset neurological disorders with common clinical features that include cerebellar neurodegeneration, ataxia and epilepsy. Of these proteins, oxidation resistance 1 (OXR1) has been implicated in multiple cellular pathways related to antioxidant function, transcriptional regulation and cellular survival; yet how this relates to the specific neuropathological features in disease remains unclear. Here, we investigate a range of loss-of-function mouse model systems and reveal that constitutive deletion of Oxr1 leads to a rapid and striking neuroinflammatory response prior to neurodegeneration that is associated with lysosomal pathology. We go on to show that neuroinflammation and cell death in Oxr1 knockouts can be completely rescued by the neuronal expression of Oxr1, suggesting that the phenotype is driven by the cell-intrinsic defects of neuronal cells lacking the gene. Next, we generate a ubiquitous, adult inducible knockout of Oxr1 that surprisingly displays rapid-onset ataxia and cerebellar neurodegeneration, establishing for the first time that the distinctive pathology associated with the loss of Oxr1 occurs irrespective of developmental stage. Finally, we describe two new homozygous human pathogenic variants in OXR1 that cause neurodevelopmental delay, including a novel stop-gain mutation. We also compare functionally two missense human pathogenic mutations in OXR1, including one newly described here, that cause different clinical phenotypes but demonstrate partially retained neuroprotective activity against oxidative stress. Together, these data highlight the essential role of Oxr1 in modulating neuroinflammatory and lysosomal pathways in the mammalian brain and support the hypothesis that OXR1 protein dosage may be critical for pathological outcomes in disease.

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

GTDrift: a resource for exploring the interplay between genetic drift, genomic and transcriptomic characteristics in eukaryotes

We present GTDrift, a comprehensive data resource that enables explorations of genomic and transcriptomic characteristics alongside proxies of the intensity of genetic drift in individual species. This resource encompasses data for 1506 eukaryotic species, including 1413 animals and 93 green plants, and is organized in three components. The first two components contain approximations of the effective population size, which serve as indicators of the extent of random genetic drift within each species. In the first component, we meticulously investigated public databases to assemble data on life history traits such as longevity, adult body length and body mass for a set of 979 species. The second component includes estimations of the ratio between the rate of non-synonymous substitutions and the rate of synonymous substitutions (dN/dS) in protein-coding sequences for 1324 species. This ratio provides an estimate of the efficiency of natural selection in purging deleterious substitutions. Additionally, we present polymorphism-derived N e estimates for 66 species. The third component encompasses various genomic and transcriptomic characteristics. With this component, we aim to facilitate comparative transcriptomics analyses across species, by providing easy-to-use processed data for more than 16 000 RNA-seq samples across 491 species. These data include intron-centered alternative splicing frequencies, gene expression levels and sequencing depth statistics for each species, obtained with a homogeneous analysis protocol. To enable cross-species comparisons, we provide orthology predictions for conserved single-copy genes based on BUSCO gene sets. To illustrate the possible uses of this database, we identify the most frequently used introns for each gene and we assess how the sequencing depth available for each species affects our power to identify major and minor splice variants.

Expression of transforming growth factor β signalling molecules and their correlations with genes in loci linked to polycystic ovary syndrome in human foetal and adult tissues

Context Altered signalling of androgens, anti-Müllerian hormone or transforming growth factor beta (TGFβ) during foetal development have been implicated in the predisposition to polycystic ovary syndrome (PCOS) in later life, aside from its genetic predisposition. In foetal ovarian fibroblasts, TGFβ1 has been shown to regulate androgen signalling and seven genes located in loci associated with PCOS. Since PCOS exhibits a myriad of symptoms, it likely involves many different organs. Aims To identify the relationships between TGFβ signalling molecules and PCOS candidate genes in different tissues associated with PCOS. Methods Using RNA sequencing data, we examined the expression patterns of TGFβ signalling molecules in the human ovary, testis, heart, liver, kidney, brain tissue, and cerebellum from 4 to 20weeks of gestation and postnatally. We also examined the correlations between gene expression of TGFβ signalling molecules and PCOS candidate genes. Key results TGFβ signalling molecules were dynamically expressed in most tissues prenatally and/or postnatally. FBN3 , a PCOS candidate gene involved in TGFβ signalling, was expressed during foetal development in all tissues. The PCOS candidate genes HMGA2, YAP1 , and RAD50 correlated significantly (P TGFBR1 in six out of the seven tissues examined. Conclusions This study suggests that possible crosstalk occurs between genes in loci associated with PCOS and TGFβ signalling molecules in multiple tissues, particularly during foetal development. Implications Thus, alteration in TGFβ signalling during foetal development could affect many tissues contributing to the multiple phenotypes of PCOS in later life.

Decoding the epigenetic mechanism of mammalian sex determination

Sex determination embodies a dynamic and intricate developmental process wielding significant influence over the destiny of bipotential gonads, steering them towards male or female gonads. Gonadal differentiation and the postnatal manifestation of the gonadal phenotype involve a sophisticated interplay of transcription factors such as SOX9 and FOXL2. Central to this interplay are chromatin modifiers regulating the mutual antagonism during this interplay. In this review, the key findings and knowledge gaps in DNA methylation, histone modification, and non-coding RNA-mediated control throughout mammalian gonadal development are covered. Furthermore, it explores the role of the developing brain in playing a pivotal role in the initiation of gonadogenesis and the subsequent involvement of gonadal hormone/hormone receptor in fine-tuning sexual differentiation. Based on promising facts, the role of the developing brain through the hypothalamic pituitary gonadal axis is explained and suggested as a novel hypothesis. The article also discusses the potential impact of ecological factors on the human epigenome in relation to sex determination and trans-generational epigenetics in uncovering novel genes and mechanisms involved in sex determination and gonadal differentiation. We have subtly emphasized the disruptions in epigenetic regulations contributing to sexual disorders, which further allows us to raise certain questions, decipher approaches for handling these questions and setting up the direction of future research.

Evolution of tissue-specific expression of ancestral genes across vertebrates and insects

Regulation of gene expression is arguably the main mechanism underlying the phenotypic diversity of tissues within and between species. Here we assembled an extensive transcriptomic dataset covering 8 tissues across 20 bilaterian species and performed analyses using a symmetric phylogeny that allowed the combined and parallel investigation of gene expression evolution between vertebrates and insects. We specifically focused on widely conserved ancestral genes, identifying strong cores of pan-bilaterian tissue-specific genes and even larger groups that diverged to define vertebrate and insect tissues. Systematic inferences of tissue-specificity gains and losses show that nearly half of all ancestral genes have been recruited into tissue-specific transcriptomes. This occurred during both ancient and, especially, recent bilaterian evolution, with several gains being associated with the emergence of unique phenotypes (for example, novel cell types). Such pervasive evolution of tissue specificity was linked to gene duplication coupled with expression specialization of one of the copies, revealing an unappreciated prolonged effect of whole-genome duplications on recent vertebrate evolution.

Variant-specific pathophysiological mechanisms of AFF3 differently influence transcriptome profiles

BACKGROUND

We previously described the KINSSHIP syndrome, an autosomal dominant disorder associated with intellectual disability (ID), mesomelic dysplasia and horseshoe kidney, caused by de novo variants in the degron of AFF3. Mouse knock-ins and overexpression in zebrafish provided evidence for a dominant-negative mode of action, wherein an increased level of AFF3 resulted in pathological effects.

METHODS

Evolutionary constraints suggest that other modes-of-inheritance could be at play. We challenged this hypothesis by screening ID cohorts for individuals with predicted-to-be damaging variants in AFF3. We used both animal and cellular models to assess the deleteriousness of the identified variants.

RESULTS

We identified an individual with a KINSSHIP-like phenotype carrying a de novo partial duplication of AFF3 further strengthening the hypothesis that an increased level of AFF3 is pathological. We also detected seventeen individuals displaying a milder syndrome with either heterozygous Loss-of-Function (LoF) or biallelic missense variants in AFF3. Consistent with semi-dominance, we discovered three patients with homozygous LoF and one compound heterozygote for a LoF and a missense variant, who presented more severe phenotypes than their heterozygous parents. Matching zebrafish knockdowns exhibit neurological defects that could be rescued by expressing human AFF3 mRNA, confirming their association with the ablation of aff3. Conversely, some of the human AFF3 mRNAs carrying missense variants identified in affected individuals did not rescue these phenotypes. Overexpression of mutated AFF3 mRNAs in zebrafish embryos produced a significant increase of abnormal larvae compared to wild-type overexpression further demonstrating deleteriousness. To further assess the effect of AFF3 variation, we profiled the transcriptome of fibroblasts from affected individuals and engineered isogenic cells harboring + / + , KINSSHIP/KINSSHIP, LoF/ + , LoF/LoF or KINSSHIP/LoF AFF3 genotypes. The expression of more than a third of the AFF3 bound loci is modified in either the KINSSHIP/KINSSHIP or the LoF/LoF lines. While the same pathways are affected, only about one third of the differentially expressed genes are common to the homozygote datasets, indicating that AFF3 LoF and KINSSHIP variants largely modulate transcriptomes differently, e.g. the DNA repair pathway displayed opposite modulation.

CONCLUSIONS

Our results and the high pleiotropy shown by variation at this locus suggest that minute changes in AFF3 function are deleterious.

Applications of single-cell RNA sequencing in spermatogenesis and molecular evolution

Spermatogenic cell heterogeneity is determined by the complex process of spermatogenesis differentiation. However, effectively revealing the regulatory mechanisms underlying mammalian spermatogenic cell development and differentiation via traditional methods is difficult. Advances in technology have led to the emergence of many single-cell transcriptome sequencing protocols, which have partially addressed these challenges. In this review, we detail the principles of 10x Genomics technology and summarize the methods for downstream analysis of single-cell transcriptome sequencing data. Furthermore, we explore the role of single-cell transcriptome sequencing in revealing the heterogeneity of testicular ecological niche cells, delineating the establishment and disruption of testicular immune homeostasis during human spermatogenesis, investigating abnormal spermatogenesis in humans, and, ultimately, elucidating the molecular evolution of mammalian spermatogenesis.

Transcriptome Complexity Disentangled: A Regulatory Molecules Approach

Transcription factors (TFs) and microRNAs (miRNAs) are fundamental regulators of gene expression, cell state, and biological processes. This study investigated whether a small subset of TFs and miRNAs could accurately predict genome-wide gene expression. We analyzed 8895 samples across 31 cancer types from The Cancer Genome Atlas and identified 28 miRNA and 28 TF clusters using unsupervised learning. Medoids of these clusters could differentiate tissues of origin with 92.8% accuracy, demonstrating their biological relevance. We developed Tissue-Agnostic and Tissue-Aware models to predict 20,000 gene expressions using the 56 selected medoid miRNAs and TFs. The Tissue-Aware model attained an R 2 of 0.70 by incorporating tissue-specific information. Despite measuring only 1/400th of the transcriptome, the prediction accuracy was comparable to that achieved by the 1000 landmark genes. This suggests the transcriptome has an intrinsically low-dimensional structure that can be captured by a few regulatory molecules. Our approach could enable cheaper transcriptome assays and analysis of low-quality samples. It also provides insights into genes that are heavily regulated by miRNAs/TFs versus alternative mechanisms. However, model transportability was impacted by dataset discrepancies, especially in miRNA distribution. Overall, this study demonstrates the potential of a biology-guided approach for robust transcriptome representation.

Comparative Analysis of Maternal Gene Expression Patterns Unravels Evolutionary Signatures Across Reproductive Modes

Maternal genes have a pivotal role in regulating metazoan early development. As such their functions have been extensively studied since the dawn of developmental biology. The temporal and spatial dynamics of their transcripts have been thoroughly described in model organisms and their functions have been undergoing heavy investigations. Yet, less is known about the evolutionary changes shaping their presence within diverse oocytes. Due to their unique maternal inheritance pattern, a high degree is predicted to be present when it comes to their expression. Insofar only limited and conflicting results have emerged around it. Here, we set out to elucidate which evolutionary changes could be detected in the maternal gene expression patterns using phylogenetic comparative methods on RNAseq data from 43 species. Using normalized gene expression values and fold change information throughout early development we set out to find the best-fitting evolutionary model. Through modeling, we find evidence supporting both the high degree of divergence and constraint on gene expression values, together with their temporal dynamics. Furthermore, we find that maternal gene expression alone can be used to explain the reproductive modes of different species. Together, these results suggest a highly dynamic evolutionary landscape of maternal gene expression. We also propose a possible functional dichotomy of maternal genes which is influenced by the reproductive strategy undertaken by examined species.

Unlocking Genetic Mysteries during the Epic Sperm Journey toward Fertilization: Further Expanding Cre Mouse Lines

The spatiotemporal expression patterns of genes are crucial for maintaining normal physiological functions in animals. Conditional gene knockout using the cyclization recombination enzyme (Cre)/locus of crossover of P1 (Cre/LoxP) strategy has been extensively employed for functional assays at specific tissue or developmental stages. This approach aids in uncovering the associations between phenotypes and gene regulation while minimizing interference among distinct tissues. Various Cre-engineered mouse models have been utilized in the male reproductive system, including Dppa3-MERCre for primordial germ cells, Ddx4-Cre and Stra8-Cre for spermatogonia, Prm1-Cre and Acrv1-iCre for haploid spermatids, Cyp17a1-iCre for the Leydig cell, Sox9-Cre for the Sertoli cell, and Lcn5/8/9-Cre for differentiated segments of the epididymis. Notably, the specificity and functioning stage of Cre recombinases vary, and the efficiency of recombination driven by Cre depends on endogenous promoters with different sequences as well as the constructed Cre vectors, even when controlled by an identical promoter. Cre mouse models generated via traditional recombination or CRISPR/Cas9 also exhibit distinct knockout properties. This review focuses on Cre-engineered mouse models applied to the male reproductive system, including Cre-targeting strategies, mouse model screening, and practical challenges encountered, particularly with novel mouse strains over the past decade. It aims to provide valuable references for studies conducted on the male reproductive system.

Machine-learning analysis reveals an important role for negative selection in shaping cancer aneuploidy landscapes

BACKGROUND

Aneuploidy, an abnormal number of chromosomes within a cell, is a hallmark of cancer. Patterns of aneuploidy differ across cancers, yet are similar in cancers affecting closely related tissues. The selection pressures underlying aneuploidy patterns are not fully understood, hindering our understanding of cancer development and progression.

RESULTS

Here, we apply interpretable machine learning methods to study tissue-selective aneuploidy patterns. We define 20 types of features corresponding to genomic attributes of chromosome-arms, normal tissues, primary tumors, and cancer cell lines (CCLs), and use them to model gains and losses of chromosome arms in 24 cancer types. To reveal the factors that shape the tissue-specific cancer aneuploidy landscapes, we interpret the machine learning models by estimating the relative contribution of each feature to the models. While confirming known drivers of positive selection, our quantitative analysis highlights the importance of negative selection for shaping aneuploidy landscapes. This is exemplified by tumor suppressor gene density being a better predictor of gain patterns than oncogene density, and vice versa for loss patterns. We also identify the importance of tissue-selective features and demonstrate them experimentally, revealing KLF5 as an important driver for chr13q gain in colon cancer. Further supporting an important role for negative selection in shaping the aneuploidy landscapes, we find compensation by paralogs to be among the top predictors of chromosome arm loss prevalence and demonstrate this relationship for one paralog interaction. Similar factors shape aneuploidy patterns in human CCLs, demonstrating their relevance for aneuploidy research.

CONCLUSIONS

Our quantitative, interpretable machine learning models improve the understanding of the genomic properties that shape cancer aneuploidy landscapes.

SRRM2 splicing factor modulates cell fate in early development

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Adaptive functions of structural variants in human brain development

Quantifying the structural variants (SVs) in nonhuman primates could provide a niche to clarify the genetic backgrounds underlying human-specific traits, but such resource is largely lacking. Here, we report an accurate SV map in a population of 562 rhesus macaques, verified by in-house benchmarks of eight macaque genomes with long-read sequencing and another one with genome assembly. This map indicates stronger selective constrains on inversions at regulatory regions, suggesting a strategy for prioritizing them with the most important functions. Accordingly, we identified 75 human-specific inversions and prioritized them. The top-ranked inversions have substantially shaped the human transcriptome, through their dual effects of reconfiguring the ancestral genomic architecture and introducing regional mutation hotspots at the inverted regions. As a proof of concept, we linked APCDD1, located on one of these inversions and down-regulated specifically in humans, to neuronal maturation and cognitive ability. We thus highlight inversions in shaping the human uniqueness in brain development.

Molecular functions of ANKLE2 and its implications in human disease

Ankyrin repeat and LEM domain-containing 2 (ANKLE2) is a scaffolding protein with established roles in cell division and development, the dysfunction of which is increasingly implicated in human disease. ANKLE2 regulates nuclear envelope disassembly at the onset of mitosis and its reassembly after chromosome segregation. ANKLE2 dysfunction is associated with abnormal nuclear morphology and cell division. It regulates the nuclear envelope by mediating protein-protein interactions with barrier to autointegration factor (BANF1; also known as BAF) and with the kinase and phosphatase that modulate the phosphorylation state of BAF. In brain development, ANKLE2 is crucial for proper asymmetric division of neural progenitor cells. In humans, pathogenic loss-of-function mutations in ANKLE2 are associated with primary congenital microcephaly, a condition in which the brain is not properly developed at birth. ANKLE2 is also linked to other disease pathologies, including congenital Zika syndrome, cancer and tauopathy. Here, we review the molecular roles of ANKLE2 and the recent literature on human diseases caused by its dysfunction.

Genome-Wide Identification of the SPP/SPPL Gene Family and BnaSPPL4 Regulating Male Fertility in Rapeseed (Brassica napus L.)

Signal peptide peptidase (SPP) and its homologs, signal peptide peptidase-like (SPPL) proteases, are members of the GxGD-type aspartyl protease family, which is widespread in plants and animals and is a class of transmembrane proteins with significant biological functions. SPP/SPPLs have been identified; however, the functions of SPP/SPPL in rapeseed (Brassica napus L.) have not been reported. In this study, 26 SPP/SPPLs were identified in rapeseed and categorized into three groups: SPP, SPPL2, and SPPL3. These members mainly contained the Peptidase_A22 and PA domains, which were distributed on 17 out of 19 chromosomes. Evolutionary analyses indicated that BnaSPP/SPPLs evolved with a large number of whole-genome duplication (WGD) events and strong purifying selection. Members are widely expressed and play a key role in the growth and development of rapeseed. The regulation of rapeseed pollen fertility by the BnaSPPL4 gene was further validated through experiments based on bioinformatics analysis, concluding that BnaSPPL4 silencing causes male sterility. Cytological observation showed that male infertility caused by loss of BnaSPPL4 gene function occurs late in the mononucleate stage due to microspore dysplasia.

The journey of a generation: advances and promises in the study of primordial germ cell migration

The germline provides the genetic and non-genetic information that passes from one generation to the next. Given this important role in species propagation, egg and sperm precursors, called primordial germ cells (PGCs), are one of the first cell types specified during embryogenesis. In fact, PGCs form well before the bipotential somatic gonad is specified. This common feature of germline development necessitates that PGCs migrate through many tissues to reach the somatic gonad. During their journey, PGCs must respond to select environmental cues while ignoring others in a dynamically developing embryo. The complex multi-tissue, combinatorial nature of PGC migration is an excellent model for understanding how cells navigate complex environments in vivo. Here, we discuss recent findings on the migratory path, the somatic cells that shepherd PGCs, the guidance cues somatic cells provide, and the PGC response to these cues to reach the gonad and establish the germline pool for future generations. We end by discussing the fate of wayward PGCs that fail to reach the gonad in diverse species. Collectively, this field is poised to yield important insights into emerging reproductive technologies.

A Single-Cell Landscape of Spermioteleosis in Mice and Pigs

(1) Background: Spermatozoa acquired motility and matured in epididymis after production in the testis. However, there is still limited understanding of the specific characteristics of sperm development across different species. In this study, we employed a comprehensive approach to analyze cell compositions in both testicular and epididymal tissues, providing valuable insights into the changes occurring during meiosis and spermiogenesis in mouse and pig models. Additionally, we identified distinct gene expression signatures associated with various spermatogenic cell types. (2) Methods: To investigate the differences in spermatogenesis between mice and pigs, we constructed a single-cell RNA dataset. (3) Results: Our findings revealed notable differences in testicular cell clusters between these two species. Furthermore, distinct gene expression patterns were observed among epithelial cells from different regions of the epididymis. Interestingly, regional gene expression patterns were also identified within principal cell clusters of the mouse epididymis. Moreover, through analysing differentially expressed genes related to the epididymis in both mouse and pig models, we successfully identified potential marker genes associated with sperm development and maturation for each species studied. (4) Conclusions: This research presented a comprehensive single-cell landscape analysis of both testicular and epididymal tissues, shedding light on the intricate processes involved in spermatogenesis and sperm maturation, specifically within mouse and pig models.

RNA-seq analysis reveals changes in mRNA expression during development in Daphnia mitsukuri

Temporal transcriptional variation is a major contributor to functional evolution and the developmental process. Parthenogenetic water fleas of the genus Daphnia (Cladocera) provide an ideal model to characterize gene expression patterns across distinct developmental stages. Herein, we report RNA-seq data for female Daphnia mitsukuri at three developmental stages: the embryo, juvenile (three timepoints) and adult. Comparisons of gene expression patterns among these three developmental stages and weighted gene co-expression network analysis based on expression data across developmental stages identified sets of genes underpinning each of the developmental stages of D. mitsukuri. Specifically, highly expressed genes (HEGs) at the embryonic developmental stage were associated with cell proliferation, ensuring the necessary foundation for subsequent development; HEGs at the juvenile stages were associated with chemosensory perception, visual perception and neurotransmission, allowing individuals to enhance detection of potential environmental risks; HEGs at the adult stage were associated with antioxidative defensive systems, enabling adults to mount an efficient response to perceived environmental risks. Additionally, we found a significant overlap between expanded gene families of Daphnia species and HEGs at the juvenile stages, and these genes were associated with visual perception and neurotransmission. Our work provides a resource of developmental transcriptomes, and comparative analyses that characterize gene expression dynamics throughout development of Daphnia.

betAS: intuitive analysis and visualization of differential alternative splicing using beta distributions

Next-generation RNA sequencing allows alternative splicing (AS) quantification with unprecedented resolution, with the relative inclusion of an alternative sequence in transcripts being commonly quantified by the proportion of reads supporting it as percent spliced-in (PSI). However, PSI values do not incorporate information about precision, proportional to the respective AS events' read coverage. Beta distributions are suitable to quantify inclusion levels of alternative sequences, using reads supporting their inclusion and exclusion as surrogates for the two distribution shape parameters. Each such beta distribution has the PSI as its mean value and is narrower when the read coverage is higher, facilitating the interpretability of its precision when plotted. We herein introduce a computational pipeline, based on beta distributions accurately modeling PSI values and their precision, to quantitatively and visually compare AS between groups of samples. Our methodology includes a differential splicing significance metric that compromises the magnitude of intergroup differences, the estimation uncertainty in individual samples, and the intragroup variability, being therefore suitable for multiple-group comparisons. To make our approach accessible and clear to both noncomputational and computational biologists, we developed betAS, an interactive web app and user-friendly R package for visual and intuitive differential splicing analysis from read count data.

Comparative transcriptome analysis between rhesus macaques ( Macaca mulatta) and crab-eating macaques ( M. fascicularis)

Understanding gene expression variations between species is pivotal for deciphering the evolutionary diversity in phenotypes. Rhesus macaques ( Macaca mulatta, MMU) and crab-eating macaques ( M. fascicularis, MFA) serve as crucial nonhuman primate biomedical models with different phenotypes. To date, however, large-scale comparative transcriptome research between these two species has not yet been fully explored. Here, we conducted systematic comparisons utilizing newly sequenced RNA-seq data from 84 samples (41 MFA samples and 43 MMU samples) encompassing 14 common tissues. Our findings revealed a small fraction of genes (3.7%) with differential expression between the two species, as well as 36.5% of genes with tissue-specific expression in both macaques. Comparison of gene expression between macaques and humans indicated that 22.6% of orthologous genes displayed differential expression in at least two tissues. Moreover, 19.41% of genes that overlapped with macaque-specific structural variants showed differential expression between humans and macaques. Of these, the FAM220A gene exhibited elevated expression in humans compared to macaques due to lineage-specific duplication. In summary, this study presents a large-scale transcriptomic comparison between MMU and MFA and between macaques and humans. The discovery of gene expression variations not only enhances the biomedical utility of macaque models but also contributes to the wider field of primate genomics.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that 'complex' organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this 'drift barrier' theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

Haploinsufficiency of ZFHX3, encoding a key player in neuronal development, causes syndromic intellectual disability

Neurodevelopmental disorders (NDDs) result from impaired development and functioning of the brain. Here, we identify loss-of-function (LoF) variation in ZFHX3 as a cause for syndromic intellectual disability (ID). ZFHX3 is a zinc-finger homeodomain transcription factor involved in various biological processes, including cell differentiation and tumorigenesis. We describe 42 individuals with protein-truncating variants (PTVs) or (partial) deletions of ZFHX3, exhibiting variable intellectual disability and autism spectrum disorder, recurrent facial features, relative short stature, brachydactyly, and, rarely, cleft palate. ZFHX3 LoF associates with a specific methylation profile in whole blood extracted DNA. Nuclear abundance of ZFHX3 increases during human brain development and neuronal differentiation. ZFHX3 was found to interact with the chromatin remodeling BRG1/Brm-associated factor complex and the cleavage and polyadenylation complex, suggesting a function in chromatin remodeling and mRNA processing. Furthermore, ChIP-seq for ZFHX3 revealed that it predominantly binds promoters of genes involved in nervous system development. We conclude that loss-of-function variants in ZFHX3 are a cause of syndromic ID associating with a specific DNA methylation profile.