Transcriptomic landscape of human induced pluripotent stem cell-derived osteogenic differentiation identifies a regulatory role of KLF16.

Osteogenic differentiation is essential for bone development and metabolism, but the underlying gene regulatory networks have not been well investigated. We differentiated mesenchymal stem cells, derived from 20 human induced pluripotent stem cell lines, into preosteoblasts and osteoblasts, and performed systematic RNA-seq analyses of 60 samples for differential gene expression. We noted a highly significant correlation in expression patterns and genomic proximity among transcription factor (TF) and long noncoding RNA (lncRNA) genes. We identified TF-TF regulatory networks, regulatory roles of lncRNAs on their neighboring coding genes for TFs and splicing factors, and differential splicing of TF, lncRNA, and splicing factor genes. TF-TF regulatory and gene co-expression network analyses suggested an inhibitory role of TF KLF16 in osteogenic differentiation. We demonstrate that in vitro overexpression of human KLF16 inhibits osteogenic differentiation and mineralization, and in vivo Klf16+/- mice exhibit increased bone mineral density, trabecular number, and cortical bone area. Thus, our model system highlights the regulatory complexity of osteogenic differentiation and identifies novel osteogenic genes.

Next-generation primate genomics: New genome assemblies unlock new questions.

Nonhuman primates provide unique evolutionary and comparative insight into the human phenotype. Genome assemblies are now available for nearly half of the species in the primate order, expanding our understanding of genetic variation within and between species and making important contributions to evolutionary biology, evolutionary anthropology, and human genetics.

Paired bone marrow and peripheral blood samples demonstrate lack of widespread dissemination of some CH clones.

Clonal hematopoiesis (CH) represents clonal expansion of mutated hematopoietic stem cells detectable in the peripheral blood or bone marrow through next generation sequencing. The current prevailing model posits that CH mutations detected in the peripheral blood mirror bone marrow mutations with clones widely disseminated across hematopoietic compartments. We sought to test the hypothesis that all clones are disseminated throughout hematopoietic tissues by comparing CH in hip vs peripheral blood specimens collected at the time of hip replacement surgery. Here, we show that patients with osteoarthritis have a high prevalence of CH, which involve genes encoding epigenetic modifiers and DNA damage repair pathway proteins. Importantly, we illustrate that CH, including clones with variant allele frequencies >10%, can be confined to specific bone marrow spaces and may be eliminated through surgical excision. Future work will define whether clones with somatic mutations in particular genes or clonal fractions of certain sizes are either more likely to be localized or are slower to disseminate into the peripheral blood and other bony sites.

Evolutionary insights into primate skeletal gene regulation using a comparative cell culture model.

The evolution of complex skeletal traits in primates was likely influenced by both genetic and environmental factors. Because skeletal tissues are notoriously challenging to study using functional genomic approaches, they remain poorly characterized even in humans, let alone across multiple species. The challenges involved in obtaining functional genomic data from the skeleton, combined with the difficulty of obtaining such tissues from nonhuman apes, motivated us to consider an alternative in vitro system with which to comparatively study gene regulation in skeletal cell types. Specifically, we differentiated six human (Homo sapiens) and six chimpanzee (Pan troglodytes) induced pluripotent stem cell lines (iPSCs) into mesenchymal stem cells (MSCs) and subsequently into osteogenic cells (bone cells). We validated differentiation using standard methods and collected single-cell RNA sequencing data from over 100,000 cells across multiple samples and replicates at each stage of differentiation. While most genes that we examined display conserved patterns of expression across species, hundreds of genes are differentially expressed (DE) between humans and chimpanzees within and across stages of osteogenic differentiation. Some of these interspecific DE genes show functional enrichments relevant in skeletal tissue trait development. Moreover, topic modeling indicates that interspecific gene programs become more pronounced as cells mature. Overall, we propose that this in vitro model can be used to identify interspecific regulatory differences that may have contributed to skeletal trait differences between species.

Evolutionary Genetic Signatures of Selection on Bone-Related Variation within Human and Chimpanzee Populations.

