Search

Abstract Background Rhesus macaques (Macaca mulatta) are the primate most used for biomedical research, but phenotypic differences between Indian-origin and Chinese rhesus macaques have encouraged genetic methods for identifying genetic differences between these two populations. The completion of the rhesus genome has led to the identification of many single nucleotide polymorphisms (SNPs) in this species. These single nucleotide polymorphisms have many advantages over the short tandem repeat (STR) loci currently used to assay genetic variation. However, the number of currently identified polymorphisms is too small for whole genome analysis or studies of quantitative trait loci. To that end, we tested a combination of methods to identify large numbers of high-confidence SNPs, and screen those with high minor allele frequencies (MAF). Results By testing our previously reported single nucleotide polymorphisms, we identified a subset of high-confidence, high-MAF polymorphisms. Resequencing revealed a large number of regionally specific SNPs not identified through a single pyrosequencing run. By resequencing a pooled sample of four individuals, we reliably identified loci with a MAF of at least 12.5%. Finally, we found that when applied to a larger, geographically variable sample of rhesus, a large proportion of our loci were variable in both populations, and very few loci were ancestry informative. Despite this fact, the SNP loci were more effective at discriminating Indian and Chinese rhesus than STR loci. Conclusion Pyrosequencing and pooled resequencing are viable methods for the identification of high-MAF SNP loci in rhesus macaques. These SNP loci are appropriate for screening both the inter- and intra-population genetic variation.

Although significant progress has been achieved in elucidating the genetic architecture of systemic lupus erythematosus (SLE), identifying genes underlying the pathogenesis has been challenging. The NZM2410-derived lupus susceptibility Sle3 locus is associated with T cell hyperactivity and activated myeloid cells. However, candidate genes associated with these phenotypes have not been identified. Here, we narrow the Sle3 locus to a smaller genomic segment (Sle3k) and show that mice carrying Sle3k and Sle1 loci developed lupus nephritis. We identify Klf13 as the primary candidate gene that is associated with genome-wide transcription changes resulting in higher levels of proinflammatory cytokines, enhanced T cell activation, and hyperresponsiveness of myeloid cells. Correspondingly, Klf13 -/- mice display repression of genes involved in mediating immune activation, including key proinflammatory cytokines/chemokines in T cells and dysregulation in cytokine signaling pathways in myeloid cells in response to toll receptor ligands. Klf13 upregulation is associated with increased production of RANTES, a key chemokine in lupus nephritis, in activated T cells and the kidneys of lupus-prone mice. In sum, our findings reveal Klf13 as a key gene in the Sle3 interval in mediating lupus pathogenesis that may have implications in the rational design of new therapies for SLE., (© 2024. The Author(s).)

Cancer cachexia (CC), a wasting syndrome of muscle and adipose tissue resulting in weight loss, is observed in 50% of patients with solid tumors. Management of CC is limited by the absence of biomarkers and knowledge of molecules that drive its phenotype. To identify such molecules, we injected 54 human non-small cell lung cancer (NSCLC) lines into immunodeficient mice, 17 of which produced an unambiguous phenotype of cachexia or non-cachexia. Whole-exome sequencing revealed that 8 of 10 cachexia lines, but none of the non-cachexia lines, possessed mutations in serine/threonine kinase 11 (STK11/LKB1), a regulator of nutrient sensor AMPK. Silencing of STK11/LKB1 in human NSCLC and murine colorectal carcinoma lines conferred a cachexia phenotype after cell transplantation into immunodeficient (human NSCLC) and immunocompetent (murine colorectal carcinoma) models. This host wasting was associated with an alteration in the immune cell repertoire of the tumor microenvironments that led to increases in local mRNA expression and serum levels of CC-associated cytokines. Mutational analysis of circulating tumor DNA from patients with NSCLC identified 89% concordance between STK11/LKB1 mutations and weight loss at cancer diagnosis. The current data provide evidence that tumor STK11/LKB1 loss of function is a driver of CC, simultaneously serving as a genetic biomarker for this wasting syndrome.

Inflammation plays a central role in immune response and macrophage activation. Emerging studies demonstrate that along with proteins and genomic factors, noncoding RNA are potentially involved in regulation of immune response and inflammation. Our recent study demonstrated that lncRNA HOTAIR plays key roles in cytokine expression and inflammation in macrophages. The primary goal of this study is to discover novel lncRNAs that are crucial players in inflammation, macrophage activation, and immune response in humans. Towards this, we have stimulated THP1-derived macrophages (THP1-MΦ) with lipopolysaccharides (LPS) and performed the whole transcriptome RNA-seq analysis. Based on this analysis, we discovered that along with well-known marker for inflammation (such as cytokines), a series of long noncoding RNAs (lncRNAs) expression were highly induced upon LPS-stimulation of macrophages, suggesting their potential roles in inflammation and macrophage activation. We termed these family of lncRNAs as Long-noncoding Inflammation Associated RNA (LinfRNA). Dose and time dependent analysis demonstrated that many human LinfRNA (hLinfRNAs) expressions follow similar patterns as cytokine expressions. Inhibition of NF-κB suppressed the expression of most hLinfRNAs suggesting their potential regulation via NF-κB activation during inflammation and macrophage activation. Antisense-mediated knockdown of hLinfRNA1 suppressed the LPS-induced expression of cytokines and pro-inflammatory genes such as IL6, IL1β, and TNFα expression, suggesting potential functionality of the hLinfRNAs in cytokine regulation and inflammation. Overall, we discovered a series of novel hLinfRNAs that are potential regulators of inflammation and macrophage activation and may be linked to inflammatory and metabolic diseases., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Acquired mutations in the ligand-binding domain (LBD) of the gene encoding estrogen receptor α (ESR1) are common mechanisms of endocrine therapy resistance in patients with metastatic ER+ breast cancer. The ESR1 Y537S mutation, in particular, is associated with development of resistance to most endocrine therapies used to treat breast cancer. Employing a high-throughput screen of nearly 1,200 Federal Drug Administration-approved (FDA-approved) drugs, we show that OTX015, a bromodomain and extraterminal domain (BET) inhibitor, is one of the top suppressors of ESR1 mutant cell growth. OTX015 was more efficacious than fulvestrant, a selective ER degrader, in inhibiting ESR1 mutant xenograft growth. When combined with abemaciclib, a CDK4/6 inhibitor, OTX015 induced more potent tumor regression than current standard-of-care treatment of abemaciclib + fulvestrant. OTX015 has preferential activity against Y537S mutant breast cancer cells and blocks their clonal selection in competition studies with WT cells. Thus, BET inhibition has the potential to both prevent and overcome ESR1 mutant-induced endocrine therapy resistance in breast cancer.

