Large scale assessment of regulatory evolution and transcriptome complexity in mammals

AbstractIn addition to genetic changes affecting the function of gene products, changes in gene expression have been suggested to underlie many or even most of the phenotypic differences among mammals. However, detailed gene expression comparisons were, until recently, restricted to closely related species, owing to technological limitations. Thus, we took advantage of the latest technologies (RNA-Seq) to generate extensive qualitative and quantitative transcriptome data for a unique collection of somatic and germline tissues from representatives of all major mammalian lineages (placental mammals, marsupials and monotremes) and birds, the evolutionary outgroup.In the first major project of my thesis, we performed global comparative analyses of gene expression levels based on these data. Our analyses provided fundamental insights into the dynamics of transcriptome change during mammalian evolution (e.g., the rate of expression change across species, tissues and chromosomes) and allowed the exploration of the functional relevance and phenotypic implications of transcription changes at a genome-wide scale (e.g., we identified numerous potentially selectively driven expression switches).In a second project of my thesis, which was also based on the unique transcriptome data generated in the context of the first project we focused on the evolution of alternative splicing in mammals. Alternative splicing contributes to transcriptome complexity by generating several transcript isoforms from a single gene, which can, thus, perform various functions. To complete the global comparative analysis of gene expression changes, we explored patterns of alternative splicing evolution. This work uncovered several general and unexpected patterns of alternative splicing evolution (e.g., we found that alternative splicing evolves extremely rapidly) as well as a large number of conserved alternative isoforms that may be crucial for the functioning of mammalian organs.Finally, the third and final project of my PhD consisted in analyzing in detail the unique functional and evolutionary properties of the testis by exploring the extent of its transcriptome complexity. This organ was previously shown to evolve rapidly both at the phenotypic and molecular level, apparently because of the specific pressures that act on this organ and are associated with its reproductive function. Moreover, my analyses of the amniote tissue transcriptome data described above, revealed strikingly widespread transcriptional activity of both functional and nonfunctional genomic elements in the testis compared to the other organs. To elucidate the cellular source and mechanisms underlying this promiscuous transcription in the testis, we generated deep coverage RNA-Seq data for all major testis cell types as well as epigenetic data (DNA and histone methylation) using the mouse as model system. The integration of these complete dataset revealed that meiotic and especially post-meiotic germ cells are the major contributors to the widespread functional and nonfunctional transcriptome complexity of the testis, and that this "promiscuous" spermatogenic transcription is resulting, at least partially, from an overall transcriptionally permissive chromatin state. We hypothesize that this particular open state of the chromatin results from the extensive chromatin remodeling that occurs during spermatogenesis which ultimately leads to the replacement of histones by protamines in the mature spermatozoa. Our results have important functional and evolutionary implications (e.g., regarding new gene birth and testicular gene expression evolution).Generally, these three large-scale projects of my thesis provide complete and massive datasets that constitute valuables resources for further functional and evolutionary analyses of mammalian genomes.

MicroRNAs shape circadian hepatic gene expression on a transcriptome-wide scale.

A considerable proportion of mammalian gene expression undergoes circadian oscillations. Post-transcriptional mechanisms likely make important contributions to mRNA abundance rhythms. We have investigated how microRNAs (miRNAs) contribute to core clock and clock-controlled gene expression using mice in which miRNA biogenesis can be inactivated in the liver. While the hepatic core clock was surprisingly resilient to miRNA loss, whole transcriptome sequencing uncovered widespread effects on clock output gene expression. Cyclic transcription paired with miRNA-mediated regulation was thus identified as a frequent phenomenon that affected up to 30% of the rhythmic transcriptome and served to post-transcriptionally adjust the phases and amplitudes of rhythmic mRNA accumulation. However, only few mRNA rhythms were actually generated by miRNAs. Overall, our study suggests that miRNAs function to adapt clock-driven gene expression to tissue-specific requirements. Finally, we pinpoint several miRNAs predicted to act as modulators of rhythmic transcripts, and identify rhythmic pathways particularly prone to miRNA regulation.DOI: http://dx.doi.org/10.7554/eLife.02510.001.

Magnetic susceptibility to identify landscape segments on a detailed scale in the region of Jaboticabal, São Paulo, Brazil