Bone strength and the incidence and severity of skeletal disorders vary significantly among human populations, due in part to underlying genetic differentiation. While clinical models predict that this variation is largely deleterious, natural population variation unrelated to disease can go unnoticed, altering our perception of how natural selection has shaped bone morphologies over deep and recent time periods. Here, we conduct the first comparative population-based genetic analysis of the main bone structural protein gene, collagen type I α 1 (COL1A1), in clinical and 1000 Genomes Project datasets in humans, and in natural populations of chimpanzees. Contrary to predictions from clinical studies, we reveal abundant COL1A1 amino acid variation, predicted to have little association with disease in the natural population. We also find signatures of positive selection associated with intron haplotype structure, linkage disequilibrium, and population differentiation in regions of known gene expression regulation in humans and chimpanzees. These results recall how recent and deep evolutionary regimes can be linked, in that bone morphology differences that developed among vertebrates over 450 million years of evolution are the result of positive selection on subtle type I collagen functional variation segregating within populations over time.

An evolutionary perspective of DNA methylation patterns in skeletal tissues using a baboon model of osteoarthritis.

Epigenetic factors, such as DNA methylation, play an influential role in the development of the degenerative joint disease osteoarthritis (OA). These molecular mechanisms have been heavily studied in humans, and although OA affects several other animals in addition to humans, few efforts have taken an evolutionary perspective. This study explores the evolution of OA epigenetics by assessing the relationship between DNA methylation variation and knee OA development in baboons (Papio spp.) and by comparing these findings to human OA epigenetic associations. Genome-wide DNA methylation patterns were identified in bone and cartilage of the right distal femora from 56 pedigreed, adult baboons (28 with and 28 without knee OA) using the Illumina Infinium MethylationEPIC BeadChip. Several significantly differentially methylated positions (DMPs) and regions were found between tissue types. Substantial OA-related differential methylation was also identified in cartilage, but not in bone, suggesting that cartilage epigenetics may be more influential in OA than bone epigenetics. Additionally, some genes containing OA-related DMPs overlap with and display methylation patterns similar to those previously identified in human OA, revealing a mixture of evolutionarily conserved and divergent OA-related methylation patterns in primates. Overall, these findings reinforce the current etiological perspectives of OA and enhance our evolutionary understanding of epigenetic mechanisms associated with OA. This study further establishes baboons as a valuable nonhuman primate model of OA, and continued investigations in baboons will help to disentangle the molecular mechanisms contributing to OA and their evolutionary histories.

Prime time for primate functional genomics.

Functional genomics research is continually improving our understanding of genotype-phenotype relationships in humans, and comparative genomics perspectives can provide additional insight into the evolutionary histories of such relationships. To specifically identify conservation or species-specific divergence in humans, we must look to our closest extant evolutionary relatives. Primate functional genomics research has been steadily advancing and expanding, in spite of several limitations and challenges that this field faces. New technologies and cheaper sequencing provide a unique opportunity to enhance and expand primate comparative studies, and we outline possible paths going forward. The potential human-specific insights that can be gained from primate functional genomics research are substantial, and we propose that now is a prime time to expand such endeavors.

Intraspecific and interspecific investigations of skeletal DNA methylation and femur morphology in primates.

OBJECTIVES: Epigenetic mechanisms influence the development and maintenance of complex phenotypes and may also contribute to the evolution of species-specific phenotypes. With respect to skeletal traits, little is known about the gene regulation underlying these hard tissues or how tissue-specific patterns are associated with bone morphology or vary among species. To begin exploring these topics, this study evaluates one epigenetic mechanism, DNA methylation, in skeletal tissues from five nonhuman primate species which display anatomical and locomotor differences representative of their phylogenetic groups. MATERIALS AND METHODS: First, we test whether intraspecific variation in skeletal DNA methylation is associated with intraspecific variation in femur morphology. Second, we identify interspecific differences in DNA methylation and assess whether these lineage-specific patterns may have contributed to species-specific morphologies. Specifically, we use the Illumina Infinium MethylationEPIC BeadChip to identify DNA methylation patterns in femur trabecular bone from baboons (n = 28), macaques (n = 10), vervets (n = 10), chimpanzees (n = 4), and marmosets (n = 6). RESULTS: Significant differentially methylated positions (DMPs) were associated with a subset of morphological variants, but these likely have small biological effects and may be confounded by other variables associated with morphological variation. Conversely, several species-specific DMPs were identified, and these are found in genes enriched for functions associated with complex skeletal traits. DISCUSSION: Overall, these findings reveal that while intraspecific epigenetic variation is not readily associated with skeletal morphology differences, some interspecific epigenetic differences in skeletal tissues exist and may contribute to evolutionarily distinct phenotypes. This work forms a foundation for future explorations of gene regulation and skeletal trait evolution in primates.

Differential DNA methylation of vocal and facial anatomy genes in modern humans.

