A multi-stem cell basis for craniosynostosis and calvarial mineralization.

Craniosynostosis is a group of disorders of premature calvarial suture fusion. The identity of the calvarial stem cells (CSCs) that produce fusion-driving osteoblasts in craniosynostosis remains poorly understood. Here we show that both physiologic calvarial mineralization and pathologic calvarial fusion in craniosynostosis reflect the interaction of two separate stem cell lineages; a previously identified cathepsin K (CTSK) lineage CSC1 (CTSK+ CSC) and a separate discoidin domain-containing receptor 2 (DDR2) lineage stem cell (DDR2+ CSC) that we identified in this study. Deletion of Twist1, a gene associated with craniosynostosis in humans2,3, solely in CTSK+ CSCs is sufficient to drive craniosynostosis in mice, but the sites that are destined to fuse exhibit an unexpected depletion of CTSK+ CSCs and a corresponding expansion of DDR2+ CSCs, with DDR2+ CSC expansion being a direct maladaptive response to CTSK+ CSC depletion. DDR2+ CSCs display full stemness features, and our results establish the presence of two distinct stem cell lineages in the sutures, with both populations contributing to physiologic calvarial mineralization. DDR2+ CSCs mediate a distinct form of endochondral ossification without the typical haematopoietic marrow formation. Implantation of DDR2+ CSCs into suture sites is sufficient to induce fusion, and this phenotype was prevented by co-transplantation of CTSK+ CSCs. Finally, the human counterparts of DDR2+ CSCs and CTSK+ CSCs display conserved functional properties in xenograft assays. The interaction between these two stem cell populations provides a new biologic interface for the modulation of calvarial mineralization and suture patency.

Identification of SARS-CoV-2 inhibitors using lung and colonic organoids.

There is an urgent need to create novel models using human disease-relevant cells to study severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) biology and to facilitate drug screening. Here, as SARS-CoV-2 primarily infects the respiratory tract, we developed a lung organoid model using human pluripotent stem cells (hPSC-LOs). The hPSC-LOs (particularly alveolar type-II-like cells) are permissive to SARS-CoV-2 infection, and showed robust induction of chemokines following SARS-CoV-2 infection, similar to what is seen in patients with COVID-19. Nearly 25% of these patients also have gastrointestinal manifestations, which are associated with worse COVID-19 outcomes1. We therefore also generated complementary hPSC-derived colonic organoids (hPSC-COs) to explore the response of colonic cells to SARS-CoV-2 infection. We found that multiple colonic cell types, especially enterocytes, express ACE2 and are permissive to SARS-CoV-2 infection. Using hPSC-LOs, we performed a high-throughput screen of drugs approved by the FDA (US Food and Drug Administration) and identified entry inhibitors of SARS-CoV-2, including imatinib, mycophenolic acid and quinacrine dihydrochloride. Treatment at physiologically relevant levels of these drugs significantly inhibited SARS-CoV-2 infection of both hPSC-LOs and hPSC-COs. Together, these data demonstrate that hPSC-LOs and hPSC-COs infected by SARS-CoV-2 can serve as disease models to study SARS-CoV-2 infection and provide a valuable resource for drug screening to identify candidate COVID-19 therapeutics.

Adaptable haemodynamic endothelial cells for organogenesis and tumorigenesis.

