An Unbiased Molecular Approach Using 3'-UTRs Resolves the Avian Family-Level Tree of Life.

Presumably, due to a rapid early diversification, major parts of the higher-level phylogeny of birds are still resolved controversially in different analyses or are considered unresolvable. To address this problem, we produced an avian tree of life, which includes molecular sequences of one or several species of ∼90% of the currently recognized family-level taxa (429 species, 379 genera) including all 106 family-level taxa of the nonpasserines and 115 of the passerines (Passeriformes). The unconstrained analyses of noncoding 3-prime untranslated region (3'-UTR) sequences and those of coding sequences yielded different trees. In contrast to the coding sequences, the 3'-UTR sequences resulted in a well-resolved and stable tree topology. The 3'-UTR contained, unexpectedly, transcription factor binding motifs that were specific for different higher-level taxa. In this tree, grebes and flamingos are the sister clade of all other Neoaves, which are subdivided into five major clades. All nonpasserine taxa were placed with robust statistical support including the long-time enigmatic hoatzin (Opisthocomiformes), which was found being the sister taxon of the Caprimulgiformes. The comparatively late radiation of family-level clades of the songbirds (oscine Passeriformes) contrasts with the attenuated diversification of nonpasseriform taxa since the early Miocene. This correlates with the evolution of vocal production learning, an important speciation factor, which is ancestral for songbirds and evolved convergent only in hummingbirds and parrots. As 3'-UTR-based phylotranscriptomics resolved the avian family-level tree of life, we suggest that this procedure will also resolve the all-species avian tree of life.

Bridge helix arginines play a critical role in Cas9 sensitivity to mismatches.

The RNA-programmable DNA-endonuclease Cas9 is widely used for genome engineering, where a high degree of specificity is required. To investigate which features of Cas9 determine the sensitivity to mismatches along the target DNA, we performed in vitro biochemical assays and bacterial survival assays in Escherichia coli. We demonstrate that arginines in the Cas9 bridge helix influence guide RNA, and target DNA binding and cleavage. They cluster in two groups that either increase or decrease the Cas9 sensitivity to mismatches. We show that the bridge helix is essential for R-loop formation and that R63 and R66 reduce Cas9 specificity by stabilizing the R-loop in the presence of mismatches. Additionally, we identify Q768 that reduces sensitivity of Cas9 to protospacer adjacent motif-distal mismatches. The Cas9_R63A/Q768A variant showed increased specificity in human cells. Our results provide a firm basis for function- and structure-guided mutagenesis to increase Cas9 specificity for genome engineering.

Genes acting in synapses and neuron projections are early targets of selection during urban colonization.

When a species colonizes an urban habitat, differences in the environment can create novel selection pressures. Successful colonization will further lead to demographic perturbations and genetic drift, which can interfere with selection. Here, we test for consistent urban selection signals in multiple populations of the burrowing owl (Athene cunicularia), a species that colonized South American cities just a few decades ago. We sequenced 213 owls from three urban-rural population pairs and performed a genome-wide comparison of urban against rural birds. We further studied genome-wide associations with flight initiation distance, a measure of harm avoidance in which urban and rural birds are known to differ. Based on four samples taken over nine years from one of the urban populations, we investigated temporal allele frequency changes. The genomic data were also used to identify urban-specific signatures of selective sweeps. Single genomic sites did not reach genome-wide significance for any association. However, a gene-set analysis on the strongest signals from these four selection scans suggests a significant enrichment of genes with known functions related to synapses and neuron projections. We identified 98 genes predominantly expressed in the brain, of which many may play a role in the modulation of brain connectivity and consequently in cognitive function and motivational behaviour during urbanization. Furthermore, polymorphisms in the promoter region of the synaptic SERT gene - one of the two candidates known to correlate with urban colonization in birds - associated with the habitat in which individuals lived (urban vs. rural).

Evolution of genomic variation in the burrowing owl in response to recent colonization of urban areas.

When a species successfully colonizes an urban habitat it can be expected that its population rapidly adapts to the new environment but also experiences demographic perturbations. It is, therefore, essential to gain an understanding of the population structure and the demographic history of the urban and neighbouring rural populations before studying adaptation at the genome level. Here, we investigate populations of the burrowing owl (Athene cunicularia), a species that colonized South American cities just a few decades ago. We assembled a high-quality genome of the burrowing owl and re-sequenced 137 owls from three urban-rural population pairs at 17-fold median sequencing coverage per individual. Our data indicate that each city was independently colonized by a limited number of founders and that restricted gene flow occurred between neighbouring urban and rural populations, but not between urban populations of different cities. Using long-range linkage disequilibrium statistics in an approximate Bayesian computation approach, we estimated consistently lower population sizes in the recent past for the urban populations in comparison to the rural ones. The current urban populations all show reduced standing variation in rare single nucleotide polymorphisms (SNPs), but with different subsets of rare SNPs in different cities. This lowers the potential for local adaptation based on rare variants and makes it harder to detect consistent signals of selection in the genome.

Hsp90 inhibition differentially destabilises MAP kinase and TGF-beta signalling components in cancer cells revealed by kinase-targeted chemoproteomics.

BACKGROUND: The heat shock protein 90 (Hsp90) is required for the stability of many signalling kinases. As a target for cancer therapy it allows the simultaneous inhibition of several signalling pathways. However, its inhibition in healthy cells could also lead to severe side effects. This is the first comprehensive analysis of the response to Hsp90 inhibition at the kinome level. METHODS: We quantitatively profiled the effects of Hsp90 inhibition by geldanamycin on the kinome of one primary (Hs68) and three tumour cell lines (SW480, U2OS, A549) by affinity proteomics based on immobilized broad spectrum kinase inhibitors ("kinobeads"). To identify affected pathways we used the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway classification. We combined Hsp90 and proteasome inhibition to identify Hsp90 substrates in Hs68 and SW480 cells. The mutational status of kinases from the used cell lines was determined using next-generation sequencing. A mutation of Hsp90 candidate client RIPK2 was mapped onto its structure. RESULTS: We measured relative abundances of > 140 protein kinases from the four cell lines in response to geldanamycin treatment and identified many new potential Hsp90 substrates. These kinases represent diverse families and cellular functions, with a strong representation of pathways involved in tumour progression like the BMP, MAPK and TGF-beta signalling cascades. Co-treatment with the proteasome inhibitor MG132 enabled us to classify 64 kinases as true Hsp90 clients. Finally, mutations in 7 kinases correlate with an altered response to Hsp90 inhibition. Structural modelling of the candidate client RIPK2 suggests an impact of the mutation on a proposed Hsp90 binding domain. CONCLUSIONS: We propose a high confidence list of Hsp90 kinase clients, which provides new opportunities for targeted and combinatorial cancer treatment and diagnostic applications.

Generation of glycan-specific nanobodies.

The development of antibodies that target specific glycan structures on cancer cells or human pathogens poses a significant challenge due to the immense complexity of naturally occurring glycans. Automated glycan assembly enables the production of structurally homogeneous glycans in amounts that are difficult to derive from natural sources. Nanobodies (Nbs) are the smallest antigen-binding domains of heavy-chain-only antibodies (hcAbs) found in camelids. To date, the development of glycan-specific Nbs using synthetic glycans has not been reported. Here, we use defined synthetic glycans for alpaca immunization to elicit glycan-specific hcAbs, and describe the identification, isolation, and production of a Nb specific for the tumor-associated carbohydrate antigen Globo-H. The Nb binds the terminal fucose of Globo-H and recognizes synthetic Globo-H in solution and native Globo-H on breast cancer cells with high specificity. These results demonstrate the potential of our approach for generating glycan-targeting Nbs to be used in biomedical and biotechnological applications.

Differential Impact of Fluid Shear Stress and YAP/TAZ on BMP/TGF-ß Induced Osteogenic Target Genes.

Bone is a remarkable dynamic structure, which integrates mechanical and biochemical signaling inputs. Interstitial fluid in the intramedullary space transmits signals derived from compression-induced fluid shear stress (FSS) to stimulate osteoblasts for bone formation. Using a flow system and human osteoblasts, this study demonstrates how BMP/TGF-ß  signaling integrates stimuli derived from FSS and YAP/TAZ and confirms these findings by transcriptome analyses. Here, FSS positively affects the phosphorylation of both SMAD1/5 and SMAD2/3, the respective BMP- and TGFß-R-SMADs. Increase in phosphorylated SMAD1/5 levels affects distinct target genes, which are susceptible to low levels of phosphorylated SMADs (such as ID1-3) or dependent on high levels of phosphorylated SMAD1/5 (NOG, noggin). Thus, FSS lowers the threshold for genes dependent on high levels of phosphorylated SMAD1/5 when less BMP is available. While the impact of FSS on direct BMP target genes is independent of YAP/TAZ, FSS acts cooperatively with YAP/TAZ on TGF-ß  target genes, which are shared by both pathways (such as CTGF). As mechanical stimuli are key in bone regeneration, their crosstalk to biochemical signaling pathways such as BMP and TGF-ß and YAP/TAZ acts on different levels, which allows now to think about new and more specified intervention strategies for age-related bone loss.

Genomic Evidence for Sensorial Adaptations to a Nocturnal Predatory Lifestyle in Owls.

Owls (Strigiformes) evolved specific adaptations to their nocturnal predatory lifestyle, such as asymmetrical ears, a facial disk, and a feather structure allowing silent flight. Owls also share some traits with diurnal raptors and other nocturnal birds, such as cryptic plumage patterns, reversed sexual size dimorphism, and acute vision and hearing. The genetic basis of some of these adaptations to a nocturnal predatory lifestyle has been studied by candidate gene approaches but rarely with genome-wide scans. Here, we used a genome-wide comparative analysis to test for selection in the early history of the owls. We estimated the substitution rates in the coding regions of 20 bird genomes, including 11 owls of which five were newly sequenced. Then, we tested for functional overrepresentation across the genes that showed signals of selection. In the ancestral branch of the owls, we found traces of positive selection in the evolution of genes functionally related to visual perception, especially to phototransduction, and to chromosome packaging. Several genes that have been previously linked to acoustic perception, circadian rhythm, and feather structure also showed signals of an accelerated evolution in the origin of the owls. We discuss the functions of the genes under positive selection and their putative association with the adaptation to the nocturnal predatory lifestyle of the owls.

Network integration and modelling of dynamic drug responses at multi-omics levels.

Uncovering cellular responses from heterogeneous genomic data is crucial for molecular medicine in particular for drug safety. This can be realized by integrating the molecular activities in networks of interacting proteins. As proof-of-concept we challenge network modeling with time-resolved proteome, transcriptome and methylome measurements in iPSC-derived human 3D cardiac microtissues to elucidate adverse mechanisms of anthracycline cardiotoxicity measured with four different drugs (doxorubicin, epirubicin, idarubicin and daunorubicin). Dynamic molecular analysis at in vivo drug exposure levels reveal a network of 175 disease-associated proteins and identify common modules of anthracycline cardiotoxicity in vitro, related to mitochondrial and sarcomere function as well as remodeling of extracellular matrix. These in vitro-identified modules are transferable and are evaluated with biopsies of cardiomyopathy patients. This to our knowledge most comprehensive study on anthracycline cardiotoxicity demonstrates a reproducible workflow for molecular medicine and serves as a template for detecting adverse drug responses from complex omics data.

Esrrb Unlocks Silenced Enhancers for Reprogramming to Naive Pluripotency.

Transcription factor (TF)-mediated reprogramming to pluripotency is a slow and inefficient process, because most pluripotency TFs fail to access relevant target sites in a refractory chromatin environment. It is still unclear how TFs actually orchestrate the opening of repressive chromatin during the long latency period of reprogramming. Here, we show that the orphan nuclear receptor Esrrb plays a pioneering role in recruiting the core pluripotency factors Oct4, Sox2, and Nanog to inactive enhancers in closed chromatin during the reprogramming of epiblast stem cells. Esrrb binds to silenced enhancers containing stable nucleosomes and hypermethylated DNA, which are inaccessible to the core factors. Esrrb binding is accompanied by local loss of DNA methylation, LIF-dependent engagement of p300, and nucleosome displacement, leading to the recruitment of core factors within approximately 2 days. These results suggest that TFs can drive rapid remodeling of the local chromatin structure, highlighting the remarkable plasticity of stable epigenetic information.

Epigenomic profiling of non-small cell lung cancer xenografts uncover LRP12 DNA methylation as predictive biomarker for carboplatin resistance.

BACKGROUND: Non-small cell lung cancer (NSCLC) is the most common cause of cancer-related deaths worldwide and is primarily treated with radiation, surgery, and platinum-based drugs like cisplatin and carboplatin. The major challenge in the treatment of NSCLC patients is intrinsic or acquired resistance to chemotherapy. Molecular markers predicting the outcome of the patients are urgently needed. METHODS: Here, we employed patient-derived xenografts (PDXs) to detect predictive methylation biomarkers for platin-based therapies. We used MeDIP-Seq to generate genome-wide DNA methylation profiles of 22 PDXs, their parental primary NSCLC, and their corresponding normal tissues and complemented the data with gene expression analyses of the same tissues. Candidate biomarkers were validated with quantitative methylation-specific PCRs (qMSP) in an independent cohort. RESULTS: Comprehensive analyses revealed that differential methylation patterns are highly similar, enriched in PDXs and lung tumor-specific when comparing differences in methylation between PDXs versus primary NSCLC. We identified a set of 40 candidate regions with methylation correlated to carboplatin response and corresponding inverse gene expression pattern even before therapy. This analysis led to the identification of a promoter CpG island methylation of LDL receptor-related protein 12 (LRP12) associated with increased resistance to carboplatin. Validation in an independent patient cohort (n = 35) confirmed that LRP12 methylation status is predictive for therapeutic response of NSCLC patients to platin therapy with a sensitivity of 80% and a specificity of 84% (p < 0.01). Similarly, we find a shorter survival time for patients with LRP12 hypermethylation in the TCGA data set for NSCLC (lung adenocarcinoma). CONCLUSIONS: Using an epigenome-wide sequencing approach, we find differential methylation patterns from primary lung cancer and PDX-derived cancers to be very similar, albeit with a lower degree of differential methylation in primary tumors. We identify LRP12 DNA methylation as a powerful predictive marker for carboplatin resistance. These findings outline a platform for the identification of epigenetic therapy resistance biomarkers based on PDX NSCLC models.

Haplotype-resolved sweet potato genome traces back its hexaploidization history.

Here we present the 15 pseudochromosomes of sweet potato, Ipomoea batatas, the seventh most important crop in the world and the fourth most significant in China. By using a novel haplotyping method based on genome assembly, we have produced a half haplotype-resolved genome from ~296 Gb of paired-end sequence reads amounting to roughly 67-fold coverage. By phylogenetic tree analysis of homologous chromosomes, it was possible to estimate the time of two recent whole-genome duplication events as occurring about 0.8 and 0.5 million years ago. This half haplotype-resolved hexaploid genome represents the first successful attempt to investigate the complexity of chromosome sequence composition directly in a polyploid genome, using sequencing of the polyploid organism itself rather than any of its simplified proxy relatives. Adaptation and application of our approach should provide higher resolution in future genomic structure investigations, especially for similarly complex genomes.

A cancer specific hypermethylation signature of the TERT promoter predicts biochemical relapse in prostate cancer: a retrospective cohort study.

The identification of new biomarkers to differentiate between indolent and aggressive prostate tumors is an important unmet need. We examined the role of THOR (TERT Hypermethylated Oncological Region) as a diagnostic and prognostic biomarker in prostate cancer (PCa).We analyzed THOR in common cancers using genome-wide methylation arrays. Methylation status of the whole TERT gene in benign and malignant prostate samples was determined by MeDIP-Seq. The prognostic role of THOR in PCa was assessed by pyrosequencing on discovery and validation cohorts from patients who underwent radical prostatectomy with long-term follow-up data.Most cancers (n = 3056) including PCa (n = 300) exhibited hypermethylation of THOR. THOR was the only region within the TERT gene that is differentially methylated between normal and malignant prostate tissue (p < 0.0001). Also, THOR was significantly hypermethylated in PCa when compared to paired benign tissues (n = 164, p < 0.0001). THOR hypermethylation correlated with Gleason scores and was associated with tumor invasiveness (p = 0.0147). Five years biochemical progression free survival (BPFS) for PCa patients in the discovery cohort was 87% (95% CI 73-100) and 65% (95% CI 52-78) for THOR non-hypermethylated and hypermethylated cancers respectively (p = 0.01). Similar differences in BPFS were noted in the validation cohort (p = 0.03). Importantly, THOR was able to predict outcome in the challenging (Gleason 6 and 7 (3 + 4)) PCa (p = 0.007). For this group, THOR was an independent risk factor for BPFS with a hazard-ratio of 3.685 (p = 0.0247). Finally, THOR hypermethylation more than doubled the risk of recurrence across all PSA levels (OR 2.5, p = 0.02).

Next-generation sequencing technologies for DNA methylation analyses in cancer genomics.

For the first time, the development of next-generation sequencing technologies has brought about tools to investigate epigenetic alterations in an unbiased, yet genome-wide approach. The importance of this innovative technology is undeniable since it has already been established that changes in DNA methylation play an important role in cancer initiation and progression. The first methylation maps have already been created, and it is only a matter of time until the complete epigenetic maps of healthy and diseased human genomes are available. In this review, we summarize the use of next-generation sequencing for diverse epigenetic technologies, give an overview of the status quo and outline future perspectives for its application in oncology and basic research.

Somatic mutation profiles of MSI and MSS colorectal cancer identified by whole exome next generation sequencing and bioinformatics analysis.

BACKGROUND: Colorectal cancer (CRC) is with approximately 1 million cases the third most common cancer worldwide. Extensive research is ongoing to decipher the underlying genetic patterns with the hope to improve early cancer diagnosis and treatment. In this direction, the recent progress in next generation sequencing technologies has revolutionized the field of cancer genomics. However, one caveat of these studies remains the large amount of genetic variations identified and their interpretation. METHODOLOGY/PRINCIPAL FINDINGS: Here we present the first work on whole exome NGS of primary colon cancers. We performed 454 whole exome pyrosequencing of tumor as well as adjacent not affected normal colonic tissue from microsatellite stable (MSS) and microsatellite instable (MSI) colon cancer patients and identified more than 50,000 small nucleotide variations for each tissue. According to predictions based on MSS and MSI pathomechanisms we identified eight times more somatic non-synonymous variations in MSI cancers than in MSS and we were able to reproduce the result in four additional CRCs. Our bioinformatics filtering approach narrowed down the rate of most significant mutations to 359 for MSI and 45 for MSS CRCs with predicted altered protein functions. In both CRCs, MSI and MSS, we found somatic mutations in the intracellular kinase domain of bone morphogenetic protein receptor 1A, BMPR1A, a gene where so far germline mutations are associated with juvenile polyposis syndrome, and show that the mutations functionally impair the protein function. CONCLUSIONS/SIGNIFICANCE: We conclude that with deep sequencing of tumor exomes one may be able to predict the microsatellite status of CRC and in addition identify potentially clinically relevant mutations.