Developmental dynamics of lncRNAs across mammalian organs and species.

Although many long noncoding RNAs (lncRNAs) have been identified in human and other mammalian genomes, there has been limited systematic functional characterization of these elements. In particular, the contribution of lncRNAs to organ development remains largely unexplored. Here we analyse the expression patterns of lncRNAs across developmental time points in seven major organs, from early organogenesis to adulthood, in seven species (human, rhesus macaque, mouse, rat, rabbit, opossum and chicken). Our analyses identified approximately 15,000 to 35,000 candidate lncRNAs in each species, most of which show species specificity. We characterized the expression patterns of lncRNAs across developmental stages, and found many with dynamic expression patterns across time that show signatures of enrichment for functionality. During development, there is a transition from broadly expressed and conserved lncRNAs towards an increasing number of lineage- and organ-specific lncRNAs. Our study provides a resource of candidate lncRNAs and their patterns of expression and evolutionary conservation across mammalian organ development.

Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage.

Sex chromosomes differentiated from different ancestral autosomes in various vertebrate lineages. Here, we trace the functional evolution of the XY Chromosomes of the green anole lizard (Anolis carolinensis), on the basis of extensive high-throughput genome, transcriptome and histone modification sequencing data and revisit dosage compensation evolution in representative mammals and birds with substantial new expression data. Our analyses show that Anolis sex chromosomes represent an ancient XY system that originated at least ≈160 million years ago in the ancestor of Iguania lizards, shortly after the separation from the snake lineage. The age of this system approximately coincides with the ages of the avian and two mammalian sex chromosomes systems. To compensate for the almost complete Y Chromosome degeneration, X-linked genes have become twofold up-regulated, restoring ancestral expression levels. The highly efficient dosage compensation mechanism of Anolis represents the only vertebrate case identified so far to fully support Ohno's original dosage compensation hypothesis. Further analyses reveal that X up-regulation occurs only in males and is mediated by a male-specific chromatin machinery that leads to global hyperacetylation of histone H4 at lysine 16 specifically on the X Chromosome. The green anole dosage compensation mechanism is highly reminiscent of that of the fruit fly, Drosophila melanogaster Altogether, our work unveils the convergent emergence of a Drosophila-like dosage compensation mechanism in an ancient reptilian sex chromosome system and highlights that the evolutionary pressures imposed by sex chromosome dosage reductions in different amniotes were resolved in fundamentally different ways.

Origins and functional evolution of Y chromosomes across mammals.

Y chromosomes underlie sex determination in mammals, but their repeat-rich nature has hampered sequencing and associated evolutionary studies. Here we trace Y evolution across 15 representative mammals on the basis of high-throughput genome and transcriptome sequencing. We uncover three independent sex chromosome originations in mammals and birds (the outgroup). The original placental and marsupial (therian) Y, containing the sex-determining gene SRY, emerged in the therian ancestor approximately 180 million years ago, in parallel with the first of five monotreme Y chromosomes, carrying the probable sex-determining gene AMH. The avian W chromosome arose approximately 140 million years ago in the bird ancestor. The small Y/W gene repertoires, enriched in regulatory functions, were rapidly defined following stratification (recombination arrest) and erosion events and have remained considerably stable. Despite expression decreases in therians, Y/W genes show notable conservation of proto-sex chromosome expression patterns, although various Y genes evolved testis-specificities through differential regulatory decay. Thus, although some genes evolved novel functions through spatial/temporal expression shifts, most Y genes probably endured, at least initially, because of dosage constraints.

PACCMIT/PACCMIT-CDS: identifying microRNA targets in 3' UTRs and coding sequences.

