Engineering anisotropic biomimetic fibrocartilage microenvironment by bioprinting mesenchymal stem cells in nanoliter gel droplets.

Over the past decade, bioprinting has emerged as a promising patterning strategy to organize cells and extracellular components both in two and three dimensions (2D and 3D) to engineer functional tissue mimicking constructs. So far, tissue printing has neither been used for 3D patterning of mesenchymal stem cells (MSCs) in multiphase growth factor embedded 3D hydrogels nor been investigated phenotypically in terms of simultaneous differentiation into different cell types within the same micropatterned 3D tissue constructs. Accordingly, we demonstrated a biochemical gradient by bioprinting nanoliter droplets encapsulating human MSCs, bone morphogenetic protein 2 (BMP-2), and transforming growth factor ß1 (TGF- ß1), engineering an anisotropic biomimetic fibrocartilage microenvironment. Assessment of the model tissue construct displayed multiphasic anisotropy of the incorporated biochemical factors after patterning. Quantitative real time polymerase chain reaction (qRT-PCR) results suggested genomic expression patterns leading to simultaneous differentiation of MSC populations into osteogenic and chondrogenic phenotype within the multiphasic construct, evidenced by upregulation of osteogenesis and condrogenesis related genes during in vitro culture. Comprehensive phenotypic network and pathway analysis results, which were based on genomic expression data, indicated activation of differentiation related mechanisms, via signaling pathways, including TGF, BMP, and vascular endothelial growth factor.

Two developmental modules establish 3D beak-shape variation in Darwin's finches.

Bird beaks display tremendous variation in shape and size, which is closely associated with the exploitation of multiple ecological niches and likely played a key role in the diversification of thousands of avian species. Previous studies have demonstrated some of the molecular mechanisms that regulate morphogenesis of the prenasal cartilage, which forms the initial beak skeleton. However, much of the beak diversity in birds depends on variation in the premaxillary bone. It forms later in development and becomes the most prominent functional and structural component of the adult upper beak/jaw, yet its regulation is unknown. Here, we studied a group of Darwin's finch species with different beak shapes. We found that TGFßIIr, ß-catenin, and Dickkopf-3, the top candidate genes from a cDNA microarray screen, are differentially expressed in the developing premaxillary bone of embryos of species with different beak shapes. Furthermore, our functional experiments demonstrate that these molecules form a regulatory network governing the morphology of the premaxillary bone, which differs from the network controlling the prenasal cartilage, but has the same species-specific domains of expression. These results offer potential mechanisms that may explain how the tightly coupled depth and width dimensions can evolve independently. The two-module program of development involving independent regulating molecules offers unique insights into how different developmental pathways may be modified and combined to induce multidimensional shifts in beak morphology. Similar modularity in development may characterize complex traits in other organisms to a greater extent than is currently appreciated.

Comparison of in-gel protein separation techniques commonly used for fractionation in mass spectrometry-based proteomic profiling.

Fractionation of complex samples at the cellular, subcellular, protein, or peptide level is an indispensable strategy to improve the sensitivity in mass spectrometry-based proteomic profiling. This study revisits, evaluates, and compares the most common gel-based protein separation techniques i.e. 1D SDS-PAGE, 1D preparative SDS-PAGE, IEF-IPG, and 2D-PAGE in their performance as fractionation approaches in nano LC-ESI-MS/MS analysis of a mixture of protein standards and mitochondrial extracts isolated from rat liver. This work demonstrates that all the above techniques provide complementary protein identification results, but 1D SDS-PAGE and IEF-IPG had the highest number of identifications. The IEF-IPG technique resulted in the highest average number of detected peptides per protein. The 2D-PAGE was evaluated as a protein fractionation approach. This work shows that the recovery of proteins and resulting proteolytic digests is highly dependent on the total volume of the gel matrix. The performed comparison of the fractionation techniques demonstrates the potential of a combination of orthogonal 1D SDS-PAGE and IEF-IPG for the improved sensitivity of profiling without significant decrease in throughput.

The calmodulin pathway and evolution of elongated beak morphology in Darwin's finches.

