Tissue-specific epigenetics in gene neighborhoods: myogenic transcription factor genes.

Myogenic regulatory factor (MRF) genes, MYOD1, MYOG, MYF6 and MYF5, are critical for the skeletal muscle lineage. Here, we used various epigenome profiles from human myoblasts (Mb), myotubes (Mt), muscle and diverse non-muscle samples to elucidate the involvement of multigene neighborhoods in the regulation of MRF genes. We found more far-distal enhancer chromatin associated with MRF genes in Mb and Mt than previously reported from studies in mice. For the MYF5/MYF6 gene-pair, regions of Mb-associated enhancer chromatin were located throughout the adjacent 236-kb PTPRQ gene even though Mb expressed negligible amounts of PTPRQ mRNA. Some enhancer chromatin regions inside PTPRQ in Mb were also seen in PTPRQ mRNA-expressing non-myogenic cells. This suggests dual-purpose PTPRQ enhancers that upregulate expression of PTPRQ in non-myogenic cells and MYF5/MYF6 in myogenic cells. In contrast, the myogenic enhancer chromatin regions distal to MYOD1 were intergenic and up to 19 kb long. Two of them contain small, known MYOD1 enhancers, and one displayed an unusually high level of 5-hydroxymethylcytosine in a quantitative DNA hydroxymethylation assay. Unexpectedly, three regions of MYOD1-distal enhancer chromatin in Mb and Mt overlapped enhancer chromatin in umbilical vein endothelial cells, which might upregulate a distant gene (PIK3C2A). Lastly, genes surrounding MYOG were preferentially transcribed in Mt, like MYOG itself, and exhibited nearby myogenic enhancer chromatin. These neighboring chromatin regions may be enhancers acting in concert to regulate myogenic expression of multiple adjacent genes. Our findings reveal the very different and complex organization of gene neighborhoods containing closely related transcription factor genes.

DNA Hypomethylation in Intragenic and Intergenic Enhancer Chromatin of Muscle-Specific Genes Usually Correlates with their Expression.

Tissue-specific enhancers are critical for gene regulation. In this study, we help elucidate the contribution of muscle-associated differential DNA methylation to the enhancer activity of highly muscle-specific genes. By bioinformatic analysis of 44 muscle-associated genes, we show that preferential gene expression in skeletal muscle (SkM) correlates with SkM-specific intragenic and intergenic enhancer chromatin and overlapping foci of DNA hypomethylation. Some genes, e.g., CASQ1 and FBXO32, displayed broad regions of both SkM- and heart-specific enhancer chromatin but exhibited focal SkM-specific DNA hypomethylation. Half of the genes had SkM-specific super-enhancers. In contrast to simple enhancer/gene-expression correlations, a super-enhancer was associated with the myogenic MYOD1 gene in both SkM and myoblasts even though SkM has < 1 percent as much MYOD1 expression. Local chromatin differences in this super-enhancer probably contribute to the SkM/myoblast differential expression. Transfection assays confirmed the tissue-specificity of the 0.3-kb core enhancer within MYOD1's super-enhancer and demonstrated its repression by methylation of its three CG dinucleotides. Our study suggests that DNA hypomethylation increases enhancer tissue-specificity and that SkM super-enhancers sometimes are poised for physiologically important, rapid up-regulation.

Differences in gastric carcinoma microenvironment stratify according to EBV infection intensity: implications for possible immune adjuvant therapy.

Epstein-Barr virus (EBV) is associated with roughly 10% of gastric carcinomas worldwide (EBVaGC). Although previous investigations provide a strong link between EBV and gastric carcinomas, these studies were performed using selected EBV gene probes. Using a cohort of gastric carcinoma RNA-seq data sets from The Cancer Genome Atlas (TCGA), we performed a quantitative and global assessment of EBV gene expression in gastric carcinomas and assessed EBV associated cellular pathway alterations. EBV transcripts were detected in 17% of samples but these samples varied significantly in EBV coverage depth. In four samples with the highest EBV coverage (hiEBVaGC - high EBV associated gastric carcinoma), transcripts from the BamHI A region comprised the majority of EBV reads. Expression of LMP2, and to a lesser extent, LMP1 were also observed as was evidence of abortive lytic replication. Analysis of cellular gene expression indicated significant immune cell infiltration and a predominant IFNG response in samples expressing high levels of EBV transcripts relative to samples expressing low or no EBV transcripts. Despite the apparent immune cell infiltration, high levels of the cytotoxic T-cell (CTL) and natural killer (NK) cell inhibitor, IDO1, was observed in the hiEBVaGCs samples suggesting an active tolerance inducing pathway in this subgroup. These results were confirmed in a separate cohort of 21 Vietnamese gastric carcinoma samples using qRT-PCR and on tissue samples using in situ hybridization and immunohistochemistry. Lastly, a panel of tumor suppressors and candidate oncogenes were expressed at lower levels in hiEBVaGC versus EBV-low and EBV-negative gastric cancers suggesting the direct regulation of tumor pathways by EBV.

