Structure of biofilm-forming functional amyloid PSMα1 from Staphylococcus aureus.

Biofilm-protected pathogenic Staphylococcus aureus causes chronic infections that are difficult to treat. An essential building block of these biofilms are functional amyloid fibrils that assemble from phenol-soluble modulins (PSMs). PSMα1 cross-seeds other PSMs into cross-ß amyloid folds and is therefore a key element in initiating biofilm formation. However, the paucity of high-resolution structures hinders efforts to prevent amyloid assembly and biofilm formation. Here, we present a 3.5 Å resolution density map of the major PSMα1 fibril form revealing a left-handed cross-ß fibril composed of two C2-symmetric U-shaped protofilaments whose subunits are unusually tilted out-of-plane. Monomeric α-helical PSMα1 is extremely cytotoxic to cells, despite the moderate toxicity of the cross-ß fibril. We suggest mechanistic insights into the PSM functional amyloid formation and conformation transformation on the path from monomer-to-fibril formation. Details of PSMα1 assembly and fibril polymorphism suggest how S. aureus utilizes functional amyloids to form biofilms and establish a framework for developing therapeutics against infection and antimicrobial resistance.

Multi-sample non-negative spatial factorization.

It is important to model biological variation when analyzing spatial transcriptomics data from multiple samples. One approach to multi-sample analysis is to spatially align samples, but this is a challenging problem. Here, we provide an alignment-free framework for generalizing a one-sample spatial factorization model to multi-sample data. Using this framework, we develop a method, called multi-sample non-negative spatial factorization (mNSF) that extends the one-sample non-negative spatial factorization (NSF) framework to a multi-sample dataset. Our model allows for a sample-specific model for the spatial correlation structure and extracts a low-dimensional representation of the data. We illustrate the performance of mNSF by simulation studies and real data. mNSF identifies true factors in simulated data, identifies shared anatomical regions across samples in real data and reveals region-specific biological functions. mNSFs performance is similar to alignment based methods when alignment is possible, but extends analysis to situations where spatial alignment is impossible. We expect multi-sample factorization methods to be a powerful class of methods for analyzing spatially resolved transcriptomics data.

Growth deficiency in a mouse model of Kabuki syndrome 2 bears mechanistic similarities to Kabuki syndrome 1.

Growth deficiency is a characteristic feature of both Kabuki syndrome 1 (KS1) and Kabuki syndrome 2 (KS2), Mendelian disorders of the epigenetic machinery with similar phenotypes but distinct genetic etiologies. We previously described skeletal growth deficiency in a mouse model of KS1 and further established that a Kmt2d-/- chondrocyte model of KS1 exhibits precocious differentiation. Here we characterized growth deficiency in a mouse model of KS2, Kdm6atm1d/+. We show that Kdm6atm1d/+ mice have decreased femur and tibia length compared to controls and exhibit abnormalities in cortical and trabecular bone structure. Kdm6atm1d/+ growth plates are also shorter, due to decreases in hypertrophic chondrocyte size and hypertrophic zone height. Given these disturbances in the growth plate, we generated Kdm6a-/- chondrogenic cell lines. Similar to our prior in vitro model of KS1, we found that Kdm6a-/- cells undergo premature, enhanced differentiation towards chondrocytes compared to Kdm6a+/+ controls. RNA-seq showed that Kdm6a-/- cells have a distinct transcriptomic profile that indicates dysregulation of cartilage development. Finally, we performed RNA-seq simultaneously on Kmt2d-/-, Kdm6a-/-, and control lines at Days 7 and 14 of differentiation. This revealed surprising resemblance in gene expression between Kmt2d-/- and Kdm6a-/- at both time points and indicates that the similarity in phenotype between KS1 and KS2 also exists at the transcriptional level.

Towards the diagnosis of osteoporosis - contributions from coincidental diagnostic imaging findings in chiropractors' practice.