The MYC oncogene is frequently amplified in triple-negative breast cancer (TNBC). Here, we show that MYC suppression induces immune-related hallmark gene set expression and tumor-infiltrating T cells in MYC-hyperactivated TNBCs. Mechanistically, MYC repressed stimulator of interferon genes (STING) expression via direct binding to the STING1 enhancer region, resulting in downregulation of the T-cell chemokines CCL5, CXCL10, and CXCL11. In primary and metastatic TNBC cohorts, tumors with high MYC expression or activity exhibited low STING expression. Using a CRISPR-mediated enhancer perturbation approach, we demonstrated that MYC-driven immune evasion is mediated by STING repression. STING repression induced resistance to PD-L1 blockade in mouse models of TNBC. Finally, a small-molecule inhibitor of MYC combined with PD-L1 blockade elicited a durable response in immune-cold TNBC with high MYC expression, suggesting a strategy to restore PD-L1 inhibitor sensitivity in MYC-overexpressing TNBC., (©2022 American Association for Cancer Research.)

Low-grade oncocytic tumor (LOT) of the kidney is a recently described entity with poorly understood pathogenesis. Using next-generation sequencing (NGS) and complementary approaches, we provide insight into its biology. We describe 22 LOT corresponding to 7 patients presenting with a median age of 75 years (range 63-86 years) and male to female ratio 2:5. All 22 tumors demonstrated prototypical microscopic features. Tumors were well-circumscribed and solid. They were composed of sheets of tumor cells in compact nests. Tumor cells had eosinophilic cytoplasm, round to oval nuclei (without nuclear membrane irregularities), focal subtle perinuclear halos, and occasional binucleation. Sharply delineated edematous stromal islands were often observed. Tumor cells were positive for PAX8, negative for CD117, and exhibited diffuse and strong cytokeratin-7 expression. Six patients presented with pT1 tumors. At a median follow-up of 29 months, four patients were alive without recurrence (three patients had died from unrelated causes). All tumors were originally classified as chromophobe renal cell carcinoma, eosinophilic variant (chRCC-eo). While none of the patients presented with known syndromic features, one patient with multiple bilateral LOTs was subsequently found to have a likely pathogenic germline TSC1 mutation. Somatic, likely activating, mutations in MTOR and RHEB were identified in all other evaluable LOTs. As assessed by phospho-S6 and phospho-4E-BP1, mTOR complex 1 (mTORC1) was activated across all cases but to different extent. MTOR mutant LOT exhibited lower levels of mTORC1 activation, possibly related to mTORC1 dimerization and the preservation of a wild-type MTOR copy (retained chromosome 1). Supporting its distinction from related entities, gene expression analyses showed that LOT clustered separately from classic chRCC, chRCC-eo, and RO. In summary, converging mTORC1 pathway mutations, mTORC1 complex activation, and a distinctive gene expression signature along with characteristic phenotypic features support LOT designation as a distinct entity with both syndromic and non-syndromic cases associated with an indolent course., (© 2021. The Author(s), under exclusive licence to United States & Canadian Academy of Pathology.)

Fundamental principles of HIV-1 integration into the human genome have been revealed in the past 2 decades. However, the impact of the integration site on proviral transcription and expression remains poorly understood. Solving this problem requires the analysis of multiple genomic datasets for thousands of proviral integration sites. Here, we generated and combined large-scale datasets, including epigenetics, transcriptome, and 3-dimensional genome architecture to interrogate the chromatin states, transcription activity, and nuclear sub-compartments around HIV-1 integrations in Jurkat CD4 + T cells to decipher human genome regulatory features shaping the transcription of proviral classes based on their position and orientation in the genome. Through a Hidden Markov Model and ranked informative values prior to a machine learning logistic regression model, we defined nuclear sub-compartments and chromatin states contributing to genomic architecture, transcriptional activity, and nucleosome density of regions neighboring the integration site, as additive features influencing HIV-1 expression. Our integrated genomics approach also allows for a robust experimental design, in which HIV-1 can be genetically introduced into precise genomic locations with known regulatory features to assess the relationship of integration positions to viral transcription and fate., Competing Interests: Declaration of Conflicting Interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2022.)

Acetyl-CoA is a key intermediate situated at the intersection of many metabolic pathways. The reliance of histone acetylation on acetyl-CoA enables the coordination of gene expression with metabolic state. Abundant acetyl-CoA has been linked to the activation of genes involved in cell growth or tumorigenesis through histone acetylation. However, the role of histone acetylation in transcription under low levels of acetyl-CoA remains poorly understood. Here, we use a yeast starvation model to observe the dramatic alteration in the global occupancy of histone acetylation following carbon starvation; the location of histone acetylation marks shifts from growth-promoting genes to gluconeogenic and fat metabolism genes. This reallocation is mediated by both the histone deacetylase Rpd3p and the acetyltransferase Gcn5p, a component of the SAGA transcriptional coactivator. Our findings reveal an unexpected switch in the specificity of histone acetylation to promote pathways that generate acetyl-CoA for oxidation when acetyl-CoA is limiting., Competing Interests: Declaration of interests B.P.T. is a member of the advisory board for Molecular Cell., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Context.—: Active surveillance of small renal masses highlights the need for accurate prognostication of biopsies., Objective.—: To comprehensively evaluate the accuracy of biopsies in assessing known prognostic parameters including histologic subtype by comparison with subsequent nephrectomy samples., Design.—: We retrospectively identified patients at University of Texas Southwestern Medical Center, Dallas, Texas, who had a biopsy for a renal mass between 2004-2018. Biopsy samples were evaluated for known prognostic factors such as tumor grade, necrosis, sarcomatoid/rhabdoid change, and BRCA1-associated protein-1 (BAP1) status, which we previously showed is an independent prognostic factor for clear cell renal cell carcinoma. Accuracy was determined by comparison with subsequent analyses of nephrectomy specimens. Statistical analyses were performed to assess biopsy accuracy and correlation with tumor size and pathologic stage., Results.—: From 805 biopsies with a diagnosis of renal neoplasm, 178 had subsequent resection of the biopsied tumor. Concordance rate for histologic subtype was 96.9% (κ [w], 0.90; 95% CI, 0.82-0.99) and excellent for small renal masses (98.8%; κ [w], 0.97; 95% CI, 0.90-1). Amongst the prognostic variables evaluated, BAP1 immunohistochemistry in clear cell renal cell carcinoma had the highest agreement (94.8%; κ [w], 0.83; 95% CI, 0.66-0.99). The presence of 1 or more aggressive features (grade 3-4, tumor necrosis, BAP1 loss, sarcomatoid/rhabdoid change) in a biopsy significantly correlated with pT stage (P = .004)., Conclusions.—: Biopsy analyses showed high accuracy for subtyping renal tumors, but it underestimated several poor prognostic features. Addition of BAP1 for clear cell renal cell carcinoma may increase prognostic accuracy. If validated, routine incorporation of BAP1 immunohistochemistry in clear cell renal cell carcinoma biopsies may refine prognosis and aid in the selection of patients for active surveillance.

