Microenvironmental reorganization in brain tumors following radiotherapy and recurrence revealed by hyperplexed immunofluorescence imaging

The tumor microenvironment plays a crucial role in determining response to treatment. This involves a series of interconnected changes in the cellular landscape, spatial organization, and extracellular matrix composition. However, assessing these alterations simultaneously is challenging from a spatial perspective, due to the limitations of current high-dimensional imaging techniques and the extent of intratumoral heterogeneity over large lesion areas. In this study, we introduce a spatial proteomic workflow termed Hyperplexed Immunofluorescence Imaging (HIFI) that overcomes these limitations. HIFI allows for the simultaneous analysis of > 45 markers in fragile tissue sections at high magnification, using a cost-effective high-throughput workflow. We integrate HIFI with machine learning feature detection, graph-based network analysis, and cluster-based neighborhood analysis to analyze the microenvironment response to radiation therapy in a preclinical model of glioblastoma, and compare this response to a mouse model of breast-to-brain metastasis. Here we show that glioblastomas undergo extensive spatial reorganization of immune cell populations and structural architecture in response to treatment, while brain metastases show no comparable reorganization. Our integrated spatial analyses reveal highly divergent responses to radiation therapy between brain tumor models, despite equivalent radiotherapy benefit.

NIPMAP: niche-phenotype mapping of multiplex histology data by community ecology

Advances in multiplex histology allow surveying millions of cells, dozens of cell types, and up to thousands of phenotypes within the spatial context of tissue sections. This leads to a combinatorial challenge in (a) summarizing the cellular and phenotypic architecture of tissues and (b) identifying phenotypes with interesting spatial architecture. To address this, we combine ideas from community ecology and machine learning into niche-phenotype mapping (NIPMAP). NIPMAP takes advantage of geometric constraints on local cellular composition imposed by the niche structure of tissues in order to automatically segment tissue sections into niches and their interfaces. Projecting phenotypes on niches and their interfaces identifies previously-reported and previously-unreported spatially-driven phenotypes, concisely summarizes the phenotypic architecture of tissues, and reveals fundamental properties of tissue architecture. NIPMAP is applicable to both protein and RNA multiplex histology of healthy and diseased tissue. An open-source R/Python package implements NIPMAP.

Microscopy-based phenotypic monitoring of MDA-MB-231 spheroids allows the evaluation of phenotype-directed therapy

Breast cancer (BC) is the most commonly diagnosed cancer among women. Prognosis has improved over the years, to a large extent, owing to personalized therapy informed by molecular profiling of hormone receptors. However, there is a need for new therapeutic approaches for a subgroup of BCs lacking molecular markers, the Triple Negative Breast Cancer (TNBC) subgroup. TNBC is the most aggressive type of BC, lacks an effective standard of care, shows high levels of resistance and relapse is often inevitable. High resistance to therapy has been hypothesized to be associated with high intratumoral phenotypic heterogeneity. To characterize and treat this phenotypic heterogeneity, we optimized a whole-mount staining and image analysis protocol for three-dimensions (3D) spheroids. Applying this protocol to TNBC spheroids located in the outer region of the spheroid the cells with selected phenotypes: dividing, migrating, and high mitochondrial mass phenotypes. To evaluate the relevance of phenotype-based targeting these cell populations were targeted with Paclitaxel, Trametinib, and Everolimus, respectively, in a dose-dependent manner. Single agents cannot specifically target all phenotypes at the same time. Therefore, we combined drugs that should target independent phenotype. With this rationale we observed that combining Trametinib and Everolimus achieves the highest cytotoxicity at lower doses from all the tested combinations. These findings suggest a rational approach to design treatments can be evaluated in spheroids prior to pre-clinical models and potentially reduce adverse effects.

Controls for Phylogeny and Robust Analysis in Pareto Task Inference

Understanding the tradeoffs faced by organisms is a major goal of evolutionary biology. One of the main approaches for identifying these tradeoffs is Pareto task inference (ParTI). Two recent papers claim that results obtained in ParTI studies are spurious due to phylogenetic dependence (Mikami T, Iwasaki W. 2021. The flipping t-ratio test: phylogenetically informed assessment of the Pareto theory for phenotypic evolution. Methods Ecol Evol. 12(4):696-706) or hypothetical p-hacking and population-structure concerns (Sun M, Zhang J. 2021. Rampant false detection of adaptive phenotypic optimization by ParTI-based Pareto front inference. Mol Biol Evol. 38(4):1653-1664). Here, we show that these claims are baseless. We present a new method to control for phylogenetic dependence, called SibSwap, and show that published ParTI inference is robust to phylogenetic dependence. We show how researchers avoided p-hacking by testing for the robustness of preprocessing choices. We also provide new methods to control for population structure and detail the experimental tests of ParTI in systems ranging from ammonites to cancer gene expression. The methods presented here may help to improve future ParTI studies.

Tumour heterogeneity and the evolutionary trade-offs of cancer

Tumours vary in gene expression programmes and genetic alterations. Understanding this diversity and its biological meaning requires a theoretical framework, which could in turn guide the development of more accurate prognosis and therapy. Here, we review the theory of multi-task evolution of cancer, which is based upon the premise that tumours evolve in the host and face selection trade-offs between multiple biological functions. This theory can help identify the major biological tasks that cancer cells perform and the trade-offs between these tasks. It introduces the concept of specialist tumours, which focus on one task, and generalist tumours, which perform several tasks. Specialist tumours are suggested to be sensitive to therapy targeting their main task. Driver mutations tune gene expression towards specific tasks in a tissue-dependent manner and thus help to determine whether a tumour is specialist or generalist. We discuss potential applications of the theory of multi-task evolution to interpret the spatial organization of tumours and intratumour heterogeneity.

