A human fetal lung cell atlas uncovers proximal-distal gradients of differentiation and key regulators of epithelial fates.

We present a multiomic cell atlas of human lung development that combines single-cell RNA and ATAC sequencing, high-throughput spatial transcriptomics, and single-cell imaging. Coupling single-cell methods with spatial analysis has allowed a comprehensive cellular survey of the epithelial, mesenchymal, endothelial, and erythrocyte/leukocyte compartments from 5-22 post-conception weeks. We identify previously uncharacterized cell states in all compartments. These include developmental-specific secretory progenitors and a subtype of neuroendocrine cell related to human small cell lung cancer. Our datasets are available through our web interface (https://lungcellatlas.org). To illustrate its general utility, we use our cell atlas to generate predictions about cell-cell signaling and transcription factor hierarchies which we rigorously test using organoid models.

Human SARS-CoV-2 challenge uncovers local and systemic response dynamics.

The COVID-19 pandemic is an ongoing global health threat, yet our understanding of the dynamics of early cellular responses to this disease remains limited1. Here in our SARS-CoV-2 human challenge study, we used single-cell multi-omics profiling of nasopharyngeal swabs and blood to temporally resolve abortive, transient and sustained infections in seronegative individuals challenged with pre-Alpha SARS-CoV-2. Our analyses revealed rapid changes in cell-type proportions and dozens of highly dynamic cellular response states in epithelial and immune cells associated with specific time points and infection status. We observed that the interferon response in blood preceded the nasopharyngeal response. Moreover, nasopharyngeal immune infiltration occurred early in samples from individuals with only transient infection and later in samples from individuals with sustained infection. High expression of HLA-DQA2 before inoculation was associated with preventing sustained infection. Ciliated cells showed multiple immune responses and were most permissive for viral replication, whereas nasopharyngeal T cells and macrophages were infected non-productively. We resolved 54 T cell states, including acutely activated T cells that clonally expanded while carrying convergent SARS-CoV-2 motifs. Our new computational pipeline Cell2TCR identifies activated antigen-responding T cells based on a gene expression signature and clusters these into clonotype groups and motifs. Overall, our detailed time series data can serve as a Rosetta stone for epithelial and immune cell responses and reveals early dynamic responses associated with protection against infection.

Spatially resolved multiomics of human cardiac niches.

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system1. The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG+ and IgA+ plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs.

Local and systemic responses to SARS-CoV-2 infection in children and adults.

It is not fully understood why COVID-19 is typically milder in children1-3. Here, to examine the differences between children and adults in their response to SARS-CoV-2 infection, we analysed paediatric and adult patients with COVID-19 as well as healthy control individuals (total n = 93) using single-cell multi-omic profiling of matched nasal, tracheal, bronchial and blood samples. In the airways of healthy paediatric individuals, we observed cells that were already in an interferon-activated state, which after SARS-CoV-2 infection was further induced especially in airway immune cells. We postulate that higher paediatric innate interferon responses restrict viral replication and disease progression. The systemic response in children was characterized by increases in naive lymphocytes and a depletion of natural killer cells, whereas, in adults, cytotoxic T cells and interferon-stimulated subpopulations were significantly increased. We provide evidence that dendritic cells initiate interferon signalling in early infection, and identify epithelial cell states associated with COVID-19 and age. Our matching nasal and blood data show a strong interferon response in the airways with the induction of systemic interferon-stimulated populations, which were substantially reduced in paediatric patients. Together, we provide several mechanisms that explain the milder clinical syndrome observed in children.

Single-cell reconstruction of the early maternal-fetal interface in humans.

During early human pregnancy the uterine mucosa transforms into the decidua, into which the fetal placenta implants and where placental trophoblast cells intermingle and communicate with maternal cells. Trophoblast-decidual interactions underlie common diseases of pregnancy, including pre-eclampsia and stillbirth. Here we profile the transcriptomes of about 70,000 single cells from first-trimester placentas with matched maternal blood and decidual cells. The cellular composition of human decidua reveals subsets of perivascular and stromal cells that are located in distinct decidual layers. There are three major subsets of decidual natural killer cells that have distinctive immunomodulatory and chemokine profiles. We develop a repository of ligand-receptor complexes and a statistical tool to predict the cell-type specificity of cell-cell communication via these molecular interactions. Our data identify many regulatory interactions that prevent harmful innate or adaptive immune responses in this environment. Our single-cell atlas of the maternal-fetal interface reveals the cellular organization of the decidua and placenta, and the interactions that are critical for placentation and reproductive success.

Dosimetric impact of contour editing on CT and MRI deep-learning autosegmentation for brain OARs.

PURPOSE: To establish the clinical applicability of deep-learning organ-at-risk autocontouring models (DL-AC) for brain radiotherapy. The dosimetric impact of contour editing, prior to model training, on performance was evaluated for both CT and MRI-based models. The correlation between geometric and dosimetric measures was also investigated to establish whether dosimetric assessment is required for clinical validation. METHOD: CT and MRI-based deep learning autosegmentation models were trained using edited and unedited clinical contours. Autosegmentations were dosimetrically compared to gold standard contours for a test cohort. D1%, D5%, D50%, and maximum dose were used as clinically relevant dosimetric measures. The statistical significance of dosimetric differences between the gold standard and autocontours was established using paired Student's t-tests. Clinically significant cases were identified via dosimetric headroom to the OAR tolerance. Pearson's Correlations were used to investigate the relationship between geometric measures and absolute percentage dose changes for each autosegmentation model. RESULTS: Except for the right orbit, when delineated using MRI models, the dosimetric statistical analysis revealed no superior model in terms of the dosimetric accuracy between the CT DL-AC models or between the MRI DL-AC for any investigated brain OARs. The number of patients where the clinical significance threshold was exceeded was higher for the optic chiasm D1% than other OARs, for all autosegmentation models. A weak correlation was consistently observed between the outcomes of dosimetric and geometric evaluations. CONCLUSIONS: Editing contours before training the DL-AC model had no significant impact on dosimetry. The geometric test metrics were inadequate to estimate the impact of contour inaccuracies on dose. Accordingly, dosimetric analysis is needed to evaluate the clinical applicability of DL-AC models in the brain.

The impact of break duration, time of break onset, and prior shift duration on the amount of sleep between shifts in heavy vehicle drivers.

This study aimed to examine the impact of break duration between consecutive shifts, time of break onset, and prior shift duration on total sleep time (TST) between shifts in heavy vehicle drivers (HVDs), and to assess the interaction between break duration and time of break onset. The sleep (actigraphy and sleep diaries) and work shifts (work diaries) of 27 HVDs were monitored during their usual work schedule for up to 9 weeks. Differences in TST between consecutive shifts and days off were assessed. Linear mixed models (followed by pairwise comparisons) assessed whether break duration, prior shift duration, time of break onset, and the interaction between break duration and break onset were related to TST between shifts. Investigators found TST between consecutive shifts (mean [SD] 6.38 [1.38] h) was significantly less than on days off (mean [SD] 7.63 [1.93] h; p < 0.001). Breaks starting between 12:01 and 8:00 a.m. led to shorter sleep (p < 0.05) compared to breaks starting between 4:01 and 8:00 p.m. Break durations up to 7, 9, and 11 h (Australian and European minimum break durations) resulted in a mean (SD) of 4.76 (1.06), 5.66 (0.77), and 6.41 (1.06) h of sleep, respectively. The impact of shift duration prior to the break and the interaction between break duration and time of break were not significant. HVDs' sleep between workdays is influenced independently by break duration and time of break onset. This naturalistic study provides evidence that current break regulations prevent sufficient sleep duration in this industry. Work regulations should evaluate appropriate break durations and break onset times to allow longer sleep opportunities for HVDs.

Intensity standardization of MRI prior to radiomic feature extraction for artificial intelligence research in glioma-a systematic review.

OBJECTIVES: Radiomics is a promising avenue in non-invasive characterisation of diffuse glioma. Clinical translation is hampered by lack of reproducibility across centres and difficulty in standardising image intensity in MRI datasets. The study aim was to perform a systematic review of different methods of MRI intensity standardisation prior to radiomic feature extraction. METHODS: MEDLINE, EMBASE, and SCOPUS were searched for articles meeting the following eligibility criteria: MRI radiomic studies where one method of intensity normalisation was compared with another or no normalisation, and original research concerning patients diagnosed with diffuse gliomas. Using PRISMA criteria, data were extracted from short-listed studies including number of patients, MRI sequences, validation status, radiomics software, method of segmentation, and intensity standardisation. QUADAS-2 was used for quality appraisal. RESULTS: After duplicate removal, 741 results were returned from database and reference searches and, from these, 12 papers were eligible. Due to a lack of common pre-processing and different analyses, a narrative synthesis was sought. Three different intensity standardisation techniques have been studied: histogram matching (5/12), limiting or rescaling signal intensity (8/12), and deep learning (1/12)-only two papers compared different methods. From these studies, histogram matching produced the more reliable features compared to other methods of altering MRI signal intensity. CONCLUSION: Multiple methods of intensity standardisation have been described in the literature without clear consensus. Further research that directly compares different methods of intensity standardisation on glioma MRI datasets is required. KEY POINTS: • Intensity standardisation is a key pre-processing step in the development of robust radiomic signatures to evaluate diffuse glioma. • A minority of studies compared the impact of two or more methods. • Further research is required to directly compare multiple methods of MRI intensity standardisation on glioma datasets.

Cell Atlas technologies and insights into tissue architecture.

Since Robert Hooke first described the existence of 'cells' in 1665, scientists have sought to identify and further characterise these fundamental units of life. While our understanding of cell location, morphology and function has expanded greatly; our understanding of cell types and states at the molecular level, and how these function within tissue architecture, is still limited. A greater understanding of our cells could revolutionise basic biology and medicine. Atlasing initiatives like the Human Cell Atlas aim to identify all cell types at the molecular level, including their physical locations, and to make this reference data openly available to the scientific community. This is made possible by a recent technology revolution: both in single-cell molecular profiling, particularly single-cell RNA sequencing, and in spatially resolved methods for assessing gene and protein expression. Here, we review available and upcoming atlasing technologies, the biological insights gained to date and the promise of this field for the future.

Inorganic nitrate, hypoxia, and the regulation of cardiac mitochondrial respiration-probing the role of PPARα.

Dietary inorganic nitrate prevents aspects of cardiac mitochondrial dysfunction induced by hypoxia, although the mechanism is not completely understood. In both heart and skeletal muscle, nitrate increases fatty acid oxidation capacity, and in the latter case, this involves up-regulation of peroxisome proliferator-activated receptor (PPAR)α expression. Here, we investigated whether dietary nitrate modifies mitochondrial function in the hypoxic heart in a PPARα-dependent manner. Wild-type (WT) mice and mice without PPARα (Ppara-/-) were given water containing 0.7 mM NaCl (control) or 0.7 mM NaNO3 for 35 d. After 7 d, mice were exposed to normoxia or hypoxia (10% O2) for the remainder of the study. Mitochondrial respiratory function and metabolism were assessed in saponin-permeabilized cardiac muscle fibers. Environmental hypoxia suppressed mass-specific mitochondrial respiration and additionally lowered the proportion of respiration supported by fatty acid oxidation by 18% (P < 0.001). This switch away from fatty acid oxidation was reversed by nitrate treatment in hypoxic WT but not Ppara-/- mice, indicating a PPARα-dependent effect. Hypoxia increased hexokinase activity by 33% in all mice, whereas lactate dehydrogenase activity increased by 71% in hypoxic WT but not Ppara-/- mice. Our findings indicate that PPARα plays a key role in mediating cardiac metabolic remodeling in response to both hypoxia and dietary nitrate supplementation.-Horscroft, J. A., O'Brien, K. A., Clark, A. D., Lindsay, R. T., Steel, A. S., Procter, N. E. K., Devaux, J., Frenneaux, M., Harridge, S. D. R., Murray, A. J. Inorganic nitrate, hypoxia, and the regulation of cardiac mitochondrial respiration-probing the role of PPARα.

C8ORF37 Is Required for Photoreceptor Outer Segment Disc Morphogenesis by Maintaining Outer Segment Membrane Protein Homeostasis.

C8ORF37 is a causative gene for three different clinical forms of incurable retinal degeneration. However, the completely unknown function of C8ORF37 limits our understanding of the pathogenicity of C8ORF37 mutations. Here, we performed a comprehensive phenotypic characterization of a C8orf37 KO mouse line, generated using CRISPR/Cas9 technology. Both C8orf37 KO male and female mice exhibited progressive and simultaneous degeneration of rod and cone photoreceptors but no non-ocular phenotypes. The major ultrastructural feature of C8orf37 KO photoreceptors was massive disorganization of the outer segment (OS) membrane discs starting from the onset of disc morphogenesis during development. At the molecular level, the amounts of multiple OS-specific membrane proteins, including proteins involved in membrane disc organization, were reduced, although these proteins were targeted normally to the OS. Considering the distribution of C8ORF37 throughout the photoreceptor cell body, the normal structure of the KO photoreceptor connecting cilium, and the absence of defects in other ciliary organs of the KO mice, our findings do not support the previous notion that C8ORF37 was a ciliary protein. Because C8ORF37 is absent in the photoreceptor OS, C8ORF37 may participate in the secretory pathway of OS membrane proteins in the photoreceptor cell body and thus maintain the homeostasis of these proteins. This study established a valid animal model for future therapeutic studies of C8ORF37-associated retinal degeneration. This study also shed new light on the role of C8ORF37 in photoreceptors and on the pathogenic mechanism underlying retinal degeneration caused by C8ORF37 mutations.SIGNIFICANCE STATEMENT Inherited retinal degeneration is a group of incurable conditions with poorly understood underlying molecular mechanisms. We investigated C8ORF37, a causative gene for three retinal degenerative conditions: retinitis pigmentosa, cone-rod dystrophy, and Bardet-Biedl syndrome. C8ORF37 encodes a protein with no known functional domains and thus its biological function is unpredictable. We knocked out the C8ORF37 ortholog in mice, which resulted in a retinal phenotype similar to that observed in patients. We further demonstrated that C8ORF37 is required for photoreceptor outer segment disc formation and alignment, a process that is critical for photoreceptor function and survival. This study advances our understanding of the pathogenesis of retinal degeneration and establishes a valuable mouse model for future therapeutic development.

Adjusting to life in a residential aged care facility: Perspectives of people with dementia, family members and facility care staff.

AIMS AND OBJECTIVES: To determine factors that facilitate or impede adjustment to residential aged care (RAC) from the perspectives of residents with dementia, families of residents with dementia and facility staff. BACKGROUND: The transition to a RAC facility can be highly stressful for people with dementia and their families, but we lack an understanding of how people with dementia experience this transition. Knowledge on adjustment to the new environment is essential in order to develop procedures and interventions that better support residents. DESIGN AND METHODS: This study consisted of interviews with 12 residents with dementia who had resided at a RAC facility for six months or less; 14 family members of RAC residents with dementia; and 12 RAC facility staff members. Parallel interview schedules were constructed, with questions on the experience of relocating to RAC for a person with dementia and views on enablers and barriers to successful adjustment. Thematic analysis guided the analysis of data. The study adhered to the consolidated criteria for reporting qualitative research (COREQ) guidelines (see Supplementary File S1). RESULTS: Adjustment to RAC appeared variable, with several residents reporting poor acceptance of their circumstances several months after the relocation. The three groups were largely congruent regarding the importance of support from families and staff, and the development of new relationships with other residents, but not all residents had succeeded in forming friendships. The provision of meaningful activities and opportunities to exert autonomy day-to-day were seen as critical, but staff experienced challenges in providing individualised care due to lack of dedicated time to engage with residents. CONCLUSION AND RELEVANCE TO CLINICAL PRACTICE: There is a need for evaluated interventions to help people with dementia to successfully transition to RAC. Attention should be paid to the way in which care is coordinated within the RAC sector, to enable staff to provide individualised approaches to facilitate adjustment.

Selective activation of microglia facilitates synaptic strength.

Synaptic plasticity is thought to be initiated by neurons only, with the prevailing view assigning glial cells mere specify supportive functions for synaptic transmission and plasticity. We now demonstrate that glial cells can control synaptic strength independent of neuronal activity. Here we show that selective activation of microglia in the rat is sufficient to rapidly facilitate synaptic strength between primary afferent C-fibers and lamina I neurons, the first synaptic relay in the nociceptive pathway. Specifically, the activation of the CX3CR1 receptor by fractalkine induces the release of interleukin-1ß from microglia, which modulates NMDA signaling in postsynaptic neurons, leading to the release of an eicosanoid messenger, which ultimately enhances presynaptic neurotransmitter release. In contrast to the conventional view, this form of plasticity does not require enhanced neuronal activity to trigger the events leading to synaptic facilitation. Augmentation of synaptic strength in nociceptive pathways represents a cellular model of pain amplification. The present data thus suggest that, under chronic pain states, CX3CR1-mediated activation of microglia drives the facilitation of excitatory synaptic transmission in the dorsal horn, which contributes to pain hypersensitivity in chronic pain states.

The role of G-protein receptor 84 in experimental neuropathic pain.

G-protein receptor 84 (GPR84) is an orphan receptor that is induced markedly in monocytes/macrophages and microglia during inflammation, but its pathophysiological function is unknown. Here, we investigate the role of GPR84 in a murine model of traumatic nerve injury. Naive GPR84 knock-out (KO) mice exhibited normal behavioral responses to acute noxious stimuli, but subsequent to partial sciatic nerve ligation (PNL), KOs did not develop mechanical or thermal hypersensitivity, in contrast to wild-type (WT) littermates. Nerve injury increased ionized calcium binding adapter molecule 1 (Iba1) and phosphorylated p38 MAPK immunoreactivity in the dorsal horn and Iba1 and cluster of differentiation 45 expression in the sciatic nerve, with no difference between genotypes. PCR array analysis revealed that Gpr84 expression was upregulated in the spinal cord and sciatic nerve of WT mice. In addition, the expression of arginase-1, a marker for anti-inflammatory macrophages, was upregulated in KO sciatic nerve. Based on this evidence, we investigated whether peripheral macrophages behave differently in the absence of GPR84. We found that lipopolysaccharide-stimulated KO macrophages exhibited attenuated expression of several proinflammatory mediators, including IL-1ß, IL-6, and TNF-α. Forskolin-stimulated KO macrophages also showed greater cAMP induction, a second messenger associated with immunosuppression. In summary, our results demonstrate that GPR84 is a proinflammatory receptor that contributes to nociceptive signaling via the modulation of macrophages, whereas in its absence the response of these cells to an inflammatory insult is impaired.

Neuron-immune mechanisms contribute to pain in early stages of arthritis.

BACKGROUND: Rheumatoid arthritis (RA) patients frequently show weak correlations between the magnitude of pain and inflammation suggesting that mechanisms other than overt peripheral inflammation contribute to pain in RA. We assessed changes in microglial reactivity and spinal excitability and their contribution to pain-like behaviour in the early stages of collagen-induced arthritis (CIA) model. METHODS: Mechanically evoked hypersensitivity, spinal nociceptive withdrawal reflexes (NWRs) and hind paw swelling were evaluated in female Lewis rats before and until 13 days following collagen immunization. In the spinal dorsal horn, microgliosis was assayed using immunohistochemistry (Iba-1/p-p38) and cyto(chemo)kine levels in the cerebrospinal fluid (CSF). Intrathecal administration of microglia-targeting drugs A-438079 (P2X7 antagonist) and LHVS (cathepsin S inhibitor) were examined upon hypersensitivity, NWRs, microgliosis and cyto(chemo)kine levels in the early phase of CIA. RESULTS: The early phase of CIA was associated with mechanical allodynia and exaggerated mechanically evoked spinal NWRs, evident before hind paw swelling, and exacerbated with the development of swelling. Concomitant with the development of hypersensitivity was the presence of reactive spinal microgliosis and an increase of IL-1ß levels in CSF (just detectable in plasma). Prolonged intrathecal administration of microglial inhibitors attenuated the development of mechanical allodynia, reduced microgliosis and attenuated IL-1ß increments. Acute spinal application of either microglial inhibitor significantly diminished the sensitization of the spinal NWRs. CONCLUSIONS: Mechanical hypersensitivity in the early phase of CIA is associated with central sensitization that is dependent upon microglial-mediated release of IL-1ß in the spinal cord. Blockade of these spinal events may provide pain relief in RA patients.

The role of glia in the spinal cord in neuropathic and inflammatory pain.

Chronic pain, both inflammatory and neuropathic, is a debilitating condition in which the pain experience persists after the painful stimulus has resolved. The efficacy of current treatment strategies using opioids, NSAIDS and anticonvulsants is limited by the extensive side effects observed in patients, underlining the necessity for novel therapeutic targets. Preclinical models of chronic pain have recently provided evidence for a critical role played by glial cells in the mechanisms underlying the chronicity of pain, both at the site of damage in the periphery and in the dorsal horn of the spinal cord. Here microglia and astrocytes respond to the increased input from the periphery and change morphology, increase in number and release pro-nociceptive mediators such as ATP, cytokines and chemokines. These gliotransmitters can sensitise neurons by activation of their cognate receptors thereby contributing to central sensitization which is fundamental for the generation of allodynia, hyperalgesia and spontaneous pain.

Lhx2 balances progenitor maintenance with neurogenic output and promotes competence state progression in the developing retina.

The LIM-Homeodomain transcription factor Lhx2 is an essential organizer of early eye development and is subsequently expressed in retinal progenitor cells (RPCs). To determine its requirement in RPCs, we performed a temporal series of conditional inactivations in mice with the early RPC driver Pax6 α-Cre and the tamoxifen-inducible Hes1(CreERT2) driver. Deletion of Lhx2 caused a significant reduction of the progenitor population and a corresponding increase in neurogenesis. Precursor fate choice correlated with the time of inactivation; early and late inactivation led to the overproduction of retinal ganglion cells (RGCs) and rod photoreceptors, respectively. In each case, however, the overproduction was selective, occurring at the expense of other cell types and indicating a role for Lhx2 in generating cell type diversity. RPCs that persisted in the absence of Lhx2 continued to generate RGC precursors beyond their normal production window, suggesting that Lhx2 facilitates a transition in competence state. These results identify Lhx2 as a key regulator of RPC properties that contribute to the ordered production of multiple cell types during retinal tissue formation.

Coastal groynes reduce beach litter accumulation along the East coast of England.

Anthropogenic marine litter (AML) is a global environmental concern. One of the most conspicuous effects of AML is beach litter accumulation, the distribution of which is typically heterogenous. Little information is available on the potential effects of coastal topographic features on litter dispersal. We analysed the abundance, composition, and sources of beach litter on the East coast of England in relation to the presence of coastal groyne structures. Six beaches were surveyed in autumn and winter 2021 using the OSPAR methodology for monitoring beach litter. Litter abundance was lower on beaches with groynes present, which could infer that groynes deflect or bury AML. The presence of groynes had no significant effect on the composition/sources of beach litter. Single-use plastic packaging, fishing waste, and sewage-related debris were the largest contributors of beach litter in this region. Our findings indicate that man-made topographic features may affect marine litter dispersal and coastal accumulation.

Fetal face shape analysis from prenatal 3D ultrasound images.

3D ultrasound imaging of fetal faces has been predominantly confined to qualitative assessment. Many genetic conditions evade diagnosis and identification could assist with parental counselling, pregnancy management and neonatal care planning. We describe a methodology to build a shape model of the third trimester fetal face from 3D ultrasound and show how it can objectively describe morphological features and gestational-age related changes of normal fetal faces. 135 fetal face 3D ultrasound volumes (117 appropriately grown, 18 growth-restricted) of 24-34 weeks gestation were included. A 3D surface model of each face was obtained using a semi-automatic segmentation workflow. Size normalisation and rescaling was performed using a growth model giving the average size at every gestation. The model demonstrated a similar growth rate to standard head circumference reference charts. A landmark-free morphometry model was estimated to characterize shape differences using non-linear deformations of an idealized template face. Advancing gestation is associated with widening/fullness of the cheeks, contraction of the chin and deepening of the eyes. Fetal growth restriction is associated with a smaller average facial size but no morphological differences. This model may eventually be used as a reference to assist in the prenatal diagnosis of congenital anomalies with characteristic facial dysmorphisms.