Regenerating zebrafish scales express a subset of evolutionary conserved genes involved in human skeletal disease

BACKGROUND: Scales are mineralised exoskeletal structures that are part of the dermal skeleton. Scales have been mostly lost during evolution of terrestrial vertebrates whilst bony fish have retained a mineralised dermal skeleton in the form of fin rays and scales. Each scale is a mineralised collagen plate that is decorated with both matrix-building and resorbing cells. When removed, an ontogenetic scale is quickly replaced following differentiation of the scale pocket-lining cells that regenerate a scale. Processes promoting de novo matrix formation and mineralisation initiated during scale regeneration are poorly understood. Therefore, we performed transcriptomic analysis to determine gene networks and their pathways involved in dermal scale regeneration. RESULTS: We defined the transcriptomic profiles of ontogenetic and regenerating scales of zebrafish and identified 604 differentially expressed genes (DEGs). These were enriched for extracellular matrix, ossification, and cell adhesion pathways, but not in enamel or dentin formation processes indicating that scales are reminiscent to bone. Hypergeometric tests involving monogenetic skeletal disorders showed that DEGs were strongly enriched for human orthologues that are mutated in low bone mass and abnormal bone mineralisation diseases (P< 2× 10(−3)). The DEGs were also enriched for human orthologues associated with polygenetic skeletal traits, including height (P< 6× 10(−4)), and estimated bone mineral density (eBMD, P< 2× 10(−5)). Zebrafish mutants of two human orthologues that were robustly associated with height (COL11A2, P=6× 10(−24)) or eBMD (SPP1, P=6× 10(−20)) showed both exo- and endo- skeletal abnormalities as predicted by our genetic association analyses; col11a2(Y228X/Y228X) mutants showed exoskeletal and endoskeletal features consistent with abnormal growth, whereas spp1(P160X/P160X) mutants predominantly showed mineralisation defects. CONCLUSION: We show that scales have a strong osteogenic expression profile comparable to other elements of the dermal skeleton, enriched in genes that favour collagen matrix growth. Despite the many differences between scale and endoskeletal developmental processes, we also show that zebrafish scales express an evolutionarily conserved sub-population of genes that are relevant to human skeletal disease. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12915-021-01209-8

Intravital FRAP imaging using an E-cadherin-GFP mouse reveals disease- and drug-dependent dynamic regulation of cell-cell junctions in live tissue

E-cadherin-mediated cell-cell junctions play a prominent role in maintaining the epithelial architecture. The disruption or deregulation of these adhesions in cancer can lead to the collapse of tumor epithelia that precedes invasion and subsequent metastasis. Here we generated an E-cadherin-GFP mouse that enables intravital photobleaching and quantification of E-cadherin mobility in live tissue without affecting normal biology. We demonstrate the broad applications of this mouse by examining E-cadherin regulation in multiple tissues, including mammary, brain, liver, and kidney tissue, while specifically monitoring E-cadherin mobility during disease progression in the pancreas. We assess E-cadherin stability in native pancreatic tissue upon genetic manipulation involving Kras and p53 or in response to anti-invasive drug treatment and gain insights into the dynamic remodeling of E-cadherin during in situ cancer progression. FRAP in the E-cadherin-GFP mouse, therefore, promises to be a valuable tool to fundamentally expand our understanding of E-cadherin-mediated events in native microenvironments

Osteocyte transcriptome mapping identifies a molecular landscape controlling skeletal homeostasis and susceptibility to skeletal disease.

Osteocytes are master regulators of the skeleton. We mapped the transcriptome of osteocytes from different skeletal sites, across age and sexes in mice to reveal genes and molecular programs that control this complex cellular-network. We define an osteocyte transcriptome signature of 1239 genes that distinguishes osteocytes from other cells. 77% have no previously known role in the skeleton and are enriched for genes regulating neuronal network formation, suggesting this programme is important in osteocyte communication. We evaluated 19 skeletal parameters in 733 knockout mouse lines and reveal 26 osteocyte transcriptome signature genes that control bone structure and function. We showed osteocyte transcriptome signature genes are enriched for human orthologs that cause monogenic skeletal disorders (P = 2.4 × 10-22) and are associated with the polygenic diseases osteoporosis (P = 1.8 × 10-13) and osteoarthritis (P = 1.6 × 10-7). Thus, we reveal the molecular landscape that regulates osteocyte network formation and function and establish the importance of osteocytes in human skeletal disease

Species Differentiation on a Dynamic Landscape: Shifts in Metapopulation Genetic Structure Using the Chronology of the Hawaiian Archipelago

Species formation during adaptive radiation often occurs in the context of a changing environment. The establishment and arrangement of populations, in space and time, sets up ecological and genetic processes that dictate the rate and pattern of differentiation. Here, we focus on how a dynamic habitat can affect genetic structure, and ultimately, differentiation among populations. We make use of the chronology and geographical history provided by the Hawaiian archipelago to examine the initial stages of population establishment and genetic divergence. We use data from a set of 6 spider lineages that differ in habitat affinities, some preferring low elevation habitats with a longer history of connection, others being more specialized for high elevation and/or wet forest, some with more general habitat affinities. We show that habitat preferences associated with lineages are important in ecological and genetic structuring. Lineages that have more restricted habitat preferences are subject to repeated episodes of isolation and fragmentation as a result of lava flows and vegetation succession. The initial dynamic set up by the landscape translates over time into discrete lineages. Further work is needed to understand how genetic changes interact with a changing set of ecological interactions amongst a shifting mosaic of landscapes to achieve species formation