Modeling, simulation and analysis of methylation profiles from reduced representation bisulfite sequencing experiments.

The ENCODE project has funded the generation of a diverse collection of methylation profiles using reduced representation bisulfite sequencing (RRBS) technology, enabling the analysis of epigenetic variation on a genomic scale at single-site resolution. A standard application of RRBS experiments is in the location of differentially methylated regions (DMRs) between two sets of samples. Despite numerous publications reporting DMRs identified from RRBS datasets, there have been no formal analyses of the effects of experimental and biological factors on the performance of existing or newly developed analytical methods. These factors include variable read coverage, differing group sample sizes across genomic regions, uneven spacing between CpG dinucleotide sites, and correlation in methylation levels among sites in close proximity. To better understand the interplay among technical and biological variables in the analysis of RRBS methylation profiles, we have developed an algorithm for the generation of experimentally realistic RRBS datasets. Applying insights derived from our simulation studies, we present a novel procedure that can identify DMRs spanning as few as three CpG sites with both high sensitivity and specificity. Using RRBS data from muscle vs. non-muscle cell cultures as an example, we demonstrate that our method reveals many more DMRs that are likely to be of biological significance than previous methods.

Epigenetics of Genes Preferentially Expressed in Dissimilar Cell Populations: Myoblasts and Cerebellum.

While studying myoblast methylomes and transcriptomes, we found that CDH15 had a remarkable preference for expression in both myoblasts and cerebellum. To understand how widespread such a relationship was and its epigenetic and biological correlates, we systematically looked for genes with similar transcription profiles and analyzed their DNA methylation and chromatin state and accessibility profiles in many different cell populations. Twenty genes were expressed preferentially in myoblasts and cerebellum (Myob/Cbl genes). Some shared DNA hypo- or hypermethylated regions in myoblasts and cerebellum. Particularly striking was ZNF556, whose promoter is hypomethylated in expressing cells but highly methylated in the many cell populations that do not express the gene. In reporter gene assays, we demonstrated that its promoter's activity is methylation sensitive. The atypical epigenetics of ZNF556 may have originated from its promoter's hypomethylation and selective activation in sperm progenitors and oocytes. Five of the Myob/Cbl genes (KCNJ12, ST8SIA5, ZIC1, VAX2, and EN2) have much higher RNA levels in cerebellum than in myoblasts and displayed myoblast-specific hypermethylation upstream and/or downstream of their promoters that may downmodulate expression. Differential DNA methylation was associated with alternative promoter usage for Myob/Cbl genes MCF2L, DOK7, CNPY1, and ANK1. Myob/Cbl genes PAX3, LBX1, ZNF556, ZIC1, EN2, and VAX2 encode sequence-specific transcription factors, which likely help drive the myoblast and cerebellum specificity of other Myob/Cbl genes. This study extends our understanding of epigenetic/transcription associations related to differentiation and may help elucidate relationships between epigenetic signatures and muscular dystrophies or cerebellar-linked neuropathologies.

A comparison of clinical and pathologic assessments for the prediction of occult nipple involvement in nipple-sparing mastectomies.

BACKGROUND: Nipple-sparing mastectomy (NSM) for both risk reduction and cancer is increasing. In the cancer setting, most studies suggest the use of both clinical and intraoperative biopsy criteria in patient selection. This study examines the use of both biopsy and clinical criteria in women undergoing total nipple-removing mastectomy. METHODS: The study consisted of 58 patients undergoing total mastectomy without nipple sparing. Biopsies of the subareola tissue (SA), proximal nipple (NC) contents and radial sections of the residual nipple (NR) were examined microscopically. Tumor size and distance from the nipple were also noted. RESULTS: Using clinical criteria alone, the false negative rate was 53.8% and a false positive rate of 44.4%. When adding subareola and nipple core biopsies to clinical criteria the false negative rate fell to 7.7% but the false positive rate remained at 44.4%. When using only SA and NC biopsies to predict occult nipple involvement, the false negative rate was 11.8%. In 4 cases the NC was positive while the SA was negative for cancer and in 6 cases the SA was positive and NC negative. In 2 cases both the NC and SA biopsies were negative while the NR was positive. CONCLUSIONS: This study supports a more limited role in the use of clinical criteria for evaluating patients for NSM. This maximizes the number of patients who are candidates for NSM with minimal risk of nipple involvement. It was also noted that intraoperative biopsies are not totally reliable in predicting occult nipple involvement.