The agricultural potential is generally assessed and managed based on a one-dimensional vision of the soil profile, however, the increased appreciation of sustainable production has stimulated studies on faster and more accurate evaluation techniques and methods of the agricultural potential on detailed scales. The objective of this study was to investigate the possibility of using soil magnetic susceptibility for the identification of landscape segments on a detailed scale in the region of Jaboticabal, São Paulo State. The studied area has two slope curvatures: linear and concave, subdivided into three landscape segments: upper slope (US, concave), middle slope (MS, linear) and lower slope (LS, linear). In each of these segments, 20 points were randomly sampled from a database with 207 samples forming a regular grid installed in each landscape segment. The soil physical and chemical properties, CO2 emissions (FCO2) and magnetic susceptibility (MS) of the samples were evaluated represented by: magnetic susceptibility of air-dried fine earth (MS ADFE), magnetic susceptibility of the total sand fraction (MS TS) and magnetic susceptibility of the clay fraction (MS Cl) in the 0.00 - 0.15 m layer. The principal component analysis showed that MS is an important property that can be used to identify landscape segments, because the correlation of this property within the first principal component was high. The hierarchical cluster analysis method identified two groups based on the variables selected by principal component analysis; of the six selected variables, three were related to magnetic susceptibility. The landscape segments were differentiated similarly by the principal component analysis and by the cluster analysis using only the properties with higher discriminatory power. The cluster analysis of MS ADFE, MS TS and MS Cl allowed the formation of three groups that agree with the segment division established in the field. The grouping by cluster analysis indicated MS as a tool that could facilitate the identification of landscape segments and enable the mapping of more homogeneous areas at similar locations.

A large-scale genetic validation study coupled with in vitro analyses reveal a role for vitamin Dsignaling in the pathogenesis and response to treatment of hepatitis C virus infection

Background and Aims: Vitamin D is an important modulatorof numerous cellular processes. Some of us recently observedan association of the 1a-hydroxylase promoter polymorphismCYP27B1-1260 rs10877012 with sustained virologic response (SVR)in a relatively small number of German patients with chronichepatitis C. In the present study, we aimed to validate thisassociation in a large and well characterized patient cohort, theSwiss Hepatitis C Cohort Study (SCCS). In addition, we examinedthe effect of vitamin D on the hepatitis C virus (HCV) life cyclein vitro.Methods: CYP27B1-1260 rs10877012 and IL28B rs12979860 singlenucleotide polymorphisms (SNPs) were genotyped in 1049 patientswith chronic hepatitis C from the SCCS, of whom 698 were treatedwith pegylated interferon-a (PEG-IFN-a) and ribavirin. In addition,112 patients with spontaneous clearance of HCV were examined.SNPs were correlated with variables reflecting the natural courseand treatment outcome of chronic hepatitis C. The effect of1,25-(OH)2D3 (calcitriol) on HCV replication and viral particleproduction was investigated in vitro using human hepatoma celllines (Huh-7.5) harbouring subgenomic replicons and cell culturederivedHCV.Results: The CYP27B1-1260 rs10877012 genotype was notassociated with SVR in patients with the good-response IL28Brs1279860 CC genotype. However, in patients with poor-responseIL28B rs1279860 genotype CT and TT, CYP27B1-1260 rs10877012was a significant independent predictor of SVR (15% difference inSVR between rs10877012 genotype AA vs. CC, p = 0.030, OR = 1.495,95% CI = 1.038-2.152). The CYPB27-1260 rs10877012 genotype wasneither associated with spontaneous clearance of HCV, nor withliver fibrosis progression rate, inflammatory activity of chronichepatitis C, or HCV viral load. Physiological doses of 1,25-(OH)2D3did not significantly affect HCVRNA replication or infectiousparticle production in vitro.Conclusions: The results of this large-scale genetic validationstudy reveal a role of vitamin D metabolism in the responseto treatment in chronic hepatitis C, but 1,25-(OH)2D3 does notexhibit a significant direct inhibitory antiviral effect. Thus, theability of vitamin D to modulate immunity against HCV shouldbe investigated.

Hydrological processes obtained on the plot scale under four simulated rainfall tests during the cycle of different crop systems

The cropping system influences the interception of water by plants, water storage in depressions on the soil surface, water infiltration into the soil and runoff. The aim of this study was to quantify some hydrological processes under no tillage cropping systems at the edge of a slope, in 2009 and 2010, in a Humic Dystrudept soil, with the following treatments: corn, soybeans, and common beans alone; and intercropped corn and common bean. Treatments consisted of four simulated rainfall tests at different times, with a planned intensity of 64 mm h-1 and 90 min duration. The first test was applied 18 days after sowing, and the others at 39, 75 and 120 days after the first test. Different times of the simulated rainfall and stages of the crop cycle affected soil water content prior to the rain, and the time runoff began and its peak flow and, thus, the surface hydrological processes. The depth of the runoff and the depth of the water intercepted by the crop + soil infiltration + soil surface storage were affected by the crop systems and the rainfall applied at different times. The corn crop was the most effective treatment for controlling runoff, with a water loss ratio of 0.38, equivalent to 75 % of the water loss ratio exhibited by common bean (0.51), the least effective treatment in relation to the others. Total water loss by runoff decreased linearly with an increase in the time that runoff began, regardless of the treatment; however, soil water content on the gravimetric basis increased linearly from the beginning to the end of the rainfall.

