Human bone marrow contains a population of non-hematopoietic stromal stem/progenitor cells (BMSCs), which play a central role for bone marrow stroma and the hematopoietic microenvironment. However, the precise characteristics and potential stem cell properties of defined BMSC populations have not yet been thoroughly investigated. Using standard adherent colony-forming unit fibroblast (CFU-F) assays, we have previously shown that BMSCs were highly enriched in the nonhematopoietic CD271pos/CD140alow/neg fraction of normal adult human bone marrow. In this study, we demonstrate that prospectively isolated CD271pos/CD140alow/neg BMSCs expressed high levels of hematopoiesis supporting genes and signature mesenchymal and multipotency genes on a single cell basis. Furthermore, CD271pos/CD140alow/neg BMSCs gave rise to non-adherent sphere colonies (mesenspheres) with typical surface marker profile and trilineage in vitro differentiation potential. Importantly, serial transplantations of CD271pos/CD140alow/neg BMSC-derived mesenspheres (single cell and bulk) into immunodeficient NOD scid gamma (NSG) mice showed increased mesensphere numbers and full differentiation potential after both primary and secondary transplantations. In contrast, BMSC self-renewal potential decreased under standard adherent culture conditions. These data therefore indicate that CD271pos/CD140alow/neg BMSCs represent a population of primary stem cells with MSC phenotype and sphere-forming capacity that fulfill stringent functional stem cell criteria in vivo in a serial transplantation setting.

Coxsackievirus A16 (CA16) is a member of the Picornaviridae family and causes mild and self-limiting hand, foot, and mouth disease (HFMD) in infants and young children. CA16 infection can also progress to central nervous system (CNS) complications; however, the underlying mechanism by which CA16 penetrates the blood-brain barrier (BBB) and then causes CNS damage remains unclear. This study aimed to explore the mechanism of CA16 neurotropic tropism by establishing an in vitro BBB model with CA16 infection and an in vivo CA16 rhesus monkey infant infection model. The results showed that CA16 infection induced increased permeability of the BBB accompanied by upregulation of matrix metalloproteinase 9 (MMP9) expression. Subsequently, high-throughput miRNA sequencing technology and bioinformatics analysis revealed that miR-1303 may regulate BBB permeability by targeting MMP9. Next, we used dual-luciferase, qRT-PCR, and western blot assays to provide evidence of MMP9 targeting by miR-1303. Further experiments revealed that CA16 infection promoted the degradation of junctional complexes (Claudin4, Claudin5, VE-Cadherin, and ZO-1), likely by downregulating miR-1303 and upregulating MMP9. Finally, EGFP-CA16 infection could enter the CNS by facilitating the degradation of junctional complexes, eventually causing neuroinflammation and injury to the CNS, which was confirmed using the in vivo rhesus monkey model. Our results indicate that CA16 might penetrate the BBB and then enter the CNS by downregulating miR-1303, which disrupts junctional complexes by directly regulating MMP9 and ultimately causing pathological CNS changes. These results provide new therapeutic targets in HFMD patients following CA16 infection.

Chronic kidney disease (CKD) is a complex gene-environmental disease affecting close to 10% of the US population. Genome-wide association studies (GWASs) have identified sequence variants, localized to non-coding genomic regions, associated with kidney function. Despite these robust observations, the mechanism by which variants lead to CKD remains a critical unanswered question. Expression quantitative trait loci (eQTL) analysis is a method to identify genetic variation associated with gene expression changes in specific tissue types. We hypothesized that an integrative analysis combining CKD GWAS and kidney eQTL results can identify candidate genes for CKD. We performed eQTL analysis by correlating genotype with RNA-seq-based gene expression levels in 96 human kidney samples. Applying stringent statistical criteria, we detected 1,886 genes whose expression differs with the sequence variants. Using direct overlap and Bayesian methods, we identified new potential target genes for CKD. With respect to one of the target genes, lysosomal beta A mannosidase (MANBA), we observed that genetic variants associated with MANBA expression in the kidney showed statistically significant colocalization with variants identified in CKD GWASs, indicating that MANBA is a potential target gene for CKD. The expression of MANBA was significantly lower in kidneys of subjects with risk alleles. Suppressing manba expression in zebrafish resulted in renal tubule defects and pericardial edema, phenotypes typically induced by kidney dysfunction. Our analysis shows that gene-expression changes driven by genetic variation in the kidney can highlight potential new target genes for CKD development.

