Primary oxidative phosphorylation defects lead to perturbations in the human B cell repertoire.

Introduction: The majority of studies on oxidative phosphorylation in immune cells have been performed in mouse models, necessitating human translation. To understand the impact of oxidative phosphorylation (OXPHOS) deficiency on human immunity, we studied children with primary mitochondrial disease (MtD). Methods: scRNAseq analysis of peripheral blood mononuclear cells was performed on matched children with MtD (N = 4) and controls (N = 4). To define B cell function we performed phage display immunoprecipitation sequencing on a cohort of children with MtD (N = 19) and controls (N = 16). Results: Via scRNAseq, we found marked reductions in select populations involved in the humoral immune response, especially antigen presenting cells, B cell and plasma populations, with sparing of T cell populations. MTRNR2L8, a marker of bioenergetic stress, was significantly elevated in populations that were most depleted. mir4485, a miRNA contained in the intron of MTRNR2L8, was co-expressed. Knockdown studies of mir4485 demonstrated its role in promoting survival by modulating apoptosis. To determine the functional consequences of our findings on humoral immunity, we studied the antiviral antibody repertoire in children with MtD and controls using phage display and immunoprecipitation sequencing. Despite similar viral exposomes, MtD displayed antiviral antibodies with less robust fold changes and limited polyclonality. Discussion: Overall, we show that children with MtD display perturbations in the B cell repertoire which may impact humoral immunity and the ability to clear viral infections.

Mapping the cortico-striatal transcriptome in attention deficit hyperactivity disorder.

Despite advances in identifying rare and common genetic variants conferring risk for ADHD, the lack of a transcriptomic understanding of cortico-striatal brain circuitry has stymied a molecular mechanistic understanding of this disorder. To address this gap, we mapped the transcriptome of the caudate nucleus and anterior cingulate cortex in post-mortem tissue from 60 individuals with and without ADHD. Significant differential expression of genes was found in the anterior cingulate cortex and, to a lesser extent, the caudate. Significant downregulation emerged of neurotransmitter gene pathways, particularly glutamatergic, in keeping with models that implicate these neurotransmitters in ADHD. Consistent with the genetic overlap between mental disorders, correlations were found between the cortico-striatal transcriptomic changes seen in ADHD and those seen in other neurodevelopmental and mood disorders. This transcriptomic evidence points to cortico-striatal neurotransmitter anomalies in the pathogenesis of ADHD, consistent with current models of the disorder.

Single-cell transcriptomics from human pancreatic islets: sample preparation matters.

Single-cell RNA sequencing (scRNA-seq) of human primary tissues is a rapidly emerging tool for investigating human health and disease at the molecular level. However, optimal processing of solid tissues presents a number of technical and logistical challenges, especially for tissues that are only available at autopsy, which includes pancreatic islets, a tissue that is highly relevant to diabetes. To assess the possible effects of different sample preparation protocols on fresh islet samples, we performed a detailed comparison of scRNA-seq data generated with islets isolated from a human donor but processed according to four treatment strategies, including fixation and cryopreservation. We found significant and reproducible differences in the proportion of cell types identified, and more minor effects on cell-specific patterns of gene expression. Fresh islets from a second donor confirmed gene expression signatures of alpha and beta subclusters. These findings may well apply to other tissues, emphasizing the need for careful consideration when choosing processing methods, comparing results between different studies, and/or interpreting data in the context of multiple cell types from preserved tissue.

MEK inhibition remodels the active chromatin landscape and induces SOX10 genomic recruitment in BRAF(V600E) mutant melanoma cells.