A classic textbook example of adaptive radiation under natural selection is the evolution of 14 closely related species of Darwin's finches (Fringillidae, Passeriformes), whose primary diversity lies in the size and shape of their beaks. Thus, ground finches have deep and wide beaks, cactus finches have long and pointed beaks (low depth and narrower width), and warbler finches have slender and pointed beaks, reflecting differences in their respective diets. Previous work has shown that even small differences in any of the three major dimensions (depth, width and length) of the beak have major consequences for the overall fitness of the birds. Recently we used a candidate gene approach to explain one pathway involved in Darwin's finch beak morphogenesis. However, this type of analysis is limited to molecules with a known association with craniofacial and/or skeletogenic development. Here we use a less constrained, complementary DNA microarray analysis of the transcripts expressed in the beak primordia to find previously unknown genes and pathways whose expression correlates with specific beak morphologies. We show that calmodulin (CaM), a molecule involved in mediating Ca2+ signalling, is expressed at higher levels in the long and pointed beaks of cactus finches than in more robust beak types of other species. We validated this observation with in situ hybridizations. When this upregulation of the CaM-dependent pathway is artificially replicated in the chick frontonasal prominence, it causes an elongation of the upper beak, recapitulating the beak morphology of the cactus finches. Our results indicate that local upregulation of the CaM-dependent pathway is likely to have been a component of the evolution of Darwin's finch species with elongated beak morphology and provide a mechanistic explanation for the independence of beak evolution along different axes. More generally, our results implicate the CaM-dependent pathway in the developmental regulation of craniofacial skeletal structures.

AnyExpress: integrated toolkit for analysis of cross-platform gene expression data using a fast interval matching algorithm.

BACKGROUND: Cross-platform analysis of gene express data requires multiple, intricate processes at different layers with various platforms. However, existing tools handle only a single platform and are not flexible enough to support custom changes, which arise from the new statistical methods, updated versions of reference data, and better platforms released every month or year. Current tools are so tightly coupled with reference information, such as reference genome, transcriptome database, and SNP, which are often erroneous or outdated, that the output results are incorrect and misleading. RESULTS: We developed AnyExpress, a software package that combines cross-platform gene expression data using a fast interval-matching algorithm. Supported platforms include next-generation-sequencing technology, microarray, SAGE, MPSS, and more. Users can define custom target transcriptome database references for probe/read mapping in any species, as well as criteria to remove undesirable probes/reads. AnyExpress offers scalable processing features such as binding, normalization, and summarization that are not present in existing software tools. As a case study, we applied AnyExpress to published Affymetrix microarray and Illumina NGS RNA-Seq data from human kidney and liver. The mean of within-platform correlation coefficient was 0.98 for within-platform samples in kidney and liver, respectively. The mean of cross-platform correlation coefficients was 0.73. These results confirmed those of the original and secondary studies. Applying filtering produced higher agreement between microarray and NGS, according to an agreement index calculated from differentially expressed genes. CONCLUSION: AnyExpress can combine cross-platform gene expression data, process data from both open- and closed-platforms, select a custom target reference, filter out undesirable probes or reads based on custom-defined biological features, and perform quantile-normalization with a large number of microarray samples. AnyExpress is fast, comprehensive, flexible, and freely available at http://anyexpress.sourceforge.net.

Thermal stabilization of tissues and the preservation of protein phosphorylation states for two-dimensional gel electrophoresis.

2-DE is typically capable of discriminating proteins differing by a single phosphorylation or dephosphorylation event. However, a reliable representation of protein phosphorylation states as they occur in vivo requires that both phosphatases and kinases are rapidly and completely inactivated. Thermal stabilization of mouse cerebral cortex homogenates effectively inactivated these enzymes, as evidenced by comparison with unstabilized tissues where abscissal pI shifts were a common feature in 2-D gels. Of the 588 matched proteins separated on 2-D gels comparing stabilized and unstabilized tissues, 53 proteins exhibited greater than twofold differences in spot volume (ANOVA, p<0.05). Phosphoprotein-specific staining was corroborated by the identification of 16 phosphoproteins by nano-LC MS/MS and phosphotyrosine kinase activity assay.

Error margin analysis for feature gene extraction.