Copy number variations (CNVs) are associated with many complex diseases. Next generation sequencing data enable one to identify precise CNV breakpoints to better under the underlying molecular mechanisms and to design more efficient assays. Using the CIGAR strings of the reads, we develop a method that can identify the exact CNV breakpoints, and in cases when the breakpoints are in a repeated region, the method reports a range where the breakpoints can slide. Our method identifies the breakpoints of a CNV using both the positions and CIGAR strings of the reads that cover breakpoints of a CNV. A read with a long soft clipped part (denoted as S in CIGAR) at its 3'(right) end can be used to identify the 5'(left)-side of the breakpoints, and a read with a long S part at the 5' end can be used to identify the breakpoint at the 3'-side. To ensure both types of reads cover the same CNV, we require the overlapped common string to include both of the soft clipped parts. When a CNV starts and ends in the same repeated regions, its breakpoints are not unique, in which case our method reports the left most positions for the breakpoints and a range within which the breakpoints can be incremented without changing the variant sequence. We have implemented the methods in a C++ package intended for the current Illumina Miseq and Hiseq platforms for both whole genome and exon-sequencing. Our simulation studies have shown that our method compares favorably with other similar methods in terms of true discovery rate, false positive rate and breakpoint accuracy. Our results from a real application have shown that the detected CNVs are consistent with zygosity and read depth information. The software package is available at http://statgene.med.upenn.edu/softprog.html.

Intense interest centers on the role of the human gut microbiome in health and disease, but optimal methods for analysis are still under development. Here we present a study of methods for surveying bacterial communities in human feces using 454/Roche pyrosequencing of 16S rRNA gene tags. We analyzed fecal samples from 10 individuals and compared methods for storage, DNA purification and sequence acquisition. To assess reproducibility, we compared samples one cm apart on a single stool specimen for each individual. To analyze storage methods, we compared 1) immediate freezing at -80 degrees C, 2) storage on ice for 24 or 3) 48 hours. For DNA purification methods, we tested three commercial kits and bead beating in hot phenol. Variations due to the different methodologies were compared to variation among individuals using two approaches--one based on presence-absence information for bacterial taxa (unweighted UniFrac) and the other taking into account their relative abundance (weighted UniFrac). In the unweighted analysis relatively little variation was associated with the different analytical procedures, and variation between individuals predominated. In the weighted analysis considerable variation was associated with the purification methods. Particularly notable was improved recovery of Firmicutes sequences using the hot phenol method. We also carried out surveys of the effects of different 454 sequencing methods (FLX versus Titanium) and amplification of different 16S rRNA variable gene segments. Based on our findings we present recommendations for protocols to collect, process and sequence bacterial 16S rDNA from fecal samples--some major points are 1) if feasible, bead-beating in hot phenol or use of the PSP kit improves recovery; 2) storage methods can be adjusted based on experimental convenience; 3) unweighted (presence-absence) comparisons are less affected by lysis method.

Estrogen and insulin-like growth factor-1 (IGF-1) play important roles in mammary gland development and breast cancer. SirT1 is a highly conserved protein deacetylase that can regulate the insulin/IGF-1 signaling in lower organisms, as well as a growing number of transcription factors, including NF-kappaB, in mammalian cells. Whether SirT1 regulates the IGF-1 signaling for mammary gland development and function, however, is not clear. In the present study, this role of SirT1 was examined by studying SirT1-deficient mice.

The gut microbiome in humans is associated with geography, diet, lifestyles and so on, but its relationship with some isolated populations is not clear. We used the 16sRNA technique to sequence the fecal microbiome in the Chinese isolated Yao population and compared it with the major minority Zhuang and the major ethnic Han populations living in the same rural area. Information about diet frequency and health status and routine serum measurements were collected. The unweighted UniFrac principal coordinates analysis showed significant structural differences in fecal microbiota among the three ethnic groups. Statistically significant differences were observed in the community richness estimator (chaos) and the diversity estimator (Shannon) among the three groups. At the genus level, the fecal samples of the isolated Yao population presented the lowest relative abundance of the Megamonas genus, which was potentially related to the high frequency of bean consumption in the diet. Two enterotypes were identified in the overall fecal microbiota in the three populations. In the isolated Yao population, a higher Bacteroides abundance was observed, but the Prevotella abundance decreased with increased alcohol consumption.

