Single-Cell RNA Sequencing of Donor-Reactive T Cells Reveals Role of Apoptosis in Donor-Specific Hyporesponsiveness of Kidney Transplant Recipients

After kidney transplantation (KT), donor-specific hyporesponsiveness (DSH) of recipient T cells develops over time. Recently, apoptosis was identified as a possible underlying mechanism. In this study, both transcriptomic profiles and complete V(D)J variable regions of TR transcripts from individual alloreactive T cells of kidney transplant recipients were determined with single-cell RNA sequencing. Alloreactive T cells were identified by CD137 expression after stimulation of peripheral blood mononuclear cells obtained from KT recipients (N = 7) prior to and 3-5 years after transplantation with cells of their donor or a third party control. The alloreactive T cells were sorted, sequenced and the transcriptome and T cell receptor profiles were analyzed using unsupervised clustering. Alloreactive T cells retain a highly polyclonal T Cell Receptor Alpha/Beta repertoire over time. Post transplantation, donor-reactive CD4+ T cells had a specific downregulation of genes involved in T cell cytokine-mediated pathways and apoptosis. The CD8+ donor-reactive T cell profile did not change significantly over time. Single-cell expression profiling shows that activated and pro-apoptotic donor-reactive CD4+ T cell clones are preferentially lost after transplantation in stable kidney transplant recipients.

PDGFRβ+ cells play a dual role as hematopoietic precursors and niche cells during mouse ontogeny

Hematopoietic stem cell (HSC) generation in the aorta-gonad-mesonephros region requires HSC specification signals from the surrounding microenvironment. In zebrafish, PDGF-B/PDGFRβ signaling controls hematopoietic stem/progenitor cell (HSPC) generation and is required in the HSC specification niche. Little is known about murine HSPC specification in vivo and whether PDGF-B/PDGFRβ is involved. Here, we show that PDGFRβ is expressed in distinct perivascular stromal cell layers surrounding the mid-gestation dorsal aorta, and its deletion impairs hematopoiesis. We demonstrate that PDGFRβ⁺ cells play a dual role in murine hematopoiesis. They act in the aortic niche to support HSPCs, and in addition, PDGFRβ⁺ embryonic precursors give rise to a subset of HSPCs that persist into adulthood. These findings provide crucial information for the controlled production of HSPCs in vitro.

•

PDGFRβ deletion affects hematopoietic development in the AGM in vivo

•

The transcriptome and hematopoietic support of the PDGFRβ-KO niche are altered

•

The osteogenic gene profile and differentiation of KO AGM MSCs are affected

•

PDGFRβ⁺ early embryonic precursors contribute to EC and HSPC lineages in vivo

JMJD3 intrinsically disordered region links the 3D-genome structure to TGFβ-dependent transcription activation

Enhancers are key regulatory elements that govern gene expression programs in response to developmental signals. However, how multiple enhancers arrange in the 3D-space to control the activation of a specific promoter remains unclear. To address this question, we exploited our previously characterized TGFβ-response model, the neural stem cells, focusing on a ~374 kb locus where enhancers abound. Our 4C-seq experiments reveal that the TGFβ pathway drives the assembly of an enhancer-cluster and precise gene activation. We discover that the TGFβ pathway coactivator JMJD3 is essential to maintain these structures. Using live-cell imaging techniques, we demonstrate that an intrinsically disordered region contained in JMJD3 is involved in the formation of phase-separated biomolecular condensates, which are found in the enhancer-cluster. Overall, in this work we uncover novel functions for the coactivator JMJD3, and we shed light on the relationships between the 3D-conformation of the chromatin and the TGFβ-driven response during mammalian neurogenesis.

Here the authors demonstrate that TGFβ drives multi-enhancer contacts and ultimately gene activation during neuronal commitment, and that this requires the intrinsically disordered region (IDR) of the histone demethylase JMJD3 likely through its role in promoting phase-separated biomolecular condensates.

The Somatic Mutation Paradigm in Congenital Malformations: Hirschsprung Disease as a Model

Patients with Hirschsprung disease (HSCR) do not always receive a genetic diagnosis after routine screening in clinical practice. One of the reasons for this could be that the causal mutation is not present in the cell types that are usually tested—whole blood, dermal fibroblasts or saliva—but is only in the affected tissue. Such mutations are called somatic, and can occur in a given cell at any stage of development after conception. They will then be present in all subsequent daughter cells. Here, we investigated the presence of somatic mutations in HSCR patients. For this, whole-exome sequencing and copy number analysis were performed in DNA isolated from purified enteric neural crest cells (ENCCs) and blood or fibroblasts of the same patient. Variants identified were subsequently validated by Sanger sequencing. Several somatic variants were identified in all patients, but causative mutations for HSCR were not specifically identified in the ENCCs of these patients. Larger copy number variants were also not found to be specific to ENCCs. Therefore, we believe that somatic mutations are unlikely to be identified, if causative for HSCR. Here, we postulate various modes of development following the occurrence of a somatic mutation, to describe the challenges in detecting such mutations, and hypothesize how somatic mutations may contribute to ‘missing heritability’ in developmental defects.

