Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial transcriptomics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Advances in mass spectrometry-enabled multiomics at single-cell resolution

We report biological organisms are multifaceted, intricate systems where slight perturbations can result in extensive changes in gene expression, protein abundance and/or activity, and metabolic flux. These changes occur at different timescales, spatially across cells of heterogeneous origins, and within single-cells. Hence, multimodal measurements at the smallest biological scales are necessary to capture dynamic changes in heterogeneous biological systems. Of the analytical techniques used to measure biomolecules, mass spectrometry (MS) has proven to be a powerful option due to its sensitivity, robustness, and flexibility with regard to the breadth of biomolecules that can be analyzed. Recently, many studies have coupled MS to other analytical techniques with the goal of measuring multiple modalities from the same single-cell. It is with these concepts in mind that we focus this review on MS-enabled multiomic measurements at single-cell or near-single- cell resolution.

47 OTHER INSTRUMENTATION↗

Protein N -Glycans in Healthy and Sclerotic Glomeruli in Diabetic Kidney Disease

Diabetes is expected to directly affect renal glycosylation; yet to date, there has not been a comprehensive evaluation of alterations in N-glycan composition in the glomeruli of patients with diabetic kidney disease (DKD). Here, we used untargeted mass spectrometry imaging to identify N-glycan structures in healthy and sclerotic glomeruli in formalin-fixed paraffin-embedded sections from needle biopsies of five patients with DKD and three healthy kidney samples. Regional proteomics was performed on glomeruli from additional biopsies from the same patients to compare the abundances of enzymes involved in glycosylation. Secondary analysis of single-nucleus RNA sequencing (snRNAseq) data were used to inform on transcript levels of glycosylation machinery in different cell types and states. We detected 120 N-glycans, and among them, we identified 12 of these protein post-translated modifications that were significantly increased in glomeruli. All glomeruli-specific N-glycans contained an N-acetyllactosamine epitope. Five N-glycan structures were highly discriminant between sclerotic and healthy glomeruli. Sclerotic glomeruli had an additional set of glycans lacking fucose linked to their core, and they did not show tetra-antennary structures that were common in healthy glomeruli. Orthogonal omics analyses revealed lower protein abundance and lower gene expression involved in synthesizing fucosylated and branched N-glycans in sclerotic podocytes. In snRNAseq and regional proteomics analyses, we observed that genes and/or proteins involved in sialylation and N-acetyllactosamine synthesis were also downregulated in DKD glomeruli, but this alteration remained undetectable by our spatial N-glycomics assay. Integrative spatial glycomics, proteomics, and transcriptomics revealed protein N-glycosylation characteristic of sclerotic glomeruli in DKD.

60 APPLIED LIFE SCIENCES↗

Challenges in spatial metabolomics and proteomics for functional tissue unit and single-cell resolution

While transcriptomics is the most broadly applied technology for global spatial and single cell measurements in healthy and diseased tissues. Transcripts are often used as a proxy for protein and even metabolite measurements, but it has become commonly accepted that extrapolating this kind of information is a poor proxy and not a substitute for direct measurement. Within the last decade advanced developments of mass spectrometry-based assays have made these direct measurements not only possible, but routine. Where mass spectrometry has become an enabling technology, and various methods can now detect hundreds of metabolites and thousands of proteins from samples. Not only can this be performed within bulk measurements, but much effort has been directed into translating these measurements to single cells and tissues at cellular resolution. The information obtained from mass spectrometry is now able to trace metabolic events and decipher feedback loops across anatomical regions, connecting genetic and metabolic networks that define phenotypes. Herein, we will broadly overview developments in the field over the past decade, leading into several case studies which highlight the direct measurement of metabolites, proteins, and proteoforms from thinly sliced tissues. Much of this work is feasible due to multidisciplinary team science, and we offer brief perspective on paths forward and the challenges that persist with adoption and application of spatial omics.

59 BASIC BIOLOGICAL SCIENCES↗

Single-cell and spatial omics in plants: from cellular atlases to regulatory mechanisms