The association of HLA mismatching with kidney allograft survival has been well established. We examined whether amino acid (AA) mismatches (MMs) at the antigen recognition site of HLA molecules represent independent and incremental risk factors for kidney graft failure (GF) beyond those MMs assessed at the antigenic (2-digit) specificity.

Endobronchial stents are increasingly used to treat airway complications in multiple conditions including lung transplantation but little is known about the biofilms that form on these devices. We applied deep sequencing to profile luminal biofilms of 46 endobronchial stents removed from 20 subjects primarily with lung transplantation-associated airway compromise. Microbial communities were analyzed by bacterial 16S rRNA and fungal ITS marker gene sequencing. Corynebacterium was the most common bacterial taxa across biofilm communities. Clustering analysis revealed three bacterial biofilm types: one low diversity and dominated by Corynebacterium; another was polymicrobial and characterized by Staphylococcus; and the third was polymicrobial and associated with Pseudomonas, Streptococcus, and Prevotella. Biofilm type was significantly correlated with stent material: covered metal with the Staphylococcus-type biofilm, silicone with the Corynebacterium-dominated biofilm, and uncovered metal with the polymicrobial biofilm. Subjects with sequential stents had frequent transitions between community types. Fungal analysis found Candida was most prevalent, Aspergillus was common and highly enriched in two of three stents associated with airway anastomotic dehiscence, and fungal taxa not typically considered pathogens were highly enriched in some stents. Thus, molecular analysis revealed a complex and dynamic endobronchial stent biofilm with three bacterial types that associate with stent material, a central role for Corynebacterium, and that both expected and unexpected fungi inhabit this unique niche. The current work provides a foundation for studies to investigate the relationship between stent biofilm composition and clinical outcomes, mechanisms of biofilm establishment, and strategies for improved stent technology and use in airway compromise.

To gain insight into the aetiology of posterior subcapsular congenital cataract from the perspective of transcriptional changes, we conducted an mRNA sequencing analysis of the lenses in posterior subcapsular congenital cataract patients and in normal children. There were 1533 differentially expressed genes from 19,072 genes in the lens epithelial cells of the posterior subcapsular congenital cataract patients compared to in the normal controls at a cut-off criteria of |log2 fold change| of >1 and a p-value of <0.05, including 847 downregulated genes and 686 upregulated genes. To further narrow down the DEGs, we utilised the stricter criteria of |log2 fold change| of >1 and an FDR value of <0.05, and we identified 551 DEGs, including 97 upregulated genes and 454 downregulated genes. This study also identified 1263 differentially expressed genes of the 18,755 genes in lens cortex and nuclear fibres, including 646 downregulated genes and 617 upregulated genes. The downregulated genes in epithelial cells were significantly enriched in the structural constituent of lenses, lens development and lens fibre cell differentiation. After filtering the DEGs using the databases iSyTE and Cat-Map, several high-priority candidate genes related to posterior subcapsular congenital cataract such as GRIFIN, HTRA1 and DAPL1 were identified. The findings of our study may provide a deeper understanding of the mechanisms of posterior subcapsular congenital cataract and help in the prevention and treatment of this disease.

Congenital cataract is the leading cause of blindness among children worldwide. Patients with posterior subcapsular congenital cataract (PSC) in the central visual axis can result in worsening vision and stimulus deprivation amblyopia. However, the pathogenesis of PSC remains unclear. This study aims to explore the functional regulation and mechanism of HTRA1 in human lens epithelial cells (HLECs). HTRA1 was significantly downregulated in the lens capsules of children with PSC compared to normal controls. HTRA1 is a suppression factor of transforming growth factor-β (TGF-β) signalling pathway, which plays a key role in cataract formation. The results showed that the TGF-β/Smad signalling pathway was activated in the lens tissue of PSC. The effect of HTRA1 on cell proliferation, migration and apoptosis was measured in HLECs. In primary HLECs, the downregulation of HTRA1 can promote the proliferation and migration of HLECs by activating the TGF-β/Smad signalling pathway and can significantly upregulate the TGF-β/Smad downstream target genes FN1 and α-SMA. HTRA1 was also knocked out in the eyes of C57BL/6J mice via adeno-associated virus-mediated RNA interference. The results showed that HTRA1 knockout can significantly upregulate p-Smad2/3 and activate the TGF-β/Smad signalling pathway, resulting in abnormal proliferation and irregular arrangement of lens epithelial cells and leading to the occurrence of subcapsular cataract. To conclude, HTRA1 was significantly downregulated in children with PSC, and the downregulation of HTRA1 enhanced the proliferation and migration of HLECs by activating the TGF-β/Smad signalling pathway, which led to the occurrence of PSC.