Tumor diversity and the trade-off between universal cancer tasks

Recent advances have enabled powerful methods to sort tumors into prognosis and treatment groups. We are still missing, however, a general theoretical framework to understand the vast diversity of tumor gene expression and mutations. Here we present a framework based on multi-task evolution theory, using the fact that tumors need to perform multiple tasks that contribute to their fitness. We find that trade-offs between tasks constrain tumor gene-expression to a continuum bounded by a polyhedron whose vertices are gene-expression profiles, each specializing in one task. We find five universal cancer tasks across tissue-types: cell-division, biomass and energy, lipogenesis, immune-interaction and invasion and tissue-remodeling. Tumors that specialize in a task are sensitive to drugs that interfere with this task. Driver, but not passenger, mutations tune gene-expression towards specialization in specific tasks. This approach can integrate additional types of molecular data into a framework of tumor diversity grounded in evolutionary theory.

An active β-lactamase is a part of an orchestrated cell wall stress resistance network of Bacillus subtilis and related rhizosphere species

A hallmark of the Gram-positive bacteria, such as the soil-dwelling bacterium Bacillus subtilis, is their cell wall. Here, we report that d-leucine and flavomycin, biofilm inhibitors targeting the cell wall, activate the β-lactamase PenP. This β-lactamase contributes to ampicillin resistance in B. subtilis under all conditions tested. In contrast, both Spo0A, a master regulator of nutritional stress, and the general cell wall stress response, differentially contribute to β-lactam resistance under different conditions. To test whether β-lactam resistance and β-lactamase genes are widespread in other Bacilli, we isolated Bacillus species from undisturbed soils, and found that their genomes can encode up to five β-lactamases with differentiated activity spectra. Surprisingly, the activity of environmental β-lactamases and PenP, as well as the general stress response, resulted in a similarly reduced lag phase of the culture in the presence of β-lactam antibiotics, with little or no impact on the logarithmic growth rate. The length of the lag phase may determine the outcome of the competition between β-lactams and β-lactamases producers. Overall, our work suggests that antibiotic resistance genes in B. subtilis and related species are ancient and widespread, and could be selected by interspecies competition in undisturbed soils.

Central dogma rates and the trade-off between precision and economy in gene expression

Steady-state protein abundance is set by four rates: transcription, translation, mRNA decay and protein decay. A given protein abundance can be obtained from infinitely many combinations of these rates. This raises the question of whether the natural rates for each gene result from historical accidents, or are there rules that give certain combinations a selective advantage? We address this question using high-throughput measurements in rapidly growing cells from diverse organisms to find that about half of the rate combinations do not exist: genes that combine high transcription with low translation are strongly depleted. This depletion is due to a trade-off between precision and economy: high transcription decreases stochastic fluctuations but increases transcription costs. Our theory quantitatively explains which rate combinations are missing, and predicts the curvature of the fitness function for each gene. It may guide the design of gene circuits with desired expression levels and noise.

Massively Parallel Interrogation of the Effects of Gene Expression Levels on Fitness

Data of gene expression levels across individuals, cell types, and disease states is expanding, yet our understanding of how expression levels impact phenotype is limited. Here, we present a massively parallel system for assaying the effect of gene expression levels on fitness in Saccharomyces cerevisiae by systematically altering the expression level of ∼100 genes at ∼100 distinct levels spanning a 500-fold range at high resolution. We show that the relationship between expression levels and growth is gene and environment specific and provides information on the function, stoichiometry, and interactions of genes. Wild-type expression levels in some conditions are not optimal for growth, and genes whose fitness is greatly affected by small changes in expression level tend to exhibit lower cell-to-cell variability in expression. Our study addresses a fundamental gap in understanding the functional significance of gene expression regulation and offers a framework for evaluating the phenotypic effects of expression variation.

Geometry of the Gene Expression Space of Individual Cells

There is a revolution in the ability to analyze gene expression of single cells in a tissue. To understand this data we must comprehend how cells are distributed in a high-dimensional gene expression space. One open question is whether cell types form discrete clusters or whether gene expression forms a continuum of states. If such a continuum exists, what is its geometry? Recent theory on evolutionary trade-offs suggests that cells that need to perform multiple tasks are arranged in a polygon or polyhedron (line, triangle, tetrahedron and so on, generally called polytopes) in gene expression space, whose vertices are the expression profiles optimal for each task. Here, we analyze single-cell data from human and mouse tissues profiled using a variety of single-cell technologies. We fit the data to shapes with different numbers of vertices, compute their statistical significance, and infer their tasks. We find cases in which single cells fill out a continuum of expression states within a polyhedron. This occurs in intestinal progenitor cells, which fill out a tetrahedron in gene expression space. The four vertices of this tetrahedron are each enriched with genes for a specific task related to stemness and early differentiation. A polyhedral continuum of states is also found in spleen dendritic cells, known to perform multiple immune tasks: cells fill out a tetrahedron whose vertices correspond to key tasks related to maturation, pathogen sensing and communication with lymphocytes. A mixture of continuum-like distributions and discrete clusters is found in other cell types, including bone marrow and differentiated intestinal crypt cells. This approach can be used to understand the geometry and biological tasks of a wide range of single-cell datasets. The present results suggest that the concept of cell type may be expanded. In addition to discreet clusters in gene-expression space, we suggest a new possibility: a continuum of states within a polyhedron, in which the vertices represent specialists at key tasks.

