Using a Combination of Phylogenetic Comparative and Paleobotanical Methods to Elucidate Patterns of Lineage Selection in Rosales (Plantae: Anthophyta)

An increasing focus in evolutionary biology is on the interplay between mesoscale ecological and evolutionary processes such as population demographics, habitat tolerance, and especially geographic distribution, as potential drivers responsible for patterns of diversification and extinction over geologic time. However, few studies to date connect organismal processes such as survival and reproduction through mesoscale patterns to long-term macroevolutionary trends. In my dissertation, I investigate how mechanism of seed dispersal, mediated through geographic range size, influences diversification rates in the Rosales (Plantae: Anthophyta). In my first chapter, I validate the phylogenetic comparative methods that I use in my second and third chapters. Available state speciation and extinction (SSE) models assumptions about evolution known to be false through fossil data. I show, however, that as long as net diversification rates remain positive – a condition likely true for the Rosales – these violations of SSE’s assumptions do not cause significantly biased results. With SSE methods validated, my second chapter reconstructs three associations that appear to increase diversification rate for Rosalean genera: (1) herbaceous habit; (2) a three-way interaction combining animal dispersal, high within-genus species richness, and geographic range on multiple continents; (3) a four-way interaction combining woody habit with the other three characteristics of (2). I suggest that the three- and four-way interactions represent colonization ability and resulting extinction resistance in the face of late Cenozoic climate change; however, there are other possibilities as well that I hope to investigate in future research. My third chapter reconstructs the phylogeographic history of the Rosales using both non-fossil-assisted SSE methods as well as fossil-informed traditional phylogeographic analysis. Ancestral state reconstructions indicate that the Rosaceae diversified in North America while the other Rosalean families diversified elsewhere, possibly in Eurasia. SSE is able to successfully identify groups of genera that were likely to have been ancestrally widespread, but has poorer taxonomic resolution than methods that use fossil data. In conclusion, these chapters together suggest several potential causal links between organismal, mesoscale, and geologic scale processes, but further work will be needed to test the hypotheses that I raise here.

Expression of HOXB2, a retinoic acid signaling target in pancreatic cancer and pancreatic intraepithelial neoplasia

Purpose: Despite significant progress in understanding the molecular pathology of pancreatic cancer and its precursor lesion: pancreatic intraepithelial neoplasia (PanIN), there remain no molecules with proven clinical utility as prognostic or therapeutic markers. Here, we used oligonucleotide microarrays to interrogate mRNA expression of pancreatic cancer tissue and normal pancreas to identify novel molecular pathways dysregulated in the development and progression of pancreatic cancer. Experimental Design: RNA was hybridized to Affymetrix Genechip HG-U133 oligonucleotide microarrays. A relational database integrating data from publicly available resources was created to identify candidate genes potentially relevant to pancreatic cancer. The protein expression of one candidate, homeobox B2 (HOXB2), in PanIN and pancreatic cancer was assessed using immunohistochemistry. Results: We identified aberrant expression of several components of the retinoic acid (RA) signaling pathway (RARa, MUC4, Id-1, MMP9, uPAR, HB-EGF, HOXB6, and HOXB2), many of which are known to be aberrantly expressed in pancreatic cancer and Pan IN. HOXB2, a downstream target of RA, was up-regulated 6.7-fold in pancreatic cancer compared with normal pancreas. Immunohistochemistry revealed ectopic expression of HOXB2 in 15% of early Pan IN lesions and 48 of 128 (38%) pancreatic cancer specimens. Expression of HOXB2 was associated with nonresectable tumors and was an independent predictor of poor survival in resected tumors. Conclusions: We identified aberrant expression of RA signaling components in pancreatic cancer, including HOXB2, which was expressed in a proportion of PanIN lesions. Ectopic expression of HOXB2 was associated with a poor prognosis for all patients with pancreatic cancer and was an independent predictor of survival in patients who underwent resection.

Thymic lineage commitment rather than selection causes genetic variations in size of CD4 and CD8 compartments.

During their development, immature CD4+ CD8+ thymocytes become committed to either the CD4 or CD8 lineage. Subsequent complete maturation of CD4+ and CD8+ cells requires a molecular match of the expressed coreceptor and the MHC specificity of the TCR. The final size of the mature CD4+ and CD8+ thymic compartments is therefore determined by a combination of lineage commitment and TCR-mediated selection. In humans and mice, the relative size of CD4+ and CD8+ peripheral T cell compartments shows marked genetic variability. We show here that genetic variations in thymic lineage commitment, rather than TCR-mediated selection processes, are responsible for the distinct CD4/CD8 ratios observed in common inbred mouse strains. Genetic variations in the regulation of lineage commitment open new ways to analyze this process and to identify the molecules involved.

