Sinogram domain angular upsampling of sparse-view micro-CT with dense residual hierarchical transformer and attention-weighted loss.

Reduced angular sampling is a key strategy for increasing scanning efficiency of micron-scale computed tomography (micro-CT). Despite boosting throughput, this strategy introduces noise and extrapolation artifacts due to undersampling. In this work, we present a solution to this issue, by proposing a novel Dense Residual Hierarchical Transformer (DRHT) network to recover high-quality sinograms from 2×, 4× and 8× undersampled scans. DRHT is trained to utilize limited information available from sparsely angular sampled scans and once trained, it can be applied to recover higher-resolution sinograms from shorter scan sessions. Our proposed DRHT model aggregates the benefits of a hierarchical- multi-scale structure along with the combination of local and global feature extraction through dense residual convolutional blocks and non-overlapping window transformer blocks respectively. We also propose a novel noise-aware loss function named KL-L1 to improve sinogram restoration to full resolution. KL-L1, a weighted combination of pixel-level and distribution-level cost functions, leverages inconsistencies in noise distribution and uses learnable spatial weight maps to improve the training of the DRHT model. We present ablation studies and evaluations of our method against other state-of-the-art (SOTA) models over multiple datasets. Our proposed DRHT network achieves an average increase in peak signal to noise ratio (PSNR) of 17.73 dB and a structural similarity index (SSIM) of 0.161, for 8× upsampling, across the three diverse datasets, compared to their respective Bicubic interpolated versions. This novel approach can be utilized to decrease radiation exposure to patients and reduce imaging time for large-scale CT imaging projects.

Quantitative Geometric Modeling of Blood Cells from X-ray Histotomograms of Whole Zebrafish Larvae.

Tissue phenotyping is foundational to understanding and assessing the cellular aspects of disease in organismal context and an important adjunct to molecular studies in the dissection of gene function, chemical effects, and disease. As a first step toward computational tissue phenotyping, we explore the potential of cellular phenotyping from 3-Dimensional (3D), 0.74 µm isotropic voxel resolution, whole zebrafish larval images derived from X-ray histotomography, a form of micro-CT customized for histopathology. As proof of principle towards computational tissue phenotyping of cells, we created a semi-automated mechanism for the segmentation of blood cells in the vascular spaces of zebrafish larvae, followed by modeling and extraction of quantitative geometric parameters. Manually segmented cells were used to train a random forest classifier for blood cells, enabling the use of a generalized cellular segmentation algorithm for the accurate segmentation of blood cells. These models were used to create an automated data segmentation and analysis pipeline to guide the steps in a 3D workflow including blood cell region prediction, cell boundary extraction, and statistical characterization of 3D geometric and cytological features. We were able to distinguish blood cells at two stages in development (4- and 5-days-post-fertilization) and wild-type vs. polA2 huli hutu ( hht ) mutants. The application of geometric modeling across cell types to and across organisms and sample types may comprise a valuable foundation for computational phenotyping that is more open, informative, rapid, objective, and reproducible.

Native American genetic ancestry and pigmentation allele contributions to skin color in a Caribbean population.

Our interest in the genetic basis of skin color variation between populations led us to seek a Native American population with genetically African admixture but low frequency of European light skin alleles. Analysis of 458 genomes from individuals residing in the Kalinago Territory of the Commonwealth of Dominica showed approximately 55% Native American, 32% African, and 12% European genetic ancestry, the highest Native American genetic ancestry among Caribbean populations to date. Skin pigmentation ranged from 20 to 80 melanin units, averaging 46. Three albino individuals were determined to be homozygous for a causative multi-nucleotide polymorphism OCA2NW273KV contained within a haplotype of African origin; its allele frequency was 0.03 and single allele effect size was -8 melanin units. Derived allele frequencies of SLC24A5A111T and SLC45A2L374F were 0.14 and 0.06, with single allele effect sizes of -6 and -4, respectively. Native American genetic ancestry by itself reduced pigmentation by more than 20 melanin units (range 24-29). The responsible hypopigmenting genetic variants remain to be identified, since none of the published polymorphisms predicted in prior literature to affect skin color in Native Americans caused detectable hypopigmentation in the Kalinago.