Endothelial cells adopt tissue-specific characteristics to instruct organ development and regeneration1,2. This adaptability is lost in cultured adult endothelial cells, which do not vascularize tissues in an organotypic manner. Here, we show that transient reactivation of the embryonic-restricted ETS variant transcription factor 2 (ETV2)3 in mature human endothelial cells cultured in a serum-free three-dimensional matrix composed of a mixture of laminin, entactin and type-IV collagen (LEC matrix) 'resets' these endothelial cells to adaptable, vasculogenic cells, which form perfusable and plastic vascular plexi. Through chromatin remodelling, ETV2 induces tubulogenic pathways, including the activation of RAP1, which promotes the formation of durable lumens4,5. In three-dimensional matrices-which do not have the constraints of bioprinted scaffolds-the 'reset' vascular endothelial cells (R-VECs) self-assemble into stable, multilayered and branching vascular networks within scalable microfluidic chambers, which are capable of transporting human blood. In vivo, R-VECs implanted subcutaneously in mice self-organize into durable pericyte-coated vessels that functionally anastomose to the host circulation and exhibit long-lasting patterning, with no evidence of malformations or angiomas. R-VECs directly interact with cells within three-dimensional co-cultured organoids, removing the need for the restrictive synthetic semipermeable membranes that are required for organ-on-chip systems, therefore providing a physiological platform for vascularization, which we call 'Organ-On-VascularNet'. R-VECs enable perfusion of glucose-responsive insulin-secreting human pancreatic islets, vascularize decellularized rat intestines and arborize healthy or cancerous human colon organoids. Using single-cell RNA sequencing and epigenetic profiling, we demonstrate that R-VECs establish an adaptive vascular niche that differentially adjusts and conforms to organoids and tumoroids in a tissue-specific manner. Our Organ-On-VascularNet model will permit metabolic, immunological and physiochemical studies and screens to decipher the crosstalk between organotypic endothelial cells and parenchymal cells for identification of determinants of endothelial cell heterogeneity, and could lead to advances in therapeutic organ repair and tumour targeting.

Molecular profiling of the stroke-induced alterations in the cerebral microvasculature reveals promising therapeutic candidates.

Stroke-induced cerebral microvascular dysfunction contributes to aggravation of neuronal injury and compromises the efficacy of current reperfusion therapies. Understanding the molecular alterations in cerebral microvessels in stroke will provide original opportunities for scientific investigation of novel therapeutic strategies. Toward this goal, using a recently optimized method which minimizes cell activation and preserves endothelial cell interactions and RNA integrity, we conducted a genome-wide transcriptomic analysis of cerebral microvessels in a mouse model of stroke and compared these transcriptomic alterations with the ones observed in human, nonfatal, brain stroke lesions. Results from these unbiased comparative analyses have revealed the common alterations in mouse stroke microvessels and human stroke lesions and identified shared molecular features associated with vascular disease (e.g., Serpine1/Plasminogen Activator Inhibitor-1, Hemoxygenase-1), endothelial activation (e.g., Angiopoietin-2), and alterations in sphingolipid metabolism and signaling (e.g., Sphigosine-1-Phosphate Receptor 2). Sphingolipid profiling of mouse cerebral microvessels validated the transcript data and revealed the enrichment of sphingomyelin and sphingoid species in the cerebral microvasculature compared to brain and the stroke-induced increase in ceramide species. In summary, our study has identified novel molecular alterations in several microvessel-enriched, translationally relevant, and druggable targets, which are potent modulators of endothelial function. Our comparative analyses have revealed the presence of molecular features associated with cerebral microvascular dysfunction in human chronic stroke lesions. The results shared here provide a detailed resource for therapeutic discovery of candidates for neurovascular protection in stroke and potentially, other pathologies exhibiting cerebral microvascular dysfunction.

Checkpoint kinase 2 controls insulin secretion and glucose homeostasis.

After the discovery of insulin, a century ago, extensive work has been done to unravel the molecular network regulating insulin secretion. Here we performed a chemical screen and identified AZD7762, a compound that potentiates glucose-stimulated insulin secretion (GSIS) of a human ß cell line, healthy and type 2 diabetic (T2D) human islets and primary cynomolgus macaque islets. In vivo studies in diabetic mouse models and cynomolgus macaques demonstrated that AZD7762 enhances GSIS and improves glucose tolerance. Furthermore, genetic manipulation confirmed that ablation of CHEK2 in human ß cells results in increased insulin secretion. Consistently, high-fat-diet-fed Chk2-/- mice show elevated insulin secretion and improved glucose clearance. Finally, untargeted metabolic profiling demonstrated the key role of the CHEK2-PP2A-PLK1-G6PD-PPP pathway in insulin secretion. This study successfully identifies a previously unknown insulin secretion regulating pathway that is conserved across rodents, cynomolgus macaques and human ß cells in both healthy and T2D conditions.