RMplex: An efficient method for analyzing 30 Y-STRs with high mutation rates

Y-chromosomal short tandem repeats (Y-STRs) with high mutation rates are recognized as valuable genetic markers for differentiating paternally related men, who typically cannot be separated with standard Y-STRs, and were shown to provide paternal lineage differentiation on a higher resolution level than standard Y-STRs. Both features make Y-STRs with high mutation rates relevant in criminal casework, particularly in sexual assault cases involving highly unbalanced male-female DNA mixtures that often fail autosomal forensic STR profiling for the male donor. Previously, the number of known Y-STRs with mutation rates higher than 10^-2 per locus per generation termed rapidly mutating Y-STRs (RM Y-STRs) was limited to 13, which has recently been overcome by the discovery and characterization of 12 additional RM Y-STRs. Here, we present the development and validation of RMplex, an efficient genotyping system for analyzing 30 Y-STRs with high mutation rates, including all currently known RM Y-STRs, using multiplex PCR with capillary electrophoresis (CE) or massively parallel sequencing (MPS), overall targeting a total of 44 male-specific loci. If previously unavailable, repeat number assignations were provided based on newly generated MPS data. Validation tests based on the CE method demonstrated that the results were both repeatable and reproducible, full profiles were achieved with minimal input DNA of 250 pg for RMplex 1 and 100 pg for RMplex 2, and in the presence of inhibitors, or with a surplus of female DNA, the assays performed reasonably well. Application of RMplex to differentiate between paternally related men was exemplified in 32 males belonging to five different paternal pedigrees. Given further successful forensic validation testing, we envision the future application of RMplex in criminal cases where it is suspected, or cannot be excluded, that the crime scene trace originated from a male relatives of the suspect who is highlighted with standard Y-STR matching. Other applications of RMplex are in criminal cases without known suspects to differentiate between male relatives highlighted in familial searching based on standard Y-STR matching.

Unraveling the Genetics of Congenital Diaphragmatic Hernia: An Ongoing Challenge

Congenital diaphragmatic hernia (CDH) is a congenital structural anomaly in which the diaphragm has not developed properly. It may occur either as an isolated anomaly or with additional anomalies. It is thought to be a multifactorial disease in which genetic factors could either substantially contribute to or directly result in the developmental defect. Patients with aneuploidies, pathogenic variants or de novo Copy Number Variations (CNVs) impacting specific genes and loci develop CDH typically in the form of a monogenetic syndrome. These patients often have other associated anatomical malformations. In patients without a known monogenetic syndrome, an increased genetic burden of de novo coding variants contributes to disease development. In early years, genetic evaluation was based on karyotyping and SNP-array. Today, genomes are commonly analyzed with next generation sequencing (NGS) based approaches. While more potential pathogenic variants are being detected, analysis of the data presents a bottleneck-largely due to the lack of full appreciation of the functional consequence and/or relevance of the detected variant. The exact heritability of CDH is still unknown. Damaging de novo alterations are associated with the more severe and complex phenotypes and worse clinical outcome. Phenotypic, genetic-and likely mechanistic-variability hampers individual patient diagnosis, short and long-term morbidity prediction and subsequent care strategies. Detailed phenotyping, clinical follow-up at regular intervals and detailed registries are needed to find associations between long-term morbidity, genetic alterations, and clinical parameters. Since CDH is a relatively rare disorder with only a few recurrent changes large cohorts of patients are needed to identify genetic associations. Retrospective whole genome sequencing of historical patient cohorts using will yield valuable data from which today's patients and parents will profit Trio whole genome sequencing has an excellent potential for future re-analysis and data-sharing increasing the chance to provide a genetic diagnosis and predict clinical prognosis. In this review, we explore the pitfalls and challenges in the analysis and interpretation of genetic information, present what is currently known and what still needs further study, and propose strategies to reap the benefits of genetic screening.

Pathogen-induced activation of disease-suppressive functions in the endophytic root microbiome

Microorganisms living inside plants can promote plant growth and health, but their genomic and functional diversity remain largely elusive. Here, metagenomics and network inference show that fungal infection of plant roots enriched for Chitinophagaceae and Flavobacteriaceae in the root endosphere and for chitinase genes and various unknown biosynthetic gene clusters encoding the production of nonribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs). After strain-level genome reconstruction, a consortium of Chitinophaga and Flavobacterium was designed that consistently suppressed fungal root disease. Site-directed mutagenesis then revealed that a previously unidentified NRPS-PKS gene cluster from Flavobacterium was essential for disease suppression by the endophytic consortium. Our results highlight that endophytic root microbiomes harbor a wealth of as yet unknown functional traits that, in concert, can protect the plant inside out.

G6PD genetic variations in neonatal Hyperbilirubinemia in Indonesian Deutromalay population

BACKGROUND

