Co-aggregation with Apolipoprotein E modulates the function of Amyloid-ß in Alzheimer's disease.

Which isoforms of apolipoprotein E (apoE) we inherit determine our risk of developing late-onset Alzheimer's Disease (AD), but the mechanism underlying this link is poorly understood. In particular, the relevance of direct interactions between apoE and amyloid-ß (Aß) remains controversial. Here, single-molecule imaging shows that all isoforms of apoE associate with Aß in the early stages of aggregation and then fall away as fibrillation happens. ApoE-Aß co-aggregates account for ~50% of the mass of diffusible Aß aggregates detected in the frontal cortices of homozygotes with the higher-risk APOE4 gene. We show how dynamic interactions between apoE and Aß tune disease-related functions of Aß aggregates throughout the course of aggregation. Our results connect inherited APOE genotype with the risk of developing AD by demonstrating how, in an isoform- and lipidation-specific way, apoE modulates the aggregation, clearance and toxicity of Aß. Selectively removing non-lipidated apoE4-Aß co-aggregates enhances clearance of toxic Aß by glial cells, and reduces secretion of inflammatory markers and membrane damage, demonstrating a clear path to AD therapeutics.

RPS24 alternative splicing is a marker of cancer progression and epithelial-mesenchymal transition.

Although alternative splicing (AS) is a major mechanism that adds diversity to gene expression patterns, its precise role in generating variability in ribosomal proteins, known as ribosomal heterogeneity, remains unclear. The ribosomal protein S24 (RPS24) gene, encoding a ribosomal component, undergoes AS; however, in-depth studies have been challenging because of three microexons between exons 4 and 6. We conducted a detailed analysis of RPS24 AS isoforms using a direct approach to investigate the splicing junctions related to these microexons, focusing on four AS isoforms. Each of these isoforms showed tissue specificity and relative differences in expression among cancer types. Significant differences in the proportions of these RPS24 AS isoforms between cancerous and normal tissues across diverse cancer types were also observed. Our study highlighted a significant correlation between the expression levels of a specific RPS24 AS isoform and the epithelial-mesenchymal transition process in lung and breast cancers. Our research contributes to a better understanding of the intricate regulatory mechanisms governing AS of ribosomal protein genes and highlights the biological implications of RPS24 AS isoforms in tissue development and tumorigenesis.

Canonical and non-canonical functions of p53 isoforms: potentiating the complexity of tumor development and therapy resistance.

Full-length p53 (p53α) plays a pivotal role in maintaining genomic integrity and preventing tumor development. Over the years, p53 was found to exist in various isoforms, which are generated through alternative splicing, alternative initiation of translation, and internal ribosome entry site. p53 isoforms, either C-terminally altered or N-terminally truncated, exhibit distinct biological roles compared to p53α, and have significant implications for tumor development and therapy resistance. Due to a lack of part and/or complete C- or N-terminal domains, ectopic expression of some p53 isoforms failed to induce expression of canonical transcriptional targets of p53α like CDKN1A or MDM2, even though they may bind their promoters. Yet, p53 isoforms like Δ40p53α still activate subsets of targets including MDM2 and BAX. Furthermore, certain p53 isoforms transactivate even novel targets compared to p53α. More recently, non-canonical functions of p53α in DNA repair and of different isoforms in DNA replication unrelated to transcriptional activities were discovered, amplifying the potential of p53 as a master regulator of physiological and tumor suppressor functions in human cells. Both regarding canonical and non-canonical functions, alternative p53 isoforms frequently exert dominant negative effects on p53α and its partners, which is modified by the relative isoform levels. Underlying mechanisms include hetero-oligomerization, changes in subcellular localization, and aggregation. These processes ultimately influence the net activities of p53α and give rise to diverse cellular outcomes. Biological roles of p53 isoforms have implications for tumor development and cancer therapy resistance. Dysregulated expression of isoforms has been observed in various cancer types and is associated with different clinical outcomes. In conclusion, p53 isoforms have expanded our understanding of the complex regulatory network involving p53 in tumors. Unraveling the mechanisms underlying the biological roles of p53 isoforms provides new avenues for studies aiming at a better understanding of tumor development and developing therapeutic interventions to overcome resistance.