Hierarchical functional differences between gyri and sulci at different scales.

Gyri and sulci are 2 fundamental cortical folding patterns of the human brain. Recent studies have suggested that gyri and sulci may play different functional roles given their structural and functional heterogeneity. However, our understanding of the functional differences between gyri and sulci remains limited due to several factors. Firstly, previous studies have typically focused on either the spatial or temporal domain, neglecting the inherently spatiotemporal nature of brain functions. Secondly, analyses have often been restricted to either local or global scales, leaving the question of hierarchical functional differences unresolved. Lastly, there has been a lack of appropriate analytical tools for interpreting the hierarchical spatiotemporal features that could provide insights into these differences. To overcome these limitations, in this paper, we proposed a novel hierarchical interpretable autoencoder (HIAE) to explore the hierarchical functional difference between gyri and sulci. Central to our approach is its capability to extract hierarchical features via a deep convolutional autoencoder and then to map these features into an embedding vector using a carefully designed feature interpreter. This process transforms the features into interpretable spatiotemporal patterns, which are pivotal in investigating the functional disparities between gyri and sulci. We evaluate the proposed framework on Human Connectome Project task functional magnetic resonance imaging dataset. The experiments demonstrate that the HIAE model can effectively extract and interpret hierarchical spatiotemporal features that are neuroscientifically meaningful. The analyses based on the interpreted features suggest that gyri are more globally activated, whereas sulci are more locally activated, demonstrating a distinct transition in activation patterns as the scale shifts from local to global. Overall, our study provides novel insights into the brain's anatomy-function relationship.

P62 promotes FSH-induced antral follicle formation by directing degradation of ubiquitinated WT1.

In females, the pathophysiological mechanism of poor ovarian response (POR) is not fully understood. Considering the expression level of p62 was significantly reduced in the granulosa cells (GCs) of POR patients, this study focused on identifying the role of the selective autophagy receptor p62 in conducting the effect of follicle-stimulating hormone (FSH) on antral follicles (AFs) formation in female mice. The results showed that p62 in GCs was FSH responsive and that its level increased to a peak and then decreased time-dependently either in ovaries or in GCs after gonadotropin induction in vivo. GC-specific deletion of p62 resulted in subfertility, a significantly reduced number of AFs and irregular estrous cycles, which were same as pathophysiological symptom of POR. By conducting mass spectrum analysis, we found the ubiquitination of proteins was decreased, and autophagic flux was blocked in GCs. Specifically, the level of nonubiquitinated Wilms tumor 1 homolog (WT1), a transcription factor and negative controller of GC differentiation, increased steadily. Co-IP results showed that p62 deletion increased the level of ubiquitin-specific peptidase 5 (USP5), which blocked the ubiquitination of WT1. Furthermore, a joint analysis of RNA-seq and the spatial transcriptome sequencing data showed the expression of steroid metabolic genes and FSH receptors pivotal for GCs differentiation decreased unanimously. Accordingly, the accumulation of WT1 in GCs deficient of p62 decreased steroid hormone levels and reduced FSH responsiveness, while the availability of p62 in GCs simultaneously ensured the degradation of WT1 through the ubiquitin‒proteasome system and autophagolysosomal system. Therefore, p62 in GCs participates in GC differentiation and AF formation in FSH induction by dynamically controlling the degradation of WT1. The findings of the study contributes to further study the pathology of POR.

MEN1 is a regulator of alternative splicing and prevents R-loop-induced genome instability through suppression of RNA polymerase II elongation.