Neonatal jaundice is a common finding in newborns in Asia, including Indonesia. In some cases, the serum total bilirubin levels exceeds the 95th percentile for hours of life (neonatal hyperbilirubinemia). Severe neonatal hyperbilirubinemia (NH) could lead to kernicterus and neonatal death. Glucose-6-Phosphage Dehydrogenase (G6PD) genetic variations and deficiency have been reported in several studies to be associated with NH. This study aimed to analyze the G6PD genetic variations and its activity in neonates with and without hyperbilirubinemia in the Deutromalay Indonesian population.

METHODS

Deoxyribose Nucleic Acid (DNA) was isolated from peripheral blood of 116 and 115 healthy term neonates with and without hyperbilirubinemia. All infants underwent the following laboratory examinations: routine hematologic evaluation, Coombs test, G6PD activity measurement using the Randox kit method, and serum total bilirubin level. All exons of the G6PD gene were targeted for deep sequencing using MiSeq (Illumina). An association study of G6PD polymorphisms with NH was performed using PLINK.

RESULTS

The prevalence of G6PD deficiency in neonates with and without hyperbilirubinemia in Indonesian Deutromalay population were 1.72% (95% Confidence Interval (CI): 0.6-4.1%) and 1.74% (95% CI: 0.7-4.1%), respectively. The most common G6PD polymorphisms, i.e. rs1050757/c.* + 357A > G, rs2230037/c.1311C > T, and rs2071429/c.1365-13 T/IVS11, were identified. However, none of those polymorphisms and their haplotype were associated with NH (p > 0.05, Odds Ratio (OR) ~1.00). The prevalence of G6PD mutations in neonates with and without hyperbilirubinemia were 6.8% (95% CI: 2.3-11.5%) and 6.9% (95% CI: 2.3-11.6%), respectively. The most frequently identified G6PD mutation was the Viangchan variant (p.V291 M), which was followed by the Canton (p.R459L) and Vanua Lava (p.L128P) variants. Two novel mutations were identified both in case (p.V369A, p.I167F) and control (p.L474=, p.I36T) groups.

CONCLUSION

The prevalence of G6PD deficiency is low in neonates with or without hyperbilirubinemia in Deutromalay Indonesian population. The majority of G6PD mutations identified among Indonesian Deutromalay population in this study are Viangchan, Canton and Vanua Lava variants.

Colony-Stimulating Factor 1 Receptor (CSF1R) Regulates Microglia Density and Distribution, but Not Microglia Differentiation In Vivo

Microglia are brain-resident macrophages with trophic and phagocytic functions. Dominant loss-of-function mutations in a key microglia regulator, colony-stimulating factor 1 receptor (CSF1R), cause adult-onset leukoencephalopathy with axonal spheroids and pigmented glia (ALSP), a progressive white matter disorder. Because it remains unclear precisely how CSF1R mutations affect microglia, we generated an allelic series of csf1r mutants in zebrafish to identify csf1r-dependent microglia changes. We found that csf1r mutations led to aberrant microglia density and distribution and regional loss of microglia. The remaining microglia still had a microglia-specific gene expression signature, indicating that they had differentiated normally. Strikingly, we also observed lower microglia numbers and widespread microglia depletion in postmortem brain tissue of ALSP patients. Both in zebrafish and in human disease, local microglia loss also presented in regions without obvious pathology. Together, this implies that CSF1R mainly regulates microglia density and that early loss of microglia may contribute to ALSP pathogenesis.

Modeling the Pathological Long-Range Regulatory Effects of Human Structural Variation with Patient-Specific hiPSCs

The pathological consequences of structural variants disrupting 3D genome organization can be difficult to elucidate in vivo due to differences in gene dosage sensitivity between mice and humans. This is illustrated by branchiooculofacial syndrome (BOFS), a rare congenital disorder caused by heterozygous mutations within TFAP2A, a neural crest regulator for which humans, but not mice, are haploinsufficient. Here, we present a BOFS patient carrying a heterozygous inversion with one breakpoint located within a topologically associating domain (TAD) containing enhancers essential for TFAP2A expression in human neural crest cells (hNCCs). Using patient-specific hiPSCs, we show that, although the inversion shuffles the TFAP2A hNCC enhancers with novel genes within the same TAD, this does not result in enhancer adoption. Instead, the inversion disconnects one TFAP2A allele from its cognate enhancers, leading to monoallelic and haploinsufficient TFAP2A expression in patient hNCCs. Our work illustrates the power of hiPSC differentiation to unveil long-range pathomechanisms.

Whole-Genome Linkage Scan Combined With Exome Sequencing Identifies Novel Candidate Genes for Carotid Intima-Media Thickness

Carotid intima-media thickness (cIMT) is an established heritable marker for subclinical atherosclerosis. In this study, we aim to identify rare variants with large effects driving differences in cIMT by performing genome-wide linkage analysis of individuals in the extremes of cIMT trait distribution (>90th percentile) in a large family-based study from a genetically isolated population in the Netherlands. Linked regions were subsequently explored by fine-mapping using exome sequencing. We observed significant evidence of linkage on chromosomes 2p16.3 [rs1017418, heterogeneity LOD (HLOD) = 3.35], 19q13.43 (rs3499, HLOD = 9.09), 20p13 (rs1434789, HLOD = 4.10), and 21q22.12 (rs2834949, HLOD = 3.59). Fine-mapping using exome sequencing data identified a non-coding variant (rs62165235) in PNPT1 gene under the linkage peak at chromosome 2 that is likely to have a regulatory function. The variant was associated with quantitative cIMT in the family-based study population (effect = 0.27, p-value = 0.013). Furthermore, we identified several genes under the linkage peak at chromosome 21 highly expressed in tissues relevant for atherosclerosis. To conclude, our linkage analysis identified four genomic regions significantly linked to cIMT. Further analyses are needed to demonstrate involvement of identified candidate genes in development of atherosclerosis.