The variation in skin colour of modern humans is a product of thousands of years of natural selection. All human ancestry can be traced back to African populations, which were dark-skinned to protect them from the intense UV rays of the sun. Over time, humans spread to other parts of the world, and people in the northern latitudes with lower UV developed lighter skin through natural selection. This was likely driven by a need for vitamin D, which requires UV rays for production. Separate genetic mechanisms were involved in the evolution of lighter skin in each of the two main branches of human migration: the European branch (which includes peoples on the Indian subcontinent and Europe) and the East Asian branch (which includes East Asia and the Americas). A variant of the gene SLC24A5 is the primary contributor to lighter skin colour in the European branch, but a corresponding variant driving light skin colour evolution in the East Asian branch remains to be identified. One obstacle to finding such variants is the high prevalence of European ancestry in most people groups, which makes it difficult to separate the influence of European genes from those of other populations. To overcome this issue, Ang et al. studied a population that had a high proportion of Native American and African ancestors, but a relatively small proportion of European ancestors, the Kalinago people. The Kalinago live on the island of Dominica, one of the last Caribbean islands to be colonised by Europeans. Ang et al. were able to collect hundreds of skin pigmentation measurements and DNA samples of the Kalinago, to trace the effect of Native American ancestry on skin colour. Genetic analysis confirmed their oral history records of primarily Native American (55%) ­ one of the highest of any Caribbean population studied to date ­ compared with African (32%) and European (12%) ancestries. Native American ancestry had the highest effect on pigmentation and reduced it by more than 20 melanin units, while the European mutations in the genes SLC24A5 and SLC45A2 and an African gene variant for albinism only contributed 5, 4 and 8 melanin units, respectively. However, none of the so far published gene candidates responsible for skin lightening in Native Americans caused a detectable effect. Therefore, the gene responsible for lighter skin in Native Americans/East Asians has yet to be identified. The work of Ang et al. represents an important step in deciphering the genetic basis of lighter skin colour in Native Americans or East Asians. A better understanding of the genetics of skin pigmentation may help to identify why, for example, East Asians are less susceptible to melanoma than Europeans, despite both having a lighter skin colour. It may also further acceptance of how variations in human skin tones are the result of human migration, random genetic variation, and natural selection for pigmentation in different solar environments.

Prehistoric human migration between Sundaland and South Asia was driven by sea-level rise.

Rapid sea-level rise between the Last Glacial Maximum (LGM) and the mid-Holocene transformed the Southeast Asian coastal landscape, but the impact on human demography remains unclear. Here, we create a paleogeographic map, focusing on sea-level changes during the period spanning the LGM to the present-day and infer the human population history in Southeast and South Asia using 763 high-coverage whole-genome sequencing datasets from 59 ethnic groups. We show that sea-level rise, in particular meltwater pulses 1 A (MWP1A, ~14,500-14,000 years ago) and 1B (MWP1B, ~11,500-11,000 years ago), reduced land area by over 50% since the LGM, resulting in segregation of local human populations. Following periods of rapid sea-level rises, population pressure drove the migration of Malaysian Negritos into South Asia. Integrated paleogeographic and population genomic analysis demonstrates the earliest documented instance of forced human migration driven by sea-level rise.

Promoting validation and cross-phylogenetic integration in model organism research.

Model organism (MO) research provides a basic understanding of biology and disease due to the evolutionary conservation of the molecular and cellular language of life. MOs have been used to identify and understand the function of orthologous genes, proteins, cells and tissues involved in biological processes, to develop and evaluate techniques and methods, and to perform whole-organism-based chemical screens to test drug efficacy and toxicity. However, a growing richness of datasets and the rising power of computation raise an important question: How do we maximize the value of MOs? In-depth discussions in over 50 virtual presentations organized by the National Institutes of Health across more than 10 weeks yielded important suggestions for improving the rigor, validation, reproducibility and translatability of MO research. The effort clarified challenges and opportunities for developing and integrating tools and resources. Maintenance of critical existing infrastructure and the implementation of suggested improvements will play important roles in maintaining productivity and facilitating the validation of animal models of human biology and disease.