In the past, biological experiments usually pooled together millions of cells, masking the differences between individual cells. Current technology takes a big step forward by measuring gene expression from individual cells. Interpreting this data is challenging because we need to understand how cells are arranged in a high dimensional gene expression space. Here we test recent theory that suggests that cells facing multiple tasks should be arranged in simple low dimensional polygons or polyhedra (generally called polytopes). The vertices of the polytopes are gene expression profiles optimal for each of the tasks. We find evidence for such simplicity in a variety of tissues—spleen, bone marrow, intestine—analyzed by different single-cell technologies. We find that cells are distributed inside polytopes, such as tetrahedrons or four-dimensional simplexes, with cells closest to each vertex responsible for a different key task. For example, intestinal progenitor cells that give rise to the other cell types show a continuous distribution in a tetrahedron whose vertices correspond to several key sub-tasks. Immune dendritic cells likewise are continuously distributed between key immune tasks. This approach of testing whether data falls in polytopes may be useful for interpreting a variety of single-cell datasets in terms of biological tasks.

miR-184 Regulates Pancreatic β-Cell Function According to Glucose Metabolism

In response to fasting or hyperglycemia, the pancreatic β-cell alters its output of secreted insulin; however, the pathways governing this adaptive response are not entirely established. Although the precise role of microRNAs (miRNAs) is also unclear, a recurring theme emphasizes their function in cellular stress responses. We recently showed that miR-184, an abundant miRNA in the β-cell, regulates compensatory proliferation and secretion during insulin resistance. Consistent with previous studies showing miR-184 suppresses insulin release, expression of this miRNA was increased in islets after fasting, demonstrating an active role in the β-cell as glucose levels lower and the insulin demand ceases. Additionally, miR-184 was negatively regulated upon the administration of a sucrose-rich diet in Drosophila, demonstrating strong conservation of this pathway through evolution. Furthermore, miR-184 and its target Argonaute2 remained inversely correlated as concentrations of extracellular glucose increased, underlining a functional relationship between this miRNA and its targets. Lastly, restoration of Argonaute2 in the presence of miR-184 rescued suppression of miR-375-targeted genes, suggesting these genes act in a coordinated manner during changes in the metabolic context. Together, these results highlight the adaptive role of miR-184 according to glucose metabolism and suggest the regulatory role of this miRNA in energy homeostasis is highly conserved.

Identification and consequences of miRNA-target interactions--beyond repression of gene expression

Comparative genomics analyses and high-throughput experimental studies indicate that a microRNA (miRNA) binds to hundreds of sites across the transcriptome. Although the knockout of components of the miRNA biogenesis pathway has profound phenotypic consequences, most predicted miRNA targets undergo small changes at the mRNA and protein levels when the expression of the miRNA is perturbed. Alternatively, miRNAs can establish thresholds in and increase the coherence of the expression of their target genes, as well as reduce the cell-to-cell variability in target gene expression. Here, we review the recent progress in identifying miRNA targets and the emerging paradigms of how miRNAs shape the dynamics of target gene expression.

Linear superposition and prediction of bacterial promoter activity dynamics in complex conditions

Bacteria often face complex environments. We asked how gene expression in complex conditions relates to expression in simpler conditions. To address this, we obtained accurate promoter activity dynamical measurements on 94 genes in E. coli in environments made up of all possible combinations of four nutrients and stresses. We find that the dynamics across conditions is well described by two principal component curves specific to each promoter. As a result, the promoter activity dynamics in a combination of conditions is a weighted average of the dynamics in each condition alone. The weights tend to sum up to approximately one. This weighted-average property, called linear superposition, allows predicting the promoter activity dynamics in a combination of conditions based on measurements of pairs of conditions. If these findings apply more generally, they can vastly reduce the number of experiments needed to understand how E. coli responds to the combinatorially huge space of possible environments.

Bacteria face complex conditions in important settings such as our body and in biotechnological applications such as biofuel production. Understanding how bacteria respond to complex conditions is a hard problem: the number of conditions that need to be tested grows exponentially with the number of nutrients, stresses and other factors that make up the environment. To overcome this exponential explosion, we present an approach that allows computing the dynamics of gene expression in a complex condition based on measurements in simple conditions. This is based on the main discovery in this paper: using accurate promoter activity measurements, we find that promoter activity dynamics in a cocktail of media is a weighted average of the dynamics in each medium alone. The weights in the average are constant across time, and can be used to predict the dynamics in arbitrary cocktails based only on measurements on pairs of conditions. Thus, dynamics in complex conditions is, for the vast majority of genes, much simpler than it might have been; this simplicity allows new mathematical formula for accurate prediction in new conditions.

MicroRNA-7a regulates pancreatic β cell function

Dysfunctional microRNA (miRNA) networks contribute to inappropriate responses following pathological stress and are the underlying cause of several disease conditions. In pancreatic β cells, miRNAs have been largely unstudied and little is known about how specific miRNAs regulate glucose-stimulated insulin secretion (GSIS) or impact the adaptation of β cell function to metabolic stress. In this study, we determined that miR-7 is a negative regulator of GSIS in β cells. Using Mir7a2 deficient mice, we revealed that miR-7a2 regulates β cell function by directly regulating genes that control late stages of insulin granule fusion with the plasma membrane and ternary SNARE complex activity. Transgenic mice overexpressing miR-7a in β cells developed diabetes due to impaired insulin secretion and β cell dedifferentiation. Interestingly, perturbation of miR-7a expression in β cells did not affect proliferation and apoptosis, indicating that miR-7 is dispensable for the maintenance of endocrine β cell mass. Furthermore, we found that miR-7a levels are decreased in obese/diabetic mouse models and human islets from obese and moderately diabetic individuals with compensated β cell function. Our results reveal an interconnecting miR-7 genomic circuit that regulates insulin granule exocytosis in pancreatic β cells and support a role for miR-7 in the adaptation of pancreatic β cell function in obesity and type 2 diabetes.