Deconvolution of regulatory mechanisms that drive transcriptional programs in cancer cells is key to understanding tumor biology. Herein, we present matched transcriptome (scRNA-seq) and chromatin accessibility (scATAC-seq) profiles at single-cell resolution from human ovarian and endometrial tumors processed immediately following surgical resection. This dataset reveals the complex cellular heterogeneity of these tumors and enabled us to quantitatively link variation in chromatin accessibility to gene expression. We show that malignant cells acquire previously unannotated regulatory elements to drive hallmark cancer pathways. Moreover, malignant cells from within the same patients show substantial variation in chromatin accessibility linked to transcriptional output, highlighting the importance of intratumoral heterogeneity. Finally, we infer the malignant cell type-specific activity of transcription factors. By defining the regulatory logic of cancer cells, this work reveals an important reliance on oncogenic regulatory elements and highlights the ability of matched scRNA-seq/scATAC-seq to uncover clinically relevant mechanisms of tumorigenesis in gynecologic cancers., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Stromal-epithelial interactions are critical to the morphogenesis, differentiation, and homeostasis of the prostate, but the molecular identity and anatomy of discrete stromal cell types is poorly understood. Using single-cell RNA sequencing, we identified and validated the in situ localization of three smooth muscle subtypes (prostate smooth muscle, pericytes, and vascular smooth muscle) and two novel fibroblast subtypes in human prostate. Peri-epithelial fibroblasts (APOD+) wrap around epithelial structures, whereas interstitial fibroblasts (C7+) are interspersed in extracellular matrix. In contrast, the mouse displayed three fibroblast subtypes with distinct proximal-distal and lobe-specific distribution patterns. Statistical analysis of mouse and human fibroblasts showed transcriptional correlation between mouse prostate (C3+) and urethral (Lgr5+) fibroblasts and the human interstitial fibroblast subtype. Both urethral fibroblasts (Lgr5+) and ductal fibroblasts (Wnt2+) in the mouse contribute to a proximal Wnt/Tgfb signaling niche that is absent in human prostate. Instead, human peri-epithelial fibroblasts express secreted WNT inhibitors SFRPs and DKK1, which could serve as a buffer against stromal WNT ligands by creating a localized signaling niche around individual prostate glands. We also identified proximal-distal fibroblast density differences in human prostate that could amplify stromal signaling around proximal prostate ducts. In human benign prostatic hyperplasia, fibroblast subtypes upregulate critical immunoregulatory pathways and show distinct distributions in stromal and glandular phenotypes. A detailed taxonomy of leukocytes in benign prostatic hyperplasia reveals an influx of myeloid dendritic cells, T cells and B cells, resembling a mucosal inflammatory disorder. A receptor-ligand interaction analysis of all cell types revealed a central role for fibroblasts in growth factor, morphogen, and chemokine signaling to endothelia, epithelia, and leukocytes. These data are foundational to the development of new therapeutic targets in benign prostatic hyperplasia. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd., (© 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.)

Chronic low-grade white adipose tissue (WAT) inflammation is a hallmark of metabolic syndrome in obesity. Here, we demonstrate that a subpopulation of mouse WAT perivascular (PDGFRβ + ) cells, termed fibro-inflammatory progenitors (FIPs), activate proinflammatory signalling cascades shortly after the onset of high-fat diet feeding and regulate proinflammatory macrophage accumulation in WAT in a TLR4-dependent manner. FIPs activation in obesity is mediated by the downregulation of zinc-finger protein 423 (ZFP423), identified here as a transcriptional corepressor of NF-κB. ZFP423 suppresses the DNA-binding capacity of the p65 subunit of NF-κB by inducing a p300-to-NuRD coregulator switch. Doxycycline-inducible expression of Zfp423 in PDGFRβ +  cells suppresses inflammatory signalling in FIPs and attenuates metabolic inflammation of visceral WAT in obesity. Inducible inactivation of Zfp423 in PDGFRβ +  cells increases FIP activity, exacerbates adipose macrophage accrual and promotes WAT dysfunction. These studies implicate perivascular mesenchymal cells as important regulators of chronic adipose-tissue inflammation in obesity and identify ZFP423 as a transcriptional break on NF-κB signalling.

Background: Castration-insensitive epithelial progenitors capable of regenerating the prostate have been proposed to be concentrated in the proximal region based on facultative assays. Functional characterization of prostate epithelial populations isolated with individual cell surface markers has failed to provide a consensus on the anatomical and transcriptional identity of proximal prostate progenitors., Methods: Here, we use single-cell RNA sequencing to obtain a complete transcriptomic profile of all epithelial cells in the mouse prostate and urethra to objectively identify cellular subtypes. Pan-transcriptomic comparison to human prostate cell types identified a mouse equivalent of human urethral luminal cells, which highly expressed putative prostate progenitor markers. Validation of the urethral luminal cell cluster was performed using immunostaining and flow cytometry., Results: Our data reveal that previously identified facultative progenitors marked by Trop2, Sca-1, KRT4, and PSCA are actually luminal epithelial cells of the urethra that extend into the proximal region of the prostate, and are resistant to castration-induced androgen deprivation. Mouse urethral luminal cells were identified to be the equivalent of previously identified human club and hillock cells that similarly extend into proximal prostate ducts. Benign prostatic hyperplasia (BPH) has long been considered an "embryonic reawakening," but the cellular origin of the hyperplastic growth concentrated in the periurethral region is unclear. We demonstrate an increase in urethral luminal cells within glandular nodules from BPH patients. Urethral luminal cells are further increased in patients treated with a 5-α reductase inhibitor., Conclusions: Our data demonstrate that cells of the proximal prostate that express putative progenitor markers, and are enriched by castration in the proximal prostate, are urethral luminal cells and that these cells may play an important role in the etiology of human BPH., (© 2020 Wiley Periodicals LLC.)

The differentiation of embryonic stem cells into various lineages is highly dependent on the chromatin state of the genome and patterns of gene expression. To identify lineage-specific enhancers driving the differentiation of progenitors into pancreatic cells, we used a previously described computational framework called Total Functional Score of Enhancer Elements (TFSEE), which integrates multiple genomic assays that probe both transcriptional and epigenomic states. First, we evaluated and compared TFSEE as an enhancer-calling algorithm with enhancers called using GRO-seq-defined enhancer transcripts (method 1) versus enhancers called using histone modification ChIP-seq data (method 2). Second, we used TFSEE to define the enhancer landscape and identify transcription factors (TFs) that maintain the multipotency of a subpopulation of endodermal stem cells during differentiation into pancreatic lineages. Collectively, our results demonstrate that TFSEE is a robust enhancer-calling algorithm that can be used to perform multilayer genomic data integration to uncover cell type-specific TFs that control lineage-specific enhancers., Competing Interests: Declaration of conflicting interests:The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2020.)