Quantitative relationship between cerebrovascular network and neuronal cell types in mice.

The cerebrovasculature and its mural cells must meet brain regional energy demands, but how their spatial relationship with different neuronal cell types varies across the brain remains largely unknown. Here we apply brain-wide mapping methods to comprehensively define the quantitative relationships between the cerebrovasculature, capillary pericytes, and glutamatergic and GABAergic neurons, including neuronal nitric oxide synthase-positive (nNOS+) neurons and their subtypes in adult mice. Our results show high densities of vasculature with high fluid conductance and capillary pericytes in primary motor sensory cortices compared with association cortices that show significant positive and negative correlations with energy-demanding parvalbumin+ and vasomotor nNOS+ neurons, respectively. Thalamo-striatal areas that are connected to primary motor sensory cortices also show high densities of vasculature and pericytes, suggesting dense energy support for motor sensory processing areas. Our cellular-resolution resource offers opportunities to examine spatial relationships between the cerebrovascular network and neuronal cell composition in largely understudied subcortical areas.

Whole-Brain Wiring Diagram of Oxytocin System in Adult Mice.

Oxytocin (Oxt) neurons regulate diverse physiological responses via direct connections with different neural circuits. However, the lack of comprehensive input-output wiring diagrams of Oxt neurons and their quantitative relationship with Oxt receptor (Oxtr) expression presents challenges to understanding circuit-specific Oxt functions. Here, we establish a whole-brain distribution and anatomic connectivity map of Oxt neurons, and their relationship with Oxtr expression using high-resolution 3D mapping methods in adult male and female mice. We use a flatmap to describe Oxt neuronal expression in four hypothalamic domains including under-characterized Oxt neurons in the tuberal nucleus (TU). Oxt neurons in the paraventricular hypothalamus (PVH) broadly project to nine functional circuits that control cognition, brain state, and somatic visceral response. In contrast, Oxt neurons in the supraoptic (SO) and accessory (AN) nuclei have limited central projection to a small subset of the nine circuits. Surprisingly, quantitative comparison between Oxt output and Oxtr expression showed no significant correlation across the whole brain, suggesting abundant indirect Oxt signaling in Oxtr-expressing areas. Unlike output, Oxt neurons in both the PVH and SO receive similar monosynaptic inputs from a subset of the nine circuits mainly in the thalamic, hypothalamic, and cerebral nuclei areas. Our results suggest that PVH-Oxt neurons serve as a central modulator to integrate external and internal information via largely reciprocal connection with the nine circuits while the SO-Oxt neurons act mainly as unidirectional Oxt hormonal output. In summary, our Oxt wiring diagram provides anatomic insights about distinct behavioral functions of Oxt signaling in the brain.SIGNIFICANCE STATEMENT Oxytocin (Oxt) neurons regulate diverse physiological functions from prosocial behavior to pain sensation via central projection in the brain. Thus, understanding detailed anatomic connectivity of Oxt neurons can provide insight on circuit-specific roles of Oxt signaling in regulating different physiological functions. Here, we use high-resolution mapping methods to describe the 3D distribution, monosynaptic input and long-range output of Oxt neurons, and their relationship with Oxt receptor (Oxtr) expression across the entire mouse brain. We found Oxt connections with nine functional circuits controlling cognition, brain state, and somatic visceral response. Furthermore, we identified a quantitatively unmatched Oxt-Oxtr relationship, suggesting broad indirect Oxt signaling. Together, our comprehensive Oxt wiring diagram advances our understanding of circuit-specific roles of Oxt neurons.

A wide-field micro-computed tomography detector: micron resolution at half-centimetre scale.