Ultrafast Evolution and Loss of CRISPRs Following a Host Shift in a Novel Wildlife Pathogen, Mycoplasma gallisepticum

Measureable rates of genome evolution are well documented in human pathogens but are less well understood in bacterial pathogens in the wild, particularly during and after host switches. Mycoplasma gallisepticum (MG) is a pathogenic bacterium that has evolved predominantly in poultry and recently jumped to wild house finches (Carpodacus mexicanus), a common North American songbird. For the first time we characterize the genome and measure rates of genome evolution in House Finch isolates of MG, as well as in poultry outgroups. Using whole-genome sequences of 12 House Finch isolates across a 13-year serial sample and an additional four newly sequenced poultry strains, we estimate a nucleotide diversity in House Finch isolates of only ∼2% of ancestral poultry strains and a nucleotide substitution rate of 0.8−1.2×10−5 per site per year both in poultry and in House Finches, an exceptionally fast rate rivaling some of the highest estimates reported thus far for bacteria. We also found high diversity and complete turnover of CRISPR arrays in poultry MG strains prior to the switch to the House Finch host, but after the invasion of House Finches there is progressive loss of CRISPR repeat diversity, and recruitment of novel CRISPR repeats ceases. Recent (2007) House Finch MG strains retain only ∼50% of the CRISPR repertoire founding (1994–95) strains and have lost the CRISPR–associated genes required for CRISPR function. Our results suggest that genome evolution in bacterial pathogens of wild birds can be extremely rapid and in this case is accompanied by apparent functional loss of CRISPRs

A RhoA-FRET Biosensor Mouse for Intravital Imaging in Normal Tissue Homeostasis and Disease Contexts.

The small GTPase RhoA is involved in a variety of fundamental processes in normal tissue. Spatiotemporal control of RhoA is thought to govern mechanosensing, growth, and motility of cells, while its deregulation is associated with disease development. Here, we describe the generation of a RhoA-fluorescence resonance energy transfer (FRET) biosensor mouse and its utility for monitoring real-time activity of RhoA in a variety of native tissues in vivo. We assess changes in RhoA activity during mechanosensing of osteocytes within the bone and during neutrophil migration. We also demonstrate spatiotemporal order of RhoA activity within crypt cells of the small intestine and during different stages of mammary gestation. Subsequently, we reveal co-option of RhoA activity in both invasive breast and pancreatic cancers, and we assess drug targeting in these disease settings, illustrating the potential for utilizing this mouse to study RhoA activity in vivo in real time

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course

Phylogeographical analysis of the dominant multidrug-resistant H58 clade of Salmonella Typhi identifies inter- and intracontinental transmission events.

The emergence of multidrug-resistant (MDR) typhoid is a major global health threat affecting many countries where the disease is endemic. Here whole-genome sequence analysis of 1,832 Salmonella enterica serovar Typhi (S. Typhi) identifies a single dominant MDR lineage, H58, that has emerged and spread throughout Asia and Africa over the last 30 years. Our analysis identifies numerous transmissions of H58, including multiple transfers from Asia to Africa and an ongoing, unrecognized MDR epidemic within Africa itself. Notably, our analysis indicates that H58 lineages are displacing antibiotic-sensitive isolates, transforming the global population structure of this pathogen. H58 isolates can harbor a complex MDR element residing either on transmissible IncHI1 plasmids or within multiple chromosomal integration sites. We also identify new mutations that define the H58 lineage. This phylogeographical analysis provides a framework to facilitate global management of MDR typhoid and is applicable to similar MDR lineages emerging in other bacterial species

The frequency of osteolytic bone metastasis is determined by conditions of the soil, not the number of seeds; evidence from in vivo models of breast and prostate cancer

Background While both preclinical and clinical studies suggest that the frequency of growing skeletal metastases is elevated in individuals with higher bone turnover, it is unclear whether this is a result of increased numbers of tumour cells arriving in active sites or of higher numbers of tumour cells being induced to divide by the bone micro-environment. Here we have investigated how the differences in bone turnover affect seeding of tumour cells and/or development of overt osteolytic bone metastasis using in vivo models of hormone-independent breast and prostate cancer. Methods Cohorts of 6 (young) and 16 (mature)-week old BALB/c nude mice were culled 1, 7 and 21 days after received intracardiac injection of luciferase expressing human prostate (PC3) or breast cancer (MDA-MB-231) cell lines labelled with a fluorescent cell membrane dye (Vybrant DiD). The presence of growing bone metastases was determined by bioluminescence using an in vivo imaging system (IVIS) and followed by anatomical confirmation of tumour metastatic sites post mortem, while the presence of individual fluorescently labelled tumour cells was evaluated using two-photon microscopy ex vivo. The bone remodelling activities were compared between young and mature naïve mice (both male and female) using micro-CT analysis, ELISA and bone histomorphometry. Results Both prostate and breast cancer cells generated higher numbers of overt skeletal lesions in young mice (~80%) than in mature mice (~20%). Although mature mice presented with fewer overt bone metastases, the number of tumour cells arriving/colonizing in the tibias was comparable between young and mature animals. Young naïve mice had lower bone volume but higher bone formation and resorption activities compared to mature animals. Conclusions Our studies suggest that higher frequencies of growing osteolytic skeletal metastases in these models are linked to increased bone turnover and not to the initial number of tumour cells entering the bone microenvironment