DOC1-Dependent Recruitment of NURD Reveals Antagonism with SWI/SNF during Epithelial-Mesenchymal Transition in Oral Cancer Cells

The Nucleosome Remodeling and Deacetylase (NURD) complex is a key regulator of cell differentiation that has also been implicated in tumorigenesis. Loss of the NURD subunit Deleted in Oral Cancer 1 (DOC1) is associated with human oral squamous cell carcinomas (OSCCs). Here, we show that restoration of DOC1 expression in OSCC cells leads to a reversal of epithelial-mesenchymal transition (EMT). This is caused by the DOC1-dependent targeting of NURD to repress key transcriptional regulators of EMT. NURD recruitment drives extensive epigenetic reprogramming, including eviction of the SWI/SNF remodeler, formation of inaccessible chromatin, H3K27 deacetylation, and binding of PRC2 and KDM1A, followed by H3K27 methylation and H3K4 demethylation. Strikingly, depletion of SWI/SNF mimics the effects of DOC1 re-expression. Our results suggest that SWI/SNF and NURD function antagonistically to control chromatin state and transcription. We propose that disturbance of this dynamic equilibrium may lead to defects in gene expression that promote oncogenesis.

Regulation of the cohesin-loading factor NIPBL: Role of the lncRNA NIPBL-AS1 and identification of a distal enhancer element

Cohesin is crucial for genome stability, cell division, transcription and chromatin organization. Its functions critically depend on NIPBL, the cohesin-loader protein that is found to be mutated in >60% of the cases of Cornelia de Lange syndrome (CdLS). Other mutations are described in the cohesin subunits SMC1A, RAD21, SMC3 and the HDAC8 protein. In 25-30% of CdLS cases no mutation in the known CdLS genes is detected. Until now, functional elements in the noncoding genome were not characterized in the molecular etiology of CdLS and therefore are excluded from mutation screening, although the impact of such mutations has now been recognized for a wide range of diseases. We have identified different elements of the noncoding genome involved in regulation of the NIPBL gene. NIPBL-AS1 is a long non-coding RNA transcribed upstream and antisense to NIPBL. By knockdown and transcription blocking experiments, we could show that not the NIPBL-AS1 gene product, but its actual transcription is important to regulate NIPBL expression levels. This reveals a possibility to boost the transcriptional activity of the NIPBL gene by interfering with the NIPBL-AS1 lncRNA. Further, we have identified a novel distal enhancer regulating both NIPBL and NIPBL-AS1. Deletion of the enhancer using CRISPR genome editing in HEK293T cells reduces expression of NIPBL, NIPBL-AS1 as well as genes found to be dysregulated in CdLS.

Loss-of-Function Variants in MYLK Cause Recessive Megacystis Microcolon Intestinal Hypoperistalsis Syndrome

Megacystis microcolon intestinal hypoperistalsis syndrome (MMIHS) is a congenital disorder characterized by loss of smooth muscle contraction in the bladder and intestine. To date, three genes are known to be involved in MMIHS pathogenesis: ACTG2, MYH11, and LMOD1. However, for approximately 10% of affected individuals, the genetic cause of the disease is unknown, suggesting that other loci are most likely involved. Here, we report on three MMIHS-affected subjects from two consanguineous families with no variants in the known MMIHS-associated genes. By performing homozygosity mapping and whole-exome sequencing, we found homozygous variants in myosin light chain kinase (MYLK) in both families. We identified a 7 bp duplication (c.3838_3844dupGAAAGCG [p.Glu1282_Glyfs∗51]) in one family and a putative splice-site variant (c.3985+5C>A) in the other. Expression studies and splicing assays indicated that both variants affect normal MYLK expression. Because MYLK encodes an important kinase required for myosin activation and subsequent interaction with actin filaments, it is likely that in its absence, contraction of smooth muscle cells is impaired. The existence of a conditional-Mylk-knockout mouse model with severe gut dysmotility and abnormal function of the bladder supports the involvement of this gene in MMIHS pathogenesis. In aggregate, our findings implicate MYLK as a gene involved in the recessive form of MMIHS, confirming that this disease of the visceral organs is heterogeneous with a myopathic origin.

Complex MAX Rearrangement in a Family With Malignant Pheochromocytoma, Renal Oncocytoma, and Erythrocytosis

CONTEXT

Familial pheochromocytoma (PCC) has been associated with germline mutations in 16 genes. Here we investigated three siblings presenting with bilateral pheochromocytomas. In addition, the index patient also exhibited renal oncocytoma and erythrocytosis, whereas the second sibling presented with a lymph node metastasis.