Changes in potential regulatory elements are thought to be key drivers of phenotypic divergence. However, identifying changes to regulatory elements that underlie human-specific traits has proven very challenging. Here, we use 63 reconstructed and experimentally measured DNA methylation maps of ancient and present-day humans, as well as of six chimpanzees, to detect differentially methylated regions that likely emerged in modern humans after the split from Neanderthals and Denisovans. We show that genes associated with face and vocal tract anatomy went through particularly extensive methylation changes. Specifically, we identify widespread hypermethylation in a network of face- and voice-associated genes (SOX9, ACAN, COL2A1, NFIX and XYLT1). We propose that these repression patterns appeared after the split from Neanderthals and Denisovans, and that they might have played a key role in shaping the modern human face and vocal tract.

An Introduction to the Callithrix Genus and Overview of Recent Advances in Marmoset Research.

We provide here a current overview of marmoset (Callithrix) evolution, hybridization, species biology, basic/biomedical research, and conservation initiatives. Composed of 2 subgroups, the aurita group (C aurita and C flaviceps) and the jacchus group (C geoffroyi, C jacchus, C kuhlii, and C penicillata), this relatively young primate radiation is endemic to the Brazilian Cerrado, Caatinga, and Atlantic Forest biomes. Significant impacts on Callithrix within these biomes resulting from anthropogenic activity include (1) population declines, particularly for the aurita group; (2) widespread geographic displacement, biological invasions, and range expansions of C jacchus and C penicillata; (3) anthropogenic hybridization; and (4) epizootic Yellow Fever and Zika viral outbreaks. A number of Brazilian legal and conservation initiatives are now in place to protect the threatened aurita group and increase research about them. Due to their small size and rapid life history, marmosets are prized biomedical models. As a result, there are increasingly sophisticated genomic Callithrix resources available and burgeoning marmoset functional, immuno-, and epigenomic research. In both the laboratory and the wild, marmosets have given us insight into cognition, social group dynamics, human disease, and pregnancy. Callithrix jacchus and C penicillata are emerging neotropical primate models for arbovirus disease, including Dengue and Zika. Wild marmoset populations are helping us understand sylvatic transmission and human spillover of Zika and Yellow Fever viruses. All of these factors are positioning marmosets as preeminent models to facilitate understanding of facets of evolution, hybridization, conservation, human disease, and emerging infectious diseases.

Assessment of DNA Methylation Patterns in the Bone and Cartilage of a Nonhuman Primate Model of Osteoarthritis.

OBJECTIVE: Osteoarthritis (OA) affects humans and several other animals. Thus, the mechanisms underlying this disorder, such as specific skeletal tissue DNA methylation patterns, may be evolutionary conserved. However, associations between methylation and OA have not been readily studied in nonhuman animals. Baboons serve as important models of disease and develop OA at rates similar to those in humans. Therefore, this study investigated the associations between methylation and OA in baboons to advance the evolutionary understanding of OA. DESIGN: Trabecular bone and cartilage was collected from the medial condyles of adult female baboon femora, 5 with and 5 without knee OA. The Infinium HumanMethylation450 BeadChip (450K array) was used to identify DNA methylation patterns in these tissues. RESULTS: Approximately 44% of the 450K array probes reliably align to the baboon genome, contain a CpG site of interest, and maintain a wide distribution throughout the genome. Of the 2 filtering methods tested, both identified significantly differentially methylated positions (DMPs) between healthy and OA individuals in cartilage tissues, and some of these patterns overlap with those previously identified in humans. Conversely, no DMPs were found between tissue types or between disease states in bone tissues. CONCLUSIONS: Overall, the 450K array can be used to measure genome-wide DNA methylation in baboon tissues and identify significant associations with complex traits. The results of this study indicate that some DNA methylation patterns associated with OA are evolutionarily conserved, while others are not. This warrants further investigation in a larger and more phylogenetically diverse sample set.

Mycobacterium leprae genomes from naturally infected nonhuman primates.