The purpose of the proposed web server, publicly available at http://paccmit.epfl.ch, is to provide a user-friendly interface to two algorithms for predicting messenger RNA (mRNA) molecules regulated by microRNAs: (i) PACCMIT (Prediction of ACcessible and/or Conserved MIcroRNA Targets), which identifies primarily mRNA transcripts targeted in their 3' untranslated regions (3' UTRs), and (ii) PACCMIT-CDS, designed to find mRNAs targeted within their coding sequences (CDSs). While PACCMIT belongs among the accurate algorithms for predicting conserved microRNA targets in the 3' UTRs, the main contribution of the web server is 2-fold: PACCMIT provides an accurate tool for predicting targets also of weakly conserved or non-conserved microRNAs, whereas PACCMIT-CDS addresses the lack of similar portals adapted specifically for targets in CDS. The web server asks the user for microRNAs and mRNAs to be analyzed, accesses the precomputed P-values for all microRNA-mRNA pairs from a database for all mRNAs and microRNAs in a given species, ranks the predicted microRNA-mRNA pairs, evaluates their significance according to the false discovery rate and finally displays the predictions in a tabular form. The results are also available for download in several standard formats.

Searching the coding region for microRNA targets.

Finding microRNA targets in the coding region is difficult due to the overwhelming signal encoding the amino acid sequence. Here, we introduce an algorithm (called PACCMIT-CDS) that finds potential microRNA targets within coding sequences by searching for conserved motifs that are complementary to the microRNA seed region and also overrepresented in comparison with a background model preserving both codon usage and amino acid sequence. Precision and sensitivity of PACCMIT-CDS are evaluated using PAR-CLIP and proteomics data sets. Thanks to the properly constructed background, the new algorithm achieves a lower rate of false positives and better ranking of predictions than do currently available algorithms, which were designed to find microRNA targets within 3' UTRs.

Analysis of the accessibility of CLIP bound sites reveals that nucleation of the miRNA:mRNA pairing occurs preferentially at the 3'-end of the seed match.

To find out whether the AGO-miRNA complex is more sensitive to the accessibility of a particular region inside the seed match, we analyze in detail the accessibility of a wide set of miRNA binding sites validated by PAR-CLIP and HITS-CLIP experiments. Our analysis reveals that nucleotides at the 3'-end of bound seed matches are significantly more accessible than nucleotides at the 5'-end as well as nucleotides at any positions in the unbound seed matches. We show that the accessibility of a single nucleotide at the 3'-end is more effective than the accessibility of several nucleotides at the 5'-end in discriminating between functional and nonfunctional binding sites. Analysis of mRNA and protein fold changes induced by miRNA overexpression demonstrates that genes with accessible nucleation regions at the 3'-end are down-regulated more strongly than genes whose accessible nucleation regions are located elsewhere within the seed match. We also observed an increase in the precision of the miRNA target prediction algorithm PACMIT when accessibility toward the 3'-end of the seed match was required. The pronounced sensitivity of the AGO-miRNA complex to the accessibility of the 3'-end of the seed match suggests that, in most cases, nucleation occurs in this region. We show that this conclusion is consistent with previous experimental studies.

Efficient use of accessibility in microRNA target prediction.

Considering accessibility of the 3'UTR is believed to increase the precision of microRNA target predictions. We show that, contrary to common belief, ranking by the hybridization energy or by the sum of the opening and hybridization energies, used in currently available algorithms, is not an efficient way to rank predictions. Instead, we describe an algorithm which also considers only the accessible binding sites but which ranks predictions according to over-representation. When compared with experimentally validated and refuted targets in the fruit fly and human, our algorithm shows a remarkable improvement in precision while significantly reducing the computational cost in comparison with other free energy based methods. In the human genome, our algorithm has at least twice higher precision than other methods with their default parameters. In the fruit fly, we find five times more validated targets among the top 500 predictions than other methods with their default parameters. Furthermore, using a common statistical framework we demonstrate explicitly the advantages of using the canonical ensemble instead of using the minimum free energy structure alone. We also find that 'naïve' global folding sometimes outperforms the local folding approach.

GIInger predicts homologous recombination deficiency and patient response to PARPi treatment from shallow genomic profiles.