Metastatic prostate cancer is characterized by recurrent genomic copy number alterations that are presumed to contribute to resistance to hormone therapy. We identified CHD1 loss as a cause of antiandrogen resistance in an in vivo small hairpin RNA (shRNA) screen of 730 genes deleted in prostate cancer. ATAC-seq and RNA-seq analyses showed that CHD1 loss resulted in global changes in open and closed chromatin with associated transcriptomic changes. Integrative analysis of this data, together with CRISPR-based functional screening, identified four transcription factors (NR3C1, POU3F2, NR2F1, and TBX2) that contribute to antiandrogen resistance, with associated activation of non-luminal lineage programs. Thus, CHD1 loss results in chromatin dysregulation, thereby establishing a state of transcriptional plasticity that enables the emergence of antiandrogen resistance through heterogeneous mechanisms., Competing Interests: Declaration of Interests C.L.S. and J.W. are co-inventors of enzalutamide and apalutamide and may be entitled to royalties. C.L.S. serves on the Board of Directors of Novartis and is a co-founder of ORIC Pharm. He is a science advisor to Agios, Beigene, Blueprint, Column Group, Foghorn, Housey Pharma, Nextech, KSQ, Petra, and PMV. He was a co-founder of Seragon, purchased by Genentech/Roche in 2014. S.W.L. is a founder and member of the scientific advisory board of ORIC Pharmaceuticals, Blueprint Medicines, and Mirimus, Inc.; he is also on the scientific advisory board of PMV Pharmaceuticals, Constellation Pharmaceuticals, and Petra Pharmaceuticals. W.A. reports consulting for Clovis Oncology, Janssen, MORE Health, and ORIC Pharmaceuticals. He received honoraria from CARET and travel accommodations from GlaxoSmith Kline, Clovis Oncology, and ORIC Pharmaceuticals. C.E.M is a co-founder and board member for Biotia and Onegevity Health, as well as an advisor for Genpro and Karius., (Copyright © 2020 Elsevier Inc. All rights reserved.)

PARP inhibitors (PARPi) prevent cancer cell growth by inducing synthetic lethality with DNA repair defects (e.g., in BRCA1/2 mutant cells). We have identified an alternative pathway for PARPi-mediated growth control in BRCA1/2-intact breast cancer cells involving rDNA transcription and ribosome biogenesis. PARP-1 binds to snoRNAs, which stimulate PARP-1 catalytic activity in the nucleolus independent of DNA damage. Activated PARP-1 ADP-ribosylates DDX21, an RNA helicase that localizes to nucleoli and promotes rDNA transcription when ADP-ribosylated. Treatment with PARPi or mutation of the ADP-ribosylation sites reduces DDX21 nucleolar localization, rDNA transcription, ribosome biogenesis, protein translation, and cell growth. The salient features of this pathway are evident in xenografts in mice and human breast cancer patient samples. Elevated levels of PARP-1 and nucleolar DDX21 are associated with cancer-related outcomes. Our studies provide a mechanistic rationale for efficacy of PARPi in cancer cells lacking defects in DNA repair whose growth is inhibited by PARPi., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The cardiogenic transcription factors (TFs) Mef2c, Gata4, and Tbx5 can directly reprogram fibroblasts to induced cardiac-like myocytes (iCLMs), presenting a potential source of cells for cardiac repair. While activity of these TFs is enhanced by Hand2 and Akt1, their genomic targets and interactions during reprogramming are not well studied. We performed genome-wide analyses of cardiogenic TF binding and enhancer profiling during cardiac reprogramming. We found that these TFs synergistically activate enhancers highlighted by Mef2c binding sites and that Hand2 and Akt1 coordinately recruit other TFs to enhancer elements. Intriguingly, these enhancer landscapes collectively resemble patterns of enhancer activation during embryonic cardiogenesis. We further constructed a cardiac reprogramming gene regulatory network and found repression of EGFR signaling pathway genes. Consistently, chemical inhibition of EGFR signaling augmented reprogramming. Thus, by defining epigenetic landscapes these findings reveal synergistic transcriptional activation across a broad landscape of cardiac enhancers and key signaling pathways that govern iCLM reprogramming., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Open, collaborative research is a powerful paradigm that can immensely strengthen the scientific process by integrating broad and diverse expertise. However, traditional research and multi-author writing processes break down at scale. We present new software named Manubot, available at https://manubot.org, to address the challenges of open scholarly writing. Manubot adopts the contribution workflow used by many large-scale open source software projects to enable collaborative authoring of scholarly manuscripts. With Manubot, manuscripts are written in Markdown and stored in a Git repository to precisely track changes over time. By hosting manuscript repositories publicly, such as on GitHub, multiple authors can simultaneously propose and review changes. A cloud service automatically evaluates proposed changes to catch errors. Publication with Manubot is continuous: When a manuscript's source changes, the rendered outputs are rebuilt and republished to a web page. Manubot automates bibliographic tasks by implementing citation by identifier, where users cite persistent identifiers (e.g. DOIs, PubMed IDs, ISBNs, URLs), whose metadata is then retrieved and converted to a user-specified style. Manubot modernizes publishing to align with the ideals of open science by making it transparent, reproducible, immediate, versioned, collaborative, and free of charge., Competing Interests: The authors have declared that no competing interests exist.

Rhabdomyosarcoma (RMS) is an aggressive pediatric cancer composed of myoblast-like cells. Recently, we discovered a unique muscle progenitor marked by the expression of the Twist2 transcription factor. Genomic analyses of 258 RMS patient tumors uncovered prevalent copy number amplification events and increased expression of TWIST2 in fusion-negative RMS. Knockdown of TWIST2 in RMS cells results in up-regulation of MYOGENIN and a decrease in proliferation, implicating TWIST2 as an oncogene in RMS. Through an inducible Twist2 expression system, we identified Twist2 as a reversible inhibitor of myogenic differentiation with the remarkable ability to promote myotube dedifferentiation in vitro. Integrated analysis of genome-wide ChIP-seq and RNA-seq data revealed the first dynamic chromatin and transcriptional landscape of Twist2 binding during myogenic differentiation. During differentiation, Twist2 competes with MyoD at shared DNA motifs to direct global gene transcription and repression of the myogenic program. Additionally, Twist2 shapes the epigenetic landscape to drive chromatin opening at oncogenic loci and chromatin closing at myogenic loci. These epigenetic changes redirect MyoD binding from myogenic genes toward oncogenic, metabolic, and growth genes. Our study reveals the dynamic interplay between two opposing transcriptional regulators that control the fate of RMS and provides insight into the molecular etiology of this aggressive form of cancer., (© 2019 Li et al.; Published by Cold Spring Harbor Laboratory Press.)