Shaping human brain development and vulnerability through alternative splicing.

Alternative splicing contributes to shaping lineage-specific gene expression and phenotypes. In this issue of Cell Genomics, Recinos, Bao, Wang, et al.1 report that the balance between splicing isoforms of the microtubule-associated protein Tau in the brain is differentially regulated among primates by the RNA-binding protein MBNL2, with consequences for protein aggregation and neurodegeneration in humans.

Alternative isoform expression of key thermogenic genes in human beige adipocytes.

Background: The beneficial effect of thermogenic adipocytes in maintaining body weight and protecting against metabolic disorders has raised interest in understanding the regulatory mechanisms defining white and beige adipocyte identity. Although alternative splicing has been shown to propagate adipose browning signals in mice, this has yet to be thoroughly investigated in human adipocytes. Methods: We performed parallel white and beige adipogenic differentiation using primary adipose stem cells from 6 unrelated healthy subjects and assessed differential gene and isoform expression in mature adipocytes by RNA sequencing. Results: We find 777 exon junctions with robust differential usage between white and beige adipocytes in all 6 subjects, mapping to 562 genes. Importantly, only 10% of these differentially spliced genes are also differentially expressed, indicating that alternative splicing constitutes an additional layer of gene expression regulation during beige adipocyte differentiation. Functional classification of alternative isoforms points to a gain of function for key thermogenic transcription factors such as PPARG and CITED1, and enzymes such as PEMT, or LPIN1. We find that a large majority of the splice variants arise from differential TSS usage, with beige-specific TSSs being enriched for PPARγ and MED1 binding compared to white-specific TSSs. Finally, we validate beige specific isoform expression at the protein level for two thermogenic regulators, PPARγ and PEMT. Discussion: These results suggest that differential isoform expression through alternative TSS usage is an important regulatory mechanism for human adipocyte thermogenic specification.

Ribosome rescue factor PELOTA modulates translation start site choice for C/EBPα protein isoforms.

Translation initiation at alternative start sites can dynamically control the synthesis of two or more functionally distinct protein isoforms from a single mRNA. Alternate isoforms of the developmental transcription factor CCAAT/enhancer-binding protein α (C/EBPα) produced from different start sites exert opposing effects during myeloid cell development. This choice between alternative start sites depends on sequence features of the CEBPA transcript, including a regulatory uORF, but the molecular basis is not fully understood. Here, we identify the factors that affect C/EBPα isoform choice using a sensitive and quantitative two-color fluorescent reporter coupled with CRISPRi screening. Our screen uncovered a role of the ribosome rescue factor PELOTA (PELO) in promoting the expression of the longer C/EBPα isoform by directly removing inhibitory unrecycled ribosomes and through indirect effects mediated by the mechanistic target of rapamycin kinase. Our work uncovers further links between ribosome recycling and translation reinitiation that regulate a key transcription factor, with implications for normal hematopoiesis and leukemogenesis.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differential isoform in skin tissues of different hair length Yak.

BACKGROUND: The hair follicle development process is regulated by sophisticated genes and signaling networks, and the hair grows from the hair follicle. The Tianzhu white yak population exhibits differences in hair length, especially on the forehead and shoulder region. However, the genetic mechanism is still unclear. Isoform sequencing (Iso-seq) technology with advantages in long reads sequencing. Hence, we combined the Iso-seq and RNA-seq methods to investigate the transcript complexity and difference between long-haired yak (LHY) and normal-haired yak (NHY). RESULTS: The hair length measurement result showed a significant difference between LHY and NHY on the forehead and the shoulder (P-value < 0.001). The skin samples from the forehead and the shoulder of LHY and NHY were pooled for isoform sequencing (Iso-seq). We obtained numerous long transcripts, including novel isoforms, long non-coding RNA, alternative splicing events, and alternative polyadenylation events. Combined with RNA-seq data, we performed differential isoforms (DEIs) analysis between LHY and NHY. We found that some hair follicle and skin development-related DEIs, like BMP4, KRT2, IGF2R, and COL1A2 in the forehead skin; BMP1, KRT1, FGF5, COL2A1, and IGFBP5 in the shoulder skin. Enrichment analysis revealed that DEIs in both two comparable groups significantly participated in skin and hair follicle development-related pathways, like ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathways. The results indicated that the hair follicle development of Tianzhu white yak may influence the hair length difference. Besides, the protein-protein interaction (PPI) network of DEIs showed COL2A1 and COL3A1 exhibited a high degree of centrality, and these two genes were suggested as potential candidates for the hair length growth of Tianzhu white yak. CONCLUSIONS: The results provided a comprehensive analysis of the transcriptome complexity and identified differential transcripts that enhance our understanding of the molecular mechanisms underlying the variation in hair length growth in Tianzhu white yak.