Single-cell RNA sequencing (scRNA-seq) has transformed transcriptomic studies by enabling gene expression profiling at the resolution of individual cells within and across a broad range of tissue types, revealing cellular heterogeneity that is obscured in bulk tissue transcriptomes. Over the past decade, improvements in microfluidics and library preparation have drastically increased throughput, allowing tens of thousands of cells to be assayed in a single experiment. Although initially developed in animal systems, scRNA-seq has rapidly emerged as a powerful and widely adopted approach in plant biology. Beyond transcriptomics, the integration of single-cell data with chromatin accessibility, proteomics, metabolomics, and spatial omics is enabling a system-level understanding of plant gene regulation and cellular organization. Network-based analytical frameworks further support the reconstruction of gene regulatory networks and the interpretation of complex single-cell data. In this review, we summarize the current technological landscape of plant single-cell studies, discuss key experimental and analytical challenges, and review emerging strategies for validating single-cell discoveries. We also discuss future directions in applying single-cell technologies to woody perennials plants and bioenergy-relevant crops, emphasizing their potential to accelerate the discovery of cell type-specific regulatory mechanisms underlying growth, stress resilience, and biomass production.

Li, Miaomiao [ORNL] (ORCID:0000000321326168)↗

Endogenous adenine mediates kidney injury in diabetic models and predicts diabetic kidney disease in patients

Diabetic kidney disease (DKD) can lead to end-stage kidney disease (ESKD) and mortality, however, few mechanistic biomarkers are available for high-risk patients, especially those without macroalbuminuria. Urine from participants with diabetes from Chronic Renal Insufficiency Cohort (CRIC), Singapore Study of Macro-Angiopathy and Reactivity in Type 2 Diabetes (SMART2D), and the Pima Indian Study determined if urine adenine/creatinine ratio (UAdCR) could be a mechanistic biomarker for ESKD. ESKD and mortality were associated with the highest UAdCR tertile in CRIC (HR 1.57, 1.18, 2.10) and SMART2D (HR 1.77, 1.00, 3.12). ESKD was associated with the highest UAdCR tertile in patients without macroalbuminuria in CRIC (HR 2.36, 1.26, 4.39), SMART2D (HR 2.39, 1.08, 5.29), and Pima Indian study (HR 4.57, CI 1.37-13.34). Empagliflozin lowered UAdCR in non-macroalbuminuric participants. Spatial metabolomics localized adenine to kidney pathology and transcriptomics identified ribonucleoprotein biogenesis as a top pathway in proximal tubules of patients without macroalbuminuria, implicating mammalian target of rapamycin (mTOR). Adenine stimulated matrix in tubular cells via mTOR and stimulated mTOR in mouse kidneys. A specific inhibitor of adenine production was found to reduce kidney hypertrophy and kidney injury in diabetic mice. We propose that endogenous adenine may be a causative factor in DKD.

60 APPLIED LIFE SCIENCES↗

A reference tissue atlas for the human kidney

Kidney Precision Medicine Project (KPMP) is building a spatially specified human kidney tissue atlas in health and disease with single-cell resolution. Here, we describe the construction of an integrated reference map of cells, pathways, and genes using unaffected regions of nephrectomy tissues and undiseased human biopsies from 56 adult subjects. We use single-cell/nucleus transcriptomics, subsegmental laser microdissection transcriptomics and proteomics, near-single-cell proteomics, 3D and CODEX imaging, and spatial metabolomics to hierarchically identify genes, pathways, and cells. Integrated data from these different technologies coherently identify cell types/subtypes within different nephron segments and the interstitium. These profiles describe cell-level functional organization of the kidney following its physiological functions and link cell subtypes to genes, proteins, metabolites, and pathways. They further show that messenger RNA levels along the nephron are congruent with the subsegmental physiological activity. This reference atlas provides a framework for the classification of kidney disease when multiple molecular mechanisms underlie convergent clinical phenotypes.

59 BASIC BIOLOGICAL SCIENCES↗

In-vivo neuronal dysfunction by Aβ and tau overlaps with brain-wide inflammatory mechanisms in Alzheimer’s disease

