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.

Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Molecular phylogeny analysis of fiddler crabs: test of the hypothesis of increasing behavioral complexity in evolution.

The current phylogenetic hypothesis for the evolution and biogeography of fiddler crabs relies on the assumption that complex behavioral traits are assumed to also be evolutionary derived. Indo-west Pacific fiddler crabs have simpler reproductive social behavior and are more marine and were thought to be ancestral to the more behaviorally complex and more terrestrial American species. It was also hypothesized that the evolution of more complex social and reproductive behavior was associated with the colonization of the higher intertidal zones. Our phylogenetic analysis, based upon a set of independent molecular characters, however, demonstrates how widely entrenched ideas about evolution and biogeography led to a reasonable, but apparently incorrect, conclusion about the evolutionary trends within this pantropical group of crustaceans. Species bearing the set of "derived traits" are phylogenetically ancestral, suggesting an alternative evolutionary scenario: the evolution of reproductive behavioral complexity in fiddler crabs may have arisen multiple times during their evolution. The evolution of behavioral complexity may have arisen by coopting of a series of other adaptations for high intertidal living and antipredator escape. A calibration of rates of molecular evolution from populations on either side of the Isthmus of Panama suggest a sequence divergence rate for 16S rRNA of 0.9% per million years. The divergence between the ancestral clade and derived forms is estimated to be approximately 22 million years ago, whereas the divergence between the American and Indo-west Pacific is estimated to be approximately 17 million years ago.

Pyramidal neurons of granular prefrontal cortex of the galago: Complexity in evolution of the psychic cell in primates

Typically, cognitive abilities of humans have been attributed to their greatly expanded cortical mantle, granular prefrontal cortex (gPFC) in particular. Recently we have demonstrated systematic differences in microstructure of gPFC in different species. Specifically, pyramidal cells in adult human gPFC are considerably more spinous than those in the gPFC of the macaque monkey, which are more spinous than those in the gPFC of marmoset and owl monkeys. As most cortical dendritic spines receive at least one excitatory input, pyramidal cells in these different species putatively receive different numbers of inputs. These differences in the gPFC pyramidal cell phenotype may be of fundamental importance in determining the functional characteristics of prefrontal circuitry and hence the cognitive styles of the different species. However, it remains unknown as to why the gPFC pyramidal cell phenotype differs between species. Differences could be attributed to, among other things, brain size, relative size of gPFC, or the lineage to which the species belong. Here we investigated pyramidal cells in the dorsolateral gPFC of the prosimian galago to extend the basis for comparison. We found these cells to be less spinous than those in human, macaque, and marmoset. (c) 2005 Wiley-Liss, Inc.

Cellular source and mechanisms of high transcriptome complexity in the Mammalian testis.

Understanding the extent of genomic transcription and its functional relevance is a central goal in genomics research. However, detailed genome-wide investigations of transcriptome complexity in major mammalian organs have been scarce. Here, using extensive RNA-seq data, we show that transcription of the genome is substantially more widespread in the testis than in other organs across representative mammals. Furthermore, we reveal that meiotic spermatocytes and especially postmeiotic round spermatids have remarkably diverse transcriptomes, which explains the high transcriptome complexity of the testis as a whole. The widespread transcriptional activity in spermatocytes and spermatids encompasses protein-coding and long noncoding RNA genes but also poorly conserves intergenic sequences, suggesting that it may not be of immediate functional relevance. Rather, our analyses of genome-wide epigenetic data suggest that this prevalent transcription, which most likely promoted the birth of new genes during evolution, is facilitated by an overall permissive chromatin in these germ cells that results from extensive chromatin remodeling.

W3233 -- Networks and Complexity in Social Systems (Collective Dynamics Group ISERP -- Columbia University) Course Syllabus

The Networks and Complexity in Social Systems course commences with an overview of the nascent field of complex networks, dividing it into three related but distinct strands: Statistical description of large scale networks, viewed as static objects; the dynamic evolution of networks, where now the structure of the network is understood in terms of a growth process; and dynamical processes that take place on fixed networks; that is, "networked dynamical systems". (A fourth area of potential research ties all the previous three strands together under the rubric of co-evolution of networks and dynamics, but very little research has been done in this vein and so it is omitted.) The remainder of the course treats each of the three strands in greater detail, introducing technical knowledge as required, summarizing the research papers that have introduced the principal ideas, and pointing out directions for future development. With regard to networked dynamical systems, the course treats in detail the more specific topic of information propagation in networks, in part because this topic is of great relevance to social science, and in part because it has received the most attention in the literature to date.