Cancer evolves through a multistep process that occurs by the temporal accumulation of genetic mutations. Tumor-derived exosomes are emerging contributors to tumorigenesis. To understand how exosomes might contribute to cell transformation, we utilized the classic two-step NIH/3T3 cell transformation assay and observed that exosomes isolated from pancreatic cancer cells, but not normal human cells, can initiate malignant cell transformation and these transformed cells formed tumors in vivo. However, cancer cell exosomes are unable to transform cells alone or to act as a promoter of cell transformation. Utilizing proteomics and exome sequencing, we discovered cancer cell exosomes act as an initiator by inducing random mutations in recipient cells. Cells from the pool of randomly mutated cells are driven to transformation by a classic promoter resulting in foci, each of which encode a unique genetic profile. Our studies describe a novel molecular understanding of how cancer cell exosomes contribute to cell transformation., Editorial Note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that major issues remain unresolved (see decision letter)., Competing Interests: KS, KS, Md, HG, JT, SC, MK, VM, RB No competing interests declared, KO Reviewing editor, eLife, (© 2019, Stefanius et al.)

Background: Inorganic phosphate (Pi) is used extensively as a preservative and a flavor enhancer in the Western diet. Physical inactivity, a common feature of Western societies, is associated with increased cardiovascular morbidity and mortality. It is unknown whether dietary Pi excess contributes to exercise intolerance and physical inactivity., Methods: To determine an association between Pi excess and physical activity in humans, we assessed the relationship between serum Pi and actigraphy-determined physical activity level, as well as left ventricular function by cardiac magnetic resonance imaging, in DHS-2 (Dallas Heart Study phase 2) participants after adjusting for relevant variables. To determine direct effects of dietary Pi on exercise capacity, oxygen uptake, serum nonesterified fatty acid, and glucose were measured during exercise treadmill test in C57/BL6 mice fed either a high-Pi (2%) or normal-Pi (0.6%) diet for 12 weeks. To determine the direct effect of Pi on muscle metabolism and expression of genes involved in fatty acid metabolism, additional studies in differentiated C2C12 myotubes were conducted after subjecting to media containing 1 to 3 mmol/L Pi (pH 7.0) to simulate in vivo phosphate conditions., Results: In participants of the DHS-2 (n=1603), higher serum Pi was independently associated with reduced time spent in moderate to vigorous physical activity ( P=0.01) and increased sedentary time ( P=0.004). There was no association between serum Pi and left ventricular ejection fraction or volumes. In animal studies, compared with the control diet, consumption of high-Pi diet for 12 weeks did not alter body weight or left ventricular function but reduced maximal oxygen uptake, treadmill duration, spontaneous locomotor activity, fat oxidation, and fatty acid levels and led to downregulation of genes involved in fatty acid synthesis, release, and oxidation, including Fabp4, Hsl, Fasn, and Pparγ, in muscle. Similar results were recapitulated in vitro by incubating C2C12 myotubes with high-Pi media., Conclusions: Our data demonstrate a detrimental effect of dietary Pi excess on skeletal muscle fatty acid metabolism and exercise capacity that is independent of obesity and cardiac contractile function. Dietary Pi may represent a novel and modifiable target to reduce physical inactivity associated with the Western diet.

Therapeutic blockade of the PD-1/PD-L1 axis is recognized as an effective treatment for numerous cancer types. However, only a subset of patients respond to this treatment, warranting a greater understanding of the biological mechanisms driving immune evasion via PD-1/PD-L1 signaling and other T-cell suppressive pathways. We previously identified a head and neck squamous cell carcinoma with human papillomavirus integration in the PD-L1 locus upstream of the transmembrane domain-encoding region, suggesting expression of a truncated form of PD-L1 (Parfenov et al., Proc Natl Acad Sci USA 111(43):15544-15549, 2014). In this study, we extended this observation by performing a computational analysis of 33 other cancer types as well as human cancer cell lines, and identified additional PD-L1 isoforms with an exon 4 enrichment expressed in 20 cancers and human cancer cell lines. We demonstrate that cancer cell lines with high expression levels of exon 4-enriched PD-L1 generate a secreted form of PD-L1. Further biochemical studies of exon 4-enriched PD-L1 demonstrated that this form is secreted and maintains the capacity to bind PD-1 as well as to serve as a negative regulator on T cell function, as measured by inhibition of IL-2 and IFNg secretion. Overall, we have demonstrated that truncated forms of PD-L1 exist in numerous cancer types, and have validated that truncated PD-L1 can be secreted and negatively regulate T cell function.

Background: Tumor orchestrated metabolic changes in the microenvironment limit generation of anti-tumor immune responses. Availability of arginine, a semi-essential amino acid, is critical for lymphocyte proliferation and function. Levels of arginine are regulated by the enzymes arginase 1,2 and nitric oxide synthase (NOS). However, the role of arginase activity in lung tumor maintenance has not been investigated in clinically relevant orthotopic tumor models., Methods: RNA sequencing (RNA-seq) of sorted cell populations from mouse lung adenocarcinomas derived from immunocompetent genetically engineered mouse models (GEMM)s was performed. To complement mouse studies, a patient tissue microarray consisting of 150 lung adenocarcinomas, 103 squamous tumors, and 54 matched normal tissue were stained for arginase, CD3, and CD66b by multiplex immunohistochemistry. Efficacy of a novel arginase inhibitor compound 9 in reversing arginase mediated T cell suppression was determined in splenocyte ex vivo assays. Additionally, the anti-tumor activity of this compound was determined in vitro and in an autochthonous immunocompetent Kras G12D GEMM of lung adenocarcinoma model., Results: Analysis of RNA-seq of sorted myeloid cells suggested that arginase expression is elevated in myeloid cells in the tumor as compared to the normal lung tissue. Accordingly, in the patient samples arginase 1 expression was mainly localized in the granulocytic myeloid cells and significantly elevated in both lung adenocarcinoma and squamous tumors as compared to the controls. Our ex vivo analysis demonstrated that myeloid derived suppressor cell (MDSC)s cause T cell suppression by arginine depletion, and suppression of arginase activity by a novel ARG1/2 inhibitor, compound 9, led to restoration of T cell function by increasing arginine. Treatment of Kras G12D GEMM of lung cancer model with compound 9 led to a significant tumor regression associated with increased T cell numbers and function, while it had no activity across several murine and human non-small cell (NSCLC) lung cancer lines in vitro., Conclusions: We show that arginase expression is elevated in mouse and patient lung tumors. In a KRAS G12D GEMM arginase inhibition diminished growth of established tumors. Our data suggest arginase as an immunomodulatory target that should further be investigated in lung tumors with high arginase activity.