Ideal three-dimensional imaging of complex samples made up of micron-scale structures extending over mm to cm, such as biological tissues, requires both wide field of view and high resolution. For existing optics and detectors used for micro-CT (computed tomography) imaging, sub-micron pixel resolution can only be achieved for fields of view of <2 mm. This article presents a unique detector system with a 6 mm field-of-view image circle and 0.5 µm pixel size that can be used in micro-CT units utilizing both synchrotron and commercial X-ray sources. A resolution-test pattern with linear microstructures and whole adult Daphnia magna were imaged at beamline 8.3.2 of the Berkeley Advanced Light Source. Volumes of 10000 × 10000 × 7096 isotropic 0.5 µm voxels were reconstructed over a 5.0 mm × 3.5 mm field of view. Measurements in the projection domain confirmed a 0.90 µm measured spatial resolution that is largely Nyquist-limited. This unprecedented combination of field of view and resolution dramatically reduces the need for sectional scans and computational stitching for large samples, ultimately offering the means to elucidate changes in tissue and cellular morphology in the context of larger, whole, intact model organisms and specimens. This system is also anticipated to benefit micro-CT imaging in materials science, microelectronics, agricultural science and biomedical engineering.

Whole-organism 3D quantitative characterization of zebrafish melanin by silver deposition micro-CT.

We previously described X-ray histotomography, a high-resolution, non-destructive form of X-ray microtomography (micro-CT) imaging customized for three-dimensional (3D), digital histology, allowing quantitative, volumetric tissue and organismal phenotyping (Ding et al., 2019). Here, we have combined micro-CT with a novel application of ionic silver staining to characterize melanin distribution in whole zebrafish larvae. The resulting images enabled whole-body, computational analyses of regional melanin content and morphology. Normalized micro-CT reconstructions of silver-stained fish consistently reproduced pigment patterns seen by light microscopy, and further allowed direct quantitative comparisons of melanin content across wild-type and mutant samples, including subtle phenotypes not previously noticed. Silver staining of melanin for micro-CT provides proof-of-principle for whole-body, 3D computational phenomic analysis of a specific cell type at cellular resolution, with potential applications in other model organisms and melanocytic neoplasms. Advances such as this in whole-organism, high-resolution phenotyping provide superior context for studying the phenotypic effects of genetic, disease, and environmental variables.

Chemical Imaging of Retinal Pigment Epithelium in Frozen Sections of Zebrafish Larvae Using ToF-SIMS.

Variants of the SLC24A5 gene, which encodes a putative potassium-dependent sodium-calcium exchanger (NCKX5) that most likely resides in the melanosome or its precursor, affect pigmentation in both humans and zebrafish (Danio rerio). This finding suggests that genetic variations influencing human skin pigmentation alter melanosome biogenesis via ionic changes. Gaining an understanding of how changes in the ionic environment of organelles impact melanosome morphogenesis and pigmentation will require a spatially resolved way to characterize the chemical environment of melanosomes in pigmented tissue such as retinal pigment epithelium (RPE). The imaging mass spectrometry technique most suited for this type of cell and tissue analysis is time-of-flight secondary ion mass spectrometry (ToF-SIMS) because it is able to detect many biochemical species with high sensitivity and with submicron spatial resolution. Here, we describe chemical imaging of the RPE in frozen-hydrated sections of larval zebrafish using cryo-ToF-SIMS. To facilitate the data interpretation, positive and negative polarity ToF-SIMS image data were transformed into a single hyperspectral data set and analyzed using principal component analysis. The combination of a novel protocol and the use of multivariate data analysis allowed us to discover new marker ions that are attributable to leucodopachrome, a metabolite specific to the biosynthesis of eumelanin. The described methodology may be adapted for the investigation of other classes of molecules in frozen tissues from zebrafish and other organisms.

Selective synthetic augmentation with HistoGAN for improved histopathology image classification.