Leprosy is caused by the bacterial pathogens Mycobacterium leprae and Mycobacterium lepromatosis. Apart from humans, animals such as nine-banded armadillos in the Americas and red squirrels in the British Isles are naturally infected with M. leprae. Natural leprosy has also been reported in certain nonhuman primates, but it is not known whether these occurrences are due to incidental infections by human M. leprae strains or by M. leprae strains specific to nonhuman primates. In this study, complete M. leprae genomes from three naturally infected nonhuman primates (a chimpanzee from Sierra Leone, a sooty mangabey from West Africa, and a cynomolgus macaque from The Philippines) were sequenced. Phylogenetic analyses showed that the cynomolgus macaque M. leprae strain is most closely related to a human M. leprae strain from New Caledonia, whereas the chimpanzee and sooty mangabey M. leprae strains belong to a human M. leprae lineage commonly found in West Africa. Additionally, samples from ring-tailed lemurs from the Bezà Mahafaly Special Reserve, Madagascar, and chimpanzees from Ngogo, Kibale National Park, Uganda, were screened using quantitative PCR assays, to assess the prevalence of M. leprae in wild nonhuman primates. However, these samples did not show evidence of M. leprae infection. Overall, this study adds genomic data for nonhuman primate M. leprae strains to the existing M. leprae literature and finds that this pathogen can be transmitted from humans to nonhuman primates as well as between nonhuman primate species. While the prevalence of natural leprosy in nonhuman primates is likely low, nevertheless, future studies should continue to explore the prevalence of leprosy-causing pathogens in the wild.

The Effects of Histone Deacetylase Inhibitor and Calpain Inhibitor Combination Therapies on Ovarian Cancer Cells.

BACKGROUND: Ovarian cancer is difficult to treat due to absence of selective drugs and tendency of platinum drugs to promote resistance. Combination therapy using epigenetic drugs is predicted to be a beneficial alternative. MATERIALS AND METHODS: This study investigated the effects of combination therapies using two structurally different histone deacetylase (HDAC) inhibitors (HDACi), sodium butyrate and suberanilohydroxamic acid (SAHA), with the calpain inhibitor calpeptin on two characteristically different ovarian cancer cell lines, CAOV-3 and SKOV-3. RESULTS: Suboptimal doses of HDACi and calpeptin produced several effects. Growth inhibition was enhanced and the epigenetically silenced tumor suppressor genes ARHI, p21 and RARß2 were re-expressed. Methylation of specific CpG residues in ARHI were reduced. Cell-cycle progression was inhibited and apoptosis, as well as autophagy, were induced. The phosphorylation of ERK and Akt were differentially effected by these inhibitors. CONCLUSION: The re-expression of tumor suppressors may sensitize ovarian cancer cells, which then undergo apoptosis and autophagy for cell death.

Correction: Validation of qPCR Methods for the Detection of Mycobacterium in New World Animal Reservoirs.

[This corrects the article DOI: 10.1371/journal.pntd.0004198.].

A Comparative Analysis of Genetic and Epigenetic Events of Breast and Ovarian Cancer Related to Tumorigenesis.

Breast cancer persists as the most common cause of cancer death in women worldwide. Ovarian cancer is also a significant source of morbidity and mortality, as the fifth leading cause of cancer death among women. This reflects the continued need for further understanding and innovation in cancer treatment. Though breast and ovarian cancer usually present as distinct clinical entities, the recent explosion of large-scale -omics research has uncovered many overlaps, particularly with respect to genetic and epigenetic alterations. We compared genetic, microenvironmental, stromal, and epigenetic changes common between breast and ovarian cancer cells, as well as the clinical relevance of these changes. Some of the most striking commonalities include genetic alterations of BRCA1 and 2, TP53, RB1, NF1, FAT3, MYC, PTEN, and PIK3CA; down regulation of miRNAs 9, 100, 125a, 125b, and 214; and epigenetic alterations such as H3K27me3, H3K9me2, H3K9me3, H4K20me3, and H3K4me. These parallels suggest shared features of pathogenesis. Furthermore, preliminary evidence suggests a shared epigenetic mechanism of oncogenesis. These similarities, warrant further investigation in order to ultimately inform development of more effective chemotherapeutics, as well as strategies to circumvent drug resistance.

Validation of qPCR Methods for the Detection of Mycobacterium in New World Animal Reservoirs.

Zoonotic pathogens that cause leprosy (Mycobacterium leprae) and tuberculosis (Mycobacterium tuberculosis complex, MTBC) continue to impact modern human populations. Therefore, methods able to survey mycobacterial infection in potential animal hosts are necessary for proper evaluation of human exposure threats. Here we tested for mycobacterial-specific single- and multi-copy loci using qPCR. In a trial study in which armadillos were artificially infected with M. leprae, these techniques were specific and sensitive to pathogen detection, while more traditional ELISAs were only specific. These assays were then employed in a case study to detect M. leprae as well as MTBC in wild marmosets. All marmosets were negative for M. leprae DNA, but 14 were positive for the mycobacterial rpoB gene assay. Targeted capture and sequencing of rpoB and other MTBC genes validated the presence of mycobacterial DNA in these samples and revealed that qPCR is useful for identifying mycobacterial-infected animal hosts.