Homologous recombination deficiency (HRD) is a predictive biomarker for poly(ADP-ribose) polymerase 1 inhibitor (PARPi) sensitivity. Routine HRD testing relies on identifying BRCA mutations, but additional HRD-positive patients can be identified by measuring genomic instability (GI), a consequence of HRD. However, the cost and complexity of available solutions hamper GI testing. We introduce a deep learning framework, GIInger, that identifies GI from HRD-induced scarring observed in low-pass whole-genome sequencing data. GIInger seamlessly integrates into standard BRCA testing workflows and yields reproducible results concordant with a reference method in a multisite study of 327 ovarian cancer samples. Applied to a BRCA wild-type enriched subgroup of 195 PAOLA-1 clinical trial patients, GIInger identified HRD-positive patients who experienced significantly extended progression-free survival when treated with PARPi. GIInger is, therefore, a cost-effective and easy-to-implement method for accurately stratifying patients with ovarian cancer for first-line PARPi treatment.

Genome sequence of the Drosophila melanogaster male-killing Spiroplasma strain MSRO endosymbiont.

UNLABELLED: Spiroplasmas are helical and motile members of a cell wall-less eubacterial group called Mollicutes. Although all spiroplasmas are associated with arthropods, they exhibit great diversity with respect to both their modes of transmission and their effects on their hosts; ranging from horizontally transmitted pathogens and commensals to endosymbionts that are transmitted transovarially (i.e., from mother to offspring). Here we provide the first genome sequence, along with proteomic validation, of an endosymbiotic inherited Spiroplasma bacterium, the Spiroplasma poulsonii MSRO strain harbored by Drosophila melanogaster. Comparison of the genome content of S. poulsonii with that of horizontally transmitted spiroplasmas indicates that S. poulsonii has lost many metabolic pathways and transporters, demonstrating a high level of interdependence with its insect host. Consistent with genome analysis, experimental studies showed that S. poulsonii metabolizes glucose but not trehalose. Notably, trehalose is more abundant than glucose in Drosophila hemolymph, and the inability to metabolize trehalose may prevent S. poulsonii from overproliferating. Our study identifies putative virulence genes, notably, those for a chitinase, the H2O2-producing glycerol-3-phosphate oxidase, and enzymes involved in the synthesis of the eukaryote-toxic lipid cardiolipin. S. poulsonii also expresses on the cell membrane one functional adhesion-related protein and two divergent spiralin proteins that have been implicated in insect cell invasion in other spiroplasmas. These lipoproteins may be involved in the colonization of the Drosophila germ line, ensuring S. poulsonii vertical transmission. The S. poulsonii genome is a valuable resource to explore the mechanisms of male killing and symbiont-mediated protection, two cardinal features of many facultative endosymbionts. IMPORTANCE: Most insect species, including important disease vectors and crop pests, harbor vertically transmitted endosymbiotic bacteria. These endosymbionts play key roles in their hosts' fitness, including protecting them against natural enemies and manipulating their reproduction in ways that increase the frequency of symbiont infection. Little is known about the molecular mechanisms that underlie these processes. Here, we provide the first genome draft of a vertically transmitted male-killing Spiroplasma bacterium, the S. poulsonii MSRO strain harbored by D. melanogaster. Analysis of the S. poulsonii genome was complemented by proteomics and ex vivo metabolic experiments. Our results indicate that S. poulsonii has reduced metabolic capabilities and expresses divergent membrane lipoproteins and potential virulence factors that likely participate in Spiroplasma-host interactions. This work fills a gap in our knowledge of insect endosymbionts and provides tools with which to decipher the interaction between Spiroplasma bacteria and their well-characterized host D. melanogaster, which is emerging as a model of endosymbiosis.

Optimal use of conservation and accessibility filters in microRNA target prediction.