Histopathological analysis is the present gold standard for precancerous lesion diagnosis. The goal of automated histopathological classification from digital images requires supervised training, which requires a large number of expert annotations that can be expensive and time-consuming to collect. Meanwhile, accurate classification of image patches cropped from whole-slide images is essential for standard sliding window based histopathology slide classification methods. To mitigate these issues, we propose a carefully designed conditional GAN model, namely HistoGAN, for synthesizing realistic histopathology image patches conditioned on class labels. We also investigate a novel synthetic augmentation framework that selectively adds new synthetic image patches generated by our proposed HistoGAN, rather than expanding directly the training set with synthetic images. By selecting synthetic images based on the confidence of their assigned labels and their feature similarity to real labeled images, our framework provides quality assurance to synthetic augmentation. Our models are evaluated on two datasets: a cervical histopathology image dataset with limited annotations, and another dataset of lymph node histopathology images with metastatic cancer. Here, we show that leveraging HistoGAN generated images with selective augmentation results in significant and consistent improvements of classification performance (6.7% and 2.8% higher accuracy, respectively) for cervical histopathology and metastatic cancer datasets.

A map of cis-regulatory elements and 3D genome structures in zebrafish.

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species1. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Enhanced and unified anatomical labeling for a common mouse brain atlas.

Anatomical atlases in standard coordinates are necessary for the interpretation and integration of research findings in a common spatial context. However, the two most-used mouse brain atlases, the Franklin-Paxinos (FP) and the common coordinate framework (CCF) from the Allen Institute for Brain Science, have accumulated inconsistencies in anatomical delineations and nomenclature, creating confusion among neuroscientists. To overcome these issues, we adopt here the FP labels into the CCF to merge the labels in the single atlas framework. We use cell type-specific transgenic mice and an MRI atlas to adjust and further segment our labels. Moreover, detailed segmentations are added to the dorsal striatum using cortico-striatal connectivity data. Lastly, we digitize our anatomical labels based on the Allen ontology, create a web-interface for visualization, and provide tools for comprehensive comparisons between the CCF and FP labels. Our open-source labels signify a key step towards a unified mouse brain atlas.

Computational 3D histological phenotyping of whole zebrafish by X-ray histotomography.

Organismal phenotypes frequently involve multiple organ systems. Histology is a powerful way to detect cellular and tissue phenotypes, but is largely descriptive and subjective. To determine how synchrotron-based X-ray micro-tomography (micro-CT) can yield 3-dimensional whole-organism images suitable for quantitative histological phenotyping, we scanned whole zebrafish, a small vertebrate model with diverse tissues, at ~1 micron voxel resolutions. Micro-CT optimized for cellular characterization (histotomography) allows brain nuclei to be computationally segmented and assigned to brain regions, and cell shapes and volumes to be computed for motor neurons and red blood cells. Striking individual phenotypic variation was apparent from color maps of computed densities of brain nuclei. Unlike histology, the histotomography also allows the study of 3-dimensional structures of millimeter scale that cross multiple tissue planes. We expect the computational and visual insights into 3D cell and tissue architecture provided by histotomography to be useful for reference atlases, hypothesis generation, comprehensive organismal screens, and diagnostics.

Rigid Embedding of Fixed and Stained, Whole, Millimeter-Scale Specimens for Section-free 3D Histology by Micro-Computed Tomography.