MicroRNA binding sites in the coding region of mRNAs: extending the repertoire of post-transcriptional gene regulation

It is well established that microRNAs (miRNAs) induce mRNA degradation by binding to 3' untranslated regions (UTRs). The functionality of sites in the coding domain sequence (CDS), on the other hand, remains under discussion. Such sites have limited impact on target mRNA abundance and recent work suggests that miRNAs bind in the CDS to inhibit translation. What then could be the regulatory benefits of translation inhibition through CDS targeting compared to mRNA degradation following 3' UTR binding? We propose that these domain-dependent effects serve to diversify the functional repertoire of post-transcriptional gene expression control. Possible regulatory benefits may include tuning the time-scale and magnitude of post-transcriptional regulation, regulating protein abundance depending on or independently of the cellular state, and regulation of the protein abundance of alternative splice variants. Finally, we review emerging evidence that these ideas may generalize to RNA-binding proteins beyond miRNAs and Argonaute proteins.

Characterization of 13 polymorphic microsatellite loci for two land snail species, Trochulus villosus and T. sericeus (Gastropoda: Pulmonata: Hygromiidae)

Thirteen new microsatellite loci were isolated and tested on two land snail species, Trochulus villosus and T. sericeus (Pulmonata: Hygromiidae), resulting in a set of eight polymorphic markers for each species. The expected heterozygosity was high for all loci and species (between 0.616 and 0.944). Such levels of variability will allow detailed insights into the population genetic structure of some Trochulus species.

A biophysical miRNA-mRNA interaction model infers canonical and noncanonical targets

We introduce a biophysical model of miRNA-target interaction and infer its parameters from Argonaute 2 cross-linking and immunoprecipitation data. We show that a substantial fraction of human miRNA target sites are noncanonical and that predicted target-site affinity correlates well with the extent of target destabilization. Our model provides a rigorous biophysical approach to miRNA target identification beyond ad hoc miRNA seed-based methods.

Analysis of CDS-located miRNA target sites suggests that they can effectively inhibit translation

Most of what is presently known about how miRNAs regulate gene expression comes from studies that characterized the regulatory effect of miRNA binding sites located in the 3′ untranslated regions (UTR) of mRNAs. In recent years, there has been increasing evidence that miRNAs also bind in the coding region (CDS), but the implication of these interactions remains obscure because they have a smaller impact on mRNA stability compared with miRNA-target interactions that involve 3′ UTRs. Here we show that miRNA-complementary sites that are located in both CDS and 3′-UTRs are under selection pressure and share the same sequence and structure properties. Analyzing recently published data of ribosome-protected fragment profiles upon miRNA transfection from the perspective of the location of miRNA-complementary sites, we find that sites located in the CDS are most potent in inhibiting translation, while sites located in the 3′ UTR are more efficient at triggering mRNA degradation. Our study suggests that miRNAs may combine targeting of CDS and 3′ UTR to flexibly tune the time scale and magnitude of their post-transcriptional regulatory effects.

Argonaute CLIP--a method to identify in vivo targets of miRNAs

microRNAs are important regulators of gene expression that guide translational repression and degradation of target mRNAs. Only relatively few miRNA targets have been characterized, and computational prediction is hampered by the relatively small number of nucleotides that seem to be involved in target recognition. Argonaute (Ago) crosslinking and immunoprecipitation (CLIP) in combination with next-generation sequencing proved to be a successful method for identifying targets of endogenous cellular miRNAs on a transcriptome-wide scale. Here we review various approaches to Ago CLIP, describe in detail the PAR-CLIP method and provide an outline of the necessary computational analysis for identification of in vivo miRNA binding sites.

MicroRNA-194 is a target of transcription factor 1 (Tcf1, HNF1α) in adult liver and controls expression of frizzled-6

Transcription factor 1 (Tcf1; hepatocyte nuclear factor 1α [HNF1α]) is critical for hepatocyte development and function. Whether Tcf1 also regulates hepatic microRNAs (miRNAs) has not been investigated yet. Here we analyzed Tcf1-dependent miRNA expression in adult mice in which this transcription factor had been genetically deleted (Tcf1(-/-) ) using miRNA microarray analysis. The miR-192/-194 cluster was markedly down-regulated in liver of Tcf1(-/-) mice. MiR-192/-194 levels were also decreased in two other tissues that express Tcf1, kidney and small intestine, although to a lesser extent than in liver. In order to identify targets of miR-192/-194 in vivo we combined Affymetrix gene analysis of liver in which miR-192/-194 had been silenced or overexpressed, respectively, and tested regulated messenger RNAs (mRNAs) with multiple binding sites for these miRNAs. This approach revealed frizzled-6 (Fzd6) as a robust endogenous target of miR-194. MiR-194 also targets human FZD6 and expression of miR-194 and Fzd6 are inversely correlated in a mouse model of hepatocellular carcinoma (Dgcr8(flox/flox) p53(flox/flox) × Alb-Cre).

CONCLUSION

Our results support a role of miR-194 in liver tumorigenesis through its endogenous target Fzd6. These results may have important implications for Tcf1-mediated liver proliferation.

Kaposi's sarcoma herpesvirus microRNAs target caspase 3 and regulate apoptosis

Kaposi's sarcoma herpesvirus (KSHV) encodes a cluster of twelve micro (mi)RNAs, which are abundantly expressed during both latent and lytic infection. Previous studies reported that KSHV is able to inhibit apoptosis during latent infection; we thus tested the involvement of viral miRNAs in this process. We found that both HEK293 epithelial cells and DG75 cells stably expressing KSHV miRNAs were protected from apoptosis. Potential cellular targets that were significantly down-regulated upon KSHV miRNAs expression were identified by microarray profiling. Among them, we validated by luciferase reporter assays, quantitative PCR and western blotting caspase 3 (Casp3), a critical factor for the control of apoptosis. Using site-directed mutagenesis, we found that three KSHV miRNAs, miR-K12-1, 3 and 4-3p, were responsible for the targeting of Casp3. Specific inhibition of these miRNAs in KSHV-infected cells resulted in increased expression levels of endogenous Casp3 and enhanced apoptosis. Altogether, our results suggest that KSHV miRNAs directly participate in the previously reported inhibition of apoptosis by the virus, and are thus likely to play a role in KSHV-induced oncogenesis.