K-r/K-c but not d(N)/d(S) correlates positively with body mass in birds, raising implications for inferring lineage-specific selection

Background: The ratio of the rates of non-synonymous and synonymous substitution (d(N)/d(S)) is commonly used to estimate selection in coding sequences. It is often suggested that, all else being equal, d(N)/d(S) should be lower in populations with large effective size (Ne) due to increased efficacy of purifying selection. As N-e is difficult to measure directly, life history traits such as body mass, which is typically negatively associated with population size, have commonly been used as proxies in empirical tests of this hypothesis. However, evidence of whether the expected positive correlation between body mass and d(N)/d(S) is consistently observed is conflicting. Results: Employing whole genome sequence data from 48 avian species, we assess the relationship between rates of molecular evolution and life history in birds. We find a negative correlation between dN/dS and body mass, contrary to nearly neutral expectation. This raises the question whether the correlation might be a method artefact. We therefore in turn consider non-stationary base composition, divergence time and saturation as possible explanations, but find no clear patterns. However, in striking contrast to d(N)/d(S), the ratio of radical to conservative amino acid substitutions (K-r/K-c) correlates positively with body mass. Conclusions: Our results in principle accord with the notion that non-synonymous substitutions causing radical amino acid changes are more efficiently removed by selection in large populations, consistent with nearly neutral theory. These findings have implications for the use of d(N)/d(S) and suggest that caution is warranted when drawing conclusions about lineage-specific modes of protein evolution using this metric.

Alpha beta lineage-committed thymocytes can be rescued by the gamma delta T cell receptor (TCR) in the absence of TCR beta chain.

Commitment of the alpha beta and gamma delta T cell lineages within the thymus has been studied in T cell receptor (TCR)-transgenic and TCR mutant murine strains. TCR gamma delta-transgenic or TCR beta knockout mice, both of which are unable to generate TCR alpha beta-positive T cells, develop phenotypically alpha beta-like thymocytes in significant proportions. We provide evidence that in the absence of functional TCR beta protein, the gamma delta TCR can promote the development of alpha beta-like thymocytes, which, however, do not expand significantly and do not mature into gamma delta T cells. These results show that commitment to the alpha beta lineage can be determined independently of the isotype of the TCR, and suggest that alpha beta versus gamma delta T cell lineage commitment is principally regulated by mechanisms distinct from TCR-mediated selection. To accommodate our data and those reported previously on the effect of TCR gamma and delta gene rearrangements on alpha beta T cell development, we propose a model in which lineage commitment occurs independently of TCR gene rearrangement.

Development and selection of NKT cells.

NKT cells utilize a restricted alphabeta TCR repertoire that recognizes glycolipids in association with CD1d. The recent development of fluorescent CD1d tetramers loaded with the synthetic glycolipid alpha-galactosyl-ceramide has led to a clearer definition of NKT-cell subsets as well as important insights into their developmental origin. As many as four subsets may exist, differing in NK1.1 expression, TCR repertoire and dependence on CD1d and various glycolipids for development. Two different lineage-commitment models have been proposed, with most evidence favoring a byproduct of conventional-T-cell development.

A critical role for the T cell receptor alpha-chain connecting peptide domain in positive selection of CD1-independent NKT cells.

Natural killer T (NKT) cells are a subset of mature alpha beta TCR(+) cells that co-express NK lineage markers. Whereas most NKT cells express a canonical Valpha14/Vbeta8.2 TCR and are selected by CD1d, a minority of NKT cells express a diverse TCR repertoire and develop independently of CD1d. Little is known about the selection requirements of CD1d-independent NKT cells. We show here that NKT cells develop in RAG-deficient mice expressing an MHC class II-restricted transgenic TCR (Valpha2/Vbeta8.1) but only under conditions that lead to negative selection of conventional T cells. Moreover development of NKT cells in these mice is absolutely dependent upon an intact TCR alpha-chain connecting peptide domain, which is required for positive selection of conventional T cells via recruitment of the ERK signaling pathway. Collectively our data demonstrate that NKT cells can develop as a result of high avidity TCR/MHC class II interactions and suggest that common signaling pathways are involved in the positive selection of CD1d-independent NKT cells and conventional T cells.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

The role of the T-cell receptor (TCR) in alpha beta/gamma delta lineage commitment: clues from intracellular TCR staining.