BACKGROUND: Osteoporosis is significantly associated with fractures and burdens the health of especially older people. Osteoporotic fractures cause pain, disability, and increased mortality. Early diagnosis of osteoporosis allows earlier initiation of treatment, thereby reducing the risk of osteoporotic fractures. Chiropractors encounter potential osteoporotic patients daily, and perform radiological evaluation of these and other patients, including evaluation of X-rays done for other purposes than osteoporosis. Therefore, chiropractors may identify vertebral fractures, vertebral deformity or osteopenia not otherwise suspected or recorded. METHODS: This study examines procedures available to the chiropractor to describe conventional X-rays with the focus of osteoporosis related findings. We review the indications for radiological examination in chiropractic practice, and in the realm of osteoporosis we describe radiological methods available for examination of conventional radiographs, and the necessity of inter-disciplinary communication. RESULTS: National guidelines are available regarding referral for X-rays in chiropractic practice. Standardized protocols ensure image acquisition of the highest quality in the chiropractors' radiological department. Conventional X-ray examination is not indicated on clinical suspicion of osteoporosis alone, as bone mineral density testing is the diagnostic test. Radiological assessment of all available X-rays of patients above the age of 50 years should include evaluation of the bone quality, and hip and vertebral fracture assessment. The Singh index, Genant Semi-Quantitative tool (GSQ), and Algorithm-Based Qualitative method (ABQ) should be used consistently during interpretation. Referral for additional imaging and evaluation should be prompt and systematic when needed. CONCLUSIONS: This article presents an overview of evidence-based radiological procedures for the purpose of promoting early diagnosis of osteoporosis. We present recommendations to the clinicians where we propose an opportunistic evaluation of X-rays, done for any reason, which include systematic evaluation of bone quality, presence of hip and vertebral fractures, and vertebral deformation of all patients above the age of 50 years. Detailed referral to healthcare professionals for further diagnostic evaluation is performed when needed. Consistent, high-quality radiological procedures in chiropractic practices could feasibly contribute to the timely diagnosis of osteoporosis, ultimately minimizing the impact of osteoporosis-related complications on patients' health.

Removing unwanted variation between samples in Hi-C experiments.