A comprehensive cellular anatomy of normal human prostate is essential for solving the cellular origins of benign prostatic hyperplasia and prostate cancer. The tools used to analyze the contribution of individual cell types are not robust. We provide a cellular atlas of the young adult human prostate and prostatic urethra using an iterative process of single-cell RNA sequencing (scRNA-seq) and flow cytometry on ∼98,000 cells taken from different anatomical regions. Immunohistochemistry with newly derived cell type-specific markers revealed the distribution of each epithelial and stromal cell type on whole mounts, revising our understanding of zonal anatomy. Based on discovered cell surface markers, flow cytometry antibody panels were designed to improve the purification of each cell type, with each gate confirmed by scRNA-seq. The molecular classification, anatomical distribution, and purification tools for each cell type in the human prostate create a powerful resource for experimental design in human prostate disease., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

Twist transcription factors function as ancestral regulators of mesodermal cell fates in organisms ranging from Drosophila to mammals. Through lineage tracing of Twist2 (Tw2)-expressing cells with tamoxifen-inducible Tw2-CreERT2 and tdTomato (tdTO) reporter mice, we discovered a unique cell population that progressively contributes to cardiomyocytes (CMs), endothelial cells, and fibroblasts in the adult heart. Clonal analysis confirmed the ability of Tw2-derived tdTO + (Tw2-tdTO + ) cells to form CMs in vitro. Within the adult heart, Tw2-tdTO + CMs accounted for ∼13% of total CMs, the majority of which resulted from fusion of Tw2-tdTO + cells with existing CMs. Tw2-tdTO + cells also contribute to cardiac remodeling after injury. We conclude that Tw2-tdTO + cells participate in lifelong maintenance of cardiac function, at least in part through de novo formation of CMs and fusion with preexisting CMs, as well as in the genesis of other cellular components of the adult heart., Competing Interests: The authors declare no conflict of interest.

Background: Epigenetic regulators are frequently mutated or aberrantly expressed in a variety of cancers, leading to altered transcription states that result in changes in cell identity, behavior, and response to therapy., Results: To define alterations in epigenetic landscapes in breast cancers, we profiled the distributions of 8 key histone modifications by ChIP-Seq, as well as primary (GRO-seq) and steady state (RNA-Seq) transcriptomes, across 13 distinct cell lines that represent 5 molecular subtypes of breast cancer and immortalized human mammary epithelial cells., Discussion: Using combinatorial patterns of distinct histone modification signals, we defined subtype-specific chromatin signatures to nominate potential biomarkers. This approach identified AFAP1-AS1 as a triple negative breast cancer-specific gene associated with cell proliferation and epithelial-mesenchymal-transition. In addition, our chromatin mapping data in basal TNBC cell lines are consistent with gene expression patterns in TCGA that indicate decreased activity of the androgen receptor pathway but increased activity of the vitamin D biosynthesis pathway., Conclusions: Together, these datasets provide a comprehensive resource for histone modification profiles that define epigenetic landscapes and reveal key chromatin signatures in breast cancer cell line subtypes with potential to identify novel and actionable targets for treatment.

Noncoding transcription is a defining feature of active enhancers, linking transcription factor (TF) binding to the molecular mechanisms controlling gene expression. To determine the relationship between enhancer activity and biological outcomes in breast cancers, we profiled the transcriptomes (using GRO-seq and RNA-seq) and epigenomes (using ChIP-seq) of 11 different human breast cancer cell lines representing five major molecular subtypes of breast cancer, as well as two immortalized ("normal") human breast cell lines. In addition, we developed a robust and unbiased computational pipeline that simultaneously identifies putative subtype-specific enhancers and their cognate TFs by integrating the magnitude of enhancer transcription, TF mRNA expression levels, TF motif P -values, and enrichment of H3K4me1 and H3K27ac. When applied across the 13 different cell lines noted above, the Total Functional Score of Enhancer Elements (TFSEE) identified key breast cancer subtype-specific TFs that act at transcribed enhancers to dictate gene expression patterns determining growth outcomes, including Forkhead TFs, FOSL1, and PLAG1. FOSL1, a Fos family TF, (1) is highly enriched at the enhancers of triple negative breast cancer (TNBC) cells, (2) acts as a key regulator of the proliferation and viability of TNBC cells, but not Luminal A cells, and (3) is associated with a poor prognosis in TNBC breast cancer patients. Taken together, our results validate our enhancer identification pipeline and reveal that enhancers transcribed in breast cancer cells direct critical gene regulatory networks that promote pathogenesis., (© 2018 Franco et al.; Published by Cold Spring Harbor Laboratory Press.)

The pluripotency of embryonic stem cells (ESCs) relies on appropriate responsiveness to developmental cues. Promoter-proximal pausing of RNA polymerase II (Pol II) has been suggested to play a role in keeping genes poised for future activation. To identify the role of Pol II pausing in regulating ESC pluripotency, we have generated mouse ESCs carrying a mutation in the pause-inducing factor SPT5. Genomic studies reveal genome-wide reduction of paused Pol II caused by mutant SPT5 and further identify a tight correlation between pausing-mediated transcription effect and local chromatin environment. Functionally, this pausing-deficient SPT5 disrupts ESC differentiation upon removal of self-renewal signals. Thus, our study uncovers an important role of Pol II pausing in regulating ESC differentiation and suggests a model that Pol II pausing coordinates with epigenetic modification to influence transcription during mESC differentiation., (Copyright © 2017 The Authors. Published by Elsevier B.V. All rights reserved.)

The NELF complex is a metazoan-specific factor essential for establishing transcription pausing. Although NELF has been implicated in cell fate regulation, the cellular regulation of NELF and its intrinsic role in specific lineage differentiation remains largely unknown. Using mammalian hematopoietic differentiation as a model system, here we identified a dynamic change of NELF-mediated transcription pausing as a novel mechanism regulating hematopoietic differentiation. We found a sharp decrease of NELF protein abundance upon granulocytic differentiation and a subsequent genome-wide reduction of transcription pausing. This loss of pausing coincides with activation of granulocyte-affiliated genes and diminished expression of progenitor markers. Functional studies revealed that sustained expression of NELF inhibits granulocytic differentiation, whereas NELF depletion in progenitor cells leads to premature differentiation towards the granulocytic lineage. Our results thus uncover a previously unrecognized regulation of transcription pausing by modulating NELF protein abundance to control cellular differentiation., Competing Interests: Conflict-of-interest disclosure: the authors declare no competing financial interests.

The Encyclopedia of DNA elements (ENCODE) project is an ongoing collaborative effort to create a comprehensive catalog of functional elements initiated shortly after the completion of the Human Genome Project. The current database exceeds 6500 experiments across more than 450 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the H. sapiens and M. musculus genomes. All ENCODE experimental data, metadata, and associated computational analyses are submitted to the ENCODE Data Coordination Center (DCC) for validation, tracking, storage, unified processing, and distribution to community resources and the scientific community. As the volume of data increases, the identification and organization of experimental details becomes increasingly intricate and demands careful curation. The ENCODE DCC has created a general purpose software system, known as SnoVault, that supports metadata and file submission, a database used for metadata storage, web pages for displaying the metadata and a robust API for querying the metadata. The software is fully open-source, code and installation instructions can be found at: http://github.com/ENCODE-DCC/snovault/ (for the generic database) and http://github.com/ENCODE-DCC/encoded/ to store genomic data in the manner of ENCODE. The core database engine, SnoVault (which is completely independent of ENCODE, genomic data, or bioinformatic data) has been released as a separate Python package.