Epigenetic changes might provide the biological explanation for the long-lasting impact of metabolic alterations of diabetic kidney disease development. Here we examined cytosine methylation of human kidney tubules using Illumina Infinium 450 K arrays from 91 subjects with and without diabetes and varying degrees of kidney disease using a cross-sectional design. We identify cytosine methylation changes associated with kidney structural damage and build a model for kidney function decline. We find that the methylation levels of 65 probes are associated with the degree of kidney fibrosis at genome wide significance. In total 471 probes improve the model for kidney function decline. Methylation probes associated with kidney damage and functional decline enrich on kidney regulatory regions and associate with gene expression changes, including epidermal growth factor (EGF). Altogether, our work shows that kidney methylation differences can be detected in patients with diabetic kidney disease and improve kidney function decline models indicating that they are potentially functionally important.

Initial microbial colonization and later succession in the gut of human infants are linked to health and disease later in life. The timing of the appearance of the first gut microbiome, and the consequences for the early life metabolome, are just starting to be defined. Here, we evaluated the gut microbiome, proteome and metabolome in 88 African-American newborns using faecal samples collected in the first few days of life. Gut bacteria became detectable using molecular methods by 16 h after birth. Detailed analysis of the three most common species, Escherichia coli, Enterococcus faecalis and Bacteroides vulgatus, did not suggest a genomic signature for neonatal gut colonization. The appearance of bacteria was associated with reduced abundance of approximately 50 human proteins, decreased levels of free amino acids and an increase in products of bacterial fermentation, including acetate and succinate. Using flux balance modelling and in vitro experiments, we provide evidence that fermentation of amino acids provides a mechanism for the initial growth of E. coli, the most common early colonizer, under anaerobic conditions. These results provide a deep characterization of the first microbes in the human gut and show how the biochemical environment is altered by their appearance.

Major depressive disorder (MDD) and anxiety were recognized in treating pancreatic ductal adenocarcinoma (PDAC). This longitudinal study identified risk factors for MDD and anxiety and established associations with patients' quality of life (QoL) and survival outcomes.

Epigenetic changes may fill a critical gap in our understanding of kidney disease development, as they not only reflect metabolic changes but are also preserved and transmitted during cell division. We conducted a genome-wide cytosine methylation analysis of 399 human kidney samples, along with single-nuclear open chromatin analysis on over 60,000 cells from 14 subjects, including controls, and diabetes and hypertension attributed chronic kidney disease (CKD) patients. We identified and validated differentially methylated positions associated with disease states, and discovered that nearly 30% of these alterations were influenced by underlying genetic variations, including variants known to be associated with kidney disease in genome-wide association studies. We also identified regions showing both methylation and open chromatin changes. These changes in methylation and open chromatin significantly associated gene expression changes, most notably those playing role in metabolism and expressed in proximal tubules. Our study further demonstrated that methylation risk scores (MRS) can improve disease state annotation and prediction of kidney disease development. Collectively, our results suggest a causal relationship between epigenetic changes and kidney disease pathogenesis, thereby providing potential pathways for the development of novel risk stratification methods.

In addition to cochleotoxicity, systemic aminoglycoside pharmacotherapy causes vestibulotoxicity resulting in imbalance and visual dysfunction. The underlying trafficking routes of systemically-administered aminoglycosides from the vasculature to the vestibular sensory hair cells are largely unknown. We investigated the trafficking of systemically-administered gentamicin into the peripheral vestibular system in C56Bl/6 mice using fluorescence-tagged gentamicin (gentamicin-Texas-Red, GTTR) imaged by scanning laser confocal microscopy to determine the cellular distribution and intensity of GTTR fluorescence in the three semicircular canal cristae, utricular, and saccular maculae at 5 time points over 4 hours.