DESIGN

First, single-nucleotide polymorphism array and exome sequencing were performed on germline and PCC-derived DNA to identify genomic alterations in the index patient. Second, alterations were confirmed and validated by Sanger sequencing, analyzed by (multiplexed) PCR to determine the loss of the wild-type allele, and investigated by immunohistochemistry in the tumors of the three siblings.

RESULTS

The index patient's germline DNA revealed a large complex genomic alteration encompassing the intragenic and promoter regions of Myc-associated factor X (MAX) and alpha-(1,6)-fucosyltransferase (FUT8). In all three siblings the MAX alteration was confirmed, and the loss of the wild-type MAX and FUT8 alleles was demonstrated in all tumors. Uniparental disomy of chromosome 14q, previously demonstrated as a hallmark for MAX-related PCC, was shown in the index patient's PCC by single-nucleotide polymorphism array. Loss of MAX and FUT8 protein expression was demonstrated by immunohistochemistry in the tumors from the three siblings.

CONCLUSIONS

Our results indicate that large genomic deletions of MAX should be considered in familial and bilateral PCC with prior negative testing for gene mutations. In addition, our results confirm that MAX is a tumor suppressor gene for renal oncocytomas.

In vivo murine hepatic microRNA and mRNA expression signatures predicting the (non-)genotoxic carcinogenic potential of chemicals

There is a high need to improve the assessment of, especially non-genotoxic, carcinogenic features of chemicals. We therefore explored a toxicogenomics-based approach using genome-wide microRNA and mRNA expression profiles upon short-term exposure in mice. For this, wild-type mice were exposed for seven days to three different classes of chemicals, i.e., four genotoxic carcinogens (GTXC), seven non-genotoxic carcinogens (NGTXC), and five toxic non-carcinogens. Hepatic expression patterns of mRNA and microRNA transcripts were determined after exposure and used to assess the discriminative power of the in vivo transcriptome for GTXC and NGTXC. A final classifier set, discriminative for GTXC and NGTXC, was generated from the transcriptomic data using a tiered approach. This appeared to be a valid approach, since the predictive power of the final classifier set in three different classifier algorithms was very high for the original training set of chemicals. Subsequent validation in an additional set of chemicals revealed that the predictive power for GTXC remained high, in contrast to NGTXC, which appeared to be more troublesome. Our study demonstrated that the in vivo microRNA-ome has less discriminative power to correctly identify (non-)genotoxic carcinogen classes. The results generally indicate that single mRNA transcripts do have the potential to be applied in risk assessment, but that additional (genomic) strategies are necessary to correctly predict the non-genotoxic carcinogenic potential of a chemical.

DNA damage responsive microRNAs misexpressed in human cancer modulate therapy sensitivity

The DNA damage response (DDR) is activated upon DNA damage and prevents accumulation of mutations and chromosomal rearrangements, both driving carcinogenesis. Tumor cells often have defects in the DDR, which in combination with continuous cell proliferation are exploited by genotoxic cancer therapies. Most cancers, overcome initial sensitivity and develop drug resistance, e.g. by modulation of the DDR. Not much is known, however, about DNA damage responsive microRNAs in cancer therapy resistance. Therefore, we mapped temporal microRNA expression changes in primary breast epithelial cells upon low and high dose exposure to the DNA damaging agents ionizing radiation and cisplatin. A third of all DDR microRNAs commonly regulated across all treatments was also misexpressed in breast cancer, indicating a DDR defect. We repeated this approach in primary lung epithelial cells and non-small cell lung cancer samples and found that more than 40% of all DDR microRNAs was deregulated in non-small cell lung cancer. Strikingly, the microRNA response upon genotoxic stress in primary breast and lung epithelial cells was markedly different, although the biological outcome of DNA damage signaling (cell death/senescence or survival) was similar. Several DDR microRNAs deregulated in cancer modulated sensitivity to anti-cancer agents. In addition we were able to distinguish between microRNAs that induced resistance by potentially inducing quiescence (miR-296-5p and miR-382) or enhancing DNA repair or increased DNA damage tolerance (miR-21). In conclusion, we provide evidence that DNA damage responsive microRNAs are frequently misexpressed in human cancer and can modulate chemotherapy sensitivity.

MicroRNA-Mediated Down-Regulation of M-CSF Receptor Contributes to Maturation of Mouse Monocyte-Derived Dendritic Cells

Dendritic cell (DC) maturation is a tightly regulated process that requires coordinated and timed developmental cues. Here we investigate whether microRNAs are involved in this process. We identify microRNAs in mouse GM-CSF-generated, monocyte-related DC (GM-DC) that are differentially expressed during both spontaneous and LPS-induced maturation and characterize M-CSF receptor (M-CSFR), encoded by the Csf1r gene, as a key target for microRNA-mediated regulation in the final step toward mature DC. MicroRNA-22, -34a, and -155 are up-regulated in mature MHCII^hi CD86^hi DC and mediate Csf1r mRNA and protein down-regulation. Experimental inhibition of Csf1r-targeting microRNAs in vitro results not only in sustained high level M-CSFR protein expression but also in impaired DC maturation upon stimulation by LPS. Accordingly, over-expression of Csf1r in GM-DC inhibits terminal differentiation. Taken together, these results show that developmentally regulated microRNAs control Csf1r expression, supplementing previously identified mechanisms that regulate its transcription and protein surface expression. Furthermore, our data indicate a novel function for Csf1r in mouse monocyte-derived DC, showing that down-regulation of M-CSFR expression is essential for final DC maturation.