BACKGROUND: The MAPK/ERK signaling pathway is an essential regulator of numerous cell processes that are crucial for normal development as well as cancer progression. While much is known regarding MAPK/ERK signal conveyance from the cell membrane to the nucleus, the transcriptional and epigenetic mechanisms that govern gene expression downstream of MAPK signaling are not fully elucidated. RESULTS: This study employed an integrated epigenome analysis approach to interrogate the effects of MAPK/ERK pathway inhibition on the global transcriptome, the active chromatin landscape, and protein-DNA interactions in 501mel melanoma cells. Treatment of these cells with the small-molecule MEK inhibitor AZD6244 induces hyperpigmentation, widespread gene expression changes including alteration of genes linked to pigmentation, and extensive epigenomic reprogramming of transcriptionally distinct regulatory regions associated with the active chromatin mark H3K27ac. Regulatory regions with differentially acetylated H3K27ac regions following AZD6244 treatment are enriched in transcription factor binding motifs of ETV/ETS and ATF family members as well as the lineage-determining factors MITF and SOX10. H3K27ac-dense enhancer clusters known as super-enhancers show similar transcription factor motif enrichment, and furthermore, these super-enhancers are associated with genes encoding MITF, SOX10, and ETV/ETS proteins. Along with genome-wide resetting of the active enhancer landscape, MEK inhibition also results in widespread SOX10 recruitment throughout the genome, including increased SOX10 binding density at H3K27ac-marked enhancers. Importantly, these MEK inhibitor-responsive enhancers marked by H3K27ac and occupied by SOX10 are located near melanocyte lineage-specific and pigmentation genes and overlap numerous human SNPs associated with pigmentation and melanoma phenotypes, highlighting the variants located within these regions for prioritization in future studies. CONCLUSIONS: These results reveal the epigenetic reprogramming underlying the re-activation of melanocyte pigmentation and developmental transcriptional programs in 501mel cells in response to MEK inhibition and suggest extensive involvement of a MEK-SOX10 axis in the regulation of these processes. The dynamic chromatin changes identified here provide a rich genomic resource for further analyses of the molecular mechanisms governing the MAPK pathway in pigmentation- and melanocyte-associated diseases.

Prostate cancer susceptibility gene HIST1H1A is a modulator of androgen receptor signaling and epithelial to mesenchymal transition.

In 2018, approximately 165,000 new prostate cancer (PC) cases will be diagnosed, and over 29,000 men will succumb to PC in the U.S. alone. The means of assessing outcome in the clinic are inaccurate, and there is a pressing need to more precisely identify men at risk of aggressive PC. We previously identified HIST1H1A as a susceptibility gene for aggressive PC. HIST1H1A encodes H1.1, a member of the linker histone family that is involved in chromatin organization and compaction. To understand the molecular basis of aggressive PC, we have characterized how germline variation modulates susceptibility to neuroendocrine differentiation, which is a form of aggressive PC. Immunohistochemistry studies revealed that HIST1H1A is over-expressed in normal human prostate tissue compared to prostate adenocarcinoma. Functional characterization of HIST1H1A in prostate LNCaP cells indicated that HIST1HA over-expression increased cell growth, as well as the expression of neuroendocrine and epithelial-to-mesenchymal markers in vitro. Assay for Transposase-Accessible Chromatin (ATAC-seq), which is used to assess chromatin compaction and thus the transcriptional availability of individual genomic regions, demonstrated that H1.1 plays a prominent role in modulating Wnt signaling pathway genes, which are implicated in prostate tumorigenesis. These results demonstrate that HIST1H1A is a modulator of aggressive PC susceptibility.

A direct link between MITF, innate immunity, and hair graying.

Melanocyte stem cells (McSCs) and mouse models of hair graying serve as useful systems to uncover mechanisms involved in stem cell self-renewal and the maintenance of regenerating tissues. Interested in assessing genetic variants that influence McSC maintenance, we found previously that heterozygosity for the melanogenesis associated transcription factor, Mitf, exacerbates McSC differentiation and hair graying in mice that are predisposed for this phenotype. Based on transcriptome and molecular analyses of Mitfmi-vga9/+ mice, we report a novel role for MITF in the regulation of systemic innate immune gene expression. We also demonstrate that the viral mimic poly(I:C) is sufficient to expose genetic susceptibility to hair graying. These observations point to a critical suppressor of innate immunity, the consequences of innate immune dysregulation on pigmentation, both of which may have implications in the autoimmune, depigmenting disease, vitiligo.

TFAP2 paralogs regulate melanocyte differentiation in parallel with MITF.