T cells belong to two mutually exclusive lineages expressing either alpha beta or gamma delta T-cell receptors (TCR). Although alpha beta and gamma delta cells are known to share a common precursor the role of TCR rearrangement and specificity in the lineage commitment process is controversial. Instructive lineage commitment models endow the alpha beta or gamma delta TCR with a deterministic role in lineage choice, whereas separate lineage models invoke TCR-independent lineage commitment followed by TCR-dependent selection and maturation of alpha beta and gamma delta cells. Here we review the published data pertaining to the role of the TCR in alpha beta/gamma delta lineage commitment and provide some additional information obtained from recent intracellular TCR staining studies. We conclude that a variant of the separate lineage model is best able to accommodate all of the available experimental results.

A limited role for beta-selection during gamma delta T cell development.

T cells belong to two distinct lineages expressing either alpha beta or gamma delta TCR. During alpha beta T cell development, it is clearly established that productive rearrangement at the TCR beta locus in immature precursor cells leads to the expression of a pre-TCR complex. Signaling through the pre-TCR results in the selective proliferation and maturation of TCR beta+ cells, a process that is known as beta-selection. However, the potential role of beta-selection during gamma delta T cell development is controversial. Whereas PCR-RFLP and sequencing techniques have provided evidence for a bias toward in-frame VDJ beta rearrangements in gamma delta cells (consistent with beta-selection), gamma delta cells apparently develop normally in mice that are unable to assemble a pre-TCR complex due to a deficiency in TCR beta or pT alpha genes. In this report, we have directly addressed the physiologic significance of beta-selection during gamma delta cell development in normal mice by quantitating intracellular TCR beta protein in gamma delta cells and correlating its presence with cell cycle status. Our results indicate that beta-selection plays a significant (although limited) role in gamma delta cell development by selectively amplifying a minor subset of gamma delta precursor cells with productively rearranged TCR beta genes.

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection.

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species' history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice-free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus, we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum-likelihood migration rates. The 'two refugia' scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine-scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.

Patterns of positive selection in seven ant genomes.

The evolution of ants is marked by remarkable adaptations that allowed the development of very complex social systems. To identify how ant-specific adaptations are associated with patterns of molecular evolution, we searched for signs of positive selection on amino-acid changes in proteins. We identified 24 functional categories of genes which were enriched for positively selected genes in the ant lineage. We also reanalyzed genome-wide data sets in bees and flies with the same methodology to check whether positive selection was specific to ants or also present in other insects. Notably, genes implicated in immunity were enriched for positively selected genes in the three lineages, ruling out the hypothesis that the evolution of hygienic behaviors in social insects caused a major relaxation of selective pressure on immune genes. Our scan also indicated that genes implicated in neurogenesis and olfaction started to undergo increased positive selection before the evolution of sociality in Hymenoptera. Finally, the comparison between these three lineages allowed us to pinpoint molecular evolution patterns that were specific to the ant lineage. In particular, there was ant-specific recurrent positive selection on genes with mitochondrial functions, suggesting that mitochondrial activity was improved during the evolution of this lineage. This might have been an important step toward the evolution of extreme lifespan that is a hallmark of ants.

Development and selection of Valpha l4i NKT cells.

Valpha14 invariant natural killer T (Valpha14i NKT) cells are a unique lineage of mouse T cells that share properties with both NK cells and memory T cells. Valpha14i NKT cells recognize CDld-associated glycolipids via a semi-invariant T cell receptor (TCR) composed of an invariant Valpha14-Jalpha 18 chain paired preferentially with a restricted set of TCRbeta chains. During development in the thymus, rare CD4+ CD8+ (DP) cortical thymocytes that successfully rearrange the semi-invariant TCR are directed to the Valpha14i NKT cell lineage via interactions with CD d-associated endogenous glycolipids expressed by other DP thymocytes. As they mature, Valphal4i NKT lineage cells upregulate activation markers such as CD44 and subsequently express NK-related molecules such as NKI.1 and members of the Ly-49 inhibitory receptor family. The developmental program of Valpha l4i NKT cells is critically regulated by a number of signaling cues that have little or no effect on conventional T cell development, such as the Fyn/SAP/SLAM pathway, NFkappaB and T-bet transcription factors, and the cytokine IL-15. The unique developmental requirements of Valphal4i NKT cells may represent a paradigm for other unconventional T cell subsets that are positively selected by agonist ligands expressed on hematopoietic cells.

Interrogating eleven fast-evolving genes for signatures of recent positive selection in worldwide human populations