Local dimension and finite time prediction in spatiotemporal chaotic systems

Predictability is related to the uncertainty in the outcome of future events during the evolution of the state of a system. The cluster weighted modeling (CWM) is interpreted as a tool to detect such an uncertainty and used it in spatially distributed systems. As such, the simple prediction algorithm in conjunction with the CWM forms a powerful set of methods to relate predictability and dimension.

Intermittent pluri-sink model and the emergence of complex heterogeneity patterns: a simple paradigm for explaining complexity in soil chemical distributions

The spatial complexity of the distribution of organic matter, chemicals, nutrients, pollutants has been demonstrated to have multifractal nature (Kravchenco et al. [1]). This fact supports the possibility of existence of some emergent heterogeneity structure built under the evolution of the system. The aim of this note is providing a consistent explanation to the mentioned results via an extremely simple model.

Slip complexity in earthquake fault models.

We summarize studies of earthquake fault models that give rise to slip complexities like those in natural earthquakes. For models of smooth faults between elastically deformable continua, it is critical that the friction laws involve a characteristic distance for slip weakening or evolution of surface state. That results in a finite nucleation size, or coherent slip patch size, h*. Models of smooth faults, using numerical cell size properly small compared to h*, show periodic response or complex and apparently chaotic histories of large events but have not been found to show small event complexity like the self-similar (power law) Gutenberg-Richter frequency-size statistics. This conclusion is supported in the present paper by fully inertial elastodynamic modeling of earthquake sequences. In contrast, some models of locally heterogeneous faults with quasi-independent fault segments, represented approximately by simulations with cell size larger than h* so that the model becomes "inherently discrete," do show small event complexity of the Gutenberg-Richter type. Models based on classical friction laws without a weakening length scale or for which the numerical procedure imposes an abrupt strength drop at the onset of slip have h* = 0 and hence always fall into the inherently discrete class. We suggest that the small-event complexity that some such models show will not survive regularization of the constitutive description, by inclusion of an appropriate length scale leading to a finite h*, and a corresponding reduction of numerical grid size.

Areal specialization of pyramidal cell structure in the visual cortex of the tree shrew: a new twist revealed in the evolution of cortical circuitry

Cortical pyramidal cells, while having a characteristic morphology, show marked phenotypic variation in primates. Differences have been reported in their size, branching structure and spine density between cortical areas. In particular, there is a systematic increase in the complexity of the structure of pyramidal cells with anterior progression through occipito-temporal cortical visual areas. These differences reflect area-specific specializations in cortical circuitry, which are believed to be important for visual processing. However, it remains unknown as to whether these regional specializations in pyramidal cell structure are restricted to primates. Here we investigated pyramidal cell structure in the visual cortex of the tree shrew, including the primary (V1), second (V2) and temporal dorsal (TD) areas. As in primates, there was a trend for more complex branching structure with anterior progression through visual areas in the tree shrew. However, contrary to the trend reported in primates, cells in the tree shrew tended to become smaller with anterior progression through V1, V2 and TD. In addition, pyramidal cells in V1 of the tree shrew are more than twice as spinous as those in primates. These data suggest that variables that shape the structure of adult cortical pyramidal cells differ among species.

Phylogenetic and ecological significance in the evolution of Cetacean tonal sounds