The molecular mechanisms underlying neuronal dysfunction in Alzheimer’s disease (AD) remain uncharacterized. Here, we identify genes, molecular pathways and cellular components associated with whole-brain dysregulation caused by amyloid-beta (Aβ) and tau deposits in the living human brain. We obtained in-vivo resting-state functional MRI (rs-fMRI), Aβ- and tau-PET for 47 cognitively unimpaired and 16 AD participants from the Translational Biomarkers in Aging and Dementia cohort. Adverse neuronal activity impacts by Aβ and tau were quantified with personalized dynamical models by fitting pathology-mediated computational signals to the participant’s real rs-fMRIs. Then, we detected robust brain-wide associations between the spatial profiles of Aβ-tau impacts and gene expression in the neurotypical transcriptome (Allen Human Brain Atlas). Within the obtained distinctive signature of in-vivo neuronal dysfunction, several genes have prominent roles in microglial activation and in interactions with Aβ and tau. Moreover, cellular vulnerability estimations revealed strong association of microglial expression patterns with Aβ and tau’s synergistic impact on neuronal activity (q < 0.001). These results further support the central role of the immune system and neuroinflammatory pathways in AD pathogenesis. Neuronal dysregulation by AD pathologies also associated with neurotypical synaptic and developmental processes. In addition, we identified drug candidates from the vast LINCS library to halt or reduce the observed Aβ-tau effects on neuronal activity. Top-ranked pharmacological interventions target inflammatory, cancer and cardiovascular pathways, including specific medications undergoing clinical evaluation in AD. Our findings, based on the examination of molecular-pathological-functional interactions in humans, may accelerate the process of bringing effective therapies into clinical practice.

60 APPLIED LIFE SCIENCES↗

Capturing an Early Gene Induction Event during Wood Decay by the Brown Rot Fungus Rhodonia placenta

Brown rot fungi dominate wood decomposition in coniferous forests, and their carbohydrate-selective mechanisms are of commercial interest. Brown rot was recently described as a two-step, sequential mechanism orchestrated by fungi using differentially expressed genes (DEGs) and consisting of oxidation via reactive oxygen species (ROS) followed by enzymatic saccharification. There have been indications, however, that the initial oxidation step itself might require induction. To capture this early gene regulation event, here, we integrated fine-scale cryosectioning with whole-transcriptome sequencing to dissect gene expression at the single-hyphal-cell scale (tens of micrometers). This improved the spatial resolution 50-fold, relative to previous work, and we were able to capture the activity of the first 100 μm of hyphal front growth by Rhodonia placenta in aspen wood. This early decay period was dominated by delayed gene expression patterns as the fungus ramped up its mechanism. These delayed DEGs included many genes implicated in ROS pathways (lignocellulose oxidation [LOX]) that were previously and incorrectly assumed to be constitutively expressed. These delayed DEGs, which include those with and without predicted functions, also create a focused subset of target genes for functional genomics. However, this delayed pattern was not universal, with a few genes being upregulated immediately at the hyphal front. Most notably, this included a gene commonly implicated in hydroquinone and iron redox cycling: benzoquinone reductase.

59 BASIC BIOLOGICAL SCIENCES↗

Spatially Resolved Mass Spectrometry at the Single Cell: Recent Innovations in Proteomics and Metabolomics

Biological systems are composed of heterogeneous populations of cells that intercommunicate to form a functional living tissue. Biological function varies greatly across populations of cells, as each single cell has a unique transcriptome, proteome, and metabolome that translates to functional differences within single species and across kingdoms. Over the past decade, substantial advancements in our ability to characterize omic profiles on a single cell level have occurred, including in multiple spectroscopic and mass spectrometry (MS)-based techniques. Of these technologies, spatially resolved mass spectrometry approaches, including mass spectrometry imaging (MSI), have shown the most progress for single cell proteomics and metabolomics. For example, reporter-based methods using heavy metal tags have allowed for targeted MS investigation of the proteome at the subcellular level, and development of technologies such as laser ablation electrospray ionization mass spectrometry (LAESI-MS) now mean that dynamic metabolomics can be performed in situ. In this Perspective, we showcase advancements in single cell spatial metabolomics and proteomics over the past decade and highlight important aspects related to high-throughput screening, data analysis, and more which are vital to the success of achieving proteomic and metabolomic profiling at the single cell scale. Finally, using this broad literature summary, we provide a perspective on how the next decade may unfold in the area of single cell MS-based proteomics and metabolomics.

59 BASIC BIOLOGICAL SCIENCES↗

Spatial metabolomics and lipidomics in kidney disease