Collagen XVIII regulates extracellular matrix integrity in the developing nephrons and impacts nephron progenitor cell behavior.

Renal development is a complex process in which two major processes, tubular branching and nephron development, regulate each other reciprocally. Our previous findings have indicated that collagen XVIII (ColXVIII), an extracellular matrix protein, affects the renal branching morphogenesis. We investigate here the role of ColXVIII in nephron formation and the behavior of nephron progenitor cells (NPCs) using isoform-specific ColXVIII knockout mice. The results show that the short ColXVIII isoform predominates in the early epithelialized nephron structures whereas the two longer isoforms are expressed only in the later phases of glomerular formation. Meanwhile, electron microscopy showed that the ColXVIII mutant embryonic kidneys have ultrastructural defects at least from embryonic day 16.5 onwards. Similar structural defects had previously been observed in adult ColXVIII-deficient mice, indicating a congenital origin. The lack of ColXVIII led to a reduced NPC population in which changes in NPC proliferation and maintenance and in macrophage influx were perceived to play a role. The changes in NPC behavior in turn led to notably reduced overall nephron formation. In conclusion, the results show that ColXVIII has multiple roles in renal development, both in ureteric branching and in NPC behavior.

Molecular and biochemical characterizations of a Fasciola gigantica retinoid X receptor-α isoform A (FgRXRα-A).

Fascioliasis is a parasitic infection in animals and humans caused by the parasitic flatworm genus Fasciola, which has two major species, F. hepatica and F. gigantica. A major concern regarding this disease is drug resistance, which is increasingly reported worldwide. Hence, the discovery of a novel drug as well as drug targets is crucially required. Therefore, this study aims to characterize the novel drug target in the adult F. gigantica. In the beginning, we hypothesized that the parasite might interact with some host molecules when it lives inside the liver parenchyma or bile ducts, specifically hormones and hormone-like molecules, through the specific receptors, primarily nuclear receptors (NRs), which are recognized as a major drug target in various diseases. The retinoid X receptor (RXR) is a member of subfamily 2 NRs that plays multitudinous roles in organisms by forming homodimers or heterodimers with other NRs. We obtained the full-length amino acid sequences of F. gigantica retinoid X receptor-alpha (FgRXRα-A) from the transcriptome of F. gigantica that existed in the NCBI database. The FgRXRα-A were computationally predicted for the basic properties, multiple aligned, phylogeny analyzed, and generated of 2D and 3D models. Moreover, FgRXRα-A was molecular cloned and expressed as a recombinant protein (rFgRXRα-A), then used for immunization for specific polyclonal antibodies. The native FgRXRα-A was detected in the parasite extracts and tissues, and the function was investigated by in vitro binding assay. The results demonstrated the conservation of FgRXRα-A to the other RXRs, especially RXRs from the trematodes. Interestingly, the native FgRXRα-A could be detected in the testes of the parasite, where the sex hormones are accumulated. Moreover, the binding assay revealed the interaction of 9-cis retinoic acid and FgRXRα-A, suggesting the function of FgRXRα-A. Our findings suggested that FgRXRα-A will be involved with the sexual reproduction of the parasite by forming heterodimers with other NRs, and it could be the potential target for further drug development of fascioliasis.

Alzheimer's disease associated isoforms of human CD33 distinctively modulate microglial cell responses in 5XFAD mice.