Cetaceans are aquatic mammals that rely primarily on sound for most daily tasks. A compendium of sounds is emitted for orientation, prey detection, and predator avoidance, and to communicate. Communicative sounds are among the most studied Cetacean signals, particularly those referred to as tonal sounds. Because tonal sounds have been studied especially well in social dolphins, it has been assumed these sounds evolved as a social adaptation. However, whistles have been reported in ‘solitary’ species and have been secondarily lost three times in social lineages. Clearly, therefore, it is necessary to examine closely the association, if any, between whistles and sociality instead of merely assuming it. Several hypotheses have been proposed to explain the evolutionary history of Cetacean tonal sounds. The main goal of this dissertation is to cast light on the evolutionary history of tonal sounds by testing these hypotheses by combining comparative phylogenetic and field methods. This dissertation provides the first species-level phylogeny of Cetacea and phylogenetic tests of evolutionary hypotheses of cetacean communicative signals. Tonal sounds evolution is complex in that has likely been shaped by a combination of factors that may influence different aspects of their acoustical structure. At the inter-specific level, these results suggest that only tonal sound minimum frequency is constrained by body size. Group size also influences tonal sound minimum frequency. Species that live in large groups tend to produce higher frequency tonal sounds. The evolutionary history of tonal sounds and sociality may be intertwined, but in a complex manner rejecting simplistic views such as the hypothesis that tonal sounds evolved ‘for’ social communication in dolphins. Levels of social and tonal sound complexity nevertheless correlate indicating the importance of tonal sounds in social communication. At the intraspecific level, tonal sound variation in frequency and temporal parameters may be product of genetic isolation and local levels of underwater noise. This dissertation provides one of the first insights into the evolution of Cetacean tonal sounds in a phylogenetic context, and points out key species where future studies would be valuable to enrich our understanding of other factors also playing a role in tonal sound evolution. ^

Phylogenetic and Ecological Significance in the Evolution of Cetacean Tonal Sounds

Cetaceans are aquatic mammals that rely primarily on sound for most daily tasks. A compendium of sounds is emitted for orientation, prey detection, and predator avoidance, and to communicate. Communicative sounds are among the most studied Cetacean signals, particularly those referred to as tonal sounds. Because tonal sounds have been studied especially well in social dolphins, it has been assumed these sounds evolved as a social adaptation. However, whistles have been reported in ‘solitary’ species and have been secondarily lost three times in social lineages. Clearly, therefore, it is necessary to examine closely the association, if any, between whistles and sociality instead of merely assuming it. Several hypotheses have been proposed to explain the evolutionary history of Cetacean tonal sounds. The main goal of this dissertation is to cast light on the evolutionary history of tonal sounds by testing these hypotheses by combining comparative phylogenetic and field methods. This dissertation provides the first species-level phylogeny of Cetacea and phylogenetic tests of evolutionary hypotheses of cetacean communicative signals. Tonal sounds evolution is complex in that has likely been shaped by a combination of factors that may influence different aspects of their acoustical structure. At the inter-specific level, these results suggest that only tonal sound minimum frequency is constrained by body size. Group size also influences tonal sound minimum frequency. Species that live in large groups tend to produce higher frequency tonal sounds. The evolutionary history of tonal sounds and sociality may be intertwined, but in a complex manner rejecting simplistic views such as the hypothesis that tonal sounds evolved ‘for’ social communication in dolphins. Levels of social and tonal sound complexity nevertheless correlate indicating the importance of tonal sounds in social communication. At the intraspecific level, tonal sound variation in frequency and temporal parameters may be product of genetic isolation and local levels of underwater noise. This dissertation provides one of the first insights into the evolution of Cetacean tonal sounds in a phylogenetic context, and points out key species where future studies would be valuable to enrich our understanding of other factors also playing a role in tonal sound evolution.

Physicochemical complexity in complex chemical systems