MicroRNAs 103 and 107 regulate insulin sensitivity

Defects in insulin signalling are among the most common and earliest defects that predispose an individual to the development of type 2 diabetes. MicroRNAs have been identified as a new class of regulatory molecules that influence many biological functions, including metabolism. However, the direct regulation of insulin sensitivity by microRNAs in vivo has not been demonstrated. Here we show that the expression of microRNAs 103 and 107 (miR-103/107) is upregulated in obese mice. Silencing of miR-103/107 leads to improved glucose homeostasis and insulin sensitivity. In contrast, gain of miR-103/107 function in either liver or fat is sufficient to induce impaired glucose homeostasis. We identify caveolin-1, a critical regulator of the insulin receptor, as a direct target gene of miR-103/107. We demonstrate that caveolin-1 is upregulated upon miR-103/107 inactivation in adipocytes and that this is concomitant with stabilization of the insulin receptor, enhanced insulin signalling, decreased adipocyte size and enhanced insulin-stimulated glucose uptake. These findings demonstrate the central importance of miR-103/107 to insulin sensitivity and identify a new target for the treatment of type 2 diabetes and obesity.

PAR-CliP--a method to identify transcriptome-wide the binding sites of RNA binding proteins

RNA transcripts are subjected to post-transcriptional gene regulation by interacting with hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) that are often expressed in a cell-type dependently. To understand how the interplay of these RNA-binding factors affects the regulation of individual transcripts, high resolution maps of in vivo protein-RNA interactions are necessary. A combination of genetic, biochemical and computational approaches are typically applied to identify RNA-RBP or RNA-RNP interactions. Microarray profiling of RNAs associated with immunopurified RBPs (RIP-Chip) defines targets at a transcriptome level, but its application is limited to the characterization of kinetically stable interactions and only in rare cases allows to identify the RBP recognition element (RRE) within the long target RNA. More direct RBP target site information is obtained by combining in vivo UV crosslinking with immunoprecipitation followed by the isolation of crosslinked RNA segments and cDNA sequencing (CLIP). CLIP was used to identify targets of a number of RBPs. However, CLIP is limited by the low efficiency of UV 254 nm RNA-protein crosslinking, and the location of the crosslink is not readily identifiable within the sequenced crosslinked fragments, making it difficult to separate UV-crosslinked target RNA segments from background non-crosslinked RNA fragments also present in the sample. We developed a powerful cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs that we term PAR-CliP (Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation) (see Fig. 1A for an outline of the method). The method relies on the incorporation of photoreactive ribonucleoside analogs, such as 4-thiouridine (4-SU) and 6-thioguanosine (6-SG) into nascent RNA transcripts by living cells. Irradiation of the cells by UV light of 365 nm induces efficient crosslinking of photoreactive nucleoside-labeled cellular RNAs to interacting RBPs. Immunoprecipitation of the RBP of interest is followed by isolation of the crosslinked and coimmunoprecipitated RNA. The isolated RNA is converted into a cDNA library and deep sequenced using Solexa technology. One characteristic feature of cDNA libraries prepared by PAR-CliP is that the precise position of crosslinking can be identified by mutations residing in the sequenced cDNA. When using 4-SU, crosslinked sequences thymidine to cytidine transition, whereas using 6-SG results in guanosine to adenosine mutations. The presence of the mutations in crosslinked sequences makes it possible to separate them from the background of sequences derived from abundant cellular RNAs. Application of the method to a number of diverse RNA binding proteins was reported in Hafner et al.

Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP

RNA transcripts are subject to posttranscriptional gene regulation involving hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) expressed in a cell-type dependent fashion. We developed a cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs. The crosslinked sites are revealed by thymidine to cytidine transitions in the cDNAs prepared from immunopurified RNPs of 4-thiouridine-treated cells. We determined the binding sites and regulatory consequences for several intensely studied RBPs and miRNPs, including PUM2, QKI, IGF2BP1-3, AGO/EIF2C1-4 and TNRC6A-C. Our study revealed that these factors bind thousands of sites containing defined sequence motifs and have distinct preferences for exonic versus intronic or coding versus untranslated transcript regions. The precise mapping of binding sites across the transcriptome will be critical to the interpretation of the rapidly emerging data on genetic variation between individuals and how these variations contribute to complex genetic diseases.

Relative contribution of sequence and structure features to the mRNA binding of Argonaute/EIF2C-miRNA complexes and the degradation of miRNA targets

How miRNAs recognize their target sites is a puzzle that many experimental and computational studies aimed to solve. Several features, such as perfect pairing of the miRNA seed, additional pairing in the 3' region of the miRNA, relative position in the 3' UTR, and the A/U content of the environment of the putative site, have been found to be relevant. Here we have used a large number of previously published data sets to assess the power that various sequence and structure features have in distinguishing between putative sites that do and those that do not appear to be functional. We found that although different data sets give widely different answers when it comes to ranking the relative importance of these features, the sites inferred from most transcriptomics experiments, as well as from comparative genomics, appear similar at this level. This suggests that miRNA target sites have been selected in evolution on their ability to trigger mRNA degradation. To understand at what step in the miRNA-induced response individual features play a role, we transfected human HEK293 cells with miRNAs and analyzed the association of Argonaute/EIF2C-miRNA complexes with target mRNAs and the degradation of these messages. We found that structural features of the target site are only important for Argonaute/EIF2C binding, while sequence features such as the A/U content of the 3' UTR are important for mRNA degradation.