The fidelity of alternative splicing (AS) patterns is essential for growth development and cell fate determination. However, the scope of the molecular switches that regulate AS remains largely unexplored. Here we show that MEN1 is a previously unknown splicing regulatory factor. MEN1 deletion resulted in reprogramming of AS patterns in mouse lung tissue and human lung cancer cells, suggesting that MEN1 has a general function in regulating alternative precursor mRNA splicing. MEN1 altered exon skipping and the abundance of mRNA splicing isoforms of certain genes with suboptimal splice sites. Chromatin immunoprecipitation and chromosome walking assays revealed that MEN1 favored the accumulation of RNA polymerase II (Pol II) in regions encoding variant exons. Our data suggest that MEN1 regulates AS by slowing the Pol II elongation rate and that defects in these processes trigger R-loop formation, DNA damage accumulation and genome instability. Furthermore, we identified 28 MEN1-regulated exon-skipping events in lung cancer cells that were closely correlated with survival in patients with lung adenocarcinoma, and MEN1 deficiency sensitized lung cancer cells to splicing inhibitors. Collectively, these findings led to the identification of a novel biological role for menin in maintaining AS homeostasis and link this role to the regulation of cancer cell behavior.

ROCK1 is a multifunctional factor maintaining the primordial follicle reserve and follicular development in mice.

The follicle is the basic structural and functional unit of the ovary in female mammals. The excessive depletion of follicles will lead to diminished ovarian reserve or even premature ovarian failure, resulting in diminished ovarian oogenesis and endocrine function. Excessive follicular depletion is mainly due to loss of primordial follicles. Our analysis of published human ovarian single-cell sequencing results by others revealed a significant increase in rho-associated protein kinase 1 (ROCK1) expression during primordial follicle development. However, the role of ROCK1 in primordial follicle development and maintenance is not clear. This study revealed a gradual increase in ROCK1 expression during primordial follicle activation. Inhibition of ROCK1 resulted in reduced primordial follicle activation, decreased follicular reserve, and delayed development of growing follicles. This effect may be achieved through the HIPPO pathway. The present study indicates that ROCK1 is a key molecule for primordial follicular reserve and follicular development.NEW & NOTEWORTHY ROCK1, one of the Rho GTPases, plays an important role in primordial follicle reserve and follicular development. ROCK1 was primarily expressed in the cytoplasm of oocytes and granulosa cell in mice. Inhibition of ROCK1 significantly reduced the primordial follicle reserve and delayed growing follicle development. ROCK1 regulates primordial follicular reserve and follicle development through the HIPPO signaling pathway. These findings shed new lights on the physiology of sustaining female reproduction.

Pregranulosa cell-derived FGF23 protects oocytes from premature apoptosis during primordial follicle formation by inhibiting p38 MAPK in mice.

A large number of oocytes in the perinatal ovary in rodents get lost for unknown reasons. The granulosa cell-oocyte mutual communication is pivotal for directing formation of the primordial follicle; however, little is known if paracrine factors participate in modulating programmed oocyte death perinatally. We report here that pregranulosa cell-derived fibroblast growth factor 23 (FGF23) functioned in preventing oocyte apoptosis in the perinatal mouse ovary. Our results showed that FGF23 was exclusively expressed in pregranulosa cells, while fibroblast growth factor receptors (FGFRs) were specifically expressed in the oocytes in perinatal ovaries. FGFR1 was one of the representative receptors in mediating FGF23 signaling during the formation of the primordial follicle. In cultured ovaries, the number of live oocytes declines significantly, accompanied by the activation of the p38 mitogen-activated protein kinase signaling pathway, under the condition of FGFR1 disruption by specific inhibitors of FGFR1 or silencing of Fgf23. As a result, oocyte apoptosis increased and eventually led to a decrease in the number of germ cells in perinatal ovaries following the treatments. In the perinatal mouse ovary, pregranulosa cell-derived FGF23 binds to FGFR1 and activates at least the p38 mitogen-activated protein kinase signaling pathway, thereby regulating the level of apoptosis during primordial follicle formation. This study reemphasizes the importance of granulosa cell-oocyte mutual communication in modulating primordial follicle formation and supporting oocyte survival under physiological conditions.