BACKGROUND: Feature gene extraction is a fundamental issue in microarray-based biomarker discovery. It is normally treated as an optimization problem of finding the best predictive feature genes that can effectively and stably discriminate distinct types of disease conditions, e.g. tumors and normals. Since gene microarray data normally involves thousands of genes at, tens or hundreds of samples, the gene extraction process may fall into local optimums if the gene set is optimized according to the maximization of classification accuracy of the classifier built from it. RESULTS: In this paper, we propose a novel gene extraction method of error margin analysis to optimize the feature genes. The proposed algorithm has been tested upon one synthetic dataset and two real microarray datasets. Meanwhile, it has been compared with five existing gene extraction algorithms on each dataset. On the synthetic dataset, the results show that the feature set extracted by our algorithm is the closest to the actual gene set. For the two real datasets, our algorithm is superior in terms of balancing the size and the validation accuracy of the resultant gene set when comparing to other algorithms. CONCLUSION: Because of its distinct features, error margin analysis method can stably extract the relevant feature genes from microarray data for high-performance classification.

A potential role for intragenic miRNAs on their hosts' interactome.

BACKGROUND: miRNAs are small, non-coding RNA molecules that mainly act as negative regulators of target gene messages. Due to their regulatory functions, they have lately been implicated in several diseases, including malignancies. Roughly half of known miRNA genes are located within previously annotated protein-coding regions ("intragenic miRNAs"). Although a role of intragenic miRNAs as negative feedback regulators has been speculated, to the best of our knowledge there have been no conclusive large-scale studies investigating the relationship between intragenic miRNAs and host genes and their pathways. RESULTS: miRNA-containing host genes were three times longer, contained more introns and had longer 5' introns compared to a randomly sampled gene cohort. These results are consistent with the observation that more than 60% of intronic miRNAs are found within the first five 5' introns. Host gene 3'-untranslated regions (3'-UTRs) were 40% longer and contained significantly more adenylate/uridylate-rich elements (AREs) compared to a randomly sampled gene cohort. Coincidentally, recent literature suggests that several components of the miRNA biogenesis pathway are required for the rapid decay of mRNAs containing AREs. A high-confidence set of predicted mRNA targets of intragenic miRNAs also shared many of these features with the host genes. Approximately 20% of intragenic miRNAs were predicted to target their host mRNA transcript. Further, KEGG pathway analysis demonstrated that 22 of the 74 pathways in which host genes were associated showed significant overrepresentation of proteins encoded by the mRNA targets of associated intragenic miRNAs. CONCLUSIONS: Our findings suggest that both host genes and intragenic miRNA targets may potentially be subject to multiple layers of regulation. Tight regulatory control of these genes is likely critical for cellular homeostasis and absence of disease. To this end, we examined the potential for negative feedback loops between intragenic miRNAs, host genes, and miRNA target genes. We describe, how higher-order miRNA feedback on hosts' interactomes may at least in part explain correlation patterns observed between expression of host genes and intragenic miRNA targets in healthy and tumor tissue.

Effects of nicotine on gene expression and osseointegration in rats.

BACKGROUND: While many studies have focused on the hazardous effects of smoking, there is little direct evidence regarding the specific detrimental effects of the nicotine on the osseointegration of implants. OBJECTIVE: To understand the effects of nicotine on gene expression and osseointegration of titanium implants in rats. MATERIAL AND METHODS: Forty-four rats were administered with nicotine or saline for a period of 8 weeks. The femurs were then harvested and analyzed using a three-point bending test. Osseointegration level was determined using bone/implant contact ratio at 2 or 4 weeks after implants were placed. Expression levels of bone matrix-related genes were measured by quantitative real-time polymerase chain reaction. RESULTS: The results of the three-point bending showed that there was no significant difference detected in stiffness between control and nicotine groups at 8 weeks post-saline/nicotine delivery (P=0.705). The bone/implant contact ratio in nicotine-delivered group was significantly decreased compared with those in the control group at 4 weeks (P<0.05). Also, expression levels of osteopontin, type II collagen, bone morphogenic protein-2, bone sialoprotein, and core-binding factor α-1 were significantly down-regulated in the nicotine-delivered group compared with the control. CONCLUSIONS: Although systemic exposure to nicotine did not affect rat bone development, bone wound healing around the implant after placement was affected. These findings suggest that nicotine might inhibit the bone matrix-related gene expressions required for wound healing and thereby diminish implant osseointegration at late stage.