Microglia play diverse pathophysiological roles in Alzheimer's disease (AD), with genetic susceptibility factors skewing microglial cell function to influence AD risk. CD33 is an immunomodulatory receptor associated with AD susceptibility through a single nucleotide polymorphism that modulates mRNA splicing, skewing protein expression from a long protein isoform (CD33M) to a short isoform (CD33m). Understanding how human CD33 isoforms differentially impact microglial cell function in vivo has been challenging due to functional divergence of CD33 between mice and humans. We address this challenge by studying transgenic mice expressing either of the human CD33 isoforms crossed with the 5XFAD mouse model of amyloidosis and find that human CD33 isoforms have opposing effects on the response of microglia to amyloid-ß (Aß) deposition. Mice expressing CD33M have increased Aß levels, more diffuse plaques, fewer disease-associated microglia, and more dystrophic neurites compared to 5XFAD control mice. Conversely, CD33m promotes plaque compaction and microglia-plaque contacts, and minimizes neuritic plaque pathology, highlighting an AD protective role for this isoform. Protective phenotypes driven by CD33m are detected at an earlier timepoint compared to the more aggressive pathology in CD33M mice that appears at a later timepoint, suggesting that CD33m has a more prominent impact on microglia cell function at earlier stages of disease progression. In addition to divergent roles in modulating phagocytosis, scRNAseq and proteomics analyses demonstrate that CD33m+ microglia upregulate nestin, an intermediate filament involved in cell migration, at plaque contact sites. Overall, our work provides new functional insights into how CD33, as a top genetic susceptibility factor for AD, modulates microglial cell function.

Targeting an Old Foe for Cancer: A Molecular Dynamics Perspective to Unravel the Specific Binding Nature of 2-Methoxy Estradiol to Human ß-Tubulin Isotypes.

Microtubules, composed of α- and ß-tubulin subunits are crucial for cell division with their dynamic tissue-specificity which is dictated by expression of isotypes. These isotypes differ in carboxy-terminal tails (CTTs), rich in negatively charged acidic residues in addition to the differences in the composition of active site residues. 2-Methoxy estradiol (2-ME) is the first antimicrotubule agent that showed less affinity toward hemopoietic-specific ß1 isotype consequently preventing myelosuppression toxicity. The present study focuses on the MD-directed conformational analysis of 2-ME and estimation of its binding affinity in the colchicine binding pocket of various ß-tubulin isotypes combined with the α-tubulin isotype, α1B. AlphaFold 2.0 was used to predict the 3D structure of phylogenetically divergent human ß-tubulin isotypes in dimer form with α1B. The dimeric complexes were subjected to induced-fit docking with 2-ME. The statistical analysis of docking showed differences in the binding characteristics of 2-ME with different isotypes. The replicas of atom-based molecular dynamic simulations of the best conformation of 2-ME provided insights into the molecular-level details of its binding pattern across the isotypes. Furthermore, the MM/GBSA analyses revealed the specific binding energy profile of 2-ME in ß-tubulin isotypes. It also highlighed, 2-ME exhibits the lowest binding affinity toward the ß1 isotype as supported by experimental study. The present study may offer useful information for designing next-generation antimicrotubule agents that are more specific and less toxic.

Dichotomous transactivation domains contribute to growth inhibitory and promotion functions of TAp73.

The transcription factor p73, a member of the p53 tumor-suppressor family, regulates cell death and also supports tumorigenesis, although the mechanistic basis for the dichotomous functions is poorly understood. We report here the identification of an alternate transactivation domain (TAD) located at the extreme carboxyl (C) terminus of TAp73ß, a commonly expressed p73 isoform. Mutational disruption of this TAD significantly reduced TAp73ß's transactivation activity, to a level observed when the amino (N)-TAD that is similar to p53's TAD, is mutated. Mutation of both TADs almost completely abolished TAp73ß's transactivation activity. Expression profiling highlighted a unique set of targets involved in extracellular matrix-receptor interaction and focal adhesion regulated by the C-TAD, resulting in FAK phosphorylation, distinct from the N-TAD targets that are common to p53 and are involved in growth inhibition. Interestingly, the C-TAD targets are also regulated by the oncogenic, amino-terminal-deficient DNp73ß isoform. Consistently, mutation of C-TAD reduces cellular migration and proliferation. Mechanistically, selective binding of TAp73ß to DNAJA1 is required for the transactivation of C-TAD target genes, and silencing DNAJA1 expression abrogated all C-TAD-mediated effects. Taken together, our results provide a mechanistic basis for the dichotomous functions of TAp73 in the regulation of cellular growth through its distinct TADs.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms.

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