Quantitative proteomic analysis of the Anopheles gambiae (Diptera: Culicidae) midgut infected with o'nyong-nyong virus.

Alphaviruses are arthropod-borne pathogens that infect a range of hosts. In humans and other mammals, alphavirus infection can cause severe disease. In mosquito hosts, however, there are generally few symptoms. Little is known about the cellular responses of mosquitoes that allow them to cope with infection. In this investigation, a six-plex tandem mass tagging proteomic approach was used to study protein accumulation changes in the midgut of Anopheles gambiae (Giles) (Diptera: Culicidae) mosquitoes infected with o'nyong-nyong virus (Togaviridae, Alphavirus). Five hundred thirty-six nonredundant proteins were identified. Twenty-two were found in significantly different quantities in infected midguts compared with controls. Of interest, analysis revealed molecular pathways possibly targeted by virus proteins, such as those involving TAF4 and DNA polymerase phi proteins. Also identified was an FK506-binding protein. FK506-binding protein orthologs have been described as conserved host resistance factors, which suppress dengue and West Nile virus infection in human HeLa cells. This investigation constitutes the first study of the midgut-specific proteome of An. gambiae in relation to alphavirus infection. Our findings offer insight into mosquito immunity, including factors that possibly contribute to the different pathological outcomes observed in vertebrate and insect hosts.

TgMAPK1 is a Toxoplasma gondii MAP kinase that hijacks host MKK3 signals to regulate virulence and interferon-γ-mediated nitric oxide production.

The parasite Toxoplasma gondii controls tissue-specific nitric oxide (NO), thereby augmenting virulence and immunopathology through poorly-understood mechanisms. We now identify TgMAPK1, a Toxoplasma mitogen-activated protein kinase (MAPK), as a virulence factor regulating tissue-specific parasite burden by manipulating host interferon (IFN)-γ-mediated inducible nitric oxide synthase (iNOS). Toxoplasma with reduced TgMAPK1 expression (TgMAPK1(lo)) demonstrated that TgMAPK1 facilitates IFN-γ-driven p38 MAPK activation, reducing IFN-γ-generated NO in an MKK3-dependent manner, blunting IFN-γ-mediated parasite control. TgMAPK1(lo) infection in wild type mice produced ≥ten-fold lower parasite burden versus control parasites with normal TgMAPK1 expression (TgMAPK1(con)). Reduced parasite burdens persisted in IFN-γ KO mice, but equalized in normally iNOS-replete organs from iNOS KO mice. Parasite MAPKs are far less studied than other parasite kinases, but deserve additional attention as targets for immunotherapy and drug discovery.

DNA hypomethylation and hemimethylation in cancer.

In contrast to earlier views that there was much compartmentalization of the types of sequences subject to cancer-linked changes in DNA epigenetics, it is now clear that both cancer-associated DNA hypomethylation and hypermethylation are found throughout the genome. The hypermethylation includes promoters of tumor suppressor genes whose expression becomes repressed, thereby facilitating cancer formation. How hypomethylation contributes to carcinogenesis has been less clear. Recent insights into tissue-specific intra- and intergenic methylation and into cancer methylomes suggest that some of the DNA hypomethylation associated with cancers is likely to aid in tumor formation and progression by many different pathways, including effects on transcription in cis. Cancer-associated loss of DNA methylation from intergenic enhancers, promoter regions, silencers, and chromatin boundary elements may alter transcription rates. In -addition, cancer-associated intragenic DNA hypomethylation might modulate -alternative promoter usage, -production of intragenic noncoding RNA transcripts, cotranscriptional splicing, and transcription initiation or elongation. Initial studies of hemimethylation of DNA in cancer and many new studies of DNA demethylation in normal tissues suggest that active demethylation with spreading of hypomethylation can explain much of the cancer-associated DNA hypomethylation. The new discoveries that genomic 5-hydroxymethylcytosine is an intermediate in DNA demethylation, a base with its own functionality, and a modified base that, like 5-methylcytosine, exhibits cancer-associated losses, suggest that both decreased hydroxymethylation and decreased methylation of DNA play important roles in carcinogenesis.

