Gene Promoter Evolution Targets the Center of the Human Protein Interaction Network

Assessing the contribution of promoters and coding sequences to gene evolution is an important step toward discovering the major genetic determinants of human evolution. Many specific examples have revealed the evolutionary importance of cis-regulatory regions. However, the relative contribution of regulatory and coding regions to the evolutionary process and whether systemic factors differentially influence their evolution remains unclear. To address these questions, we carried out an analysis at the genome scale to identify signatures of positive selection in human proximal promoters. Next, we examined whether genes with positively selected promoters (Prom+ genes) show systemic differences with respect to a set of genes with positively selected protein-coding regions (Cod+ genes). We found that the number of genes in each set was not significantly different (8.1% and 8.5%, respectively). Furthermore, a functional analysis showed that, in both cases, positive selection affects almost all biological processes and only a few genes of each group are located in enriched categories, indicating that promoters and coding regions are not evolutionarily specialized with respect to gene function. On the other hand, we show that the topology of the human protein network has a different influence on the molecular evolution of proximal promoters and coding regions. Notably, Prom+ genes have an unexpectedly high centrality when compared with a reference distribution (P = 0.008, for Eigenvalue centrality). Moreover, the frequency of Prom+ genes increases from the periphery to the center of the protein network (P = 0.02, for the logistic regression coefficient). This means that gene centrality does not constrain the evolution of proximal promoters, unlike the case with coding regions, and further indicates that the evolution of proximal promoters is more efficient in the center of the protein network than in the periphery. These results show that proximal promoters have had a systemic contribution to human evolution by increasing the participation of central genes in the evolutionary process

Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smokin

Pan-Cancer Analysis of lncRNA Regulation Supports Their Targeting of Cancer Genes in Each Tumor Context

Long noncoding RNAs (lncRNAs) are commonly dys-regulated in tumors, but only a handful are known toplay pathophysiological roles in cancer. We inferredlncRNAs that dysregulate cancer pathways, onco-genes, and tumor suppressors (cancer genes) bymodeling their effects on the activity of transcriptionfactors, RNA-binding proteins, and microRNAs in5,185 TCGA tumors and 1,019 ENCODE assays.Our predictions included hundreds of candidateonco- and tumor-suppressor lncRNAs (cancerlncRNAs) whose somatic alterations account for thedysregulation of dozens of cancer genes and path-ways in each of 14 tumor contexts. To demonstrateproof of concept, we showed that perturbations tar-geting OIP5-AS1 (an inferred tumor suppressor) andTUG1 and WT1-AS (inferred onco-lncRNAs) dysre-gulated cancer genes and altered proliferation ofbreast and gynecologic cancer cells. Our analysis in-dicates that, although most lncRNAs are dysregu-lated in a tumor-specific manner, some, includingOIP5-AS1, TUG1, NEAT1, MEG3, and TSIX, synergis-tically dysregulate cancer pathways in multiple tumorcontexts

Spatial Organization and Molecular Correlation of Tumor-Infiltrating Lymphocytes Using Deep Learning on Pathology Images

Beyond sample curation and basic pathologic characterization, the digitized H&E-stained images of TCGA samples remain underutilized. To highlight this resource, we present mappings of tumorinfiltrating lymphocytes (TILs) based on H&E images from 13 TCGA tumor types. These TIL maps are derived through computational staining using a convolutional neural network trained to classify patches of images. Affinity propagation revealed local spatial structure in TIL patterns and correlation with overall survival. TIL map structural patterns were grouped using standard histopathological parameters. These patterns are enriched in particular T cell subpopulations derived from molecular measures. TIL densities and spatial structure were differentially enriched among tumor types, immune subtypes, and tumor molecular subtypes, implying that spatial infiltrate state could reflect particular tumor cell aberration states. Obtaining spatial lymphocytic patterns linked to the rich genomic characterization of TCGA samples demonstrates one use for the TCGA image archives with insights into the tumor-immune microenvironment

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although theMYConcogene has been implicated incancer, a systematic assessment of alterations ofMYC, related transcription factors, and co-regulatoryproteins, forming the proximal MYC network (PMN),across human cancers is lacking. Using computa-tional approaches, we define genomic and proteo-mic features associated with MYC and the PMNacross the 33 cancers of The Cancer Genome Atlas.Pan-cancer, 28% of all samples had at least one ofthe MYC paralogs amplified. In contrast, the MYCantagonists MGA and MNT were the most frequentlymutated or deleted members, proposing a roleas tumor suppressors.MYCalterations were mutu-ally exclusive withPIK3CA,PTEN,APC,orBRAFalterations, suggesting that MYC is a distinct onco-genic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such asimmune response and growth factor signaling; chro-matin, translation, and DNA replication/repair wereconserved pan-cancer. This analysis reveals insightsinto MYC biology and is a reference for biomarkersand therapeutics for cancers with alterations ofMYC or the PMN