Histone chaperone NAP1 mediates sister chromatid resolution by counteracting protein phosphatase 2A

Chromosome duplication and transmission into daughter cells requires the precisely orchestrated binding and release of cohesin. We found that the Drosophila histone chaperone NAP1 is required for cohesin release and sister chromatid resolution during mitosis. Genome-wide surveys revealed that NAP1 and cohesin co-localize at multiple genomic loci. Proteomic and biochemical analysis established that NAP1 associates with the full cohesin complex, but it also forms a separate complex with the cohesin subunit stromalin (SA). NAP1 binding to cohesin is cell-cycle regulated and increases during G2/M phase. This causes the dissociation of protein phosphatase 2A (PP2A) from cohesin, increased phosphorylation of SA and cohesin removal in early mitosis. PP2A depletion led to a loss of centromeric cohesion. The distinct mitotic phenotypes caused by the loss of either PP2A or NAP1, were both rescued by their concomitant depletion. We conclude that the balanced antagonism between NAP1 and PP2A controls cohesin dissociation during mitosis.

Dynamic regulation of the transcription initiation landscape at single nucleotide resolution during vertebrate embryogenesis

Spatiotemporal control of gene expression is central to animal development. Core promoters represent a previously unanticipated regulatory level by interacting with cis-regulatory elements and transcription initiation in different physiological and developmental contexts. Here, we provide a first and comprehensive description of the core promoter repertoire and its dynamic use during the development of a vertebrate embryo. By using cap analysis of gene expression (CAGE), we mapped transcription initiation events at single nucleotide resolution across 12 stages of zebrafish development. These CAGE-based transcriptome maps reveal genome-wide rules of core promoter usage, structure, and dynamics, key to understanding the control of gene regulation during vertebrate ontogeny. They revealed the existence of multiple classes of pervasive intra- and intergenic post-transcriptionally processed RNA products and their developmental dynamics. Among these RNAs, we report splice donor site-associated intronic RNA (sRNA) to be specific to genes of the splicing machinery. For the identification of conserved features, we compared the zebrafish data sets to the first CAGE promoter map of Tetraodon and the existing human CAGE data. We show that a number of features, such as promoter type, newly discovered promoter properties such as a specialized purine-rich initiator motif, as well as sRNAs and the genes in which they are detected, are conserved in mammalian and Tetraodon CAGE-defined promoter maps. The zebrafish developmental promoterome represents a powerful resource for studying developmental gene regulation and revealing promoter features shared across vertebrates.

Contribution of rare and common variants determine complex diseases-Hirschsprung disease as a model

Finding genes for complex diseases has been the goal of many genetic studies. Most of these studies have been successful by searching for genes and mutations in rare familial cases, by screening candidate genes and by performing genome wide association studies. However, only a small fraction of the total genetic risk for these complex genetic diseases can be explained by the identified mutations and associated genetic loci. In this review we focus on Hirschsprung disease (HSCR) as an example of a complex genetic disorder. We describe the genes identified in this congenital malformation and postulate that both common 'low penetrant' variants in combination with rare or private 'high penetrant' variants determine the risk on HSCR, and likely, on other complex diseases. We also discuss how new technological advances can be used to gain further insights in the genetic background of complex diseases. Finally, we outline a few steps to develop functional assays in order to determine the involvement of these variants in disease development.

Early divergence of Th1 and Th2 transcriptomes involves a small core response and sets of transiently expressed genes

Th cells can adopt a number of different phenotypes. We performed microarray-assisted mRNA profiling on antigen-stimulated, TCR transgenic murine splenocytes that were cultured in the presence of cytokines. Transcriptome snapshots of Th cells differentiating into Th1 and Th2 phenotypes were obtained at various time points. Principal component analysis shows that time since activation and Th skewing are the largest sources of variance (i.e. the largest contributing factors) in our profiling experiments. Divergence between the Th1 and Th2 phenotypes is established early and does not increase in terms of number of differential genes from day 1 to day 4 after stimulation. Notwithstanding the lack of further divergence between the Th1 and Th2 lineages, we show that gene expression is best described by a 'turnover' rather than a 'core response' model, although we find evidence for both. We identify clusters of skewed genes associated with early persistent ('core response') and late ('turnover') Th1 and Th2 gene expression. In addition to the classical Th genes, members of the Batf transcription factor family are differentially expressed in particular helper phenotypes, suggesting an important role for this family in Th-cell phenotype differentiation.

Metabolic enzyme IMPDH is also a transcription factor regulated by cellular state