EMT and tumor metastasis.

EMT and MET comprise the processes by which cells transit between epithelial and mesenchymal states, and they play integral roles in both normal development and cancer metastasis. This article reviews these processes and the molecular pathways that contribute to them. First, we compare embryogenesis and development with cancer metastasis. We then discuss the signaling pathways and the differential expression and down-regulation of receptors in both tumor cells and stromal cells, which play a role in EMT and metastasis. We further delve into the clinical implications of EMT and MET in several types of tumors, and lastly, we discuss the role of epigenetic events that regulate EMT/MET processes. We hypothesize that reversible epigenetic events regulate both EMT and MET, and thus, also regulate the development of different types of metastatic cancers.

Drug resistance in cancer: an overview.

Cancers have the ability to develop resistance to traditional therapies, and the increasing prevalence of these drug resistant cancers necessitates further research and treatment development. This paper outlines the current knowledge of mechanisms that promote or enable drug resistance, such as drug inactivation, drug target alteration, drug efflux, DNA damage repair, cell death inhibition, and the epithelial-mesenchymal transition, as well as how inherent tumor cell heterogeneity plays a role in drug resistance. It also describes the epigenetic modifications that can induce drug resistance and considers how such epigenetic factors may contribute to the development of cancer progenitor cells, which are not killed by conventional cancer therapies. Lastly, this review concludes with a discussion on the best treatment options for existing drug resistant cancers, ways to prevent the formation of drug resistant cancers and cancer progenitor cells, and future directions of study.

Genetic and epigenetic aspects of breast cancer progression and therapy.

Although breast cancer is a heterogeneous disease that is challenging to characterize and treat, the recent explosion of genetic and epigenetic research may help improve these endeavors. In the present review, we use genetic diversity to characterize and classify different types of breast cancer. We also discuss genetic and epigenetic changes that are involved in the development of different breast cancer types and examine how these changes affect prognosis. It appears that while a combination of mutations and copy number changes determine the type of breast cancer, epigenetic alterations may be the primary initiators of cancer development. Understanding these critical biomarkers and molecular changes will advance our ability to effectively treat breast cancer. Next, we examine potential drug therapies directed at epigenetic changes, as such epigenetic drug treatments may prove useful for treating patient-specific tumors, breast cancer progenitor cells, and drug-resistant cells. Lastly, we discuss on mechanisms of carcinogenesis, including a novel hypothesis outlining the role of epigenetics in the development of cancer progenitor cells and metastasis. Overall, breast cancer subtypes may have a similar epigenetic signal that promotes cancer development, and treatment may be most effective if both epigenetic and genetic differences are targeted.

Population genetic structure of traditional populations in the Peruvian Central Andes and implications for South American population history.

Molecular-based characterizations of Andean peoples are traditionally conducted in the service of elucidating continent-level evolutionary processes in South America. Consequently, genetic variation among "western" Andean populations is often represented in relation to variation among "eastern" Amazon and Orinoco River Basin populations. This west-east contrast in patterns of population genetic variation is typically attributed to large-scale phenomena, such as dual founder colonization events or differing long-term microevolutionary histories. However, alternative explanations that consider the nature and causes of population genetic diversity within the Andean region remain underexplored. Here we examine population genetic diversity in the Peruvian Central Andes using data from the mtDNA first hypervariable region and Y-chromosome short tandem repeats among 17 newly sampled populations and 15 published samples. Using this geographically comprehensive data set, we first reassessed the currently accepted pattern of western versus eastern population genetic structure, which our results ultimately reject: mtDNA population diversities were lower, rather than higher, within Andean versus eastern populations, and only highland Y-chromosomes exhibited significantly higher within-population diversities compared with eastern groups. Multiple populations, including several highland samples, exhibited low genetic diversities for both genetic systems. Second, we explored whether the implementation of Inca state and Spanish colonial policies starting at about ad 1400 could have substantially restructured population genetic variation and consequently constitute a primary explanation for the extant pattern of population diversity in the Peruvian Central Andes. Our results suggest that Peruvian Central Andean population structure cannot be parsimoniously explained as the sole outcome of combined Inca and Spanish policies on the region's population demography: highland populations differed from coastal and lowland populations in mtDNA genetic structure only; highland groups also showed strong evidence of female-biased gene flow and/or effective sizes relative to other Peruvian ecozones. Taken together, these findings indicate that population genetic structure in the Peruvian Central Andes is considerably more complex than previously reported and that characterizations of and explanations for genetic variation may be best pursued within more localized regions and defined time periods.