Human bone marrow (BM) contains a rare population of nonhematopoietic mesenchymal stromal cells (MSCs), which are of central importance for the hematopoietic microenvironment. However, the precise phenotypic definition of these cells in adult BM has not yet been reported. In this study, we show that low/negative expression of CD140a (PDGFR-α) on lin(-)/CD45(-)/CD271(+) BM cells identified a cell population with very high MSC activity, measured as fibroblastic colony-forming unit frequency and typical in vitro and in vivo stroma formation and differentiation capacities. Furthermore, these cells exhibited high levels of genes associated with mesenchymal lineages and HSC supportive function. Moreover, lin(-)/CD45(-)/CD271(+)/CD140a(low/-) cells effectively mediated the ex vivo expansion of transplantable CD34(+) hematopoietic stem cells. Taken together, these data indicate that CD140a is a key negative selection marker for adult human BM-MSCs, which enables to prospectively isolate a close to pure population of candidate human adult stroma stem/progenitor cells with potent hematopoiesis-supporting capacity.

Mesenchymal stem cells (MSC) have not only been implicated in the development of lung diseases, but they have also been proposed as a future cell-based therapy for lung diseases. However, the cellular identity of the primary MSC in human lung tissues has not yet been reported. This study therefore aimed to identify and characterise the 'bona fide' MSC in human lungs and to investigate if the MSC numbers correlate with the development of bronchiolitis obliterans syndrome in lung-transplanted patients.

Enterovirus 71 (EV71) and Coxsackievirus A16 (CA16) remain the predominant pathogens in hand, foot, and mouth disease (HFMD), but the factors underlying the pathogenesis of EV71 and CA16 infections have not been elucidated. Recently, the functions of microRNAs (miRNAs) in pathogen-host interactions have been highlighted. In the present study, we performed comprehensive miRNA profiling in EV71- and CA16-infected human umbilical vein endothelial cells (HUVECs) at multiple time points using high-throughput sequencing. The results showed that 135 known miRNAs exhibited remarkable differences in expression. Of these, 30 differentially expressed miRNAs presented opposite trends in EV71- and CA16-infected samples. Subsequently, we mainly focused on the 30 key differentially expressed miRNAs through further screening to predict targets. Gene ontology (GO) and pathway analysis of the predicted targets showed the enrichment of 14 biological processes, 9 molecular functions, 8 cellular components, and 85 pathways. The regulatory networks of these miRNAs with predicted targets, GOs, pathways, and co-expression genes were determined, suggesting that miRNAs display intricate regulatory mechanisms during the infection phase. Consequently, we specifically analyzed the hierarchical GO categories of the predicted targets involved in biological adhesion. The results indicated that the distinct changes induced by EV71 and CA16 infection may be partly linked to the function of the blood-brain barrier. Taken together, this is the first report describing miRNA expression profiles in HUVECs with EV71 and CA16 infections using high-throughput sequencing. Our data provide useful insights that may help to elucidate the different host-pathogen interactions following EV71 and CA16 infection and offer novel therapeutic targets for these infections.

Bone marrow mesenchymal stromal cells (BM-MSCs) are a rare population of cells that gives rise to skeletal tissues and the hematopoietic stroma in vivo. Recently, we have demonstrated that BM-MSCs fulfill stringent in vivo stem cell criteria when propagated as non-adherent mesenspheres but not as adherent-cultured cells. Motivated by these profound functional differences, the current study aimed to identify potential important MSC regulators by investigating global gene expression profiles of adherent and non-adherent culture-derived BM-MSCs in comparison with primary BM-MSCs. A substantial number of genes were differentially expressed between primary and culture-expanded cells already early upon culture, and numerous genes were found to be different when comparing adherent and non-adherent BM-MSCs. Cluster analysis identified 16 sets of genes of which two displayed comparable gene expression levels in primary and non-adherent cultured cells, but not in adherent cultured cells. This pattern suggested that these clusters contained candidate regulators of BM-MSCs. Gene expression differences were confirmed for selected genes and BM-MSC transcription factors by protein analysis and RT-PCR, respectively. Taken together, these data demonstrated profound gene expression changes upon culture of primary BM-MSCs. Moreover, gene cluster differences provide the basis to uncover the regulatory mechanisms that control primary and cultured BM-MSCs.