For over a hundred years, the histological study of tissues has been the gold standard for medical diagnosis because histology allows all cell types in every tissue to be identified and characterized. Our laboratory is actively working to make technological advances in X-ray micro-computed tomography (micro-CT) that will bring the diagnostic power of histology to the study of full tissue volumes at cellular resolution (i.e., an X-ray Histo-tomography modality). Toward this end, we have made targeted improvements to the sample preparation pipeline. One key optimization, and the focus of the present work, is a straightforward method for rigid embedding of fixed and stained millimeter-scale samples. Many of the published methods for sample immobilization and correlative micro-CT imaging rely on placing the samples in paraffin wax, agarose, or liquids such as alcohol. Our approach extends this work with custom procedures and the design of a 3-dimensional printable apparatus to embed the samples in an acrylic resin directly into polyimide tubing, which is relatively transparent to X-rays. Herein, sample preparation procedures are described for the samples from 0.5 to 10 mm in diameter, which would be suitable for whole zebrafish larvae and juveniles, or other animals and tissue samples of similar dimensions. As proof of concept, we have embedded the specimens from Danio, Drosophila, Daphnia, and a mouse embryo; representative images from 3-dimensional scans for three of these samples are shown. Importantly, our methodology leads to multiple benefits including rigid immobilization, long-term preservation of laboriously-created resources, and the ability to re-interrogate samples.

Dietary lipophilic iron alters amyloidogenesis and microglial morphology in Alzheimer's disease knock-in APP mice.

Alzheimer's disease (AD) is a progressive neurodegenerative disorder characterized pathologically by amyloid beta (Aß) deposition, microgliosis, and iron dyshomeostasis. Increased labile iron due to homeostatic dysregulation is believed to facilitate amyloidogenesis. Free iron is incorporated into aggregating amyloid peptides during Aß plaque formation and increases potential for oxidative stress surrounding plaques. The goal of this work was to observe how brain iron levels temporally influence Aß plaque formation, plaque iron concentration, and microgliosis. We fed humanized APPNL-F and APPNL-G-F knock-in mice lipophilic iron compound 3,5,5-trimethylhexanoyl ferrocene (TMHF) and iron deficient diets for twelve months. TMHF elevated brain iron by 22% and iron deficiency decreased brain iron 21% relative to control diet. Increasing brain iron with TMHF accelerated plaque formation, increased Aß staining, and increased senile morphology of amyloid plaques. Increased brain iron was associated with increased plaque-iron loading and microglial iron inclusions. TMHF decreased IBA1+ microglia branch length while increasing roundness indicative of microglial activation. This body of work suggests that increasing mouse brain iron with TMHF potentiates a more human-like Alzheimer's disease phenotype with iron integration into Aß plaques and associated microgliosis.

Comparative analysis of fixation and embedding techniques for optimized histological preparation of zebrafish.