Kidney disease is a global health issue that affects over 850 million people, and early detection is key to preventing severe disease and complications. Kidney diseases are associated with complex dysregulation of lipid metabolism. Spatial metabolomics through mass spectrometry imaging (MSI) enables spatial mapping of the lipids in tissue and includes a variety of techniques that can be used to image lipids. In the kidney MSI studies often seek to resolve individual functional tissue units such as glomeruli and proximal tubules. Several different MSI techniques such as matrix-assisted laser desorption/ionization and desorption electrospray ionization have been used to characterize lipids and small molecules in chronic kidney disease, acute kidney injury, genetic kidney disease, and cancer. In this review we provide several examples of how spatial metabolomics data can provide critical information concerning localization of changes in disease states. Additionally, when combined with pathology, measurements, transcriptomics, or proteomics, the metabolomic changes can illuminate underlying mechanisms and provide new clinical insights.

59 BASIC BIOLOGICAL SCIENCES↗

Cell Population–resolved Multiomics Atlas of the Developing Lung

The lung is a vital organ that undergoes extensive morphological and functional changes during postnatal development. To disambiguate how different cell populations contribute to organ development, we performed proteomic and transcriptomic analyses of four sorted cell populations from the lung of human subjects aged 0 to 8 years-old with a focus on early life. The cell populations analyzed included epithelial, endothelial, mesenchymal, and immune cells. Our results revealed distinct molecular signatures for each of the sorted cell populations that enable the description of molecular shifts occurring in these populations during post-natal development. Here, we confirmed that the proteome of the different cell populations was distinct regardless of age and identified functions specific to each population. We identified a series of cell population protein markers, including those located at the cell surface, that show differential expression and distribution on RNA in situ hybridization and immunofluorescence imaging. We validated the spatial distribution of AT1 and endothelial cell surface markers. Temporal analyses of the proteome of each of the four populations revealed processes modulated during postnatal development and disambiguating results obtained on whole tissue proteome. Finally, the proteome was compared to a transcriptomics survey performed on the same lung samples to evaluate processes under post-transcriptional control.

59 BASIC BIOLOGICAL SCIENCES↗

A multi-scale pipeline linking drug transcriptomics with pharmacokinetics predicts in vivo interactions of tuberculosis drugs

Tuberculosis (TB) is the deadliest infectious disease worldwide. The design of new treatments for TB is hindered by the large number of candidate drugs, drug combinations, dosing choices, and complex pharmaco-kinetics/dynamics (PK/PD). Here we study the interplay of these factors in designing combination therapies by linking a machine-learning model, INDIGO-MTB, which predicts in vitro drug interactions using drug transcriptomics, with a multi-scale model of drug PK/PD and pathogen-immune interactions called GranSim. We calculate an in vivo drug interaction score (iDIS) from dynamics of drug diffusion, spatial distribution, and activity within lesions against various pathogen sub-populations. The iDIS of drug regimens evaluated against non-replicating bacteria significantly correlates with efficacy metrics from clinical trials. Our approach identifies mechanisms that can amplify synergistic or mitigate antagonistic drug interactions in vivo by modulating the relative distribution of drugs. Our mechanistic framework enables efficient evaluation of in vivo drug interactions and optimization of combination therapies.

59 BASIC BIOLOGICAL SCIENCES↗

Autonomous adaptive data acquisition for scanning hyperspectral imaging

Non-invasive and label-free spectral microscopy (spectromicroscopy) techniques can provide quantitative biochemical information complementary to genomic sequencing, transcriptomic profiling, and proteomic analyses. However, spectromicroscopy techniques generate high-dimensional data; acquisition of a single spectral image can range from tens of minutes to hours, depending on the desired spatial resolution and the image size. This substantially limits the timescales of observable transient biological processes. To address this challenge and move spectromicroscopy towards efficient real-time spatiochemical imaging, we developed a grid-less autonomous adaptive sampling method. Our method substantially decreases image acquisition time while increasing sampling density in regions of steeper physico-chemical gradients. When implemented with scanning Fourier Transform infrared spectromicroscopy experiments, this grid-less adaptive sampling approach outperformed standard uniform grid sampling in a two-component chemical model system and in a complex biological sample, Caenorhabditis elegans. We quantitatively and qualitatively assess the efficiency of data acquisition using performance metrics and multivariate infrared spectral analysis, respectively.