Cells need to coordinate gene expression and metabolic state. Inosine monophosphate dehydrogenase (IMPDH) controls the guanine nucleotide pool and, thereby, cell proliferation. We found that Drosophila IMPDH is also a DNA-binding transcriptional repressor. IMPDH attenuates expression of histone genes and E2f, a key driver of cell proliferation. Nuclear IMPDH accumulates during the G2 phase of the cell cycle or following replicative or oxidative stress. Thus, IMPDH can couple the expression of histones and E2F to cellular state. Genome-wide profiling and in vitro binding assays established that IMPDH binds sequence specifically to single-stranded, CT-rich DNA elements. Surprisingly, this DNA-binding function is conserved in E. coli IMPDH. The catalytic function of IMPDH is not required for DNA binding. Yet substitutions that correspond to human retinitis pigmentosa mutations disrupt IMPDH binding to CT-rich, single-stranded DNA elements. By doubling as nucleotide biosynthetic enzyme or transcription factor, IMPDH can either enable or restrict cell proliferation.

Progesterone inhibits epithelial-to-mesenchymal transition in endometrial cancer

BACKGROUND

Every year approximately 74,000 women die of endometrial cancer, mainly due to recurrent or metastatic disease. The presence of tumor infiltrating lymphocytes (TILs) as well as progesterone receptor (PR) positivity has been correlated with improved prognosis. This study describes two mechanisms by which progesterone inhibits metastatic spread of endometrial cancer: by stimulating T-cell infiltration and by inhibiting epithelial-to-mesenchymal cell transition (EMT).

METHODOLOGY AND PRINCIPAL FINDINGS

Paraffin sections from patients with (n = 9) or without (n = 9) progressive endometrial cancer (recurrent or metastatic disease) were assessed for the presence of CD4+ (helper), CD8+ (cytotoxic) and Foxp3+ (regulatory) T-lymphocytes and PR expression. Progressive disease was observed to be associated with significant loss of TILs and loss of PR expression. Frozen tumor samples, used for genome-wide expression analysis, showed significant regulation of pathways involved in immunesurveillance, EMT and metastasis. For a number of genes, such as CXCL14, DKK1, DKK4, PEG10 and WIF1, quantitive RT-PCR was performed to verify up- or downregulation in progressive disease. To corroborate the role of progesterone in regulating invasion, Ishikawa (IK) endometrial cancer cell lines stably transfected with PRA (IKPRA), PRB (IKPRB) and PRA+PRB (IKPRAB) were cultured in presence/absence of progesterone (MPA) and used for genome-wide expression analysis, Boyden- and wound healing migration assays, and IHC for known EMT markers. IKPRB and IKPRAB cell lines showed MPA induced inhibition of migration and loss of the mesenchymal marker vimentin at the invasive front of the wound healing assay. Furthermore, pathway analysis of significantly MPA regulated genes showed significant down regulation of important pathways involved in EMT, immunesuppression and metastasis: such as IL6-, TGF-β and Wnt/β-catenin signaling.

CONCLUSION

Intact progesterone signaling in non-progressive endometrial cancer seems to be an important factor stimulating immunosurveilance and inhibiting transition from an epithelial to a more mesenchymal, more invasive phenotype.

Transcriptional profiling of fibroblasts from patients with mutations in MCT8 and comparative analysis with the human brain transcriptome

Thyroid hormone (TH) is crucial for normal brain development. TH transporters control TH homeostasis in brain as evidenced by the complex endocrine and neurological phenotype of patients with mutations in monocarboxylate transporter 8 (MCT8). We investigated the mechanisms of disease by analyzing gene expression profiles in fibroblasts from patients with MCT8 mutations. Studying MCT8 and its transcriptional context in different comprehensive spatial and temporal human brain transcriptome data sets revealed distinct region-specific MCT8 expression. Furthermore, MCT8 demonstrated a clear age-dependent decrease, suggesting its importance in early brain development. Performing comparative transcriptome analysis, we linked the genes differentially expressed (DE) in patient fibroblasts to the human brain transcriptome. DE genes in patient fibroblasts were strongly over-represented among genes highly correlated with MCT8 expression in brain. Furthermore, using the same approach we identified which genes in the classical TH signaling pathway are affected in patients. Finally, we provide evidence that the TRα2 receptor variant is closely connected to MCT8. The present study provides a molecular basis for understanding which pathways are likely affected in the brains of patients with mutations in MCT8. Our data regarding a functional relationship between MCT8 and TRα2 suggest an unanticipated role for TRα2 in the (patho)physiology of TH signaling in the brain. This study demonstrates how genome-wide expression data from patient-derived non-neuronal tissue related to the human brain transcriptome may be successfully employed to improve our understanding of neurological disease.

Unacylated ghrelin rapidly modulates lipogenic and insulin signaling pathway gene expression in metabolically active tissues of GHSR deleted mice

Background

There is increasing evidence that unacylated ghrelin (UAG) improves insulin sensitivity and glucose homeostasis; however, the mechanism for this activity is not fully understood since a UAG receptor has not been discovered.

Methodology/Principal Findings