Validation of oligoarrays for quantitative exploration of the transcriptome.

BACKGROUND: Oligoarrays have become an accessible technique for exploring the transcriptome, but it is presently unclear how absolute transcript data from this technique compare to the data achieved with tag-based quantitative techniques, such as massively parallel signature sequencing (MPSS) and serial analysis of gene expression (SAGE). By use of the TransCount method we calculated absolute transcript concentrations from spotted oligoarray intensities, enabling direct comparisons with tag counts obtained with MPSS and SAGE. The tag counts were converted to number of transcripts per cell by assuming that the sum of all transcripts in a single cell was 5.105. Our aim was to investigate whether the less resource demanding and more widespread oligoarray technique could provide data that were correlated to and had the same absolute scale as those obtained with MPSS and SAGE. RESULTS: A number of 1,777 unique transcripts were detected in common for the three technologies and served as the basis for our analyses. The correlations involving the oligoarray data were not weaker than, but, similar to the correlation between the MPSS and SAGE data, both when the entire concentration range was considered and at high concentrations. The data sets were more strongly correlated at high transcript concentrations than at low concentrations. On an absolute scale, the number of transcripts per cell and gene was generally higher based on oligoarrays than on MPSS and SAGE, and ranged from 1.6 to 9,705 for the 1,777 overlapping genes. The MPSS data were on same scale as the SAGE data, ranging from 0.5 to 3,180 (MPSS) and 9 to1,268 (SAGE) transcripts per cell and gene. The sum of all transcripts per cell for these genes was 3.8.105 (oligoarrays), 1.1.105 (MPSS) and 7.6.104 (SAGE), whereas the corresponding sum for all detected transcripts was 1.1.106 (oligoarrays), 2.8.105 (MPSS) and 3.8.105 (SAGE). CONCLUSION: The oligoarrays and TransCount provide quantitative transcript concentrations that are correlated to MPSS and SAGE data, but, the absolute scale of the measurements differs across the technologies. The discrepancy questions whether the sum of all transcripts within a single cell might be higher than the number of 5.105 suggested in the literature and used to convert tag counts to transcripts per cell. If so, this may explain the apparent higher transcript detection efficiency of the oligoarrays, and has to be clarified before absolute transcript concentrations can be interchanged across the technologies. The ability to obtain transcript concentrations from oligoarrays opens up the possibility of efficient generation of universal transcript databases with low resource demands.

Automatic correspondence of tags and genes (ACTG): a tool for the analysis of SAGE, MPSS and SBS data.

UNLABELLED: A critical step in any SAGE, MPSS and SBS data analysis is tag-to-gene assignment. Current available tools are limited by a tag-by-tag annotation process and/or do not provide the dataset that is used to produce a complete tag-to-gene mapping. We developed ACTG, a web-based application that allows a large-scale tag-to-gene mapping using several reference datasets. ACTG can annotate SAGE (14 or 21 bp), MPSS (17 or 20 bp) and SBS (16 bp) data for both human and mouse organisms. AVAILABILITY: http://retina.med.harvard.edu/ACTG/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genomic analysis of mouse retinal development.

The vertebrate retina is comprised of seven major cell types that are generated in overlapping but well-defined intervals. To identify genes that might regulate retinal development, gene expression in the developing retina was profiled at multiple time points using serial analysis of gene expression (SAGE). The expression patterns of 1,051 genes that showed developmentally dynamic expression by SAGE were investigated using in situ hybridization. A molecular atlas of gene expression in the developing and mature retina was thereby constructed, along with a taxonomic classification of developmental gene expression patterns. Genes were identified that label both temporal and spatial subsets of mitotic progenitor cells. For each developing and mature major retinal cell type, genes selectively expressed in that cell type were identified. The gene expression profiles of retinal Müller glia and mitotic progenitor cells were found to be highly similar, suggesting that Müller glia might serve to produce multiple retinal cell types under the right conditions. In addition, multiple transcripts that were evolutionarily conserved that did not appear to encode open reading frames of more than 100 amino acids in length ("noncoding RNAs") were found to be dynamically and specifically expressed in developing and mature retinal cell types. Finally, many photoreceptor-enriched genes that mapped to chromosomal intervals containing retinal disease genes were identified. These data serve as a starting point for functional investigations of the roles of these genes in retinal development and physiology.