It is generally accepted that filtering microRNA (miRNA) target predictions by conservation or by accessibility can reduce the false discovery rate. However, these two strategies are usually not exploited in a combined and flexible manner. Here, we introduce PACCMIT, a flexible method that filters miRNA binding sites by their conservation, accessibility, or both. The improvement in performance obtained with each of these three filters is demonstrated on the prediction of targets for both i) highly and ii) weakly conserved miRNAs, i.e., in two scenarios in which the miRNA-target interactions are subjected to different evolutionary pressures. We show that in the first scenario conservation is a better filter than accessibility (as both sensitivity and precision are higher among the top predictions) and that the combined filter improves performance of PACCMIT even further. In the second scenario, on the other hand, the accessibility filter performs better than both the conservation and combined filters, suggesting that the site conservation is not equally effective in rejecting false positive predictions for all miRNAs. Regarding the quality of the ranking criterion proposed by Robins and Press and used in PACCMIT, it is shown that top ranking interactions correspond to more downregulated proteins than do the lower ranking interactions. Comparison with several other target prediction algorithms shows that the ranking of predictions provided by PACCMIT is at least as good as the ranking generated by other conservation-based methods and considerably better than the energy-based ranking used in other accessibility-based methods.

Characterization of activity landscapes using 2D and 3D similarity methods: consensus activity cliffs.

Activity landscape characterization has been demonstrated to be a valuable tool in lead optimization, virtual screening, and computational modeling of active compounds. In this work, we present a general protocol to explore systematically the activity landscape of a lead series using 11 2D and 3D structural representations. As a test case we employed a set of 48 bicyclic guanidines (BCGs) with kappa-opioid receptor binding affinity, identified in our group. MACCS keys, graph-based three point pharmacophores, circular fingerprints, ROCS shape descriptors, and the TARIS approach, that compares structures based on molecular electrostatic potentials, were employed as orthogonal descriptors. Based on 'activity cliffs' common to a series of descriptors, we introduce the concept of consensus activity cliffs. Results for the current test case suggest that the presence or absence of a methoxybenzyl group may lead to different modes of binding for the active BCGs with the kappa-opioid receptor. The most active compound (IC50 = 37 nM) is involved in a number of consensus cliffs making it a more challenge query for future virtual screening than would be expected from affinity alone. Results also reveal the importance of screening high density combinatorial libraries, especially in the "cliff-rich" regions of activity landscapes. The protocol presented here can be applied to other data sets to develop a consensus model of the activity landscape.

Graph theoretical similarity approach to compare molecular electrostatic potentials.

In this work we introduce a graph theoretical method to compare MEPs, which is independent of molecular alignment. It is based on the edit distance of weighted rooted trees, which encode the geometrical and topological information of Negative Molecular Isopotential Surfaces. A meaningful chemical classification of a set of 46 molecules with different functional groups was achieved. Structure--activity relationships for the corticosteroid binding affinity (CBG) of 31 steroids by means of hierarchical clustering resulted in a clear partitioning in high, intermediate, and low activity groups, whereas the results from quantitative structure--activity relationships, obtained from a partial least-squares analysis, showed comparable or better cross-validated correlation coefficients than the ones reported for previous methods based solely in the MEP.

Radial scan of the molecular electrostatic potential of RNA double helices: an application to the enzyme-tRNA recognition.

We introduced a method to characterize quantitatively the molecular electrostatic potential (MEP) of the minor and major grooves of base pairs located at nucleic acid double helices. By means of a radial MEP scan, we obtained a n-tuple of potential values corresponding to each groove, which can be analyzed by plotting the MEP values as a function of the angle in the radial scan. We studied base pairs of two different tRNAs, relevant in the recognition process with their cognate aminoacyl tRNA synthetases (aaRSs), in order to correlate their electrostatic behavior with the corresponding aminoacylation activity. We analyzed the first three base pairs of the Escherichia coli tRNA(Ala) acceptor stem, finding several cases where the MEP profiles obtained from the plots are in agreement with the reported aminoacylation activities. Additionally, a non-hierarchical clustering performed over the MEP n-tuples resulted in meaningful classifications that correlate with the activity and with the predicted stereochemistry of the reaction. We also studied the first two base pairs of the E. coli tRNA(Thr) acceptor stem but constraining the analysis to the angle intervals that seem relevant for the binding sites of the enzyme. These intervals were deduced from the ThrRS-tRNA(Thr) complex crystal structure. In this case, we also found a good agreement between the MEP profiles and the activity, supporting the idea that the tRNA identity elements function is to allow an optimal electrostatic complementarity between the aminoacyl-tRNA synthetase and the tRNA.