Absence of epidemic threshold in scale free networks with degree correlations

Random scale-free networks have the peculiar property of being prone to the spreading of infections. Here we provide for the susceptible-infected-susceptible model an exact result showing that a scale-free degree distribution with diverging second moment is a sufficient condition to have null epidemic threshold in unstructured networks with either assortative or disassortative mixing. Degree correlations result therefore irrelevant for the epidemic spreading picture in these scale-free networks. The present result is related to the divergence of the average nearest neighbors degree, enforced by the degree detailed balance condition.

Generation of uncorrelated random scale-free networks

Uncorrelated random scale-free networks are useful null models to check the accuracy and the analytical solutions of dynamical processes defined on complex networks. We propose and analyze a model capable of generating random uncorrelated scale-free networks with no multiple and self-connections. The model is based on the classical configuration model, with an additional restriction on the maximum possible degree of the vertices. We check numerically that the proposed model indeed generates scale-free networks with no two- and three-vertex correlations, as measured by the average degree of the nearest neighbors and the clustering coefficient of the vertices of degree k, respectively.

Scale decisions can reverse conclusions on community assembly processes.

AIM: Phylogenetic diversity patterns are increasingly being used to better understand the role of ecological and evolutionary processes in community assembly. Here, we quantify how these patterns are influenced by scale choices in terms of spatial and environmental extent and organismic scales. LOCATION: European Alps. METHODS: We applied 42 sampling strategies differing in their combination of focal scales. For each resulting sub-dataset, we estimated the phylogenetic diversity of the species pools, phylogenetic α-diversities of local communities, and statistics commonly used together with null models in order to infer non-random diversity patterns (i.e. phylogenetic clustering versus over-dispersion). Finally, we studied the effects of scale choices on these measures using regression analyses. RESULTS: Scale choices were decisive for revealing signals in diversity patterns. Notably, changes in focal scales sometimes reversed a pattern of over-dispersion into clustering. Organismic scale had a stronger effect than spatial and environmental extent. However, we did not find general rules for the direction of change from over-dispersion to clustering with changing scales. Importantly, these scale issues had only a weak influence when focusing on regional diversity patterns that change along abiotic gradients. MAIN CONCLUSIONS: Our results call for caution when combining phylogenetic data with distributional data to study how and why communities differ from random expectations of phylogenetic relatedness. These analyses seem to be robust when the focus is on relating community diversity patterns to variation in habitat conditions, such as abiotic gradients. However, if the focus is on identifying relevant assembly rules for local communities, the uncertainty arising from a certain scale choice can be immense. In the latter case, it becomes necessary to test whether emerging patterns are robust to alternative scale choices.

Density-based hierarchical clustering of pyro-sequences on a large scale--the case of fungal ITS1.

MOTIVATION: Analysis of millions of pyro-sequences is currently playing a crucial role in the advance of environmental microbiology. Taxonomy-independent, i.e. unsupervised, clustering of these sequences is essential for the definition of Operational Taxonomic Units. For this application, reproducibility and robustness should be the most sought after qualities, but have thus far largely been overlooked. RESULTS: More than 1 million hyper-variable internal transcribed spacer 1 (ITS1) sequences of fungal origin have been analyzed. The ITS1 sequences were first properly extracted from 454 reads using generalized profiles. Then, otupipe, cd-hit-454, ESPRIT-Tree and DBC454, a new algorithm presented here, were used to analyze the sequences. A numerical assay was developed to measure the reproducibility and robustness of these algorithms. DBC454 was the most robust, closely followed by ESPRIT-Tree. DBC454 features density-based hierarchical clustering, which complements the other methods by providing insights into the structure of the data. AVAILABILITY: An executable is freely available for non-commercial users at ftp://ftp.vital-it.ch/tools/dbc454. It is designed to run under MPI on a cluster of 64-bit Linux machines running Red Hat 4.x, or on a multi-core OSX system. CONTACT: dbc454@vital-it.ch or nicolas.guex@isb-sib.ch.

Glasgow Coma Score für den Patienten mit Schädel-Hirn-Trauma [Glasgow Coma Scale in traumatic brain injury].

Even 30 years after its first publication the Glasgow Coma Scale (GCS) is still used worldwide to describe and assess coma. The GCS consists of three components, the ocular, motor and verbal response to standardized stimulation, and is used as a severity of illness indicator for coma of various origins. The GCS facilitates information transfer and monitoring changes in coma. In addition, it is used as a triage tool in patients with traumatic brain injury. Its prognostic value regarding the outcome after a traumatic brain injury still lacks evidence. One of the main problems is the evaluation of the GCS in sedated, paralysed and/or intubated patients. A multitude of pseudoscores exists but a universal definition has yet to be defined.