Analysis of repeatability in spotted cDNA microarrays.

We report a strategy for analysis of data quality in cDNA microarrays based on the repeatability of repeatedly spotted clones. We describe how repeatability can be used to control data quality by developing adaptive filtering criteria for microarray data containing clones spotted in multiple spots. We have applied the method on five publicly available cDNA microarray data sets and one previously unpublished data set from our own laboratory. The results demonstrate the feasibility of the approach as a foundation for data filtering, and indicate a high degree of variation in data quality, both across the data sets and between arrays within data sets.

Endothelial and circulating C19MC microRNAs are biomarkers of infantile hemangioma.

Infantile hemangioma (IH) is the most common vascular tumor of infancy, and it uniquely regresses in response to oral propranolol. MicroRNAs (miRNAs) have emerged as key regulators of vascular development and are dysregulated in many disease processes, but the role of miRNAs in IH growth has not been investigated. We report expression of C19MC, a primate-specific megacluster of miRNAs expressed in placenta with rare expression in postnatal tissues, in glucose transporter 1-expressing (GLUT-1-expressing) IH endothelial cells and in the plasma of children with IH. Tissue or circulating C19MC miRNAs were not detectable in patients having 9 other types of vascular anomalies or unaffected children, identifying C19MC miRNAs as the first circulating biomarkers of IH. Levels of circulating C19MC miRNAs correlated with IH tumor size and propranolol treatment response, and IH tissue from children treated with propranolol or from children with partially involuted tumors contained lower levels of C19MC miRNAs than untreated, proliferative tumors, implicating C19MC miRNAs as potential drivers of IH pathogenesis. Detection of C19MC miRNAs in the circulation of infants with IH may provide a specific and noninvasive means of IH diagnosis and identification of candidates for propranolol therapy as well as a means to monitor treatment response.

Epithelialization of mouse ovarian tumor cells originating in the fallopian tube stroma.

Epithelial ovarian carcinoma accounts for 90% of all ovarian cancer and is the most deadly gynecologic malignancy. Recent studies have suggested that fallopian tube fimbriae can be the origin of cells for high-grade serous subtype of epithelial ovarian carcinoma (HGSOC). A mouse HGSOC model with conditional Dicer-Pten double knockout (Dicer-Pten DKO) developed primary tumors, intriguingly, from the fallopian tube stroma. We examined the growth and epithelial phenotypes of the Dicer-Pten DKO mouse tumor cells contributable by each gene knockout. Unlike human ovarian epithelial cancer cells that expressed full-length E-cadherin, the Dicer-Pten DKO stromal tumor cells expressed cleaved E-cadherin fragments and metalloproteinase 2, a mixture of epithelial and mesenchymal markers. Although the Dicer-Pten DKO tumor cells lost the expression of mature microRNAs as expected, they showed high levels of tRNA fragment expression and enhanced AKT activation due to the loss of PTEN function. Introduction of a Dicer1-expressing construct into the DKO mouse tumor cells significantly reduced DNA synthesis and the cell growth rate, with concurrent diminished adhesion and ZO1 epithelial staining. Hence, it is likely that the loss of Dicer promoted mesenchymal-epithelial transition in fallopian tube stromal cells, and in conjunction with Pten loss, further promoted cell proliferation and epithelial-like tumorigenesis.

Infected erythrocyte-derived extracellular vesicles alter vascular function via regulatory Ago2-miRNA complexes in malaria.

Malaria remains one of the greatest public health challenges worldwide, particularly in sub-Saharan Africa. The clinical outcome of individuals infected with Plasmodium falciparum parasites depends on many factors including host systemic inflammatory responses, parasite sequestration in tissues and vascular dysfunction. Production of pro-inflammatory cytokines and chemokines promotes endothelial activation as well as recruitment and infiltration of inflammatory cells, which in turn triggers further endothelial cell activation and parasite sequestration. Inflammatory responses are triggered in part by bioactive parasite products such as hemozoin and infected red blood cell-derived extracellular vesicles (iRBC-derived EVs). Here we demonstrate that such EVs contain functional miRNA-Argonaute 2 complexes that are derived from the host RBC. Moreover, we show that EVs are efficiently internalized by endothelial cells, where the miRNA-Argonaute 2 complexes modulate target gene expression and barrier properties. Altogether, these findings provide a mechanistic link between EVs and vascular dysfunction during malaria infection.