Mutations in the gene encoding transcription factor TFAP2A result in pigmentation anomalies in model organisms and premature hair graying in humans. However, the pleiotropic functions of TFAP2A and its redundantly-acting paralogs have made the precise contribution of TFAP2-type activity to melanocyte differentiation unclear. Defining this contribution may help to explain why TFAP2A expression is reduced in advanced-stage melanoma compared to benign nevi. To identify genes with TFAP2A-dependent expression in melanocytes, we profile zebrafish tissue and mouse melanocytes deficient in Tfap2a, and find that expression of a small subset of genes underlying pigmentation phenotypes is TFAP2A-dependent, including Dct, Mc1r, Mlph, and Pmel. We then conduct TFAP2A ChIP-seq in mouse and human melanocytes and find that a much larger subset of pigmentation genes is associated with active regulatory elements bound by TFAP2A. These elements are also frequently bound by MITF, which is considered the "master regulator" of melanocyte development. For example, the promoter of TRPM1 is bound by both TFAP2A and MITF, and we show that the activity of a minimal TRPM1 promoter is lost upon deletion of the TFAP2A binding sites. However, the expression of Trpm1 is not TFAP2A-dependent, implying that additional TFAP2 paralogs function redundantly to drive melanocyte differentiation, which is consistent with previous results from zebrafish. Paralogs Tfap2a and Tfap2b are both expressed in mouse melanocytes, and we show that mouse embryos with Wnt1-Cre-mediated deletion of Tfap2a and Tfap2b in the neural crest almost completely lack melanocytes but retain neural crest-derived sensory ganglia. These results suggest that TFAP2 paralogs, like MITF, are also necessary for induction of the melanocyte lineage. Finally, we observe a genetic interaction between tfap2a and mitfa in zebrafish, but find that artificially elevating expression of tfap2a does not increase levels of melanin in mitfa hypomorphic or loss-of-function mutants. Collectively, these results show that TFAP2 paralogs, operating alongside lineage-specific transcription factors such as MITF, directly regulate effectors of terminal differentiation in melanocytes. In addition, they suggest that TFAP2A activity, like MITF activity, has the potential to modulate the phenotype of melanoma cells.

Mapping Complex Traits in a Diversity Outbred F1 Mouse Population Identifies Germline Modifiers of Metastasis in Human Prostate Cancer.

It is unclear how standing genetic variation affects the prognosis of prostate cancer patients. To provide one controlled answer to this problem, we crossed a dominant, penetrant mouse model of prostate cancer to Diversity Outbred mice, a collection of animals that carries over 40 million SNPs. Integration of disease phenotype and SNP variation data in 493 F1 males identified a metastasis modifier locus on Chromosome 8 (LOD = 8.42); further analysis identified the genes Rwdd4, Cenpu, and Casp3 as functional effectors of this locus. Accordingly, analysis of over 5,300 prostate cancer patient samples revealed correlations between the presence of genetic variants at these loci, their expression levels, cancer aggressiveness, and patient survival. We also observed that ectopic overexpression of RWDD4 and CENPU increased the aggressiveness of two human prostate cancer cell lines. In aggregate, our approach demonstrates how well-characterized genetic variation in mice can be harnessed in conjunction with systems genetics approaches to identify and characterize germline modifiers of human disease processes.

Genomic analysis reveals distinct mechanisms and functional classes of SOX10-regulated genes in melanocytes.

SOX10 is required for melanocyte development and maintenance, and has been linked to melanoma initiation and progression. However, the molecular mechanisms by which SOX10 guides the appropriate gene expression programs necessary to promote the melanocyte lineage are not fully understood. Here we employ genetic and epigenomic analysis approaches to uncover novel genomic targets and previously unappreciated molecular roles of SOX10 in melanocytes. Through global analysis of SOX10-binding sites and epigenetic characteristics of chromatin states, we uncover an extensive catalog of SOX10 targets genome-wide. Our findings reveal that SOX10 predominantly engages 'open' chromatin regions and binds to distal regulatory elements, including novel and previously known melanocyte enhancers. Integrated chromatin occupancy and transcriptome analysis suggest a role for SOX10 in both transcriptional activation and repression to regulate functionally distinct classes of genes. We demonstrate that distinct epigenetic signatures and cis-regulatory sequence motifs predicted to bind putative co-regulatory transcription factors define SOX10-activated and SOX10-repressed target genes. Collectively, these findings uncover a central role of SOX10 as a global regulator of gene expression in the melanocyte lineage by targeting diverse regulatory pathways.

MLV integration site selection is driven by strong enhancers and active promoters.

Retroviruses integrate into the host genome in patterns specific to each virus. Understanding the causes of these patterns can provide insight into viral integration mechanisms, pathology and genome evolution, and is critical to the development of safe gene therapy vectors. We generated murine leukemia virus integrations in human HepG2 and K562 cells and subjected them to second-generation sequencing, using a DNA barcoding technique that allowed us to quantify independent integration events. We characterized >3,700,000 unique integration events in two ENCODE-characterized cell lines. We find that integrations were most highly enriched in a subset of strong enhancers and active promoters. In both cell types, approximately half the integrations were found in <2% of the genome, demonstrating genomic influences even narrower than previously believed. The integration pattern of murine leukemia virus appears to be largely driven by regions that have high enrichment for multiple marks of active chromatin; the combination of histone marks present was sufficient to explain why some strong enhancers were more prone to integration than others. The approach we used is applicable to analyzing the integration pattern of any exogenous element and could be a valuable preclinical screen to evaluate the safety of gene therapy vectors.