Self-replication and compartmentalization are two central properties thought to be essential for minimal life, and understanding how such processes interact in the emergence of complex reaction networks is crucial to exploring the development of complexity in chemistry and biology. Autocatalysis can emerge from multiple different mechanisms such as formation of an initiator, template self-replication and physical autocatalysis (where micelles formed from the reaction product solubilize the reactants, leading to higher local concentrations and therefore higher rates). Amphiphiles are also used in artificial life studies to create protocell models such as micelles, vesicles and oil-in-water droplets, and can increase reaction rates by encapsulation of reactants. So far, no template self-replicator exists which is capable of compartmentalization, or transferring this molecular scale phenomenon to micro or macro-scale assemblies. Here a system is demonstrated where an amphiphilic imine catalyses its own formation by joining a non-polar alkyl tail group with a polar carboxylic acid head group to form a template, which was shown to form reverse micelles by Dynamic Light Scattering (DLS). The kinetics of this system were investigated by 1H NMR spectroscopy, showing clearly that a template self-replication mechanism operates, though there was no evidence that the reverse micelles participated in physical autocatalysis. Active oil droplets, composed from a mixture of insoluble organic compounds in an aqueous sub-phase, can undergo processes such as division, self-propulsion and chemotaxis, and are studied as models for minimal cells, or protocells. Although in most cases the Marangoni effect is responsible for the forces on the droplet, the behaviour of the droplet depends heavily on the exact composition. Though theoretical models are able to calculate the forces on a droplet, to model a mixture of oils on an aqueous surface where compounds from the oil phase are dissolving and diffusing through the aqueous phase is beyond current computational capability. The behaviour of a droplet in an aqueous phase can only be discovered through experiment, though it is determined by the droplet's composition. By using an evolutionary algorithm and a liquid handling robot to conduct droplet experiments and decide which compositions to test next, entirely autonomously, the composition of the droplet becomes a chemical genome capable of evolution. The selection is carried out according to a fitness function, which ranks the formulation based on how well it conforms to the chosen fitness criteria (e.g. movement or division). Over successive generations, significant increases in fitness are achieved, and this increase is higher with more components (i.e. greater complexity). Other chemical processes such as chemiluminescence and gelation were investigated in active oil droplets, demonstrating the possibility of controlling chemical reactions by selective droplet fusion. Potential future applications for this might include combinatorial chemistry, or additional fitness goals for the genetic algorithm. Combining the self-replication and the droplet protocells research, it was demonstrated that the presence of the amphiphilic replicator lowers the interfacial tension between droplets of a reaction mixture in organic solution and the alkaline aqueous phase, causing them to divide. Periodic sampling by a liquid handling robot revealed that the extent of droplet fission increased as the reaction progressed, producing more individual protocells with increased self-replication. This demonstrates coupling of the molecular scale phenomenon of template self-replication to a macroscale physicochemical effect.

Prognostic Value Of Dna And Mrna E6/e7 Of Human Papillomavirus In The Evolution Of Cervical Intraepithelial Neoplasia Grade 2.

This study aimed at evaluating whether human papillomavirus (HPV) groups and E6/E7 mRNA of HPV 16, 18, 31, 33, and 45 are prognostic of cervical intraepithelial neoplasia (CIN) 2 outcome in women with a cervical smear showing a low-grade squamous intraepithelial lesion (LSIL). This cohort study included women with biopsy-confirmed CIN 2 who were followed up for 12 months, with cervical smear and colposcopy performed every three months. Women with a negative or low-risk HPV status showed 100% CIN 2 regression. The CIN 2 regression rates at the 12-month follow-up were 69.4% for women with alpha-9 HPV versus 91.7% for other HPV species or HPV-negative status (P < 0.05). For women with HPV 16, the CIN 2 regression rate at the 12-month follow-up was 61.4% versus 89.5% for other HPV types or HPV-negative status (P < 0.05). The CIN 2 regression rate was 68.3% for women who tested positive for HPV E6/E7 mRNA versus 82.0% for the negative results, but this difference was not statistically significant. The expectant management for women with biopsy-confirmed CIN 2 and previous cytological tests showing LSIL exhibited a very high rate of spontaneous regression. HPV 16 is associated with a higher CIN 2 progression rate than other HPV infections. HPV E6/E7 mRNA is not a prognostic marker of the CIN 2 clinical outcome, although this analysis cannot be considered conclusive. Given the small sample size, this study could be considered a pilot for future larger studies on the role of predictive markers of CIN 2 evolution.

Relative frequency of knowledge of results and task complexity in the motor skill acqusition

The aim of this study was to investigate the effects of knowledge of results (KR) frequency and task complexity on motor skill acquisition. The task consisted of throwing a bocha ball to place it as close as possible to the target ball. 120 students ages 11 to 73 years were assigned to one of eight experimental groups according to knowledge of results frequency (25, 50, 75, and 100%) and task complexity (simple and complex). Subjects performed 90 trials in the acquisition phase and 10 trials in the transfer test. The results showed that knowledge of results given at a frequency of 25% resulted in an inferior absolute error than 50% and inferior variable error than 50, 75, and 100 I frequencies, but no effect of task complexity was found.