Multivariate analysis of food consumption profiles in crisis settings.

Preventing malnutrition is one of the primary objectives of many humanitarian agencies, and household surveys are regularly employed to monitor food insecurity caused by political, economic, or environmental crises. Consumption frequencies for standard food groups are often collected to characterize the depth of food insecurity in a community and measure the impact of food assistance programs, producing a vector of bounded, correlated counts for each household. While aggregate indicators are typically used to summarize these results with a single statistic, they can be difficult to interpret and provide insufficient detail to judge the effectiveness of assistance programs. To address these limitations, we have developed a multivariate modeling framework for consumption frequency data. We introduce methods to update baseline models for the analysis of the smaller and more variable surveys typically collected in crisis settings, and we present an application of our approach to national consumption data collected in Yemen in 2014 and 2016 by the World Food Programme. The approach provides more nuanced and interpretable information about consumption changes in response to shocks and the effectiveness of humanitarian assistance.

Tumor necrosis factor receptor 1 functions as a tumor suppressor.

Tumor necrosis factor (TNF) is a key player in inflammatory bowel disease and has been variably associated with carcinogenesis, but details of the cross talk between inflammatory and tumorigenic pathways remain incompletely understood. It has been shown that, in C57BL/6 mice, signaling via TNF receptor 1 (TNFR1) is protective from injury and inflammation in experimental colitis. Therefore, we hypothesized that loss of TNFR1 signaling would confer increased risk of developing colitis-associated carcinoma. Using three models of murine tumorigenesis based on repeated bouts of inflammation or systemic tumor initiator, we sought to determine the roles of TNF and TNFR1 with regard to neoplastic transformation in the colon in wild-type (WT), TNFR1 knockout (R1KO), and TNF knockout (TNFKO) mice. We found R1KO animals to have more severe disease, as defined by weight loss, hematochezia, and histology. TNFKO mice demonstrated less weight loss but were consistently smaller, and rates and duration of hematochezia were comparable to WT mice. Histological inflammation scores were higher and neoplastic lesions occurred more frequently and earlier in R1KO mice. Apoptosis is not affected in R1KO mice although epithelial proliferation following injury is more ardent even before tumorigenesis is apparent. Lastly, there is earlier and more intense expression of activated ß-catenin in these mice, implying a connection between TNFR1 and Wnt signaling. Taken together, these findings show that in the context of colitis-associated carcinogenesis TNFR1 functions as a tumor suppressor, exerting this effect not via apoptosis but by modulating activation of ß-catenin and controlling epithelial proliferation.

Transcriptome and targetome analysis in MIR155 expressing cells using RNA-seq.

Previous studies have demonstrated the utility of microarray expression analysis to identify potential microRNA targets. Nevertheless, technical limitations intrinsic to this platform constrain its ability to fully exploit the potential of assessing transcript level changes to explore microRNA targetomes. High-throughput multiplexed Illumina-based next-generation sequencing (NGS) provides a digital readout of absolute transcript levels and imparts a higher level of accuracy and dynamic range than microarray platforms. We used Illumina NGS to analyze transcriptome changes induced by the human microRNA MIR155. This analysis resulted in a larger inferred targetome than similar studies carried out using microarray platforms. A comparison with 3' UTR reporter data demonstrated general concordance between NGS and corresponding 3' UTR reporter results. Nonharmonious results were investigated more deeply using transcript structure information assembled from the NGS data. This analysis revealed that transcript structure plays a substantial role in mitigated targeting and in frank targeting failures. With its high level of accuracy, its broad dynamic range, its utility in assessing transcript structure, and its capacity to accurately interrogate global direct and indirect transcriptome changes, NGS is a useful tool for investigating the biology and mechanisms of action of microRNAs.

Deciphering transcription dysregulation in FSH muscular dystrophy.