CARACTERIZACIÓN DEL TALLO ACEPTOR DEL tARN MEDIANTE DESCRIPTORES LOCALES BASADOS EN CARGAS PARCIALES / tRNA ACCEPTOR STEM CHARACTERIZATION BY MEANS OF LOCAL DESCRIPTORS BASED ON PARTIAL CHARGES / CARACTERIZAÇÃO DO TALO ACEITADOR DO tRNA MEDIANTE DESCRITORES LOCAIS BASEADOS EM CARGAS PARCIAIS

En este trabajo se caracteriza la distribución de carga del tallo aceptor del tARN, considerando todas las posibles combinaciones de pares Watson-Crick. El estudio se realizó con 256 fragmentos moleculares de 10 nucleótidos que modelan los tres primeros pares del tallo aceptor, la base diferenciadora y el extremo CCA. Para caracterizar los nucleótidos se proponen dos descriptores locales basados en la distribución de carga de la base nitrogenada de cada nucleótido, los cuales se calculan a partir de las cargas parciales de Mulliken obtenidas de cálculos HF/6-31G. La caracterización y clasificación de los tallos según estos descriptores mostró cómo la base diferenciadora tiene un comportamiento particular respecto a los demás nucleótidos del tallo y una fuerte influencia sobre el extremo CCA. La clasificación de nueve variaciones del tallo aceptor del tARNAla mostró una buena relación estructura-actividad que pone en evidencia la bondad de los descriptores propuestos para caracterizar de manera local la distribución de carga de estas biomoléculas.

In this work the charge distribution of the tRNA acceptor stem is characterized, considering all the possible Watson- Crick base pair combinations. 256 RNA fragments modeled by 10 nucleotides were used in order to model the first three pairs of the acceptor stem, the discriminator base and the CCA end. We propose two local charge descriptors based on the charge distribution of the nitrogenated base to characterize each nucleotide. These descriptors were computed from atomic partial charges derived from HF/6-31G calculations. From the characterization and classification of the stems according to the proposed descriptors, we found a special behavior for the discriminator base (in contrast to the other positions) and a strong effect of this position on the CCA end. The classification of nine variations of the tRNAAla acceptor stem showed a good structure-activity relationship that makes evident the usefulness of the proposed descriptors to characterize the local charge distributions of these biomolecules.

Nesse estudo é caracterizada a distribuição da carga do talo aceitador considerando- se todas as combinações possíveis dos pares Watson-Crick. O estudo realizouse com 256 fragmentos moleculares dos 10 nucleotídeos que modelam os três primeiros pares do talo aceitador, a base diferenciadora e o extremo CCA. Com o intuito de caracterizar cada nucleotídeo, foram propostos dois descritores locais baseados na distribuição de carga da base nitrogenada de cada nucleotídeo, os quais se calculam a partir das cargas parciais de Mulliken obtidas de cálculos HF/6-31G. A caracterização e classificação dos talos segundo esses descritores demonstrou um particular comportamento da base diferenciadora em relação aos demais nucleotídeos do talo e uma forte influência sobre o extremo CCA. A classificação de nove variações do talo aceitador do tRNA mostrou uma boa relação estrutura-atividade que colocam em evidência a utilidade dos descritores propostos para caracterizar de maneira local a distribuição de carga dessas bio-moléculas.