Molecular changes in endometriosis-associated ovarian clear cell carcinoma.

BACKGROUND: Endometriosis is frequently associated with and thought of having propensity to develop into ovarian clear cell carcinoma (OCCC), although the molecular transformation mechanism is not completely understood. METHODS: We employed immunohistochemical (IHC) staining for marker expression along the potential progression continuum. Expression profiling of microdissected endometriotic and OCCC cells from patient-matched formalin-fixed, paraffin-embedded samples was performed to explore the carcinogenic pathways. Function of novel biomarkers was confirmed by knockdown experiments. RESULTS: PTEN was significantly lost in both endometriosis and invasive tumour tissues, while oestrogen receptor (ER) expression was lost in OCCC relative to endometriosis. XRCC5, PTCH2, eEF1A2 and PPP1R14B were significantly overexpressed in OCCC and associated endometriosis, but not in benign endometriosis (p ⩽ 0.004). Knockdown experiments with XRCC5 and PTCH2 in a clear cell cancer cell line resulted in significant growth inhibition. There was also significant silencing of a panel of target genes with histone H3 lysine 27 trimethylation, a signature of polycomb chromatin-remodelling complex in OCCC. IHC confirmed the loss of expression of one such polycomb target gene, the serous ovarian cancer lineage marker Wilms' tumour protein 1 (WT1) in OCCC, while endometriotic tissues showed significant co-expression of WT1 and ER. CONCLUSIONS: Loss of PTEN expression is proposed as an early and permissive event in endometriosis development, while the loss of ER and polycomb-mediated transcriptional reprogramming for pluripotency may play an important role in the ultimate transformation process. Our study provides new evidence to redefine the pathogenic programme for lineage-specific transformation of endometriosis to OCCC.

Gene expression profiling using laser capture microdissection.

Human disease is governed by a complex array of cellular populations and sub-population. Gene expression profiling is proving an important means of understanding and classifying pathophysiologic processes by identifying genes, gene pathways and pathway networks not previously known to be associated with particular diseases. However, disease-associated gene expression can be obscured by surrounding 'normal' tissue. Laser capture microdissection allows gene expression analysis of pooled single cells, cell subpopulations and cell populations. Analysis of laser capture microdissection-procured cells will allow a better understanding of the cellular components of disease.

The Goodman-Kruskal coefficient and its applications in genetic diagnosis of cancer.

Increasing interest in new pattern recognition methods has been motivated by bioinformatics research. The analysis of gene expression data originated from microarrays constitutes an important application area for classification algorithms and illustrates the need for identifying important predictors. We show that the Goodman-Kruskal coefficient can be used for constructing minimal classifiers for tabular data, and we give an algorithm that can construct such classifiers.

MicroRNA-146a and microRNA-155 show tissue-dependent expression in dental pulp, gingival and periodontal ligament fibroblasts in vitro.

MicroRNAs (miRNAs) are small non-coding RNAs showing a tissue-specific expression pattern, and whose function is to suppress protein synthesis. In this study, we hypothesized that expression of miRNAs would differ among fibroblasts from dental pulp (DPF), gingiva (GF) and periodontal ligament (PLF) in vitro. Once established by an explant technique, DPF, GF and PLF were collected for RNA isolation and subjected to a miRNA microarray. Next, cells were stimulated with E. coli lipopolysaccharide (LPS) for 24 h and then collected for RNA isolation. Expression of miR-146a and miR-155 was investigated by qPCR. Microarray screening revealed several miRNAs that showed specifically high expression in at least one of the fibroblast subtypes. These molecules are potentially involved in the regulation of extracellular matrix turnover and production of inflammatory mediators. Microarray analysis showed that both miR-146a and miR-155 were among the miRNAs expressed exclusively in GF. qPCR demonstrated significant upregulation of miR-146a only in GF after LPS stimulation, whereas basal expression of miR-155 was higher in GF than in the other cell subtypes. LPS downregulated the expression of miR-155 only in GF. Our results suggest that the expression and regulation of miR-146a and miR-155 are more pronounced in GF than in DPF and PLF.