Transcriptional enhancers are DNA regulatory elements that are bound by transcription factors and act to positively regulate the expression of nearby or distally located target genes. Enhancers have many features that have been discovered using genomic analyses. Recent studies have shown that active enhancers recruit RNA polymerase II (Pol II) and are transcribed, producing enhancer RNAs (eRNAs). GRO-seq, a method for identifying the location and orientation of all actively transcribing RNA polymerases across the genome, is a powerful approach for monitoring nascent enhancer transcription. Furthermore, the unique pattern of enhancer transcription can be used to identify enhancers in the absence of any information about the underlying transcription factors. Here, we describe the computational approaches required to identify and analyze active enhancers using GRO-seq data, including data pre-processing, alignment, and transcript calling. In addition, we describe protocols and computational pipelines for mining GRO-seq data to identify active enhancers, as well as known transcription factor binding sites that are transcribed. Furthermore, we discuss approaches for integrating GRO-seq-based enhancer data with other genomic data, including target gene expression and function. Finally, we describe molecular biology assays that can be used to confirm and explore further the function of enhancers that have been identified using genomic assays. Together, these approaches should allow the user to identify and explore the features and biological functions of new cell type-specific enhancers.

The Encyclopedia of DNA Elements (ENCODE) Data Coordinating Center (DCC) is responsible for organizing, describing and providing access to the diverse data generated by the ENCODE project. The description of these data, known as metadata, includes the biological sample used as input, the protocols and assays performed on these samples, the data files generated from the results and the computational methods used to analyze the data. Here, we outline the principles and philosophy used to define the ENCODE metadata in order to create a metadata standard that can be applied to diverse assays and multiple genomic projects. In addition, we present how the data are validated and used by the ENCODE DCC in creating the ENCODE Portal (https://www.encodeproject.org/). Database URL: www.encodeproject.org., (© The Author(s) 2016. Published by Oxford University Press.)

The Encyclopedia of DNA Elements (ENCODE) Project is in its third phase of creating a comprehensive catalog of functional elements in the human genome. This phase of the project includes an expansion of assays that measure diverse RNA populations, identify proteins that interact with RNA and DNA, probe regions of DNA hypersensitivity, and measure levels of DNA methylation in a wide range of cell and tissue types to identify putative regulatory elements. To date, results for almost 5000 experiments have been released for use by the scientific community. These data are available for searching, visualization and download at the new ENCODE Portal (www.encodeproject.org). The revamped ENCODE Portal provides new ways to browse and search the ENCODE data based on the metadata that describe the assays as well as summaries of the assays that focus on data provenance. In addition, it is a flexible platform that allows integration of genomic data from multiple projects. The portal experience was designed to improve access to ENCODE data by relying on metadata that allow reusability and reproducibility of the experiments., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The Encyclopedia of DNA elements (ENCODE) project is an ongoing collaborative effort to create a catalog of genomic annotations. To date, the project has generated over 4000 experiments across more than 350 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory network and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All ENCODE experimental data, metadata and associated computational analyses are submitted to the ENCODE Data Coordination Center (DCC) for validation, tracking, storage and distribution to community resources and the scientific community. As the volume of data increases, the organization of experimental details becomes increasingly complicated and demands careful curation to identify related experiments. Here, we describe the ENCODE DCC's use of ontologies to standardize experimental metadata. We discuss how ontologies, when used to annotate metadata, provide improved searching capabilities and facilitate the ability to find connections within a set of experiments. Additionally, we provide examples of how ontologies are used to annotate ENCODE metadata and how the annotations can be identified via ontology-driven searches at the ENCODE portal. As genomic datasets grow larger and more interconnected, standardization of metadata becomes increasingly vital to allow for exploration and comparison of data between different scientific projects., (© The Author(s) 2015. Published by Oxford University Press.)

The laboratory mouse shares the majority of its protein-coding genes with humans, making it the premier model organism in biomedical research, yet the two mammals differ in significant ways. To gain greater insights into both shared and species-specific transcriptional and cellular regulatory programs in the mouse, the Mouse ENCODE Consortium has mapped transcription, DNase I hypersensitivity, transcription factor binding, chromatin modifications and replication domains throughout the mouse genome in diverse cell and tissue types. By comparing with the human genome, we not only confirm substantial conservation in the newly annotated potential functional sequences, but also find a large degree of divergence of sequences involved in transcriptional regulation, chromatin state and higher order chromatin organization. Our results illuminate the wide range of evolutionary forces acting on genes and their regulatory regions, and provide a general resource for research into mammalian biology and mechanisms of human diseases.

The University of California Santa Cruz (UCSC) Genome Browser (http://genome.ucsc.edu) offers online public access to a growing database of genomic sequence and annotations for a wide variety of organisms. The Browser is an integrated tool set for visualizing, comparing, analysing and sharing both publicly available and user-generated genomic datasets. As of September 2012, genomic sequence and a basic set of annotation 'tracks' are provided for 63 organisms, including 26 mammals, 13 non-mammal vertebrates, 3 invertebrate deuterostomes, 13 insects, 6 worms, yeast and sea hare. In the past year 19 new genome assemblies have been added, and we anticipate releasing another 28 in early 2013. Further, a large number of annotation tracks have been either added, updated by contributors or remapped to the latest human reference genome. Among these are an updated UCSC Genes track for human and mouse assemblies. We have also introduced several features to improve usability, including new navigation menus. This article provides an update to the UCSC Genome Browser database, which has been previously featured in the Database issue of this journal.

The Encyclopedia of DNA Elements (ENCODE), http://encodeproject.org, has completed its fifth year of scientific collaboration to create a comprehensive catalog of functional elements in the human genome, and its third year of investigations in the mouse genome. Since the last report in this journal, the ENCODE human data repertoire has grown by 898 new experiments (totaling 2886), accompanied by a major integrative analysis. In the mouse genome, results from 404 new experiments became available this year, increasing the total to 583, collected during the course of the project. The University of California, Santa Cruz, makes this data available on the public Genome Browser http://genome.ucsc.edu for visual browsing and data mining. Download of raw and processed data files are all supported. The ENCODE portal provides specialized tools and information about the ENCODE data sets.

To complement the human Encyclopedia of DNA Elements (ENCODE) project and to enable a broad range of mouse genomics efforts, the Mouse ENCODE Consortium is applying the same experimental pipelines developed for human ENCODE to annotate the mouse genome.

The Encyclopedia of DNA Elements (ENCODE) Consortium is entering its 5th year of production-level effort generating high-quality whole-genome functional annotations of the human genome. The past year has brought the ENCODE compendium of functional elements to critical mass, with a diverse set of 27 biochemical assays now covering 200 distinct human cell types. Within the mouse genome, which has been under study by ENCODE groups for the past 2 years, 37 cell types have been assayed. Over 2000 individual experiments have been completed and submitted to the Data Coordination Center for public use. UCSC makes this data available on the quality-reviewed public Genome Browser (http://genome.ucsc.edu) and on an early-access Preview Browser (http://genome-preview.ucsc.edu). Visual browsing, data mining and download of raw and processed data files are all supported. An ENCODE portal (http://encodeproject.org) provides specialized tools and information about the ENCODE data sets.

The University of California Santa Cruz Genome Browser (http://genome.ucsc.edu) offers online public access to a growing database of genomic sequence and annotations for a wide variety of organisms. The Browser is an integrated tool set for visualizing, comparing, analyzing and sharing both publicly available and user-generated genomic data sets. In the past year, the local database has been updated with four new species assemblies, and we anticipate another four will be released by the end of 2011. Further, a large number of annotation tracks have been either added, updated by contributors, or remapped to the latest human reference genome. Among these are new phenotype and disease annotations, UCSC genes, and a major dbSNP update, which required new visualization methods. Growing beyond the local database, this year we have introduced 'track data hubs', which allow the Genome Browser to provide access to remotely located sets of annotations. This feature is designed to significantly extend the number and variety of annotation tracks that are publicly available for visualization and analysis from within our site. We have also introduced several usability features including track search and a context-sensitive menu of options available with a right-click anywhere on the Browser's image.

This study was designed to address issues regarding sample size and marker location that have arisen from the discovery of SNPs in the genomes of poorly characterized primate species and the application of these markers to the study of primate population genetics. We predict the effect of discovery sample size on the probability of discovering both rare and common SNPs and then compare this prediction with the proportion of common and rare SNPs discovered when different numbers of individuals are sequenced. Second, we examine the effect of genomic region on estimates of common population genetic data, comparing markers from both coding and non-coding regions of the rhesus macaque genome and the population genetic data calculated from these markers, to measure the degree and direction of bias introduced by SNPs located in coding versus non-coding regions of the genome. We found that both discovery sample size and genomic region surveyed affect SNP marker attributes and population genetic estimates, even when these are calculated from an expanded data set containing more individuals than the original discovery data set. Although none of the SNP detection methods or genomic regions tested in this study was completely uninformative, these results show that each has a different kind of genetic variation that is suitable for different purposes, and each introduces specific types of bias. Given that each SNP marker has an individual evolutionary history, we calculated that the most complete and unbiased representation of the genetic diversity present in the individual can be obtained by incorporating at least 10 individuals into the discovery sample set, to ensure the discovery of both common and rare polymorphisms.

The ENCODE project is an international consortium with a goal of cataloguing all the functional elements in the human genome. The ENCODE Data Coordination Center (DCC) at the University of California, Santa Cruz serves as the central repository for ENCODE data. In this role, the DCC offers a collection of high-throughput, genome-wide data generated with technologies such as ChIP-Seq, RNA-Seq, DNA digestion and others. This data helps illuminate transcription factor-binding sites, histone marks, chromatin accessibility, DNA methylation, RNA expression, RNA binding and other cell-state indicators. It includes sequences with quality scores, alignments, signals calculated from the alignments, and in most cases, element or peak calls calculated from the signal data. Each data set is available for visualization and download via the UCSC Genome Browser (http://genome.ucsc.edu/). ENCODE data can also be retrieved using a metadata system that captures the experimental parameters of each assay. The ENCODE web portal at UCSC (http://encodeproject.org/) provides information about the ENCODE data and links for access.

A 60 bp sequence variation hotspot in the canine mitochondrial DNA hypervariable region 1 was evaluated for its use in forensic investigations. Nineteen haplotypes containing 18 single nucleotide polymorphisms were observed among laboratory-generated and GenBank-derived domestic dog sequences representing five regional localities in the U.S. Samples from the different localities were highly variable with the levels of intra-population variability being similar among the populations studied. AMOVA further confirmed that there was no significant genetic structuring of the populations. Assays using these haplotypes were robust, canid specific and portend a rapid method for correctly excluding individual dogs as noncontributors of forensic evidence. Species-specificity of the primers was confirmed by means of in-tube polymerase chain reaction of human and cat DNA and in-silico assessment of the genomes of several animal species. Breed-specific fragments were not detected among the common haplotypes but there is evidence that this assay may be capable of differentiating domestic dog, wolf, and coyote sequences.

Background: Rhesus macaques (Macaca mulatta) are the primate most used for biomedical research, but phenotypic differences between Indian-origin and Chinese rhesus macaques have encouraged genetic methods for identifying genetic differences between these two populations. The completion of the rhesus genome has led to the identification of many single nucleotide polymorphisms (SNPs) in this species. These single nucleotide polymorphisms have many advantages over the short tandem repeat (STR) loci currently used to assay genetic variation. However, the number of currently identified polymorphisms is too small for whole genome analysis or studies of quantitative trait loci. To that end, we tested a combination of methods to identify large numbers of high-confidence SNPs, and screen those with high minor allele frequencies (MAF)., Results: By testing our previously reported single nucleotide polymorphisms, we identified a subset of high-confidence, high-MAF polymorphisms. Resequencing revealed a large number of regionally specific SNPs not identified through a single pyrosequencing run. By resequencing a pooled sample of four individuals, we reliably identified loci with a MAF of at least 12.5%. Finally, we found that when applied to a larger, geographically variable sample of rhesus, a large proportion of our loci were variable in both populations, and very few loci were ancestry informative. Despite this fact, the SNP loci were more effective at discriminating Indian and Chinese rhesus than STR loci., Conclusion: Pyrosequencing and pooled resequencing are viable methods for the identification of high-MAF SNP loci in rhesus macaques. These SNP loci are appropriate for screening both the inter- and intra-population genetic variation.

The 608-bp hypervariable region 1 (HV1) sequences from 36 local dogs were analyzed to characterize the population genetic structure of canid mitochondrial DNA (mtDNA). Sixteen haplotypes were identified. A 417-bp segment of this sequence was compared with GenBank sequences from a geographically representative sample of 201 dogs, two coyotes, and two wolves. Sixty-six haplotypes were identified including 62 found only in domestic dogs. Fourteen of these correspond to the 16 local haplotypes and were among the most frequent haplotypes. The local sample was judged to be representative of the much broader geographic sample. No correlation was observed between local haplotypes and the owner's characterization of dog breed. A 60-bp variation "hotspot" within the canid HV1 was identified as a potentially valuable molecular tool, particularly for assaying limited or degraded DNA samples.