47 OTHER INSTRUMENTATION↗

Elucidation of Aromatic Catabolic Pathways in White-Rot Fungi

This project aims to investigate the hypothesis that white-rot fungi can simultaneously depolymerize lignin extracellularly and catabolize depolymerization products intracellularly as carbon and energy sources. Evaluating this hypothesis will provide deeper understanding of the role of white-rot fungi in facilitating carbon sequestration in Nature. Additionally, identifying the most promising fungal strains for lignin turnover and catabolism will catalyze future efforts in genetic tool development to enable metabolic engineering in white-rot fungi for lignin bioconversion to bioproducts. Lignin is the second most abundant plant-based biopolymer on Earth and represents up to 40% of the energy density of lignocellulosic biomass. Even though lignin is a massive natural carbon and energy reservoir, only a small group of basidiomycete fungi, namely white-rot fungi (WRF), have evolved the ability to efficiently depolymerize and mineralize lignin to CO2 and H2O. Considerable research efforts have been undertaken to understand how WRF depolymerize lignin but the biochemical reactions that convert lignin into CO2 have been largely neglected. In fact, it is unclear if WRF intracellularly catabolize lignin-derived aromatic compounds to utilize them as a carbon and energy source, or rather if lignin is depolymerized and mineralized extracellularly merely to facilitate access to cellulose and hemicellulose for use as a primary carbon source. To date, we have employed 13C-isotope labeling, systems biology approaches, and in vitro enzyme assays to definitively demonstrate that two WRF, Trametes versicolor and Gelatoporia (Ceriporiopsis) subvermispora, funnel carbon from lignin-derived aromatic compounds into central carbon metabolism via intracellular catabolic pathways [1]. Specifically, 13C-isotopic labeling approaches showed that these WRF utilize poplar-derived aromatic compounds (e.g. 4-hydroxybenzoic acid (4-HBA)) as a carbon source. In silico genome analysis led us to hypothesize a complete catabolic pathway for 4-HBA and identify multiple homologous sequences for enzymes with putative oxidative decarboxylase, hydroxylase, and ring-opening dioxygenase activities, which are among the main biochemical reactions acting on aromatic compounds. Spatial and differential proteomic and metabolomic analyses supported the proposed catabolic pathways and showed alternative catabolic steps in T. versicolor that were not present in G. subvermispora. Based on the in silico, proteomics, and transcriptomics results, we down-selected enzymes for further in vitro characterization, and we have assigned a function to six fungal enzymes (including oxidative decarboxylases, hydroxylases, and ring-opening dioxygenases). Interestingly, even though we selected homologous enzyme pairs from both WRF with similar -omics trends, in a few cases only one of the studied fungi showed activity for the proposed substrate. Based on all the observations from this study [1], we hypothesized that 4-HBA preferentially undergoes oxidative decarboxylation to hydroquinone and subsequent hydroxylation to 1,2,4-benzenetriol in G. subvermispora before ring cleavage, whereas 4-HBA would preferentially undergo hydroxylation to protocatechuate and further oxidative decarboxylation to 1,2,4-benzenetriol in T. versicolor. Examining additional decarboxylases and hydroxylases as well as enzymes from other protein families that can perform the same or similar oxidative reactions, such as cytochromes P450 with aromatic hydroxylation activity, will also be key for elucidating enzyme preferences for specific substrates. Overall, this work forms the foundation of a new research area based on lignin catabolism by WRF, which could be further exploited to convert the undervalued biopolymer lignin into value-added compounds.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Single-Cell Metabolomics with Rapid Determination of Chemical Formulas from Isotopic Fine Structures

Metabolomic measurements can provide functional readouts of cellular states and phenotypes. Here, we present a protocol for single-cell metabolomics that permits direct untargeted detection of a broad number of metabolites under ambient conditions, without the need for sample processing, and with high confidence in the discovery and identification of the molecular formulas for detected metabolites. This protocol describes combining fiber-based laser ablation electrospray ionization (f-LAESI) with a 21 Tesla Fourier transform ion cyclotron resonance mass spectrometer (21T-FTICR-MS) to obtain high confidence molecular formula information about detected metabolites. The f-LAESI source utilizes mid-infrared laser ablation through a sharp optical fiber tip, affording direct ambient analysis of cells without the need for sample processing. Using the 21T-FTICR-MS as a mass analyzer enabled measurement of the isotopic fine structure (IFS) for numerous metabolites simultaneously from single cells, and the IFSs were in turn computationally processed to rapidly determine the corresponding elemental compositions. This metabolomics technique complements other single cell omics measurement methods, helping to resolve complex molecular interactions that take place within cells unattainable from single cell transcriptomic and proteomics methods.