Hi-C data are commonly normalized using single sample processing methods, with focus on comparisons between regions within a given contact map. Here, we aim to compare contact maps across different samples. We demonstrate that unwanted variation, of likely technical origin, is present in Hi-C data with replicates from different individuals, and that properties of this unwanted variation change across the contact map. We present band-wise normalization and batch correction, a method for normalization and batch correction of Hi-C data and show that it substantially improves comparisons across samples, including in a quantitative trait loci analysis as well as differential enrichment across cell types.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice.

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically shows phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption-the root of the pathogenesis-is similar in the different disease-relevant cell types. This is possible in principle, because all these cell types are subject to effects from the same causative gene, which has the same kind of function (e.g., methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2 and find that the chromatin accessibility changes in neurons are mostly distinct from changes in B or T cells. This is not because the neuronal accessibility changes occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that regulatory elements disrupted in B/T cells do show chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators and suggest that blood-derived episignatures, although useful diagnostically, may not be well suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Selective Cell-Surface Expression of Triheteromeric NMDA Receptors.

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

Genotype prediction of 336,463 samples from public expression data.

Tens of thousands of RNA-sequencing experiments comprising hundreds of thousands of individual samples have now been performed. These data represent a broad range of experimental conditions, sequencing technologies, and hypotheses under study. The Recount project has aggregated and uniformly processed hundreds of thousands of publicly available RNA-seq samples. Most of these samples only include RNA expression measurements; genotype data for these same samples would enable a wide range of analyses including variant prioritization, eQTL analysis, and studies of allele specific expression. Here, we developed a statistical model based on the existing reference and alternative read counts from the RNA-seq experiments available through Recount3 to predict genotypes at autosomal biallelic loci in coding regions. We demonstrate the accuracy of our model using large-scale studies that measured both gene expression and genotype genome-wide. We show that our predictive model is highly accurate with 99.5% overall accuracy, 99.6% major allele accuracy, and 90.4% minor allele accuracy. Our model is robust to tissue and study effects, provided the coverage is high enough. We applied this model to genotype all the samples in Recount 3 and provide the largest ready-to-use expression repository containing genotype information. We illustrate that the predicted genotype from RNA-seq data is sufficient to unravel the underlying population structure of samples in Recount3 using Principal Component Analysis.

Dynamic multi-pinhole collimated brain SPECT of Parkinson's disease by [123I]FP-CIT: a feasibility study of fSPECT.

We investigated the feasibility of using a dopamine transporter (DaT) tracer ligand ([123I]FP-CIT) along with novel multi-pinhole brain collimators for dynamic brain single photon emission computed tomography (SPECT) in suspected Parkinson's disease patients. Thirteen patients underwent dynamic tracer acquisitions before standard imaging. Uptake values were corrected for partial volume effects. Specific binding ratio (SBRcalc) was calculated, reflecting binding potential relative to non-displaceable binding (BPND) in the cortex. Additional pharmacokinetic parameters (BPND, R1, k2) were estimated using the simplified reference tissue model, revealing differences between Kahraman low-score (LS) and high-score (HS) groups. Results showed increasing striatal tracer uptake until 100 min post-injection, with consistent values afterward. Uptake and SBRcalc ratios matched visual assessment. LS patients had lower putamen than caudate nucleus tracer uptake, decreased BPND values, while R1 and k2 values were comparable to HS patients. In conclusion, dynamic multi-pinhole SPECT using DaT tracer with the extraction of pharmacokinetic parameters is feasible and could help enable early differentiation of reduced and normal DaT values.

Intrinsically disordered Pseudomonas chaperone FapA slows down the fibrillation of major biofilm-forming functional amyloid FapC.

Functional bacterial amyloids play a crucial role in the formation of biofilms, which mediate chronic infections and contribute to antimicrobial resistance. This study focuses on the FapC amyloid fibrillar protein from Pseudomonas, a major contributor to biofilm formation. We investigate the initial steps of FapC amyloid formation and the impact of the chaperone-like protein FapA on this process. Using solution nuclear magnetic resonance (NMR), we recently showed that both FapC and FapA are intrinsically disordered proteins (IDPs). Here, the secondary structure propensities (SSPs) are compared to alphafold (DeepMind, protein structure prediction tool/algorithm: https://alphafold.ebi.ac.uk/) models. We further demonstrate that the FapA chaperone interacts with FapC and significantly slows down the formation of FapC fibrils. Our NMR titrations reveal ~ 18% of the resonances show FapA-induced chemical shift perturbations (CSPs), which has not been previously observed, the largest being for A82, N201, C237, C240, A241, and G245. These sites may suggest a specific interaction site and/or hotspots of fibrillation inhibition/control interface at the repeat-1 (R1)/loop-2 (L2) and L2/R3 transition areas and at the C-terminus of FapC. Remarkably, ~ 90% of FapA NMR signals exhibit substantial CSPs upon titration with FapC, the largest being for S63, A69, A80, and I92. A temperature-dependent effect of FapA was observed on FapC by thioflavin T (ThT) and NMR experiments. This study provides a detailed understanding of the interaction between the FapA and FapC, shedding light on the regulation and slowing down of amyloid formation, and has important implications for the development of therapeutic strategies targeting biofilms and associated infections.

A mouse model of Weaver syndrome displays overgrowth and excess osteogenesis reversible with KDM6A/6B inhibition.

Weaver syndrome is a Mendelian disorder of the epigenetic machinery (MDEM) caused by germline pathogenic variants in EZH2, which encodes the predominant H3K27 methyltransferase and key enzymatic component of Polycomb repressive complex 2 (PRC2). Weaver syndrome is characterized by striking overgrowth and advanced bone age, intellectual disability, and distinctive facies. We generated a mouse model for the most common Weaver syndrome missense variant, EZH2 p.R684C. Ezh2R684C/R684C mouse embryonic fibroblasts (MEFs) showed global depletion of H3K27me3. Ezh2R684C/+ mice had abnormal bone parameters, indicative of skeletal overgrowth, and Ezh2R684C/+ osteoblasts showed increased osteogenic activity. RNA-Seq comparing osteoblasts differentiated from Ezh2R684C/+, and Ezh2+/+ BM-mesenchymal stem cells (BM-MSCs) indicated collective dysregulation of the BMP pathway and osteoblast differentiation. Inhibition of the opposing H3K27 demethylases KDM6A and KDM6B substantially reversed the excessive osteogenesis in Ezh2R684C/+ cells both at the transcriptional and phenotypic levels. This supports both the ideas that writers and erasers of histone marks exist in a fine balance to maintain epigenome state and that epigenetic modulating agents have therapeutic potential for the treatment of MDEMs.

Precision pharmacological reversal of strain-specific diet-induced metabolic syndrome in mice informed by epigenetic and transcriptional regulation.

Diet-related metabolic syndrome is the largest contributor to adverse health in the United States. However, the study of gene-environment interactions and their epigenomic and transcriptomic integration is complicated by the lack of environmental and genetic control in humans that is possible in mouse models. Here we exposed three mouse strains, C57BL/6J (BL6), A/J, and NOD/ShiLtJ (NOD), to a high-fat, high-carbohydrate diet, leading to varying degrees of metabolic syndrome. We then performed transcriptomic and genome-wide DNA methylation analyses for each strain and found overlapping but also highly divergent changes in gene expression and methylation upstream of the discordant metabolic phenotypes. Strain-specific pathway analysis of dietary effects revealed a dysregulation of cholesterol biosynthesis common to all three strains but distinct regulatory networks driving this dysregulation. This suggests a strategy for strain-specific targeted pharmacologic intervention of these upstream regulators informed by epigenetic and transcriptional regulation. As a pilot study, we administered the drug GW4064 to target one of these genotype-dependent networks, the farnesoid X receptor pathway, and found that GW4064 exerts strain-specific protection against dietary effects in BL6, as predicted by our transcriptomic analysis. Furthermore, GW4064 treatment induced inflammatory-related gene expression changes in NOD, indicating a strain-specific effect in its associated toxicities as well as its therapeutic efficacy. This pilot study demonstrates the potential efficacy of precision therapeutics for genotype-informed dietary metabolic intervention and a mouse platform for guiding this approach.

Pumping the brakes on RNA velocity by understanding and interpreting RNA velocity estimates.

BACKGROUND: RNA velocity analysis of single cells offers the potential to predict temporal dynamics from gene expression. In many systems, RNA velocity has been observed to produce a vector field that qualitatively reflects known features of the system. However, the limitations of RNA velocity estimates are still not well understood. RESULTS: We analyze the impact of different steps in the RNA velocity workflow on direction and speed. We consider both high-dimensional velocity estimates and low-dimensional velocity vector fields mapped onto an embedding. We conclude the transition probability method for mapping velocity estimates onto an embedding is effectively interpolating in the embedding space. Our findings reveal a significant dependence of the RNA velocity workflow on smoothing via the k-nearest-neighbors (k-NN) graph of the observed data. This reliance results in considerable estimation errors for both direction and speed in both high- and low-dimensional settings when the k-NN graph fails to accurately represent the true data structure; this is an unknown feature of real data. RNA velocity performs poorly at estimating speed in both low- and high-dimensional spaces, except in very low noise settings. We introduce a novel quality measure that can identify when RNA velocity should not be used. CONCLUSIONS: Our findings emphasize the importance of choices in the RNA velocity workflow and highlight critical limitations of data analysis. We advise against over-interpreting expression dynamics using RNA velocity, particularly in terms of speed. Finally, we emphasize that the use of RNA velocity in assessing the correctness of a low-dimensional embedding is circular.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice.

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically exhibits phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption - the root of the pathogenesis - is similar in the different disease-relevant cell types. This is possible in principle, since all these cell-types are subject to effects from the same causative gene, that has the same kind of function (e.g. methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2, and find that the chromatin accessibility abnormalities in neurons are mostly distinct from those in B or T cells. This is not because the neuronal abnormalities occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that the regions disrupted in B/T cells do exhibit chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators, and suggest that blood-derived "episignatures" may not be well-suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Novel mouse model of Weaver syndrome displays overgrowth and excess osteogenesis reversible with KDM6A/6B inhibition.

Weaver syndrome is a Mendelian disorder of the epigenetic machinery (MDEM) caused by germline pathogenic variants in EZH2, which encodes the predominant H3K27 methyltransferase and key enzymatic component of Polycomb repressive complex 2 (PRC2). Weaver syndrome is characterized by striking overgrowth and advanced bone age, intellectual disability, and distinctive facies. We generated a mouse model for the most common Weaver syndrome missense variant, EZH2 p.R684C. Ezh2R684C/R684C mouse embryonic fibroblasts (MEFs) showed global depletion of H3K27me3. Ezh2R684C/+ mice had abnormal bone parameters indicative of skeletal overgrowth, and Ezh2R684C/+ osteoblasts showed increased osteogenic activity. RNA-seq comparing osteoblasts differentiated from Ezh2R684C/+ and Ezh2+/+ bone marrow mesenchymal stem cells (BM-MSCs) indicated collective dysregulation of the BMP pathway and osteoblast differentiation. Inhibition of the opposing H3K27 demethylases Kdm6a/6b substantially reversed the excessive osteogenesis in Ezh2R684C/+ cells both at the transcriptional and phenotypic levels. This supports both the ideas that writers and erasers of histone marks exist in a fine balance to maintain epigenome state, and that epigenetic modulating agents have therapeutic potential for the treatment of MDEMs.

nnSVG for the scalable identification of spatially variable genes using nearest-neighbor Gaussian processes.

Feature selection to identify spatially variable genes or other biologically informative genes is a key step during analyses of spatially-resolved transcriptomics data. Here, we propose nnSVG, a scalable approach to identify spatially variable genes based on nearest-neighbor Gaussian processes. Our method (i) identifies genes that vary in expression continuously across the entire tissue or within a priori defined spatial domains, (ii) uses gene-specific estimates of length scale parameters within the Gaussian process models, and (iii) scales linearly with the number of spatial locations. We demonstrate the performance of our method using experimental data from several technological platforms and simulations. A software implementation is available at https://bioconductor.org/packages/nnSVG .

Association of serum and fecal microRNA profiles in cats with gastrointestinal cancer and chronic inflammatory enteropathy.

BACKGROUND: Differentiation of gastrointestinal cancer (GIC) from chronic inflammatory enteropathies (CIE) in cats can be challenging and often requires extensive diagnostic testing. MicroRNAs (miRNAs) have promise as non-invasive biomarkers in serum and feces for diagnosis of GIC. HYPOTHESIS/OBJECTIVES: Cats with GIC will have serum and fecal miRNA profiles that differ significantly from healthy cats and cats with CIE. Identify serum and fecal miRNAs with diagnostic potential for differentiation between cats with GIC and CIE as compared to healthy cats. ANIMALS: Ten healthy cats, 9 cats with CIE, and 10 cats with GIC; all client-owned. METHODS: Cats were recruited for an international multicenter observational prospective case-control study. Serum and feces were screened using small RNA sequencing for miRNAs that differed in abundance between cats with GIC and CIE, and healthy cats. Diagnostic biomarker potential of relevant miRNAs from small RNA sequencing and the literature was confirmed using reverse transcription quantitative real-time PCR (RT-qPCR). RESULTS: Serum miR-223-3p was found to distinguish between cats with GIC and CIE with an area under the curve (AUC) of 0.9 (95% confidence interval [CI], 0.760-1.0), sensitivity of 90% (95% CI, 59.6-99.5%), and specificity of 77.8% (95% CI, 45.3-96.1%). Serum miR-223-3p likewise showed promise in differentiating a subgroup of cats with small cell lymphoma (SCL) from those with CIE. No fecal miRNAs could distinguish between cats with GIC and CIE. CONCLUSION AND CLINICAL IMPORTANCE: Serum miR-223-3p potentially may serve as a noninvasive diagnostic biomarker of GIC in cats, in addition to providing a much needed tool for the differentiation of CIE and SCL.

Special Challenges in PET Imaging of Ectothermic Vertebrates.

The bulk of biomedical positron emission tomography (PET)-scanning experiments are performed on mammals (ie, rodents, pigs, and dogs), and the technique is only infrequently applied to answer research questions in ectothermic vertebrates such as fish, amphibians, and reptiles. Nevertheless, many unique and interesting physiological characteristics in these ectothermic vertebrates could be addressed in detail through PET. The low metabolic rate of ectothermic animals, however, may compromise the validity of physiological and biochemical parameters derived from the images created by PET and other scanning modalities. Here, we review some of the considerations that should be taken into account when PET scanning fish, amphibians, and reptiles. We present specific results from our own experiments, many of which remain previously unpublished, and we draw on examples from the literature. We conclude that knowledge on the natural history and physiology of the species studied and an understanding of the limitations of the PET scanning techniques are necessary to avoid the design of faulty experiments and erroneous conclusions.

3D printed skulls in court - a benefit to stakeholders?

Forensic pathologists may use 3D prints as demonstrative aids when providing expert testimony in court of law, but the effects remain unclear despite many assumed benefits. In this qualitative study, the effects of using a 3D print, demonstrating a blunt force skull fracture, in court were explored by thematic analysis of interviews with judges, prosecutors, defence counsels, and forensic pathologists with the aim of improving the expert testimony. Five semi-structured focus groups and eight one-to-one interviews with a total of 29 stakeholders were transcribed ad verbatim and analysed using thematic analysis. The study found that a highly accurate 3D print of a skull demonstrated autopsy findings in detail and provided a quick overview, but sense of touch was of little benefit as the 3D print had different material characteristics than the human skull. Virtual 3D models were expected to provide all the benefits of 3D prints, be less emotionally confronting, and be logistically feasible. Both 3D prints and virtual 3D models were expected to be less emotionally confronting than autopsy photos. Regardless of fidelity, an expert witness was necessary to translate technical language and explain autopsy findings, and low-fidelity models may be equally suited as demonstrative aids. The court infrequently challenged the expert witnesses' conclusions and, therefore, rarely had a need for viewing autopsy findings in detail, therefore rarely needing a 3D print.