DUX4, a homeobox-containing gene present in a tandem array, is implicated in facioscapulohumeral muscular dystrophy (FSHD), a dominant autosomal disease. New findings about DUX4 have raised as many fundamental questions about the molecular pathology of this unique disease as they have answered. This review discusses recent studies addressing the question of whether there is extensive FSHD-related transcription dysregulation in adult-derived myoblasts and myotubes, the precursors for muscle repair. Two models for the role of DUX4 in FSHD are presented. One involves transient pathogenic expression of DUX4 in many cells in the muscle lineage before the myoblast stage resulting in a persistent, disease-related transcription profile ('Majority Rules'), which might be enhanced by subsequent oscillatory expression of DUX4. The other model emphasizes the toxic effects of inappropriate expression of DUX4 in only an extremely small percentage of FSHD myoblasts or myotube nuclei ('Minority Rules'). The currently favored Minority Rules model is not supported by recent studies of transcription dysregulation in FSHD myoblasts and myotubes. It also presents other difficulties, for example, explaining the expression of full-length DUX4 transcripts in FSHD fibroblasts. The Majority Rules model is the simpler explanation of findings about FSHD-associated gene expression and the DUX4-encoded homeodomain-type protein.

Endocrine disrupting chemicals in New Orleans surface waters and Mississippi Sound sediments.

Endocrine disrupting compounds (EDCs), represented by steroid hormones, organochlorine pesticides (OCPs), polychlorinated biphenyls (PCBs), and bisphenol A have been determined in four sediment cores from the Gulf of Mexico, from New Orleans surface water (Lake Pontchartrain and Mississippi River), and from the influent and effluent of a New Orleans municipal sewage treatment plant. During the five-month monitoring of selected EDCs in the Mississippi River (MR) and Lake Pontchartrain (LP) in 2008, 21 of 29 OCPs in MR and 17 of 29 OCPs in LP were detected; bisphenol A was detected in all of the samples. Steroid hormones (estrone, 17ß-estradiol and 17α-ethinylestradiol) were detected occasionally. Total EDC (OCPs + PCBs + steroid hormones + bisphenol A) concentrations in the two surface water samples were found to vary from 148 to 1112 ng L(-1). Strong correlation of the distribution of total OCPs, total PCBs and total EDCs between solid and water phases was found in LP, while moderate or no correlation existed in MR. OCPs, PCBs, steroid hormones, and bisphenol A were all detected in the ocean sediments, and total EDCs were measured in the range of 77 to 1796 ng g(-1) dry sediment weight. The EDCs were also found in untreated and treated municipal sewage samples with a removal efficiency of 83% for OCPs but no removal efficiency for 17α-ethinylestradiol.

Evaluating differences in Young's Modulus of regenerated and uninjured mouse digit bone through microCT density-based calculation and nanoindentation testing.

The mouse digit tip amputation model is an excellent model of bone regeneration, but its size and shape present an obstacle for biomechanical testing. As a result, assessing the structural quality of the regenerated bone in this model has focused on mineral density and bone architecture analysis. Here we describe an image-processing based method for assessment of mechanical properties in the regenerated digit by using micro-computed tomography mineral density data to calculate spatially discrete Young's modulus values throughout the entire distal third phalange. Further, we validate this method through comparison to nanoindentation-measured values for Young's modulus. Application to a set of regenerated and unamputated digits shows that regenerated bone has a lower Young's modulus compared to the uninjured digit, with a similar trend for experimental hardness values. Importantly, this method heightens the utility of the digit regeneration model, allows for more impactful treatment evaluation using the model, and introduces an analysis platform that can be used for other bones that do not conform to a standard long-bone model.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage.

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

Spatial transcriptomics reveals metabolic changes underly age-dependent declines in digit regeneration.

De novo limb regeneration after amputation is restricted in mammals to the distal digit tip. Central to this regenerative process is the blastema, a heterogeneous population of lineage-restricted, dedifferentiated cells that ultimately orchestrates regeneration of the amputated bone and surrounding soft tissue. To investigate skeletal regeneration, we made use of spatial transcriptomics to characterize the transcriptional profile specifically within the blastema. Using this technique, we generated a gene signature with high specificity for the blastema in both our spatial data, as well as other previously published single-cell RNA-sequencing transcriptomic studies. To elucidate potential mechanisms distinguishing regenerative from non-regenerative healing, we applied spatial transcriptomics to an aging model. Consistent with other forms of repair, our digit amputation mouse model showed a significant impairment in regeneration in aged mice. Contrasting young and aged mice, spatial analysis revealed a metabolic shift in aged blastema associated with an increased bioenergetic requirement. This enhanced metabolic turnover was associated with increased hypoxia and angiogenic signaling, leading to excessive vascularization and altered regenerated bone architecture in aged mice. Administration of the metabolite oxaloacetate decreased the oxygen consumption rate of the aged blastema and increased WNT signaling, leading to enhanced in vivo bone regeneration. Thus, targeting cell metabolism may be a promising strategy to mitigate aging-induced declines in tissue regeneration.