A large-scale zebrafish gene knockout resource for the genome-wide study of gene function.

With the completion of the zebrafish genome sequencing project, it becomes possible to analyze the function of zebrafish genes in a systematic way. The first step in such an analysis is to inactivate each protein-coding gene by targeted or random mutation. Here we describe a streamlined pipeline using proviral insertions coupled with high-throughput sequencing and mapping technologies to widely mutagenize genes in the zebrafish genome. We also report the first 6144 mutagenized and archived F1's predicted to carry up to 3776 mutations in annotated genes. Using in vitro fertilization, we have rescued and characterized ~0.5% of the predicted mutations, showing mutation efficacy and a variety of phenotypes relevant to both developmental processes and human genetic diseases. Mutagenized fish lines are being made freely available to the public through the Zebrafish International Resource Center. These fish lines establish an important milestone for zebrafish genetics research and should greatly facilitate systematic functional studies of the vertebrate genome.

The Zebrafish Insertion Collection (ZInC): a web based, searchable collection of zebrafish mutations generated by DNA insertion.

ZInC (Zebrafish Insertional Collection, http://research.nhgri.nih.gov/ZInC/) is a web-searchable interface of insertional mutants in zebrafish. Over the last two decades, the zebrafish has become a popular model organism for studying vertebrate development as well as for modeling human diseases. To facilitate such studies, we are generating a genome-wide knockout resource that targets every zebrafish protein-coding gene. All mutant fish are freely available to the scientific community through the Zebrafish International Resource Center (ZIRC). To assist researchers in finding mutant and insertion information, we developed a comprehensive database with a web front-end, the ZInC. It can be queried using multiple types of input such as ZFIN (Zebrafish Information Network) IDs, UniGene accession numbers and gene symbols from zebrafish, human and mouse. In the future, ZInC may include data from other insertional mutation projects as well. ZInC cross-references all integration data with the ZFIN (http://zfin.org/).

The pleiotropic mouse phenotype extra-toes spotting is caused by translation initiation factor Eif3c mutations and is associated with disrupted sonic hedgehog signaling.

Polydactyly is a common malformation and can be an isolated anomaly or part of a pleiotropic syndrome. The elucidation of the mutated genes that cause polydactyly provides insight into limb development pathways. The extra-toes spotting (Xs) mouse phenotype manifests anterior polydactyly, predominantly in the forelimbs, with ventral hypopigmenation. The mapping of Xs(J) to chromosome 7 was confirmed, and the interval was narrowed to 322 kb using intersubspecific crosses. Two mutations were identified in eukaryotic translation initiation factor 3 subunit C (Eif3c). An Eif3c c.907C>T mutation (p.Arg303X) was identified in Xs(J), and a c.1702_1758del mutation (p.Leu568_Leu586del) was identified in extra-toes spotting-like (Xsl), an allele of Xs(J). The effect of the Xs(J) mutation on the SHH/GLI3 pathway was analyzed by in situ hybridization analysis, and we show that Xs mouse embryos have ectopic Shh and Ptch1 expression in the anterior limb. In addition, anterior limb buds show aberrant Gli3 processing, consistent with perturbed SHH/GLI3 signaling. Based on the occurrence of Eif3c mutations in 2 Xs lines and haploinsufficiency of the Xs(J) allele, we conclude that the Xs phenotype is caused by a mutation in Eif3c, a component of the translation initiation complex, and that the phenotype is associated with aberrant SHH/GLI3 signaling.

Genome-wide scan of Swedish families with hereditary prostate cancer: suggestive evidence of linkage at 5q11.2 and 19p13.3.

BACKGROUND: Prostate cancer (CaP) is a common disorder with multiple genetic and environmental factors contributing to the disease. CaP susceptibility loci can be identified through genome-wide scans of high-risk families. METHODS: Allele sharing at 405 markers, distributed across the genome, among 50 families with hereditary prostate cancer, ascertained throughout Sweden, was evaluated through linkage analyses. Genotype data were analyzed utilizing multipoint parametric and non-parametric methods. RESULTS: Two regions provided suggestive evidence for linkage: 19p13.3 (marker D19S209, LOD = 2.91, P = 0.0001) and 5q11.2 (marker D5S407, LOD = 2.24, P = 0.0007). Additional regions with moderate evidence for linkage in the complete set of families, or stratified subsets, were observed on chromosome 1, 4, 6, 7, 8, and X. CONCLUSIONS: Our results provide strong confirmatory evidence of linkage at 19q13.3 and 5q11.2. The lack of confirmation of linkage at several loci identified in other genome-wide scans emphasizes the need to combine linkage data between research groups.