Different signatures of natural selection persist over varying time scales in our genome, revealing possible episodes of adaptative evolution during human history. Here, we identify genes showing signatures of ancestral positive selection in the human lineage and investigate whether some of those genes have been evolving adaptatively in extant human populations. Specifically, we compared more than 11,000 human genes with their orthologs inchimpanzee, mouse, rat and dog and applied a branch-site likelihood method to test for positive selection on the human lineage. Among the significant cases, a robust set of 11 genes were then further explored for signatures of recent positive selection using SNP data. We genotyped 223 SNPs in 39 worldwide populations from the HGDP Diversity panel and supplemented this information with available genotypes for up to 4,814 SNPs distributed along 2 Mb centered on each gene. After exploring the allele frequency spectrum, population differentiation and the maintainance of long unbroken haplotypes, we found signals of recent adaptative phenomena in only one of the 11 candidate gene regions. However, the signal ofrecent selection in this region may come from a different, neighbouring gene (CD5) ratherthan from the candidate gene itself (VPS37C). For this set of positively-selected genes in thehuman lineage, we find no indication that these genes maintained their rapid evolutionarypace among human populations. Based on these data, it therefore appears that adaptation forhuman-specific and for population-specific traits may have involved different genes.

Application of RNAi to silence tartrate-resistant acid phosphatase: unexpected effects on the monocyte-macrophage lineage

Tartraatti-resistentin happaman fosfataasin hiljentäminen RNAi menetelmällä: odottamaton vaikutus monosyytti-makrofagi linjan soluissa RNA interferenssi (RNAi) eli RNA:n hiljentyminen löydettiin ensimmäisenä kasveissa, ja 2000-luvulla RNAi menetelmä on otettu käyttöön myös nisäkässoluissa. RNAi on mekanismi, jossa lyhyet kaksi juosteiset RNA molekyylit eli siRNA:t sitoutuvat proteiinikompleksiin ja sitoutuvat komplementaarisesti proteiinia koodaavaan lähetti RNA:han katalysoiden lähetti RNA:n hajoamisen. Tällöin RNA:n koodaamaa proteiinia ei solussa tuoteta. Tässä työssä on RNA interferenssi menetelmän avuksi kehitetty uusi siRNA molekyylien suunnittelualgoritmi siRNA_profile, joka etsii lähetti RNA:sta geenin hiljentämiseen sopivia kohdealueita. Optimaalisesti suunnitellulla siRNA molekyylillä voi olla mahdollista saavuttaa pitkäaikainen geenin hiljeneminen ja spesifinen kohdeproteiinin määrän aleneminen solussa. Erilaiset kemialliset modifikaatiot, mm. 2´-Fluoro-modifikaatio, siRNA molekyylin riboosirenkaassa lisäsivät siRNA molekyylin stabiilisuutta veren plasmassa sekä siRNA molekyylin tehokkuutta. Nämä ovat tärkeitä siRNA molekyylien ominaisuuksia kun RNAi menetelmää sovelletaan lääketieteellisiin tarkoituksiin. Tartraatti-resistentti hapan fosfataasi (TRACP) on entsyymi, joka esiintyy luunsyöjäsoluissa eli osteoklasteissa, antigeenejä esittelevissä dendiriittisissä soluissa sekä eri kudosten makrofageissa, jotka ovat syöjäsoluja. TRACP entsyymin biologista tehtävää ei ole saatu selville, mutta oletetaan että TRACP entsyymin kyvyllä tuottaa reaktiivisia happiradikaaleja on tehtävä sekä luuta hajoittavissa osteoklasteissa sekä antigeenia esittelevissä dendriittisissä soluissa. Makrofageilla, jotka yliekpressoivat TRACP entsyymiä, on myös solunsisäinen reaktiivisten happiradikaalien tuotanto sekä bakteerin tappokyky lisääntynyt. TRACP-geenin hiljentämiseen tarkoitetut spesifiset DNA ja siRNA molekyylit aiheuttivat monosyytti-makrofagilinjan soluviljelymallissa TRACP entsyymin tuoton lisääntymistä odotusten vastaisesti. DNA ja RNA molekyylien vaikutusta TRACP entsyymin tuoton lisääntymiseen tutkittiin myös Tolllike reseptori 9 (TLR9) poistogeenisestä hiirestä eristetyissä monosyyttimakrofaagisoluissa. TRACP entsyymin tuoton lisääntyminen todettiin sekvenssistä ja TLR9:stä riippumattomaksi vasteeksi solun ulkopuolisia DNA ja RNA molekyylejä vastaan. Havainto TRACP entsyymin tuoton lisääntymisestä viittaa siihen, että TRACP entsyymillä on tehtävä solun immuunipuolustusjärjestelmässä.