Emerging concepts of myosin 18A isoform mechanobiology in organismal and immune system physiology, development, and function.

Alternative and combinatorial splicing of myosin 18A (MYO18A) gene transcripts results in expression of MYO18A protein isoforms and isoform variants with different membrane and subcellular localizations, and functional properties. MYO18A proteins are members of the myosin superfamily consisting of a myosin-like motor domain, an IQ motif, and a coiled-coil domain. MYO18A isoforms, however, lack the ability to hydrolyze ATP and do not perform ATP-dependent motor activity. MYO18A isoforms are distinguished by different amino- and carboxy-terminal extensions and domains. The domain organization and functions of MYO18Aα, MYO18Aß, and MYO18Aγ have been studied experimentally. MYO18Aα and MYO18Aß have a common carboxy-terminal extension but differ by the presence or absence of an amino-terminal KE repeat and PDZ domain, respectively. The amino- and carboxy-terminal extensions of MYO18Aγ contain unique proline and serine-rich domains. Computationally predicted MYO18Aε and MYO18Aδ isoforms contain the carboxy-terminal serine-rich extension but differ by the presence or absence of the amino-terminal KE/PDZ extension. Additional isoform variants within each category arise by alternative utilization or inclusion/exclusion of small exons. MYO18Aα variants are expressed in somatic cells and mature immune cells, whereas MYO18Aß variants occur mainly in myeloid and natural killer cells. MYO18Aγ expression is selective to cardiac and skeletal muscle. In the present review perspective, we discuss current and emerging concepts of the functional specialization of MYO18A proteins in membrane and cytoskeletal dynamics, cellular communication and signaling, endocytic and exocytic organelle movement, viral infection, and as the SP-R210 receptor for surfactant protein A.

Selective Cell-Surface Expression of Triheteromeric NMDA Receptors.

NMDA-type ionotropic glutamate receptors are critically involved in many brain functions and are implicated in a variety of brain disorders. Seven NMDA receptor subunits exist (GluN1, GluN2A-D, and GluN3A-B) that assemble into tetrameric receptor subtypes with distinct functional properties and physiological roles. The majority NMDA receptors are composed of two GluN1 and two GluN2 subunits, which can assemble into four diheteromeric receptors subtypes composed of GluN1 and one type of GluN2 subunit (e.g., GluN1/2A), and presumably also six triheteromeric receptor subtypes composed of GluN1 and two different GluN2 subunits (e.g., GluN1/2A/2B). Furthermore, the GluN1 subunit exists as eight splice variants (e.g., GluN1-1a and GluN1-1b isoforms), and two different GluN1 isoforms can co-assemble to also form triheteromeric NMDA receptors (e.g., GluN1-1a/1b/2A). Here, we describe a method to faithfully express triheteromeric NMDA receptors in heterologous expression systems by controlling the identity of two of the four subunits. This method overcomes the problem that co-expression of three different NMDA receptor subunits generates two distinct diheteromeric receptor subtypes as well as one triheteromeric receptor subtype, thereby confounding studies that require a homogenous population of triheteromeric NMDA receptors. The method has been applied to selectively express recombinant triheteromeric GluN1/2A/2B, GluN1/2A/2C, GluN1/2B/2D, GluN1-1a/GluN1-1b/2A, GluN1-1a/GluN1-1b/2B receptors with negligible co-expression of the respective diheteromeric receptor subtypes. This method therefore enables quantitative evaluation of functional and pharmacological properties of triheteromeric NMDA receptors, some of which are abundant NMDA receptor subtypes in the adult brain.

A ZBP1 isoform blocks ZBP1-mediated cell death.