Structural Insights into the Reaction between Hydrogen Peroxide and Di-iron Complexes at the Ferroxidase Center of Ferritin.

The Fe(II) oxidation mechanism in the ferroxidase center of heavy chain ferritin has been studied extensively. However, the actual production of H2O2 was found to be substantially lower than expected at low flux of Fe(II) to ferritin subunits. Here, we demonstrated that H2O2 could interact with the di-iron nuclear center, leading to the production of hydroxyl radicals and oxygen. Two reaction intermediates were captured in the ferroxidase center by using the time-lapse crystallographic techniques in a shellfish ferritin. The crystal structures revealed the binding of H2O2 as a µ -1,2-peroxo-diferric species and the binding of O2 to the diferric structure. This investigation sheds light on the reaction between the di-iron nuclear center and H2O2 and provides insights for the exploitation of metalloenzymes.

Metal-Binding Protein TaGlo1 Improves Fungal Resistance to Arsenite (AsIII) and Methylarsenite (MAsIII) in Paddy Soil.

Trivalent arsenicals such as arsenite (AsIII) and methylarsenite (MAsIII) are thought to be ubiquitous in flooded paddy soils and have higher toxicity than pentavalent forms. Fungi are widely prevalent in the rice rhizosphere, and the latter is considered a hotspot for As uptake. However, few studies have focused on alleviating As toxicity in paddy soils using fungi. In this study, we investigated the mechanism by which the protein TaGlo1, derived from the As-resistant fungal strain Trichoderma asperellum SM-12F1, mitigates AsIII and MAsIII toxicity in paddy soils. Taglo1 gene expression in Escherichia coli BL21 conferred strong resistance to AsIII and MAsIII, while purified TaGlo1 showed a high affinity for AsIII and MAsIII. Three cysteine residues (Cys13, Cys18, and Cys71) play crucial roles in binding with AsIII, while only two (Cys13 and Cys18) play crucial roles for MAsIII binding. TaGlo1 had a stronger binding strength for MAsIII than AsIII. Importantly, up to 90.2% of the homologous TaGlo1 proteins originate from fungi by GenBank searching. In the rhizospheres of 14 Chinese paddy soils, Taglo1 was widely distributed and its gene abundance increased with porewater As. This study highlights the potential of fungi to mitigate As toxicity and availability in the soil-rice continuum and suggests future microbial strategies for bioremediation.

Discovery of a periosteal stem cell mediating intramembranous bone formation.

Bone consists of separate inner endosteal and outer periosteal compartments, each with distinct contributions to bone physiology and each maintaining separate pools of cells owing to physical separation by the bone cortex. The skeletal stem cell that gives rise to endosteal osteoblasts has been extensively studied; however, the identity of periosteal stem cells remains unclear1-5. Here we identify a periosteal stem cell (PSC) that is present in the long bones and calvarium of mice, displays clonal multipotency and self-renewal, and sits at the apex of a differentiation hierarchy. Single-cell and bulk transcriptional profiling show that PSCs display transcriptional signatures that are distinct from those of other skeletal stem cells and mature mesenchymal cells. Whereas other skeletal stem cells form bone via an initial cartilage template using the endochondral pathway4, PSCs form bone via a direct intramembranous route, providing a cellular basis for the divergence between intramembranous versus endochondral developmental pathways. However, there is plasticity in this division, as PSCs acquire endochondral bone formation capacity in response to injury. Genetic blockade of the ability of PSCs to give rise to bone-forming osteoblasts results in selective impairments in cortical bone architecture and defects in fracture healing. A cell analogous to mouse PSCs is present in the human periosteum, raising the possibility that PSCs are attractive targets for drug and cellular therapy for skeletal disorders. The identification of PSCs provides evidence that bone contains multiple pools of stem cells, each with distinct physiologic functions.