MirZ: an integrated microRNA expression atlas and target prediction resource

MicroRNAs (miRNAs) are short RNAs that act as guides for the degradation and translational repression of protein-coding mRNAs. A large body of work showed that miRNAs are involved in the regulation of a broad range of biological functions, from development to cardiac and immune system function, to metabolism, to cancer. For most of the over 500 miRNAs that are encoded in the human genome the functions still remain to be uncovered. Identifying miRNAs whose expression changes between cell types or between normal and pathological conditions is an important step towards characterizing their function as is the prediction of mRNAs that could be targeted by these miRNAs. To provide the community the possibility of exploring interactively miRNA expression patterns and the candidate targets of miRNAs in an integrated environment, we developed the MirZ web server, which is accessible at www.mirz.unibas.ch. The server provides experimental and computational biologists with statistical analysis and data mining tools operating on up-to-date databases of sequencing-based miRNA expression profiles and of predicted miRNA target sites in species ranging from Caenorhabditis elegans to Homo sapiens.

Chromosomal rearrangements and gene flow over time in an inter-specific hybrid zone of the Sorex araneus group

Most hybrid zones have existed for hundreds or thousands of years but have generally been observed for only a short time period. Studies extending over periods long enough to track evolutionary changes in the zones or assess the ultimate outcome of hybridization are scarce. Here, we describe the evolution over time of the level of genetic isolation between two karyotypically different species of shrews (Sorex araneus and Sorex antinorii) at a hybrid zone located in the Swiss Alps. We first evaluated hybrid zone movement by contrasting patterns of gene flow and changes in cline parameters (centre and width) using 24 microsatellite loci, between two periods separated by 10 years apart. Additionally, we tested the role of chromosomal rearrangements on gene flow by analysing microsatellite loci located on both rearranged and common chromosomes to both species. We did not detect any movement of the hybrid zone during the period analysed, suggesting that the zone is a typical tension zone. However, the gene flow was significantly lower among the rearranged than the common chromosomes for the second period, whereas the difference was only marginally significant for the first period. This further supports the role of chromosomal rearrangements on gene flow between these taxa.

Testis size, sperm characteristics and testosterone concentrations in four species of shrews (Mammalia, Soricidae)

The aim of this study was to establish and compare the sperm characteristics in four shrew species in the context of the sperm competition hypothesis. As expected, the large relative testis size in promiscuous species was associated with a high number of cauda epididymal spermatozoa and a high concentration of circulating testosterone. In addition, in Sorex and Neomys, species with high intensity of sperm competition, the spermatozoa stored in cauda epididymis were characterized by high percentage of progressive motility whereas in Crocidura and Suncus, the cauda epididymal spermatozoa were motile but with very low percentage of progressive motility. This capability is achieved only following the passage through the vas gland, a specialized region for sperm storage located along the vas deferens in these shrew species. The hypothesis that sperm competition is positively correlated with spermatozoa length could not be confirmed. In Crocidura and Suncus, the total sperm length is increased by the large sperm head due to a big acrosome. This trait, specific to the subfamily Crocidurinae, may results from a selective pressure independent of the context of sperm competition, related to a specific, but as yet unclear role, for the acrosome during the fertilization.

Chromosomal rearrangements and genetic structure at different evolutionary levels of the Sorex araneus group

Robertsonian (Rb) fusions received large theoretical support for their role in speciation, but empirical evidence is often lacking. Here, we address the role of Rb rearrangements on the genetic differentiation of the karyotypically diversified group of shrews, Sorex araneus. We compared genetic structure between 'rearranged' and 'common' chromosomes in pairwise comparisons of five karyotypic taxa of the group. Considering all possible comparisons, we found a significantly greater differentiation at rearranged chromosomes, supporting the role of chromosomal rearrangements in the general genetic diversification of this group. Intertaxa structure and distance were larger across rearranged chromosomes for most of the comparisons, although these differences were not significant. This last result could be explained by the large variance observed among microsatellite-based estimates. The differences observed among the pairs of taxa analysed support the role of both the hybrid karyotypic complexity and the level of evolutionary divergence.

A new perspective on the evolutionary history of western European Sorex araneus group revealed by paternal and maternal molecular markers

The species of the common shrew (Sorex araneus) group are morphologically very similar, but have undergone a spectacular chromosomal evolution. We investigate here the evolutionary history of the Sorex araneus group distributed in western Europe. In particular, we clarify the position of a difficult species, S. granarius, using sex-specific (mtDNA and Y-chromosome) markers. The karyotype of S. granarius is generally considered similar to the common ancestor of the restricted group considered here. The mtDNA data (1.4 kb) confirms the close relationship between S. granarius and S. araneus sensu stricto (hereafter S. araneus s.s.), but the Y-chromosome (3.4 kb) produces a quite different picture: S. granarius is closely related to another species, S. coronatus. Comparison of mtDNA and Y-chromosome phylogenies suggests that the genetic and chromosomal evolution in this group are disconnected processes. The evolutionary history of the south-western European populations of the S. araneus group can only be understood considering secondary contacts between taxa after their divergence, implying genetic exchanges by means of hybridization and/or introgression.

[Drowning]

Drowning is always a serious event, with high morbidity and mortality. Definitions were recently revised and simplified according to Utstein-style. The circumstances and prognosis differ in different age groups. Consequences of drowning are hypoxemia, acidosis and hypothermia. The management of drowning victims should be directed to the correction of these 3 elements. An early basic life support and an efficient oxygenation are keys of good prognosis, this one being mainly related to hypoxic cerebral damages. The most frequent complications except anoxic encephalopathy are bacterial pneumonia and ARDS.