In recognition of the importance of zebrafish as a model organism for studying human disease, we have created zebrafish content for a web-based reference atlas of microanatomy for comparing histology and histopathology between model systems and with humans (http://bio-atlas.psu.edu). Fixation, decalcification, embedding, and sectioning of zebrafish were optimized to maximize section quality. A comparison of protocols involving six fixatives showed that 10% Neutral Buffered Formalin at 21°C for 24h yielded excellent results. Sectioning of juveniles and adults requires bone decalcification; EDTA at 0.35M produced effective decalcification in 21-day-old juveniles through adults (≥~3Months). To improve section plane consistency in sets of larvae, we have developed new array casting molds based on the outside contours of larvae derived from 3D microCT images. Tissue discontinuity in sections, a common barrier to creating quality sections of zebrafish, was minimized by processing and embedding the formalin-fixed zebrafish tissues in plasticized forms of paraffin wax, and by periodic hydration of the block surface in ice water between sets of sections. Optimal H&E (Hematoxylin and Eosin) staining was achieved through refinement of standard protocols. High quality slide scans produced from glass histology slides were digitally processed to maximize image quality, and experimental replicates posted as full slides as part of this publication. Modifications to tissue processing are still needed to eliminate the need for block surface hydration. The further addition of slide collections from other model systems and 3D tools for visualizing tissue architecture would greatly increase the utility of the digital atlas.

Synchrotron microCT imaging of soft tissue in juvenile zebrafish reveals retinotectal projections.

Biomedical research and clinical diagnosis would benefit greatly from full volume determinations of anatomical phenotype. Comprehensive tools for morphological phenotyping are central for the emerging field of phenomics, which requires high-throughput, systematic, accurate, and reproducible data collection from organisms affected by genetic, disease, or environmental variables. Theoretically, complete anatomical phenotyping requires the assessment of every cell type in the whole organism, but this ideal is presently untenable due to the lack of an unbiased 3D imaging method that allows histopathological assessment of any cell type despite optical opacity. Histopathology, the current clinical standard for diagnostic phenotyping, involves the microscopic study of tissue sections to assess qualitative aspects of tissue architecture, disease mechanisms, and physiological state. However, quantitative features of tissue architecture such as cellular composition and cell counting in tissue volumes can only be approximated due to characteristics of tissue sectioning, including incomplete sampling and the constraints of 2D imaging of 5 micron thick tissue slabs. We have used a small, vertebrate organism, the zebrafish, to test the potential of microCT for systematic macroscopic and microscopic morphological phenotyping. While cell resolution is routinely achieved using methods such as light sheet fluorescence microscopy and optical tomography, these methods do not provide the pancellular perspective characteristic of histology, and are constrained by the limited penetration of visible light through pigmented and opaque specimens, as characterizes zebrafish juveniles. Here, we provide an example of neuroanatomy that can be studied by microCT of stained soft tissue at 1.43 micron isotropic voxel resolution. We conclude that synchrotron microCT is a form of 3D imaging that may potentially be adopted towards more reproducible, large-scale, morphological phenotyping of optically opaque tissues. Further development of soft tissue microCT, visualization and quantitative tools will enhance its utility.

Molecular phylogeography of a human autosomal skin color locus under natural selection.

Divergent natural selection caused by differences in solar exposure has resulted in distinctive variations in skin color between human populations. The derived light skin color allele of the SLC24A5 gene, A111T, predominates in populations of Western Eurasian ancestry. To gain insight into when and where this mutation arose, we defined common haplotypes in the genomic region around SLC24A5 across diverse human populations and deduced phylogenetic relationships between them. Virtually all chromosomes carrying the A111T allele share a single 78-kb haplotype that we call C11, indicating that all instances of this mutation in human populations share a common origin. The C11 haplotype was most likely created by a crossover between two haplotypes, followed by the A111T mutation. The two parental precursor haplotypes are found from East Asia to the Americas but are nearly absent in Africa. The distributions of C11 and its parental haplotypes make it most likely that these two last steps occurred between the Middle East and the Indian subcontinent, with the A111T mutation occurring after the split between the ancestors of Europeans and East Asians.

Functional assessment of human coding mutations affecting skin pigmentation using zebrafish.

A major challenge in personalized medicine is the lack of a standard way to define the functional significance of the numerous nonsynonymous, single nucleotide coding variants that are present in each human individual. To begin to address this problem, we have used pigmentation as a model polygenic trait, three common human polymorphisms thought to influence pigmentation, and the zebrafish as a model system. The approach is based on the rescue of embryonic zebrafish mutant phenotypes by "humanized" zebrafish orthologous mRNA. Two hypomorphic polymorphisms, L374F in SLC45A2, and A111T in SLC24A5, have been linked to lighter skin color in Europeans. The phenotypic effect of a second coding polymorphism in SLC45A2, E272K, is unclear. None of these polymorphisms had been tested in the context of a model organism. We have confirmed that zebrafish albino fish are mutant in slc45a2; wild-type slc45a2 mRNA rescued the albino mutant phenotype. Introduction of the L374F polymorphism into albino or the A111T polymorphism into slc24a5 (golden) abolished mRNA rescue of the respective mutant phenotypes, consistent with their known contributions to European skin color. In contrast, the E272K polymorphism had no effect on phenotypic rescue. The experimental conclusion that E272K is unlikely to affect pigmentation is consistent with a lack of correlation between this polymorphism and quantitatively measured skin color in 59 East Asian humans. A survey of mutations causing human oculocutaneous albinism yielded 257 missense mutations, 82% of which are theoretically testable in zebrafish. The developed approach may be extended to other model systems and may potentially contribute to our understanding the functional relationships between DNA sequence variation, human biology, and disease.