Modeling functional difference between gyri and sulci within intrinsic connectivity networks.

Recently, the functional roles of the human cortical folding patterns have attracted increasing interest in the neuroimaging community. However, most existing studies have focused on the gyro-sulcal functional relationship on a whole-brain scale but possibly overlooked the localized and subtle functional differences of brain networks. Actually, accumulating evidences suggest that functional brain networks are the basic unit to realize the brain function; thus, the functional relationships between gyri and sulci still need to be further explored within different functional brain networks. Inspired by these evidences, we proposed a novel intrinsic connectivity network (ICN)-guided pooling-trimmed convolutional neural network (I-ptFCN) to revisit the functional difference between gyri and sulci. By testing the proposed model on the task functional magnetic resonance imaging (fMRI) datasets of the Human Connectome Project, we found that the classification accuracy of gyral and sulcal fMRI signals varied significantly for different ICNs, indicating functional heterogeneity of cortical folding patterns in different brain networks. The heterogeneity may be contributed by sulci, as only sulcal signals show heterogeneous frequency features across different ICNs, whereas the frequency features of gyri are homogeneous. These results offer novel insights into the functional difference between gyri and sulci and enlighten the functional roles of cortical folding patterns.

Cerebral cortical regions always connect with each other via the shortest paths.

In human society, the choice of transportation mode between two cities is largely influenced by the distance between the regions. Similarly, when neurons communicate with each other within the cerebral cortex, do they establish their connections based on their physical distance? In this study, we employed a data-driven approach to explore the relationships between fiber length and corresponding geodesic distance between the fiber's two endpoints on brain surface. Diffusion-MRI-derived fiber streamlines were used to represent extra-cortical axonal connections between neurons or cortical regions, while geodesic paths between cortical points were employed to simulate intra-cortical connections. The results demonstrated that the geodesic distance between two cortical regions connected by a fiber streamline was greater than the fiber length most of the time, indicating that cortical regions tend to choose the shortest path for connection; whether it be an intra-cortical or extra-cortical route, especially when intra-cortical routes within cortical regions are longer than potential extrinsic fiber routes, there is an increased probability to establish fiber routes to connect the both regions. These findings were validated in a group of human brains and may provide insights into the underlying mechanisms of neuronal growth, connection, and wiring.

Gyral peaks and patterns in human brains.

Cortical folding patterns are related to brain function, cognition, and behavior. Since the relationship has not been fully explained on a coarse scale, many efforts have been devoted to the identification of finer grained cortical landmarks, such as sulcal pits and gyral peaks, which were found to remain invariant across subjects and ages and the invariance may be related to gene mediated proto-map. However, gyral peaks were only investigated on macaque monkey brains, but not on human brains where the investigation is challenged due to high inter-individual variabilities. To this end, in this work, we successfully identified 96 gyral peaks both on the left and right hemispheres of human brains, respectively. These peaks are spatially consistent across individuals. Higher or sharper peaks are more consistent across subjects. Both structural and functional graph metrics of peaks are significantly different from other cortical regions, and more importantly, these nodal graph metrics are anti-correlated with the spatial consistency metrics within peaks. In addition, the distribution of peaks and various cortical anatomical, structural/functional connective features show hemispheric symmetry. These findings provide new clues to understanding the cortical landmarks, as well as their relationship with brain functions, cognition, behavior in both healthy and aberrant brains.

Interactions between paramyosin and actin greatly improve their thermostability and gel properties.