To assess potential mechanisms of UAG action in vivo, we examined rapid effects of UAG on genome-wide expression patterns in fat, muscle and liver of growth hormone secretagogue receptor (GHSR)-ablated mice using microarrays. Expression data were analyzed using Ingenuity Pathways Analysis and Gene Set Enrichment Analysis. Regulation of subsets of these genes was verified by quantitative PCR in an independent experiment. UAG acutely regulated clusters of genes involved in glucose and lipid metabolism in all three tissues, consistent with enhancement of insulin sensitivity.

Conclusions/Significance

Fat, muscle and liver are central to the control of lipid and glucose homeostasis. UAG rapidly modulates the expression of metabolically important genes in these tissues in GHSR-deleted mice indicating a direct, GHSR-independent, action of UAG to improve insulin sensitivity and metabolic profile.

Physiological thyroid hormone levels regulate numerous skeletal muscle transcripts

CONTEXT

Skeletal muscle is an important target tissue for thyroid hormone (TH). It is currently unknown which genes are regulated by physiological TH levels.

OBJECTIVE

We examined the effects of l-thyroxine on human skeletal muscle transcriptome.

DESIGN

Microarray analysis of transcript levels was performed using skeletal muscle biopsies from patients under euthyroid and hypothyroid conditions.

SETTING

The study was conducted in a university hospital laboratory.

PATIENTS

We studied skeletal muscle obtained from 10 thyroidectomized patients with differentiated thyroid carcinoma on and after 4 wk off L-thyroxine replacement.

MEAN OUTCOME MEASURES

Gene expression changes were measured using microarrays. Results were analyzed using dedicated statistical methods.

RESULTS

We detected 607 differentially expressed genes on L-thyroxine treatment, of which approximately 60% were positively and approximately 40% were negatively regulated. Representative genes were validated by quantitative PCR. Genes involved in energy and fuel metabolism were overrepresented among the up-regulated genes, of which a large number were newly associated with thyroid state. L-thyroxine therapy induced a large down-regulation of the primary transcripts of the noncoding microRNA pair miR-206/miR-133b.

CONCLUSION

We demonstrated that physiological levels of TH regulate a myriad of genes in human skeletal muscle. The identification of novel putatively TH-responsive genes may provide the molecular basis of clinical effects in subjects with different TH status. The observation that TH regulates microRNAs reveals a new layer of complexity by which TH influences cellular processes.

HPV related VIN: highly proliferative and diminished responsiveness to extracellular signals

Vulvar intraepithelial neoplasia (VIN) is a premalignant disorder caused by human papillomaviruses. Basic knowledge about the molecular pathogenesis of VIN is sparse. Therefore, we have analyzed the gene expression profile of 9 VIN samples in comparison to 10 control samples by using genome wide Affymetrix Human U133A plus2 GeneChips. Results were validated by quantitative real-time RT-PCR analysis and immunostaining of a few representative genes (TACSTD1, CCNE2, AR and ESR1). Significance analysis of microarrays (SAM) showed that 1,497 genes were differentially expressed in VIN compared to controls. By analyzing the biological processes affected by the observed differences, we found that VIN appears to be a highly proliferative disease; many cyclins (CCNA, CCNB and CCNE) and almost all prereplication complex proteins are upregulated. Thereby, VIN does not seem to depend for its proliferation on paracrine or endocrine signals. Many receptors (for example ESR1 and AR) and ligands are downregulated. Furthermore, although VIN is not an invasive disease, the inhibition of expression of a marked number of cell-cell adhesion molecules seems to indicate development towards invasion. Upon reviewing apoptosis and angiogenesis, it was observed that these processes have not become significantly disregulated in VIN.

IN CONCLUSION

although VIN is still a premalignant disease, it already displays several hallmarks of cancer.

Functional differentiation of SWI/SNF remodelers in transcription and cell cycle control

Drosophila BAP and PBAP represent two evolutionarily conserved subclasses of SWI/SNF chromatin remodelers. The two complexes share the same core subunits, including the BRM ATPase, but differ in a few signature subunits: OSA defines BAP, whereas Polybromo (PB) and BAP170 specify PBAP. Here, we present a comprehensive structure-function analysis of BAP and PBAP. An RNA interference knockdown survey revealed that the core subunits BRM and MOR are critical for the structural integrity of both complexes. Whole-genome expression profiling suggested that the SWI/SNF core complex is largely dysfunctional in cells. Regulation of the majority of target genes required the signature subunit OSA, PB, or BAP170, suggesting that SWI/SNF remodelers function mostly as holoenzymes. BAP and PBAP execute similar, independent, or antagonistic functions in transcription control and appear to direct mostly distinct biological processes. BAP, but not PBAP, is required for cell cycle progression through mitosis. Because in yeast the PBAP-homologous complex, RSC, controls cell cycle progression, our finding reveals a functional switch during evolution. BAP mediates G(2)/M transition through direct regulation of string/cdc25. Its signature subunit, OSA, is required for directing BAP to the string/cdc25 promoter. Our results suggest that the core subunits play architectural and enzymatic roles but that the signature subunits determine most of the functional specificity of SWI/SNF holoenzymes in general gene control.