mass spectrometry, FTICR, FTMS, ultrahigh mass res↗

Imaging and spatially resolved mass spectrometry applications in nephrology

The application of spatially resolved mass spectrometry (MS) and MS imaging approaches for studying biomolecular processes in the kidney is rapidly growing. These powerful methods, which enable label-free and multiplexed detection of many molecular classes across omics domains (including metabolites, drugs, proteins and protein post-translational modifications), are beginning to reveal new molecular insights related to kidney health and disease. Further, the complexity of the kidney often necessitates multiple scales of analysis for interrogating biofluids, whole organs, functional tissue units, single cells and subcellular compartments. Various MS methods can generate omics data across these spatial domains and facilitate both basic science and pathological assessment of the kidney. Optimal processes related to sample preparation and handling for different MS applications are rapidly evolving. Emerging technology and methods, improvement of spatial resolution, broader molecular characterization, multimodal and multiomics approaches and the use of machine learning and artificial intelligence approaches promise to make these applications even more valuable in the field of nephology. Overall, spatially resolved MS and MS imaging methods have the potential to fill much of the omics gap in systems biology analysis of the kidney and provide functional outputs that cannot be obtained using genomics and transcriptomic methods.

60 APPLIED LIFE SCIENCES↗

Single molecule-based fliFISH validates radial and heterogeneous gene expression patterns in pancreatic islet ß cells

Single cell RNA sequencing (scRNA-Seq) technologies have greatly enhanced our understanding of islet cell transcriptomes and have revealed the existence of ß cell heterogeneity. However, comparison of scRNA-Seq datasets from different groups have highlighted inconsistencies in gene expression patterns, primarily due to variable detection of lower abundance transcripts. Furthermore, such analyses are unable to uncover the spatial organization of heterogeneous gene expression. Here we used fluctuation localization imaging-based fluorescence in situ hybridization (fliFISH) to quantify transcripts in single cells in mouse pancreatic islet sections. We compared the expression patterns of Insulin 2 (Ins2) with Mafa and Ucn3 – two genes expressed in ß cells as they mature, as well as Rgs4 – a factor with variably reported expression in the islet. This approach accurately quantified transcripts across a wide range of expression levels - from single copies to over hundred copies per cell in one islet. Importantly, fliFISH allowed evaluation of transcript heterogeneity in the spatial context of an intact islet. These studies confirm the existence of a high degree of heterogeneous gene expression levels within the islet and highlight relative and radial expression patterns that likely reflect distinct ß cell maturation states along the radial axis of the islet.

Li, Fangjia↗

Cis-regulatory evolution shapes facial diversity in birds and mammals

Birds and mammals exhibit extraordinary facial diversity, reflecting adaptations to distinct ecological niches and feeding strategies. While core face-building developmental programs are conserved and orchestrated by interactions between ectodermal organizers and the underlying mesenchyme, mechanisms driving facial shape variation remain poorly understood. Here, we integrate single-cell transcriptomic and chromatin accessibility profiling of mouse and chicken developing face to construct a comparative regulatory map. Although both ectodermal and mesenchymal populations display distinct regulatory features in each species, the mesenchyme exhibits markedly greater divergence, pointing to its central role in shaping facial morphology. We further reveal unexpected molecular complexity in the main face-shaping organizer, including a mouse-specific Shh/Wnt5a expression domain. At key morphogen loci (Bmp4, Fgf8, and Wnt5a), conserved and lineage-specific enhancers exhibit spatially restricted activity patterns that mirror divergent signaling domains. These findings demonstrate how cis-regulatory evolution modulates conserved developmental programs to generate morphological novelty, providing a valuable resource for studying vertebrate facial evolution.

Kyomen, Stella [Max Planck Institute for Evolution↗