BACKGROUND: Myofibrillar proteins, the main contributors to the quality of meat products, are the main structural protein component of muscle and have functional properties such as the formation of a 3D protein gel network and water binding. The susceptibility of meat-derived proteins to heat-induced aggregation is the functional constraint that hinders their applications in industry, and so establishing an effective but simple method to improve their thermostability of the proteins is of great importance. RESULTS: In the present study, we describe an easy approach to perform high colloidal thermostability of both paramyosin and actin by mixing them at low ionic strength. The improvement in thermal stability was found to be derived from intermolecular interactions between these two different proteins through non-covalent binding with each other. Consequently, such interactions protected each of them from thermal-induced degradation compared to individual components. Notably, this binary native protein mixture rather than single paramyosin or actin component has the ability to form protein hydrogels with a shear-thinning and reversible sol-gel transformation behavior, which is markedly different from most of reported heat-induced, denatured protein hydrogels. CONCLUSION: The present study not only presents a facile and effective strategy for improvement of the thermal stability and gel properties of a binary paramyosin and actin mixture, but also enhances our understanding of how mutual interactions of protein components affect their physicochemical and functional properties. © 2023 Society of Chemical Industry.

The complex world of kefir: Structural insights and symbiotic relationships.

Kefir milk, known for its high nutritional value and health benefits, is traditionally produced by fermenting milk with kefir grains. These grains are a complex symbiotic community of lactic acid bacteria, acetic acid bacteria, yeasts, and other microorganisms. However, the intricate coexistence mechanisms within these microbial colonies remain a mystery, posing challenges in predicting their biological and functional traits. This uncertainty often leads to variability in kefir milk's quality and safety. This review delves into the unique structural characteristics of kefir grains, particularly their distinctive hollow structure. We propose hypotheses on their formation, which appears to be influenced by the aggregation behaviors of the community members and their alliances. In kefir milk, a systematic colonization process is driven by metabolite release, orchestrating the spatiotemporal rearrangement of ecological niches. We place special emphasis on the dynamic spatiotemporal changes within the kefir microbial community. Spatially, we observe variations in species morphology and distribution across different locations within the grain structure. Temporally, the review highlights the succession patterns of the microbial community, shedding light on their evolving interactions.Furthermore, we explore the ecological mechanisms underpinning the formation of a stable community composition. The interplay of cooperative and competitive species within these microorganisms ensures a dynamic balance, contributing to the community's richness and stability. In kefir community, competitive species foster diversity and stability, whereas cooperative species bolster mutualistic symbiosis. By deepening our understanding of the behaviors of these complex microbial communities, we can pave the way for future advancements in the development and diversification of starter cultures for food fermentation processes.

Genetic Influence on Gyral Peaks.

Genetic mechanisms have been hypothesized to be a major determinant in the formation of cortical folding. Although there is an increasing number of studies examining the heritability of cortical folding, most of them focus on sulcal pits rather than gyral peaks. Gyral peaks, which reflect the highest local foci on gyri and are consistent across individuals, remain unstudied in terms of heritability. To address this knowledge gap, we used high-resolution data from the Human Connectome Project (HCP) to perform classical twin analysis and estimate the heritability of gyral peaks across various brain regions. Our results showed that the heritability of gyral peaks was heterogeneous across different cortical regions, but relatively symmetric between hemispheres. We also found that pits and peaks are different in a variety of anatomic and functional measures. Further, we explored the relationship between the levels of heritability and the formation of cortical folding by utilizing the evolutionary timeline of gyrification. Our findings indicate that the heritability estimates of both gyral peaks and sulcal pits decrease linearly with the evolution timeline of gyrification. This suggests that the cortical folds which formed earlier during gyrification are subject to stronger genetic influences than the later ones. Moreover, the pits and peaks coupled by their time of appearance are also positively correlated in respect of their heritability estimates. These results fill the knowledge gap regarding genetic influences on gyral peaks and significantly advance our understanding of how genetic factors shape the formation of cortical folding. The comparison between peaks and pits suggests that peaks are not a simple morphological mirror of pits but could help complete the understanding of folding patterns.