SHIV-162P3 Infection of Rhesus Macaques Given Maraviroc Gel Vaginally Does Not Involve Resistant Viruses

Maraviroc (MVC) gels are effective at protecting rhesus macaques from vaginal SHIV transmission, but breakthrough infections can occur. To determine the effects of a vaginal MVC gel on infecting SHIV populations in a macaque model, we analyzed plasma samples from three rhesus macaques that received a MVC vaginal gel (day 0) but became infected after high-dose SHIV-162P3 vaginal challenge. Two infected macaques that received a placebo gel served as controls. The infecting SHIV-162P3 stock had an overall mean genetic distance of 0.294±0.027%; limited entropy changes were noted across the envelope (gp160). No envelope mutations were observed consistently in viruses isolated from infected macaques at days 14–21, the time of first detectable viremia, nor selected at later time points, days 42–70. No statistically significant differences in MVC susceptibilities were observed between the SHIV inoculum (50% inhibitory concentration [IC50] 1.87 nM) and virus isolated from the three MVC-treated macaques (MVC IC50 1.18 nM, 1.69 nM, and 1.53 nM, respectively). Highlighter plot analyses suggested that infection was established in each MVC-treated animal by one founder virus genotype. The expected Poisson distribution of pairwise Hamming Distance frequency counts was observed and a phylogenetic analysis did not identify infections with distinct lineages from the challenge stock. These data suggest that breakthrough infections most likely result from incomplete viral inhibition and not the selection of MVC-resistant variants

Next-generation mitogenomics: A comparison of approaches applied to caecilian amphibian phylogeny

Mitochondrial genome (mitogenome) sequences are being generated with increasing speed due to the advances of next-generation sequencing (NGS) technology and associated analytical tools. However, detailed comparisons to explore the utility of alternative NGS approaches applied to the same taxa have not been undertaken. We compared a 'traditional' Sanger sequencing method with two NGS approaches (shotgun sequencing and non-indexed, multiplex amplicon sequencing) on four different sequencing platforms (Illumina's HiSeq and MiSeq, Roche's 454 GS FLX, and Life Technologies' Ion Torrent) to produce seven (near-) complete mitogenomes from six species that form a small radiation of caecilian amphibians from the Seychelles. The fastest, most accurate method of obtaining mitogenome sequences that we tested was direct sequencing of genomic DNA (shotgun sequencing) using the MiSeq platform. Bayesian inference and maximum likelihood analyses using seven different partitioning strategies were unable to resolve compellingly all phylogenetic relationships among the Seychelles caecilian species, indicating the need for additional data in this case

The Felix-trial. Double-blind randomization of interspinous implant or bony decompression for treatment of spinal stenosis related intermittent neurogenic claudication

Abstract. Background. Decompressive laminotomy is the standard surgical procedure in the treatment of patients with canal stenosis related intermittent neurogenic claudication. New techniques, such as interspinous process implants, claim a shorter hospital stay, less post-operative pain and equal long-term functional outcome. A comparative (cost-) effectiveness study has not been performed yet. This protocol describes the design of a randomized controlled trial (RCT) on (cost-) effectiveness of the use of interspinous process implants versus conventi

Biallelic and monoallelic ESR2 variants associated with 46,XY disorders of sex development

Purpose: Disorders or differences of sex development (DSDs) are rare congenital conditions characterized by atypical sex development. Despite advances in genomic technologies, the molecular cause remains unknown in 50% of cases. Methods: Homozygosity mapping and whole-exome sequencing revealed an ESR2 variant in an individual with syndromic 46, XY DSD. Additional cases with 46, XY DSD underwent whole-exome sequencing and targeted next-generation sequencing of ESR2. Functional characterization of the identified variants included luciferase assays and protein structure analysis. Gonadal ESR2 expression was assessed in human embryonic data sets and immunostaining of estrogen receptor-beta (ER-beta) was performed in an 8-week-old human male embryo. Results: We identified a homozygous ESR2 variant, c.541_543del p. (Asn181del), located in the highly conserved DNA-binding domain of ER-beta, in an individual with syndromic 46, XY DSD. Two additional heterozygous missense variants, c.251G>T p.(Gly84Val) and c.1277T>G p.(Leu426Arg), located in the N-terminus and the ligand-binding domain of ER-beta, were found in unrelated, nonsyndromic 46, XY DSD cases. Significantly increased transcriptional activation and an impact on protein conformation were shown for the p.(Asn181del) and p.(Leu426Arg) variants. Testicular ESR2 expression was previously documented and ER-beta immunostaining was positive in the developing intestine and eyes. Conclusion: Our study supports a role for ESR2 as a novel candidate gene for 46, XY DSD