ZBP1 is an interferon (IFN)-induced nucleic acid (NA) sensor that senses unusual Z-form NA (Z-NA) to promote cell death and inflammation. However, the mechanisms that dampen ZBP1 activation to fine-tune inflammatory responses are unclear. Here, we characterize a short isoform of ZBP1 (referred to as ZBP1-S) as an intrinsic suppressor of the inflammatory signaling mediated by full-length ZBP1. Mechanistically, ZBP1-S depresses ZBP1-mediated cell death by competitive binding with Z-NA for Zα domains of ZBP1. Cells from mice (Ripk1D325A/D325A) with cleavage-resistant RIPK1-induced autoinflammatory (CRIA) syndrome are alive but sensitive to IFN-induced and ZBP1-dependent cell death. Intriguingly, Ripk1D325A/D325A cells die spontaneously when ZBP1-S is deleted, indicating that cell death driven by ZBP1 is under the control of ZBP1-S. Thus, our findings reveal that alternative splicing of Zbp1 represents autogenic inhibition for regulating ZBP1 signaling and indicate that uncoupling of Z-NA with ZBP1 could be an effective strategy against autoinflammations.

Transcription factor 4 promotes increased corneal endothelial cellular migration by altering microtubules in Fuchs endothelial corneal dystrophy.

Fuchs endothelial corneal dystrophy (FECD) is a complex corneal disease characterized by the progressive decline and morphological changes of corneal endothelial cells (CECs) that leads to corneal edema and vision loss. The most common mutation in FECD is an intronic CTG repeat expansion in transcription factor 4 (TCF4) that leads to its altered expression. Corneal endothelial wound healing occurs primarily through cell enlargement and migration, and FECD CECs have been shown to display increased migration speeds. In this study, we aim to determine whether TCF4 can promote cellular migration in FECD CECs. We generated stable CEC lines derived from FECD patients that overexpressed different TCF4 isoforms and investigated epithelial-to-mesenchymal (EMT) expression, morphological analysis and cellular migration speeds. We found that full length TCF4-B isoform overexpression promotes cellular migration in FECD CECs in an EMT-independent manner. RNA-sequencing identified several pathways including the negative regulation of microtubules, with TUBB4A (tubulin beta 4A class IVa) as the top upregulated gene. TUBB4A expression was increased in FECD ex vivo specimens, and there was altered expression of cytoskeleton proteins, tubulin and actin, compared to normal healthy donor ex vivo specimens. Additionally, there was increased acetylation and detyrosination of microtubules in FECD supporting that microtubule stability is altered in FECD and could promote cellular migration. Future studies could be aimed at investigating if targeting the cytoskeleton and microtubules would have therapeutic potential for FECD by promoting cellular migration and regeneration.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms.

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

Gap junctions in Turing-type periodic feather pattern formation.

Periodic patterning requires coordinated cell-cell interactions at the tissue level. Turing showed, using mathematical modeling, how spatial patterns could arise from the reactions of a diffusive activator-inhibitor pair in an initially homogeneous 2D field. Most activators and inhibitors studied in biological systems are proteins, and the roles of cell-cell interaction, ions, bioelectricity, etc. are only now being identified. Gap junctions (GJs) mediate direct exchanges of ions or small molecules between cells, enabling rapid long-distance communications in a cell collective. They are therefore good candidates for propagating nonprotein-based patterning signals that may act according to the Turing principles. Here, we explore the possible roles of GJs in Turing-type patterning using feather pattern formation as a model. We found 7 of the 12 investigated GJ isoforms are highly dynamically expressed in the developing chicken skin. In ovo functional perturbations of the GJ isoform, connexin 30, by siRNA and the dominant-negative mutant applied before placode development led to disrupted primary feather bud formation. Interestingly, inhibition of gap junctional intercellular communication (GJIC) in the ex vivo skin explant culture allowed the sequential emergence of new feather buds at specific spatial locations relative to the existing primary buds. The results suggest that GJIC may facilitate the propagation of long-distance inhibitory signals. Thus, inhibition of GJs may stimulate Turing-type periodic feather pattern formation during chick skin development, and the removal of GJ activity would enable the emergence of new feather buds if the local environment were competent and the threshold to form buds was reached. We further propose Turing-based computational simulations that can predict the sequential appearance of these ectopic buds. Our models demonstrate how a Turing activator-inhibitor system can continue to generate patterns in the competent morphogenetic field when the level of intercellular communication at the tissue scale is modulated.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain.

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.