Acceptability and response to a postal survey using self-taken samples for HPV vaccine impact monitoring.

OBJECTIVE: To assess the feasibility and acceptance of a postal survey to measure human papillomavirus (HPV) prevalence and monitor vaccine impact, using self-taken specimens from young women who do not attend their first cervical screening appointment. METHODS: Focus groups informed the survey design identifying factors that would influence acceptability. Postal testing kits were sent to a nationally representative sample of unscreened women. Overall response rate, the influence of different specimen types (urine or vaginal swab) and the receipt of a reminder letter on participation were calculated. Specimens were tested anonymously for HPV. Individual test results were not provided. RESULTS: Of 5500 kits sent, 725 were returned (13.2%). Fifty-two women actively opted out. There was a higher return rate for urine kits (13.7% vs 12%) and from those who received a reminder letter (15.5% vs 12.2%). Response was influenced by deprivation (10.3% in the most deprived quintile vs 16.2% in the least). Overall weighted HPV prevalence was 35.9% (40.0% from swab specimens and 31.9% from urine). CONCLUSIONS: Some women were willing to participate in anonymised postal testing. However, the low uptake means that HPV prevalence results are difficult to interpret for ongoing surveillance. Monitoring HPV vaccine impact outwith the cervical screening programme remains challenging.

In vitro high-capacity assay to quantify the clonal heterogeneity in trilineage potential of mesenchymal stem cells reveals a complex hierarchy of lineage commitment.

In regenerative medicine, bone marrow is a promising source of mesenchymal stem cells (MSCs) for a broad range of cellular therapies. This research addresses a basic prerequisite to realize the therapeutic potential of MSCs by developing a novel high-capacity assay to quantify the clonal heterogeneity in potency that is inherent to MSC preparations. The assay utilizes a 96-well format to (1) classify MSCs according to colony-forming efficiency as a measure of proliferation capacity and trilineage potential to exhibit adipo-, chondro-, and osteogenesis as a measure of multipotency and (2) preserve a frozen template of MSC clones of known potency for future use. The heterogeneity in trilineage potential of normal bone marrow MSCs is more complex than previously reported: all eight possible categories of trilineage potential were detected. In this study, the average colony-forming efficiency of MSC preparations was 55-62%, and tripotent MSCs accounted for nearly 50% of the colony-forming cells. The multiple phenotypes detected in this study infer a more convoluted hierarchy of lineage commitment than described in the literature. Greater cell amplification, colony-forming efficiency, and colony diameter for tri- versus unipotent clones suggest that MSC proliferation may be a function of potency. CD146 may be a marker of multipotency, with approximately 2-fold difference in mean fluorescence intensity between tri- and unipotent clones. The significance of these findings is discussed in the context of the efficacy of MSC therapies. The in vitro assay described herein will likely have numerous applications given the importance of heterogeneity to the therapeutic potential of MSCs.

Clonal analysis of the proliferation potential of human bone marrow mesenchymal stem cells as a function of potency.

Human mesenchymal stem cells (MSCs) from bone marrow are a heterogeneous ensemble of progenitors and lineage-committed cells, with a broad range of regenerative properties. Ex vivo expansion to produce sufficient quantities of MSCs is essential for most therapeutic applications. The present study resolves the relationship between proliferation potential of MSCs and their potency. Clonal analysis generated single-cell derived colonies of MSCs that were classified according to their trilineage potential to exhibit adipo- (A), chondro- (C), and osteogenesis (O) as a measure of potency. Multipotent OAC clones were highly proliferative with colony-forming efficiencies that ranged from 35% to 90%; whereas, O clones formed colonies with an efficiency of 5% or less (P < 0.01). Similar trends were evident during ex vivo expansion: for example, the median specific growth rate was 0.8 day(-1) (20 h doubling time) for cultures inoculated with OAC clones and was 5-fold less for inocula of O clones (P < 0.01). OA and OC clones had similar proliferation potentials. More than 75% of cells in subconfluent cultures inoculated with O clones stained positive for senescence-associated ß-galactosidase activity vs. less than 10% for OAC clones (P < 0.001). Apoptotic cells were in the minority for all potency groups. Preliminary data generated during clonal analysis suggest that osteogenic potential of MSCs to produce mineralized matrix is a function of potency, as well. These results are discussed in the context of the preparation of efficacious MSC therapies by ex vivo expansion.