Using a Bayesian method to assign individuals to karyotypic taxa in shrew hybrid zones

Individuals sampled in hybrid zones are usually analysed according to their sampling locality, morphology, behaviour or karyotype. But the increasing availability of genetic information more and more favours its use for individual sorting purposes and numerous assignment methods based on the genetic composition of individuals have been developed. The shrews of the Sorex araneus group offer good opportunities to test the genetic assignment on individuals identified by their karyotype. Here we explored the potential and efficiency of a Bayesian assignment method combined or not with a reference dataset to study admixture and individual assignment in the difficult context of two hybrid zones between karyotypic species of the Sorex araneus group. As a whole, we assigned more than 80% of the individuals to their respective karyotypic categories (i.e. 'pure' species or hybrids). This assignment level is comparable to what was obtained for the same species away from hybrid zones. Additionally, we showed that the assignment result for several individuals was strongly affected by the inclusion or not of a reference dataset. This highlights the importance of such comparisons when analysing hybrid zones. Finally, differences between the admixture levels detected in both hybrid zones support the hypothesis of an impact of chromosomal rearrangements on gene flow.

Genetic and karyotypic structure in the shrews of the Sorex araneus group: are they independent?

The species of the common shrew (Sorex araneus) group are morphologically very similar but exhibit high levels of karyotypic variation. Here we used genetic variation at 10 microsatellite markers in a data set of 212 individuals mostly sampled in the western Alps and composed of five karyotypic taxa (Sorex coronatus, Sorex antinorii and the S. araneus chromosome races Cordon, Bretolet and Vaud) to investigate the concordance between genetic and karyotypic structure. Bayesian analysis confirmed the taxonomic status of the three sampled species since individuals consistently grouped according to their taxonomical status. However, introgression can still be detected between S. antinorii and the race Cordon of S. araneus. This observation is consistent with the expected low karyotypic complexity of hybrids between these two taxa. Geographically based cryptic substructure was discovered within S. antinorii, a pattern consistent with the different postglaciation recolonization routes of this species. Additionally, we detected two genetic groups within S. araneus notwithstanding the presence of three chromosome races. This pattern can be explained by the probable hybrid status of the Bretolet race but also suggests a relatively low impact of chromosomal differences on genetic structure compared to historical factors. Finally, we propose that the current data set (available at http://www.unil.ch/dee/page7010_en.html#1) could be used as a reference by those wanting to identify Sorex individuals sampled in the western Alps.

Funktionelle Ergebnisse nach autologer Rotationskeratoplastik

BACKGROUND

Autologous rotational keratoplasty is a therapeutic option in patients with central corneal scars. Its advantage compared to allogeneic procedures is the absence of immunologic destruction. In this retrospective analysis, long-term results are presented.

PATIENTS

Since 1988, autologous rotational keratoplasty was performed in 19 patients. The trephine size varied between 7.0 mm and 9.0 mm. Mean follow-up was 7.1+/-4.6 years.

RESULTS

In 16 of the patients, the graft center remained clear during the entire follow-up period. Visual acuity improved from 0.25+/-0.16 preoperatively to 0.34+/-0.27 postoperatively with spectacles and to 0.43+/-0.3 with contact lenses. The mean postoperative astigmatism was 6.41+/-3.87 dpt with irregular components in six patients.

CONCLUSION

Only a few postoperative graft complications were observed. Especially in patients with a markedly elevated risk for immune reactions (e.g., small infants), autologous rotational keratoplasty is a safe alternative to allogenic procedures.

Alps, genes, and chromosomes: their role in the formation of species in the Sorex araneus group (Mammalia, Insectivora), as inferred from two hybrid zones

During the Pleistocene glaciations, the Alps were an efficient barrier to gene flow between isolated populations, often leading to allopatric speciation. Afterwards, the Alps strongly influenced the post-glacial recolonization of Europe and represent a major suture zone between differentiated populations. Two hybrid zones in the Swiss and French Alps between genetically and chromosomally well-differentiated species-the Valais shrew, Sorex antinorii, and the common shrew, S. araneus-were studied karyotypically and by analyzing the distribution of seven microsatellite loci. In the center of the Haslital hybrid zone the two species coexist over a distance of 900 m. Hybrid karyotypes, among them the most complex known in Sorex, are rare. F-statistics based on microsatellite data revealed a strong heterozygote deficit only in the center of the zone, due to the sympatric distribution of the two species with little hybridization between them. Structuring within the species (both F(IS) and F(ST)) was low. An hierarchical analysis showed a high level of interspecific differentiation. Results were compared with those previously reported in another hybrid zone located at Les Houches in the French Alps. Genetic structuring within and between species was comparable in both hybrid zones, although chromosomal incompatibilities are more important in Haslital, where a linkage block of the race-specific chromosomes should additionally impede gene flow. Evidence for a more restricted gene flow in Haslital comes from the genetically intermediate hybrid karyotypes, whereas in Les Houches, hybrid karyotypes are genetically identical to individuals of the pure karyotypic races. Genic and chromosomal introgression was observed in Les Houches, but not in Haslital. The possible influence of a river, separating the two species at Les Houches, on gene flow is discussed.

Phylogeographical structure, postglacial recolonization and barriers to gene flow in the distinctive Valais chromosome race of the common shrew (Sorex araneus)

Using one male-inherited and eight biparentally inherited microsatellite markers, we investigate the population genetic structure of the Valais chromosome race of the common shrew (Sorex araneus) in the Central Alps of Europe. Unexpectedly, the Y-chromosome microsatellite suggests nearly complete absence of male gene flow among populations from the St-Bernard and Simplon regions (Switzerland). Autosomal markers also show significant genetic structuring among these two geographical areas. Isolation by distance is significant and possible barriers to gene flow exist in the study area. Two different approaches are used to better understand the geographical patterns and the causes of this structuring. Using a principal component analysis for which testing procedure exists, and partial Mantel tests, we show that the St-Bernard pass does not represent a significant barrier to gene flow although it culminates at 2469 m, close to the highest altitudinal record for this species. Similar results are found for the Simplon pass, indicating that both passes represented potential postglacial recolonization routes into Switzerland from Italian refugia after the last Pleistocene glaciations. In contrast with the weak effect of these mountain passes, the Rhône valley lowlands significantly reduce gene flow in this species. Natural obstacles (the large Rhône river) and unsuitable habitats (dry slopes) are both present in the valley. Moreover, anthropogenic changes to landscape structures are likely to have strongly reduced available habitats for this shrew in the lowlands, thereby promoting genetic differentiation of populations found on opposite sides of the Rhône valley.

Is the Gibraltar strait a barrier to gene flow for the bat Myotis myotis (Chiroptera: Vespertilionidae)?

Because of their role in limiting gene flow, geographical barriers like mountains or seas often coincide with intraspecific genetic discontinuities. Although the Strait of Gibraltar represents such a potential barrier for both plants and animals, few studies have been conducted on its impact on gene flow. Here we test this effect on a bat species (Myotis myotis) which is apparently distributed on both sides of the strait. Six colonies of 20 Myotis myotis each were sampled in southern Spain and northern Morocco along a linear transect of 1350 km. Results based on six nuclear microsatellite loci reveal no significant population structure within regions, but a complete isolation between bats sampled on each side of the strait. Variability at 600 bp of a mitochondrial gene (cytochrome b) confirms the existence of two genetically distinct and perfectly segregating clades, which diverged several million years ago. Despite the narrowness of the Gibraltar Strait (14 km), these molecular data suggest that neither males, nor females from either region have ever reproduced on the opposite side of the strait. Comparisons of molecular divergence with bats from a closely related species (M. blythii) suggest that the North African clade is possibly a distinct taxon warranting full species rank. We provisionally refer to it as Myotis cf punicus Felten 1977, but a definitive systematic understanding of the whole Mouse-eared bat species complex awaits further genetic sampling, especially in the Eastern Mediterranean areas.

Meiotic drive favors Robertsonian metacentric chromosomes in the common shrew (Sorex araneus, Insectivora, mammalia)

Meiotic drive has attracted much interest because it concerns the robustness of Mendelian segregation and its genetic and evolutionary stability. We studied chromosomal meiotic drive in the common shrew (Sorex araneus, Insectivora, Mammalia), which exhibits one of the most remarkable chromosomal polymorphisms within mammalian species. The open question of the evolutionary success of metacentric chromosomes (Robertsonian fusions) versus acrocentrics in the common shrew prompted us to test whether a segregation distortion in favor of metacentrics is present in female and/or male meiosis. Performing crosses under controlled laboratory conditions with animals from natural populations, we found a clear trend toward a segregation distortion in favor of metacentrics during male meiosis, two chromosome combinations (gm and jl) being significantly preferred over their acrocentric homologs. Apart for one Robertsonian fusion (hi), this trend was absent in female meiosis. We propose a model based on recombination events between twin acrocentrics to explain the difference in transmission ratios of the same metacentric in different sexes and unequal drive of particular metacentrics in the same sex. Pooled data for female and male meiosis revealed a trend toward stronger segregation distortion for larger metacentrics. This is partially in agreement with the frequency of metacentrics occurring in natural populations of a chromosome race showing a high degree of chromosomal polymorphism.

Microsatellites can be misleading: an empirical and simulation study

It has been long recognized that highly polymorphic genetic markers can lead to underestimation of divergence between populations when migration is low. Microsatellite loci, which are characterized by extremely high mutation rates, are particularly likely to be affected. Here, we report genetic differentiation estimates in a contact zone between two chromosome races of the common shrew (Sorex araneus), based on 10 autosomal microsatellites, a newly developed Y-chromosome microsatellite, and mitochondrial DNA. These results are compared to previous data on proteins and karyotypes. Estimates of genetic differentiation based on F- and R-statistics are much lower for autosomal microsatellites than for all other genetic markers. We show by simulations that this discrepancy stems mainly from the high mutation rate of microsatellite markers for F-statistics and from deviations from a single-step mutation model for R-statistics. The sex-linked genetic markers show that all gene exchange between races is mediated by females. The absence of male-mediated gene flow most likely results from male hybrid sterility.

Do riverine barriers, history or introgression shape the genetic structuring of a common shrew (Sorex araneus) population?

The common shrew (Sorex araneus) is subdivided into numerous chromosome races. The Valais and Cordon chromosome races meet and hybridize at a mountain river in Les Houches (French Alps). Significant genetic structuring was recently reported among populations found on the Valais side of this hybrid zone. In this paper, a phylogenetic analysis and partial Mantel tests are used to investigate the patterns and causes of this structuring. A total of 185 shrews were trapped at 12 localities. All individuals were typed for nine microsatellite loci. Although several mountain rivers are found in the study area, riverine barriers do not have a significant influence on gene flow. Partial Mantel tests show that our result is caused by the influence of the hybrid zone with the Cordon race. The geographical patterns of this structuring are discussed in the context of the contact zone, which appears to extend up to a group of two rivers. The glacier they originate from is known to have cut the Arve valley as recently as 1818. The recent history of this glacier, its moraine and possibly rivers, may therefore be linked to the history of this hybrid zone.

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.